Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High dimensional modeling”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗
NASA NTRS1993

Optimization for efficient structure-control systems

The efficiency of a structure-control system is a nondimensional parameter which indicates the fraction of the total control power expended usefully in controlling a finite-dimensional system. The balance of control power is wasted on the truncated dynamics serving no useful purpose towards the control objectives. Recently, it has been demonstrated that the concept of efficiency can be used to address a number of control issues encountered in the control of dynamic systems such as the spillover effects, selection of a good input configuration and obtaining reduced order control models. Reference (1) introduced the concept and presented analyses of several Linear Quadratic Regulator designs on the basis of their efficiencies. Encouraged by the results of Ref. (1), Ref. (2) introduces an efficiency modal analysis of a structure-control system which gives an internal characterization of the controller design and establishes the link between the control design and the initial disturbances to affect efficient structure-control system designs. The efficiency modal analysis leads to identification of principal controller directions (or controller modes) distinct from the structural natural modes. Thus ultimately, many issues of the structure-control system revolve around the idea of insuring compatibility of the structural modes and the controller modes with each other, the better the match the higher the efficiency. A key feature in controlling a reduced order model of a high dimensional (or infinity-dimensional distributed parameter system) structural dynamic system must be to achieve high efficiency of the control system while satisfying the control objectives and/or constraints. Formally, this can be achieved by designing the control system and structural parameters simultaneously within an optimization framework. The subject of this paper is to present such a design procedure.

Oz, Hayrani↗
DOE OSTI2025

Dynamical control in a prethermalized molecular ultracold plasma: Local dissipation drives global relaxation

Prethermalization occurs as an important phase in the dynamics of isolated many-body systems when coupled degrees of freedom relax in a subspace separated from a state of complete thermodynamic equilibrium by a gap in energy or other conserved quantity. Slow equilibration from a prethermal phase can localize the dynamics of natural and model systems despite high dimensionality and limited disorder. Here, we report the signature of an enduring prethermal regime of arrested relaxation in the molecular ultracold plasma that forms following the avalanche of a state-selected Rydberg gas of nitric oxide. For a wide range of initial conditions, this system enters a critical phase in which a density of NO + and electrons balances a population of Rydberg molecules. Electron collisions mix orbital angular momentum, scattering Rydberg electrons to states of very high ℓ . The rapid predissociation of molecules in states of low- ℓ purifies this high- ℓ ensemble. The angular momentum barrier that separates NO + ions and electrons creates an extraordinary gap between the plasma states of n ≈ ℓ , with measured n > 200 and penetrating states of ℓ = 0 , 1 , and 2. Evolution to a statistically equilibrated state of N and O atoms cannot occur without Rydberg electron penetration, and this gap blocks relaxation for a millisecond or more. A weak radiofrequency (RF) field drives ℓ -mixing electron collisions that erase this gap, causing wholesale dissipation. Remarkably, a local quantum-state transition promoting an exceedingly small fraction of the molecules in the prethermalized ensemble to a predissociative state also acts with a global effect, driving the entire system to a dissipative equilibrium. Using the Lindblad master equation, we illustrate qualitatively similar dynamics for a toy model of an open quantum system that consists of a localized set of spins on which dissipation acts at a single site as a gateway to equilibrium. Published by the American Physical Society 2025

Wang, Ruoxi↗
DOE OSTI2025

A review of low-rank methods for time-dependent kinetic simulations

Time-dependent kinetic models are ubiquitous in computational science and engineering. The underlying integro-differential equations in these models are high-dimensional, comprised of a six–dimensional phase space, making simulations of such phenomena extremely expensive. In this article we demonstrate that in many situations, the solution to kinetics problems lives on a low dimensional manifold that can be described by a low-rank matrix or tensor approximation. We then review the recent development of so-called low-rank methods that evolve the solution on this manifold. The two classes of methods we review are the dynamical low-rank (DLR) method, which derives differential equations for the low-rank factors, and a Step-and-Truncate (SAT) approach, which projects the solution onto the low-rank representation after each time step. Thorough discussions of time integrators, tensor decompositions, and method properties such as structure preservation and computational efficiency are included. We further show examples of low-rank methods as applied to particle transport and plasma dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Theoretical study of the high-latitude ionosphere's response to multicell convection patterns

A time-dependent three-dimensional model of the high-latitude ionosphere is used to study the characteristic ionospheric signatures associated with two-, three-, and four-cell plasma convection patterns. It is found that, for two-cell convection, the antisunward flow of plasma from the dayside into the polar cap acts to maintain the densities in this region in winter. For four-cell convection, the two additional convection cells in the polar cap are in darkness most of the time, and the resulting O(+) decay acts to produce twin polar holes that are separated by a sun-aligned ridge of enhanced ionization due to theta-auroral precipitation. For three-cell convection, only one polar hole forms in the total electron density, and an additional O(+) depletion region develops near noon. In this region there are strong electric fields, high ion temperatures, and an enhanced rate of O(+) - NO(+) conversion.

Sojka, J. J.↗
NASA NTRS1992

The neutral circulation in the vicinity of a stable auroral arc

The effects of asymmetric potential distributions, arc-neutral feedback, and cross-arc winds on the neutral circulation near auroral arc were examined using a sophisticated two-dimensional high-resolution model of the neutral atmosphere which was preconditioned by allowing weak diffuse aurora to exist prior to the formation of an arc. Simulations were performed with asymmetric potential structures and with an initial state which was warmed and accelerated due to a preexisting diffuse aurora, as well as with arc-zonal wind feedback and with large-scale cross-arc flow. It is shown that, with a preexisting zonal flow and asymmetrical arc structures, strong zonal winds are generated within the arc. Simulations with the arc-zonal feedback indicate that the effect of feedback is not large because of the limited duration of an arc. Simulations with a strong large-scale cross-arc wind indicate that the zonal wind is not strongly sensitive to the large-scale wind.

Walterscheid, R. L.↗
NASA NTRS2015

Assessment of Slat Noise Predictions for 30P30N High-Lift Configuration From BANC-III Workshop

This paper presents a summary of the computational predictions and measurement data contributed to Category 7 of the 3rd AIAA Workshop on Benchmark Problems for Airframe Noise Computations (BANC-III), which was held in Atlanta, GA, on June 14-15, 2014. Category 7 represents the first slat-noise configuration to be investigated under the BANC series of workshops, namely, the 30P30N two-dimensional high-lift model (with a slat contour that was slightly modified to enable unsteady pressure measurements) at an angle of attack that is relevant to approach conditions. Originally developed for a CFD challenge workshop to assess computational fluid dynamics techniques for steady high-lift predictions, the 30P30N configurations has provided a valuable opportunity for the airframe noise community to collectively assess and advance the computational and experimental techniques for slat noise. The contributed solutions are compared with each other as well as with the initial measurements that became available just prior to the BANC-III Workshop. Specific features of a number of computational solutions on the finer grids compare reasonably well with the initial measurements from FSU and JAXA facilities and/or with each other. However, no single solution (or a subset of solutions) could be identified as clearly superior to the remaining solutions. Grid sensitivity studies presented by multiple BANC-III participants demonstrated a relatively consistent trend of reduced surface pressure fluctuations, higher levels of turbulent kinetic energy in the flow, and lower levels of both narrow band peaks and the broadband component of unsteady pressure spectra in the nearfield and farfield. The lessons learned from the BANC-III contributions have been used to identify improvements to the problem statement for future Category-7 investigations.

Choudhari, Meelan↗
DOE OSTI2024

Accurate data-driven surrogates of dynamical systems for forward propagation of uncertainty

Stochastic collocation (SC) is a well-known non-intrusive method of constructing surrogate models for uncertainty quantification. In dynamical systems, SC is especially suited for full-field uncertainty propagation that characterizes the distributions of the high-dimensional solution fields of a model with stochastic input parameters. However, due to the highly nonlinear nature of the parameter-to-solution map in even the simplest dynamical systems, the constructed SC surrogates are often inaccurate. Here, this work presents an alternative approach, where we apply the SC approximation over the dynamics of the model, rather than the solution. By combining the data-driven sparse identification of nonlinear dynamics framework with SC, we construct dynamics surrogates and integrate them through time to construct the surrogate solutions. We demonstrate that the SC-over-dynamics framework leads to smaller errors, both in terms of the approximated system trajectories as well as the model state distributions, when compared against full-field SC applied to the solutions directly. We present numerical evidence of this improvement using three test problems: a chaotic ordinary differential equation, and two partial differential equations from solid mechanics.

42 ENGINEERING↗
NASA NTRS2005

The Importance of Detailed Component Simulations in the Feedsystem Development for a Two-Stage-to Orbit Reusable Launch Vehicle

To meet the requirements for the 2nd Generation Reusable Launch Vehicle (RLV), a unique propulsion feed system concept was identified using crossfeed between the booster and orbiter stages that could reduce the Two-Stage-to-Orbit (TSTO) vehicle weight and development cost by approximately 25%. A Main Propulsion System (MPS) crossfeed water demonstration test program was configured to address all the activities required to reduce the risks for the MPS crossfeed system. A transient, one-dimensional system simulation was developed for the subscale crossfeed water flow tests. To ensure accurate representation of the crossfeed valve's dynamics in the system model, a high-fidelity, three-dimensional, computational fluid-dynamics (CFD) model was employed. The results from the CFD model were used to specify the valve's flow characteristics in the system simulation. This yielded a crossfeed system model that was anchored to the specific valve hardware and achieved good agreement with the measured test data. These results allowed the transient models to be correlated and validated and used for full scale mission predictions. The full scale model simulations indicate crossfeed is ' viable with the system pressure disturbances at the crossfeed transition being less than experienced by the propulsion system during engine start and shutdown transients.

Mazurkivich, Pete↗
NASA NTRS1990

A two-dimensional kinematic dynamo model of the ionospheric magnetic field at Venus

The results of a high-resolution, two-dimensional, time dependent, kinematic dynamo model of the ionospheric magnetic field of Venus are presented. Various one-dimensional models are considered and the two-dimensional model is then detailed. In this model, the two-dimensional magnetic induction equation, the magnetic diffusion-convection equation, is numerically solved using specified plasma velocities. Origins of the vertical velocity profile and of the horizontal velocities are discussed. It is argued that the basic features of the vertical magnetic field profile remain unaltered by horizontal flow effects and also that horizontal plasma flow can strongly affect the magnetic field for altitudes above 300 km.

Cravens, T. E.↗
NASA NTRS1991

Effects of engine emissions from high-speed civil transport aircraft: A two-dimensional modeling study, part 1

The AER two-dimensional chemistry-transport model is used to study the effect on stratospheric ozone (O3) from operations of supersonic and subsonic aircraft. The study is based on six emission scenarios provided to AER. The study showed that: (1) the O3 response is dominated by the portion of the emitted nitrogen compounds that is entrained in the stratosphere; (2) the entrainment is a sensitive function of the altitude at which the material is injected; (3) the O3 removal efficiency of the emitted material depends on the concentrations of trace gases in the background atmosphere; and (4) evaluation of the impact of fleet operations in the future atmosphere must take into account the expected changes in trace gas concentrations from other activities. Areas for model improvements in future studies are also discussed.

Ko, Malcolm K. W.↗
NASA NTRS1991

Effects of engine emissions from high-speed civil transport aircraft: A two-dimensional modeling study, part 2

The AER two-dimensional chemistry-transport model is used to study the effect of supersonic and subsonic aircraft operation in the 2010 atmosphere on stratospheric ozone (O3). The results show that: (1) the calculated O3 response is smaller in the 2010 atmosphere compared to previous calculations performed in the 1980 atmosphere; (2) with the emissions provided, the calculated decrease in O3 column is less than 1 percent; and (3) the effect of model grid resolution on O3 response is small provided that the physics is not modified.

Ko, Malcolm K. W.↗
NASA NTRS1992

Two-equation turbulence modeling for 3-D hypersonic flows

An investigation to verify, incorporate and develop two-equation turbulence models for three-dimensional high speed flows is presented. The current design effort of hypersonic vehicles has led to an intensive study of turbulence models for compressible hypersonic flows. This research complements an extensive review of experimental data and the current development of 2D turbulence models. The review of experimental data on 2D and 3D flows includes complex hypersonic flows with pressure profiles, skin friction, wall heat transfer, and turbulence statistics data. In a parallel effort, turbulence models for high speed flows have been tested against flat plate boundary layers, and are being tested against the 2D database. In the present paper, we present the results of 3D Navier-Stokes numerical simulations with an improved k-omega two-equation turbulence model against experimental data and empirical correlations of an adiabatic flat plate boundary layer, a cold wall flat plate boundary layer, and a 3D database flow, the interaction of an oblique shock wave and a thick turbulent boundary layer with a free stream Mach number = 8.18 and Reynolds number = 5 x 10 to the 6th.

Bardina, J. E.↗
NASA NTRS1996

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗
DOE OSTI2026

A procedure for rule extraction from a Self-Organising plasma disruption predictor for JET

In a previous paper, a Self-Organizing Map had proven to be able to identify the regions of the plasma operative space characterizing the pre-disruptive phase at JET without relying on any a priori information. One of the strengths of this disruption predictor lies in its inherent self-organization capability. The Self-Organizing Map discovers non-trivial relationships and captures the complicated interplay of device diagnostics on the internal plasma states directly from the experimental data. Moreover, the provided model allows the visualization of high-dimensional plasma parameters and facilitates easy interrogation of the model to understand the reasons behind its correlations. In this paper, an additional step is taken towards the interpretability of models for predicting disruptions by training a Decision Tree to classify the plasma states according to the interpretation provided by the Self-Organizing Map (stable or at high risk of disruptions). The Decision tree provides a set of rules which describe the transition of the plasma towards the pre-disruptive phase as visualized in the Self-Organizing Map. The obtained rules for the database explored in the study identify four regions in the map, two of which are at risk of disruption. These regions correspond to partitions of a 3D space based on the peaking factors of the core and divertor radiation, as well as the Locked Mode. The agreement between the Self-Organizing Map answers and the rules supplied by the Decision Tree is confirmed by the comparison of the performance exhibited by the two models in the prediction of disruptions.

Setzu, Samuele [Univ. of Cagliari, Monserrato, Cag↗
NASA NTRS1982

The baroclinic instability of highly structured one-dimensional basic states

A frictionless, adiabatic, quasi-geostrophic model on a beta-plane was employed to examine the baroclinic instability of highly structured one-dimensional basic states. Definitions of the model eigenvalue equation, the boundary conditions, and the fluxes and energetics are provided, and highly structured average wind speed and temperature profiles were selected to represent a wide range of latitudes. Attention was given to the instability source, wavelength, growth rate, phase velocity, steering levels, and the vertical structure of their amplitude, phase, meridional entropy transport, and conversion of potential-to-kinetic energy. Green modes, Charney-Eady modes, and the modes defined in the paper are examined in connection with the highly structured profiles. The eigenvalue approach is shown to be superior to the numerical time integration approach and necessary conditions for instabilities are defined.

Fullmer, J. W. A.↗
DOE OSTI2025

Persistent Sampling: Enhancing the Efficiency of Sequential Monte Carlo

Sequential Monte Carlo (SMC) samplers are powerful tools for Bayesian inference but suffer from high computational costs due to their reliance on large particle ensembles for accurate estimates. We introduce persistent sampling (PS), an extension of SMC that systematically retains and reuses particles from all prior iterations to construct a growing, weighted ensemble. By leveraging multiple importance sampling and resampling from a mixture of historical distributions, PS mitigates the need for excessively large particle counts, directly addressing key limitations of SMC such as particle impoverishment and mode collapse. Crucially, PS achieves this without additional likelihood evaluations-weights for persistent particles are computed using cached likelihood values. This framework not only yields more accurate posterior approximations but also produces marginal likelihood estimates with significantly lower variance, enhancing reliability in model comparison. Furthermore, the persistent ensemble enables efficient adaptation of transition kernels by leveraging a larger, decorrelated particle pool. Experiments on high-dimensional Gaussian mixtures, hierarchical models, and non-convex targets demonstrate that PS consistently outperforms standard SMC and related variants, including recycled and waste-free SMC, achieving substantial reductions in mean squared error for posterior expectations and evidence estimates, all at reduced computational cost. PS thus establishes itself as a robust, scalable, and efficient alternative for complex Bayesian inference tasks.

Karamanis, Minas↗