Engineering Papers⌕ Search

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.

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

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Control Design and EMT Simulation Toward the Restoration of Faulted Wind-Dominant Grids

Present power engineers have been challenged to employ electromagnetic transient (EMT) simulations to study the restoration of wind-dominant grids. This paper addresses this hurdle by developing strategies to: (i) reliably control the dc-link voltage of grid-forming wind turbines with dc-coupled batteries, (ii) tune controllers to stably interconnect wind power plants into a grid under recovery, and (iii) autonomously make restoration decisions if transmission faults occur during restoration. These contributions are assessed via EMT simulations of wind-dominant versions of the WSCC 9-bus grid and a 22-bus power system. Furthermore, this paper is significant to address recommendations by the North American Electric Reliability Corporation.

17 WIND ENERGY↗
DOE OSTI2021

Crystallization of the P 3 Sn 4 Phase upon Cooling P 2 Sn 5 Liquid by Molecular Dynamics Simulation Using a Machine Learning Interatomic Potential

We performed molecular dynamics simulations to study the crystallization of the P3Sn4 phase from P 2 Sn 5 liquid using a machine learning (ML) interatomic potential with desirable efficiency and accuracy. Our results capture the liquid properties of P 2 Sn 5 at 1300 K, which is well above the melting temperature. The phase separation and crystallization are observed when P 2 Sn 5 liquid is cooled down below 832 and 505 K, respectively. The simulation results are in good agreement with the experimentally observed phase transformation behaviors and provide useful insights into the complex nucleation and crystallization process at the details of atomistic scale. Our work also demonstrated that ML interatomic potentials based on neural network deep learning are robust and capable of accurately describing the energetics and kinetics of complex materials through molecular dynamics simulations.

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

Gyrokinetic Particle Simulations of Core Turbulence in Fusion Plasmas

In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Impact of Regional and Seasonal Characteristics on Battery Electric Vehicle Operational Costs in the U.S.

This study investigates the operational cost competitiveness of battery electric vehicles (BEVs) in the United States, considering regional climates, energy prices, and driving patterns. By comparing BEVs with plug-in hybrid electric vehicles (PHEVs), hybrid electric vehicles (HEVs), and the alternative use of BEVs and conventional vehicles (Convs), the analysis incorporates thermal dynamometer tests, real-world vehicle miles traveled (VMT), and state-specific energy prices. Using detailed simulations, the study evaluates energy consumption across varying temperatures and driving distances. The findings reveal that, while BEVs remain cost-effective for short trips in moderate climates, PHEVs are more economical for long-range trips and cold environments, due to the excessive cost of using external direct current fast chargers (DCFCs) and reduced BEV efficiency at low temperatures. HEVs are identified as the most cost-efficient option in regions like New England, characterized by high residential electricity prices. These insights are critical for shaping vehicle electrification strategies, particularly under diverse regional and seasonal conditions, and for advancing policies on alternative energy and fuels.

Kim, Kyung-Ho (ORCID:0009000450851732)↗
DOE OSTI2020

Quantifying CanESM5 and EAMv1 sensitivities to Mt. Pinatubo volcanic forcing for the CMIP6 historical experiment

Large volcanic eruptions reaching the stratosphere have caused marked perturbations to the global climate including cooling at the Earth's surface, changes in large-scale circulation and precipitation patterns and marked temporary reductions in global ocean heat content. Many studies have investigated these effects using climate models; however, uncertainties remain in the modelled response to these eruptions. This is due in part to the diversity of forcing datasets that are used to prescribe the distribution of stratospheric aerosols resulting from these volcanic eruptions, as well as uncertainties in optical property derivations from these datasets. To improve this situation for the sixth phase of the Coupled Model Intercomparison Project (CMIP6), a two-step process was undertaken. First, a combined stratospheric aerosol dataset, the Global Space-based Stratospheric Aerosol Climatology (GloSSAC; 1979–2016), was constructed. Next, GloSSAC, along with information from ice cores and Sun photometers, was used to generate aerosol distributions, characteristics and optical properties to construct a more consistent stratospheric aerosol forcing dataset for models participating in CMIP6. This “version 3” of the stratospheric aerosol forcing has been endorsed for use in all contributing CMIP6 simulations. Recent updates to the underlying GloSSAC from version 1 to version 1.1 affected the 1991–1994 period and necessitated an update to the stratospheric aerosol forcing from version 3 to version 4. As version 3 remains the official CMIP6 input, quantification of the impact on radiative forcing and climate is both relevant and timely for interpreting results from experiments such as the CMIP6 historical simulations. This study uses two models, the Canadian Earth System Model version 5 (CanESM5) and the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1), to estimate the difference in instantaneous radiative forcing in simulated post-Pinatubo climate response when using version 4 instead of version 3. Differences in temperature, precipitation and radiative forcings are generally found to be small compared to internal variability. An exception to this is differences in monthly temperature anomalies near 24 km altitude in the tropics, which can be as large as 3 °C following the eruption of Mt. Pinatubo.

58 GEOSCIENCES↗
DOE OSTI2020

Assessment of Methane Generation from the Thermolysis of Organics in Caustic Waste

Methane generation has been detected during thermolytic Hydrogen Generation Rate (HGR) testing of SRS tank waste and in simulant samples. Simulant testing has shown methane generation in the presence of trimethylsilanol (TMS) and radioactive sample testing has shown methane generation in presence of methylmercury, both of which are found within SRS waste tanks. The 2019 report Evaluation of the Current State of Knowledge for Thermolysis of Organics within SRS Waste Forming Volatile Organic Compounds (VOCs) reviewed historical data and assessed literature pertaining to thermolytic methane generation. That review recommended undertaking thermolysis studies on simulated waste to detect and quantify VOCs. The purpose of this assessment is to review chemistry and current knowledge of organics in SRS waste to recommend which organics in SRS waste tanks, if any, should be tested for thermolytic methane generation. An evaluation of observed methane generation rates has been performed to assess the implication for current CSTF controls. Based on the inputs used in this report, it was concluded that no methane generation data observed to date has challenged the existing 5% organic contribution to the Composite Lower Flammability Limit (CLFL) assumption in the CSTF tank farm at temperatures below 100 °C. Evidence shows that thermolytic pathways to methane generation become significant at temperatures >100 °C in caustic waste. Currently, the influences of thermolytic methane production on the flammability of vessels at >100 °C (e.g., evaporators) is unclear. It is recommended that testing should be undertaken using the HGR experimental approach for methane generation. The initial simulant testing should obtain quantitative methane generation rates at elevated temperatures for methylated siloxanes/silanes and methylmercury with the goal of determining the current mechanism for thermolytic methane generation seen in SRS tank waste. The results of this testing should then be compiled and analyzed to ascertain the risk of challenging the current projected organic contribution of 5% to the CLFL at elevated temperatures. Following this testing, an evaluation of possible flammable gas generation in the 2H and 3H Evaporators should be performed to assess the risk of methane generation at the high temperatures employed by the evaporators.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗
DOE OSTI2026

Passive Confirmation of the Presence Of High-Explosive Material Via Neutron Transmission Spectroscopy

The goal of this project was to explore a novel approach for identifying presence and potential types of high explosives (HE) in treaty-controlled items with passive neutron sources by using passive neutron transmission spectroscopy. We predominantly look to measure relative elemental abundance of Carbon, Hydrogen, Oxygen, and Nitrogen (CHON) since most relevant materials are certain mix of these elements, as shown in Table 1. Previously, most studies exploring potential identification of CHON elemental content of targets in the vicinity of passive neutron source were focused solely on gamma spectroscopy using high resolution gamma detectors. Using neutron transmission spectroscopy with pulse-shape discrimination (PSD) capable organic scintillators is not only a novel idea in of itself, but arguably necessary for accurate identification of the CHON elemental content of an interrogated target. While high resolution gamma spectroscopy in principle can determine a presence of hydrogen and nitrogen, it is effectively blind to a quantitatively measuring areal densities of C, O. The initially proposed approach, described in detail in the project proposal, was to leverage previous Monte Carlo study on neutron transmission spectroscopy using slab shaped target interrogated with well-collimated (“pencil”) neutron beam. In this study, a simulated test CO 2 target was interrogated by a pencil neutron beam and transmitted neutron spectra were measured using PSD capable organic scintillator using MLEM based unfolding technique. To establish relative elemental content of carbon and oxygen, the unfolded neutron spectrum was fit with a parametrized combination of their respective elemental spectral templates. The individual spectral template was calculated by convoluting ENDF neutron crosssection with the known resolution of the PSD detector used in the study. This approach in the studied configuration successfully quantitatively established relative elemental composition of carbon and oxygen.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Blending as an In-Tank Waste Processing Strategy: Year 1 Progress Report

The Hanford Site stores approximately 56 million gallons of radioactive legacy defense waste in underground tanks. Direct feed of southeast quadrant tank waste to a high-level waste treatment facility is currently being explored to accelerate tank treatment activities. This study seeks to de-risk waste transfers by favorably altering the resuspension and critical velocities of wastes through strategic waste blending. Wastes rich in aluminum- and zirconium-bearing solids pose significant risks due to their high settled strength and density. A series of tests explores how solid phase particle size and suspending phase chemistry alter waste transport properties, quantified through the just-suspended mixing speed (NJS), for pure and blended aluminum- and zirconium-bearing waste solid simulants. The study utilizes a range of gibbsites (three sizes spanning 10 to 90 µm) and zirconia solids (three sizes spanning 10 to 100 µm). As expected, the NJS of pure solids largely follows predictions based on literature correlations; blends of gibbsite and zirconium-bearing solids of differing particle sizes exhibit NJS maxima attributed to dense packing of the mixed-size settled solids. The findings of this study will assist in developing strategies to optimize the transport of waste solids for the safe and efficient delivery of tank waste.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Molecular dynamics simulation of metallic Al–Ce liquids using a neural network machine learning interatomic potential

Al-rich Al-Ce alloys have the possibility of replacing heavier steel and cast-irons for use in high-temperature applications. Knowledge about the structures and properties of Al-Ce alloys at liquid state is vital for optimizing the manufacture process to produce desired allows. However, reliable molecular dynamics simulation of Al-Ce alloy systems remains a great challenge due to the lack of accurate Al-Ce interatomic potential. In this work, an artificial neural network (ANN) deep machine learning (ML) method is used to develop a reliable interatomic potential for Al-Ce alloy. Ab initio molecular dynamics (AIMD) simulation data on Al-Ce liquid with small unit cell (~200 atoms) and on the known Al-Ce crystalline compounds are collected to train the interatomic potential using ANN-ML. The obtained ANN-ML model reproduces well the energies, forces, and atomic structure of Al 90 Ce 10 liquid and crystalline phases of Al-Ce compounds in comparison with ab initio results. The developed ANN-ML potential is applied in molecular dynamics simulations to study the structures and properties of metallic Al 90 Ce 10 liquid, which would provide useful insight for guiding experimental process to produce desired Al-Ce allows.

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

Probe Chiral Magnetic Effect with Signed Balance Function

In this paper a pair of observables are proposed as alternative ways, by examining the fluctuation of net momentum-ordering of charged pairs, to study the charge separation induced by the Chiral Magnetic Effect (CME) in relativistic heavy ion collisions. They are, the out-of-plane to in-plane ratio of fluctuation of the difference between signed balance functions measured in pair’s rest frame, and the ratio of it to similar measurement made in the laboratory frame. In this work, both observables have been studied with simulations including flow-related backgrounds, and for the first time, backgrounds that are related to resonance's global spin alignment. The two observables have similar positive responses to signal, and opposite, limited responses to identifiable backgrounds arising from resonance flow and spin alignment. Both observables have also been tested with two realistic models, namely, a multi-phase transport (AMPT) model and the anomalous-viscous fluid dynamics (AVFD) model. These two observables, when cross examined, will provide useful insights in the study of CME-induced charge separation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Self‐Healable Poly(ionic liquid) Copolymers Driven by Polar and Dipolar Forces

Abstract Commodity aliphatic and aromatic acrylic‐based copolymers self‐heal due to ubiquitouskey‐and‐lock,ring‐and‐lock, andfluorophilic‐σ‐lockvan der Waals (vdW) interactions. However, the role of these interactions in the presence of covalently copolymerized ionic liquid (IL) is not known. This study is driven by the hypothesis that covalently incorporated cation–anion pairs to form poly(ionic liquid) copolymers (PILCs) can perturb inter‐ or intra‐chain vdW interactions reflected in mechanical and electrical responses. To test this hypothesis, we synthesized a series of PILCs comprising of pentafluorostyrene (PFS) and imidazolium‐based IL monomers with variable‐length aliphatic tails (methyl and butyl). Using a combination of 2D 1 H‐ 1 H and 19 F ‐ 19 F NOESY NMR and FTIR measurements supplemented by molecular dynamic (MD) simulations, these studies demonstrate that preferentially alternating/random PILCs topologies facilitate self‐healing. The introduction of cation–anion moieties modifies thefluorophilic‐σ‐lockinteractions and, along with longer aliphatic tails ─(CH 2 ) 3 CH 3 covalently attached to the imidazolium cation, enhances cation‐anion mobility, thus faster recovery from mechanical damage occurs. These findings underline how precise control over dipolar and ionic interactions through copolymer composition enables self‐healing in PILCs. These insights may open pathways for designing sustainable, mechanically resilient materials for applications in energy storage and energy harvesting.

Chemistry↗
DOE OSTI2020

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

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

Carbon nanotube porin diffusion in mixed composition supported lipid bilayers

Carbon nanotube porins (CNTPs), short pieces of carbon nanotubes capable of self-inserting into a lipid bilayer, represent a simplified model of biological membrane channels. We have used high-speed atomic force microscopy (HS-AFM) and all-atom molecular dynamics (MD) simulations to study the behavior of CNTPs in a mixed lipid membrane consisting of DOPC lipid with a variable percentage of DMPC lipid added to it. HS-AFM data reveal that the CNTPs undergo diffusive motion in the bilayer plane. Motion trajectories extracted from the HS-AFM movies indicate that CNTPs exhibit diffusion coefficient values broadly similar to values reported for membrane proteins in supported lipid bilayers. The data also indicate that increasing the percentage of DMPC leads to a marked slowing of CNTP diffusion. MD simulations reveal a CNTP-lipid assembly that diffuses in the membrane and show trends that are consistent with the experimental observations.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

The structure of dissipative dark matter halos

Dissipative dark matter self-interactions can affect halo evolution and change its structure. We perform a series of controlled N-body simulations to study impacts of the dissipative interactions on halo properties. The interplay between gravitational contraction and collisional dissipation can significantly speed up the onset of gravothermal collapse, resulting in a steep inner density profile. For reasonable choices of model parameters controlling the dissipation, the collapse timescale can be a factor of 10–100 shorter than that predicted in purely elastic self-interacting dark matter. The effect is maximized when energy loss per collision is comparable to characteristic kinetic energy of dark matter particles in the halo. Our simulations provide guidance for testing the dissipative nature of dark matter with astrophysical observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Decadal Thermal Variability of the Upper Southern Ocean: Zonal Asymmetry

Abstract As the major sink of anthropogenic heat, the Southern Ocean has shown quasi-symmetric, deep-reaching warming since the mid-twentieth century. In comparison, the shorter-term heat storage pattern of the Southern Ocean is more complex and has notable impacts on regional climate and marine ecosystems. By analyzing observational datasets and climate model simulations, this study reveals that the Southern Ocean exhibits prominent decadal (>8 years) variability extending to ∼700-m depth and is characterized by out-of-phase changes in the Pacific and Atlantic–Indian Ocean sectors. Changes in the Pacific sector are larger in magnitude than those in the Atlantic–Indian Ocean sectors and dominate the total heat storage of the Southern Ocean on decadal time scales. Instead of heat uptake through surface heat fluxes, these asymmetric variations arise primarily from wind-driven heat redistribution. Pacemaker and preindustrial simulations of the Community Earth System Model version 1 (CESM1) suggest that these variations in Southern Ocean winds arise primarily from natural variability of the tropical Pacific, as represented by the interdecadal Pacific oscillation (IPO). Through atmospheric teleconnection, the positive phase of the IPO gives rise to higher-than-normal sea level pressure and anticyclonic wind anomalies in the 50°–70°S band of the Pacific sector. These winds lead to warming of 0–700 m by driving the convergence of warm water. The opposite processes, involving cyclonic winds and upper-layer divergence, occur in the Atlantic–Indian Ocean sector. These findings aid our understanding of the time-varying heat storage of the Southern Ocean and provide useful implications on initialized decadal climate prediction.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Focusing Optics for Axion Detection: Simulating Sensing Enhancements of Photons in InfraBREAD

The Broadband Reflector Experiment for Axion Detection (BREAD) will search for axions, a promising particle candidate for dark matter, in the high mass range of [$10^{-3}$,1]eV. BREAD acts like a telescope for dark matter, using a parabolic reflector to direct axion-induced photons onto a precise focal spot. InfraBREAD is the experimental version designed to detect photons resulting from axions with infrared frequency using a 1mm × 1mm Superconducting Nanowire Singe Photon Detector (SNSPD). However, at this frequency scale two problems arise: 1) the focal spot becomes smeared across an area larger than the size of the SNSPD detector; 2) any millimeter-scale misalignment shifts in the SNSPD can greatly diminish photon signal. This study performed simulations of novel optical configurations of lenses and reflectors placed around the SNSPD to achieve a greater efficiency than the SNSPD alone through both of these effects. One optical setup was identified called the “Parabolo id to Winston Cone (PTW) Configuration,” which exploits the optics of parabolic reflectors to focus incoming photons from the larger area of the smeared focal spot onto the smaller 1mm 2 area of the detector. The PTW configuration demonstrated overall improvements in efficiency by around 10% through misalignment shifts, with a maximal improvement of 55% in one region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗