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

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗

Thermal-hydraulic performance of a high-temperature zigzag-channel printed circuit heat exchanger - 18563

Advanced nuclear reactors, such as high-temperature gas-cooled reactors (HTGRs) from the Generation IV nuclear systems, are designed with the capability of delivering high-pressure, high-temperature helium to power conversion units for electricity generation and industrial plants for process heat applications. The efficiency of the electricity generation and process heat applications of HTGRs is critically dependent upon an intermediate heat exchanger (IHX), which is a key component in transferring thermal energy from the primary coolant to a secondary coolant. Printed circuit heat exchangers (PCHEs) are promising to be employed in HTGR designs due to their capability for high-temperature, high-pressure applications and compactness. In the current study, a laboratory-scale zigzag-channel PCHE was fabricated and its heat transfer and pressure drop characteristics were investigated experimentally in a high-temperature helium test facility (HTHF). A commercial computational fluid dynamics (CFD) code, STARCCM+, was used to simulate the thermal-hydraulic performance of the fabricated PCHE with a simplified geometry model. Comparisons showed some differences between experimental data available for the zigzag-channel PCHE and the numerical results. Local thermal-hydraulic performance analyses indicated that fully-developed flow conditions and periodic changes were not observed, which could be attributed to the temperature variations along the flow channels, resulting in large fluid property variations. In addition, the effects of sharp-edged or rounded zigzag channel geometry were numerically evaluated. The use of rounded corners at each bend of the channels is recommended since this method can reduce the pressure drop on both the hot and cold sides of the heat exchanger while maintaining similar heat transfer capabilities. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Quantitative modeling of gate-tunable cavity exciton-polaritons in transition-metal dichalcogenides

We performed a comprehensive real-space simulation study of propagative exciton-polaritons (EPs) of tungsten disulfide (WS 2 ) inside a Fabry–Pérot cavity. Through our simulations, we were able to construct the dispersion relation and determine the propagation length (L p ) of EPs, both of which show sensitive dependence on the applied gate voltage. We found that L p increases at higher gate voltages when excited at the exciton energy but shows opposite gate dependence at the trion energy. These observations are direct evidence of the coupling between the cavity mode with excitons and trions of WS 2 . Furthermore, we demonstrate a proof-of-concept device based on the tunable cavity EPs of WS 2 , which is promising for application as photonic transistors, modulators, or electrical-to-optical converters. Furthermore, our work unveils the real-space transport properties of gate-tunable cavity EPs and paves the way for future applications of EPs in nanophotonic devices and circuits.

42 ENGINEERING↗

Multiscale-Informed Modeling of High Temperature Component Response with Uncertainty Quantification

This report summarizes a joint effort between Argonne National Laboratory, Idaho National Laboratory, and Los Alamos National Laboratory to develop and deploy constitutive models targeted at predicting the life of Grade 91 alloy components subjected to high temperature environments typical of those that structural components in advanced nuclear reactors would experience. Two distinct, but complementary constitutive modeling approaches have been taken here. The first employs a phenomenological viscoplastic model for which parameters have been calibrated based on experimental data for a wide range of Grade 91 alloy that has undergone a variety of processing. A Bayesian approach was used to derive distributions of uncertain parameters for this model based on this data set. The second approach is a reduced order model suitable for engineering-scale analysis that is based on the results of a large set of mesoscale simulations. Mesoscale models allow for the microstructure and composition of a particular alloy to be directly taken into account in the computation of the viscoplastic response, but are computationally expensive, which makes it impractical to directly call those models for the material constitutive response in an engineering-scale simulation. The reduced-order representation of the response of the underlying model used here allows for an engineering-scale model to take into account the characteristics of the underlying microstructure, while only incurring a reasonable computational expense. Both of these approaches have different strengths, and are applicable for different parts of the design/analysis process. The phenomenological models can be readily parameterized based on a set of experimental data for a given class of materials and used for scoping calculations. Once a specific material is chosen and adequately characterized, the reduced order models can accurately predict the response of that specific alloy, and because the models are based on predictive models of the underlying microstructure, they can be used to more confidently predict the response under conditions in regions where there is limited experimental data. Both of these models have been integrated in the Grizzly code, which is used here to perform proof-of-concept uncertainty quantification analyses of a simple component under prototypical conditions. The built- in stochastic analysis capabilities in the MOOSE framework that Grizzly is built on are used here to run large sets of simulations for this uncertainty quantification analysis. As would be expected, because the reduced order models are developed for a much more tightly defined alloy, they predict tighter distributions of the time to failure than the phenomenological models, which are calibrated to a broader set of data. Also important is that these simulations demonstrate that a reduced order modeling approach can be successfully deployed to propagate uncertainties from the material scale to practical engineering-scale component simulations.

42 ENGINEERING↗

Simultaneous optimization of simulated moving bed adsorption and distillation for 2,3‐butanediol recovery

Abstract A combined simulated moving bed (SMB) and distillation separation scheme is developed to recover 2,3‐butanediol (BDO) from a dilute fermentation broth. The scheme was integrated into a lignocellulosic biorefinery that produces hydrocarbon fuels from corn stover with BDO as an intermediate. BDO recovery is one of the most challenging processes in this biorefinery; and given the high associated energy duties, direct distillation is considered cost‐prohibitive. An alternative separation is SMB adsorption in nanoporous materials, which can reject 90% of the water and reduce subsequent distillation costs. Rigorous models were used to optimize the SMB and distillation simultaneously. The separation can be added to the biorefinery while keeping the projected minimum fuel selling price (MFSP) below $0.66 USD (US dollars) per liter gasoline‐equivalent ($2.50/GGE, gallon gasoline equivalent). Finally, sensitivity analyses were conducted to assess the effects of cost and lifetime of the adsorbent, titer concentration, and BDO purity.

09 BIOMASS FUELS↗

ORNL/grit

Grit is a software library designed specifically for conducting particle-based Lagrangian simulations with CPU/GPUperformance portability. This library allows researchers and engineers to perform a wide range of simulations, including multiphase flow, and more. Grit employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs.

Ge, Wenjun [Oak Ridge National Laboratory (ORNL), ↗

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗

Molten pool dynamics and humping suppression in high-speed laser welding via tailored beam configurations

High-speed laser welding is essential for increasing the production rate of fuel cell fabrication. However, when the welding speed exceeds a critical limit, humping occurs and reduces the weld quality. In this study, two tailored beam configurations, including an adjustable ring mode and a dual-beam configuration, were employed to suppress humping. Computational fluid dynamics simulations were performed to elucidate the underlying suppression mechanisms. Here, the results show that, in the adjustable ring mode, humping mitigation arises from a reduced backward cross-sectional melt flow rate and a more stable molten pool. In the dual-beam configuration, humping suppression is attributed to the deceleration of melt flow, the conduction-mode behavior of the trailing beam, and the widening of the molten pool induced by the trailing laser. Furthermore, because the dual-beam configuration directly modifies the trailing molten pool dynamics, it achieves more effective humping suppression, extending the welding speed limit to 1.50 m/s, compared with 1.00 m/s for the adjustable ring mode.

08 HYDROGEN↗

Thermal-hydraulic performance of printed circuit heat exchangers: CFD analysis with experimental validation - 18434

Results of a computational fluid dynamic (CFD) study of flow and heat transfer in a printed circuit heat exchanger (PCHE) geometry are presented. CFD results obtained from a two-plate model are compared to corresponding experimental results for the validation. This process provides the basis for further application of the CFD code to PCHE design and performance analysis in a variety of internal flow geometries. As a part of the code verification and validation (V and V) process, CFD simulation of a single semicircular straight channel under laminar isothermal conditions was also performed and compared to theoretical results. This comparison yielded good agreement with the theoretical values. The two-plate CFD model based on the experimental PCHE design overestimated the effectiveness and underestimated the pressure drop. However, it is found that the discrepancy between the CFD result and experimental data was a result of the uncertainty in the geometry of heat exchanger during the fabrication. The CFD results obtained using a slightly smaller channel diameter yielded good agreement with the experimental data. A separate investigation revealed that the average channel diameter of the OSU PCHE after the diffusion-bonding was 1.93 mm on the cold fluid side and 1.90 mm on the hot fluid side which are both smaller than the nominal design value. Consequently, the CFD code was shown to have sufficient capability to evaluate the heat exchanger thermal-hydraulic performance. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Cycle Design and Optimization of Pressure Swing Adsorption Cycles for Pre-Combustion CO 2 Capture

Novel pressure-swing adsorption (PSA) cycles were developed based on patented TDA AMS- 19 (activated carbon) adsorbent for pre-combustion CO 2 capture in integrated gasi cation combined cycle (IGCC) power plants. A variety of cycles comprising of counter-current blowdown, pressure equalization, steam purge and light product pressurization steps were designed and simulated using an in-house one dimensional detailed model. Full process opti- mization studies were performed for all cycles to evaluate their feasibility for pre-combustion CO 2 capture. The CO 2 purity and recovery Pareto fronts obtained using the multi-objective optimization were used to assess their ability to simultaneously achieve high CO 2 purity (>95%) and recovery (>90%). The cycles that achieved the purity-recovery (95%-90%) requirements were subjected to energy-productivity optimizations under the constraints of CO 2 purity and recovery. Three cycle designs were ranked in terms of lowest energy con- sumption at 95% CO 2 purities and 90% CO 2 recoveries. It was found that a 10-step cycle with three pressure equalization steps achieved a minimum energy consumption of 95.7 kWhe/tonne of CO 2 captured at a productivity of 3.3 mol CO 2 captured/m 3 adsorbent/s.

36 MATERIALS SCIENCE↗

Resolving the dynamic correlated disorder in KTa 1- x Nb x O 3

Understanding the complex temporal and spatial correlations of ions in disordered perovskite oxides is critical to rationalize their functional properties. Here, in this study, we provide new insights into the longstanding controversy regarding the off-centering of transition metal ions in the archetypal ferroelectric alloy KTa 1-x Nb x O 3 (KTN). By mapping the full energy (E) and wavevector (Q) dependence of the dynamical structure factor S(Q, E) using neutron scattering, and rationalizing our observations with atomistic simulations leveraging machine learning, we fully resolve the static vs dynamic nature of diffuse scattering sheets, as well as their composition (x) and temperature dependence. Our first-principles simulations, extended with machine-learning molecular dynamics, reproduce both inelastic neutron spectra and diffuse features, and establish how dynamically-correlated transition metal off-centerings couple to phonons, unifying local and collective viewpoints. This study sheds new light into an exemplary ferroelectric systems and shows the importance of mapping the full S(Q, E) to reveal critical spatio-temporal correlations of atomic disorder from which functional properties emerge.

42 ENGINEERING↗

A streamline coordinate analysis of a turbulent boundary layer subject to pressure gradients and curvature on the windward side of a bump

Direct numerical simulation (DNS) of a turbulent boundary layer over the Gaussian (Boeing) bump is performed. This boundary layer exhibits a series of adverse and favourable pressure gradients and convex and concave curvature effects before separating. These effects on turbulent boundary layers are characterised and compared with a lower-Reynolds-number flow over the same geometry. The momentum budgets are analysed in the streamline-aligned coordinate system upstream of the separation region. These momentum budgets allow the simplification of equations to facilitate an integral analysis. Integral-analysis-based approximations for Reynolds stresses in the inner and outer regions of the boundary layer are also formulated. The shear and wall-normal Reynolds stress profiles normalised by these approximations exhibit a better collapse compared with friction velocity and Zagarola–Smits normalisations in the strong favourable pressure gradient region and in the mild adverse pressure region that precedes it in this flow. Simplification of these Reynolds stress approximations along with results from the DNS are used to obtain semi-empirical approximations that are able to provide stress closure in terms of wall solution fields for the turbulent boundary layer under consideration.

42 ENGINEERING↗

Modeling of internal and near-nozzle flow for a GDI fuel injector

A numerical study of two-phase flow inside the nozzle holes and the issuing spray jets for a multi-hole direct injection gasoline injector has been presented in this work. The injector geometry is representative of the Spray G nozzle, an eight-hole counterbore injector, from the Engine Combustion Network (ECN). Simulations have been carried out for a fixed needle lift. Effects of turbulence, compressibility and non-condensable gases have been considered in this work. Standard k -ε turbulence model has been used to model the turbulence. Homogeneous Relaxation Model (HRM) coupled with Volume of Fluid (VOF) approach has been utilized to capture the phase change phenomena inside and outside the injector nozzle. Three different boundary conditions for the outlet domain have been imposed to examine non-flashing and evaporative, non-flashing and non-evaporative and flashing conditions. Noticeable hole-to-hole variations have been observed in terms of mass flow rates for all the holes under all the operating conditions considered in this study. Inside the nozzle holes mild cavitation-like and in the near-nozzle region flash boiling phenomena have been predicted when liquid fuel is subjected to superheated ambiance. Under favorable conditions considerable flashing has been observed in the near-nozzle regions. An enormous volume is occupied by the gasoline vapor, stantial computational cost. Volume-averaging instead of mass-averaging is observed to be more effective, especially for finer mesh resolutions.

42 ENGINEERING↗

CFD simulations of Rapid Compression Machines using detailed chemistry: Evaluation of the ‘crevice containment’ concept

The use of creviced pistons in Rapid Compression Machines (RCMs) has proven to be very effcient in making the temperature homogeneous inside the reaction chamber but has the disadvantage of inducing a mass transfer from the reaction chamber to the crevice, especially during the preliminary heat release of two- stage ignition processes. Aiming to mitigate this mass transfer, the technique of `crevice containment' (CC) has been proposed. It consists of a physical separa- tion between the reaction chamber and the crevice region that is engaged at the end of the compression, physically preventing any mass transfer between both parts of the geometry. In order to numerically assess this novel design concept across a broader range of conditions than previously investigated, reactive simu- lations using detailed chemical kinetic mechanisms are performed for n-heptane and iso-octane. For compressed temperatures outside of the NTC (negative temperature coeffcient) region, the CC approach is very effective in suppress- ing the influence of the crevice mass transfer and thus increases the validity of the widely-used 0-D model based on the adiabatic core assumption. Still, for most of the temperature cases inside the NTC region, the ignition appears to be precipitately initiated in the residual vortex region formed after the seal engagement, possibly inducing very significant differences with the 0-D model. The benefits of eliminating the post-compression crevice mass transfer appear to be counter-balanced by effects that have not been previously investigated in detail.

42 ENGINEERING↗

Summary of BISON Development and Validation Activities - NEAMS FY15 Report

The objective of the NEAMS ToolKit is to develop a “pellet-to-plant” simulation capability useful for predicting performance and safety for a broad range of nuclear reactor power systems. The NEAMS ToolKit has been organized into a Fuels Product Line (FPL) and a Reactor Product Line (RPL) and is modular in design. Within the FPL, a multiscale approach has been adopted in which simulations of fuel performance at the engineering scale are informed by material property and irradiation performance models developed from mesoscale simulations of microstructural evolution. The focus in this report is on development and validation of the engineering-scale fuel performance analysis tool within the FPL, which is BISON. This summary report contains an overview of work performed under the work package enti- tled “FY2015 NEAMS INL-Engineering Scale Fuel Performance & Interface with RPL Tools.” A first chapter identifies the specific FY-15 milestones, providing a basic description of the asso- ciated work and references to related detailed documentation. Where applicable, a representative technical result is provided. A second chapter summarizes major additional accomplishments including 1) a BISON validation paper, 2) a multiphysics coupling paper, 3) application of BI- SON to Reactivity Insertion Accident analysis and 4) use of BISON for uncertainity analysis in fuel modeling. A final chapter outlines FY-16 future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantifying the Impacts of Grain-scale Heterogeneity on Mechanical Response

Driven by the exceedingly high computational demands of simulating mechanical response in complex engineered systems with finely resolved finite element models, there is a critical need to optimally reduce the fidelity of such simulations. The minimum required fidelity is constrained by error tolerances on the simulation results, but error bounds are often impossible to obtain a priori. One such source of error is the variability of material properties within a body due to spatially non-uniform processing conditions and inherent stochasticity in material microstructure. This study seeks to quantify the effects of microstructural heterogeneity on component- and system-scale performance to aid in the choice of an appropriate material model and spatial resolution for finite element analysis.

36 MATERIALS SCIENCE↗

A cut cell method for the 3D simulation of Crookes radiometer

Devices involved in engineering applications, such as vacuum pumps or MEMS, may be made of several moving parts. This raise the issue of the simulation of rarefied gas flow around moving boundaries. We propose a simple process, known as cut cell method, to treat the motion of a solid body in the framework of the deterministic solving of a kinetic equation. Up to our knowledge, this is the first time that this approach has been used for this kind of simulations. The method is illustrated by the 2D and 3D simulations of a Crookes radiometer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗