SEARCH · Search NASA

Results for “numerics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗

DOE OSTI2020

A Digital Twin Approach to Study Additive Manufacturing Processing Using Embedded Optical Fiber Sensors and Numerical Modeling

One of the major challenges for metal-powder-based additive manufacturing is measuring and mitigating residual strain induced during the manufacturing processes. This article reports distributed fiber optic sensors embedded in Inconel alloy components as experimental means to validate numerical models of additive manufacturing process. Electroplating was used to deposit a metal protective jacket onto standard telecom single-mode fibers for strain measurements, Fiber sensors were embedded in an Inconel alloy substrate using the laser engineered net shaping (LENS) process. Here by using a Rayleigh-scattering optical frequency domain reflectometer (OFDR), temperature changes, and residual strain in the metal substrate were monitored with 5 mm spatial resolution during the LENS process. Using finite element analysis, temperature and strain profiles induced by the LENS deposition processes were also numerically studied. Discrepancies between the simulated temperature and strain profiles and those measured directly were less than 10%. Results presented in this article demonstrates a digital twin approach to fuse modeling results with distributed fiber sensor measurement data to study additive manufacturing process toward design and fabrication process optimization.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The Numerical Computation of Clear-Sky Transmittance and Reflectance of Solar Radiation for the Modeling of PV System Performance

Modeling of photovoltaics (PV) system performance requires numerical computation of plane-of-array (POA) irradiances that are influenced by atmospheric absorption and scattering. Clear-sky transmittance and reflectance of diffuse irradiance are essential for the understanding of solar radiation related to multiple reflection between land surface and clouds. This study analytically solves the radiative transfer equation for clear-sky conditions and numerically computes the transmittance and reflectance of diffuse irradiance as functions of the optical thickness of the atmosphere. The computation will be implemented in all-sky radiative transfer models to improve their computational efficiency.

14 SOLAR ENERGY↗

DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

DOE OSTI2022

Experimental and Numerical Study on the Effect of Annular Combustor Design on Thermal Uniformity Jet in Crossflow

Abstract Non-uniformity of the exit flow temperature represents one of the significant damages to gas-turbine components, particularly turbine blades. This may occur in the course of gas-turbine operation. This paper aims to provide passive techniques by modifying the combustor design rather than changing the flow parameters to improve the thermal uniformity and turbine blades to reduce thermal stresses and increase turbine blades’ life span. An acceptable agreement between the numerical and experimental results has been achieved, and the agreement includes the velocity and temperature profile. Four different angles have been tested numerically and experimentally with a maximum error of 5% at two different Reynolds numbers. Designing the outer combustor surface with a 45-deg angle bend can give a more uniform temperature distribution of 37% higher than the basic design with only a 0.5% higher pressure drop.

Energy & Fuels↗

DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗

DOE OSTI2024

Direct Numerical Simulation of Involute Channel Turbulence

A direct numerical simulation (DNS) study was performed on turbulent flow in the high flux isotope reactor involute channel geometry to develop a numerical database and determine the differences compared with a flat parallel channel. The varying channel curvature along the walls was studied for differences in mean profiles. Parameters of interest include streamwise velocity, turbulent kinetic energy (TKE), and turbulence dissipation rate, as well as Reynolds stresses and turbulence transport terms. Profile sampling was carried out at 10 locations along the span of the involute. Additional DNS studies were performed on smaller domains of comparable curvature to the involute domain: a high curvature channel (high circular), a low curvature channel (low circular), and a flat channel (flat). Here, each of these four cases was compared against each other and to other DNS studies performed on parallel flows. The results indicate that the bulk involute channel flow does not differ significantly from a flat parallel channel flow and that the curvature of the walls does not significantly alter the mean flow parameters. However, the regions of the involute channel near the side walls exhibit relatively low magnitude twin recirculation structures driven toward the side walls from the centerline of the channel, which warrants further study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗

DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis under Boosted Spark-ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel's performance irrespective of its chemical composition, was numerically investigated In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel's knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (= 98) and MON (= 90), were first generated using a non-linear regression model based on artificial neural network (ANN). Three unconventional bioblendstocks were included in the analysis: Di-isobutylene (DIB), Isobutanol and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed was tabulated for each blend as a function of pressure, temperature and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON and engine operating conditions, the TPRF-Anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (DI) expectation was shown to be caused by differences in laminar flame speed (LFS). However; it was found that relatively large fuel-specific differences in LFS (> 20%) would have to be present to cause any appreciable deviation from the 01 framework. Otherwise, RON and MON would still be robust enough to predict a fuel's knock-limited performance.

Pal, Pinaki↗

DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗

DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part I: the Linear Setting

Fully implicit Runge--Kutta (IRK) methods have many desirable properties as time integration schemes in terms of accuracy and stability, but high-order IRK methods are not commonly used in practice with numerical PDEs due to the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic preconditioning framework for solving the systems of equations that arise from IRK methods applied to linear numerical PDEs (without algebraic constraints). Additionally, this framework also naturally applies to discontinuous Galerkin discretizations in time. Under quite general assumptions on the spatial discretization that yield stable time integration, the preconditioned operator is proven to have condition number bounded by a small, order-one constant, independent of the spatial mesh and time-step size, and with only weak dependence on number of stages/polynomial order; for example, the preconditioned operator for 10th-order Gauss IRK has condition number less than two, independent of the spatial discretization and time step. The new method can be used with arbitrary existing preconditioners for backward Euler-type time-stepping schemes and is amenable to the use of three-term recursion Krylov methods when the underlying spatial discretization is symmetric. The new method is demonstrated to be effective on various high-order finite-difference and finite element discretizations of linear parabolic and hyperbolic problems, demonstrating fast, scalable solution of up to 10th-order accuracy. The new method consistently outperforms existing block preconditioning approaches, and in several cases, the new method can achieve 4th-order accuracy using Gauss integration with roughly half the number of preconditioner applications and wallclock time as required using standard diagonally IRK methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Numerical methods for fractional Fokker–Planck equation with multiplicative Marcus Lévy noises

The Fokker–Planck equation (FPE) is an important deterministic tool for investigating stochastic dynamical systems. In this paper, we consider the space-time fractional FPE driven by multiplicative Marcus Lévy noises. Efficient numerical schemes are presented to solve the equations. Stability and convergence of the methods are also discussed. We give some numerical experiments to validate our schemes, and examine the effects of parameters on solutions. Additionally, we analyze the maximal likely trajectories and the critical time for the change of the most probability location.

Mathematics↗

DOE OSTI2022

A Comprehensive Numerical and Experimental Study for the Passive Thermal Management in Battery Modules and Packs

Cooling plates in battery packs of electric vehicles play critical roles in passive thermal management systems to reduce risks of catastrophic thermal runaway. In this work, a series of numerical simulations and experiments are carried out to unveil the role of cooling plates (both between cells and a bottom plate parallel to the cell stack) on the thermal behavior of battery modules and packs under nail penetrations. First, we investigated the role of side cooling plates on the thermal runaway propagation mitigation in battery modules (1S3P) and packs (3S3P) by varying the key parameters of the side cooling plates, such as plate thicknesses, thermal contact resistances, and materials. Then, three important factors for passive thermal management systems are identified: (i) thermal mass of side cooling plates, (ii) interfacial thermal contact resistances, and (iii) the effective heat transfer coefficients at exterior surfaces. The roles of bottom cooling plates on thermal runaway propagation mitigation in 1S3P and 1S5P battery modules are numerically investigated by comparing the thermal behavior of the modules with only side cooling plates and with both side and bottom cooling plates.

25 ENERGY STORAGE↗

DOE OSTI2021

Reduced Order Modeling for Accelerating Numerical Simulations (ROMANS)

SAND2021-15064 O ROMANS, or Reduced Order Modeling for Accelerating Numerical Simulations, is an open source software for accelerating numerical simulations. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Martin, Shawn↗

DOE OSTI2024

FTTN: Feature-Targeted Testing for Numerical Properties of NVIDIA & AMD Matrix Accelerators

FTTN is a test suite to evaluate the numerical behaviors of matrix accelerators of GPUs (NVIDIA Tensor Cores and AMD Matrix Cores) in a quick and simple setting. Matrix accelerators are heavily used in today's computationally intense applications to speed up matrix multiplications. This test suite provides a comprehensive study on the numerical behaviors of these accelerators, including support for subnormals, rounding modes, extra precision bits and FMA features. Is there

Laguna Peralta, Ignacio↗

DOE OSTI2024

Efficient general method for numerically modeling laser pulse propagation, overlap, and lifetime effects in amplifiers

An efficient numerical time-dependent general method is developed to address incoherent pulse overlap and lifetime effects in laser amplifiers. The alternating propagation-population laser energetics method (APPLE) has been validated against a semi-discrete coupled rate equation numerical method (SDRE) and analytic formalisms in bounding cases. APPLE is based on decoupled rates applied to a time-dependent framework where both space-time-dependent populations and pulse energetics are consistently updated in each time step. A significant advantage of APPLE lies in its conceptual simplicity, ease of implementation, and relatively small computational cost. SDRE tracks the populations through coupled rates and uses the method of lines to discretize the hyperbolic partial differential transport equations allowing for use of ordinary differential equation solvers. With reasonably sized mesh, we report both energetic and power pulse shape relative differences on the order of one percent between the models over a large range of initial conditions.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Experimental and numerical evaluation of laser shinethrough in polystyrene at 3 ω and 4 ω laser frequencies

In laser-direct-drive inertial confinement fusion (ICF), laser shinethrough is the process associated with laser propagation through the target plastic ablator, polystyrene is often used in ICF experiments, and leads to energy deposition in the target until the critical density is established at the ablator front surface. This process causes material modifications that may change the shock propagation and lead to density and pressure perturbations at the interface between the polystyrene and cryogenic deuterium-tritium (DT) layer. This work reports on experimental observations and numerical investigations of laser shinethrough at the third and the fourth harmonic of ICF class laser systems (1 ω corresponds to a wavelength of ∼1.06 μm ). The experiments evidence that, at 4 ω , the laser is absorbed over a distance of about a few microns and laser shinethrough is completely eliminated in contrast with the behavior at 3 ω . These experimental observations are reproduced by numerical simulations. The laser fluence required to reach the critical density at the ablator front surface is also found to be more than one order of magnitude higher at the third harmonic. Finally, modeling suggests that the back surface of the ablator exhibits ion temperatures greater than the melting temperature of cryogenic DT and electron pressure larger than the bulk and shear moduli of cryogenic DT for the case of using a 3 ω laser.

Pineau, A. [Univ. of Rochester, NY (United States)↗

DOE OSTI2022

NUMERICAL INVESTIGATION OF AIR DEHUMIDIFICATION THROUGH WATER DROPLETS DIELECTROPHORESIS

This paper numerically investigated the dehumidification of moist airflow in a converging rectangular duct with electrically enhanced condensation. The aim was to develop a condensation model that predicted water vapor condensation by employing the dielectrophoresis principle. The charged droplets were injected into the computational domain to produce a gradient electric field. The dielectrophoretic interaction between the droplets and vapor molecules of the humid air resulted in local water vapor condensation at the vapor and liquid droplet interface. This phenomenon is described in detail in this paper. A hybrid Eulerian-Lagrangian solver SprayFoam was developed based on OpenFOAM® to simulate condensation in two-phase gas-liquid mixtures. The following developments were made based on the standard compressible multiphase flow solver SprayFoam in OpenFOAM®: (1) Eulerian solver for gas phase (2) multi-component species transport, and (3) Lagrangian solver for gas-droplet two-phase flows and sub-models for liquid droplets. The newly developed numerical model was experimentally validated with data from a series of tests conducted in the authors' laboratory. The results showed that the simulations followed the same trends as the data, and the model predicted the condensation due to the electro-spray injection process in the air stream. The injected water droplets, which are electrically charged to their Rayleigh limits, increased in size while wiping out the humidity from the air. The simulations indicated that the dehumidification was about 1% for 0.5 cubic feet per minute airflow rate. Scaling up to larger flows is a future follow-up work.

Yel Mahi, *Maliha↗