Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

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,621 records · Page 90

Simulation of the 1979 spring bloom in the Mid-Atlantic Bight - A coupled physical/biological/optical model

A coupled physical/biological/optical model is developed for studies of phytoplankton variability in the spring 1979 Mid-Atlantic Bight, as observed by CZCS imagery. The model incorporates advection, mixing, sinking, growth as a function of light, temperature, nutrient availability, and death as a function of ingestion. It produced chlorophyll concentrations within the first attenuated depth within 1 standard deviation of CZCS imagery on large scale. The primary production estimates obtained using this model were within reasonable agreement with those measured in situ.

Gregg, Watson W.↗

Multiscale CFD simulation of biomass fast pyrolysis with a machine learning derived intra-particle model and detailed pyrolysis kinetics

Coupling particle and reactor scale models is as essential as reactor fluid dynamics and particle motion for accurate Computational Fluid Dynamic (CFD) simulations of biomass fast pyrolysis reactors due to intraparticle heat transfer and chemical reactions controlling conversion time and product distributions. Direct online coupling of a particle model with a reactor model is computationally expensive, while offline coupling is case-dependent. In this research, solutions from a series of particle pyrolysis simulations were regressed with Artificial Neural Network (ANN). Furthermoer, this machine learning-derived model predicted the same temperature and conversion profiles compared with particle resolved simulation while the isothermal approach overpredicted the temperature by 130 K and underpredicted the conversion time by 30 s. The ANN model was then integrated into CFD simulations of fluidized bed biomass fast pyrolysis with varied feedstocks via coupling PyTorch and MFiX. The averaged error of simulation predicted bio-oil yields with four feedstocks is 6.4%. This multi-scale approach provides an efficient tool for the coupled particle and reactor scale simulations of biomass pyrolysis.

09 BIOMASS FUELS↗

Simulating Sulfate Aerosol Size in the Hunga Tonga Plume: Comparisons with Observations

Using the NASA GEOS Chemistry Climate Model (CCM) coupled to a detailed sectional aerosol microphysical model, we have simulated the January 2022 Hunga Tonga-Hunga Ha’apai (HTHH) plume. We have developed a version of the Global Modeling Initiative (GMI) chemistry mechanism to include the requisite sulfur chemistry and coupled it to the Community Aerosol and Radiation Model for Atmospheres (CARMA) sectional aerosol microphysics within the GEOS CCM. Our simulations focus on the lifecycle of stratospheric sulfur injected by the eruption–from the oxidation of the injected sulfur dioxide to the nucleation, growth, coagulation, hydration, and eventual settling of the volcanic sulfate aerosol particles. In order to constrain the simulated microphysics, we focus on comparing our model to observations with aerosol size information, namely in-situ particle counter observations and multi-wavelength satellite observations from OMPS-LP and SAGE-III. This eruption provides an interesting test case for our model; both the modeled gas-phase chemistry as well as the modeled aerosol microphysics are perturbed by the massive injection of water that accompanied the primary sulfur emissions from the eruption.

aerosols↗

ToPolyAgent: AI agents for coarse-grained bead-spring topological polymer simulations

We introduce ToPolyAgent, a multi-agent AI framework for performing coarse-grained molecular dynamics (MD) simulations of topological polymers through natural language instructions. By integrating large language models (LLMs) with domain-specific computational tools, ToPolyAgent supports both interactive and autonomous simulation workflows across diverse polymer architectures, including linear, ring, brush, and star polymers, as well as dendrimers. The system consists of four LLM-powered agents: a Config Agent for generating initial polymer–solvent configurations, a Simulation Agent for executing LAMMPS-based MD simulations and conformational analyses, a Report Agent for compiling markdown reports, and a Workflow Agent for streamlined autonomous operations. Interactive mode incorporates user feedback loops for iterative refinements, while autonomous mode enables end-to-end task execution from detailed prompts. We demonstrate ToPolyAgent's versatility through case studies involving diverse polymer architectures under varying solvent conditions, thermostats, and simulation lengths. Furthermore, we highlight its potential as a research assistant by directing it to investigate the effect of interaction parameters on the linear polymer conformation, and the influence of grafting density on the persistence length of the brush polymer. By coupling natural language interfaces with rigorous simulation tools, ToPolyAgent lowers barriers to complex computational workflows and advances AI-driven materials discovery in polymer science. It lays the foundation for autonomous and extensible multi-agent scientific research ecosystems.

Ding, Lijie [Oak Ridge National Laboratory (ORNL),↗

Octave-spanning terahertz quarter-wave plates based on over-coupled Fabry–Pérot resonances in reflective metal–dielectric–metal metasurfaces

Compact devices for broadband polarization control in the terahertz (THz) regime remain challenging to realize due to the intrinsic phase dispersion associated with birefringent materials and resonant structures. Here, we demonstrate achromatic THz quarter-wave plates with near-unity efficiency and operating under a wide range of incidence angles, based on over-coupled metal–dielectric–metal reflective metasurfaces. The devices operate as single-port anisotropic Fabry–Pérot cavities, in which the phase dispersion of the over-coupled resonances is engineered to maintain an approximately constant relative phase delay between orthogonal field components over a broad frequency range. Four complementary metasurface designs operating at an incidence angle of 45° collectively cover the 0.25–3 THz frequency range accessible to a typical THz time-domain spectroscopy system. Each device exhibits an approximately octave-wide operational bandwidth, with axial ratios below 3 dB and polarization-conversion efficiencies exceeding 80% across most of the operating band. Systematic structural optimization suppresses coupling to diffraction and guided-wave modes, thereby further extending the usable bandwidth while preserving the required phase relationship between orthogonal polarizations. The demonstrated metasurfaces are compatible with wafer-scale fabrication processes, and experimental results show excellent agreement with simulations. These findings establish over-coupled reflective metasurfaces as a robust, scalable, and versatile platform for broadband THz polarization control.

36 MATERIALS SCIENCE↗

ANOPP/VMS HSCT ground contour system

This viewgraph shows the integration of the Visual Motion Simulator with ANOPP. ANOPP is an acronym for the Aircraft NOise Prediction Program. It is a computer code consisting of dedicated noise prediction modules for jet, propeller, and rotor powered aircraft along with flight support and noise propagation modules, all executed under the control of an executive system. The Visual Motion Simulator (VMS) is a ground based motion simulator with six degrees of freedom. The transport-type cockpit is equipped with conventional flight and engine-thrust controls and with flight instrument displays. Control forces on the wheel, column, and rudder pedals are provided by a hydraulic system coupled with an analog computer. The simulator provides variable-feel characteristics of stiffness, damping, coulomb friction, breakout forces, and inertia. The VMS provides a wide range of realistic flight trajectories necessary for computing accurate ground contours. The NASA VMS will be discussed in detail later in this presentation. An equally important part of the system for both ANOPP and VMS is the engine performance. This will also be discussed in the presentation.

Rawls, John, Jr.↗

Information Management Workflow and Tools Enabling Multiscale Modeling Within ICME Paradigm

With the increased emphasis on reducing the cost and time to market of new materials, the need for analytical tools that enable the virtual design and optimization of materials throughout their processing - internal structure - property - performance envelope, along with the capturing and storing of the associated material and model information across its lifecycle, has become critical. This need is also fueled by the demands for higher efficiency in material testing; consistency, quality and traceability of data; product design; engineering analysis; as well as control of access to proprietary or sensitive information. Fortunately, material information management systems and physics-based multiscale modeling methods have kept pace with the growing user demands. Herein, recent efforts to establish workflow for and demonstrate a unique set of web application tools for linking NASA GRC's Integrated Computational Materials Engineering (ICME) Granta MI database schema and NASA GRC's Integrated multiscale Micromechanics Analysis Code (ImMAC) software toolset are presented. The goal is to enable seamless coupling between both test data and simulation data, which is captured and tracked automatically within Granta MI®, with full model pedigree information. These tools, and this type of linkage, are foundational to realizing the full potential of ICME, in which materials processing, microstructure, properties, and performance are coupled to enable application-driven design and optimization of materials and structures.

Materials Engineering↗

Integration of Information Management System, Workflow and Computational Tools Enabling Multiscale Modeling Within an ICME Paradigm

With the increased emphasis on reducing the cost and time to market of new materials, the need for analytical tools that enable the virtual design and optimization of materials throughout their processing - internal structure - property - performance envelope, along with the capturing and storing of the associated material and model information across its lifecycle, has become critical. This need is also fueled by the demands for higher efficiency in material testing; consistency, quality and traceability of data; product design; engineering analysis; as well as control of access to proprietary or sensitive information. Fortunately, material information management systems and physics-based multiscale modeling methods have kept pace with the growing user demands. Herein, recent efforts to develop a set of Python functions that exchange information between NASA GRC's Integrated multiscale Micromechanics Analysis Code (ImMAC) software toolset and its Integrated Computational Materials Engineering (ICME), Granta MI® database schema is presented. The goal is to enable seamless coupling between both test data and simulation data, which is captured and tracked automatically within Granta MI®, with full model pedigree information. These tools, and this type of linkage, are foundational to realizing the full potential of ICME, in which materials processing, microstructure, properties, and performance are coupled to enable application-driven design and optimization of materials and structures.

multiscale modeling; Micromechanics; Computational↗

Droplet-Gas Interactions in Nanoelectrospray Multiphase Flow

Electrospray (ES) is production of a charged droplet plume that transits from an electrically biased supply capillary to a counter electrode under the influence of an electric field. Drag interactions with the surrounding gas decelerate the droplets while transferring the momentum from droplets to accelerate the surrounding gas, resulting in an induced gas flow. This work seeks to characterize the structure of the resulting gas jet and identify the key mechanisms defining the flow structure. The phenomenon of gas jets produced by momentum transfer from nano-electrospray (nES) plumes is explored with schlieren visualization, thermal anemometry, and numerical simulations. Schlieren visualization experiments provide information on the flow structure in support of simulation predictions, and the hot thermistor anemometry measurements of gas velocities outside the spray demonstrate quantitatively validated simulation results. The study reveals the formation of a moderately high velocity coaxial gas jet within the nES plume and provides insight into the overall flow structure of the induced flow. The multiphase electrohydrodynamic simulations enable numerical experimentation to explore the fundamental physics of coupled droplet-gas transport and the resulting flow structure. The simulations, confirmed by the experiments, reveal gas jetting induced by nES with a narrow (sub hundred micrometers in diameter) core originating from the nES liquid-jet breakup region and a surrounding larger-in-extent zone (several hundred micrometers in diameter) of lower velocity gas flow within the electrospray plume. This behavior is due to nES ejecting a stream of droplets from the tip of a narrow liquid cone-jet, where the combined effect of many small droplets transferring momentum to a confined region of gas yields a narrow but high-velocity gas stream. Furthermore, the micro-jets produced by nES have practical utility for mass spectrometry, 3D printing and fabrication, and thermal management.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

GRid Analysis and Visualization Interface (GRAVI) [SWR-24-16]

GRAVI (GRid Analysis and Visualization Interface) is a web application for viewing and analyzing nodal Production Cost Model (PCM) and Capacity Expansion Model (CEM) simulations. The web application provides the ability to animate geospatially coupled timeseries data in an agnostic way regardless of the underlying simulation tool used to generate the data. GRAVI also provides capabilities to animate non-geospatial data relevant to a PCM or CEM model. Furthermore, this web application can be tailored as an real-time operational tool to better understand a live grid.

Webb, Micah↗

Increasing CO2 Coupled with Other Anthropogenic Perturbations: Effects on Ozone and Other Trace Gases

The GSFC 2D interactive chemistry-radiation-dynamics model has been used to study the effects on stratospheric trace gases of past and future CO2 increases coupled with changes in CFC'S, methane, and nitrous oxide. Previous simulations with the GSFC model showed that the stratospheric cooling calculated to result from doubling atmospheric CO2 would lead, in the absence of a growth of other anthropogenic gases, to a decrease in upper stratospheric NO(y) of roughly 15%. This work has been extended to simulate changes in stratospheric chemistry and dynamics occurring between the years 1960 and 2050. The simulations have been carried out with and without changes in CO2. In the low latitude upper stratosphere ozone is predicted to be 10% greater in 2050 than in 1990 when increased CO2 is included, compared with an increase of only 2% without the inclusion of CO2. In the low latitude lower stratosphere, ozone is predicted to decrease by about 1% between 1990 and 2050 when CO2 changes are taken into account, in contrast to an approximate 3% increase when they are not. The simulated behavior of water vapor is another example of the coupled responses. Between 1990 and 2050 low latitude water vapor is predicted to increase by 4% and 2% in the upper and lower stratosphere, respectively, without the inclusion of CO2 increases. with the inclusion of CO2 changes, the water vapor increases are predicted to be roughly 12% and 8%, for the upper and lower stratosphere, respectively.

Rosenfield, J. E.↗

New Insights in Tropospheric Ozone and its Variability

We have produced time-slice simulations using the Goddard Earth Observing System Version 5 (GEOS-5) coupled to a comprehensive stratospheric and tropospheric chemical mechanism. These simulations are forced with observed sea surface temperatures over the past 25 years and use constant specified surface emissions, thereby providing a measure of the dynamically controlled ozone response. We examine the model performance in simulating tropospheric ozone and its variability. Here we show targeted comparisons results from our simulations with a multi-decadal tropical tropospheric column ozone dataset obtained from satellite observations of total column ozone. We use SHADOZ ozonesondes to gain insight into the observed vertical response and compare with the simulated vertical structure. This work includes but is not limited to ENSO related variability.

Oman, Luke D.↗

LOKI

LOKI is a simulation code that models weakly-coupled plasma evolution using Eulerian kinetic models in four phase space dimensions. The code is useful to understand basic plasma nonlinear plasma physics phenomena - in particular, wave-particle and wave-wave resonances that occur in laser plasma interactions. LOKI can simulate fundamental physical process in low-density plasmas on a scale of hundreds of light wavelengths. It does so by solving the relativistic Vlasov and Valsov-Fokker-Planck differential transport equations describing the collision less and low-collisionality evolution of charged particles in electromagnetic fields in two configuration space and two velocity space dimensions. Two principle electromagnetic models are provided: electrostatic (Poisson's Equation) and electromagnetic (Maxwell's equations).

DIMITS, ANDRIS↗

3D Numerical Simulation of Kink-driven Rayleigh–Taylor Instability Leading to Fast Magnetic Reconnection

Fast magnetic reconnection involving non-MHD microscale physics is thought to underlie both solar eruptions and laboratory plasma current disruptions. While there is extensive research on both the MHD macroscale physics and the non-MHD microscale physics, the process by which large-scale MHD couples to the microscale physics is not well understood. An MHD instability cascade from a kink to a secondary Rayleigh–Taylor instability in the Caltech astrophysical jet laboratory experiment provides new insights into this coupling and motivates a 3D numerical simulation of this transition from large to small scale. A critical finding from the simulation is that the axial magnetic field inside the current-carrying dense plasma must exceed the field outside. In addition, the simulation verifies a theoretical prediction and experimental observation that, depending on the strength of the effective gravity produced by the primary kink instability, the secondary instability can be Rayleigh–Taylor or mini-kink. Finally, it is shown that the kink-driven Rayleigh–Taylor instability generates a localized electric field sufficiently strong to accelerate electrons to very high energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data accompanying the publication "Sensitivity of grass fires burning in marginal conditions to atmospheric turbulence" in the Journal of Geophysical Research-Atmospheres

Data accompanying Jonko, Yedinak, Conley, and Linn (2021): "Sensitivity of grass fires burning in marginal conditions to atmospheric turbulence", JGR-Atmospheres. LA-UR-21-24832. The following files were generated from simulations performed with Los Alamos National Laboratory coupled atmosphere-fire behavior model FIRETEC v1.3: -area_burned.csv contains the area burned per second in m^2 for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -fire_perimeter.csv containing the fire perimeters in m for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -consumption_per_second.csv containing consumption in kg/s for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -rate_of_spread.csv containing the distance of the fire front from the ignition line in m for 45 CORE simulations, 20 INIT simulations and 20 BOUND simulations (85 columns), and 600 seconds of simulation (600 rows). -CORE1_u_10m_300_sec.hdf, INIT3_u_10m_300_sec.hdf, and BOUND3_u_10m_300_sec.hdf contain u wind speeds at 10 m height for 300 seconds (3D: time, x, y) of simulation for simulations CORE1, INIT3, and BOUND3, respectively. -CORE.final.fuel.density.hdf contains final fuel densities in the surface layer of the model (3D: ensemble member, x, y) for all 45 CORE simulations. -CORE{1-45}_3D_u_120_sec.hdf: 45 files containing vertically and horizontally resolved u winds for 120 seconds of simulation (4D: time, x, y, z) for the 45 CORE ensemble members. For additional information about this data, please contact the corresponding author at ajonko@lanl.gov.

54 ENVIRONMENTAL SCIENCES↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

A comparative study of experimental and simulated ultrasonic pulse-echo signals from multilayered structures

The effect on the system acoustic response of variations in the adhesive thickness, coupling thickness, and paint thickness is considered. Both simulations and experimental measurements are used to characterize and classify A-scans from test regions, and to study the effects of various parameters such as paint thickness and epoxy thickness on the variations in the reflected signals. A 1D model of sound propagation in multilayered structures is used to verify the validity of the measured signals, and is also used to computationally generate signals for a class of test locations with gradually varying parameters. This approach exploits the ability of numerical simulations to provide a good understanding of the ultrasonic pulses reflected at disbonds.

Abedin, M. N.↗

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗