Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Numerical simulation of narrow bipolar electromagnetic pulses generated by thunderstorm discharges

Using the concept of avalanche relativistic runaway electrons (REs), we perform numerical simulations of compact intracloud discharge (CID) as a generator of powerful natural electromagnetic pulses (EMPs) in the HF-VHF range, called narrow bipolar pulses (NBPs). For several values of the field overvoltage and altitude at which the discharge develops, the numbers of seed electrons initiating the avalanche are evaluated, with which the calculated EMP characteristics are consistent with the measured NBP parameters. We note shortcomings in the hypothesis assuming participation of cosmic ray air showers in avalanche initiation. The discharge capable of generating NBPs produces REs in numbers close to those in the source of terrestrial {gamma}-ray flashes (TGFs), which can be an argument in favor of a unified NBP and TGF source.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Numerical Simulations of Geologic Storage Reservoir Management to Support Risk Mitigation Evaluation

This report provides a detailed description of a set of numerical simulations that represent reservoir behavior over time in response to different operational decision scenarios for detection of potential leakage and reduction or avoidance of leakage impact at a hypothetical geological carbon storage (GCS) site. These simulations serve as the basis for a series of GCS leakage risk forecasts that are to be developed using the National Risk Assessment Partnership’s Open-Source Integrated Assessment Model (NRAP-Open-IAM), and a demonstration of a simple decision support workflow for evaluation of mitigation strategies based on results of those system model forecasts. This risk assessment and decision support study is forthcoming. Four injection scenarios were considered: a constant rate carbon dioxide (CO 2 ) injection case (base case), a case with CO 2 injection rate adjustment, a case with early termination of injection operations, and a case with brine extraction. CO 2 injection operations were controlled to ensure that the pressure transient remains below the defined manageable reservoir fracture pressure, with consideration shown to hypothetical locations within the modeled spatial domain where the overburden was weaker and lower transient pressure increases were allowable. Additionally, a brine extraction alternative was considered as a reservoir management and risk mitigation option to reduce reservoir pressure, steer the plume away from any hypothetical geohazard such as fault as needed, and enhance storage capacity. Such operational actions contribute to risk management overtime. This study explores the potential utility of reservoir management for risk reduction at GCS sites. This study shows that the injection design may modify the time to CO 2 breakthrough at a legacy well; in particular, these results show that brine extraction can add value for mitigating risk both by delaying leakage and reducing pressure build-up. For the scenario considered, both pressure plots and pressure distributions demonstrate that pressure build-up was decreased by 3% with brine extraction. Additionally, extraction of brine afforded enhancement of CO 2 storage capacity by 5% compared to the base case. These findings suggest that brine extraction has substantial potential to steer the risk-related reservoir effects away from known geohazards (e.g., faults and legacy wells) by conducting pressure transient effects and CO 2 plume movement toward the production well. Injection rate adjustment scenarios considered in this study show potential value for managing both reservoir pressure transients and CO 2 plume behavior. Operational actions for reducing injection and/or early termination of injection (as compared to the base case), however, require careful design; tailoring both the extent and timing of injection rate adjustment over the injection and post-injection operational period must be thoroughly planned to balance maximizing storage and minimizing subsurface environmental risk. The study also gives preliminary consideration to the effectiveness that monitoring strategy may play in providing useful information to inform reservoir management decisions for risk reduction. Two types of monitoring were considered: 1) pressure build-up or pressure transient; and 2) potential leakage detection from a CO 2 mass or plume. Four hypothetical legacy wells, two plugged and two abandoned, were placed in the model domain. Monitoring along these wells was measured over time in individual stacked reservoir formations and shale formations to support risk mitigation decisions, especially operational decisions that assisted in risk reduction. These simulations will serve as the basis for a series of GCS leakage risk forecasts that are to be developed using NRAP’s Open-IAM, and demonstration of a simple decision support workflow for comparative assessment of mitigation alternatives based on results of those system model forecasts. This risk assessment and decision support study is forthcoming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Pore-Scale and Multiscale Numerical Simulation of Flow and Transport in a Laboratory-Scale Column

Pore-scale models are useful for studying relationships between fundamental processes and phenomena at larger (i.e., Darcy) scales. However, the size of domains that can be simulated with explicit pore-scale resolution is limited by computational and observational constraints. Direct numerical simulation of pore-scale flow and transport is typically performed on millimeter-scale volumes at which X-ray computed tomography (XCT), often used to characterize pore geometry, can achieve micrometer resolution. In contrast, the scale at which a continuum approximation of a porous medium is valid is usually larger, on the order of centimeters to decimeters. Furthermore, laboratory experiments that measure continuum properties are typically performed on decimeter-scale columns. At this scale, XCT resolution is coarse (tens to hundreds of micrometers) and prohibits characterization of small pores and grains. We performed simulations of pore-scale processes over a decimeter-scale volume of natural porous media with a wide range of grain sizes, and compared to results of column experiments using the same sample. Simulations were conducted using high-performance codes executed on a supercomputer. Two approaches to XCT image segmentation were evaluated, a binary (pores and solids) segmentation and a ternary segmentation that resolved a third category (porous solids with pores smaller than the imaged resolution). We used a mixed Stokes-Darcy simulation method to simulate the combination of Stokes flow in large open pores and Darcy-like flow in porous solid regions. Simulations based on the ternary segmentation provided results that were consistent with experimental observations, demonstrating our ability to successfully model pore-scale flow over a column-scale domain.

CFD↗
DOE OSTI2022

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗
DOE OSTI2016

The application of the mesh-free method in the numerical simulations of the higher-order continuum structures

This paper attempts to study the application of mesh-free method in the numerical simulations of the higher-order continuum structures. A high-order bending beam considers the effect of the third-order derivative of deflections, and can be viewed as a one-dimensional higher-order continuum structure. The moving least-squares method is used to construct the shape function with the high-order continuum property, the curvature and the third-order derivative of deflections are directly interpolated with nodal variables and the second- and third-order derivative of the shape function, and the mesh-free computational scheme is establish for beams. The coupled stress theory is introduced to describe the special constitutive response of the layered rock mass in which the bending effect of thin layer is considered. The strain and the curvature are directly interpolated with the nodal variables, and the mesh-free method is established for the layered rock mass. The good computational efficiency is achieved based on the developed mesh-free method, and some key issues are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Turbulent nonpremixed cool flames: Experimental measurements, Direct Numerical Simulation, and manifold-based combustion modeling

Turbulence, low-temperature chemistry, and their interactions in the form of turbulent cool flames are critical to understanding and improving advanced engines. Design of such engines requires tractable simulations which in turn necessitate turbulent combustion models that can account for cool flames. While manifold-based turbulent combustion models are an attractive option for hot flames, their applicability to cool flames is not yet fully understood. This is partially due to the lack of turbulent cool flame experiments, which has made model validation difficult. This work addresses these points with a combined experimental and computational investigation of turbulent nonpremixed dimethyl ether cool flames. First, a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner is developed and tested. Turbulent cool flames are studied using the formaldehyde planar laser-induced fluorescence (PLIF), acetone PLIF, and planar Rayleigh scattering techniques. The acetone PLIF signals are converted into mixture fraction values, and measurements of time-averaged temperature are derived from the Rayleigh scattering signals by taking advantage of the similarities between cool flames and unburned mixtures. Second, a Direct Numerical Simulation (DNS) of the CARAT burner is conducted. The flame is shown to be sensitive to thermal boundary conditions such that the stabilization method is dependent on a 10 K change in inlet temperature. Comparisons of first- and second-order statistics of temperature, mixture fraction, and formaldehyde between the DNS and experiment show good agreement. The validity of manifold-based turbulent combustion models in turbulent cool flames is then explored by first studying the effective Lewis numbers in the flame using a local differential diffusion parameter analysis. Turbulent cool flames are shown to have effective unity Lewis number transport even at a much lower Reynolds number than typically required for hot flames. Good agreement is shown for temperature and formaldehyde between DNS conditional means and one-dimensional nonpremixed flame solutions. Additionally, the motion of the cool flame in mixture fraction space is shown to be well described by one-dimensional nonpremixed flame solutions. These points indicate that the fundamental physics of turbulent cool flames are analogous to turbulent hot flames, implying that the same reduced-order manifold-based modeling approaches utilized for turbulent hot flames can be utilized for turbulent cool flames.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2015

Numerical simulation of photoexcited polaron states in water

We consider the dynamic polaron model of the hydrated electron state on the basis of a system of three nonlinear partial differential equations with appropriate initial and boundary conditions. A parallel numerical algorithm for the numerical solution of this system has been developed. Its effectiveness has been tested on a few multi-processor systems. A numerical simulation of the polaron states formation in water under the action of the ultraviolet range laser irradiation has been performed. The numerical results are shown to be in a reasonable agreement with experimental data and theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Thermal characterization of full-scale PCM products and numerical simulations, including hysteresis, to evaluate energy impacts in an envelope application

This article presents combined measurements of fatty acid-based organic PCM products and numerical simulations to evaluate the energy benefits of adding a PCM layer to an exterior wall. The thermal storage characteristics of the PCM were measured using a heat flow meter apparatus (HFMA). The PCM characterization is based on a recent ASTM International standard test method, ASTM C1784. The PCM samples were subjected to step changes in temperature and allowed to stabilize at each temperature. By measuring the heat absorbed or released by the PCM, the temperature-dependent enthalpy functions for melting and freezing were determined.Here, the simulations were done using a previously-validated two-dimensional (2D) wall model containing a PCM layer and incorporating the HFMA-measured enthalpy functions. The wall model was modified to include the hysteresis phenomenon observed in PCMs, which is reflected in different melting and freezing temperatures of the PCM. Simulations were done with a single enthalpy curve based on the PCM melting tests, both melting and freezing enthalpy curves, and with different degrees of hysteresis between the melting and freezing curves. Significant differences were observed between the thermal performances of the modeled wall with the PCM layer under the different scenarios.

2D wall model with PCM↗
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2016

Three-dimensional direct numerical simulation of turbulent lean premixed methane combustion with detailed kinetics

The interaction of maintained homogeneous isotropic turbulence with lean premixed methane flames is investigated using direct numerical simulation with detailed chemistry. The conditions are chosen to be close to those found in atmospheric laboratory experiments. As the Karlovitz number is increased from 1 to 36, the preheat zone becomes thickened, while the reaction zone remains largely unaffected. A negative correlation of fuel consumption with mean flame surface curvature is observed. With increasing turbulence intensity, the chemical composition in the preheat zone tends towards that of an idealised unity Lewis number flame, which we argue is the onset of the transition to distributed burning, and the response of the various chemical species is shown to fall into broad classes. Smaller-scale simulations are used to isolate the specific role of species diffusion at high turbulent intensities. Diffusion of atomic hydrogen is shown to be related to the observed curvature correlations, but does not have significant consequential impact on the thickening of the preheat zone. It is also shown that susceptibility of the preheat zone to thickening by turbulence is related to the 'global' Lewis number (the Lewis number of the deficient reactant); higher global Lewis number flames tend to be more prone to thickening.

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

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2021

Propagation of Input Uncertainties in Numerical Simulations of Laser Powder Bed Fusion

Laser powder bed fusion has the potential of redefining state-of-the-art processing and production methods, but defect formation and inconsistent build quality have limited the implementation of this process. Numerical models are widely used to study this process and predict the formation of these defects. Presently, the uncertainties of model input parameters and thermophysical properties used by these numerical simulations have not been investigated. In the present study, the uncertainty in these input parameters and material properties are quantified for laser powder bed fusion, with and without a simulated powder bed, to determine their influence on the predictive accuracy of an experimentally validated numerical model. Accounting for all possible sources of uncertainty quickly becomes computationally expensive on account of the curse of dimensionality. Uncertainty in laser absorption, solid, and liquid specific heat of the metal were found to have the largest effect on model prediction reliability with or without the use of a powder bed. Results also illustrate that accounting for these three uncertain parameters still captures the majority of model prediction uncertainty. Furthermore, the methodology of this study may be used to understand the uncertainty in as-built microstructure through propagation to microstructure prediction models, or applied under processing conditions where high Péclet numbers are observed and the thermal convection and fluid flow within the molten pool are substantial.

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

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Monitoring of chlorin-based photosensitiser localisation with dual-wavelength fluorescence imaging: numerical simulations

An approach to the estimation of chlorin-based photosensitiser localisation with dual-wavelength fluorescence imaging after topical application or intravenous injection is developed based on numerical simulations. Fluorescence imaging, which is commonly used for the monitoring of photosensitiser distribution within tissue, does not allow for its depth-resolved imaging. Chlorin-based photosensitisers feature two peaks in their absorption spectra corresponding to 402 and 662 nm, which allows for fluorescence excitation at two wavelengths with significantly different light penetration depths into biotissues. The ratio of fluorescence signals detected upon irradiation at different excitation wavelengths is shown to provide the information about photosensitiser localisation. The cases of exponential in-depth decay of the photosensitiser concentration and uniform photosensitiser distribution within the layer at a certain depth within the medium are analysed. These cases simulate topical application of photosensitiser and its intravenous injection and further accumulation in tumour, respectively. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Numerical simulations of quiet sun magnetism: On the contribution from a small-scale dynamo

We present a series of radiative MHD simulations addressing the origin and distribution of the mixed polarity magnetic field in the solar photosphere. To this end, we consider numerical simulations that cover the uppermost 2-6 Mm of the solar convection zone and we explore scales ranging from 2 km to 25 Mm. We study how the strength and distribution of the magnetic field in the photosphere and subsurface layers depend on resolution, domain size, and boundary conditions. We find that 50% of the magnetic energy at the τ = 1 level comes from fields with the less than 500 G strength and that 50% of the energy resides on scales smaller than about 100 km. While the probability distribution functions are essentially independent of resolution, properly describing the spectral energy distribution requires grid spacings of 8 km or smaller. The formation of flux concentrations in the photosphere exceeding 1 kG requires a mean vertical field strength greater than 30-40 G at τ = 1. The filling factor of kG flux concentrations increases with overall domain size as the magnetic field becomes organized by larger, longer-lived flow structures. A solution with a mean vertical field strength of around 85 G at τ = 1 requires a subsurface rms field strength increasing with depth at the same rate as the equipartition field strength. We consider this an upper limit for the quiet Sun field strength, which implies that most of the convection zone is magnetized close to the equipartition. We discuss these findings in view of recent high-resolution spectropolarimetric observations of quiet Sun magnetism.

79 ASTRONOMY AND ASTROPHYSICS↗