Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence modeling”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Aeroacoustics Noise Model of OpenFAST

The report describes theory and application of a newly released module of OpenFAST to simulate the aeroacoustic noise generated by the rotor of an arbitrary wind turbine. OpenFAST is a fully open-source publicly available wind turbine analysis tool actively developed at the National Renewable Energy Laboratory. The aeroacoustic module, which is also fully open-source and publicly available, is based on work performed over the past three decades. Frequency-based models for turbulent inflow, turbulent boundary layer – trailing edge, laminar boundary layer – vortex shedding, tip vortex, and trailing edge bluntness – vortex shedding noise mechanisms are included. A simple directivity model is also included. The appendices of this report describe in detail the inputs to the module and how to find them in the input files of OpenFAST. The noise models are exercised simulating the aeroacoustic noise emissions of the IEA Wind Task 37 land-based reference wind turbine. A code-to-code comparison between the implementation presented here and the implementation available at the Wind Energy Institute of the Technical University of Munich, Germany, is also presented.

17 WIND ENERGY↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2026

Quantitative analysis of zonal flow influence on turbulent plasmas driven by trapped electrons

The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.

Song, Jiheon [Hanyang University, Seoul (Korea, Re↗
DOE OSTI2026

Laser absorption measurements of temperature, pressure, CO, and CO 2 at near-MHz rates in post-detonation fireballs with comparison to synthetic measurements

A laser absorption spectroscopy (LAS) diagnostic was used to obtain measurements of temperature, pressure, CO, and CO 2 at 500 kHz or 1 MHz in post-detonation fireballs produced by hemispherical samples of pentaerythritol tetranitrate (PETN). A quantum-cascade laser was scanned over multiple CO absorption transitions near 2008.5 cm −1 at 1 MHz, while an interband-cascade laser was scanned over a CO 2 absorption transition near 2394.8 cm −1 at 500 kHz. Light from each laser was combined onto a single path and passed through a detonation chamber approximately 83 mm above the 12-mm diameter hemispherical PETN charge. The CO and CO 2 absorption signals were post-processed to obtain time histories of temperature, pressure, species column pressures (P CO L, P CO2 L), and species column mole fractions (X CO L, X CO2 L). Additionally, schlieren imaging was performed simultaneously at 500 kHz to aid interpretation of the LAS measurements. Experimental and synthetic (i.e., CFD based) LAS measurements were compared to evaluate the accuracy of the CFD model and its ability to model the turbulent afterburning of the detonation products in air. In general, the experimental measurements exhibit reasonable agreement with the synthetic measurements at early times; thereby supporting the accuracy of the CFD model. Periods of disagreement between experimental and synthetic measurements at later times are most likely due to a reflected shock and detonator cavity jetting, which are not accounted for in the CFD model.

Schwartz, Charles J. [Purdue Univ., West Lafayette↗
DOE OSTI2022

Coupling a Lagrangian–Eulerian Spark-Ignition (LESI) model with LES combustion models for engine simulations

In the United States transportation sector, Light-Duty Vehicles (LDVs) are the largest energy consumers and CO 2 emitters. Electrification of LDVs is posed as a potential solution, but SI engines can still contribute to decarbonization. Car manufacturers have turned to unconventional engine operation to increase the efficiency of Spark-Ignition (SI) engines and reduce the carbon emissions of their fleets. Dilute, lean, and stratified-charge engine operation has the potential for engine efficiency improvements at the expense of increased cyclic variability and combustion instability. At such demanding engine conditions, the spark ignition event is key for flame initiation and propagation and for enhanced combustion stability. Reliable and accurate spark ignition models can help design ignition systems that reduce cyclic variability. Multiple computational spark-ignition models exist that perform well under conventional conditions, but the underlying physics needs to be expanded, for unconventional engine operation. In this paper, a hybrid Lagrangian–Eulerian Spark-Ignition (LESI) model is coupled with different turbulent flame propagation models for engine simulations. LESI relies on Lagrangian arc tracking and Eulerian energy deposition. The LESI model is coupled with the Well-Stirred Reactor (WSR), Thickened Flame Model (TFM), and g-equation model and used to simulate several cycles of a Direct-Injection Spark-Ignition (DISI) engine using a commercial Computational Fluid Dynamics (CFD) engine solver. The results showcase the successful coupling of LESI with the combustion models. Global engine metrics, such as pressure and Apparent Heat Release Rate (AHRR), for each simulation setup are compared to experimental engine results, for validation. In addition, results highlight the successful prediction of spark channel movement by comparing simulation images to experimental optical engine images. Finally, the successful coupling of LESI to combustion models, making it a usable model in the engine modeling community, is emphasized and future development details are discussed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Laminar-to-turbulent flame transition and cycle-to-cycle variations in large eddy simulation of spark-ignition engines

This paper investigates the effect of laminar-to-turbulent flame transition modeling on the prediction of cycle-to-cycle variations (CCVs) in large eddy simulation (LES) of spark-ignition (SI) engines. A laminar-to-turbulent flame transition model that describes the non-equilibrium sub-filter flame speed evolution during an early stage of flame kernel growth is developed. In the present model, the flame transition is characterized by the flame kernel size at which the flame transition ends, defined here as the flame transition scale. The proposed model captures the effects that variations in a turbulent flow field have on the evolution of early-stage burning rates, through variations in the flame transition scale. The proposed flame transition model is combined with the front propagation formulation (FPF) method and a spark-ignition model to predict CCVs in a gasoline direct injection SI engine. It is found that multi-cycle LES with the proposed flame transition model reproduces experimentally-observed CCVs satisfactorily. When the transition model is not considered or when variations in the transition process are neglected, CCVs are significantly under-predicted for the case considered here. These results indicate the importance of modeling the laminar-to-turbulent flame transition and the effect of turbulence on the transition process, when predicting CCVs, under certain engine conditions. The LES results are also used to analyze sources for variations in the flame transition. It is found, for the present engine case, that the most important source is the cycle-to-cycle variation in the turbulence dissipation rate, which is used to measure the strength of turbulence in the proposed model, near a spark plug. The large-scale velocity field and the variations of the laminar flame speed due to the mixture composition and thermal stratification are also found to be important factors to contribute to the variations in the flame transition.

42 ENGINEERING↗
DOE OSTI2020

DPM: A deep learning PDE augmentation method with application to large-eddy simulation

A framework is introduced that leverages known physics to reduce overfitting in machine learning for scientific applications. The partial differential equation (PDE) that expresses the physics is augmented with a neural network that uses available data to learn a description of the corresponding unknown or unrepresented physics. Training within this combined system corrects for missing, unknown, or erroneously represented physics, including discretization errors associated with the PDE's numerical solution. For optimization of the network within the PDE, an adjoint PDE is solved to provide high-dimensional gradients, and a stochastic adjoint method (SAM) further accelerates training. Additionally, the approach is demonstrated for large-eddy simulation (LES) of turbulence. High-fidelity direct numerical simulations (DNS) of decaying isotropic turbulence provide the training data used to learn sub-filter-scale closures for the filtered Navier–Stokes equations. Out-of-sample comparisons show that the deep learning PDE method outperforms widely-used models, even for filter sizes so large that they become qualitatively incorrect. It also significantly outperforms the same neural network when a priori trained based on simple data mismatch, not accounting for the full PDE. Measures of discretization errors, which are well-known to be consequential in LES, point to the importance of the unified training formulation's design, which without modification corrects for them. For comparable accuracy, simulation runtime is significantly reduced. A relaxation of the typical discrete enforcement of the divergence-free constraint in the solver is also successful, instead allowing the DPM to approximately enforce incompressibility physics. Since the training loss function is not restricted to correspond directly to the closure to be learned, training can incorporate diverse data, including experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A passive model for the evolution of subgrid-scale instabilities in turbulent flow regimes

In this work, a zero-dimensional model for the development of turbulent mixing layers driven by Rayleigh–Taylor (RT) or Richtmyer–Meshkov (RM) instabilities is presented, for application in Reynolds-Averaged Navier–Stokes (RANS) simulations. The model approximates the behavior of the turbulent mixing layer at a passively advected tracer particle that sits on the material interface of interest in the simulation, and is used to provide initial conditions for the RANS model once the mixing layer width has developed to a size that may be resolved on the computational mesh. The model relies on approximations to the Besnard–Harlow–Rauenzahn (BHR) RANS closure equations and is thus consistent the turbulent development of the mixing layer, unlike previous models of this type which have focused on the early-time laminar development of RT or RM instabilities. These instabilities will often grow through an initially laminar state before developing into turbulence, and criteria for switching from an existing laminar model to the proposed turbulent model are considered. The proposed model reproduces the mixing layer development of RT and RM instabilities seen in mixed width resolving BHR RANS simulations to a reasonable degree of accuracy, and relaxes the need for the RANS simulation to resolve the mixed width as soon as the mixing layer becomes turbulent, which may be impractical in cases where the mixing layer grows through a wide range of scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗
DOE OSTI2021

Radiation GRMHD simulations of M87: funnel properties and prospects for gap acceleration

Here, we use the public code EBHLIGHT to carry out 3D radiative general relativistic magnetohydrodynamics (GRMHD) simulations of accretion on to the supermassive black hole in M87. The simulations self-consistently evolve a frequency-dependent Monte Carlo description of the radiation field produced by the accretion flow. We explore two limits of accumulated magnetic flux at the black hole (SANE and MAD), each coupled to several subgrid prescriptions for electron heating that are motivated by models of turbulence and magnetic reconnection. We present convergence studies for the radiation field and study its properties. We find that the near-horizon photon energy density is an order of magnitude higher than is predicted by simple isotropic estimates from the observed luminosity. The radially dependent photon momentum distribution is anisotropic and can be modeled by a set of point-sources near the equatorial plane. We draw properties of the radiation and magnetic field from the simulation and feed them into an analytic model of gap acceleration to estimate the very high energy (VHE) γ-ray luminosity from the magnetized jet funnel, assuming that a gap is able to form. We find luminosities of $\rm \sim 10^{41} \, erg \, s^{-1}$ for MAD models and $\rm \sim 2\times 10^{40} \, erg \, s^{-1}$ for SANE models, which are comparable to measurements of M87’s VHE flares. The time-dependence seen in our calculations is insufficient to explain the flaring behaviour. Our results provide a step towards bridging theoretical models of near-horizon properties seen in black hole images with the VHE activity of M87.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

The persistent shadow of the supermassive black hole of M87. II. Model comparisons and theoretical interpretations

The Event Horizon Telescope (EHT) observation of M87∗ in 2018 has revealed a ring with a diameter that is consistent with the 2017 observation. The brightest part of the ring is shifted to the southwest from the southeast. In this paper, we provide theoretical interpretations for the multi-epoch EHT observations for M87∗ by comparing a new general relativistic magnetohydrodynamics model image library with the EHT observations for M87∗ in both 2017 and 2018. The model images include aligned and tilted accretion with parameterized thermal and nonthermal synchrotron emission properties. The 2018 observation again shows that the spin vector of the M87∗ supermassive black hole is pointed away from Earth. A shift of the brightest part of the ring during the multi-epoch observations can naturally be explained by the turbulent nature of black hole accretion, which is supported by the fact that the more turbulent retrograde models can explain the multi-epoch observations better than the prograde models. The EHT data are inconsistent with the tilted models in our model image library. Assuming that the black hole spin axis and its large-scale jet direction are roughly aligned, we expect the brightest part of the ring to be most commonly observed 90 deg clockwise from the forward jet. This prediction can be statistically tested through future observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Grid resolution requirement of chemical explosive mode analysis for large eddy simulations of premixed turbulent combustion

Full Article Figures & data References Citations Metrics Reprints & Permissions Read this article Abstract The grid resolution requirement for trustworthy Chemical Explosive Mode Analysis (CEMA) in Large Eddy Simulation (LES) of premixed turbulent combustion is proposed. Explicit filtering, to emulate the effect of the LES filter, is applied to one-dimensional laminar flame and three-dimensional planar turbulent flames across a wide range of Karlovitz numbers (5 - 239). The identification of the flame front by CEMA is found relatively insensitive to the cell size (Δ), while the combustion mode identification shows more significant sensitivity. Specifically, increasing Δ falsely enhances the auto-ignition and local extinction modes and suppresses the diffusion-assisted mode. Limited dependence of the CEMA performance on the turbulent combustion regime (Karlovitz number) is observed. A simple grid size criterion for reliable CEMA mode identification in LES is proposed as Δ ≲ δ L /2; The criterion can be relaxed to Δ ≲ δ L in the laminar flame limit. Furthermore, theoretical analysis is conducted on an idealised chemistry-diffusion system. The effects of the filtering process and turbulence on the local combustion mode are demonstrated, which is consistent with the numerical observations. Further, by incorporating turbulent combustion models in CEMA, potential improvement in identifying local combustion modes can be expected.

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

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 1: Dust budget analyses and the impacts of a revised coupling scheme

Abstract. An earlier study evaluating dust life cycle in the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) has revealed that the simulated global mean dust lifetime is substantially shorter when higher vertical resolution is used, primarily due to significant strengthening of dust dry removal in source regions. This paper demonstrates that the sequential splitting of aerosol emissions, dry removal, and turbulent mixing in the model's time integration loop, especially the calculation of dry removal after surface emissions and before turbulent mixing, is the primary reason for the vertical resolution sensitivity reported in that earlier study. Based on this reasoning, we propose a revised numerical process coupling scheme that requires the least amount of code changes, in which the surface emissions are applied before turbulent mixing instead of before dry removal. The revised scheme allows newly emitted particles to be transported aloft by turbulence before being removed from the atmosphere, and hence better resembles the dust life cycle in the real world. Sensitivity experiments show that the revised process coupling substantially weakens dry removal and strengthens vertical mixing in dust source regions. It also strengthens the large-scale transport from source to non-source regions, strengthens dry removal outside the source regions, and strengthens wet removal and activation globally. In transient simulations of the years 2000–2009 conducted using 1∘ horizontal grid spacing, 72 vertical layers, and unchanged tuning parameters of emission strength, the revised process coupling leads to a 40 % increase in the global total dust burden and an increase of dust lifetime from 1.8 to 2.5 d in terms of 10-year averages. Weakened dry removal and increased mixing ratios are also seen for other aerosol species that have substantial surface emissions, although the changes in mixing ratio are considerably smaller for the submicron species than for dust and sea salt. Numerical experiments confirm that the revised coupling scheme significantly reduces the strong and non-physical sensitivities of model results to vertical resolution in the original EAMv1. This provides a motivation for adopting the revised scheme in EAM as well as for further improvements on the simple revision presented in this paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

3D Convective Urca Process in a Simmering White Dwarf

Abstract A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. In the early stages of carbon burning, called the simmering phase, energy released by the reactions in the core drive the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection, which may alter the composition and structure of the white dwarf. The convective Urca process is not well understood and requires 3D fluid simulations to properly model the turbulent convection, an inherently 3D process. Because the neutron excess of the fluid both sets and is set by the extent of the convection zone, the realistic steady state can only be determined in simulations with real 3D mixing processes. Additionally, the convection is relatively slow (Mach number less than 0.005) and thus a low Mach number method is needed to model the flow over many convective turnovers. Using the MAESTROeX low Mach number hydrodynamic software, we present the first full-star 3D simulations of the A = 23 convective Urca process, spanning hundreds of convective turnover times. Our findings on the extent of mixing across the Urca shell, the characteristic velocities of the flow, the energy-loss rates due to neutrino emission, and the structure of the convective boundary can be used to inform 1D stellar models that track the longer-timescale evolution.

Boyd, Brendan (ORCID:0000000254199751)↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗