SEARCH · Engineering Papers

Results for “simulator study”

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 1,333 records · Page 74

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Subsurface Characterization of Hydraulic Fracture Test Site-2 (HFTS-2), Delaware Basin

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp Formation of the Permian (Delaware) Basin. This paper focuses on integration, advanced geological characterization, and 3D subsurface modeling of the comprehensive HFTS-2 dataset. The study showcases a multidisciplinary reservoir characterization approach that incorporates geology, petrophysics, geochemistry, geomechanics, microseismic, and subsurface engineering analysis. Subsurface characterization of organic-rich mudstone formations requires understanding complex hydraulic fracture network growth in relation to inherent lithology, geomechanical properties, and interaction with pre-existing natural fractures. This paper presents a characterization workflow incorporating pre- and post-stimulation subsurface data, unique to the HFTS-2 dataset. The study integrated: (1) rock properties from logs, cores, and thin sections; (2) natural and hydraulic fracture descriptions from cores and image logs; (3) local and regional stresses; (4) geomechanics; (5) microseismic; (6) fiber optic (FO) and bottomhole pressure gauge (BHPG) response; and (7) produced fluids analysis. During a stimulation treatment, creation of the stimulated rock volume (SRV) is influenced by several subsurface factors. Key contributing factors include structural context, stress conditions, lithology, facies architecture, pre-existing natural fractures, and geomechanical properties. The HFTS-2 subsurface data integration indicates that the SRV is comprised of a complex juxtaposition of hydraulic fracture swarms, as evidenced by image logs analysis, core description, and microseismic monitoring. The HFTS-2 microseismic event density was used to generate 3D heat maps that serve as a representative SRV footprint, corroborated by secondary datasets. These maps were further integrated with petrophysical and geomechanical characteristics, as well as responses from FO and BHPG, to estimate the lateral and vertical dimensions of the effective fractures. The geological characterization for the HFTS-2 dataset combined with 3D modeling for petrophysical and geomechanical properties provides a strong foundation for subsurface simulation and optimization studies. Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S005R001/2477501/urtec-2021-5243-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5243 The workflow improved our understanding of HFTS-2 hydraulic fracture propagation and characteristics in relation to offset pressure depletion and interaction with pre-existing natural fractures. Analysis showed that fracture geometry varies by stage and by well, and a complex fracture network is generated with varying fracture density. The multidisciplinary workflow presented herein for integration and characterization serves as a foundation to evaluate completion efficiency and estimate areal and vertical stimulation and depletion extent for the project. Furthermore, the workflow and learnings can also be transferred to other unconventional plays.

58 GEOSCIENCES↗

DOE OSTI2020

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: The TRAPPIST-1 Case Study

Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in X-ray and ultraviolet (XUV) irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on an ultraviolet continuum flux measurement (centered around ∼1900 Å) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photodissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 × 10{sup −4} L {sub *}, while the full flux including non-flaring periods is 6.1 × 10{sup −5} L {sub *}, when L {sub *} is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass loss and tidal evolution and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio Q{sup ′}=3Q/2k{sub 2} for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator VPLanet to study ocean retention for these planets using our derived system parameters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Effect of Production Bias on Radiation-Induced Segregation in Ni-Cr Alloys

We present an in-depth investigation into the Radiation-Induced Segregation (RIS) phenomenon in Ni-Cr alloys. All the pivotal factors affecting RIS such as surface’s absorption efficiency, grain size, production bias, dose rate, temperature, and sink density were systematically studied. Through comprehensive simulations, the individual and collective impacts of these factors were analyzed, enabling a refined understanding of RIS. A notable finding was the significant influence of production bias on point defects’ interactions with grain boundaries/surfaces, thereby playing a crucial role in RIS processes. Production bias alters the neutrality of these interactions, leading to a preferential absorption of one type of point defect by the boundary and consequent establishment of distinct surface-mediated patterns of point defects. These spatial patterns further result in non-monotonic spatial profiles of solute atoms near surfaces/grain boundaries, corroborated by experimental observations. In particular, a positive production bias, signifying a higher production rate of vacancies over interstitials, drives more Cr depletion at the grain boundary. Moreover, a temperature-dependent production bias must be considered to recover the experimentally reported dependence of RIS on temperature. The severity of radiation damage and RIS becomes more pronounced with increased production bias, dose rate, and grain size, while high temperatures or sink density suppress the RIS severity. Model predictions were validated against experimental data, showcasing robust qualitative and quantitative agreements. The findings pave the way for further exploration of these spatial dependencies in subsequent studies, aiming to augment the comprehension and predictability of RIS processes in alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2021

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Response to Community-Scale Geothermal Energy Production

Geothermal heat, as renewable energy, shows great advantage with respect to its environmental impact due to its significantly lower CO2 emissions than conventional fossil fuel. Open and closed-loop geothermal heat pumps, which utilize shallow geothermal systems, are an efficient technology for cooling and heating buildings, especially in urban areas. Integrated use of geothermal energy technologies for district heating, cooling, and thermal energy storage can be applied to optimize the subsurface for communities to provide them with multiple sustainable energy and community resilience benefits. The utilization of the subsurface resources may lead to a variation in the underground environment, which might further impact local environmental conditions. However, very few simulators can handle such a highly complex set of coupled computations on a regional or city scale. We have developed high-performance computing (HPC) based hydrothermal finite element (FE) simulator that can simulate the subsurface and its hydrothermal conditions at a scale of tens of km. The HPC simulator enables us to investigate the subsurface thermal and hydrologic response to the built underground environment (such as basements and subways) at the community scale. In this study, a coupled hydrothermal simulator is developed based on the open-source finite element library deal.II. The HPC simulator was validated by comparing the results of a benchmark case study against COMSOL Multiphysics, in which Aquifer Thermal Energy Storage (ATES) is modeled and a process of heat injection into ATES is simulated. The use of an energy pile system at the Treasure Island redevelopment site (San Francisco, CA, USA) was selected as a case study to demonstrate the HPC capability of the developed simulator. The simulator is capable of modeling multiple city-scale geothermal scenarios in a reasonable amount of time.

Sun, Xiang↗

DOE OSTI2025

Kinetics and Thermodynamics of Sr Permeation in CeO 2 -Based Barrier Layers for Solid-Oxide Electrolyzer Cells

Solid-oxide electrolyzer cells (SOECs) convert steam to hydrogen efficiently at high temperatures. However, during operation, the diffusion of cations or impurities through the cells due to electrode degradation can cause unwanted secondary phases to form, which may degrade device performance. Here, in this study, we use atomistic and mesoscale simulations coupled with experimental analysis to study the diffusion of Sr through the Gd-doped CeO 2 (GDC) barrier layer used to protect the yttria-stabilized zirconia (YSZ) electrolyte in SOECs. From our atomistic calculations, we find Sr diffusion to be negligibly slow in bulk GDC; however, surface diffusion is much more favorable. Subsequent mesoscale simulations show that Sr diffusion is activated when the porosity of GDC exceeds ∼10% and significantly exceeds diffusion in bulk and grain boundary regions. We also find that SrO-based species can accumulate at GDC surfaces; however, SrO aggregation and coarsening will be limited by the large lattice mismatch between GDC and SrO. Energy-dispersive X-ray spectroscopy (EDS) and electron diffraction confirm that Sr can accumulate within GDC pores and form disperse Sr-containing secondary phases. Altogether, Sr diffusion in dense GDC is unlikely to give rise to thick SrO layers, which would severely limit device performance. The formation of Sr-containing secondary phases can largely be avoided by restricting the porosity of the GDC layer as much as possible.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Revisiting W–ZrC interfaces: A first principles study

We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Injection simulations of space charge dominated proton beams in IOTA

A 2.5 MeV proton injector is being constructed for the IOTA ring at Fermilab to study the interaction of nonlinear integrable optics (NIO) with high space charge beams. Space charge in the transport line from the RFQ to the injection location has a significant current dependent effect on the phase space. Simulation studies to support efficient injection of intense bunches into IOTA are presented, included schemes to inject directly into NIO lattices.

Accelerator Physics↗

DOE OSTI2025

Injection simulations of space charge dominated proton beams in IOTA

Wieland, J. [Fermilab; Michigan State U.] (ORCID:0↗

DOE OSTI2020

The passivity of low-temperature carburized austenitic stainless steel AISI-316L in a simulated boiling-water-reactor environment

The passive film on low-temperature-carburized AISI-316L austenitic stainless steel was studied in an oxygenated, simulated boiling-water-reactor environment with normal water chemistry. The microstructure was studied using scanning electron microscopy, grazing-incidence X-ray diffractometry, Auger electron spectroscopy, X-ray energy-dispersive spectroscopy, and transmission electron microscopy. For non-surface engineered AISI-316L, three distinct layers were observed after exposure: an outer layer of larger nickel-oxide enriched particles with the spinel structure, an intermediate layer of small Cr-oxide-enriched hematite particles, and a compact inner layer of Cr-enriched oxide also having the spinel structure. Exposure of low-temperature-carburized AISI-316L, in contrast, resulted in only two distinct layers. These layers did not include a compact Cr-rich inner layer. Instead, the innermost layer consisted of an oxide with the spinel structure that was a mixture of Fe, Cr and Ni oxide. The outer layer on the low-temperature-carburized specimen consisted of large loosely stacked particles with the spinel structure of mixed composition. In addition to passive-film characterization, the corrosion rate was studied by monitoring mass loss. Finally, after 500 h of immersion, the mass changes measured were relatively low for all specimens indicating good adhesion and low oxidation rates.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Systematics of the chemical freeze-out line in the high baryon density regime explored at SIS100

The systematic uncertainties of chemical freeze-out fits at SIS100 energies (Au+Au reactions at $\sqrt{s_{NN}}$ = 3 – 5 GeV) are studied using UrQMD simulations. Although hadron production in UrQMD does not occur on a sharp chemical freeze-out hyper-surface, the extracted fit quality is shown to be very good. The extracted chemical parameters depend on the selected hadron species as well as the underlying equation of state (EoS) of the matter. Including light nuclei and anti-protons in the fit increases the expected freeze-out temperature, while a stiffer EoS increases the obtained chemical potential. Similarly, the baryon densities extracted by the thermal fits depend on the choice of hadrons as well as the underlying equation of state. These results are important for the upcoming CBM@FAIR physics program and highlight that a degree of caution is advised when one relates the chemical freeze-out curve to features on the QCD phase diagram like the critical endpoint or a possible phase transition.

Hofmann, Emma Lilith [GSI Helmholtzzentrum für Sch↗

DOE OSTI2022

Performance Evaluation of the Meteorology and Air Quality Conditions From Multiscale WRF-CMAQ Simulations for the Long Island Sound Tropospheric Ozone Study (LISTOS)

The Long Island Sound (LIS) Tropospheric Ozone Study was a multi-agency collaborative field campaign conducted during the summer of 2018 to improve the understanding of ozone chemistry and transport from New York City to areas downstream, especially the LIS and adjacent Connecticut coastline. Measurements made during this campaign were leveraged to test and evaluate the coupled WRF-CMAQ model at 12 km, 4 and 1.33 km horizontal grid spacing. Special attention was placed on the model's representation of sea breeze circulations, low level jets, and boundary layer evolution. The evaluation suggests using higher resolutions resulted in improved surface meteorology statistics throughout the whole summer, with temperature biases seeing the biggest statistical improvements when using 1.33-km grid spacing, going from –0.12 to 0.08 K. Additionally, 4-km grid spacing provided the biggest advantage when simulating ozone over the region of interest, with biases being reduced from 2.40 to 0.57 to 0.37 ppbV with increased resolution. Case studies of two high ozone concentration events (July 10 and August 6) revealed that sound breezes and low-level jets had a critical role in transporting pollutant-rich, shallow marine air masses from the LIS inland over the Connecticut coast. Modifications were made to the representation of sea surface temperatures, which subsequently improved the simulation of surface ozone predictions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

CFD Simulation of Floating Body Motion with Mooring Dynamics: Coupling MoorDyn with OpenFOAM

It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.

floating body↗

DOE OSTI2021

Interpreting internal consistency of DES measurements

ABSTRACT Bayesian evidence ratios are widely used to quantify the statistical consistency between different experiments. However, since the evidence ratio is prior dependent, the precise translation between its value and the degree of concordance/discordance requires additional information. The most commonly adopted metric, the Jeffreys scale, can falsely suggest agreement between data sets when priors are chosen to be sufficiently wide. This work examines evidence ratios in a DES-Y1 simulated analysis, focusing on the internal consistency between weak lensing and galaxy clustering. We study two scenarios using simulated data in controlled experiments. First, we calibrate the expected evidence ratio distribution given noise realizations around the best-fitting DES-Y1 ΛCDM cosmology. Secondly, we show the behaviour of evidence ratios for noiseless fiducial data vectors simulated using a modified gravity model, which generates internal tension in the LCDM analysis. We find that the evidence ratio of noise realizations generated at all confidence levels was biased towards agreement and show, with a modified gravity model, that the choice of prior could conceal the discrepancies between weak lensing and galaxy clustering induced by prior effects in unlike cosmological models, concluding that the evidence ratio in a DES-Y1 study is, indeed, biased towards agreement. Boundary effects can also influence conclusions about the inconsistency induced by modified gravity, even in a noiseless data vector simulation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗

DOE OSTI2021

Enhancing the Operational Resilience of Advanced Reactors with Digital Twins by Recurrent Neural Networks

Because of a lack of operational data and uncertainty in evaluation model for abnormal and accident scenarios, the established operating procedures can be biased in characterizing the reactor states and ensuring operational resilience. To reduce uncertainty associated with actual plant conditions, digital twin (DT) technology is suggested to support operator’s decision-making by effectively extracting and using knowledge of the current and future plant states from the knowledge base. This study first builds a knowledge base based on the characterization of issue space and the simulation tool. Next, this study discusses diagnosis and prognosis DTs for enhancing operational resilience by recovering the complete states of reactors and by predicting the future reactor behaviors. Finally, the decision-making module of the control system can determine the optimal control strategy that meets operational goals during loss-of-flow scenarios. To demonstrate and evaluate the DTs capability for supporting the operations of nuclear reactors, this study develops and assesses both the diagnosis and prognosis DTs in a nearly autonomous management and control system for an Experimental Breeder Reactor-II simulator during different loss-of-flow scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Data driven methods to recognize patterns in EIC weak-strong simulation

Beam-Beam simulations are currently being studied in preparation for future EIC experiments to study beam-beam effects and, in turn, maximize luminosity. Weak-strong methods are studied for single-particle dynamics during collision. 1 million macro-particles for 1 million turns are typically tracked, corresponding to only 10 seconds in the EIC. The goal of this study is to predict beam properties over the scale of hours. A potential solution focuses on using data-driven methods such as machine learning methods to analyze and extend the insights of the beam properties such as long-term nonlinear effects. This would aid in long-term predictions where results would be more efficiently acquired than a typical tracking simulation. Some limitations such as inaccurate predictions and spatial complexity are also discussed. These methods can then be applied to strong-strong simulations in the future studies.

Accelerator Physics↗

DOE OSTI2025

Computational Study of a Cryopump Design for Fusion Exhaust Gas Purification Using Direct Simulation Monte Carlo

Cryopump-based direct internal recycling (DIR) of fusion fuel is an attractive prospect because it provides pumping of helium ash from a reactor along with separation of helium and other impurities from the fuel. Previous studies have demonstrated a continuously regenerating cryopump, referred to as the Snail pump, that separates helium ash from fusion fuel for reactor-relevant flow rates. Here, in this study, a conceptual cryopump was designed as a proof-of-principle study to target other impurities, besides helium ash, from the fusion exhaust that would complement the Snail pump. The impurity-removal cryopump consists of four sets of chevron fins, two sets operating at 80 K and the other two at 30 K. A computational study using direct simulation Monte Carlo (DSMC) was performed with a gas mixture of 96.5% D2, 2.0% He, and 0.5% each of CO2, N2, and CH4. In this configuration, 80 K chevron fin sets are capable of capturing CO2 and 30 K sets capable of capturing CO2, N2, and CH4. The computational study showed that the cryopump is capable of reducing the impurity content by more than two orders of magnitude from the flow.

Adhikari, Nirajan [Oak Ridge National Laboratory (↗

DOE OSTI2023

Classroom aerosol dispersion modeling: experimental assessment of a low-cost flow simulation tool

The purpose of this study was to assess the utility of a low-cost flow simulation tool for an indoor air modeling application by comparing its outputs with the results of a physical experiment, as well as those from a more advanced computational fluid dynamics (CFD) software package. In this study, five aerosol dispersion tests were performed in two different classrooms by releasing a CO 2 tracer gas from six student locations. Resultant steady-state concentrations were monitored at 13 locations around the periphery of the room. Subsequently, the experiments were modeled using both a low-cost tool (SolidWorks Flow Simulation) and a more sophisticated tool (STAR-CCM+). Models were evaluated based on their ability to predict the experimentally measured concentrations at the 13 monitoring locations by calculating four performance parameters commonly used in the evaluation of dispersion models: fractional mean bias (FB), normalized mean-square error (NMSE), fraction of predicted value within a factor of two (FAC2), and normalized absolute difference (NAD). The more sophisticated model performed better in 15 of the 20 possible cases (five tests at four parameters each), with parameters meeting acceptance criteria in 19 of 20 cases. However, the lower-cost tool was only slightly worse, with parameters meeting acceptance criteria in 18 of 20 cases, and it performed better than the other tool in 3 of 20 cases. Because it provides useful results at a fraction of the monetary and training cost and is already widely accessible to many institutions, such a tool may be worthwhile for many indoor aerosol dispersion applications, especially for students or researchers just beginning CFD modeling.

54 ENVIRONMENTAL SCIENCES↗