Search NASA⌕ Search

SEARCH · Search NASA

Results for “engine 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗
DOE OSTI2017

Detailed assessment of diesel spray atomization models using visible and X-ray extinction measurements

The physical mechanisms characterizing the breakup of a diesel spray into droplets are still unknown. This gap in knowledge has largely been due to the challenges of directly imaging this process or quantitatively measuring the outcomes of spray breakup, such as droplet size. Recent x-ray measurements by Argonne National Laboratory, utilized in this work, provide needed information about the spatial evolution of droplet sizes in selected regions of the spray under a range of injection pressures (50–150 MPa) and ambient densities (7.6–22.8 kg/m3) relevant for diesel operating conditions. Ultra-small angle x-ray scattering (USAXS) measurements performed at the Advanced Photon Source are presented, which quantify Sauter mean diameters (SMD) within optically thick regions of the spray that are inaccessible by conventional droplet sizing measurement techniques, namely in the near-nozzle region, along the spray centerline, and within the core of the spray. To quantify droplet sizes along the periphery of the spray, a complementary technique is proposed and introduced, which leverages the ratio of path-integrated x-ray and visible laser extinction (SAMR) measurements to quantify SMD. The SAMR and USAXS measurements are then utilized to evaluate current spray models used for engine computational fluid dynamic (CFD) simulations. We explore the ability of a carefully calibrated spray model, premised on aerodynamic wave growth theory, to capture the experimentally observed trends of SMD throughout the spray. The spray structure is best predicted with an aerodynamic primary and secondary breakup process that is represented with a slower time constant and larger formed droplet size than conventionally recommended for diesel spray models. Additionally, spray model predictions suggest that droplet collisions may not influence the resultant droplet size distribution along the spray centerline in downstream regions of the spray.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Predicting extreme anisotropy and shape variations in impact testing of tantalum single crystals

We report recent Taylor cylinder impact tests carried out for Ta single crystals showed strong variations in dimensional changes for different crystallographic directions aligned with the cylindrical axis. In order to capture the effect of crystallography on the deformation characteristics and final shapes of the impacted cylinders, a single crystal material subroutine is adapted and embedded in the solid mechanics/dynamics Finite Element solver Abaqus to simulate the aforementioned single crystal Ta Taylor impact experiments. Details of the coupled model implementation, and insights on the role played by single crystal anisotropic flow on the deformation behavior across a broad range of strain rates and temperatures for different single crystal orientations are presented and discussed. We demonstrate the predictive capability of the adopted crystal plasticity model to capture the significant role played by crystal orientation-induced anisotropy, as well as strain hardening and adiabatic heating, on the dynamic deformation response of crystalline materials. This re-emphasizes the need of microstructure-aware models to improve the accuracy of simulations for high-consequence engineering design.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A modified Embedded-Atom Method interatomic potential for uranium-silicide

Uranium-silicide (U-Si) fuels are being pursued as a possible accident tolerant fuel (ATF). U-Si fuels benefit from higher thermal conductivity and higher fissile density compared to UO2. In order to perform engineering scale nuclear fuel performance simulations, the material properties of the fuel must be known. Currently, the experimental data available for U-Si fuels is rather limited. Thus, multiscale modeling efforts are underway to address this gap in knowledge. In this study, a semi-empirical Modified Embedded Atom Method (MEAM) potential is presented for the description of the U-Si system. The potential accurately describes not only the primary phase of interest (U3Si2), but also a variety of U-Si phases across the composition spectrum.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Molecular Origin of Wettability Alteration of Subsurface Porous Media upon Gas Pressure Variations

Upon extraction/injection of a large quantity of gas from/into a subsurface system in shale gas production or carbon sequestration, the gas pressure varies remarkably, which may significantly change the wettability of porous media involved. Mechanistic understanding of such changes is critical for designing and optimizing a related subsurface engineering process. Using molecular dynamics simulations, in this work we have calculated the contact angle of a water droplet on various solid surfaces (kerogen, pyrophyllite, calcite, gibbsite, and montmorillonite) as a function of CO 2 or CH 4 gas pressure up to 200 atm at a temperature of 300 K. The calculation reveals a complex behavior of surface wettability alteration by gas pressure variation depending on surface chemistry and structure, and molecular interactions of fluid molecules with surfaces. As the CO 2 gas pressure increases, a partially hydrophilic kerogen surface becomes highly hydrophobic, while a calcite surface becomes more hydrophilic. Considering kerogen and calcite being the major components of a shale formation, we postulate that the wettability alteration of a solid surface induced by a gas pressure change may play an important role in fluid flows in shale gas production and geological carbon sequestration.

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

A parametrically programmable delay line for microwave photons

Delay lines that store quantum information are crucial for advancing quantum repeaters and hardware efficient quantum computers. Traditionally, they are realized as extended systems that support wave propagation but provide limited control over the propagating fields. Here, we introduce a parametrically addressed delay line for microwave photons that provides a high level of control over the stored pulses. By parametrically driving a three-wave mixing circuit element that is weakly hybridized with an ensemble of resonators, we engineer a spectral response that simulates that of a physical delay line, while providing fast control over the delay line’s properties. We demonstrate this novel degree of control by choosing which photon echo to emit, translating pulses in time, and even swapping two pulses, all with pulse energies on the order of a single photon. We also measure the noise added from our parametric interactions and find it is much less than one photon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Making metals linear super-elastic with ultralow modulus and nearly zero hysteresis

We demonstrate a novel materials design approach to achieve unprecedented properties by utilizing nanoscale chemo-mechanical coupling. In particular, by using computer simulations we demonstrate how to engineer ultralow modulus (12 GPa), nearly hysteresisfree, and linear super-elastic metals with a giant elastic strain limit (2.7%) by creating appropriate concentration modulations (CMs) at the nanoscale in the parent phase and by pre-straining to regulate the stress-induced martensitic transformation (MT). The nanoscale CMs created via spinodal decomposition produce corresponding phase stability modulations, suppress autocatalysis in nucleation, impose nano-confinements on growth, and hinder long-range ordering of transformation strain during the MT, which changes the otherwise sharp first-order transition into a smeared, macroscopically continuous transition over a large stress range. The pre-straining generates retained martensitic particles that are stable at the test temperature after unloading and act as operational nuclei in subsequent load cycles, eliminating the stress–strain hysteresis and offering an ultralow apparent Young’s modulus. Materials with a high strength and an ultralow apparent Young’s modulus have great potential for application in orthopaedic implants.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Privacy Preservation from High-Performance Computing to Autonomous Science [Industrial and Governmental Activities]

High-Performance Computing (HPC) and Leadership-Class Supercomputing are driving forces behind scientific advancements, enabling researchers to tackle complex challenges in physics, chemistry, biology, and engineering. These systems power vast simulations and data analyses, fueling discoveries in fields ranging from materials science to climate modeling. However, their use often involves processing sensitive data—such as proprietary industry simulations, biomedical records, and national security computations—posing significant privacy concerns. In conclusion, this issue is amplified in collaborative environments like Department of Energy (DOE) user facilities, where HPC resources are shared across institutions to foster innovation.

Kotevska, Olivera [Oak Ridge National Laboratory (↗
DOE OSTI2017

Modeling defect and fission gas properties in U-Si fuels

Uranium silicides, in particular U 3 Si 2 , are being explored as an advanced nuclear fuel with increased accident tolerance as well as competitive economics compared to the baseline UO2 fuel. They benefit from high thermal conductivity (metallic) compared to UO 2 fuel (insulator or semi-conductor) used in current Light Water Reactors (LWRs). The U-Si fuels also have higher fissile density. In order to perform meaningful engineering scale nuclear fuel performance simulations, the material properties of the fuel, including the response to irradiation environments, must be known. Unfortunately, the data available for USi fuels are rather limited, in particular for the temperature range where LWRs would operate. The ATF HIP is using multi-scale modeling and simulations to address this knowledge gap.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Modeling defect and fission gas properties in U-Si fuels

Uranium silicides, in particular U 3 Si 2 , are being explored as an advanced nuclear fuel with increased accident tolerance as well as competitive economics compared to the baseline UO 2 fuel. They benefit from high thermal conductivity (metallic) compared to UO 2 fuel (insulator or semi-conductor) used in current Light Water Reactors (LWRs). The U-Si fuels also have higher fissile density. In order to perform meaningful engineering scale nuclear fuel performance simulations, the material properties of the fuel, including the response to irradiation environments, must be known. Unfortunately, the data available for USi fuels are rather limited, in particular for the temperature range where LWRs would operate. The ATF HIP is using multi-scale modeling and simulations to address this knowledge gap.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Comoparison of Slit Trench and Engineered Trench 3D Vadose Zone Conceptual Models

In the next revision of the E-Area Low-Level Waste Facility (ELLWF) Performance Assessment (PA), radionuclide groundwater transport from slit (ST) and engineered trenches (ET) will be simulated using 3-dimensional PORFLOW models representing the vadose zone. In total, the conceptual framework will aim to model 29 total trenches with a degree of conservatism that provides a high level of confidence in meeting performance objectives. The implementation of generic conceptual models allows for a reasonable number of simulation scenarios, thereby reducing the chances of errors and enabling a better one-to-one comparison between the results, while still capturing the most important features of each trench. The current investigation seeks to understand if a generic trench model based on the general ET footprint can justifiably be used for representing both STs and ETs while maintaining a reasonable degree of conservatism. In this process, the groundwater transport of several key radionuclides was simulated for both ET and ST geometries and the flux-to-the-water-table profiles were compared. For all radionuclides simulated, the ET geometry is shown to be conservative (i.e., produces the highest peak flux to the water table).

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

Refinement and Modeling of a Blackbody-Based Calibration Method in the InfraBREAD Detector

The Broadband Reflector Experiment for Axion Detection (BREAD) is an ongoing experiment searching for the conversion of yet undiscovered axion-like dark matter particles to photons in the presence of a magnetic field. InfraBREAD, a sub-experiment of BREAD, uses a superconducting nanowire single photon detector (SNSPD), a high-efficiency and low-noise device, to specifically detect infrared-range photons produced by $\mathcal{O}$(eV) axion-like particles. The unique BREAD reflector setup allows for the focusing of all converted signal photons to a 1mm $\times$ 1mm SNSPD. However, when the detector is cooled to cryogenic temperatures during operation, uneven thermal contraction of reflector components may lead to a small shift in the location of the focal spot. A novel calibration method using blackbody radiation is proposed to locate the true focal spot of the detector \textit{in situ}. Through ray tracing simulations done in FRED Optical Engineering Software, this method is demonstrated to locate the focus to within 50 $\mu$m in three dimensions.

Rao, Shardul↗
DOE OSTI2022

Materials Design using an Active Subspace-based Batch Bayesian Optimization Approach

Integrated computational materials engineering (ICME) calls for integrating simulation tools and/or experiments to develop new materials and materials systems. However, implementation of ICME approaches is challenging mainly due to the considerable computational expense of such frameworks and large dimensionality of the design space. Addressing these challenges is thus critical to the success of ICME initiatives. We present here a specific Bayesian optimization framework designed to address these two challenges. In particular, we propose an active subspace batch Bayesian optimization framework. The framework makes use of dimension reduction via the active subspace method and makes use of the ability to query in parallel via the batch Bayesian optimization approach. Here, the integration of these techniques leads to significant efficiency improvements while maintaining accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

DME-Propane Blends Ignition Experiments and Modeling for Heavy-Duty Mixing Controlled Compression Engines

A blend of dimethyl ether (DME) and propane (C3H8) is being studied in a shock tube at heavy-duty engine conditions at 110 bar. Due to its intrinsic combustion properties, DME/propane blend can potentially replace diesel in mixing controlled compression ignition engines. A blend of DME/propane can reduce emissions in mixing controlled compression ignition in heavy-duty engines through modifications, which require simulations using a high-fidelity chemical kinetics model that can accurately predict the chemistry of the blend. An essential aspect of testing the chemical kinetics model is doing baseline fundamental chemistry studies on neat DME and propane, which include ignition delay time measurements. In this work, using a high-pressure shock tube, ignition delay times were gathered for DME/Propane blends at 110 bar diluted with AR to test chemical kinetic models published in the literature. These models include Aramco 3.0, NUIG V1.1, C3mech V3.3, and Dames et al. Comparisons with the experimental IDTs and models were conducted, and general agreement was observed. A sensitivity analysis was conducted, and important reactions were outlined.

Mohammed, Zuhayr Pasha [University of Central Flor↗
DOE OSTI2018

Mass imbalances in EPANET water-quality simulations

EPANET is widely employed to simulate water quality in water distribution systems. However, in general, the time-driven simulation approach used to determine concentrations of water-quality constituents provides accurate results only for short water-quality time steps. Overly long time steps can yield errors in concentration estimates and can result in situations in which constituent mass is not conserved. The use of a time step that is sufficiently short to avoid these problems may not always be feasible. The absence of EPANET errors or warnings does not ensure conservation of mass. This paper provides examples illustrating mass imbalances and explains how such imbalances can occur because of fundamental limitations in the water-quality routing algorithm used in EPANET. In general, these limitations cannot be overcome by the use of improved water-quality modeling practices. This paper also presents a preliminary event-driven approach that conserves mass with a water-quality time step that is as long as the hydraulic time step. Results obtained using the current approach converge, or tend to converge, toward those obtained using the preliminary event-driven approach as the water-quality time step decreases. Improving the water-quality routing algorithm used in EPANET could eliminate mass imbalances and related errors in estimated concentrations. Furthermore, the results presented in this paper should be of value to those who perform water-quality simulations using EPANET or use the results of such simulations, including utility managers and engineers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

DEVELOPING A TRAFFIC SIGNAL ECO-COOPERATIVE ADAPTIVE CRUISE CONTROL SYSTEM FOR BATTERY ELECTRIC VEHICLES

This study developed an eco-driving strategy for battery electric vehicles (BEVs) in the vicinity of signalized intersections, known as BEV Eco-Cooperative Adaptive Cruise Control at Intersections (Eco-CACC-I). The system computes a real-time, energy-optimized vehicle trajectory using a BEV vehicle dynamics and energy consumption model. In the proposed approach, the relationship between vehicle speed, acceleration level, and tractive/resistance forces on the vehicle follows the vehicle dynamics model. The Virginia Tech Comprehensive Power-based Electric Vehicle Energy Consumption Model (VT-CPEM) is used to compute the instantaneous battery energy consumption, with the consideration of battery energy regeneration during braking. The Eco-CACC-I energy-optimum trajectory computation is formulated as an optimization problem with constraints. Considering that the optimization solution needs to be computed at a rapid frequency (e.g., 10 Hz) for real-time applications, an A-star algorithm was used here to expedite the computational speed. Lastly, the proposed BEV Eco-CACC-I system was tested in simulation by considering the impacts of various road grades, speed limits, and signal timings. Simulations for an Internal Combustion Engine Vehicles (ICEVs) was conducted to compare with the BEV test results. The results of the comparison indicate that the energy-optimum solution for BEVs is different from that for ICEVs and thus demonstrating the need for vehicle-tailored optimum trajectories.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗