Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “criticality simulations”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Project Description for Student Symposium 2022: Modeling Approach to Critical in the Upcoming CERBERUS Experiment

The Critical Experiment Reflected By copper to BettEr Understand Scattering [CERBERUS] is an integral experiment that will be executed at the end of FY23. The experiment is designed to study the intermediate energy range, in particular elastic and inelastic scattering of neutrons in copper. Approaching criticality safely is of utmost importance, and having a computer mode is a useful tool when determining changes that can safely be made to the experiment. Multiple simulations are needed to analyze the approach to critical, and manually creating MCNP input files is both tedious and error prone. To solve this problem, a python script was created to develop and run these input files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Global micro-tearing modes in the wide pedestal of an NSTX plasma

Global linear electromagnetic gyrokinetic simulations of micro-tearing modes (MTMs) growing in the wide pedestal of the ELM-free NSTX discharge #132 588 are reported. Two micro-tearing branches growing near the surfaces of safety factor q = 4 (pedestal top) and q=5–6 (pedestal) are simulated. The frequencies of these MTMs with low toroidal mode number (n ranging from 3 to 6) are in agreement with the ones of quasi-coherent fluctuations present on the magnetic spectrogram. No kinetic ballooning mode (KBM) is observed at experimental values of β; instead, the experimental plasma lies 20% below the critical global KBM threshold. This confirms local simulation results stating that experimental NSTX profiles lie in the KBM stability region. Carbon impurities, which are abundant in this plasma, were found to play a significant role. Neglecting them leads to a ∼50% error in the growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Critical Scaling of Solid Fragmentation at Quasistatic and Finite Strain Rates

Here, using two-dimensional simulations of sheared, brittle solids, we characterize the resulting fragmentation and explore its underlying critical nature. Under quasistatic loading, a power-law distribution of fragment masses emerges after fracture which grows with increasing strain. With increasing strain rate, the maximum size of a grain decreases and a shallower distribution is produced. We propose a scaling theory for distributions based on a fractal scaling of the largest mass with system size in the quasistatic limit or with a correlation length that diverges as a power of rate in the finite-rate limit. Critical exponents are measured using finite-size scaling techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oscillatory Spreading and Inertia in Power Grids

The increase in variable renewable generators (VRGs) in power systems has altered the dynamics from a historical experience. VRGs introduce new sources of power oscillations, and the stabilizing response provided by synchronous generators (SGs, e.g., natural gas, coal, etc.), which help avoid some power fluctuations, will lessen as VRGs replace SGs. These changes have led to the need for new methods and metrics to quickly assess the likely oscillatory behavior for a particular network without performing computationally expensive simulations. This work studies the impact of a critical dynamical parameter - the inertia value - on the rest of a power system's oscillatory response to representative VRG perturbations. We use a known localization metric in a novel way to quantify the number of nodes responding to a perturbation and the magnitude of those responses. This metric allows us to relate the spread and severity of a system's power oscillations with inertia. We find that as inertia increases, the system response to node perturbations transitions from localized (only a few close nodes respond) to delocalized (many nodes across the network respond). We introduce a heuristic computed from the network Laplacian to relate this oscillatory transition to the network structure. We show that our heuristic accurately describes the spread of oscillations for a realistic power-system test case. Using a heuristic to determine the likely oscillatory behavior of a system given a set of parameters has wide applicability in power systems, and it could decrease the computational workload of planning and operation.

dynamical systems↗
NASA NTRS2009

Ultrasonic Phased Array Simulations of Welded Components at NASA

Comprehensive and accurate inspections of welded components have become of increasing importance as NASA develops new hardware such as Ares rocket segments for future exploration missions. Simulation and modeling will play an increasing role in the future for nondestructive evaluation in order to better understand the physics of the inspection process, to prove or disprove the feasibility for an inspection method or inspection scenario, for inspection optimization, for better understanding of experimental results, and for assessment of probability of detection. This study presents simulation and experimental results for an ultrasonic phased array inspection of a critical welded structure important for NASA future exploration vehicles. Keywords: nondestructive evaluation, computational simulation, ultrasonics, weld, modeling, phased array

Roth, D. J.↗
DOE OSTI2025

Acoustic Tomography of the Atmosphere: A Large-Eddy Simulation Sensitivity Study

Accurate measurement of atmospheric turbulent fluctuations is critical for understanding environmental dynamics and improving models in applications such as wind energy. Advanced remote sensing technologies are essential for capturing instantaneous velocity and temperature fluctuations. Acoustic tomography (AT) offers a promising approach that utilizes sound travel times between an array of transducers to reconstruct turbulence fields. This study presents a systematic evaluation of the time-dependent stochastic inversion (TDSI) algorithm for AT using synthetic travel-time measurements derived from large-eddy simulation (LES) fields under both neutral and convective atmospheric boundary-layer conditions. Unlike prior work that relied on field observations or idealized fields, the LES framework provides a ground-truth atmospheric state, enabling quantitative assessment of TDSI retrieval reliability, sensitivity to travel-time measurement noise, and dependence on covariance model parameters and temporal data integration. A detailed sensitivity analysis was conducted to determine the best-fit model parameters, identify the tolerance thresholds for parameter mismatch, and establish a maximum spatial resolution. The TDSI algorithm successfully reconstructed large-scale velocity and temperature fluctuations with root mean square errors ( RMSE s) below 0.35 m/s and 0.12 K, respectively. Spectral analysis established a maximum spatial resolution of approximately 1.4 m, and reconstructions remained robust for travel-time measurement uncertainties up to 0.002 s. These findings provide critical insights into the operational limits of TDSI and inform future applications of AT for atmospheric turbulence characterization and system design.

17 WIND ENERGY↗
DOE OSTI2021

Insights into Spontaneous Solid Electrolyte Interphase Formation at Magnesium Metal Anode Surface from Ab Initio Molecular Dynamics Simulations

Spontaneous chemical reactivity at multivalent (Mg, Ca, Zn, Al) electrode surfaces is critical to solid electrolyte interphase (SEI) formation, and hence, directly affects the longevity of batteries. Here, we report an investigation of the reactivity of 0.5 M Mg(TFSI)2 in 1,2-dimethoxyethane (DME) solvent at a Mg(0001) surface using ab initio molecular dynamics (AIMD) simulations and detailed Bader charge analysis. Based on the simulations, the initial degradation reactions of the electrolyte strongly depend on the structure of the Mg(TFSI)2 species near the anode surface. At the surface, the dissociation of Mg(TFSI)2 species occurs via cleavage of the N-S bond for the solvent separated ion pair (SSIP) and via cleavage of the C-S bond for the contact ion pair (CIP) configuration. In the case of the CIP, both TFSI anions undergo spontaneous bond dissociation reactions to form atomic O, C, S, F, and N species adsorbed on the surface of the Mg anode. These products indicate that the initial SEI layer formed on the surface of the pristine Mg anode consists of a complex mixture of multiple components such as oxides, carbides, sulfides, fluorides, and nitrides. We believe that the atomic level insights gained from these simulations will lay the groundwork for the rational design of tailored and functional interphases that are critical for the success of multivalent battery technology.

Agarwal, Garvit↗
NASA NTRS1997

Multi-Body Analysis of a Tiltrotor Configuration

The paper describes the aeroelastic analysis of a tiltrotor configuration. The 1/5 scale wind tunnel semispan model of the V-22 tiltrotor aircraft is considered. The analysis is performed by means of a multi-body code, based on an original formulation. The differential equilibrium problem is stated in terms of first order differential equations. The equilibrium equations of every rigid body are written, together with the definitions of the momenta. The bodies are connected by kinematic constraints, applied in form of Lagrangian multipliers. Deformable components are mainly modelled by means of beam elements, based on an original finite volume formulation. Multi-disciplinar problems can be solved by adding user-defined differential equations. In the presented analysis the equations related to the control of the swash-plate of the model are considered. Advantages of a multi-body aeroelastic code over existing comprehensive rotorcraft codes include the exact modelling of the kinematics of the hub, the detailed modelling of the flexibility of critical hub components, and the possibility to simulate steady flight conditions as well as wind-up and maneuvers. The simulations described in the paper include: 1) the analysis of the aeroelastic stability, with particular regard to the proprotor/pylon instability that is peculiar to tiltrotors, 2) the determination of the dynamic behavior of the system and of the loads due to typical maneuvers, with particular regard to the conversion from helicopter to airplane mode, and 3) the stress evaluation in critical components, such as the pitch links and the conversion downstop spring.

Ghiringhelli, G. L.↗
NASA NTRS2020

InSight Entry, Descent and Landing Pre-Flight Performance Predictions

On November 26, 2018, the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) lander successfully touched down on the surface of Mars. Over its seven-plus year development, NASA Langley Research Center’s (LaRC) Program to Optimize Simulated Trajectories II (POST2) was used to assess the mission’s Entry, Descent and Landing (EDL) vehicle system performance against related requirements across the full range of possible environmental and spacecraft conditions. Much of the simulation code was derived from the Phoenix mission, for which this vehicle is very similar. The InSight six degree-of-freedom simulation included models for Mars atmosphere, gravity and digital elevation maps of the landing location. Additionally, vehicle specific aerodynamic, parachute, engine, navigation sensor, flight software and landing radar models were also included. A set of dispersions for each model, as well as for additional simulation input parameters, were also included in order to provide a statistical, Monte Carlo prediction of the EDL system performance. An overview of the pre-flight performance assessments completed, including the various simulation campaigns used, will be provided. Ultimately, this work was critical in the assessment of readiness for InSight launch. A brief description of the use of this simulation in support of flight operations is also discussed.

Robert W Maddock↗
NASA NTRS1991

Behavioral interactions across various aircraft types - Results of systematic observations of line operations and simulations

The NASA/UT Line/LOS checklist is designed to capture critical components of crew interaction. The behaviors deemed critical to flight crew interaction include briefings, communications, inquiry, assertion/advocacy, and decisions communicated and acknowledged. Data shows significant behavioral interaction differences as a function of aircraft type, indicating that crew size and technology level were at least partly driving that difference.

Clothier, Cathy C.↗
NASA NTRS2019

Fe(exp3) - A Monte Carlo Simulation Capability for Evaluation of New Air Traffic Management Concepts

The concepts of unmanned aircraft system traffic management (UTM) and urban air mobility (UAM) introduce high-density operations in low-altitude airspace and will change the paradigm of the traditional air traffic system. The Flexible engine for Fast-time evaluation of Flight environments (Fe3) provides the capability of statistically analyzing high-density, high-fidelity, and low-altitude traffic system without conducting infeasible and cost-prohibitive flight tests that involve a large volume of aerial vehicles. With this simulation capability, stakeholders can study the impacts of critical factors, define requirements, policies, and protocols needed to support a safe yet efficient traffic system, assess operational risks, and optimize flight schedules. This work provides an introduction to this simulation tool including its architecture and various models involved. Its performance and applications in high density air traffic operations are also presented.

Paralellization↗
DOE OSTI2026

Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations

Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

An Overview of Distributed Spacecraft Autonomy at NASA Ames

Autonomous decision-making significantly increases mission effectiveness by mitigating the effects of communication constraints, like latency and bandwidth, and mission complexity on multi-spacecraft operations. To advance the state of the art in autonomous Distributed Space Systems (DSS), the Distributed Spacecraft Autonomy (DSA) team at NASA's Ames Research Center is developing within five relevant technical areas: distributed resource and task management, reactive operations, system modeling and simulation, human-swarm interaction, and ad hoc network communications. DSA is maturing these technologies - critical for future large autonomous DSS - from concept to launch via simulation studies and orbital deployments. A 100-node heterogenous Processor-in-the-Loop (PiL) testbed aids distributed autonomy capability development and verification of multi-spacecraft missions. The DSA software payload deployed to the D-Orbit SCV-004 spacecraft demonstrates multi-agent reconfigurability and reliability as part of an ESA-sponsored in-orbit technology demonstration. Finally, DSA's primary flight mission showcases collaborative resource allocation for multipoint science data collection with four small spacecraft as a payload on NASA's Starling 1.0 satellites.

Caleb Ashmore Adams↗
NASA NTRSDate not supplied

Extending ISS Life Beyond 2030

The United States On-orbit Segment (USOS) of the International Space Station (ISS) was designed to meet a 15-year on-orbit life. Since the first hardware was launched in late 1998, the ISS would have reached its end of life in 2013. With the realization that the ISS would be needed well into the next decade, and beyond, a multi-disciplinary effort was undertaken to extend the ISS’ life through 2028 and to show that further extension to 2040 and beyond is not only feasible but achievable. Currently, NASA and the ISS international partners have agreed to extend its operations through 2030. This collaborative effort ensures that the ISS will continue to serve as a hub for scientific research, international cooperation, and educational endeavors for the next decade. Maintaining a continuous human presence in Low Earth Orbit (LEO) is desirable for testing new LEO, lunar, and deep-space technologies; conducting scientific research in micro-gravity for the benefit of life on Earth; and enabling a seamless transition of capabilities to one or more commercially owned and operated destinations. This paper provides an overview of the ISS life extension project with a particular focus on the analytical approach used to assess the primary structure. This analytical approach includes future operations planning, critical location screening, on-orbit dynamic load simulation, on-orbit optical property degradation studies, on-orbit thermal analyses, spectra generation, crack model idealization, fracture analyses, and post processing. Structural life results and identification of the most critical on-orbit events are presented. Also addressed are life extension approaches for other affected sub-systems, including: 1. Secondary Structure. 2. Materials. Evaluations consider environmental exposure to atomic oxygen, ionizing and gamma radiation, fluids, etc. Life limited materials, wear, and usage effects are also considered. 3. Environmental Control and Life Support Systems, including oxygen supply and generation, water recovery and management, and regenerative hardware. Evaluations are performed to determine which hardware can be run to failure and which are assessed for life extension. 4. Electrical power system. ISS is powered by eight channels of solar arrays and an electrical energy storage system providing 357 kWh power. Power generation and balance analyses are performed considering hardware degradation and increasing power demand. 5. Logistics and maintenance. Analyses are performed to determine critical spares required to maintain functionality. Consideration is given to supply chain health, obsolescence issues, onboard stowage availability, and up-mass capability. The successful life extension results have built confidence to safely operate, maintain and enhance the ISS well beyond the current decade. Extending the operational life of the ISS maintains an international presence in LEO and serves to avoid a gap in capability necessary to fulfill the exploration and research needs of NASA, international partners, and industry without interruption until a commercial space station is operational.

design life↗
NASA NTRS2015

F-16XL Hybrid Reynolds-Averaged Navier-Stokes/Large Eddy Simulation on Unstructured Grids

This study continues the Cranked Arrow Wing Aerodynamics Program, International (CAWAPI) investigation with the FUN3D and USM3D flow solvers. CAWAPI was established to study the F-16XL, because it provides a unique opportunity to fuse fight test, wind tunnel test, and simulation to understand the aerodynamic features of swept wings. The high-lift performance of the cranked-arrow wing planform is critical for recent and past supersonic transport design concepts. Simulations of the low speed high angle of attack Flight Condition 25 are compared: Detached Eddy Simulation (DES), Modi ed Delayed Detached Eddy Simulation (MDDES), and the Spalart-Allmaras (SA) RANS model. Iso- surfaces of Q criterion show the development of coherent primary and secondary vortices on the upper surface of the wing that spiral, burst, and commingle. SA produces higher pressure peaks nearer to the leading-edge of the wing than flight test measurements. Mean DES and MDDES pressures better predict the flight test measurements, especially on the outer wing section. Vorticies and vortex-vortex interaction impact unsteady surface pressures. USM3D showed many sharp tones in volume points spectra near the wing apex with low broadband noise and FUN3D showed more broadband noise with weaker tones. Spectra of the volume points near the outer wing leading-edge was primarily broadband for both codes. Without unsteady flight measurements, the flight pressure environment can not be used to validate the simulations containing tonal or broadband spectra. Mean forces and moment are very similar between FUN3D models and between USM3D models. Spectra of the unsteady forces and moment are broadband with a few sharp peaks for USM3D.

Park, Michael A.↗
DOE OSTI2024

Influence of Real Gas Effects on Chemical Kinetics in Oxycombustion in Supercritical Carbon Dioxide

Oxycombustion in supercritical carbon dioxide is an integral part of the Allam Cycle, a technology that enables carbon-neutral use of fossil-fuels and carbon-negative use of biofuels. We simulate oxycombustion in a realistic combustor geometry for two sets of fuel conditions: pure methane, and a 40 percent methane/60 percent carbon dioxide blend, both at 343.15 K. The fuel jet mixes with a preheated swirler of 20 percent oxygen and 80 percent carbon dioxide at 1005.35 K, with a 100 percent carbon dioxide coflow at 783.15 K. The entire system operates at a pressure of 300 bar, putting the entire system above the critical temperature and pressure of carbon dioxide. Simulations are performed using PeleC, a compressible block-structured adaptive-mesh refinement (AMR) reacting flow code. Direct comparison of results using thermodynamically self-consistent implementations of the ideal gas equation of state and the Soave-Redlich-Kwong (SRK) equation of state allow quantification of the real gas impacts on combustion kinetics, flow temperatures, and pollutant formation.

carbon dioxide↗
NASA NTRS1990

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-2

Jet-A spray combustion has been evaluated in gas turbine combustion with the use of propane chemical kinetics as the first approximation for the chemical reactions. Here, the numerical solutions are obtained by using the KIVA-2 computer code. The KIVA-2 code is the most developed of the available multidimensional combustion computer programs for application of the in-cylinder combustion dynamics of internal combustion engines. The released version of KIVA-2 assumes that 12 chemical species are present; the code uses an Arrhenius kinetic-controlled combustion model governed by a four-step global chemical reaction and six equilibrium reactions. Researchers efforts involve the addition of Jet-A thermophysical properties and the implementation of detailed reaction mechanisms for propane oxidation. Three different detailed reaction mechanism models are considered. The first model consists of 131 reactions and 45 species. This is considered as the full mechanism which is developed through the study of chemical kinetics of propane combustion in an enclosed chamber. The full mechanism is evaluated by comparing calculated ignition delay times with available shock tube data. However, these detailed reactions occupy too much computer memory and CPU time for the computation. Therefore, it only serves as a benchmark case by which to evaluate other simplified models. Two possible simplified models were tested in the existing computer code KIVA-2 for the same conditions as used with the full mechanism. One model is obtained through a sensitivity analysis using LSENS, the general kinetics and sensitivity analysis program code of D. A. Bittker and K. Radhakrishnan. This model consists of 45 chemical reactions and 27 species. The other model is based on the work published by C. K. Westbrook and F. L. Dryer.

Nguyen, H. L.↗
NASA NTRS1990

Critical evaluation of Jet-A spray combustion using propane chemical kinetics in gas turbine combustion simulated by KIVA-II

Numerical solutions of the Jet-A spray combustion were obtained by means of the KIVA-II computer code after Jet-A properties were added to the 12 chemical species the program had initially contained. Three different reaction mechanism models are considered. The first model consists of 131 reactions and 45 species; it is evaluated by comparing calculated ignition delay times with available shock tube data, and it is used in the evaluation of the other two simplified models. The simplified mechanisms consider 45 reactions and 27 species and 5 reactions and 12 species, respectively. In the prediction of pollutants NOx and CO, the full mechanism of 131 reactions is considered to be more reliable. The numerical results indicate that the variation of the maximum flame temperature is within 20 percent as compared with that of the full mechanism of 131 reactions. The chemical compositions of major components such as C3H8, H2O, O2, CO2, and N2 are of the same order of magnitude. However, the concentrations of pollutants are quite different.

Nguyen, H. L.↗