SEARCH · Search NASA
Results for “Numerical”
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.
An experimental and numerical investigation of HD diesel engine DOC efficiency in oxidizing NO to NO 2
Reducing pollutant emissions from heavy-duty (HD) diesel engines is critical due to their significant environmental impact, particularly concerning NOx emissions. Understanding and optimizing modern diesel oxidation catalyst (DOC) and selective catalytic reduction (SCR) performance is essential for improving exhaust aftertreatment (EAT) system efficiency to meet stringent emissions regulations. The oxidation of nitric oxide (NO) to nitrogen dioxide (NO 2 ) in DOC plays a key role in improving SCR efficiency in reducing NO x . This study investigates the DOC performance in oxidizing NO to NO2 and its impact on the SCR efficiency of a 2021 MY Navistar E39 HD diesel engine. The influence of engine speed, load, exhaust gas temperature, and composition on DOC efficiency is experimentally investigated. The relationships between DOC inlet temperature, oxygen availability and NO 2 /NO x ratio at the DOC inlet are examined to better understand their effects on the overall DOC efficiency. The results indicate that DOC NO oxidation efficiency is highly dependent on exhaust temperature, with optimal oxidation occurring within a specific temperature range (275-350°C). Below this threshold, the chemical reactions are kinetically limited, while at higher temperatures, thermodynamic constraints reduce the efficiency of DOC in oxidizing NO to NO 2 . The experimental data further reveal that the NO 2 /NO x ratio peaks at medium loads before declining at higher loads due to reduced residence time and mass transfer effects. Additionally, the SCR NO x conversion efficiency is significantly influenced by the NO 2 /NO x ratio, achieving peak performance when the NO 2 /NO x ratio approaches 0.5. A DOC chemistry model was developed and validated against the experimental data to predict DOC oxidation behavior under various operating conditions. The findings of this study provide insights into the interdependencies between DOC and SCR performance, contributing to the optimization of SCR systems for optimized NOx reduction.
Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning
We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.
Numerical simulation of asteroid geometry variance on airburst threat
For an atmospheric airburst the primary source of concern when assessing uncertainty is the size and velocity. Determining these properties provides the basis for threat assessment, as the total energy of the asteroid may then be estimated, and the threat investigated thoroughly. Even with clarity as to how much energy an asteroid may deposit, a great deal of uncertainty still exists for the actual energy deposition process. One such source of uncertainty is the geometry of the incoming asteroid. The geometry of an asteroid will alter the stress distribution during entry, which adds uncertainty to when fracture will occur. Here, in this study, we use Smoothed Particle Hydrodynamics to model the atmospheric airburst of Tunguska-scale asteroids with varying geometric profiles, including a sphere, ellipsoid, binary and superellipsoid. Each asteroid is modeled as a homogenous structure with strength. We assess uncertainty through a series of planar 2D simulation cases for each geometry, comparing the source of stochasticity across geometries. A single 3D airburst simulation for each geometry is also analyzed. Additionally, the 3D cases are compared to the highly uncertain Tunguska event, predicting variance in burst height across geometries, but all bounded by theoretical burst heights proposed for Tunguska.
Numerical assessment of triply periodic minimal surfaces for direct air capture of carbon dioxide
Direct air capture (DAC) systems often consist of packing material wetted by a capture fluid that reacts with CO 2 in the airstream. The efficiency of the contactor is determined by a complex relationship of fluid dynamics, heat and mass transfer, contactor geometry, and chemical properties. The efficiency of the contactor must be balanced with other factors, primarily pressure drop through the system. Triply periodic minimal surfaces (TPMS) are a class of differential surfaces that have been explored in multiple engineering applications and have been shown to exhibit excellent performance when used in heat exchangers. Their tortuous path provides a high surface-to-volume ratio and favorable trade-off between contact area and pressure drop. In this work, a gyroid-type TPMS contactor was evaluated using computational fluid dynamics for a variety of geometric parameters to explore the potential benefit of TPMS shapes for DAC applications. A thin-film model was employed to model the flow and distribution of the capture solvent, allowing efficient simulations of TPMS structures at scale by eliminating the need for a computationally intensive interface capturing method. A liquid-gas mass transfer model was implemented in the commercial software STAR-CCM+ and used to predict the CO 2 capture efficiency and study the trade-off between capture performance and pressure drop through analysis of capture rates, mass transfer coefficients, and other relevant variables. TPMS contactors with a variety of geometric parameters and two capture solvent options were investigated to determine the effect of design choices on the operational performance of DAC systems. In conclusion, results showed that while contactor geometry is the dominant factor in efficiency and pressure drop, the physiochemical properties of the solvent are an important secondary influence on the contactor performance.
Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers
Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.
Experimental and numerical study of soot formation in hydrocarbon sprays under high-pressure fuel pyrolysis conditions
Not Available
Numerical methods and improvements for simulating quasi-static elastoplastic materials
Not Available
Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade
Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.
Numerical analysis of rotating scans’ effect on surface roughness in laser-powder bed fusion
Not Available
Numerical and experimental analysis of mechanically induced failure in electric vehicle battery modules
Mitigating thermal runaway and cell-to-cell propagation is essential for improving the safety of electric and hybrid vehicles. Enhancing digital twin capabilities to predict battery mechanical abuse is particularly critical for automotive and aerospace applications, where crashworthiness is a key concern. Understanding failure conditions and propagation in battery modules during mechanical abuse is complex due to interactions between structural deformation, heat transfer, electrochemical processes, exothermic reactions and mechanical fracture. While prior studies have focused on modeling cell-level behavior, extending these models to module or pack level is necessary for a system level understating of electric vehicle safety. This study develops coupled large deformation finite element models that simultaneously solve for electrochemistry, material failure, internal short circuit and thermal runaway propagation. The models account for mechanical and thermal interactions between lithium-ion cells and other battery components while the contact interfaces are evolving with time. Model-predicted voltage, temperature and force responses are compared with experimental data for validation. The results demonstrate that the approach captures key failure mechanisms, including thermal propagation through heat transfer, electrical propagation from short circuits in parallel-connected cells, and mechanical propagation via penetration and crack formation. These findings show that computational models are valuable tools for understanding battery module failure and providing insight that can reduce the need for extensive experimental testing.
Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation
We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.
Numerical implementation of the MOOSE subchannel module (SCM) algorithm
The subchannel module of moose (SCM), previously referred to as Pronghorn-SC, is a subchannel code designed to resolve single-phase flow fields and calculate the relevant flow variables, in nuclear reactor fuel-pin assemblies. The assemblies it models are: water-cooled, with bare fuel-pins in a square ducted, quadrilateral lattice and liquid-metal/water-cooled, with bare/wire-wrapped fuel pins in a hexagonal ducted, triangular lattice. Previous publications have presented the development, validation and verification of SCM. This work presents an comparative overview of the different solvers implemented within SCM. This includes the solver algorithm description, and a performance comparison. The two test cases chosen to demonstrate the solver performance were taken from the PSBT enthalpy mixing benchmark and the ORNL-19 pin benchmark.
Analytical and numerical investigation of gravity anchors for floating photovoltaic systems
Not Available
High resolution numerical simulations of methane pool fires using adaptive mesh refinement
The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.
Numerical investigations of deep learning-assisted delamination characterization using ultrasonic guided waves
Not Available