Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear processes”

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 1,441 records · Page 80

Iterative Nonlinear Tikhonov Algorithm with Constraints for Electromagnetic Tomography

Low frequency electromagnetic tomography such as the capacitance tomography (ECT) has been proposed for monitoring and mass-gauging of gas-liquid two-phase system under microgravity condition in NASA's future long-term space missions. Due to the ill-posed inverse problem of ECT, images reconstructed using conventional linear algorithms often suffer from limitations such as low resolution and blurred edges. Hence, new efficient high resolution nonlinear imaging algorithms are needed for accurate two-phase imaging. The proposed Iterative Nonlinear Tikhonov Regularized Algorithm with Constraints (INTAC) is based on an efficient finite element method (FEM) forward model of quasi-static electromagnetic problem. It iteratively minimizes the discrepancy between FEM simulated and actual measured capacitances by adjusting the reconstructed image using the Tikhonov regularized method. More importantly, it enforces the known permittivity of two phases to the unknown pixels which exceed the reasonable range of permittivity in each iteration. This strategy does not only stabilize the converging process, but also produces sharper images. Simulations show that resolution improvement of over 2 times can be achieved by INTAC with respect to conventional approaches. Strategies to further improve spatial imaging resolution are suggested, as well as techniques to accelerate nonlinear forward model and thus increase the temporal resolution.

Xu, Feng↗

Quantum kinetics of anomalous and nonlinear Hall effects in topological semimetals

Highlights: • Linear and photogalvanic anomalous Hall responses are systematically derived. • Extrinsic mechanisms of AHE include Gaussian, diffractive, hybrid skew scattering. • Diagrammatic calculations are matched to semiclassical picture of AHE. • The Pancharatnam phase of multifold fermions determines the skew scattering amplitude. • Photon-induced interband scattering are accompanied by coordinate shifts. We present a systematic microscopic derivation of the semiclassical Boltzmann equation for band structures with the finite Berry curvature based on Keldysh technique of nonequilibrium systems. In the analysis, an AC electrical driving field is kept up to quadratic order, and both cases of small and large frequencies corresponding to intra- and interband transitions are considered. In particular, this formulation is suitable for the study of nonlinear Hall effect and photogalvanic phenomena. The role of impurity scattering is carefully addressed. Specifically, in addition to previously studied side-jump and skew-scattering processes, quantum interference diffractive contributions are now explicitly incorporated within the developed framework. This theory is applied to multifold fermions in topological semimetals, for which the generic formula for the skew scattering rate from the Pancharatnam phase is obtained along with the corresponding anomalous Hall conductivity.

36 MATERIALS SCIENCE↗

An experimental study on round-trip efficiency of a preheating control with a medium office building

The interactions between building energy efficiency and demand response (DR) with buildings may complicate their impacts on the electricity grid and buildings, yet those interactions are not well understood. This work fills this gap by presenting an experimental study that quantifies the effects of DR on building energy efficiency. Unlike existing studies, this work 1) investigates the uncertainty levels when assessing the effects of DR; 2) studies the impact of rooftop heat pumps with direct expansion refrigerant-to-air coils on the round-trip efficiency. Specifically, we considered a medium office building (~1860 m 2 ), which is served by rooftop units, located in eastern Washington State. We performed preheating control tests during typical heating days in the year 2022 and compared the building’s energy consumption during the testing days to that during normal operation days under similar weather conditions. To understand how well those normal operation days represent the baseline for the assessment, the relationship between the energy performance of this building and weather conditions was also studied. Experimental results reveal that the round-trip efficiency (i.e., the ratio of discharging energy output from the building to charging energy input into the building) was between 0.16 and 0.92. The results show that significant energy loss may occur during the preheating process. Finally, it also suggests a significant correlation between the energy efficiency of the studied building under the preheating control and the outdoor temperature; and a highly nonlinear relationship between the zone temperature changing rate and the weather condition.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗

CI-MOR Final Report: Analysis and Validation of Critical Infrastructure Models using Model Order Reduction

This report summarizes the research and capabilities developed as part of the project “Analysis and Validation of Critical Infrastructure Models using Model Order Reduction” (CI-MOR) LDRD project. CI-MOR research enables the solution of large, complex optimization models that naturally arise in national security challenges involving critical infrastructures. Specifically, CI-MOR researchers developed methods to (1) rigorously approximate complex, nonlinear optimization formulations, (2) identify alternative near-optimal solutions, (3) accelerate optimization workflows used for complex applications, and (4) rigorously integrate domain knowledge in stochastic-process models. This report provides an overview of the research done in CI-MOR, and we describe application exemplars used to illustrate CI-MOR capabilities. Furthermore, we describe the software developed by CI-MOR that researchers can leverage to analyze new applications.

97 MATHEMATICS AND COMPUTING↗

Surrogate Modeling of Landau Damping with Deep Operator Networks

Kinetic simulations excel at capturing microscale plasma physics phenomena with high accuracy, but their computational demands make them impractical for modeling large-scale space and astrophysical systems. In this context, we build a surrogate model, using Deep Operator Networks (DeepONets), based upon the Vlasov–Poisson simulation data to model the dynamical evolution of plasmas, focusing on the Landau damping process—a fundamental kinetic phenomenon in space and astrophysical plasmas. The trained DeepONets are able to capture the evolution of electric field energy in both linear and nonlinear regimes under various conditions. Extensive validation highlights DeepONets’ robust performance in reproducing complex plasma behaviors with high accuracy, paving the way for large-scale modeling of space and astrophysical plasmas.

plasma astrophysics↗

On the role of steepened structures and drift waves in equatorial spread F

A description is presented of a unified model of equatorial spread F which organizes much of the observational data and theoretical results in a self-consistent way. It is suggested that an initial bottomside instability nonlinearly evolves to the point that very steep gradients develop and that these steepened structures are responsible for the power spectra observed by probe experiments. This primary process is very likely the Rayleigh-Taylor instability. Observations indicate that some of the bottomside structures rise into the topside. It is shown that the observed gradients on the edges of these rising bubbles should be unstable to drift waves with a linear growth rate of about 1 per second and a perpendicular wavelength, for maximum growth, near the ion gyroradius.

Costa, E.↗

Dynamic properties of the thermosphere inferred from Pioneer Venus mass spectrometer measurements

The data obtained by the Pioneer Venus spectrometer experiments indicate that the day-night temperature contrast on Venus is associated with wind velocities of about 200 m/s which transport oxygen, helium, and hydrogen toward the night side. A mass exchange with the mesosphere commensurate with an eddy diffusion coefficient of 3 x 10 to the 7th is required to buffer the horizontal advection so as to reproduce the observed day time bulge in oxygen and the small diurnal variations in helium. The observed time response and magnitude of the day-night density variations require transport processes to be effective over time periods between five and ten days, implying a superrotation rate or prevailing winds in excess of 50 m/s at the equator. Nonlinear mass transport results in wave steepening and contributes to the amplification of the density extrema in hydrogen and helium.

Mayr, H. G.↗

Analysis of shell-type structures subjected to time-dependent mechanical and thermal loading

A general mathematical model and solution methodologies are being developed for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic, or static thermomechanical loads. Among the system responses, which were associated with these load conditions, were thermal buckling, creep buckling, and ratcheting. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model. Furthermore, this must also be accommodated in the solution process.

Simitses, G. J.↗

Simulator evaluation of system identification with on-line control law update for the controls and astrophysics experiment in space

A procedure for optimizing the performance of large flexible spacecraft that require active vibration suppression to achieve required performance is presented. The procedure is to conduct on-orbit testing and system identification followed by a control system design. It is applied via simulation to a spacecraft configuration currently being considered for flight test by NASA - the Controls, Astrophysics, and Structures Experiment in Space (CASES). The system simulator is based on a NASTRAN finite element structural model. A finite number of modes is used to represent the structural dynamics. The system simulator also includes models of the electronics, actuators, sensors, the digital controller, and the internal and external disturbances. Nonlinearities caused by quantization are included in the study to examine tolerance of the procedure to modelling errors. Disturbance and sensor noise is modelled as a Gaussian process. For system identification, the system is excited using sinusoidal inputs at the resonant frequencies of the structure using each actuator. Mode shapes, frequencies, and damping ratios are identified from the unforced response sensor data after each excitation. Then, the excitation data is used to identify the actuator influence coefficients. The results of the individual parameter identification analyses are assembled into an aggregate system model. The control design is accomplished based only on the identified model using multi-input/output linear quadratic Gaussian theory. Its performance is evaluated based on time-to-damp as compared with the uncontrolled structure.

Montgomery, Raymond C.↗

Plasma diffusion at the magnetopause? The case of lower hybrid drift waves

The diffusion expected from the quasilinear theory of the lower hybrid drift instability at the Earth's magnetopause is recalculated. The resulting diffusion coefficient is in principle just marginally large enough to explain the thickness of the boundary layer under quiet conditions, based on observational upper limits for the wave intensities. Thus, one possible model for the boundary layer could involve equilibrium between the diffusion arising from lower hybrid waves and various low processes. However, some recent data and simulations seems to indicate that the magnetopause is not consistent with such a soft diffusive equilibrium model. Furthermore, investigation of the nonlinear equations for the lower hybrid waves for magnetopause parameters indicates that the quasilinear state may never arise because coalescence to large wavelengths, followed by collapse once a critical wavelengths is reached, occur on a time scale faster than the quasilinear diffusion. In this case, an inhomogeneous boundary layer is to be expected. More simulations are required over longer time periods to explore whether this nonlinear evolution really takes place at the magnetopause.

Treumann, R. A.↗

Coherent Pump-Probe Interactions and Terahertz Intersubband Gain in Semiconductor Quantum Wells

In recent years there has been considerable interest in intersubband-transition-based infrared semiconductor quantum well (QW) lasers because of their potential applications. In the mid-infrared range, both electrically-injected quantum cascade lasers [1] and optically-pumped multiple QW lasers [2] have been experimentally realized. In these studies, optical gain is due to population inversion between the lasing subbands. It was also proposed that stimulated Raman scattering in QW systems can produce net infrared optical gain [3j. In such a nonlinear optical scheme, the appearance of optical gain that may lead to intersubband Raman lasers does not rely on the population inversion. Since, in tile resonant Raman process (Raman gain is the largest in this case), the pump field induces population redistribution among subbands in the QW s ystem, it seems that a realistic estimate of the optical gain has to include this effect. Perturbative calculations used in the previous work [3] may overestimate the Raman gain. In this paper we present a nonperturbative calculation of terahertz gain of optically-pumped semiconductor step quantum wells. Limiting optical transitions within the conduction band of QW, we solve the pump-field-induced nonequilibrium distribution function for each subband of the QW system from a set of coupled rate equations. Both intrasubband and intersubband relaxation processes in the quantum well system are included. Taking into account the coherent interactions between pump and THz (signal) waves, we we derive the susceptibility of the QW system for the THz field. For a GaAs/AlGaAs step QW, we calculate the Thz gain spectrum for different pump frequencies and intensities. Under moderately strong pumping (approximately 0.3 MW/sq cm), a significant THz gain (approximately 300/m) is predicted. It is also shown that the coherent wave interactions (resonant stimulated Raman processes) contribute significantly to the THz gain.

Liu, Ansheng↗

The Magnetospheric Constellation Mission. Dynamic Response and Coupling Observatory (DRACO): Understanding the Global Dynamics of the Structured Magnetotail

Magnetospheric Constellation Dynamic Response and Coupling Observatory (DRACO) is the Solar Terrestrial Probe (STP) designed to understand the nonlinear dynamics, responses, and connections within the Earth's structured magnetotail, using a constellation of approximately 50 to 100 distributed vector measurement spacecraft. DRACO will reveal magnetotail processes operating within a domain extending 20 Earth radii (R(sub E)) across the tail and 40 R(sub E)down the tail, on spatial and time scales accessible to global circulation models, i.e., approximately 2 R(sub E) and 10 seconds.

Source record↗

Parametric Study on the Response of Compression-Loaded Composite Shells With Geometric and Material Imperfections

The results of a parametric study of the effects of initial imperfections on the buckling and postbuckling response of three unstiffened thinwalled compression-loaded graphite-epoxy cylindrical shells with different orthotropic and quasi-isotropic shell-wall laminates are presented. The imperfections considered include initial geometric shell-wall midsurface imperfections, shell-wall thickness variations, local shell-wall ply-gaps associated with the fabrication process, shell-end geometric imperfections, nonuniform applied end loads, and variations in the boundary conditions including the effects of elastic boundary conditions. A high-fidelity nonlinear shell analysis procedure that accurately accounts for the effects of these imperfections on the nonlinear responses and buckling loads of the shells is described. The analysis procedure includes a nonlinear static analysis that predicts stable response characteristics of the shells and a nonlinear transient analysis that predicts unstable response characteristics.

Hilburger, Mark W.↗

The Magnetospheric Constellation Mission Dynamic Response and Coupling Observatory (DRACO): Understanding the Global Dynamics of the Structured Magnetotail

Magnetospheric Constellation Dynamic Response and Coupling Observatory (DRACO) is the Solar Terrestrial Probe (STP) designed to understand the nonlinear dynamics, responses, and connections within the Earth's structured magnetotail, using a constellation of approximately 50 to 100 distributed vector measurement spacecraft. DRACO will reveal magnetotail processes operating within a domain extending 20 Earth radii (R(sub E)) across the tail and 40 R(sub E)down the tail, on spatial and time scales accessible to global circulation models, i.e., approximately 2 R(sub E) and 10 seconds

Source record↗

Novel Estimation of Pilot Performance Characteristics

Two mechanisms internal to the pilot that affect performance during a tracking task are: 1) Pilot equalization (i.e. lead/lag); and 2) Pilot gain (i.e. sensitivity to the error signal). For some applications McRuer's Crossover Model can be used to anticipate what equalization will be employed to control a vehicle's dynamics. McRuer also established approximate time delays associated with different types of equalization - the more cognitive processing that is required due to equalization difficulty, the larger the time delay. However, the Crossover Model does not predict what the pilot gain will be. A nonlinear pilot control technique, observed and coined by the authors as 'amplitude clipping', is shown to improve stability, performance, and reduce workload when employed with vehicle dynamics that require high lead compensation by the pilot. Combining linear and nonlinear methods a novel approach is used to measure the pilot control parameters when amplitude clipping is present, allowing precise measurement in real time of key pilot control parameters. Based on the results of an experiment which was designed to probe workload primary drivers, a method is developed that estimates pilot spare capacity from readily observable measures and is tested for generality using multi-axis flight data. This paper documents the initial steps to developing a novel, simple objective metric for assessing pilot workload and its variation over time across a wide variety of tasks. Additionally, it offers a tangible, easily implementable methodology for anticipating a pilot's operating parameters and workload, and an effective design tool. The model shows promise in being able to precisely predict the actual pilot settings and workload, and observed tolerance of pilot parameter variation over the course of operation. Finally, an approach is proposed for generating Cooper-Harper ratings based on the workload and parameter estimation methodology.

Human pilot models↗

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗