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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

The catalytic behavior of precisely synthesized Pt–Pd bimetallic catalysts for use as diesel oxidation catalysts

The demands of stricter diesel engine emission regulations have created challenges for current exhaust systems. With advances in low-temperature internal combustion engines and their operations, advances must also be made in vehicle exhaust catalysts. Most current diesel oxidation catalysts use heavy amounts of precious group metals (PGMs) for hydrocarbon (HC), CO, and NO oxidation. These catalysts are expensive and are most often synthesized with poor bimetallic interaction and dispersion. In this paper, the goal was to study the effect of aging on diesel emission abatement of Pt–Pd bimetallic nanoparticles precisely prepared with different morphologies: well dispersed core–shell vs. well dispersed homogeneously alloyed vs. poorly dispersed, poorly alloyed particles. Alumina and silica supports were studied. Particle morphology and dispersion were analyzed before and after hydrothermal treatments by XRD, EDX, and STEM. Reactivity as a function of aging was measured in simulated diesel engine exhaust. While carefully controlled bimetallic catalyst nanoparticle structure has a profound influence on initial or low temperature catalytic activity, the differences in behavior disappear with higher temperature aging as thermodynamic equilibrium is achieved. The metallic character of Pt-rich alumina-supported catalysts is such that behavior rather closely follows the Pt–Pd metal phase diagram. Nanoparticles disparately composed as well-dispersed core–shell (via seq-SEA), well-dispersed homogeneously alloyed (via co-SEA), and poorly dispersed, poorly alloyed (via co-DI) end up as well alloyed, large particles of almost the same size and activity. With Pd-rich systems, the oxidation of Pd also figures into the equilibrium, such that Pd-rich oxide phases appear in the high temperature forms along with alloyed metal cores. Finally, the small differences in activity after high temperature aging can be attributed to the synthesis methods, sequential SEA and co-DI which give rise, after aging, to a bimetallic surface enriched in Pd.

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

Guide to Using Sierra

Sierra is an engineering mechanics simulation code suite supporting the Nation's Nuclear Weapons mission as well as other customers. It has explicit ties to Sandia National Labs' workfow, including geometry and meshing, design and optimization, and visualization. Distinguishing strengths include "application aware" development, scalability, SQA and V&V, multiple scales, and multi-physics coupling. This document is intended to help new and existing users of Sierra as a user manual and troubleshooting guide.

42 ENGINEERING↗
DOE OSTI2018

Milestone Report: Grain growth and fission gas behaviour in doped UO 2 M2MS-18LA0201034

Here, the microstructure of UO 2 can be modified during nuclear fuel fabrication by using additives. Several dopants (e.g. Cr, Ti, V, Mg, Nb) are used to enhance grain growth and densification during sintering. In previous work, we used atomic scale simulation techniques to identify a common interstitial solution mechanism for a range of dopants that form positively charged defects at sintering temperatures. As a result, negatively/positively charged defect concentrations were enhanced/suppressed relative to undoped fuel under the same conditions (temperature and oxygen partial pressure). Importantly U and fission gas migration is vacancy mediated, such that the increase in the concentration uranium vacancies (which have negative charge), due to doping, causes an increase in U and fission gas diffusivity. This is important because U diffusivity is the rate limiting step in mass transport during sintering and bulk fission gas diffusivity is the underlying property that governs fission gas release. In this work, we implement analytical expressions for uranium vacancy concentrations in doped fuel (based on atomic scale results) into mesoscale and engineering scale simulations of sintering and fission gas release. It is shown that densification is dramatically enhanced at 1850 K and 1950 K due to dopant solution and that Mn-doping has a greater effect than the more widely used Cr-doping. The effect of both dopants is strongly temperature dependent and is negligible below approximately 1800 K. The effect of enlarged grain size and enhanced fission gas diffusivity has been included in BISON. Fuel temperatures under normal operating conditions are not sufficient to activate the enhanced fission gas diffusivity mechanism predicted by atomic scale calculations. Therefore, the benefits of enlarged grains for reduced fission gas release should not undermined by enhanced fission gas diffusivity for doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Modernization and optimization of a legacy open-source CFD code for high-performance computing architectures

Legacy codes remain a crucial element of today's simulation-based engineering ecosystem due to the extensive validation process and investment in such software. The rapid evolution of high-performance computing architectures necessitates the modernization of these codes. One approach to modernization is a complete overhaul of the code. However, this could require extensive investments, such as rewriting in modern languages, new data constructs, etc., which will necessitate systematic verification and validation to re-establish the credibility of the computational models. The current study advocates using a more incremental approach and is a culmination of several modernization efforts of the legacy code MFIX, which is an open-source computational fluid dynamics code that has evolved over several decades, widely used in multiphase flows and still being developed by the National Energy Technology Laboratory. Two different modernization approaches,‘bottom-up’ and ‘top-down’, are illustrated. Here, preliminary results show up to 8.5x improvement at the selected kernel level with the first approach, and up to 50% improvement in total simulated time with the latter were achieved for the demonstration cases and target HPC systems employed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Effect of Pressure on NO x Entitlement and Reaction Timescales in a Premixed Axial Jet-In-Crossflow

This paper investigates the pressure dependency of a lean premixed jet injected into a lean vitiated crossflow with an experimentally verified detailed chemistry computational fluid dynamics (CFD) model and 53 species considered. Experimental data were taken in an axially staged combustor with an optically accessible test section, allowing the use of particle image velocimetry (PIV) and CH* chemiluminescence techniques as well as point measurement of species concentration, temperature, and pressure. The experimental data cases at one, three, and five atmospheres were selected to describe the flame stabilization dependency on pressure and gain the required knowledge for an extrapolation to engine condition. Simulated exit nitrogen oxide levels were validated with experimental emission data, and a global emission trend for the NO reduction at elevated pressure and constant turbine inlet temperature level was defined. The nitrogen oxide benefit at elevated operating pressure was justified with the significantly smaller flame surface area: the analysis of the simulated spanwise and top-view profiles showed a relatively short receded core flame with nitrogen oxide production in the center at high pressure relative to a longer and larger shear layer flame at atmospheric condition that produced NO toward the inner and outer side of the flame. Furthermore, decomposition of the Damköhler number revealed the strong influence of the reaction timescales with higher reaction rates at elevated pressure, along with a moderate influence of the turbulent timescales, showing higher turbulence intensity in the lee-side recirculation zone at lower pressure.

03 NATURAL GAS↗
DOE OSTI2022

Analytical homogenization techniques applied to the Fickian diffusion: Effective diffusivity coefficient

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. In this work, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing TRISO particles. An extensive survey on effective thermal conductivity modeling was performed in our previous study, considering a random distribution of mono-sized spherical inclusions in a continuous matrix. Using the analogy between heat conduction and the simplified Fickian diffusion (or fission product species conservation), we can use the same analytical homogenization methods to obtain ETC as for the effective diffusivity coefficient (EDC). We performed several numerical experiments at varying conditions to assess the validity of our hypothesis for EDC calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Precipitation of calcium carbonate and calcium phosphate under diffusion controlled mixing

Multi-component mineral precipitation in porous, subsurface environments is challenging to simulate or engineer when in situ reactant mixing is controlled by diffusion. In contrast to well-mixed systems, the conditions that favor mineral precipitation in porous media are distributed along chemical gradients, which evolve spatially due to concurrent mineral precipitation and modification of solute transport in the media. The resulting physical and chemical characteristics of a mixing/precipitation zone are a consequence of coupling between transport and chemical processes, and the distinctive properties of individual chemical systems. We examined the spatial distribution of precipitates formed in “double diffusion” columns for two chemical systems, calcium carbonate and calcium phosphate. Polyacrylamide hydrogel was used as a low permeability, high porosity medium to maximize diffusive mixing and minimize pressure- and density-driven flow between reactant solutions. In the calcium phosphate system, multiple, visually dense and narrow bands of precipitates were observed that were reminiscent of previously reported Liesegang patterns. In the calcium carbonate system, wider precipitation zones characterized by more sparse distributions of precipitates and a more open channel structure were observed. In both cases, formation of precipitates inhibited, but did not necessarily eliminate, continued transport and mixing of the reactants. A reactive transport model with fully implicit coupling between diffusion, chemical speciation and precipitation kinetics, but where explicit details of nucleation processes were neglected, was able to qualitatively simulate properties of the precipitation zones. The results help to illustrate how changes in the physical properties of a precipitation zone depend on coupling between diffusion-controlled reactant mixing and chemistry-specific details of precipitation kinetics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

PREDICTION OF CYCLIC VARIABILITY AND KNOCK-LIMITED SPARK ADVANCE (KLSA) IN SPARK-IGNITION (SI) ENGINE

Engine knock remains one of the major barriers to further improve the thermal efficiency of spark-ignition (SI) engines. SI engine is usually operated at knock-limited spark advance (KLSA) to achieve possibly maximum efficiency with given engine hardware and fuel properties. Co-optimization of fuels and engines is promising to improve engine efficiency, and predictive computational fluid dynamics (CFD) models can be used to facilitate this process. However, cyclic variability of SI engine demands that multi cycle results are required to capture the extreme conditions. In addition, Mach Courant Friedrichs-Lewy (CFL) number of 1 is desired to accurately predict the knock intensity (KI), resulting in unaffordable computational cost. In this study, a new approach to numerically predict KLSA using large Mach CFL of 50 with ten consecutive cycle simulation is proposed. This approach is validated against the experimental data for a boosted SI engine at multiple loads and spark timings with good agreements in terms of cylinder pressure, combustion phasing, and cyclic variation. Engine knock is predicted with early spark timing, indicated by significant pressure oscillation and end-gas heat release. Maximum amplitude of pressure oscillation analysis is performed to quantify the KI, and the slope change point in KI extrema is used to indicate the KLSA accurately. Using a smaller Mach CFL number of 5 also results in the same conclusions, thus demonstrating that this approach is insensitive to the Mach CFL number. The use of large Mach CFL number allows us to achieve fast turn-around time for multicycle engine CFD simulations.

CCV↗
DOE OSTI2021

NEAMS-Multiphysics Technical Assistance in FY21

The Multiphysics Object Oriented Simulation Environment (MOOSE) [1] is a massively parallel finite-element/volume package for multiphysics simulation in science and engineering. The package focuses on providing rapid-development capabilities for engineering applications by leveraging well-built features from libMesh [2] and the Portable Extensible Toolkit for Scientific Computation (PETSc) [3]. Fiscal year 2021 (FY-21) was the first year with funding dedicated to supporting MOOSE-derived applications relevant to the Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program. In this report we outline the work done to support NEAMS applications such as BISON, Griffin, Pronghorn, and System Analysis Module (SAM).

42 ENGINEERING↗
DOE OSTI2015

Optimal Ge/SiGe nanofin geometries for hole mobility enhancement: Technology limit from atomic simulations

We use first principles simulations to engineer Ge nanofins for maximum hole mobility by controlling strain tri-axially through nano-patterning. Large-scale molecular dynamics predict fully relaxed, atomic structures for experimentally achievable nanofins, and orthogonal tight binding is used to obtain the corresponding electronic structure. Hole transport properties are then obtained via a linearized Boltzmann formalism. This approach explicitly accounts for free surfaces and associated strain relaxation as well as strain gradients which are critical for quantitative predictions in nanoscale structures. We show that the transverse strain relaxation resulting from the reduction in the aspect ratio of the fins leads to a significant enhancement in phonon limited hole mobility (7× over unstrained, bulk Ge, and 3.5× over biaxially strained Ge). Maximum enhancement is achieved by reducing the width to be approximately 1.5 times the height and further reduction in width does not result in additional gains. These results indicate significant room for improvement over current-generation Ge nanofins, provide geometrical guidelines to design optimized geometries and insight into the physics behind the significant mobility enhancement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

General Correlations for Laminar Flow Friction Loss and Heat Transfer in Plain Rectangular Plate-Fin Cores

In this paper, new generalized correlations for predicting the average fanning friction factor f and average Nusselt number Nu for laminar flow in plain plate-fin compact cores of rectangular cross section are presented. These are based on extended experimental data, as well as three-dimensional computational simulations, obtained for a broad range of fin density and geometrical attributes. The results indicate that while the fully developed forced convection scales only with the interfin channel cross-sectional ratio α (fin spacing by fin height), the entrance region hydrodynamic and thermal performance is additionally a function of the fin-core length L, flow Reynolds number Re, and fluid Prandtl number Pr. The developing flow and convection is further shown to scale as: (fRe)~(L/d h Re) 1/2 , and Nu ~(L/d h Re) 1/2 Pr 1/3 Φ(α), where f, Re, and Nu are all based on the hydraulic diameter dh of the interfin flow channel. Generalized correlations for both (fRe) and Nu are developed by the corresponding scaling of the forced convection behavior and asymptotic matching of the entrance or developing flow (short fin-core flow length) and the fully developed flow (large fin-core flow length) region performance. Finally, the predictions from these correlations are found to be within ±15% of all available experimental data for air, water, and glycol (0.71 ≤ Pr ≤ 10), and fin cores with 0 < α ≤ 1.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Modeling the Resilience of Nuclear Power Plant for Unexpected Reactor Trip

Resilience can be defined as the intrinsic ability of a system to adjust to its functioning prior to, during, or following changes and disturbances, so that it can sustain required operations under both expected and unexpected conditions. Resilience engineering is a relatively new paradigm for safety management that focuses on how to cope with complexity under pressure or disturbance to achieve success. A resilient system combines inherent human abilities with engineered safety features and organizational features in a controlled and flexible manner both in expected and in unexpected situations. The concept of resilience has been suggested to complement conventional safety analyses. Conventional safety analysis methods such as PSA have several limitations: 1) they primarily focus on technical dimension, 2) the analysis are linear and sequential, 3) they are dominated by static models, 4) they do not take a systemic view into account, and 5) they focus primarily on why accidents happen and not how success is achieved. The concept of resilience provides a more integrated view for the safety analysis in safety-critical systems, e.g., nuclear power plants. The concept of resilience in the organization of NPPs has been highlighted since the Fukushima Daiichi Nuclear Power Plant accident. An IAEA report addresses that a resilient organization is one that quickly realizes deviation from normal operations and has the ability to make even the toughest and least popular decisions and to manage the margins in which it can maneuver. While traditional strategy to safety in NPPs is to identify what could go or has gone wrong, the concept of resilience is focused on what the organization does well and what it does to maintain its successful operation and to capitalize on those processes in the event of an unexpected situation. This study attempts to model the resilience of NPP dynamically, based on the operating experiences in Korean NPPs. First, this study revised the characterization of NPP. Then, a survey on the literature was conducted to identify the influences between the elements in the model of resilience in NPPs. These works provide a qualitative model for the resilience in NPPs. Then, this study reviewed 168 licensee events reports submitted to the Korean regulatory body since 2004 and quantified the relation between the attributes in the model by performing a Pearson correlation analysis. Finally, a few simulations using the System Dynamics were conducted to discuss which factors are most influencing the resilience. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

HELICS-HLA (Hierarchical Engine for Large-scale Infrastructure Co-Simulation High Level Architecture) [SWR-18-46]

The Hierarchical Engine for Large-scale Infrastructure Co-Simulation (HELICS) is a co-simulation platform that is currently under development as part of a GMLC project. The High Level Architecture (HLA) is a co-simulation standard, various implementations of which have been widely used for co-simulations for several years, particularly in the defense community. The objective of developing the HELICS-HLA interface is to provide the ability of making software connected with HLA implementations to work with software connected with HELICS. Since each HLA federation is dependent on the RTI being used, the HELICS-HLA interface developed here uses an implementation of the HLA, the Universal CPS Environment for Federation (UCEF), to create this interface. Both HELICS and UCEF are open source software. UCEF-gateway that enables the integration of HELICS and UCEF is also an open source software.

Jain, Himanshu↗
DOE OSTI2017

TC2015: Life cycle analysis of co-formed coal fines and hydrochar produced in twin-screw extruder (TSE)

Life cycle analysis (LCA) is a quantitative tool that evaluates the environmental performance of a process or system. The main objective of this LCA study is to assess the greenhouse gas (GHG) emissions and life cycle energy use associated with solid fuels produced by co-forming coal fines and hydrochar. The hydrochar is derived from loblolly pine wood via hydrothermal carbonization (HTC) conducted in a reactive twin-screw extruder (TSE). The final product is a drop-in replacement for coal. The life cycle impacts of the co-formed end products are quantified from ‘cradle-to-grave’. Model input parameters are collected from life cycle databases, literature, experimental measurements and simulations by engineering software. Results show that electricity generated from the co-formed products has significantly lower GHG intensity and slightly higher life cycle energy use than coal-generated electricity. The most carbon-intensive component in the overall hydrochar system is the HTC plant. However, this could be improved in the future by plant re-design to recover and recycle thermal energy for more efficient operation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Partial relay selection in mixed RF/FSO dual-hop system over unified M-distributed fading channel with non-zero boresight pointing errors

In this paper, the performance of a partial relay selection-based amplify-and-forward mixed radio frequency/free space optical (RF/FSO) system with a direct RF link and outdated channel state information (CSI) is presented. The source-to-relay and source-to-destination RF links are subjected to Rayleigh distributions while the relay-to-destination FSO link experiences Malaga (M)-distribution with generalized non-zero boresight pointing errors. It is also assumed that selection combining diversity with heterodyne detection is considered for the signal reception at the destination. Moreover, the cumulative distribution functions of the end-to-end signal-to-noise under the influence of with/or without a direct link are then derived for the concerned system. Through this, the closed-form expressions for the system outage probability and average bit error rate are obtained. Thus, the influence of atmospheric turbulence, outdated CSI and relay selection are evaluated for the system under zero boresight and non-zero boresight displacement pointing errors scenarios. The results illustrated that the mixed RF/FSO system with a direct link offers better performance than the one without a direct link. The accuracy of the presented results is validated by Monte-Carlo simulation.

42 ENGINEERING↗
DOE OSTI2025

GPU-accelerated kinetic Hall thruster simulations in WarpX

Abstract Two-dimensional (axial-azimuthal) simulations of a Hall thruster are performed using the open-source particle-in-cell code WarpX. The simulation conditions are chosen to match those of the axial-azimuthal benchmark first reported by Charoy et al. in 2019. A range of numerical and solver parameters is investigated in order to find those which yield the best performance. It is found that WarpX completes the benchmark case in 3.8 days on an Nvidia V100 GPU, and in as low as 1.5 days on a more recent Nvidia H100 GPU. Of the numerical parameters investigated, it is determined that the field-solve tolerance and particle resampling thresholds have the largest effect on the simulation wall time and that particle resampling may artificially widen electron velocity distribution functions, leading to unphysical heating. A semi-implicit scheme for the electrostatic field solve is tested and is found to produce results consistent to within 10% of the benchmark in less than twelve hours. The scaling properties of the electrostatic solver to multiple GPUs are also assessed on a uniform plasma test problem. The results of this work are discussed in the context of advancements in GPU hardware and the suitability of kinetic Hall thruster simulations for engineering applications.

Marks, Thomas A.↗