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 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An experimental and computational analysis of combustion heat release transformation in dual fuel combustion

Dual fuel (DF) diesel-methane combustion, which employs a high-reactivity fuel (diesel) to ignite a low-reactivity fuel (methane), is a widely studied combustion strategy for internal combustion engines, with significant potential for engine-out emissions reductions without the need for major hardware modifications. A phenomenon, which has been reported in the DF literature, but not explained fully, is the transformation of the shape of the apparent heat release rate (AHRR) curve as the start of injection (SOI) of diesel is advanced beyond a certain threshold; coincidentally, this AHRR transformation is usually accompanied by a sharp decrease in engine-out emissions of oxides of nitrogen (NOx). The goal of the present work is to establish the underlying physical reason(s) that cause the AHRR transformation. The AHRR transformation was observed on a single cylinder research engine (SCRE) at an indicated mean effective pressure (IMEP) of 5 bar at a speed of 1500 rev/min. The transformation occurred over a range of SOIs from 330 to 320 crank angle degrees (CAD). While the 330 CAD SOI exhibited a typical two-stage AHRR curve, with a clearly definable first-stage peak followed by a second-stage AHRR with little-to-no low temperature heat release (LTHR) present and high engine-out NOx, the 320 CAD SOI exhibited a single-stage, Gaussian-like AHRR curve, with noticeable LTHR and at least one order-of-magnitude lower NOx emissions. Here, leveraging analysis of experimental data and three-dimensional computational fluid dynamic simulations, the authors show that the AHRR transformation is impacted mainly by differences in local equivalence ratio distributions within the cylinder at ignition onset for different diesel SOIs.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

Three-dimensional (3D) analysis of white etching bands (WEBs) in AISI M50 bearing steel using automated serial sectioning

Dark etching regions (DER) and white etching bands (WEB) form in the subsurface stressed regions of bearing raceways and rolling elements. Also, “butterfly” formations form around non-metallic inclusions and large carbides. The formation of WEBs, DERs and butterfly are dependent on operating conditions, material alloying elements and heat treatment. The DER, WEB and butterfly formations are well studied in two dimensions (2D). The objective of this study was to analyze WEB formation in three dimensions (3D) that will allow better understanding of material transformation after stress cycling. Three rolling contact fatigue (RCF) tested vacuum induction melted vacuum arc re-melted (VIM-VAR) AISI M50 samples (from a bench test, sub-scale test and commercial gas turbine engine) were analyzed using automated serial sectioning to map WEB networks and a butterfly formation in 3D. Automated serial sectioning and 2D optical microstructural data for 3D reconstructions and analysis was done using a fully automated serial sectioning system – Robo-Met.3D®. 3D analysis was performed using Image-Pro Premier 3D. The data set allowed 3D visualization and analysis of the evolution and propagation of WEB and butterfly in M50 bearing steel, including the positive microscopic identification of a ceramic inclusion that was the source of a WEB.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electronic Excitation‐Driven β‐Ga 2 O 3 Metastability Transformation and Self‐Organization Mechanism: β→κ/γ/δ Phases

Here, irradiation-driven multiphase self-organization presents emergent opportunities for the customization of nanoscale engineering properties, dynamically tuning strain-field distributions and interfacial electronic structures. Responding to intense electronic excitation-induced energy deposition, the dominant phase transformations, with varying Gibbs free energy $Δ,G^o_f$ are confirmed as β → κ → γ → δ that are located in specific microregions for Gallium (III) oxide (Ga 2 O 3 ), as follows: (i) Surface-localized interstitial accumulation under compressive stress triggers β → δ via semi-coherent interface formation. (ii) Tensile stress within latent tracks drives vacancy-mediated oxygen layer truncation (4/12 periodicity along ⟨0001⟩), stabilizing coherent 4H (ABCB) κ and 3C (ABC) β (ABC) interfaces through strain-compensated octahedral distortion. (iii) Screw dislocation-mediated lattice relaxation induces β → γ via cation disordering (Ga 3+ occupancy at β-interstitial sites), forming metastable spinel γ with mixed occupancy across 16d/8a Wyckoff sites. Irradiation-driven β-Ga 2 O 3 →κ/γ/δ transitions, as mechanistically revealed via inelastic thermal spike (i-TS) calculations and molecular dynamics simulations, induce defect-mediated nonlinear photoresponse, critical for optoelectronic engineering.

Han, Xinqing [Shandong Univ., Jinan (China)]↗
DOE OSTI2021

Reinforcement learning in discrete action space applied to inverse defect design

Abstract Reinforcement learning (RL) algorithms that include Monte Carlo Tree Search (MCTS) have found tremendous success in computer games such as Go, Shiga and Chess. Such learning algorithms have demonstrated super-human capabilities in navigating through an exhaustive discrete action search space. Motivated by their success in computer games, we demonstrate that RL can be applied to inverse materials design problems. We deploy RL for a representative case of the optimal atomic scale inverse design of extended defects via rearrangement of chalcogen (e.g. S) vacancies in 2D transition metal dichalcogenides (e.g. MoS 2 ). These defect rearrangements and their dynamics are important from the perspective of tunable phase transition in 2D materials i.e. 2H (semi-conducting) to 1T (metallic) in MoS 2 . We demonstrate the ability of MCTS interfaced with a reactive molecular dynamics simulator to efficiently sample the defect phase space and perform inverse design—starting from randomly distributed S vacancies, the optimal defect rearrangement of defects corresponds a line defect of S vacancies. We compare MCTS performance with evolutionary optimization i.e. genetic algorithms and show that MCTS converges to a better optimal solution (lower objective) and in fewer evaluations compared to GA. We also comprehensively evaluate and discuss the effect of MCTS hyperparameters on the convergence to solution. Overall, our study demonstrates the effectives of using RL approaches that operate in discrete action space for inverse defect design problems.

42 ENGINEERING↗
DOE OSTI2023

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Overview of lower length scale model development for accident tolerant fuels regarding U 3 Si 2 fuel and FeCrAl cladding

U 3 Si 2 and FeCrAl have been proposed as fuel and cladding concepts, respectively, for accident tolerance fuels with higher tolerance to accident scenarios compared to UO 2 . However, a lot of key physics and material properties regarding their in-pile performance are yet to be explored. To accelerate the understanding and reduce the cost of experimental studies, multiscale modeling and simulation are used to develop physics-based materials models to assist engineering scale fuel performance modeling. In this report, the lower-length-scale efforts in method and material model development supported by the Accident Tolerance Fuel (ATF) high-impact-problem (HIP) under the NEAMS program are summarized. Significant progresses have been made regarding interatomic potential, phase field models for phase decomposition and gas bubble formation, and thermal conductivity for U 3 Si 2 fuel, and precipitation in FeCrAl cladding. The accomplishments are very useful by providing atomistic and mesoscale tools, improving the current understanding, and delivering engineering scale models for these two ATF concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

White paper on plasma simulation

Increases in computer power combined with improved algorithms have considerably improved the ability of simulations to model plasma experiments. Higher fidelity simulations with more complete physics treatment and better resolution have improved our understanding of the physics controlling the performance of these experiments. For example, in inertial confinement fusion (ICF) simulations of indirect drive ignition targets and capsules can now include the broad range of perturbation sources, including surface imperfections, engineering features and drive asymmetries together in a fully 3D geometry. These have significantly improved the ability of simulations to match a suite of experimental observations across a range of targets.1-3 Simulations of direct drive targets and magnetically driven targets have also increased in fidelity and realism. The US is the leader in its capabilities to perform plasma simulations. This leadership has highly important implications for national security as well as future energy security and economic competitiveness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗
DOE OSTI2015

Initiation of fission chains during delayed criticality accidents - 14366

We present the simulation of the waiting time probability distribution before a prescribed reaction rate in a delayed supercritical system. The Kolmogorov backward equations for the count number probability distribution taking into account the delayed neutron precursors are solved within the time dependent point model framework. The initiation of fission chains is also studied using analog Monte-Carlo simulations of the prompt neutrons and delayed neutron precursors population. The fission rate is calculated for comparison with deterministic simulations. Numerical results are presented using both deterministic and Monte-Carlo calculations to study the effect of the source strength and reactivity level. (authors)

42 ENGINEERING↗
DOE OSTI2022

An incompressible Eulerian method for fluid–structure interaction with mixed soft and rigid solids

Here, we present a general simulation approach for incompressible fluid–structure interactions in a fully Eulerian framework using the reference map technique. The approach is suitable for modeling one or more rigid or finitely deformable objects or soft objects with rigid components interacting with the fluid and with each other. It is also extended to control the kinematics of structures in fluids. The model is based on our previous Eulerian fluid–soft solver [Rycroft et al., “Reference map technique for incompressible fluid–structure interaction,” J. Fluid Mech. 898, A9 (2020)] and generalized to rigid structures by constraining the deformation-rate tensor in a projection framework. Several numerical examples are presented to illustrate the capability of the method.

42 ENGINEERING↗
DOE OSTI2021

The IDAES process modeling framework and model library—Flexibility for process simulation and optimization

Abstract Energy systems and manufacturing processes of the 21st century are becoming increasingly dynamic and interconnected, which require new capabilities to effectively model and optimize their design and operations. Such next generation computational tools must leverage state‐of‐the‐art techniques in optimization and be able to rapidly incorporate new advances. To address these requirements, we have developed the Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform, which builds on the strengths of both process simulators (model libraries) and algebraic modeling languages (advanced solvers). This paper specifically presents the IDAES Core Modeling Framework (IDAES‐CMF), along with a case study demonstrating the application of the framework to solve process optimization problems. Capabilities provided by this framework include a flexible, modifiable, open‐source platform for optimization of process flowsheets utilizing state‐of‐the‐art solvers and solution techniques, fully open and extensible libraries of dynamic unit operations models and thermophysical property models, and integrated support for superstructure‐based conceptual design and optimization under uncertainty.

42 ENGINEERING↗
DOE OSTI2022

A Framework for Inverse Prediction Using Functional Response Data

Inverse prediction models have commonly been developed to handle scalar data from physical experiments. However, it is not uncommon for data to be collected in functional form. When data are collected in functional form, it must be aggregated to fit the form of traditional methods, which often results in a loss of information. For expensive experiments, this loss of information can be costly. In this study, we introduce the functional inverse prediction (FIP) framework, a general approach which uses the full information in functional response data to provide inverse predictions with probabilistic prediction uncertainties obtained with the bootstrap. The FIP framework is a general methodology that can be modified by practitioners to accommodate many different applications and types of data. We demonstrate the framework, highlighting points of flexibility, with a simulation example and applications to weather data and to nuclear forensics. Results show how functional models can improve the accuracy and precision of predictions.

42 ENGINEERING↗
DOE OSTI2016

Dakota Uncertainty Quantification Methods Applied to the CFD code Nek5000

This report presents the state of advancement of a Nuclear Energy Advanced Modeling and Simulation (NEAMS) project to characterize the uncertainty of the computational fluid dynamics (CFD) code Nek5000 using the Dakota package for flows encountered in the nuclear engineering industry. Nek5000 is a high-order spectral element CFD code developed at Argonne National Laboratory for high-resolution spectral-filtered large eddy simulations (LESs) and unsteady Reynolds-averaged Navier-Stokes (URANS) simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Considerations for the temperature stratification in a pre-burn constant-volume combustion chamber

In recent years, the Engine Combustion Network (ECN) has developed as a worldwide reference for understanding and describing engine combustion processes, successfully bringing together experimental and numerical efforts. Since experiments and numerical simulations both target the same boundary conditions, an accurate characterization of the stratified environment that is inevitably present in experimental facilities is required. The difference between the core-, and pressure-derived bulk-temperature of pre-burn combustion vessels has been addressed in various previous publications. Additionally, thermocouple measurements have provided initial data on the boundary layer close to the injector nozzle, showing a transition to reduced ambient temperatures. The conditions at the start of fuel injection influence physicochemical properties of a fuel spray, including near nozzle mixing, heat release computations, and combustion parameters. To address the temperature stratification in more detail, thermocouple measurements at larger distances from the spray axis have been conducted. Both the temperature field prior to the pre-combustion event that preconditions the high-temperature, high-pressure ambient, as well as the stratification at the moment of fuel injection were studied. To reveal the cold boundary layer near the injector with a better spatial resolution, Rayleigh scattering experiments and thermocouple measurements at various distances close to the nozzle have been carried out. The impact of the boundary layers and temperature stratification are illustrated and quantified using numerical simulations at Spray A conditions. Next to a reference simulation with a uniform temperature field, six different stratified temperature distributions have been generated. These distributions were based on the mean experimental temperature superimposed by a randomized variance, again derived from the experiments. The results showed that an asymmetric flame structure arises in the computed results when the temperature stratification input is used. In these predictions, first-stage ignition is advanced by 24 $\mu$s , while second-stage ignition is delayed by 11 $\mu$s. At the same time a lift-off length difference between the top and the bottom of up to 1.1 mm is observed. Furthermore, the lift-off length is less stable over time. Given the shown dependency, the temperature data is made available along with the vessel geometry data as a recommended basis for future numerical simulations.

42 ENGINEERING↗
DOE OSTI2020

Materials Development and Advanced Process Simulation for Additive Manufacturing with Fiber-Reinforced Thermoplastics (Final Technical Report)

DuPont, Local Motors, and the Composites Manufacturing and Simulation Center (CMSC) at Purdue worked together to advance process simulation capabilities for Additive Manufacturing (AM) and develop printable short fiber reinforced semi-crystalline engineering thermoplastics suitable for vehicle components. DuPont provided one candidate thermoplastic material, which was evaluated against the following three criteria: 1) Printability using extrusion additive manufacturing processes. 2) Suitability for vehicle applications based on mechanical coupon screening tests. 3) Predictive capability of process simulations. Indiana Manufacturing Institute provided modeling and simulation technology and access to Purdue’s lab scale AM printer. Local Motors provided mechanical property requirements and printed demonstration prototype components using their Thermwood Large Scale Additive Manufacturing (LSAM) printer. To meet the performance and cost requirements of the automotive industry, high fiber content semi-crystalline engineering materials are required. Printing experience and property data for a 50% by weight short glass fiber reinforced material were gathered. Further, material properties required for simulating the printing process with the glass-fiber reinforced material were characterized. Knowledge was acquired with regards to the processability and performance of this thermoplastic material. Predictions of part deformation and crystallinity level were verified against experimental observations. Any significant results other than verification of simulation, e.g., novel material identified here would be impactful.

36 MATERIALS SCIENCE↗
DOE OSTI2016

A surface ice module for wind turbine dynamic response simulation using FAST

It is a fact that developing offshore wind energy has become more and more serious worldwide in recent years. Many of the promising offshore wind farm locations are in cold regions that may have ice cover during wintertime. The challenge of possible ice loads on offshore wind turbines raises the demand of modeling capacity of dynamic wind turbine response under the joint action of ice, wind, wave, and current. The simulation software FAST is an open source computer-aided engineering (CAE) package maintained by the National Renewable Energy Laboratory. In this paper, a new module of FAST for assessing the dynamic response of offshore wind turbines subjected to ice forcing is presented. In the ice module, several models are presented which involve both prescribed forcing and coupled response. For conditions in which the ice forcing is essentially decoupled from the structural response, ice forces are established from existing models for brittle and ductile ice failure. For conditions in which the ice failure and the structural response are coupled, such as lock-in conditions, a rate-dependent ice model is described, which is developed in conjunction with a new modularization framework for FAST. In this paper, analytical ice mechanics models are presented that incorporate ice floe forcing, deformation, and failure. For lower speeds, forces slowly build until the ice strength is reached and ice fails resulting in a quasi-static condition. For intermediate speeds, the ice failure can be coupled with the structural response and resulting in coinciding periods of the ice failure and the structural response. A third regime occurs at high speeds of encounter in which brittle fracturing of the ice feature occurs in a random pattern, which results in a random vibration excitation of the structure. An example wind turbine response is simulated under ice loading of each of the presented models. This module adds to FAST the capabilities for analyzing the response of wind turbines subjected to forces resulting from ice impact on the turbine support structure. The conditions considered in this module are specifically addressed in the International Organization for Standardization (ISO) standard 19906:2010 for arctic offshore structures design consideration. Special consideration of lock-in vibrations is required due to the detrimental effects of such response with regard to fatigue and foundation/soil response. Finally, the use of FAST for transient, time domain simulation with the new ice module is well suited for such analyses.

17 WIND ENERGY↗
DOE OSTI2022

Development and verification of the pyIsoDep depletion package for nuclear thermal propulsion reactors

Nuclear thermal propulsion technology is currently being developed by the National Air and Space Administration for future manned Mars missions. There are multiple on-going activities to develop and enhance computational capabilities, single- and multi-physics, to advance the readiness level of this technology. The prediction of short-lived fission products is crucial to understand the reactivity penalty and the decay heat production during the different burnout stages of the vehicle. Generally, the solution entails coupled neutronic-depletion analysis. However, such simulations are somewhat expensive, especially if the duration of each engine burn is to be optimized. Therefore, there is a need to complement such high-fidelity techniques with time-efficient alternatives that do not require the sequential neutronic and depletion analysis. This paper details the methods and capabilities of pyIsoDep, an open-source code used for depletion analysis of nuclear thermal propulsion reactors. The package is general, in principle, and thus can be applied to any fission system, and can potentially be extended for fusion systems. The package was developed to provide a flexible and efficient alternative to pre-existing higher order depletion codes that could integrate well into mission and transient analysis of nuclear thermal propulsion reactors. Finally, this paper presents a validation and verification of pyIsoDep against analytical and Serpent-generated solutions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior

Our goal is to develop, apply and validate a holistic thermal, hydrologic, mechanical, and chemical (THMC) workflow that includes evaluation of induced seismic slip in EGS reservoirs. We will integrate experimental and modeling approaches to reduce parameter uncertainty and better predict/mitigate seismic hazard at EGS sites. Our novel approach couples 3D physics-based earthquake simulations with THMC models (THMc+E). This capability will enable improve engineering decisions at Utah-FORGE and move EGS operations toward repeatable, robust, economically viable, and socially accepted development. For example, our THMC+E models will predict circulation scenarios and related seismic hazard for a suite of flow rates and under uncertainty, thus enabling evaluation of optimal circulation strategy. Laboratory experiments will be performed to constrain key model parameters and Bayesian techniques will provide a probabilistic evaluation of parameters used in models. THMC+E simulations will enable exploration various circumstances that may hinder EGS success and develop mitigation strategies.

58 GEOSCIENCES↗