Engineering Papers⌕ Search

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.

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

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Providing Ramping Service with Wind to Enhance Power System Operational Flexibility

This project develops an innovative, integrated, and transformative approach to mitigate the impact of wind ramping by providing flexible ramping products from wind power. The project will significantly contribute to the reduction of wind integration costs by making wind power dispatchable and allowing the efficient management of wind ramping characteristics. In addition, we collaborated with industry to test and validate the methodology, considering economic and reliability goals, by integrating the proposed methodology into the simulated market frameworks of two independent system operators: California Independent System Operator and the Electric Reliability Council of Texas.

17 WIND ENERGY↗
DOE OSTI2020

A Bezier Curve Informed Melt Pool Geometry to Model Additive Manufacturing Microstructures Using SPPARKS

Additive manufacturing is a transformative technology with the potential to manufacture designs which traditional subtractive machining methods cannot. Additive manufacturing offers fast builds at near final desired geometry; however, material properties and variability from part to part remain a challenge for certification and qualification of metallic components. AM induced metallic microstructures are spatially heterogeneous and highly process dependent. Engineering properties such as strength and toughness are significantly affected by microstructure morphologies resulting from the manufacturing process Linking process parameters to microstructures and ultimately to the dynamic response of AM materials is critical to certifying and qualifying AM built parts and components and improving the performance of AM materials. The AM fabrication process is characterized by building parts layer by layer using a selective laser melt process guided by a computer. A laser selectively scans and melts metal according to a designated geometry. As the laser scans, metal melts, fuses, and solidifies forming the final geometry in a layerwise fashion. As the laser heat source moves away, the metal cools and solidifies forming metallic microstructures. This work describes a microstructure modeling application implemented in the SPPARKS kinetic Monte Carlo computational framework for simulating the resulting microstructures. The application uses Bzier curves and surfaces to model the melt pool surface and spatial temperature profile induced by moving the laser heat source; it simulates the melting and fusing of metal at the laser hot spot and microstructure formation and evolution when the laser moves away. The geometry of the melt pool is quite flexible and we explore effects of variances in model parameters on simulated microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Selecting and Implementing Resilience Metrics in Existing Energy Sector Models [Slides]

Resilience is a topic receiving much attention in relation to energy systems, with particular attention being paid to the supply of electricity. As a result of the growing interest in energy sector resilience, research communities have proposed a plethora of candidate resilience indicators and metrics, most of which remain immature at different scales and segments within the energy system. A necessary focus of the research community lies in implementing, testing, and validating resilience metrics and analysis approaches in energy sector models, which will be invaluable for informing resilience planning and investment decisions. Recognizing these challenges that need to be addressed, we explore how to effectively integrate resilience considerations into energy sector models and tools. The overarching goal of the effort was to evaluate the data needs, methodologies, and outcomes - including consequences and/or changes in investment or operational decisions due to avoided consequences - based on resilience analysis in a range of existing tools. In particular, we selected five models originally built at NREL to explore non-resilience energy research questions to implement and exercise resilience metrics and analysis approaches. To demonstrate the importance of perspective, we selected models that represent different segments of the energy sector, geographic scales, and modeling approaches. A second important aspect of our effort was the development of generalized power interruption scenarios. These scenarios were intended to help establish a framework for simulating the effects of real-world threats in terms of their impacts on system components and, in turn, power interruption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Demonstration of Advanced Experimental and Theoretical Characterization of Hydrogen Dynamics and Associated Behavior in Advanced Reactors

Advanced materials development, manufacturing, and modeling capabilities for innovative reactor designs support nuclear security and mission-focused science through enhanced technology for safer and more efficient and secure production of nuclear energy. The research in this project has established: 1) a state-of-the-art neutron-based hydrogen mapping and cross-section measurement capability as well as detailed crystallographic characterization of hydrogen atoms at LANSCE, and 2) a multi-physics framework for simulating behavior of moderator materials and other material performance in advanced nuclear reactors. Through the course of this project, we successfully developed and demonstrated measurement techniques for hydrogen distribution and atomistic-scale behavior of hydrogen atoms using pulsed neutron techniques. In parallel, advanced multi-physics simulation tools to predict the behavior of hydrogen atoms, e.g. in a moderator for a nuclear reactor, through materials performance, neutron transport, and thermal mechanical behavior were enhanced. Multi-discipline areas across the laboratory were involved in the project as the integration of improved experimental capabilities with enhanced modeling and simulation through MST, NEN, SIGMA, and XCP division subject matter experts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Bioreactor Optimization through Multi-Phase Flow Models (CRADA Final Report)

Chemical manufacturing uses 29% of energy in the United States and produces 925 million metric tons of CO2 annually. Biomanufacturing offers the potential to leverage America’s rich agricultural resources to produce critical chemicals such as lubricants, pharmaceutical precursors, and components of energetic materials that today are sourced extensively from overseas. The bioreactors used in biomanufacturing applications, such as one developed by Capra Biosciences, involve multiphase flow of biofilm-coated solid support particles that are continuously circulated in a fluidized state within the reactor along with a constant supply of oxygen via an aeration mechanism. In this project, Capra Biosciences and LBNL developed a multiscale modeling framework to simulate the multiphase flows of solid particles in a liquid-gas bubble mixture that occurs in the bioreactor using the current MFIX-Exa software, an opensource multiphase flow solver developed and maintained at LBNL and NETL. By leveraging HPC capabilities, this high-fidelity multiscale model was used to inform design decisions for bioreactor architecture to make them operationally efficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
DOE OSTI2025

The Monolithic Heat Pipe Microreactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the generic monolithic heat-pipe-cooled microreactor reference plant model for the United States Nuclear Regulatory Commission. This model serves as the foundation for the future detailed design evaluation models based on license applications. This model has been developed with the Comprehensive Reactor Analysis Bundle (BlueCRAB) and its specifications are based on open literature publications for the eVinci™ design . BlueCRAB is the U.S. Nuclear Regulatory Commission non-light-water reactor analysis system based on MOOSE, the Multiphysics Object-Oriented Simulation Environment framework, which can couple the Griffin, BISON, and Sockeye applications to resolve the various physics that are essential for the safety analysis of this type of reactor system. The core specifications include tristructural isotropic fuel, graphite monolith, graphite reflectors, and drums composed of graphite and B 4 C.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Stochastic microstructural analysis of ceramic matrix composites using a high-fidelity multiscale framework

This paper presents a high-fidelity multiscale framework to simulate the mechanical behavior of ceramic matrix composites (CMCs), accounting for the complex micromorphology captured using detailed material characterization. First, high-resolution micrographs are obtained for the specific carbon fiber silicon-carbide-nitride matrix (C/SiNC) CMC to characterize the variability of the architectural features and manufacturing-induced defects at the microscale. An image processing algorithm is then used to precisely estimate the size and distribution of all subscale features and defects from the micrographs. The information is then used to generate a three-dimensional stochastic representative volume element (SRVE) to reconstruct microscale constituents accounting for the variability. As a result, the generated SRVEs are simulated using the high-fidelity generalized method of cells (HFGMC) micromechanics theory to investigate the effects of defects on the elastic properties of C/SiNC CMCs.

Ceramic Matrix Composites↗
DOE OSTI2015

Species interactions differ in their genetic robustness

Conflict and cooperation between bacterial species drive the composition and function of microbial communities. Stability of these emergent properties will be influenced by the degree to which species' interactions are robust to genetic perturbations. We use genome-scale metabolic modeling to computationally analyze the impact of genetic changes when Escherichia coli and Salmonella enterica compete, or cooperate. We systematically knocked out in silico each reaction in the metabolic network of E. coli to construct all 2583 mutant stoichiometric models. Then, using a recently developed multi-scale computational framework, we simulated the growth of each mutant E. coli in the presence of S. enterica . The type of interaction between species was set by modulating the initial metabolites present in the environment. We found that the community was most robust to genetic perturbations when the organisms were cooperating. Species ratios were more stable in the cooperative community, and community biomass had equal variance in the two contexts. Additionally, the number of mutations that have a substantial effect is lower when the species cooperate than when they are competing. In contrast, when mutations were added to the S. enterica network the system was more robust when the bacteria were competing. These results highlight the utility of connecting metabolic mechanisms and studies of ecological stability. Cooperation and conflict alter the connection between genetic changes and properties that emerge at higher levels of biological organization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Housing-Performance Atlas of Baltimore Row Homes: Archetype-Based Multi-Hazard Baseline of Energy, Heat, Survivability, and Durability

Baltimore’s historic row-home neighborhoods face escalating risks to energy, heat, and durability under intensifying climate stress. This study develops a Housing-Performance Atlas that quantifies multi-hazard performance for eight representative archetypes using DesignBuilder/EnergyPlus Version 7.3.1.003, under Baltimore TMY3 boundary conditions. Performance is evaluated across the following four adaptation domains: energy use intensity, passive survivability during 72 h outage events, roof overheating exposure (>150 °F exceedance hours), and material service life derived from ISO 15686 and synthesized into Lean and Full Deficit Indices for comparative resilience ranking. Results show that EUI ranged from 46.7 to 67.6 kBtu ft −2 ·yr −1 , survivability from 0 to 23 h, and roof temperatures exceeded 150 °F for 150–210 h, shortening roof service life by up to 10 years. Composite Lean and Full Deficit Indices ranged 7.8–92.4, ranking Model 5 (end-unit, flat roof, two-story with basement) as the most resilient configuration and Model 8 (end-unit, pitched roof, three-story above-grade) as the least resilient due to compounded overheating and energy losses. Heat-related domains accounted for nearly 70% of overall resilience deficits, confirming thermal safety and roof reflectivity as retrofit priorities. The Housing-Performance Atlas establishes a reproducible diagnostic framework linking simulation, service life, and resilience metrics to guide cost-effective, climate-responsive retrofits in Baltimore’s aging urban housing stock.

Housing-Performance atlas↗
DOE OSTI2016

Cloud structure of galactic OB cluster-forming regions from combining ground- and space-based bolometric observations

We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground-based (CSO, JCMT, APEX, and IRAM-30 m) and space telescopes (Herschel and Planck). For the seven luminous (L>10{sup 6} L {sub ⊙}) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, and G10.3-0.1, we have performed single-component, modified blackbody fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The ∼10″ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3 kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed ∼10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function (N-PDF), the two-point correlation function of column density, and the probability distribution function of clump/core separations). In particular, the massive molecular gas clumps located at the center of G10.6-0.4 and W49A, which contribute to a considerable fraction of their overall cloud masses, may be special OB cluster-forming environments as a direct consequence of global cloud collapse. These centralized massive molecular gas clumps also uniquely occupy much higher column densities than what is determined by the overall fit of power-law N-PDF. We have made efforts to archive the derived statistical quantities of individual target sources, to permit comparisons with theoretical frameworks, numerical simulations, and other observations in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Observing Exoplanets with High Dispersion Coronagraphy. I. The Scientific Potential of Current and Next-generation Large Ground and Space Telescopes

Direct imaging of exoplanets presents a formidable technical challenge owing to the small angular separation and high contrast between exoplanets and their host stars. High Dispersion Coronagraphy (HDC) is a pathway to achieve unprecedented sensitivity to Earth-like planets in the habitable zone. Here, we present a framework to simulate HDC observations and data analyses. The goal of these simulations is to perform a detailed analysis of the trade-off between raw star light suppression and spectral resolution for various instrument configurations, target types, and science cases. We predict the performance of an HDC instrument at Keck observatory for characterizing directly imaged gas-giant planets in near-infrared bands. We also simulate HDC observations of an Earth-like planet using next-generation ground-based (TMT) and spaced-base telescopes (HabEx and LUVOIR). We conclude that ground-based ELTs are more suitable for HDC observations of an Earth-like planet than future space-based missions owing to the considerable difference in collecting area. For ground-based telescopes, HDC observations can detect an Earth-like planet in the habitable zone around an M-dwarf star at 10{sup −4} star light suppression level. Compared to the 10{sup −7} planet/star contrast, HDC relaxes the star light suppression requirement by a factor of 10{sup 3}. For space-based telescopes, detector noise will be a major limitation at spectral resolutions higher than 10{sup 4}. Considering detector noise and speckle chromatic noise, R = 400 (1600) is the optimal spectral resolutions for HabEx (LUVOIR). The corresponding star light suppression requirement to detect a planet with planet/star contrast = 6.1×10{sup −11} is relaxed by a factor of 10 (100) for HabEx (LUVOIR).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bioverse: A Comprehensive Assessment of the Capabilities of Extremely Large Telescopes to Probe Earth-like O 2 Levels in Nearby Transiting Habitable-zone Exoplanets

Molecular oxygen is a strong indicator of life on Earth and may indicate biological processes on exoplanets too. Recent studies proposed that Earth-like O 2 levels might be detectable on nearby exoplanets using high-resolution spectrographs on future extremely large telescopes (ELTs). However, these studies did not consider constraints like relative velocities, planet occurrence rates, and target observability. We expanded on past studies by creating a homogeneous catalog of 286,391 main-sequence stars within 120 pc using Gaia DR3 and used the Bioverse framework to simulate the likelihood of finding nearby transiting Earth analogs. We also simulated a survey of M dwarfs within 20 pc accounting for η ⊕ estimates, transit probabilities, relative velocities, and target observability to determine how long ELTs and theoretical 50–100 m ground-based telescopes need to observe to probe for Earth-like O 2 levels with an R = 100,000 spectrograph. This would only be possible within 50 yr for up to ~21% of nearby M-dwarf systems if a suitable transiting habitable-zone Earth analog was discovered, assuming signals from every observable partial transit from each ELT can be combined. If so, Earth-like O 2 levels could be detectable on TRAPPIST-1 d–g within 16–55 yr, respectively, and about half that time with an R = 500,000 spectrograph. These results have important implications for whether ELTs can survey nearby habitable-zone Earth analogs for O 2 via transmission spectroscopy. Our work provides the most comprehensive assessment to date of the ground-based capabilities to search for life beyond the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Planetary Accretion Shock. I. Framework for Radiation-hydrodynamical Simulations and First Results

The key aspect determining the postformation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the surface of the planet. We use one-dimensional radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain postshock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux that escapes the system (roughly the Hill sphere), taking into account the energy recycling that occurs ahead of the shock in a radiative precursor. We focus in this paper on a constant equation of state (EOS) to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will have to await a self-consistent treatment including hydrogen dissociation and ionization, the results presented here show the correct qualitative behavior and can be understood from semianalytical calculations. The shock is found to be isothermal and supercritical for a range of conditions relevant to the core accretion formation scenario (CA), with Mach numbers M≳3. Across the shock, the entropy decreases significantly by a few times k{sub B}/baryon. While nearly 100% of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40%, implying that a significant fraction of the total accretion energy is brought into the planet. However, for realistic parameter combinations in the CA scenario, we find that a nonzero fraction of the luminosity always escapes the Hill sphere. This luminosity could explain, at least in part, recent observations in the young LkCa 15 and HD 100546 systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Dust Polarization Maps from TIGRESS: E/B Power Asymmetry and TE Correlation

We present the first large set of all-sky synthetic dust polarization maps derived directly from a self-consistent magnetohydrodynamics simulation using the TIGRESS framework. Turbulence in this simulation is predominantly driven by supernova explosions, with rates that are self-consistently regulated by feedback loops. The simulation covers both the outer scale and inertial range of turbulence with uniformly high resolution. The shearing-box utilized in the simulation, in concert with resolved supernova-driven turbulence, enables the capturing of generation, growth, and saturation of both turbulent and mean magnetic fields. We construct polarization maps at 353 GHz, as seen by observers inside a model of the multiphase, turbulent, magnetized interstellar medium (ISM). To fully sample the simulated ISM state, we use 350 snapshots spanning over $\sim 350\,\mathrm{Myr}$ (more than six feedback loops) and nine representative observers. The synthetic skies show a prevalent E/B power asymmetry (${EE}\gt {BB}$) and positive TE correlation in broad agreement with observations by the Planck satellite. However, the ranges of ${EE}/{BB}\sim 1.4\mbox{-}1.7$ and ${TE}/{({TT}\cdot {EE})}^{1/2}\sim 0.2\mbox{-}0.3$ are generally lower than those measured by Planck. We find large fluctuations of E/B asymmetry and TE correlation depending on the observer's position and temporal fluctuations of ISM properties due to bursts of star formation. The synthetic maps are made publicly available to provide novel models of the microwave sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Multigroup Thermal Radiation Transport with Tensor Trains

We investigate the application of tensor-train (TT) algorithms to multigroup thermal radiation transport (i.e., photon radiation transport). The TT framework enables simulations at discretizations that might otherwise be computationally infeasible on conventional hardware. We show that solutions to certain multigroup problems possess an intrinsic low-rank structure, which the TT representation leverages effectively. This enables us to solve problems where the discretized solution size exceeds a trillion parameters on a single node. The solver is evaluated on a range of test problems with varying levels of complexity, consistently achieving compression factors greater than 100× and speedups exceeding 2×. We also investigate alternative TT topologies by analyzing the low-rank structure of the merged spatio-spectral core to assess the potential for greater compression. This analysis suggests that compression gains could increase by factors as large as 7. Our results indicate that the low-rank structure of the merged spatio-spectral core captures the spatio-spectral complexity of the solution, largely driven by the opacity structure of the medium. Beyond identifying opportunities for improved compression, this analysis highlights the types of errors that may arise in angle-integrated quantities when exploiting this low-rank structure.

79 ASTRONOMY AND ASTROPHYSICS↗