Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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

Soil Moisture Buffers the Impact of Precipitation Variability on Ecosystem Productivity

Water availability governs ecosystem productivity, yet estimates of vegetation sensitivity to water can differ greatly depending on whether the sensitivity is examined spatially or temporally. In particular, the spatial sensitivity is often reported to be much stronger than temporal sensitivities, leading to highly uncertain projections of ecosystem responses to future climate change when using space-for-time substitution. The large difference between spatial and temporal sensitivities remains unexplained. Prior research, however, primarily relied on precipitation as the water availability proxy, whereas vegetation responds to soil moisture. Here, we combined satellite estimates of vegetation productivity with soil moisture data across water-limited ecosystems of the continental United States (CONUS) to identify a convergent sensitivity of productivity to water availability. Using precipitation, we show that temporal sensitivity is 66% lower than spatial sensitivity overall. Our analysis identified the cause of the difference to be primarily driven by the seasonal variability of water availability, rooting depth, and soil properties. When using soil moisture instead of precipitation, we observed widespread convergence in the spatial and temporal sensitivities—that is, the two sensitivities became much more similar in magnitude across all water-limited ecosystems within CONUS. These results show that overlooking soil hydrology can inflate perceived discrepancies between spatial and temporal vegetation sensitivities, leading to biased projections of ecosystem dynamics under future hydro-climatic change.

Wang, Huiqi [University of California, Berkeley, C↗
NASA NTRS1980

Constrained optimum trajectories with specified range

The characteristics of optimum fixed-range trajectories whose structure is constrained to climb, steady cruise, and descent segments are derived by application of optimal control theory. The performance function consists of the sum of fuel and time costs, referred to as direct operating costs (DOC). The state variable is range-to-go and the independent variable is energy. In this formulation a cruise segment always occurs at the optimum cruise energy for sufficiently large range. At short ranges (500 n. mi. and less) a cruise segment may also occur below the optimum cruise energy. The existence of such a cruise segment depends primarily on the fuel flow vs thrust characteristics and on thrust constraints. If thrust is a free control variable along with airspeed, it is shown that such cruise segments will not generally occur. If thrust is constrained to some maximum value in climb and to some minimum in descent, such cruise segments generally will occur. The performance difference between free thrust and constrained thrust trajectories has been determined in computer calculations for an example transport aircraft.

Erzberger, H.↗
NASA NTRS1984

Horizontal energy propagation in a baroclinic atmosphere with meridional and zonal structure

A baroclinic model was developed. The model is linear, steady state, baroclinic (2 level) and based on the primitive equations. The unique features of the model are: (1) The basic state u and T fields are functions of wavelength, theta, and rho, (2) A mean meridional circulation v, w, which may be a function of wavelength, theta, and rho, and p may be included, (3) Dissipation may be a function of theta, wavelength, rho, p (e.g., differ over land and sea). Because the basic state is a function of longitude, the equations are not separable in wavelength and the usual semi-spectral solution procedure is not applicable. The model solutions can be found in a straight forward (brute force) way by discretizing the model equations directly on a wavelength, theta, rho, p grid. There are six dependent (perturbation) variables: u, v, and geopotential at levels 1 and 2. The thermal forcing is prescribed at the intermediate vertical level. The linearized equations may be written.

Youngblut, C.↗
NASA NTRS2020

Discovery of Thermonuclear (Type 1) X-ray Bursts in the X-ray Binary Swift J1858.6-0814 Observed with NICER and NuSTAR

Swift J1858.6–0814 is a recently discovered X-ray binary notable for extremely strong variability (by factors of >100 in soft X-rays) in its discovery state. We present the detection of five thermonuclear (Type I) X-ray bursts from Swift J1858.6–0814, implying that the compact object in the system is a neutron star (NS). Some of the bursts show photospheric radius expansion, so their peak flux can be used to estimate the distance to the system. The peak luminosity, and hence distance, can depend on several system parameters; for the most likely values, a high inclination and a helium atmosphere, D = 12.8+0.8−0.6 kpc, although systematic effects allow a conservative range of 9–18 kpc. Before one burst, we detect a QPO at 9.6 ± 0.5 mHz with a fractional rms amplitude of 2.2 ± 0.2 per cent (0.5–10 keV), likely due to marginally stable burning of helium; similar oscillations may be present before the other bursts but the light curves are not long enough to allow their detection. We also search for burst oscillations but do not detect any, with an upper limit in the best case of 15 per cent fractional amplitude (over 1–8 keV). Finally, we discuss the implications of the NS accretor and this distance on other inferences which have been made about the system. In particular, we find that Swift J1858.6–0814 was observed at super-Eddington luminosities at least during bright flares during the variable stage of its outburst.

D J K Buisson↗
DOE OSTI2025

Impact of grassland fire dynamics on particulate emission factors

Understanding particulate emission factors in grassland fires is critical for improving air quality assessments and refining emissions inventories. This study explores the dynamic and heterogeneous nature of fire processes and their influence on particulate emissions. Traditional emission factor inventories often rely on static values, overlooking the spatially and temporally variable dynamics of fire behavior. Through field experiments conducted at Konza Prairie Biological Station measuring emissions with uncrewed aerial vehicles, and computational modeling using HIGRAD/FIRETEC, this research highlights the dependency of emission production on localized fire dynamics, including ignition patterns, fuel properties, and atmospheric conditions. The results demonstrate a need for fire dynamic considerations to accurately capture emissions variability.

09 BIOMASS FUELS↗
DOE OSTI2020

Spectral variability of a sample of extreme variability quasars and implications for the MgII broad-line region

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed bymore than 1.5mag in the g band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey epochs (separated by a few years in the quasar rest frame). These EVQs are selected from quasars in the SDSS Stripe 82 region, covering a redshift range of 0.5 < z < 2.1. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the g band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on time-scales of only months. To increase sample statistics, we use a supplemental sample of 33 EVQs with multi-epoch spectra from SDSS that cover the broad Mg II lambda 2798 line. Leveraging on the large dynamic range in continuum variability between the multi-epoch spectra, we explore the associated variations in the broad MgII line, whose variability properties have not been well studied before. The broad MgII flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the full width at half-maximum (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad MgII to estimate the black hole mass with single epoch spectra therefore introduces a luminosity-dependent bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Clathrate BaNi 2 P 4 : An Interplay of Heat and Charge Transport Due to Strong Host–Guest Interactions

Heat and charge transport properties of the metallic unconventional clathrate BaNi 2 P 4 , hosting Ba cations in oversized Ni 8 P 16 cages, are investigated. A novel method of single-crystal growth was developed, yielding 2–3 mm sized crystals of BaNi 2 P 4 . We also developed a setup to accurately measure thermal conductivity and electrical resistivity of the synthesized single crystals in a wide temperature range avoiding crystal remounting. BaNi 2 P 4 has a metallic temperature dependence of its electrical resistivity (decreasing with decreasing temperature) and manifests an unconventional T 2 power law for 50 K < T < 300 K; below 50 K, the power-law exponent increases gradually such that below 10 K the power law is T 5 , a predicted but extremely rarely experimentally observed dependence for peculiar electron–phonon interactions. Furthermore, electronic band structure calculations, consistent with measurements of de Haas–van Alphen oscillations, show large band dispersions with significant contributions of Ba orbitals to states near the Fermi level, which is atypical for clathrates. The thermal properties of BaNi 2 P 4 were probed using a combination of variable-temperature single-crystal X-ray diffraction experiments, heat capacity measurements, first-principles phonon dispersion calculations, and inelastic neutron scattering measurements. BaNi 2 P 4 exhibits significant hybridization of the Ba-guest and Ni–P-framework vibrational modes, which may be enhanced via the detected split of the Ba position, which results in strong Ba-framework interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Defense Waste Processing Facility Nitric-Glycolic Flowsheet Chemical Process Cell Chemistry: Part 2

The conversions of nitrite to nitrate, the destruction of glycolate, and the conversion of glycolate to formate and oxalate were modeled for the Nitric-Glycolic flowsheet using data from Chemical Process Cell (CPC) simulant runs conducted by Savannah River National Laboratory (SRNL) from 2011 to 2016. The goal of this work was to develop empirical correlation models to predict these values from measurable variables from the chemical process so that these quantities could be predicted a-priori from the sludge or simulant composition and measurable processing variables. The need for these predictions arises from the need to predict the REDuction/OXidation (REDOX) state of the glass from the Defense Waste Processing Facility (DWPF) melter. This report summarizes the work on these correlations based on the aforementioned data. Previous work on these correlations was documented in a technical report covering data from 2011-2015. This current report supersedes this previous report. Further refinement of the models as additional data are collected is recommended. The glass REDOX depends on the concentrations of nitrate and manganese (oxidants), and of glycolate, formate, oxalate, carbon, and antifoam (reductants) in the melter feed. The waste sludge contains nitrite, nitrate, manganese (Mn), and oxalate. Virtually all of the nitrite is converted to nitrate or NO+NO 2 +N 2 O gases in the CPC. The portion of the nitrite converted to nitrate increases the amount of nitrate in the sludge. The amount of glycolate in the final melter feed depends on the amount of the glycolic acid feed that is destroyed. Similarly, the amounts of formate and oxalate formed during the decomposition of glycolic acid are required. The material balance on carbon was found to not close in most cases. Generally, there was less carbon at the end of testing compared to the inputs. The most uncertain product variable was glycolate, so material balances were performed where the glycolate concentration was adjusted, usually upward, to close the balance. Correlation versus the original, as-measured, data was generally poor, but correlation against the material balance adjusted values was greatly improved. It was also shown that the correlation of the measured REDOX versus the predicted REDOX was much better when the material balance adjusted glycolate values were used. Three data series were primarily used during the regressions of the data; these series were 1) Sludge Batch 9 NG flowsheet simulant runs NG51-62 (SB9-NG); 2) Scaled Runs + Bounding Hydrogen Runs (SR+BH); and 3) Runs GN43-50 and 57 (43-50,57). The glycolate destruction was found to correlate with acid stoichiometry (AS), percent reducing acid (PRA), and for some data series, headspace to simulant volume ratio (HSV), mercury (Hg), and nitrate. Although glycolate destruction for pairs of data series (e.g., [SB9-NG] and [SR+BH]) were found to depend on HSV, the combination of all three data series was not found to have significant dependence on this variable. The best model for glycolate destruction depended on AS, nitrate, and Hg. This model predicted the product glycolate compositions of the data to within 92-106%. The conversion of glycolate to formate was high when noble metals and Hg were not present, with values up to 100%. When noble metals and Hg were present, this conversion ranged from zero to 7%, and was dependent on AS. Lower AS gave higher conversions to formate. The conversion to oxalate was found to depend on the AS and the initial concentration of nitrite. An alternative fit versus AS and the form of ruthenium (Ru) used is a possible alternative. This fit was somewhat less statistically significant. This second model predicts that more oxalate is formed when Ru-nitrosyl nitrate is used rather than Ru chloride. The conversion of glycolate to oxalate ranged from zero to 6%. The conversion of nitrite to nitrate depended primarily on AS and PRA, with HSV and Hg being significant when these variables were varied. For multiple series of data, nitrite was also needed to SRNL-STI-2017-00172 5HYLVLRQ viL distinguish between data series, and the effect of HSV became insignificant. The best model for nitrite to nitrate conversion depended on AS, PRA, nitrite, and Hg. The 95% confidence intervals on the predicted values of glycolate destruction, glycolate to oxalate conversion, and nitrite to nitrate conversion were used to determine the uncertainty in the predicted REDOX when starting with only the composition of the sludge, AS, and PRA. Using the 95% confidences on an individual value (that is the confidence in getting a particular value for one single test as opposed to what the mean would be for multiple tests), the uncertainty in the predicted REDOX was calculated. The uncertainty in the actual product composition glycolate, oxalate, formate, and nitrate concentrations translated to an uncertainty in the REDOX value of ±0.1,which is approximately the uncertainty claimed in the REDOX model itself.

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

Iodate interactions with calcite: implications for natural attenuation

Solid-phase interactions and speciation are important to radioiodine transport in groundwater. At the Hanford Site in Southeastern Washington State, iodate (IO 3 - ) is the main aqueous species in dilute radioiodine groundwater plumes. Like other oxyanions, IO 3 - may be incorporated into and/or adsorbed onto calcite, a common mineral at Hanford, decreasing its mobility in the environment. A series of macroscale batch experiments combined with solid phase characterization were conducted to identify variables impacting time-dependent aqueous IO 3 - removal via calcite precipitation and determine the location of IO 3 - within the calcite crystal structure. Results demonstrated 11.5-97% aqueous IO 3 - removal during initial rapid calcite precipitation. Incorporation was apparently the main removal mechanism, although later slower precipitation and/or adsorption may have also contributed to IO 3 - removal. Using a higher concentration of the calcite forming solutions (i.e., using 1M vs. 0.1M concentrations) resulted in an increase in the amount of precipitated calcite and a greater percentage of IO 3 - removed; however, calcite formed with lower molarity solutions resulted in higher IO 3 - mass (µg/g) removal. Solubility testing of laboratory produced calcites showed only small differences in solubility for calcite with and without IO 3 - incorporated into the structure. Evidence collected from SEM/FIB and TEM/SAED suggested that the IO 3 - incorporated into calcite was present in regions close to surface (implying potential easy release upon calcite dissolution).

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

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗
DOE OSTI2021

Estimating Submicron Aerosol Mixing State at the Global Scale With Machine Learning and Earth System Modeling

Abstract This study integrates machine learning and particle‐resolved aerosol simulations to develop emulators that predict submicron aerosol mixing state indices from the Earth system model (ESM) simulations. The emulators predict aerosol mixing state using only quantities that are predicted by the ESM, including bulk aerosol species concentrations, which do not by themselves carry mixing state information. We used PartMC‐MOSAIC as the particle‐resolved model and NCAR's CESM as the ESM. We trained emulators for three different mixing state indices for submicron aerosol in terms of chemical species abundance ( χ a ), the mixing of optically absorbing and nonabsorbing species ( χ o ), and the mixing of hygroscopic and nonhygroscopic species ( χ h ). Our global mixing state maps show considerable spatial and seasonal variability unique to each mixing state index. Seasonal averages varied spatially between 13% and 94% for χ a , between 38% and 94% for χ o , and between 20% and 87% for χ h with global annual averages of 67%, 68%, and 56%, respectively. High values in one index can be consistent with low values in another index depending on the grouping of species and their relative abundance, meaning that each mixing state index captures different aspects of the population mixing state. Although a direct validation with observational data has not been possible yet, our results are consistent with mixing state index values derived from ambient observations. This work is a prototypical example of using machine learning emulators to add information to ESM simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Geometric adiabatic angle in anisotropic oscillators

We discuss a classical anisotropic oscillator and the Foucault pendulum as examples illustrating non-conservation of action variables in integrable classical mechanical systems with adiabatically slow evolution. We also emphasize the importance of the mass parameter of a harmonic oscillator, alongside its frequency, in explicitly time-dependent situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Sigma: Uncertainty Propagation for C++

Sigma is a header-only C++-17 library for uncertainty propagation, inspired by uncertainties for Python and Measurements.jl for Julia. The library tracks the functional correlation between dependent and independent variables, ensuring that the uncertainty of the independent variables is properly considered in the calculation of the dependent variables’ uncertainties. It is intended as a near drop-in replacement for the standard floating point types (aside from uncertainty specification), and aims to be easily interoperable with the existing standard types.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Time-controlled descent guidance in uncertain winds

A procedure has been developed for constructing a statistical model of the altitude-dependent mean wind profile from the historical record of wind measurements at particular locations. The model is constructed by fitting a Markov process, with altitude as the stage variable, to the historical wind data. The wind model, together with the aircraft dynamics and the error characteristics of the navigation system, are incorporated in the design of a state estimator, which gives the minimum variance estimate of the aircraft state and the wind vector. The state and wind estimates are used as inputs to a linear feedback law for guiding the aircraft along the nominal trajectory. An example design of a time-constrained (4D RNAV) descent guidance system is presented, showing tracking accuracy, control activity, and probability of arrival time with and without the wind estimator.

Menga, G.↗
DOE OSTI2025

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗
DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗
DOE OSTI2021

A Differential Measurement of Circumstellar Extinction for AA Tau’s 2011 Dimming Event

AA Tau is a classical T Tauri star with a highly inclined, warped circumstellar disk. For decades, AA Tau exhibited photometric and spectroscopic variabilities that were successfully modeled as occultations of the primary star by circumstellar material. In 2011, AA Tau entered an extended faint state, presumably due to enhanced levels of circumstellar dust. We use two sets of contemporaneous optical–near-IR (NIR) spectra of AA Tau, obtained in 2008 and 2014 December, to directly measure the wavelength-dependent extinction associated with the dust enhancement driving AA Tau’s 2011 optical fade. Taken alone, AA Tau’s apparent optical–NIR increased extinction curve cannot be fit well with standard extinction laws. At optical wavelengths, AA Tau’s dimming event is consistent with predictions of common extinction models for an increase of A {sub V} = 2, but no such model reproduces AA Tau’s color–color excess at NIR wavelengths. We show that veiling emission accounts for the apparent anomalous NIR extinction: after including this veiling flux, AA Tau’s dimming event is consistent with a standard A {sub V} = 2 extinction law across the full optical–NIR range. We also report an increase in AA Tau’s mid-IR flux since its 2011 fade and suggest that an increase in the height of AA Tau’s inner disk is the most likely explanation for both the IR brightening and the additional extinction along the line of sight. In addition to informing our understanding of AA Tau, this analysis demonstrates that caution should be exercised when inferring extinction (and stellar parameters) from the NIR color–color excess of young stars with measurable NIR veiling fluxes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗