Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic model reduction”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Comparing Delay-, Distance-, and Cordon-Based Congestion Pricing Strategies Via Large-Scale Simulation

This study compares the impacts of delay-, distance-, and cordon-based congestion pricing strategies for Austin, Texas, using the POLARIS agent-based activity-based travel demand simulation model. This approach enables agent-level heterogeneity and realistic choice options (including destination, mode, and activity scheduling) for dynamic traffic assignment and congestion feedbacks across a major metro region, which are features lacking in past work. To ensure comparability, distance-based tolls were set to generate the same revenue as delay-based tolling of $3.5 M/day, averaging $1.17/resident/day or $0.42/vehicle-trip. Delay-based pricing delivers 44% lower network delay and 13% lower VHT compared to the no-toll baseline, levels unmatched by other pricing strategies. At the height of the AM peak, drivers pay up to $0.13/mile on average, though most links in the network remain untolled. Distance-based pricing is the most effective at reducing VMT (by 4%), but VHT reductions (of 6%) primarily stem from drivers selecting closer destinations, achieving only one-fourth the delay reduction of delay-based pricing. Across various implementations of delay- and distance-based pricing, the results suggest that spatial variations of tolls are far more important than temporal variations. Cordon tolls produce minimal impacts at the network-wide level, but offer substantial delay reductions inside the cordon. Other major findings include: 1) delay-based pricing increases trip-making during the PM peak period due to backward shifts in discretionary-activity start times by higher-income residents; and 2) tolls’ spatial impacts, including changes in network flows and tolls paid by residents, vary substantially between delay- and distance-based pricing strategies.

Agent-based modeling↗
DOE OSTI2021

Townsend to glow discharge transition for a nanosecond pulse plasma in helium: space charge formation and resulting electric field dynamics

Stark polarization spectroscopy is used to investigate the temporal evolution of the electric field distribution in the cathode region of a nanosecond pulsed discharge in helium at 120 Torr. The measurements are performed on the He I transition at 492.19 nm, during the early stages of the discharge formation. The experimental results are compared with the predictions of a 1D fluid model. Time-resolved ICCD images show that the discharge develops as a diffuse, cathode-directed ionization wave with a Townsend-like feature before transitioning into a glow-like structure. Near anode instabilities characterized by filament formation were observed near the high voltage electrode. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. Additionally, the cathode electric field corresponding to this sheath with a thickness of 250 μm is about 40 kV cm -1 . A subsequent steep increase of the discharge current leads to a further reduction of the sheath width. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solving differential‐algebraic equations in power system dynamic analysis with quantum computing

Abstract Power system dynamics are generally modeled by high dimensional non‐linear differential‐algebraic equations (DAEs) given a large number of components forming the network. These DAEs' complexity can grow exponentially due to the increasing penetration of distributed energy resources, whereas their computation time becomes sensitive due to the increasing interconnection of the power grid with other energy systems. This paper demonstrates the use of quantum computing algorithms to solve DAEs for power system dynamic analysis. We leverage a symbolic programming framework to equivalently convert the power system's DAEs into ordinary differential equations (ODEs) using index reduction methods and then encode their data into qubits using amplitude encoding. The system non‐linearity is captured by Hamiltonian simulation with truncated Taylor expansion so that state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can solve the power system's DAEs accurately with a computational complexity polynomial in the logarithm of the system dimension. We also illustrate the use of recent advanced tools in scientific machine learning for implementing complex computing concepts, that is, Taylor expansion, DAEs/ODEs transformation, and quantum computing solver with abstract representation for power engineering applications.

computational complexity↗
DOE OSTI2025

Symmetry instability induced by topological phase transitions

The symmetry-topology interplay dictates how to define order parameters and classify material ordered phases. However, current understanding of this interplay has been predominately approached from a one-sided perspective, with topological states being classified within the constraints imposed by specific fixed symmetries. Here, in this study, we complete this full circle by demonstrating spontaneous symmetry breaking that results from a periodic alteration of topological phases induced by light in a centrosymmetric Dirac material ZrTe 5 . The distinguishing feature is the observation of robust correlation and striking anomalies in the fluence and temperature dependence of key transport parameters. First, both shift current 𝐽 𝑠 and displacement current 𝐽 𝑑 , arising from interband transition and infrared phonon driving, respectively, along with charge carrier pumping, exhibit similar behaviors. Second, they all peak at similar low pump fluence, followed by a subsequent reduction as the fluence further increases. This behavior cannot be explained by conventional energetically allowed, direct excitations. Third, all the three observables exhibit anomalies when they approach the topological phase transition temperature. These results highlight the unique low-energy pumping behaviors in ZrTe 5 , characterized by reversible fluence dependence and a “hinge-like” interaction that connects various electronic and lattice observables, including phonons, charge carriers, and currents. Our findings, supported by model analysis, provide key insights into the fragility of crystalline (inversion) and time-reversal symmetries during the dynamics of topological phase transitions. This fragility drives spontaneous symmetry breaking, evidenced by the synchronized emergence of off-resonant infrared phonons and broken-symmetry photocurrents.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2014

Aeroelastic Tailoring of the NASA Common Research Model via Novel Material and Structural Configurations

This work explores the use of tow steered composite laminates, functionally graded metals (FGM), thickness distributions, and curvilinear rib/spar/stringer topologies for aeroelastic tailoring. Parameterized models of the Common Research Model (CRM) wing box have been developed for passive aeroelastic tailoring trade studies. Metrics of interest include the wing weight, the onset of dynamic flutter, and the static aeroelastic stresses. Compared to a baseline structure, the lowest aggregate static wing stresses could be obtained with tow steered skins (47% improvement), and many of these designs could reduce weight as well (up to 14%). For these structures, the trade-off between flutter speed and weight is generally strong, although one case showed both a 100% flutter improvement and a 3.5% weight reduction. Material grading showed no benefit in the skins, but moderate flutter speed improvements (with no weight or stress increase) could be obtained by grading the spars (4.8%) or ribs (3.2%), where the best flutter results were obtained by grading both thickness and material. For the topology work, large weight reductions were obtained by removing an inner spar, and performance was maintained by shifting stringers forward and/or using curvilinear ribs: 5.6% weight reduction, a 13.9% improvement in flutter speed, but a 3.0% increase in stress levels. Flutter resistance was also maintained using straightrotated ribs although the design had a 4.2% lower flutter speed than the curved ribs of similar weight and stress levels were higher. These results will guide the development of a future design optimization scheme established to exploit and combine the individual attributes of these technologies.

Jutte, Christine V.↗
NASA NTRS2016

Elevated Aerosol Layers and Their Radiative Impact over Kanpur During Monsoon Onset Period

Accurate information about aerosol vertical distribution is needed to reduce uncertainties in aerosol radiative forcing and its effect on atmospheric dynamics. The present study deals with synergistic analyses of aerosol vertical distribution and aerosol optical depth (AOD) with meteorological variables using multisatellite and ground-based remote sensors over Kanpur in central Indo-Gangetic Plain (IGP). Micro-Pulse Lidar Network-derived aerosol vertical extinction (sigma) profiles are analyzed to quantify the interannual and daytime variations during monsoon onset period (May-June) for 2009-2011. The mean aerosol profile is broadly categorized into two layers viz., a surface layer (SL) extending up to 1.5 km (where sigma decreased exponentially with height) and an elevated aerosol layer (EAL) extending between 1.5 and 5.5 km. The increase in total columnar aerosol loading is associated with relatively higher increase in contribution from EAL loading than that from SL. The mean contributions of EALs are about 60%, 51%, and 50% to total columnar AOD during 2009, 2010, and 2011, respectively. We observe distinct parabolic EALs during early morning and late evening but uniformly mixed EALs during midday. The interannual and daytime variations of EALs are mainly influenced by long-range transport and convective capacity of the local emissions, respectively. Radiative flux analysis shows that clear-sky incoming solar radiation at surface is reduced with increase in AOD, which indicates significant cooling at surface. Collocated analysis of atmospheric temperature and aerosol loading reveals that increase in AOD not only resulted in surface dimming but also reduced the temperature (approximately 2-3 C) of lower troposphere (below 3 km altitude). Radiative transfer simulations indicate that the reduction of incoming solar radiation at surface is mainly due to increased absorption by EALs (with increase in total AOD). The observed cooling in lower troposphere in high aerosol loading scenario could be understood as a dynamical feedback of EAL-induced stratification of lower troposphere. Further, the observed radiative effect of EALs increases the stability of the lower troposphere, which could modulate the large-scale atmospheric dynamics during monsoon onset period. These findings encourage follow-up studies on the implication of EALs to the Indian summer monsoon dynamics using numerical models.

aerosol vertical distribution↗
NASA NTRS2012

Integration of Linear Dynamic Emission and Climate Models with Air Traffic Simulations

Future air traffic management systems are required to balance the conflicting objectives of maximizing safety and efficiency of traffic flows while minimizing the climate impact of aviation emissions and contrails. Integrating emission and climate models together with air traffic simulations improve the understanding of the complex interaction between the physical climate system, carbon and other greenhouse gas emissions and aviation activity. This paper integrates a national-level air traffic simulation and optimization capability with simple climate models and carbon cycle models, and climate metrics to assess the impact of aviation on climate. The capability can be used to make trade-offs between extra fuel cost and reduction in global surface temperature change. The parameters in the simulation can be used to evaluate the effect of various uncertainties in emission models and contrails and the impact of different decision horizons. Alternatively, the optimization results from the simulation can be used as inputs to other tools that monetize global climate impacts like the FAA s Aviation Environmental Portfolio Management Tool for Impacts.

Sridhar, Banavar↗
NASA NTRS1999

Fractal analysis of heart rate dynamics as a predictor of mortality in patients with depressed left ventricular function after acute myocardial infarction. TRACE Investigators. TRAndolapril Cardiac Evaluation

A number of new methods have been recently developed to quantify complex heart rate (HR) dynamics based on nonlinear and fractal analysis, but their value in risk stratification has not been evaluated. This study was designed to determine whether selected new dynamic analysis methods of HR variability predict mortality in patients with depressed left ventricular (LV) function after acute myocardial infarction (AMI). Traditional time- and frequency-domain HR variability indexes along with short-term fractal-like correlation properties of RR intervals (exponent alpha) and power-law scaling (exponent beta) were studied in 159 patients with depressed LV function (ejection fraction <35%) after an AMI. By the end of 4-year follow-up, 72 patients (45%) had died and 87 (55%) were still alive. Short-term scaling exponent alpha (1.07 +/- 0.26 vs 0.90 +/- 0.26, p <0.001) and power-law slope beta (-1.35 +/- 0.23 vs -1.44 +/- 0.25, p <0.05) differed between survivors and those who died, but none of the traditional HR variability measures differed between these groups. Among all analyzed variables, reduced scaling exponent alpha (<0.85) was the best univariable predictor of mortality (relative risk 3.17, 95% confidence interval 1.96 to 5.15, p <0.0001), with positive and negative predictive accuracies of 65% and 86%, respectively. In the multivariable Cox proportional hazards analysis, mortality was independently predicted by the reduced exponent alpha (p <0.001) after adjustment for several clinical variables and LV function. A short-term fractal-like scaling exponent was the most powerful HR variability index in predicting mortality in patients with depressed LV function. Reduction in fractal correlation properties implies more random short-term HR dynamics in patients with increased risk of death after AMI.

Non-NASA Center↗
DOE OSTI2025

Oscillate and renormalize: Fast phonons reshape the Kondo effect in flat-band systems

Here, we examine the interplay between electron correlations and phonons in an Anderson-Holstein impurity model with an Einstein phonon. When the phonons are slow compared to charge fluctuations (frequency 𝜔 0 ≪ 𝑈/2, the onsite Coulomb scale 𝑈/2), we demonstrate analytically that the expected phonon-mediated reduction of interactions is completely suppressed, even in the strong-coupling regime. This suppression arises from the oscillator's inability to respond to rapid charge fluctuations, manifested as a compensation effect between the polaronic cloud and the excited-state phonons associated with valence fluctuations. We identify a frozen mixed valence phase, above a threshold dimensionless electron-phonon coupling 𝛼* when the phonons are slow, where the static phonon cloud locks the impurity into specific valence configurations, potentially explaining the puzzling coexistence of mixed valence behavior and insulating properties in materials like rust. Conversely, when the phonon is fast (𝜔 0 ≳ 𝑈/2), the system exhibits conventional polaronic behavior with renormalized onsite interactions effectively 𝑈 eff due to phonon-mediated attraction, with additional satellite features in the local spectral function due to phonon excitations. Using numerical renormalization group calculations, a fully dynamic renormalization technique, we confirm these behaviors in both regimes. These findings have important implications for strongly correlated systems where phonon energy scales may be comparable to the Coulomb scale, such as in twisted bilayer graphene, necessitating careful consideration of interaction renormalizations in theoretical models.

Anderson impurity model↗
NASA NTRS1977

Electron dynamics in a plasma focus

Results are presented of a numerical integration of the three-dimensional relativistic equations of motion of electrons subject to given electric and magnetic fields deduced from experiments. Fields due to two different models are investigated. For the first model, the fields are those due to a circular distribution of axial current filaments. As the current filaments collapse toward the axis, large azimuthal magnetic and axial electric fields are induced. These fields effectively heat the electrons to a temperature of approximately 8 keV and accelerate electrons within the radius of the filaments to high axial velocities. Similar results are obtained for the current-reduction phase of focus formation. For the second model, the fields are those due to a uniform current distribution. Both the current-reduction and the compression phases were studied. These is little heating or acceleration of electrons during the compression phase because the electrons are tied to the magnetic field. However, during the current-reduction phase, electrons near the axis are accelerated toward the center electrode and reach energies of 100 keV. A criterion is obtained which limits the runaway electron current to about 400 A.

Hohl, F.↗
DOE OSTI2023

Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN

Model simulations using the E3SM Land Model (ELM) coupled to the PFLOTRAN reactive transport model via the Alquimia interface. The simulations were conducted for a tidal salt marsh at the Plum Island Ecosystems LTER near Rowley, Massachusetts, USA. Model simulations were forced using site-specific tidal cycles and salinity, and the simulations used a biogeochemical reaction network including aerobic decomposition, sulfate reduction, iron reduction, and methanogenesis. Model outputs include simulated carbon stocks, carbon dioxide and methane fluxes, and porewater concentrations of key solutes related to sulfur, iron, and carbon cycling. The model simulations included a saline simulation (with tidal sulfate inputs), a fresh simulation (with low salinity and low sulfate inputs), and a saline simulation with lower vegetation productivity to represent the effect of salinity on vegetation. These simulations were conducted to demonstrate that a new model framework incorporating subsurface redox and biogeochemical interactions into a land surface model could reproduce measured surface greenhouse gas fluxes and biogeochemical dynamics in tidal marsh ecosystems, and to test whether including redox interactions in a land surface model would allow the model to resolve contrasts in biogeochemical cycling and greenhouse gas production between saline and freshwater wetlands.The data package includes gzipped tar archives (which can be expanded using standard tar and gzip utilities) of model outputs from three model configurations: saline subsurface and reduced vegetation productivity related to salinity; saline subsurface with vegetation productivity not reduced; and freshwater. Also included are code for the modified E3SM model, Alquimia interface, and PFLOTRAN reactive transport simulator in gzipped tar format; plain text parameter and configuration files; python code files for visualizing model output and defining model configurations; and model output, tide and salinity forcing, and configuration files in netCDF format. See the README.md file in the data package for a detailed description of all files contained in the package. All files are in netCDF (.nc), gzipped tar archive (.tar.gz or .tgz), or text (all other files).Updated: May 13, 2024. Model output, E3SM code, PFLOTRAN input files, and python codes for visualizing results were updated to reflect changes made for the manuscript revision. The updated archive reflects the code and model output from the final accepted manuscript. Changes included updated reaction parameters reflecting improved parameterization and additional comparisons with field measurements. E3SM code changes included better support for multiple grid cells and improved flow and transport parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Wind Resource Data for Southeast Asia Using a Hybrid Numerical Weather Prediction with Machine Learning Super Resolution Approach

In this work we develop and present a machine learning based downscaling approach using generative adversarial networks (GANs). GANs learn to distinguish the relationships between low-resolution and high-resolution simulations and generate accurate high-resolution output from low-resolution input (Stengel, Glaws, Hettinger, & King, 2020). Low-resolution numerical weather prediction (NWP) simulations at 9-km spatial and 60-minute temporal resolution were executed over Southeast Asia to provide input to the GANs model. GANs for wind, temperature, and pressure were trained on coarsened WIND Toolkit data with a diverse sampling of terrain and meteorological conditions. After training, the NWP simulations over Southeast Asia were enhanced by 3x along each horizontal spatial dimension and 4x along the temporal dimension. This novel downscaling approach generated 15-year high-resolution wind, temperature, and pressure data from January 2007 through December 2021 at multiple hub heights over Southeast Asia at 3-km spatial and 15-minute temporal resolution with a 16x reduction in compute time over standard dynamical downscaling.

17 WIND ENERGY↗
NASA NTRS1991

A hybrid approach to test-analysis-model development for large space structures

The present FEM reduction method for the generation of test-analysis-models (TAMs) in test-analysis correlation addresses contentions that the current modal TAM is hypersensitive to differences between test model shapes and analysis mode shapes, thereby generating large off-diagonal terms within the orthogonality and cross-orthogonality matrices employed in test-analysis mode-shape correlation. A hybrid TAM methodology is accordingly developed which combines the exact representation of the FEM target modes from the modal TAM with the more accurate TAM representation of the residual modes; the superior residual dynamics representation of the hybrid TAM is demonstrated for a detailed representation of a large space structure.

Kammer, D. C.↗
NASA NTRS1992

Rolling maneuver load alleviation using active controls

Rolling Maneuver Load Alleviation (RMLA) was demonstrated on the Active Flexible Wing (AFW) wind tunnel model in the LaRC Transonic Dynamics Tunnel. The design objective was to develop a systematic approach for developing active control laws to alleviate wing incremental loads during roll maneuvers. Using linear load models for the AFW wind-tunnel model which were based on experimental measurements, two RMLA control laws were developed based on a single-degree-of-freedom roll model. The RMLA control laws utilized actuation of outboard control surface pairs to counteract incremental loads generated during rolling maneuvers and roll performance. To evaluate the RMLA control laws, roll maneuvers were performed in the wind tunnel at dynamic pressures of 150, 200, and 250 psf and Mach numbers of .33, .38, and .44, respectively. Loads obtained during these maneuvers were compared to baseline maneuver loads. For both RMLA controllers, the incremental torsion moments were reduced by up to 60 percent at all dynamic pressures and performance times. Results for bending moment load reductions during roll maneuvers varied. In addition, in a multiple function test, RMLA and flutter suppression system control laws were operated simultaneously during roll maneuvers at dynamic pressures 11 percent above the open-loop flutter dynamic pressure.

Woods-Vedeler, Jessica A.↗
NASA NTRSDate not supplied

Development of Data Reduction Methods for Dynamic Stability Tests in a Magnetic Suspension and Balance System

The NASA/ODU 6-inch Magnetic Suspension and Balance System (MSBS) is operational in a “free-to-yaw” configuration in a subsonic wind tunnel. The MSBS hardware has been heavily modified over recent years to utilize a transverse magnetization scheme, with mixed electromagnetic and motion tracking camera position sensors and a relatively sophisticated control system. The project aims to provide a unique test capability for dynamic stability testing of planetary entry capsules, free from support interference. Aerodynamic tests are being conducted, with a focus of validation of the data reduction methods used to extract dynamic stability parameters. The unique complication is the existence of magnetic, aerodynamic, and mass centers, where magnetic, aerodynamic, and inertia forces respectively are resolved. The paper will report on tests of a family of EEV-like geometries, where the conical half-angle varied between 45 and 60 degrees in steps of 3.75 degrees. The paper will present experimental results from the parametric family of EEV-like geometries. Simulation results will also be presented, based on a two or three degree-of-freedom Simulink model. The model allows the exploration of the effects of non-coincident centers and the sensitivity of results arising from uncertainties in location.

Entry Systems Modeling↗
DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗
NASA NTRS1977

Wind tunnel test and calibration techniques for the measurement of damping and dynamic cross derivatives due to pitching and yawing

The paper shows how the oscillatory apparatus described previously by Orlik-Rueckemann et al. (1974) for determining cross and cross-coupling moment derivatives due to pitching and yawing oscillation is strengthened, extensively modified, and adapted to permit operation in a large continuous-flow wind tunnel. Particular attention is given to an electromagnetic three-degrees-of-freedom calibrating system developed to verify the validity of the experimental method and the various aspects of data reduction procedure. Typical results from calibration tests and from a series of wind tunnel experiments are presented. The only significant limitation is that the apparatus can only be used for models that are relatively small and have a conically-shaped longitudinal cavity.

Hanff, E. S.↗
NASA NTRS1991

Numerical and literal aeroelastic-vehicle-model reduction for feedback control synthesis

The simplification of a high-order, literal model for large flexible aircraft is discussed. Areas of model fidelity that are critical if the model is to be used for control law synthesis are presented. Several simplification techniques, some new and some widely available, that can deliver the necessary model fidelity are presented and applied to a model from the literature. The techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory, is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle vibrational characteristics and rigid-body and aeroelastic stability derivatives, thus providing insight into the underlying causes for critical dynamic characteristics.

Newman, Brett↗