Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Hybrid Reduced Order Model for N2-N2 Interactions for Application to Dissociation and Energy Transfer Processes

Recent work in the aerothermodynamics community has focused on the development of reduced-order models for thermo-chemical non-equilibrium which avoid the restrictive assumptions of multi-temperature models and the prohibitive cost associated with State-to-State (StS) models. In the present work, this is accomplished by lumping energy states together and assuming groups of states are roughly in equilibrium. As a result, the non-equilibrium behavior of a gas can be captured at a reduced computational cost from a full StS simulation. In this work, we present a hybrid grouping model for studying energy transfer and dissociation in a mixture of nitrogen molecules due to N2-N2 reactions. This is accomplished by making use of a grouping strategy informed by data from the N2-N StS kinetic data. However, due to the massive computational cost associated with constructing StS data for the N2-N2 system, the kinetic data for the hybrid grouping model are calculated using the quasi-classical trajectory (QCT) method by sampling states for trajectory within the groups. This general framework is called the Maximum-Entropy Quasi-Classical Trajectory (ME-QCT) method. The primary challenge associated with this method is that rates for reverse grouped reactions cannot be obtained through detailed balance at a group level, due to the variation of group internal temperatures. To construct the full model for N2-N2 grouped kinetics using the ME-QCT method, detailed balance is invoked at the microscopic level, allowing for the calculation of the full kinetic data from QCT. Results will be presented using the full ME-QCT model for the N2-N2 system in an isothermal and isochoric reactor simulation. In addition, simple CFD test cases for a one-dimensional standing shock and for a quasi-one-dimensional nozzle will be used for demonstration of the ME-QCT method. This method allows for the calculation of non-equilibrium behavior for the N2-N2 system without the prohibitive cost of a full StS simulation. Moreover, it enables the construction of a unified model for the dissociating and recombining non-equilibrium flows.

nonequilibrium modeling↗
DOE OSTI2025

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗
NASA NTRS1990

NASA/MSFC Large Space Structures Ground Test Facility

The NASA/MFSC Large Space Structures Ground Test Facility (LSS GTF) is described in terms of the testing, evaluation, and implementation of control and system identification techniques for typical large space structures. The GTF comprises Control, Astrophysics, and Structures Experiment in Space (CASES) GTF which is being developed and the operational Single Structure Control (SSC) laboratory (an LSS flexible beam suspended vertically with sensor and actuator systems, a real-time computer system, a disturbance system, and an optical pointing system). The configuration of the laboratory and the systems used are set forth in terms of monitoring simulated disturbances. The Shuttle-based CASES experiment involves a 105-foot boom as part of an X-ray experiment, and the control of this structure will be tested at the CASES GTF. The actuation, measurement, sensor, and computer systems are described, and the configuration of the GTF is given.

Jones, Victoria L.↗
DOE OSTI2022

Internal Model-Based Active Damping Strategy for a Back-to-Back Modular Multilevel Converter System for Advanced Grid Support: Preprint

The proposed work focuses on the possibility of achieving the reduction in the size of the interfacing filter for the grid connection for a back-to-back modular multilevel converter system with the use of third order LCL filters. To achieve the same attenuation, it is possible to reduce the size of the interfacing filters by the usage of third order filters. However, this kind of filtering comes with the limitation of having sustained oscillations due to lack of damping especially during transient changes. Therefore, in this work, a methodology has been proposed based on the principle of internal model to cater for this unwanted oscillations. The third order system is modeled inside the microcontroller with the damping enabled. The error between the output from the model and the actual are compared and the error is utilized to accomplish the active damping strategy. The proposed architecture is verified via computer simulations based on MATLAB/Simulink domain and various case study results along with their discussion is presented in this paper.

dq control↗
DOE OSTI2025

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗
NASA NTRS1992

A hybrid flow model for charged and neutral particles around spacecraft in Low Earth Orbit

A flow simulation for charged and neutral particles around spacecraft is given that combines the fully coupled effects of neutral particle flow, plasma flow, electromagnetic field effects, and spacecraft charging. The simulation allows for chemically reacting flows (associative ionization, dissociation, and chain exchange) and for thermal accommodation at a spacecraft surface. The solution procedure was applied to several types of flows with orbital altitudes ranging from 250 km to 500 km. Excellent agreement with previously verified computational algorithms was obtained for both the neutral and charged particle cases. Results for the full flow simulation are presented for altitudes of 250 km and 500 km.

Justiz, Charles R.↗
NASA NTRS1993

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 323 simulation (Run B). Other than the grid size, the primary difference between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. The numerical results indicate that gravitational collapse can proceed in a variety of ways. In the Run A, collapse led to an elongated tube-like structure, while in the Run B, collapse led to a flatter, disklike structure.

Shebalin, John V.↗
NASA NTRS2006

Fluorescence Imaging of Underexpanded Jets and Comparison with CFD

An experimental study of underexpanded and highly underexpanded axisymmetric nitrogen free jets seeded with 0.5% nitric oxide (NO) and issuing from a sonic orifice was conducted at NASA Langley Research Center. Reynolds numbers based on nozzle exit conditions ranged from 770 to 35,700, and nozzle exit-to-ambient jet pressure ratios ranged from 2 to 35. These flows were non-intrusively visualized with a spatial resolution of approximately 0.14 mm x 0.14 mm x 1 mm thick and a temporal resolution of 1 s using planar laser-induced fluorescence (PLIF) of NO, with the laser tuned to the strongly-fluorescing UV absorption bands of the Q1 band head near 226.256 nm. Three laminar cases were selected for comparison with computational fluid dynamics (CFD). The cases were run using GASP (General Aerodynamic Simulation Program) Version 4. Comparisons of the fundamental wavelength of the jet flow showed good agreement between CFD and experiment for all three test cases, while comparisons of Mach disk location and Mach disk diameter showed good agreement at lower jet pressure ratios, with a tendency to slightly underpredict these parameters with increasing jet pressure ratio.

Wilkes, Jennifer A.↗
DOE OSTI2022

PythonFOAM: In-situ data analyses with OpenFOAM and Python

Here, we outline the development of a general-purpose Python-based data analysis tool for OpenFOAM. Our implementation relies on the construction of OpenFOAM applications that have bindings to data analysis libraries in Python. Double precision data in OpenFOAM is cast to a NumPy array using the NumPy C-API and Python modules may then be used for arbitrary data analysis and manipulation on flow-field information. We highlight how the proposed wrapper may be used for an in-situ online singular value decomposition (SVD) implemented in Python and accessed from the OpenFOAM solver PimpleFOAM. Here, 'in-situ' refers to a programming paradigm that allows for a concurrent computation of the data analysis on the same computational resources utilized for the partial differential equation solver. In addition, to demonstrate parallel deployments, we deploy a distributed SVD, which collects snapshot data across the ranks of a distributed simulation to compute the global left singular vectors. Crucially, both OpenFOAM and Python share the same message passing interface (MPI) communicator for this deployment which allows Python objects and functions to exchange NumPy arrays across ranks. Subsequently, we provide scaling assessments of this distributed SVD on multiple nodes of Intel Broadwell and KNL architectures for canonical test cases such as the large eddy simulations of a backward facing step and a channel flow at friction Reynolds number of 395. Finally, we demonstrate the deployment of a deep neural network for compressing the flow-field information using an autoencoder to demonstrate an ability to use state-of-the-art machine learning tools in the Python ecosystem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Comparison of measurement and simulation of ATLAS cavern radiation background

Sixteen Medipix2 pixel detector based (MPX) devices were operated at various positions within the ATLAS detector and cavern continuously from early 2008 up to 2013. In addition to photons, each MPX detector is capable to detect charged particles, and neutrons as it is covered with a mask of converter materials dividing its area into regions sensitive to thermal or fast neutrons. The MPX detector network was effectively used for real-time measurements of the spectral characteristics and composition of complex radiation fields in ATLAS. This article reports comparison of the results of measurements performed with MPX detectors during the LHC operation period in 2010 and 2011 with Monte Carlo simulations results from the FLUGG and GCALOR codes. For the purpose of this comparison, the MPX detectors were operated in tracking mode with low threshold (8–10 keV) allowing one to distinguish among particle categories based on the recognition of track patterns left by the particles in the MPX sensitive layer. The comparison of measurements with simulations shows that the agreement between measured and simulated data is satisfactory in most cases within a factor of two.

47 OTHER INSTRUMENTATION↗
NASA NTRS1990

Bounded free shear flows - Linear and nonlinear growth

The physics of bounded free shear flows is studied at supersonic Mach numbers. Wall locations are at 1500 momentum thicknesses away from the shear layer center line. Intrinsic to this flow are high multiplicities of supersonic modes not present in the unbounded case. It is seen that these supersonic modes are resistant to subharmonic growth, a necessary precursor for efficient mixing. Previous results (Eberhardt et al., 1988) with wall locations at less than one half momentum thickness away from the shear layer center indicate that modes resembling channel modes (maximum amplification at the walls) did not roll up in numerical simulations. The similarity between these two cases is pointed out.

Macaraeg, Michele G.↗
DOE OSTI2020

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

A Systems Approach to AI Model Integration and Performance Evaluation for the Generic UAM Simulation Framework

This paper introduces py-guam, an open-source experimentation framework developed for the NASA Generic Urban Air Mobility simulation (GUAM) environment, facilitating the integration and evaluation of advanced artificial intelligence (AI) algorithms. We present a systems approach which enables the seamless incorporation of data-driven models, including off-nominal and failure state detection, into the GUAM’s Cognitive Architecture (CA). The framework supports customizable experimentation parameters, derives Safety Performance Indicators (SPIs) from UL 4600 safety case analyses, and employs rapid UAM simulations to assess AI impacts on flight performance across diverse scenarios. Through comprehensive testing and validation experiments, we demonstrate GUAM’s capability to enhance safety and efficiency in urban air mobility operations. Additionally, the open-source nature of py-guam fosters community collaboration, ensuring continuous improvement and adaptability to evolving technological advancements. This work establishes a robust tool for developing and testing AI-driven urban air mobility (UAM) systems, advancing the safety and reliability of autonomous urban air vehicles.

Artificial Intelligence↗
DOE OSTI2025

Data-model files associated with the manuscript "Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California"

This package contains the data, simulation setups, notebooks and figures used in “Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California” (Xu et al., 2025). In this study, we selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise on coastal hydrology using the process-based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site-specific data. We designed a novel modeling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modeled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected sea level rise over the next few decades. The low-lying Yampah Marsh, the reference site, is likely to be inundated due to future sea level rise when highest tides arrive; while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.Several files can be found from this data package.1. README.md: This file describes the title, journal, co-authors, abstract, repository structure and model version.2. Simulation_Setups.zip: The file contains the model configuration files (XML format) for ATS. 3. Notebooks.zip: The file contains the Jupyter notebooks for generating the pre- and post-restoration meshes and the meshes of future scenarios. 4. Figures.zip: The file contains the figures used in the manuscript.5. Data.zip: The file contains the data used to drive the model simulations, including watershed and wetlands boundaries, mesh files and references to additional datasets (e.g., meteorological forcing, tidal dataset, DEMs, land cover, soil properties). Also, it contains water level observations at the restored wetland.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Protection Against Graph-Based False Data Injection Attacks on Power Systems

Graph signal processing (GSP) has emerged as a powerful tool for practical network applications, including power system monitoring. By representing power system voltages as smooth graph signals, recent research has focused on developing GSP-based methods for state estimation, attack detection, and topology identification. Included, efficient methods have been developed for detecting false data injection (FDI) attacks, which until now were perceived as non-smooth with respect to the graph Laplacian matrix. Consequently, these methods may not be effective against smooth FDI attacks. In this paper, we propose a graph FDI (GFDI) attack that minimizes the Laplacian-based graph total variation (TV) under practical constraints. In addition, we develop a low-complexity algorithm that solves the non-convex GDFI attack optimization problem using ell_1-norm relaxation, the projected gradient descent (PGD) algorithm, and the alternating direction method of multipliers (ADMM). We then propose a protection scheme that identifies the minimal set of measurements necessary to constrain the GFDI output to high graph TV, thereby enabling its detection by existing GSP-based detectors. Our numerical simulations on the IEEE-57 bus test case reveal the potential threat posed by well-designed GSP-based FDI attacks. Moreover, we demonstrate that integrating the proposed protection design with GSP-based detection can lead to significant hardware cost savings compared to previous designs of protection methods against FDI attacks.

Morgenstern, Gal↗
DOE OSTI2020

Moment preserving constrained resampling with applications to particle-in-cell methods

The Moment Preserving Constrained Resampling (MPCR) algorithm for particle resampling is introduced and applied to particle-in-cell (PIC) methods to increase simulation accuracy, reduce compute cost, and/or avoid numerical instabilities. The general algorithm partitions the system space into smaller subsets and resamples the distribution within each subset. Further, the algorithm is designed to conserve any number of particle and grid moments with a high degree of accuracy (i.e. machine accuracy). The effectiveness of MPCR is demonstrated with several numerical tests, including a use-case study in gyrokinetic fusion plasma simulations. Finally, the computational cost of MPCR is negligible compared to the cost of particle evolution in PIC methods, and the tests demonstrate that periodic particle resampling yields a significant improvement in the accuracy and stability of the results.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Decoupling of peeling and ballooning thresholds for pedestal stability and reduction in ELM frequency via enhanced turbulence with edge electron cyclotron heating in DIII-D

Abstract The edge localized mode (ELM) frequency ( f ELM ) decreased by 63% when electron cyclotron heating (ECH) deposition location is shifted from ρ = 0.4 to ρ = 0.8 in DIII-D discharges where the power ratio between neutral beam injection (NBI) and ECH ( P NBI / P ECH ) is kept at ∼1. The performance of the pedestal in the ECH heated case is compared with a pure NBI reference discharge while keeping the total input power constant. All these discharges are performed at balanced input torque conditions. Furthermore, in the pure NBI discharge a strong decoupling of the peeling–ballooning (PB) thresholds is observed. The PB decoupling is preserved when the ECH is deposited at ρ = 0.8 and P NBI / P ECH ∼ 1, while the thresholds manifest a closed stability boundary when the ECH is deposited at ρ = 0.4. The inter-ELM pedestal recovery time is considerably larger for the ECH at ρ = 0.8 case. Increased pedestal turbulence is observed in beam emission spectroscopy (BES), Doppler backscattering and magnetic diagnostics for the ECH at the ρ = 0.8 case. Strong growth of a TEM-like mode is observed in BES and the mode growth is correlated with the decrease in f ELM . In view of these observations, the increased pedestal turbulence seems to be the plausible reason behind the delayed pedestal recovery following an ELM event in the ECH at ρ = 0.8 case, and the preservation of PB decoupling through temperature pedestal profile widening. TRANSP interpretative simulations show that the ECH at the ρ = 0.8 case is more susceptible to ITG/TEM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗