Search NASA⌕ Search

SEARCH · Search NASA

Results for “scalar dissipation rate”

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 37 records · Page 2

Self-similar Reynolds-averaged mechanical–scalar turbulence models for reshocked Richtmyer–Meshkov instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Richtmyer–Meshkov instability-induced mixing in planar geometry derived in the small Atwood number limit are extended to construct models for reshocked Richtmyer–Meshkov mixing. In this study, the models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ modeled either differentially or algebraically. The three- and four-equation models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing. Mixing layer growth parameters and other physical observables were obtained explicitly as functions of the model coefficients and were used to calibrate the model coefficients. Here, the solutions for the singly shocked Richtmyer–Meshkov case for the mixing layer width and the turbulent fields are used to construct piecewise-continuous generalizations of these quantities for times after reshock. For generality, the post-reshock mixing layer width is not assumed to grow with the same power-law as the pre-reshock width, and an impulsive approximation applied to Rayleigh–Taylor instability growth is used to establish the expression for the post-reshock width. A four-equation model is then used to illustrate the spatiotemporal behavior of the mean and turbulent fields and late-time turbulent equation budgets across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for reshocked Richtmyer–Meshkov instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

Scalar Dissipation Modeling for Passive and Active Scalars: a priori Study Using Direct Numerical Simulation

Transitional databases from Direct Numerical Simulation (DNS) of three-dimensional mixing layers for single-phase flows and two-phase flows with evaporation are analyzed and used to examine the typical hypothesis that the scalar dissipation Probability Distribution Function (PDF) may be modeled as a Gaussian. The databases encompass a single-component fuel and four multicomponent fuels, two initial Reynolds numbers (Re), two mass loadings for two-phase flows and two free-stream gas temperatures. Using the DNS calculated moments of the scalar-dissipation PDF, it is shown, consistent with existing experimental information on single-phase flows, that the Gaussian is a modest approximation of the DNS-extracted PDF, particularly poor in the range of the high scalar-dissipation values, which are significant for turbulent reaction rate modeling in non-premixed flows using flamelet models. With the same DNS calculated moments of the scalar-dissipation PDF and making a change of variables, a model of this PDF is proposed in the form of the (beta)-PDF which is shown to approximate much better the DNS-extracted PDF, particularly in the regime of the high scalar-dissipation values. Several types of statistical measures are calculated over the ensemble of the fourteen databases. For each statistical measure, the proposed (beta)-PDF model is shown to be much superior to the Gaussian in approximating the DNS-extracted PDF. Additionally, the agreement between the DNS-extracted PDF and the (beta)-PDF even improves when the comparison is performed for higher initial Re layers, whereas the comparison with the Gaussian is independent of the initial Re values. For two-phase flows, the comparison between the DNS-extracted PDF and the (beta)-PDF also improves with increasing free-stream gas temperature and mass loading. The higher fidelity approximation of the DNS-extracted PDF by the (beta)-PDF with increasing Re, gas temperature and mass loading bodes well for turbulent reaction rate modeling.

scalar dissipation modeling↗

Conditional analysis of temperature and strain rate effects on dissipation structure in turbulent non-premixed jet flames

Here, this work presents results from simultaneous high-resolution temperature and velocity measurements in a series of turbulent non-premixed jet flames. The filtered Rayleigh scattering (FRS)-based temperature measurements demonstrate sufficient signal-to-noise (SNR) and spatial resolution to estimate the smallest scalar length scales and accurately determine dissipation rate fields. A comprehensive set of conditional statistics are used to characterize the small-scale structure, including the dependence of dissipation layer widths on Reynolds number, temperature, and dissipation magnitude. In general, the dissipation layer thickness decrease with increasing Reynolds number and increase with increasing temperature. However, dissipation layer widths show two distinct behaviors with respect to dissipation magnitude. For small dissipation values, increases in magnitude results in broadening of the dissipation layer, while for larger magnitude values of dissipation, the layer widths are thinned, highlighting the complexity of small-scale turbulent mixing. Additionally, measured ratios of the dissipation layer width to the Batchelor length scale are consistent across all Reynolds numbers and agree with previous studies in non-reacting flows. The unique aspect about the current set of measurements is the ability to examine the interaction of dissipation structure with turbulent flow parameters for the first time in turbulent non-premixed flames. Particularly, the strain rate/dissipation relationship is examined and compared to previous studies in non-reacting flows. It is found that the dissipation layers tend to align normal to the principal compressive strain axis and this tendency increases with increasing Reynolds number. For the lowest Reynolds number case, no dependence of the dissipation layer width nor dissipation rate magnitude on strain rate is found. However, for higher Reynolds numbers, a strong dependence of the dissipation layer width and dissipation rate magnitude on the principal compressive strain rate is observed. These results indicate the direct role of the compressive strain rate field on small-scale mixing structure in reacting flows.

42 ENGINEERING↗

${\tt Warm}$SPy: a numerical study of cosmological perturbations in warm inflation

We present ${\tt Warm}$SPy, a numerical code in Python designed to solve for the perturbations' equations in warm inflation models and compute the corresponding scalar power spectrum at CMB horizon crossing. In models of warm inflation, a radiation bath of temperature T during inflation induces a dissipation (friction) rate of strength Q ∝ T c /Φ m in the equation of motion for the inflaton field Φ. While for a temperature-independent dissipation rate (c = 0) an analytic expression for the scalar power spectrum exists, in the case of a non-zero value for c the set of equations can only be solved numerically. For c > 0 (c < 0), the coupling between the perturbations in the inflaton field and radiation induces a growing (decaying) mode in the scalar perturbations, generally parameterized by a multiplicative function G(Q) which we refer to as the scalar dissipation function. Using ${\tt Warm}$SPy, we provide an analytic fit for G(Q) for the cases of c = {3,1,-1}, corresponding to three cases that have been realized in physical models. Compared to previous literature results, our fits are more robust and valid over a broader range of dissipation strengths Q ϵ [10 -7 ,10 4 ]. Additionally, for the first time, we numerically assess the stability of the scalar dissipation function against various model parameters, inflationary histories as well as the effects of metric perturbations. As a whole, the results do not depend appreciably on most of the parameters in the analysis, except for the dissipation index c, providing evidence for the universal behaviour of the scalar dissipation function G(Q).

79 ASTRONOMY AND ASTROPHYSICS↗

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗

A critical comparison of second order closures with direct numerical simulation of homogeneous turbulence

Recently, several second order closure models have been proposed for closing the second moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behavior in different flows have not been fully examined because they are un-measurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the opportunity to do this kind of study. With the direct numerical simulation, we may use the solution to exactly calculate the values of these correlation terms and then directly compare them with the values from their modeled formulations (models). Here, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of forty five homogeneous flows which were done by Rogers et al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones which give the best performance. The modeling procedure, model constraints, and the various evaluated models are described. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

Direct numerical simulation of turbulent H2-O2 combustion using reduced chemistry

Results of direct numerical simulations of hydrogen-oxygen combustion using a partial-equilibrium chemistry scheme in constant density, decaying, isotropic turbulence are reported. The simulations qualitatively reproduce many features of experimental results, such as superequilibrium radical species mole fractions, with temperature and major species mole fractions closer to chemical equilibrium. It was also observed that the peak reaction rates occur in narrow zones where the stoichiometric surface intersects regions of high scalar dissipation, as might be expected for combustion conditions close to chemical equilibrium. Another finding was that high OH mole fraction correspond more closely to the stoichiometric surface than to areas of high reaction rate for conditions of the simulations. Simulation results were compared to predictions of the Conditional Moment Closure model. This model was found to give good results for all quantities of interest when the conditionally averaged scalar dissipation was used in the prediction. When the nonconditioned average dissipation was used, the predictions compared well to the simulations for most of the species and temperature, but not for the reaction rate. The comparison would be expected to improve for higher Reynolds number flows, however.

Montgomery, Christopher J.↗

'Reduced' magnetohydrodynamics and minimum dissipation rates

It is demonstrated that all solutions of the equations of 'reduced' magnetohydrodynamics approach a uniform-current, zero-flow state for long times, given a constant wall electric field, uniform scalar viscosity and resistivity, and uniform mass density. This state is the state of minimum energy dissipation rate for these boundary conditions. No steady-state turbulence is possible. The result contrasts sharply with results for full three-dimensional magnetohydrodynamics before the reduction occurs.

Montgomery, David↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗

A transport equation for the scalar dissipation in reacting flows with variable density: First results

Although the different regimes of premixed combustion are not well defined, most of the recent developments in turbulent combustion modeling are led in the so-called flamelet regime. The goal of these models is to give a realistic expression to the mean reaction rate (w). Several methods can be used to estimate (w). Bray and coworkers (Libby & Bray 1980, Bray 1985, Bray & Libby 1986) express the instantaneous reaction rate by means of a flamelet library and a frequency which describes the local interaction between the laminar flamelets and the turbulent flowfield. In another way, the mean reaction rate can be directly connected to the flame surface density (Sigma). This quantity can be given by the transport equation of the coherent flame model initially proposed by Marble & Broadwell 1977 and developed elsewhere. The mean reaction rate, (w), can also be estimated thanks to the evolution of an arbitrary scalar field G(x, t) = G(sub O) which represents the flame sheet. G(x, t) is obtained from the G-equation proposed by Williams 1985, Kerstein et al. 1988 and Peters 1993. Another possibility proposed in a recent study by Mantel & Borghi 1991, where a transport equation for the mean dissipation rate (epsilon(sub c)) of the progress variable c is used to determine (w). In their model, Mantel & Borghi 1991 considered a medium with constant density and constant diffusivity in the determination of the transport equation for (epsilon(sub c)). A comparison of different flamelet models made by Duclos et al. 1993 shows the realistic behavior of this model even in the case of constant density. Our objective in this present report is to present preliminary results on the study of this equation in the case of variable density and variable diffusivity. Assumptions of constant pressure and a Lewis number equal to unity allow us to significantly simplify the equation. A systematic order of magnitude analysis based on adequate scale relations is performed on each term of the equation. As in the case of constant density and constant diffusivity, the effects of stretching of the scalar field by the turbulent strain field, of local curvature, and of chemical reactions are predominant. In this preliminary work, we suggest closure models for certain terms, which will be validated after comparisons with DNS data.

Mantel, T.↗

Enhanced damping of Alfven waves in the solar corona by a turbulent wave spectrum

The effect of a background spectrum of Alfven waves on the rate of dissipation of a test shear Alfven wave is numerically calculated. The results demonstrate that as the classical resistivity eta and classical viscosity mu become small, the damping rate of the Alfven wave remains large and depends only on the amplitude for the scalar potential of the wave spectrum and the wavenumber of the Alfven wave. The damping rate is virtually independent of eta and mu. The wave spectrum need not be turbulent or stochastic to affect the damping rate. The dissipation rate is nonlinear enhanced by nonstochastic spectra as well as by stochastic spectra if two conditions are met. First, the perpendicular magnetic field associated with Alfven wave spectrum must exceed a certain collision-frequency threshold and second, for nonstochastic spectra only, the magnetic field must exceed a threshold proportional to the parallel wavenumber of the shear Alfven wave. These conditions can be easily satisfied in the solar corona.

Kleva, Robert G.↗

Investigation of Strain/Vorticity and Large-Scale Flow Structure in Turbulent Nonpremixed Jet Flames

Our study will use the microgravity environment to investigate the underlying flow structure of turbulent nonpremixed round jet flames. In particular, we aim to investigate the large-scale turbulent structure using planar laser Mie scattering (PLMS), and the strain rate and vorticity fields using particle image velocimetry (PIV). This work is motivated by recent studies in our laboratory that have led to several interesting observations of nominally momentum-driven turbulent nonpremixed planar flames. First of all, the organized large-scale turbulent structures that are observed in nonreacting planar jets may be substantially modified or suppressed in nonpremixed planar jet flames. Furthermore, a recent study using PIV and planar laser-induced fluorescence of OH has shown that in transitional and turbulent nonpremixed planar jet flames the presence of the flame seems to greatly influence the underlying vorticity and strain fields, as compared to nonreacting jets. For example, the reaction zones in the jet flames are strongly correlated with regions of high vorticity. A related study has demonstrated that vorticity is not correlated in the same way with either iso-scalar surfaces or scalar dissipation layers in nonreacting planar jets. Furthermore, the relationship between strain and the reaction zone appears to be modified by the presence of high levels of heat release. In particular, the strain rate field in planar jet flames exhibits a preferred direction of principal compressive strain that apparently is related to strong shear across the reaction zone. This preferred direction of strain was not observed in nonreacting jets. One of the major problems encountered when conducting these types of studies is that it is difficult to know to what extent buoyancy influences the results. Therefore, the microgravity environment provides us with an excellent opportunity to explore these issues without the complicating effects of buoyancy. This is particularly the case when studying flames that are transitional between laminar and turbulent states. For example, the strong correlation of vorticity with the reaction zone (discussed above) was observed in both transitional and turbulent planar flames, but the effect was stronger for the transitional case. To date, the reason for the presence of the vorticity-reaction zone correlation is not known, although vorticity production via baroclinic torque is a likely cause. The microgravity environment will allow us to specifically determine whether the vorticity is produced by baroclinic torque resulting from the flame density gradient acting with the hydrostatic pressure gradient.AAdditional details of the planned experiments are provided.

Clemens, N. T.↗

Adiabatic Evolution Due to the Conservative Scalar Self-force During Orbital Resonances

We calculate the scalar self-force experienced by a scalar point-charge orbiting a Kerr black holealongrθ-resonant geodesics. We use the self-force to calculate the averaged rate of change of thecharge’s orbital energy⟨ ̇E⟩, angular momentum⟨ ̇Lz⟩, and Carter constant⟨ ̇Q⟩, which togethercapture the leading-order adiabatic, secular evolution of the point-charge. Away from resonances,only the dissipative (time antisymmetric) components of the self-force contribute to⟨ ̇E⟩,⟨ ̇Lz⟩, and⟨ ̇Q⟩. We demonstrate, using a new numerical code, that duringrθresonances conservative (timesymmetric) scalar perturbations also contribute to⟨ ̇Q⟩and, thus, help drive the adiabatic evolutionof the orbit. Furthermore, we observe that the relative impact of these conservative contributionsto⟨ ̇Q⟩is particularly strong for eccentric 2:3 resonances. These results provide the first conclusivenumerical evidence that conservative scalar perturbations of Kerr spacetime are nonintegrable duringrθresonances.

Zachary Nasipak↗

Constraints on the scalar-field potential in warm inflation

In this work, we quantify the degree of fine-tuning required to achieve an observationally viable period of inflation in the strongly dissipative regime of warm inflation. The “fine-tuning” parameter λ is taken to be the ratio of the change in the height of the potential Δ V to the change in the scalar field (ΔΦ) 4 , i.e., the width of the potential, and therefore measures the requisite degree of flatness in the potential. The best motivated warm inflationary scenarios involve a dissipation rate of the kind Γ ∝ T c with c ≥ 0 , and for all such cases, the bounds on λ are tighter than those for standard cold inflation by at least 3 orders of magnitude. In other words, these models require an even flatter potential than standard inflation. On the other hand for the case of warm inflation with c < 0, we find that in a strongly dissipative regime the bound on λ can significantly weaken with respect to cold inflation. Thus, if a warm inflation model can be constructed in a strongly dissipative, negatively temperature-dependent regime, it accommodates steeper potentials otherwise ruled out in standard inflation.

79 ASTRONOMY AND ASTROPHYSICS↗

Second order modeling of boundary-free turbulent shear flows

A set of realizable second order models for boundary-free turbulent flows is presented. The constraints on second order models based on the realizability principle are re-examined. The rapid terms in the pressure correlations for both the Reynolds stress and the passive scalar flux equations are constructed to exactly satisfy the joint realizability. All other model terms (return-to-isotropy, third moments, and terms in the dissipation equations) already satisfy realizability. To correct the spreading rate of the axisymmetric jet, an extra term is added to the dissipation equation which accounts for the effect of mean vortex stretching on dissipation. The test flows used in this study are the mixing shear layer, plane jet, axisymmetric jet, and plane wake. The numerical solutions show that the unified model equations predict all these flows reasonably. It is expected that these models would be suitable for more complex and critical flows.

Shih, T.-H.↗

Second order modeling of boundary-free turbulent shear flows

A set of realizable second order models for boundary-free turbulent flows is presented. The constraints on second order models based on the realizability principle are re-examined. The rapid terms in the pressure correlations for both the Reynolds stress and the passive scalar flux equations are constructed to exactly satisfy the joint realizability. All other model terms (return-to-isotropy, third moments, and terms in the dissipation equations) already satisfy realizability. To correct the spreading rate of the axisymmetric jet, an extra term is added to the dissipation equation which accounts for the effect of mean vortex stretching on dissipation. The test flows used in this study are the mixing shear layer, plane jet, axisymmetric jet, and plane wake. The numerical solutions show that the unified model equations predict all these flows reasonably. It is expected that these models would be suitable for more complex and critical flows.

Shih, T.-H.↗