Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Effects of air chemistry and stiffened EOS of air in numerical simulations of bubble collapse in water

In this paper we study the effects of a stiffened gas equation of state for air at high pressure and air chemistry with dissociation in numerical simulations of bubble collapse in water. Two types of bubble collapse are studied, with one corresponding to spherically symmetric collapse simulating single-bubble sonoluminescence (SBSL), and the other a GPa-shock-induced axisymmetric bubble collapse. The numerical method is a finite-volume based solver with diffuse material interface model. Verification and validation of the solver are demonstrated by comparing to analytical solutions and experimental observations. We find that for sonoluminescence, air chemistry and equation of state for air at high pressures can have significant effects on the peak temperature and pressure attained during the evolution, and that the peak temperature is on the order of 1 eV, which is close to that observed in experiments. For shock-induced bubble collapse, the pressure inside the bubble will not be as high as that in sonoluminescence due to the absence of spherical symmetry. However, the air chemistry still has a significant effect on the temperature of the bubble. We also examine the effect of multiple bubble interaction in shock-induced bubble collapse and find several mechanisms that can significantly increase the pressure in the water. These mechanisms include frontal-distal-side-collision, shock focusing, upstream-traveling shocks, and compression of the water near the centerline by the vortex generated by bubble collapse. These have important implications for cavitation erosion.

42 ENGINEERING↗
DOE OSTI2021

Grid resolution requirement for resolving rare and high intensity wall-shear stress events in direct numerical simulations

Turbulent signals are intermittent with large instantaneous fluctuations. Such large fluctuations lead to small Kolmogorov scales that are hard to resolve in numerical simulations [P. K. Yeung, K. R. Sreenivasan, and S. B. Pope, Effects of finite spatial and temporal resolution in direct numerical simulations of incompressible isotropic turbulence, Phys. Rev. Fluids 3, 064603 (2018)]. Here, the present paper follows the above basic logic, but instead of dissipation events in isotropic turbulence, we study wall-shear stress events in plane channel flow. Wall-shear stress fluctuations are increasingly more intermittent as the Reynolds number increases. Hence, one has to employ higher grid resolutions as the Reynolds number increases in order to resolve a given percentage of wall-shear stress events. The objective of this paper is to quantify effects of the grid resolutions on the rare and high intensity wall-shear stress events. We find that the standard grid resolution resolves about 99% of the wall-shear stress events at Re τ = 180. A slightly higher grid resolution has to be employed in order to resolve 99% of the wall-shear stress events at higher Reynolds numbers, and if the standard grid resolution is used for, e.g., a Re τ = 10000 channel flow, one resolves about 90%–95% wall-shear stress events.

42 ENGINEERING↗
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 OSTI2019

Numerical Oscillation Prevention for PV Inverter Controllers in Quasi-Steady-State Simulators

Droop-based PV inverter control methods are popular for preventing overvoltage in low voltage networks with high installations of PV solar. Quasi-steady-state simulators (QSTS) can be used to reduce the computational burden of conducting long-term techno-economic analyses of such PV inverter controllers. If dynamic models for droop-based controllers are directly implemented in QSTS, numerical oscillation occurs and the controllers fail to converge to their steady-state operation. In this paper, a power network sensitivity-based algorithm is proposed for preventing numerical oscillation while implementing active and reactive droop-based PV inverter controllers in quasisteady state simulation tool (e.g. in GridLAB-D). Convergence of the sensitivity-based algorithm is verified using a typical radial distribution feeder. It is found that all three controllers reached converged state in less than 15 iterations.

GridLAB-D, PV inverter control, quasi-steady-state↗
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 OSTI2016

A domain-specific compiler for a parallel multiresolution adaptive numerical simulation environment

This paper describes the design and implementation of a layered domain-specific compiler to support MADNESS---Multiresolution ADaptive Numerical Environment for Scientific Simulation. MADNESS is a high-level software environment for the solution of integral and differential equations in many dimensions, using adaptive and fast harmonic analysis methods with guaranteed precision. MADNESS uses k-d trees to represent spatial functions and implements operators like addition, multiplication, differentiation, and integration on the numerical representation of functions. The MADNESS runtime system provides global namespace support and a task-based execution model including futures. MADNESS is currently deployed on massively parallel supercomputers and has enabled many science advances. Due to the highly irregular and statically unpredictable structure of the k-d trees representing the spatial functions encountered in MADNESS applications, only purely runtime approaches to optimization have previously been implemented in the MADNESS framework. This paper describes a layered domain-specific compiler developed to address some performance bottlenecks in MADNESS. The newly developed static compile-time optimizations, in conjunction with the MADNESS runtime support, enable significant performance improvement for the MADNESS framework.

Rajbhandari, Samyam↗
DOE OSTI2015

Numerical—experimental analysis of two floor platform designs for auxiliary construction truck

In an auxiliary construction truck, the permanent deformations created by impact of rocks against the floor platform during the daily work generate the premature failure and need for replacement of the platform floor. This article shows the process followed in the analysis of two floor platforms for auxiliary construction trucks. This analysis has been carried out by means of the methodology developed in collaboration with truck manufacturer (ZAMARBU S.L) and the research group “New Technologies Applied on Vehicles and Road Safety” (VEHIVIAL) of the University of Zaragoza. The methodology applied in the floor platform development has been based on the application of numerical techniques by means of the finite element method (FEM) with explicit integration of dynamic equilibrium equation and the validation of the numerical results by means of experimental test. Applied load cases correspond to the impact of a sphere of 50, 150, and 250 kg. These load cases reproduce the impact caused by stones falling on the floor platform, which are the most restrictive for the floor design, because the vehicle’s rigidity and resistance against torsion and bending are provided for the chassis. The results of the analysis served to demonstrate the efficiency of the methodology developed for the optimization of the floor platform for auxiliary construction truck, getting an innovative floor platform design that improves the current models.

42 ENGINEERING↗
DOE OSTI2018

Multidimensional Numerical Simulations of Knocking Combustion in a Cooperative Fuel Research Engine

A numerical approach was developed based on multidimensional computational fluid dynamics (CFD) to predict knocking combustion in a cooperative fuel research (CFR) engine. G-equation model was employed to track the turbulent flame front and a multizone model was used to capture auto-ignition in the end-gas. Furthermore, a novel methodology was developed wherein a lookup table generated from a chemical kinetic mechanism could be employed to provide laminar flame speed as an input to the Gequation model, instead of using empirical correlations. To account for fuel chemistry effects accurately and lower the computational cost, a compact 121-species primary reference fuel (PRF) skeletal mechanism was developed from a detailed gasoline surrogate mechanism using the directed relation graph (DRG) assisted sensitivity analysis (DRGASA) reduction technique. Extensive validation of the skeletal mechanism was performed against experimental data available from the literature on both homogeneous ignition delay and laminar flame speed. The skeletal mechanism was used to generate lookup tables for laminar flame speed as a function of pressure, temperature, and equivalence ratio. The numerical model incorporating the skeletal mechanism was employed to perform simulations under research octane number (RON) and motor octane number (MON) conditions for two different PRFs. Parametric tests were conducted at different compression ratios (CR) and the predicted values of critical CR, delineating the boundary between “no knock” and “knock,” were found to be in good agreement with available experimental data. The virtual CFR engine model was, therefore, demonstrated to be capable of adequately capturing the sensitivity of knock propensity to fuel chemistry.

CFD↗
DOE OSTI2019

Experimental and Numerical Analysis of Additively Manufactured Coupons With Parallel Channels and Inline Wall Jets

Present paper developed a series of additively manufactured parallel cooling channels with streamwise wall jets inside, which could be suitable for double-wall and near-wall cooling configurations in gas turbine hot section components. The tested coupons consisted of parallel channels, each channel further divided into small chambers using several spanwise separation walls. Height of these walls was kept less than channel height, thus forming a slot with one of the end walls. Coolant entered from one side of channel and formed streamwise wall jet while crossing through the slot over to the downstream chamber. The test coupons were additively manufactured by selective laser sintering (SLS) technique using Inconel 718 alloy. Steady-state heat transfer experiments with constant wall temperature boundary condition were performed to analyze effect of pitch between subsequent slots and blockage ratio (ratio of separation wall height to channel height) on heat transfer. The channel Reynolds number ranged from 1800 to 5000. Numerical simulations were performed using ANSYS cfx solver with SST k–ω turbulence model to obtain detailed understanding of existing flow field. Experimental results showed heat transfer enhancement of up to 6.5 times that of a smooth channel for the highest blockage ratio of 0.75. Numerical results revealed complex flow field which consisted of wall jets along with impingement, separation, and recirculation zones in each chamber. For all configurations, gain in heat transfer was accompanied with high pressure drops. However, coupled with the high heat transfer, this design could lead to potential coolant savings.

Engineering↗
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 OSTI2015

Gaussian variational ansatz in the problem of anomalous sea waves: Comparison with direct numerical simulation

The nonlinear dynamics of an obliquely oriented wave packet on a sea surface is analyzed analytically and numerically for various initial parameters of the packet in relation to the problem of the so-called rogue waves. Within the Gaussian variational ansatz applied to the corresponding (1+2)-dimensional hyperbolic nonlinear Schrödinger equation (NLSE), a simplified Lagrangian system of differential equations is derived that describes the evolution of the coefficients of the real and imaginary quadratic forms appearing in the Gaussian. This model provides a semi-quantitative description of the process of nonlinear spatiotemporal focusing, which is one of the most probable mechanisms of rogue wave formation in random wave fields. The system of equations is integrated in quadratures, which allows one to better understand the qualitative differences between linear and nonlinear focusing regimes of a wave packet. Predictions of the Gaussian model are compared with the results of direct numerical simulation of fully nonlinear long-crested waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗