SEARCH · Engineering Papers

Results for “Simulation framework”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Possible routes to a strong science case for a joint CLEoPATRA-Roman microlensing survey

The Contemporaneous LEnsing Parallax and Autonomous TRansient Array (CLEoPATRA) is a mis sion concept for a space-based microlensing survey, proposed to enhance science for the Roman Space Telescope (RST) microlensing survey. From CLEoPATRA’s planned vantage point in geosynchronous equatorial orbit (GEO), it will observe microlensing events simultaneously and from a different point of view to RST at L2. This joint L2-GEO observation baseline gives a projected angular separation of order the size of the characteristic angular lensing scale of Free-Floating-Planets (FFPs) in the Galaxy, making it powerful for constraining FFP properties. We find that a strong science case will demonstrate which parts of the FFP parameter space a joint RST-CLEoPATRA survey will give lead ing lens constraints compared with competing ground- and space- based RST observatory companions (e.g., PRIME, Euclid). Both the simulations and framework that connect telescope capabilities to physical constraints already exist for RST FFPs and can be used to demonstrate the power of joint RST-CLEoPATRA survey, with the addition of a simulated CLEoPATRA event blending distribution. Population-level FFP constraints could also be calculated by adapting ongoing black hole microlensing research at Lawrence Livermore National Laboratory. We also find that considering CLEoPATRA as a third space-based observatory (in addition to RST and Euclid/PRIME), which will break microlensing degeneracies leading to unambiguous FFP constraints, is worth investigating as part of a science case. Finally, we highlight possible auxiliary CLEoPATRA science cases including relieving tension between the RST Core Community Surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

Large-Scale Simulation of a Distributed Sensing Network Supporting Regional Urban Air Mobility Operations

Urban Air Mobility (UAM) is set to transform transportation in densely populated regions like the San Francisco Bay Area. This paper introduces an innovative simulation approach to explore large-scale UAM scenarios, emphasizing the use of distributed sensing to enhance operational efficiency and safety. The Revolutionary Vertical Lift Technology (RVLT) model is employed as the framework for simulating complex interactions among multiple vehicles within urban landscapes. Strategically deployed ground sensor nodes enable distributed sensing, enhancing situational awareness and operational effectiveness. By integrating empirical data and geographical realism, the simulations provide a systematic analysis of the feasibility, efficiency, and safety considerations associated with UAM deployment in urban environments. Factors such as air traffic density and infrastructural requirements are thoroughly examined, offering actionable insights for policymakers and industry stakeholders. This paper aims to refine the structure and scenarios for large-scale simulations based on distributed sensing, thereby contributing to the advancement of UAM operations.

UAM↗

NASA NTRSDate not supplied

Large-Scale Simulation of a Distributed Sensing Network Supporting Regional Urban Air Mobility Operations

Aircraft Mobility↗

NASA NTRS1993

An RL10A-3-3A rocket engine model using the rocket engine transient simulator (ROCETS) software

Steady-state and transient computer models of the RL10A-3-3A rocket engine have been created using the Rocket Engine Transient Simulation (ROCETS) code. These models were created for several purposes. The RL10 engine is a critical component of past, present, and future space missions; the model will give NASA an in-house capability to simulate the performance of the engine under various operating conditions and mission profiles. The RL10 simulation activity is also an opportunity to further validate the ROCETS program. The ROCETS code is an important tool for modeling rocket engine systems at NASA Lewis. ROCETS provides a modular and general framework for simulating the steady-state and transient behavior of any desired propulsion system. Although the ROCETS code is being used in a number of different analysis and design projects within NASA, it has not been extensively validated for any system using actual test data. The RL10A-3-3A has a ten year history of test and flight applications; it should provide sufficient data to validate the ROCETS program capability. The ROCETS models of the RL10 system were created using design information provided by Pratt & Whitney, the engine manufacturer. These models are in the process of being validated using test-stand and flight data. This paper includes a brief description of the models and comparison of preliminary simulation output against flight and test-stand data.

Binder, Michael↗

NASA NTRSDate not supplied

Transonic Lift and Drag Predictions using Wall Modelled Large Eddy Simulations

Wall-modelled Large Eddy Simulations (WMLES) of the NASA Common Research Model (CRM) at transonic conditions and various angles of attacks leading up to and including shock-induced flow separation are performed using the LAVA computational framework. The simulations are shown to accurately predict the lift curve slope and the onset of separation characterized by the break in the pitching moment. Furthermore, careful assessment of skin-friction drag at cruise condition and its subsequent decrease with increasing angles of attack are shown to be in agreement with viscous sublayer resolving Reynolds Averaged Navier Stokes(RANS) simulations. The small differences between WMLES and RANS appear to be of the same order as differences seen between two RANS models at the cruise-point. Some sensitivity is observed to the coefficient used in the subgrid scale model, although this can be reconciled by noting the low chord- and shock-incidence Reynolds numbers (Rex≈106) seen in the outboard regions of the wing, along with the uncertainties associated with tripping and the numerical transition that occurs near the leading edge. Among the primary drawbacks seen in the WMLES predictions is the subdued increase in wave-drag with increasing angle-of-attack and the shock intensity when compared with experimental data. Preliminary assessment suggests that additional span- and stream-wise grid resolution is likely needed in mid- and out-board portions of the wing to better resolve the shock-induced separated flow and to further investigate prediction accuracy of unsteady temporal characteristics of the problem.

ARMD↗

NASA NTRS1991

Numerical simulations of thermal detection of disbonds in lap joints

A 2D time-dependent FEM heat-transfer algorithm that is unconditionally stable and implicit is presently used to solve for the time-dependence of a structure's temperature. By examining data-reduction techniques in the framework of simulation data, the contrast enhancement afforded by a given technique can be evaluated without regard to experimental error. The present computational simulations demonstrate that the time-derivative technique can clearly delineate the substructure of lap-joint geometries; in addition, histograms of the front surface images are able to quantitatively measure contrast-enhancements for various inspection and data-reduction techniques.

Howell, P. A.↗

NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗

DOE OSTI2022

Logistics simulation of a remediation effort for a hypothetical radiological contamination scenario

To mitigate the effects following a large-scale nuclear or radiological material release in an urban environment and to expedite recovery, the Integrated Wash-Aid Treatment Emergency Reuse System (IWATERS) was developed. IWATERS consists of three operations: washing contaminated surfaces with an ionic wash solution, collecting, and treating the contaminated wash solution on-site to remove contaminants, and reusing the treated solution throughout operations to preserve the clean water resource. This study develops a framework to simulate the logistics of IWATERS deployment, thereby gaining an understanding of the timeline for decontamination operations. For this purpose, the Analysis of Mobility Platform and GoldSim were leveraged for a hypothetical contamination scenario covering 65,200 m 2 of an urban center. The framework reveals that remediation progress is limited by several resources, notably the availability of vermiculite, a reactive clay that is required to treat the contaminated wash solution. Further, this study also presents how the simulation approach can be used to characterize alternatives to reduce the influence of limited resources on operational progress. Overall, this work lays the foundation for evaluating different decontamination methods through detailed logistics simulation, i.e., by refining simulation assumptions and expanding the range of scenarios the simulation can depict.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2023

Application of community data to surface complexation modeling framework development: Iron oxide protolysis

This study presents a comprehensive community data-driven surface complexation modeling framework for simulating potentiometric titration of mineral surfaces. Compiled community data for ferrihydrite, goethite, hematite, and magnetite are fit to produce representative protolysis constants that can reproduce potentiometric titration data collected from multiple literature sources. Using this framework, the impact of surface complexation model type and surface site density (SSD) on the fit quality and protolysis constants can be readily evaluated. For example, the non-electrostatic model yielded a poor data fit compared to diffuse double layer model and constant capacitance models due to the absence of known surface charge effects. Regardless of the choice of iron oxide mineral, pK a1 decreased with increasing SSD while the opposite tendency was observed for pK a2 . This newly developed framework demonstrates a method to reconcile community data-wide potentiometric titration data using Findable, Accessible, Interoperable, Reusable data principles to produce mineral protolysis constants that improve robustness of surface complexation models for applications in metal sorption and reactive transport modeling. The framework is readily expandable (as community data increase) and extensible (as the number of minerals increase). The framework provides a path forward for developing self-consistent, comprehensive, and updateable surface complexation databases for surface complexation and reactive transport modeling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

An uncertainty visualization framework for large-scale cardiovascular flow simulations: A case study on aortic stenosis

We present a generalizable uncertainty quantification (UQ) and visualization framework for lattice Boltzmann method simulations of high Reynolds number vascular flows, demonstrated on a patient-specific stenosed aorta. The framework combines EasyVVUQ for parameter sampling with large-eddy simulation turbulence modeling in HemeLB, and executes ensembles on the Frontier exascale supercomputer. Spatially resolved metrics, including entropy and isosurface-crossing probability, are used to map uncertainty in pressure and wall shear stress fields directly onto vascular geometries. Two sources of model variability are examined: inlet peak velocity and the Smagorinsky constant. Inlet velocity variation produces high uncertainty downstream of the stenosis where turbulence develops, while upstream regions remain stable. Smagorinsky constant variation has little effect on the large-scale pressure field but increases WSS uncertainty in localized high-shear regions. In both cases, the stenotic throat manifests low entropy, indicative of robust identification of elevated WSS. By linking quantitative UQ measures to three-dimensional anatomy, the framework improves interpretability over conventional 1D UQ plots and supports clinically relevant decision-making, with broad applicability to vascular flow problems requiring both accuracy and spatial insight.

Hemodynamics↗

DOE OSTI2025

Molecular Simulation of Functionalized Covalent Organic Framework Membranes for Inorganic Salt Separation

Covalent organic frameworks (COFs) enable molecular-level design of nanochannels for selective separation in pressure-driven membrane processes. Through variation of building blocks and, subsequently, the pore structure and chemistry, membrane performance can be tailored. This study employs nonequilibrium molecular dynamics simulations to theoretically demonstrate the tunable selectivity of COF membranes through a bottom-up functionalization approach. Water and salt transport are evaluated for six β-ketoenamine-linked COFs with varying multilayer thicknesses. For the thinnest multilayer (0.64–0.72 nm), all COFs exhibit low sodium sulfate (Na 2 SO 4 ) rejection (55–66%). However, 20 stacked sheets (6.4–7.2 nm) provide 67–98% rejection, with the sulfonated COF providing the highest Na 2 SO 4 rejection. Analysis of time-resolved ion density profiles reveals that solute rejection is primarily governed by interfacial exclusion arising from pore size and functional group chemistry. Although increasing salt rejection compromises water permeance, the permeance of all COF membranes is at least two orders of magnitude greater than that of a commercially available nanofiltration membrane. Overall, this work guides the rational design of COF membranes for aqueous salt separation.

Nanofiltration↗

DOE OSTI2024

Error-correcting codes for fermionic quantum simulation

Utilizing the framework of \mathbb{Z}_2 ℤ 2 lattice gauge theories in the context of Pauli stabilizer codes, we present methodologies for simulating fermions via qubit systems on a two-dimensional square lattice. We investigate the symplectic automorphisms of the Pauli module over the Laurent polynomial ring. This enables us to systematically increase the code distances of stabilizer codes while fixing the rate between encoded logical fermions and physical qubits. We identify a family of stabilizer codes suitable for fermion simulation, achieving code distances of d=2,3,4,5,6,7, allowing correction of any \lfloor \frac{d-1}{2} \rfloor ⌊ d − 1 2 ⌋ -qubit error. In contrast to the traditional code concatenation approach, our method can increase the code distances without decreasing the (fermionic) code rate. In particular, we explicitly show all stabilizers and logical operators for codes with code distances of d=3,4,5. We provide syndromes for all Pauli errors and invent a syndrome-matching algorithm to compute code distances numerically.

Physics↗

DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗

NASA NTRS1990

Radiatively driven stratosphere-troposphere interactions near the tops of tropical cloud clusters

Results are presented of two numerical simulations of the mechanism involved in the dehydration of air, using the model of Churchill (1988) and Churchill and Houze (1990) which combines the water and ice physics parameterizations and IR and solar-radiation parameterization with a convective adjustment scheme in a kinematic nondynamic framework. One simulation, a cirrus cloud simulation, was to test the Danielsen (1982) hypothesis of a dehydration mechanism for the stratosphere; the other was to simulate the mesoscale updraft in order to test an alternative mechanism for 'freeze-drying' the air. The results show that the physical processes simulated in the mesoscale updraft differ from those in the thin-cirrus simulation. While in the thin-cirrus case, eddy fluxes occur in response to IR radiative destabilization, and, hence, no net transfer occurs between troposphere and stratosphere, the mesosphere updraft case has net upward mass transport into the lower stratosphere.

Churchill, Dean D.↗

NASA NTRS1993

Gas turbine system simulation: An object-oriented approach

A prototype gas turbine engine simulation has been developed that offers a generalized framework for the simulation of engines subject to steady-state and transient operating conditions. The prototype is in preliminary form, but it successfully demonstrates the viability of an object-oriented approach for generalized simulation applications. Although object oriented programming languages are-relative to FORTRAN-somewhat austere, it is proposed that gas turbine simulations of an interdisciplinary nature will benefit significantly in terms of code reliability, maintainability, and manageability. This report elucidates specific gas turbine simulation obstacles that an object-oriented framework can overcome and describes the opportunity for interdisciplinary simulation that the approach offers.

Drummond, Colin K.↗

NASA NTRS1990

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗