Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

Elastic Flow Modeling for Hydropower Digital Twins

This report details the elastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The elastic flow model is accurate even in cases of fast transients, long penstock length, and high gravitation head. The compressibility of the water and elasticity of the pipe walls are explicitly accounted for in the developed models so that the water hammer phenomena can be accurately captured. The elastic flow model is coupled to a mechanistic turbine model to resolve the dynamic feedbacks between the elastic water column and the turbine rotation rate. A finite volume method is employed to solve the governing equations, and the method is shown to have low numerical diffusion in sharp gradient phenomena, such as those encountered when simulating water hammer. The method is applied to single dimensional linear advection benchmark problem and numerical results are compared with analytic results. This is followed on by an application to a full hydropower system with a coupled turbine. The elastic flow model results are compared against an inelastic model.

13 HYDRO ENERGY↗

Understanding the Physics Representation of Deep Learning Models in Environmental Applications

Deep learning (DL) models have been popular in earth and environmental modeling and analysis, which exhibit huge potential in capturing and reconstructing the non-linearity of relevant environmental processes. They are extensively used as analytical tools or emulators for multiple domains (atmosphere, land surface, ocean, and biogeochemistry). Despite their success, their internal working mechanism remains largely unknown. Such a lack of knowledge hinders the identification of physically consistent models that are fully adaptive to non-stationary climate, as well as the development of physics-informed machine learning such as physics-informed neural network (PINN). To establish preliminary knowledge and framework of such physics representation evaluation, this project focuses on an improved understanding of DL models in the environmental applications. DL models are increasingly applied to environmental modeling and prediction. However, they have been evaluated mostly from a performance perspective, and there is a gap in understanding how they represent the known physics internally. Such knowledge is especially critical when applying DL models under climate change conditions, where new inputs are likely outside the ranges of the training datasets. In this project, we reveal how the known physical processes are represented within DL models from both statistical and mechanistic perspectives. Leveraging the traditional model evaluations that focus more on the accuracies of predictions, we establish a framework that examines both the accuracy and physics representation of DL models. This analysis framework can identify DL models that make the correct predictions based on correct physics, thus enhancing the existing explainable artificial intelligence (explainable-AI) portfolio. It lays a foundation for developing novel metrics to evaluate the emerging DL models in environmental applications. This knowledge also informs the development of physics-informed DL models by revealing the direct connections between the known physical processes and specific model components or structures.

54 ENVIRONMENTAL SCIENCES↗

The kinetic Sunyaev-Zel'dovich effect from reionization: simulated full-sky maps at arcminute resolution

The kinetic Sunyaev–Zel’dovich (kSZ) effect results from Thomson scattering by coherent flows in the reionized intergalactic medium. We present new results based on ray-tracing an 8 Gpc/h realization of reionization with resolution elements 2 Mpc/h (subtending ∼1′ at z = 6) on a side to create a full-sky kSZ map. The realization includes, self-consistently, the effects of reionization on scales corresponding to multipoles 10≲ℓ≲5000. We separate the kSZ map into Doppler (v), Ostriker–Vishniac (δv), patchy (xv), and third-order (xδv) components, and compute explicitly all the auto- and cross-correlations (e.g., 〈vv〉, 〈δvxv〉, etc.) that contribute to the total power. We find a complex and nonmonotonic dependence on the duration of reionization at ℓ∼300 and evidence for a non-negligible (10%–30%) contribution from connected four-point correlations, 〈xvxv〉{sub c}, usually neglected in analytical models. We also investigate the cross-correlation of linear matter and large-scale kSZ temperature fluctuations, focusing on (1) cross-power spectra with biased tracers of the matter density and (2) cold spots from infall onto large, rare H ii regions centered on peaks in the matter distribution at redshifts z>10 that are a generic non-Gaussian feature of patchy reionization. Finally, we show that the reionization history can be reconstructed at 5σ–10σ significance by correlating full-sky 21 cm maps stacked in bins with Δν = 10 MHz with existing cosmic microwave background (CMB) temperature maps at ℓ<500, raising the prospects for probing reionization by correlating CMB and LSS measurements. The resulting kSZ maps have been made publicly available at www.cita.utoronto.ca/∼malvarez/research/ksz-data/.

79 ASTRONOMY AND ASTROPHYSICS↗

KELVIN–HELMHOLTZ INSTABILITY IN SOLAR CHROMOSPHERIC JETS: THEORY AND OBSERVATION

Using data obtained by the high-resolution CRisp Imaging SpectroPolarimeter instrument on the Swedish 1 m Solar Telescope, we investigate the dynamics and stability of quiet-Sun chromospheric jets observed at the disk center. Small-scale features, such as rapid redshifted and blueshifted excursions, appearing as high-speed jets in the wings of the H α line, are characterized by short lifetimes and rapid fading without any descending behavior. To study the theoretical aspects of their stability without considering their formation mechanism, we model chromospheric jets as twisted magnetic flux tubes moving along their axis, and use the ideal linear incompressible magnetohydrodynamic approximation to derive the governing dispersion equation. Analytical solutions of the dispersion equation indicate that this type of jet is unstable to Kelvin–Helmholtz instability (KHI), with a very short (few seconds) instability growth time at high upflow speeds. The generated vortices and unresolved turbulent flows associated with the KHI could be observed as a broadening of chromospheric spectral lines. Analysis of the H α line profiles shows that the detected structures have enhanced line widths with respect to the background. We also investigate the stability of a larger-scale H α jet that was ejected along the line of sight. Vortex-like features, rapidly developing around the jet’s boundary, are considered as evidence of the KHI. The analysis of the energy equation in the partially ionized plasma shows that ion–neutral collisions may lead to fast heating of the KH vortices over timescales comparable to the lifetime of chromospheric jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal analytic multiburn trajectories

Fuel optimal analytic multiburn transfer trajectories and shuttle rendezvous, assuming gravity vector linear dependence on position vector for burn arcs

Jezewski, D. J.↗

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations. Phase 5B: Analysis of gust alleviation methods and rotor dynamic stability

The effects of various gust alleviation methods on the random blade response in flapping were studied analytically, assuming a rigid rotor support. The analytical model assumes rigid flapping blades with elastic root restraints. Linearized equations which are approximately valid at low lift conditions were used. Because of the interblade coupling from the feedback devices, the method of multiblade generalized coordinates was most convenient and was extended to include coning, differential coning, and warping of the rotor. The numerical examples cover the dynamic stability characteristics as affected by feedback gains of three- to six-bladed rotors. The number of blades had large effects on stability limits and modal time functions at these limits. The random flapping response of the blades to atmospheric turbulence was determined at 1.6 rotor advance ratio using feedback gains below the stability limit. The most effective reduction of the flapping response per unit gain was achieved with a rotor coning angle feedback.

Hohenemser, K. H.↗

Sound generated in a cascade by three-dimensional disturbances convected in a subsonic flow

Discrete tone sound generation in a subsonic fan subject to three-dimensional disturbances is investigated. The analytical model used treats the fan rotor and stator as linear cascades of thin airfoils in a rectangular duct subject to a three-dimensional gust for which a complete aerodynamic theory already exists. The sound pressure can then be cast as the sum of a finite number of discrete sound waves (modes) the magnitude of which depends on an unknown function satisfying a singular integral equation. Similarity rules are derived to reduce the problem to that of a two-dimensional gust. Three-dimensional effects on the cut-off condition, the sound pressure, and the acoustic power are first investigated for each mode. The theory is then applied to noise generated by typical rotor-wake-defect and rotor-tip-vortex disturbances interacting with a stator.

Atassi, H.↗

Energetics and the resistive tearing mode - Effects of Joule heating and radiation

The contribution of energy flux to the dynamics of magnetic field reconnection is analytically studied in order to determine the influence of Joule heating and radiation on the linear development of the tearing instability in slab geometry. A temperature-dependent Coulomb-like resistivity is used to provide the coupling between the dynamics and the energy equation. Analytical expressions are derived for the growth rates utilizing constant-psi and long-wavelength approximations. The solutions indicate the occurrence of several modes in addition to the usual tearing mode, several of which have relatively slow, complex growth rates. At large values of the magnetic Reynolds number, there are at least two modes with purely exponential growth when the radiative loss decreases with increasing temperature. If the radiation is neglected, the Joule heating alone also results in two modes with real, positive growth at large S. Below a particular value of S, all the modes are generally stabilized.

Steinolfson, R. S.↗

Tether damping in space

This paper analyzes the dynamics of a long tether connecting two spacecraft in earth orbit, one of the spacecraft having dominant mass. In particular, it considers the material damping of the tether. The nominal position of the tether is stabilized by the gravity gradient such that it is aligned with the local vertical. The tether is modeled as a viscoelastic flexible continuum. Modal frequencies are derived in an analytical approximation form. Damping ratios are estimated according to the linear model calibrated by ground measurements. The results show that, with properly chosen tether material and braiding structure, longitudinal vibrations of the tethered system are well damped.

Powell, J. David↗

Combined control-structure optimization

An approach for combined control-structure optimization keyed to enhancing early design trade-offs is outlined and illustrated by numerical examples. The approach employs a homotopic strategy and appears to be effective for generating families of designs that can be used in these early trade studies. Analytical results were obtained for classes of structure/control objectives with linear quadratic Gaussian (LQG) and linear quadratic regulator (LQR) costs. For these, researchers demonstrated that global optima can be computed for small values of the homotopy parameter. Conditions for local optima along the homotopy path were also given. Details of two numerical examples employing the LQR control cost were given showing variations of the optimal design variables along the homotopy path. The results of the second example suggest that introducing a second homotopy parameter relating the two parts of the control index in the LQG/LQR formulation might serve to enlarge the family of Pareto optima, but its effect on modifying the optimal structural shapes may be analogous to the original parameter lambda.

Salama, M.↗

Vortex breakdown incipience: Theoretical considerations

The sensitivity of the onset and the location of vortex breakdowns in concentrated vortex cores, and the pronounced tendency of the breakdowns to migrate upstream have been characteristic observations of experimental investigations; they have also been features of numerical simulations and led to questions about the validity of these simulations. This behavior seems to be inconsistent with the strong time-like axial evolution of the flow, as expressed explicitly, for example, by the quasi-cylindrical approximate equations for this flow. An order-of-magnitude analysis of the equations of motion near breakdown leads to a modified set of governing equations, analysis of which demonstrates that the interplay between radial inertial, pressure, and viscous forces gives an elliptic character to these concentrated swirling flows. Analytical, asymptotic, and numerical solutions of a simplified non-linear equation are presented; these qualitatively exhibit the features of vortex onset and location noted above.

Berger, Stanley A.↗

The Middeck Active Control Experiment (MACE)

The Middeck Active Control Experiment (MACE) is a NASA In-Step and Control Structure Interaction (CSI) Office funded Shuttle middeck experiment. The objective is to investigate the extent to which closed-loop behavior of flexible spacecraft in zero-gravity (0-g) can be predicted. This prediction becomes particularly difficult when dynamic behavior during ground testing exhibits extensive suspension and direct gravity coupling. On-orbit system identification and control reconfiguration is investigated to improve performance which would otherwise be limited due to errors in prediction. The program is presently in its preliminary design phase with launch expected in the summer of 1994. The MACE test article consists of three attitude control torque wheels, a two axis gimballing payload, inertial sensors and a flexible support structure. With the acquisition of a second payload, this will represent a multiple payload platform with significant structural flexibility. This paper presents on-going work in the areas of modelling and control of the MACE test article in the zero and one-gravity environments. Finite element models, which include suspension and gravity effects, and measurement models, derived from experimental data, are used as the basis for Linear Quadratic Gaussian controller designs. Finite element based controllers are analytically used to study the differences in closed-loop performance as the test article transitions between the 0-g and 1-g environments. Measurement based controllers are experimentally applied to the MACE test article in the 1-g environment and achieve over an order of magnitude improvement in payload pointing accuracy when disturbed by a broadband torque disturbance. The various aspects of the flight portion of the experiment are also discussed.

Miller, David W.↗

Active control of panel vibrations induced by a boundary layer flow

The problems of active and passive control of sound and vibration has been investigated by many researchers for a number of years. However, few of the articles are concerned with the sound and vibration with flow-structure interaction. Experimental and numerical studies on the coupling between panel vibration and acoustic radiation due to flow excitation have been done by Maestrello and his associates at NASA/Langley Research Center. Since the coupled system of nonlinear partial differential equations is formidable, an analytical solution to the full problem seems impossible. For this reason, we have to simplify the problem to that of the nonlinear panel vibration induced by a uniform flow or a boundary-layer flow with a given wall pressure distribution. Based on this simplified model, we have been able to consider the control and stabilization of the nonlinear panel vibration, which have not been treated satisfactorily by other authors. Although the sound radiation has not been included, the vibration suppression will clearly reduce the sound radiation power from the panel. The major research findings are presented in three sections. In section two we describe results on the boundary control of nonlinear panel vibration, with or without flow excitation. Sections three and four are concerned with some analytical and numerical results in the optimal control of the linear and nonlinear panel vibrations, respectively, excited by the flow pressure fluctuations. Finally, in section five, we draw some conclusions from research findings.

Chow, Pao-Liu↗

Control-Theory-Informed Feature Selection for Detecting Malicious Tampering in Additive Layer Manufacturing Processes

Additive layer manufacturing (ALM) is rapidly becoming an appealing solution to the low-volume manufacturing of metal, polymer, or composite parts. However, ALM’s reliance on digital part specifications, microcontrollers, and modern networking makes these devices vulnerable to malicious tampering by cyber attackers, which can negatively affect part performance and even result in catastrophic failure. We present a hybrid analytic approach to feature discovery using control theoretic techniques and linear modelling on input-output data collected from a representative controller system. Employing this approach, we design, train, and test an anomaly detection system. The preliminary results show that the proposed approach effectively discovers useful input-output relationships for anomaly detection in a simulated ALM process. Application to larger and more complex systems are discussed.

Dawson, Joel↗

Formulation of the linear model from the nonlinear simulation for the F18 HARV

The F-18 HARV is a modified F-18 Aircraft which is capable of flying in the post-stall regime in order to achieve superagility. The onset of aerodynamic stall, and continued into the post-stall region, is characterized by nonlinearities in the aerodynamic coefficients. These aerodynamic coefficients are not expressed as analytic functions, but rather in the form of tabular data. The nonlinearities in the aerodynamic coefficients yield a nonlinear model of the aircraft's dynamics. Nonlinear system theory has made many advances, but this area is not sufficiently developed to allow its application to this problem, since many of the theorems are existance theorems and that the systems are composed of analytic functions. Thus, the feedback matrices and the state estimators are obtained from linear system theory techniques. It is important, in order to obtain the correct feedback matrices and state estimators, that the linear description of the nonlinear flight dynamics be as accurate as possible. A nonlinear simulation is run under the Advanced Continuous Simulation Language (ACSL). The ACSL simulation uses FORTRAN subroutines to interface to the look-up tables for the aerodynamic data. ACSL has commands to form the linear representation for the system. Other aspects of this investigation are discussed.

Hall, Charles E., Jr.↗

TH-C-17A-07: Visualizing and Quantifying Radiation Therapy in Real-Time Using a Novel Beam Imaging Technique

Purpose: Radiation therapy uses invisible high energy X-rays to treat an invisible tumor. Proper positioning of the treatment beam relative to the patient's anatomy during dose delivery is critically important to the success of treatment. We develop and characterize a novel radiation therapy beam visualization technique for real-time monitoring of patient treatment. Methods: Custom made flexible scintillator sheets were fabricated from gadolinium oxysulfide (GOS) particles that had been doped with terbium within a silicone elastomer matrix. Sheets of several thicknesses ranging from 0.3 to 1mm were prepared and tested. Sheets were exposed to megavoltage X-ray and electron beams from a Varian linac and the resulting optical signal was collected by multiple CMOS cameras placed in the treatment room. Real-time images were collected for different beam energies and dose rates. Signal intensity and SNR were calculated by processing the acquired images. Results: All signals were detectable in the presence of full room lighting and at an integration time of 45ms. Average signal intensity and SNR increased with both sheet thickness and dose rate and decreased with beam energy and incident light. For a given sheet thickness and beam energy the correlation between dose rate and signal intensity was highly linear. Increased sheet thickness or dose rate results in a linear increase in the detected signal. All results are consistent with analytical approximations. Conclusion: The technique offers a means of accurately visualizing a radiation therapy beam shape and fluence in real time. The effects of salient parameters have been characterized and will enable further optimization of the technique as it is implemented into the clinical workflow. The project described was supported by the National Center for Advancing Translational Sciences of the National Institutes of Health through UL1 TR001085.

60 APPLIED LIFE SCIENCES↗

TLD efficiency calculations for heavy ions: an analytical approach

The use of thermoluminescent dosimeters (TLDs) in heavy charged particles’ dosimetry is limited by their non-linear dose response curve and by their response dependence on the radiation quality. Thus, in order to use TLDs with particle beams, a model that can reproduce the behavior of these detectors under different conditions is needed. Here a new, simple and completely analytical algorithm for the calculation of the relative TL-efficiency depending on the ion charge Z and energy E is presented. In addition, the detector response is evaluated starting from the single ion case, where the computed effectiveness values have been compared with experimental data as well as with predictions from a different method. The main advantage of this approach is that, being fully analytical, it is computationally fast and can be efficiently integrated into treatment planning verification tools. In conclusion, the calculated efficiency values have been then implemented in the treatment planning code TRiP98 and dose calculations on a macroscopic target irradiated with an extended carbon ion field have been performed and verified against experimental data.

43 PARTICLE ACCELERATORS↗

Stability with large step sizes for multistep discretizations of stiff ordinary differential equations

One-leg and multistep discretizations of variable-coefficient linear systems of ODEs having both slow and fast time scales are investigated analytically. The stability properties of these discretizations are obtained independent of ODE stiffness and compared. The results of numerical computations are presented in tables, and it is shown that for large step sizes the stability of one-leg methods is better than that of the corresponding linear multistep methods.

Majda, George↗