Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stratified”

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.

360 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗
DOE OSTI2021

Quantifying wind plant blockage under stable atmospheric conditions

Wind plant blockage reduces the wind velocity upstream undermining turbine performance for the first row of the plant. We assess how atmospheric stability modifies the induction zone of a wind plant in flat terrain. We also explore different approaches to quantifying the magnitude and extent of the induction zone from field-like observations. To investigate the influence from atmospheric stability, we compare simulations of two stable boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. We find a faster cooling rate at the surface, which produces a stronger stably stratified boundary layer, amplifies the induction zone of both an isolated turbine and of a large wind plant. A statistical analysis on the hub-height wind speed field shows wind slowdowns only extend far upstream (up to 15D) of a wind plant in strong stable boundary layers. To evaluate different ways of measuring wind plant blockage from field-like observations, we consider various ways of estimating the freestream velocity upstream of the plant. Sampling a large area upstream is the most accurate approach to estimating the freestream conditions, and thus of measuring the blockage effect. Also, the choice of sampling method may induce errors of the same order as the velocity deficit in the induction zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Mountain waves can impact wind power generation

Abstract. Mountains can modify the weather downstream of the terrain. In particular, when stably stratified air ascends a mountain barrier, buoyancy perturbations develop. These perturbations can trigger mountain waves downstream of the mountains that can reach deep into the atmospheric boundary layer where wind turbines operate. Several such cases of mountain waves occurred during the Second Wind Forecast Improvement Project (WFIP2) in the Columbia River basin in the lee of the Cascade Range bounding the states of Washington and Oregon in the Pacific Northwest of the United States. Signals from the mountain waves appear in boundary layer sodar and lidar observations as well as in nacelle wind speeds and power observations from wind plants. Weather Research and Forecasting (WRF) model simulations also produce mountain waves and are compared to satellite, lidar, and sodar observations. Simulated mountain wave wavelengths and wave propagation speeds (group velocities) are analyzed using the fast Fourier transform. We found that not all mountain waves exhibit the same speed and conclude that the speed of propagation, magnitudes of wind speeds, or wavelengths are important parameters for forecasters to recognize the risk for mountain waves and associated large drops or surges in power. When analyzing wind farm power output and nacelle wind speeds, we found that even small oscillations in wind speed caused by mountain waves can induce oscillations between full-rated power of a wind farm and half of the power output, depending on the position of the mountain wave's crests and troughs. For the wind plant analyzed in this paper, mountain-wave-induced fluctuations translate to approximately 11 % of the total wind farm output being influenced by mountain waves. Oscillations in measured wind speeds agree well with WRF simulations in timing and magnitude. We conclude that mountain waves can impact wind turbine and wind farm power output and, therefore, should be considered in complex terrain when designing, building, and forecasting for wind farms.

17 WIND ENERGY↗
DOE OSTI2025

Data and code from: Multivariate bayesian regression model for predicting disposed ash composition at U.S. coal fired power stations

This dataset contains the code and data files needed for implementation of a Multivariate Bayesian Regression model, described in Jin et al. (2025), for the historical prediction of the chemical composition of disposed coal ash at U.S. coal fired power plants as a function of annualized coal purchase data. The integrated coal supply data file (CoalSupplyDataset.csv) represents a compilation of monthly fuel purchase records for the period 1973-2022 at major U.S. power stations. These records were obtained from the U.S. Energy Information Administration. The CSV file also contains, for each coal purchase record, the coal region of the mine as defined by the U.S. Geological Survey. Data entry errors and data gaps in the EIA records were corrected as described in Jin et al. This CSV file represents the integrated coal supply data after corrections were made. The model structure and fitting parameters are encoded in pickle file format (Bayesian.pkl). The model was developed with the coal supply data and coal ash composition data, apportioned according to the Stratified Shuffle Split for training and testing subsets. The model was built using Python and the PyMC library. Reference Publication: Jin, Z.; Huang, J.; Hower, J.C.; Hsu-Kim, H.(2025). Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites. Environmental Science & Technology.

Coal ash composition↗
DOE OSTI2020

Spark-ignited direct-injection engine combustion systems

A direct-injection stratified charge internal combustion engine includes a combustion cylinder to receive an air-fuel mixture, and an air intake port to inlet air into the combustion cylinder. The direct-injection engine also includes a fuel injector configured to deliver fuel within the cylinder in a spray pattern substantially aligned to a cylinder central axis to create the air-fuel mixture. A spark igniter is located within a path of the spray pattern to ignite combustion of the air-fuel mixture. The direct-injection engine further includes a movable piston defining a lower boundary of the combustion cylinder to contain the combustion of the air-fuel mixture. The piston is configured to include a bowl portion having local geometric features located on an intake port side of the combustion cylinder to redirect fluid flow towards a vortex in fluid communication with a combustion location near the cylinder central axis.

42 ENGINEERING↗
DOE OSTI2020

Viscosity of the Mercurian Magma Ocean: Implications of Sulfur-Free and Sulfur-Bearing Magma Oceans for Differentiation and Crustal Petrogenesis

Mercury is believed to have had a magma ocean in its early evolution whose solidification dynamics contributed to the planet’s present day surface composition. The surface is characterized by two main geological provinces; the Northern Smooth Plains (NSP),and the Heavily Cratered Terrain (HCT) and Intercratered Terrian (IcT). These provinces are defined by the differences in morphology, composition, and age. The NSP are located in the northern hemisphere and have low Mg/Si, Ca/Si, and S/Si ratios and high Al/Si ratios, and are compositionally analogous to a magnesian alkali basalt. The HCT/IcT cover the majority of the planet and are older than the NSP, and are characterized by high Mg/Si, Ca/Si, and S/Si ratios and low Al/Si ratios. The dichotomy in surface composition motivates the study into the structure of the mantle as a product of magma ocean solidification. Two important and interrelated parameters that determine the efficiency of crystal fractionation during magma ocean solidification are convection vigor and viscosity of the magma ocean. Crystal suspension in a cooling magma ocean could be efficient or inefficient which suggests two end member possibilities for magma ocean solidification. The first is the formation of a compositionally stratified mantle reflecting a cumulate sequence produced along a fractional liquid line of descent. The second is the formation of a well-mixed, batch crystallized mantle by high viscosity and/or vigorously convecting magma ocean. Determining which of these end members, or where on the spectrum between them is most relevant to Mercury’s magma ocean can be experimentally constrained through high P-T falling-sphere viscometry experiments.

Mouser, M. D.↗
DOE OSTI2020

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Updating Critical Temperature Calculations for Several Secondary High Explosives

In designing experiments where high explosives (HEs) are heated, it is important to have an understanding of where thermal runaway may occur. This determination is often done by using the Frank-Kamenetskii (FK) equation. For several highly-studied HEs, the parameters necessary for FK calculations are typically referenced from decades-old literature, and more recent experimental data have shown that these values require adjustment. For example, some of the legacy parameters dangerously overpredict the critical temperature by tens of degrees relative to more recent observed values. Herein, we first summarize historical results and highlight the importance of insulation effects when estimating critical temperatures for small samples, which is especially relevant to small-scale aging, compatibility, and characterization experiments. Here, we then present updated FK parameters for several common secondary CHNO HEs like RDX, HMX, PETN, and TATB, as well as the first reported values for 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105). Our updated parameters produce critical-temperature curves that are consistent with both mm-scale differentical scanning calorimetry (DSC) measurements and cm-scale data from the Lawrence Livermore National Laboratory one-dimensional time-to-explosion (ODTX) experiment. In analyzing the critical temperature versus sample size curves derived from the transcendental FK equation, we discovered that these curves are remarkably well described by a simple power-law function with a universal power of 0.142 and a HE-specific scaling factor.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

On Master-Length Scale Formulations for Stable Conditions in Turbulence Closure Models

Three formulations of the turbulence-length scales used in numerical modelling of atmospheric flows are compared. The comparison is made using the Mellor–Yamada–Nakanishi–Niino turbulence closure model within the stable boundary layer local similarity framework. With an appropriate choice of model constants, the model predictions are barely discernible and compare well with the empirical data obtained from the SHEBA campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigation of ventilation-coupled high energy density sensible thermal energy storage

Low-cost energy storage will play an important role in supporting the decarbonization of the energy sector. Here, a novel approach to thermal energy storage for buildings is proposed, in which a tank of antifreeze liquid can be used to heat outdoor air in mechanical ventilation systems. If heated to a high enough temperature, the fluid can potentially undergo temperature swings of nearly 100°C during discharge in cold climates. Because the energy density of sensible storage systems scales proportionally to the system temperature change, this concept has the potential to offer higher energy densities than other liquid-based sensible storage devices used in building applications. A one-dimensional numerical model of the storage system was developed and experimentally validated using 30 wt% (wt%) potassium acetate as an aqueous antifreeze solution. Using a commercial hot water tank, energy densities of 47.0 kW-hour per cubic meter (kWh/m 3 ) were demonstrated in a laboratory setting. Material energy densities of approximately 74.7 (kWh/m 3 ) were measured. Design improvements may boost this energy storage density even further. Because of the system's relative simplicity, ventilation-coupled sensible storage has the potential to be an easily deployable, low-cost energy storage solution for building systems.

25 ENERGY STORAGE↗
DOE OSTI2023

Relationships among diurnal variations of polar night cloud, precipitation, surface temperatures, and the fair-weather return current of the global electric circuit (GEC)

Over the past five years (June 2017-current), the vertical electric field (E z ) as well as numerous cloud, precipitation and radiation properties have been monitored at the Department of Energy-Atmospheric Radiation Measurement (DOE ARM) North Slope of Alaska (NSA) field site. Comparisons between the composite diurnal averaged fair-weather E z , and composite ceilometer derived cloud base height during the polar night, reveal a significant correlation between the parameters (r = 0.62), supporting previous studies that there is high correlation between local electric field and cloud properties. With the use of extensive instrumentation at the site, such as the Micro Pulse Lidar (MPL), Ka-band Zenith Radar (KAZR), ceilometer, SKYRAD, Precipitation Imaging Package (PIP), among others, this study provides a more comprehensive examination of the diurnal cycle of cloud and precipitation properties along with the localized fair-weather return current of the larger Global Electric Circuit (GEC) system. Comparisons between the composite diurnal averaged fair-weather E z , and cloud thickness, maximum column backscatter, and precipitation particle counts all show similar diurnal variability during the polar night, indicating that during the largest magnitude fair weather E z time periods, clouds bases tend to be higher, clouds are thicker, have a larger column backscatter, and display more precipitating particles. Furthermore, a slight diurnal variability in the polar night surface temperature was found to be highly correlated (r = 0.87) to the longwave downwelling irradiance measured by SKYRAD, indicating that the variations in the physical properties of local clouds could modulate the diurnal polar night surface temperature variability on the order of 0.5 °C/day. In conclusion, these findings emphasize the importance and global nature of the GEC system, with the global aggregate of thunderstorms and electrified clouds potentially influencing polar night cloud properties as well as diurnal wintertime polar surface temperature variation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A statistical approach to screening isotopic signatures in monitoring for underground nuclear explosions

The ability to differentiate between atmospheric radionuclide signatures from underground nuclear explosions (UNEs) and signals from other sources, such as medical isotope-production facilities and nuclear reactors, can be critical to the detection and monitoring of unannounced, low-yield nuclear events. Signatures having anomalously high amplitudes, compared to background levels, remain the best indicator in screening for a UNE. However, isotopic composition can further validate a suspected UNE signature, but separation from any atmospheric background composition is first necessary. To date, evaluating the challenges of performing this separation has typically involved comparing an observed background with a highly idealized deterministic model of radioxenon signature production by a UNE that does not consider the influence of post-detonation chemical/physical processes in the detonation cavity or the subsequent gas transport mechanisms that can also affect the isotopic composition of the detected gas signature. In addition, purely deterministic models, as previously employed, overlook the uncertainty inherent in estimating critical parameters characterizing the UNE and its detonation environment. In this paper, we create detailed, multi-parameter models of radionuclide evolution using the widely accepted England and Rider post-detonation radionuclide decay-chain network coupled to detailed models simulating physical production and transport processes affecting the gas signature. Because these models are governed by uncertain parameters including barometric fluctuations, realistic ranges of variation for each of the parameters influencing isotopic composition are then defined. A Latin-Hypercube sampling approach is used to obtain a random distribution of isotopic production and gas transport results associated with a given value of each parameter. We apply these results to background histories of two stations, one providing 4-isotope background measurements and the other providing two-isotope measurements associated with the 2013 DPRK announced UNE.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2025

Multi-modal effects on indirect noise induced by turbulent entropy fields

Planar entropy waves are commonly assumed for predicting indirect combustion noise. However, the non-planar and turbulent nature of flows found in most practical combustors challenges this assumption. In the present paper, we examine the indirect noise generated by non-planar and turbulent entropy fields through subsonic nozzles. Firstly, we introduce a new transfer function framework that accounts for the contribution of non-planar Fourier modes of the entropy field to the indirect noise spectra. When applied to a turbulent flow field, this method demonstrates a significant improvement in spectral predictions compared with a conventional approach that only considers the planar mode. Secondly, simulations show that non-planar Fourier modes become significant above a threshold frequency ƒ thresh , found in the mid- to high-frequency range. This contribution of non-planar modes is explained by two-dimensional shear effects that distort the entropy waves. A scaling relation that uses residence times along streamlines is developed for ƒ thresh , showing good agreement with simulation results. Finally, we show that the indirect noise from non-planar entropy modes found in aviation combustors can be significant at frequencies below 1 kHz, which might be relevant in situations of thermo-acoustic instabilities coupled to indirect noise.

42 ENGINEERING↗
DOE OSTI2025

Particle‐In‐Cell Simulations of Starfish Prime

The Starfish Prime high altitude nuclear test created a transient diamagnetic cavity in the Earth's magnetic field above Johnston Island and launched an electromagnetic pulse (EMP) that was detected around the globe. For this work, we use the ion‐kinetic particle‐in‐cell code Topanga to simulate diamagnetic cavity evolution and the E3 EMP signal for Starfish Prime out to over a minute of physical time. The simulation domain has a longitudinal and latitudinal extent of 60° X 60° and a vertical extent of 2,000 km from the surface of the Earth. We compare our simulated results to magnetometer measurements taken in space and on the ground, finding good agreement in both cases. The diamagnetic cavity in the simulation forms in about a second, while the associated debris flux tube takes approximately 30 s to decay. The debris flux tube undergoes significant motion during this timeframe, rising upward. The measured E3 EMP signal on the ground consists of several components, all of which are present in our simulations. We discuss the physical origin of these components in relation to E3a (blast) and E3b (heave).

Belyaev, Mikhail A. [Lawrence Livermore National L↗
DOE OSTI2023

Mass and momentum transport in the Tilted Rocket Rig experiment

Here, high-fidelity large eddy simulations of the inclined Rayleigh-Taylor instability in the Tilted Rocket Rig experimental configuration [Smeeton and Youngs, AWE Report No. O 35/87 (1987)] are performed using a tenth-order compact finite difference code. These simulations are analyzed for spatial distributions of turbulent kinetic energy, turbulent mass flux velocity, species mass fraction flux, species mass fraction variance, and Favre-averaged Reynolds stresses at two time instances, t = 37 ms and t = 52 ms. Additionally, the vertical distribution of the components of the unclosed budget equations over the center of the domain for these quantities are examined. The dominant terms of these budget equations are further decomposed to examine the principal contributions to these terms in each axis. Notably, the principal contribution to the horizontal turbulent mass flux velocity budget is found to be from a term which is commonly neglected in many Reynolds-averaged Navier-Stokes models.

42 ENGINEERING↗