Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct numerical simulation (DNS)”

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 91 records · Page 5

Reynolds-number scaling of wall-pressure–velocity correlations in wall-bounded turbulence

Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure–velocity coupling requires investigation to properly design a controller with such input data so that it can actuate upon the desired turbulent structures. A comprehensive database from direct numerical simulations (DNS) of turbulent channel flow is used for this purpose, spanning a Reynolds-number range$Re_\tau \approx 550\unicode{x2013}5200$. Spectral analysis reveals that the streamwise velocity is most strongly coupled to the linear term of the wall pressure, at a Reynolds-number invariant distance-from-the-wall scaling of$\lambda _x/y \approx 14$(and$\lambda _x/y \approx 8$for the wall-normal velocity). When extending the analysis to both homogeneous directions in$x$and$y$, the peak coherence is centred at$\lambda _x/\lambda _z \approx 2$and$\lambda _x/\lambda _z \approx 1$for$p_w$and$u$, and$p_w$and$v$, respectively. A stronger coherence is retrieved when the quadratic term of the wall pressure is concerned, but there is only little evidence for a wall-attached-eddy type of scaling. An experimental dataset comprising simultaneous measurements of wall pressure and velocity complements the DNS-based findings at one value of$Re_\tau \approx 2$k, with ample evidence that the DNS-inferred correlations can be replicated with experimental pressure data subject to significant levels of (acoustic) facility noise. It is furthermore shown that velocity-state estimations can be achieved with good accuracy by including both the linear and quadratic terms of the wall pressure. An accuracy of up to 72 % in the binary state of the streamwise velocity fluctuations in the logarithmic region is achieved; this corresponds to a correlation coefficient of$\approx$0.6. This thus demonstrates that wall-pressure sensing for velocity-state estimation – e.g. for use in real-time control of wall-bounded turbulence – has merit in terms of its realization at a range of Reynolds numbers.

Mechanics↗

The effects of turbulence on the flame structure and NO formation of ammonia turbulent premixed combustion at various equivalence ratios

Ammonia is carbon-free and is regarded as a potential fuel to address global warming issues. Here, in this work, three-dimensional direct numerical simulations (DNS) of ammonia/air turbulent premixed flames were performed to explore the influence of turbulence and equivalence ratio on the flame structure and NO formation characteristics. Two equivalence ratios were considered, i.e. Φ=0.9 and Φ=1.1. The general flame structures were presented and species distributions were examined. The NO mass fraction was found to be the highest in the product for the lean case and in the reaction zone for the rich case. The conditional mean values of species mass fractions and reaction rates were compared with those of the unstrained and strained laminar premixed flames to explore how well the laminar flame structures can approach those of the turbulent flame. The budget analysis of the species transport equations showed that turbulent diffusion plays an important role in species transport. The turbulent diffusivity D T was estimated using the gradient hypothesis based on the DNS data. Various laminar flame simulations with different diffusivities were carried out. It was shown that the conditional means of the DNS agree well with those of the laminar flames including D T in the transport property calculation. The global and local NO formation characteristics were investigated. It was found that the mean NO production conditioned on the progress variable is lower compared with the corresponding laminar flame in the rich case. However, the relative contributions from various NO pathways are rarely affected by turbulence. The NO mass fraction is higher in negative curvature regions compared with positive curvature regions of the flame surface for the rich case, which is due to the preferential diffusion of H 2 and other radicals and the enhanced NO pathways in negatively curved regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗

High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility

Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.

buoyant mixing↗

Effects of gravity on sheared and nonsheared turbulent nonpremixed flames

The present numerical study is concerned with the fundamental physics of the multiway interaction between turbulence, chemical reaction, and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous, three-dimensional governing equations. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration. The effects of buoyancy include the generation of local shear, baroclinic production or destruction of vorticity, and countergradient heat and mass transport. Increased vorticity and small-scale turbulence promote further mixing and reaction. However, if the strain-rates become too high, local flame extinction can occur. Our objective is to gain an understanding of the complex interactions between the physical phenomena involved, with particular attention to the effects of buoyancy on the turbulence structure, flame behavior, and factors influencing flame extinction.

Elghobashi, Said↗

Numerical Studies of Turbulent Supersonic Plane-Channel Flows

This two-part talk begins with an overview of the `turbulence problem’, the computational strategies used to address it, and the role of direct numerical simulation (DNS). The second part presents DNS results to determine the turbulent Prandtl number Pr t above cold (isothermal) and hot (adiabatic) walls in a family of low-supersonic channel flows. A range of mean temperature/density variations, corresponding to effective/edge Mach numbers between 1.1 to 2.2, and wall-variable-based Reynolds number Reτw from 73 to 3800, is considered. The adiabatic condition is a new feature of special interest. The value of Pr t away from the wall approaches 0.85 above both the isothermal and adiabatic walls. The variations of the near-wall Pr t profiles in both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗wall scaling proposed by Huang et al. (1995), with only a weak dependence on R eτw . This leads to a rather simple proposal for a model of heat transfer, attached to an eddy-viscosity model.

Direct numerical simulation↗

Energy Transfer in Rotating Turbulence

The influence or rotation on the spectral energy transfer of homogeneous turbulence is investigated in this paper. Given the fact that linear dynamics, e.g. the inertial waves regime tackled in an RDT (Rapid Distortion Theory) fashion, cannot Affect st homogeneous isotropic turbulent flow, the study of nonlinear dynamics is of prime importance in the case of rotating flows. Previous theoretical (including both weakly nonlinear and EDQNM theories), experimental and DNS (Direct Numerical Simulation) results are gathered here and compared in order to give a self-consistent picture of the nonlinear effects of rotation on tile turbulence. The inhibition of the energy cascade, which is linked to a reduction of the dissipation rate, is shown to be related to a damping due to rotation of the energy transfer. A model for this effect is quantified by a model equation for the derivative-skewness factor, which only involves a micro-Rossby number Ro(sup omega) = omega'/(2(OMEGA))-ratio of rms vorticity and background vorticity as the relevant rotation parameter, in accordance with DNS and EDQNM results fit addition, anisotropy is shown also to develop through nonlinear interactions modified by rotation, in an intermediate range of Rossby numbers (Ro(omega) = (omega)' and Ro(omega)w greater than 1), which is characterized by a marco-Rossby number Ro(sup L) less than 1 and Ro(omega) greater than 1 which is characterized by a macro-Rossby number based on an integral lengthscale L and the micro-Rossby number previously defined. This anisotropy is mainly an angular drain of spectral energy which tends to concentrate energy in tile wave-plane normal to the rotation axis, which is exactly both the slow and the two-dimensional manifold. In Addition, a polarization of the energy distribution in this slow 2D manifold enhances horizontal (normal to the rotation axis) velocity components, and underlies the anisotropic structure of the integral lengthscales. Finally is demonstrated the ability of a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) model to predict the underlying spectral transfer structure and all the subsequent developments of classic anisotropy indicators in physical space, when compared to recent LES results. Even if the applications mainly concern developed strong turbulence, a particular emphasis is given to the strong formal analogy of this EDQNM2 model with recent weakly nonlinear approaches to wave-turbulence.

Cambon, Claude↗

Part 2 of a Computational Study of a Drop-Laden Mixing Layer

This second of three reports on a computational study of a mixing layer laden with evaporating liquid drops presents the evaluation of Large Eddy Simulation (LES) models. The LES models were evaluated on an existing database that had been generated using Direct Numerical Simulation (DNS). The DNS method and the database are described in the first report of this series, Part 1 of a Computational Study of a Drop-Laden Mixing Layer (NPO-30719), NASA Tech Briefs, Vol. 28, No.7 (July 2004), page 59. The LES equations, which are derived by applying a spatial filter to the DNS set, govern the evolution of the larger scales of the flow and can therefore be solved on a coarser grid. Consistent with the reduction in grid points, the DNS drops would be represented by fewer drops, called computational drops in the LES context. The LES equations contain terms that cannot be directly computed on the coarser grid and that must instead be modeled. Two types of models are necessary: (1) those for the filtered source terms representing the effects of drops on the filtered flow field and (2) those for the sub-grid scale (SGS) fluxes arising from filtering the convective terms in the DNS equations. All of the filtered-sourceterm models that were developed were found to overestimate the filtered source terms. For modeling the SGS fluxes, constant-coefficient Smagorinsky, gradient, and scale-similarity models were assessed and calibrated on the DNS database. The Smagorinsky model correlated poorly with the SGS fluxes, whereas the gradient and scale-similarity models were well correlated with the SGS quantities that they represented.

Okongo, Nora↗

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Direct numerical simulation of incompressible axisymmetric flows

In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Loulou, Patrick↗

Effects of Gravity on Sheared and Nonsheared Turbulent Nonpremixed Flames

The present numerical study is concerned with the fundamental physics of the multi-way interaction between turbulence, chemical reaction and buoyancy in a nonpremixed flame. The method of direct numerical simulation (DNS) is used to solve the instantaneous three-dimensional governing equations (continuity, Navier-Stokes, species mass fractions and energy) under the zero-Mach number assumption. Because of the present supercomputer limitations, we consider two simple flow geometries, namely an initially uniform flow without shear (equivalent to grid-generated turbulence) and an initially uniform shear flow. In each flow, the fuel and oxidant initially exist as two separate streams. As the reactants mix, chemical reaction takes place and exothermic energy is released causing variations in density. In the presence of a gravity field, the spatial and temporal distributions of the induced buoyancy forces depend on the local density gradients and the direction of the gravitational acceleration.

Elghobashi, Said↗

Proper orthogonal decomposition based turbulence modeling

The proper orthogonal decomposition is utilized to provide one mode estimates of various turbulence quantities in the near-wall region of the turbulent boundary layer. The experimental data of Herzog (1986) and the direct numerical simulation (DNS) results of Zang (1991) proved the two-point correlation tensor. The correlation tensor is then used as the kernel in the integral eigenvalue problem from which the proper orthogonal modes are obtained. The one mode estimates of the turbulent kinetic energy, dissipation, and triple moments all exhibit the proper asymtotic behavior near the wall. As a validation test, a model for the triple moments is developed and its near-wall behavior is studied. This information can then be used to obtain an estimate of the damping function near the wall, which is then compared to known existing behavior. Since this is an initial investigation, the streamwise and spanwise dependence is surpressed and only the wall normal direction (y+) behavior of various turbulence statistics near the wall is studied. Before examining one mode estimates of the triple moments and the dissipation (on which the model is based), various second order moments are examined. If these do not exhibit the proper behavior near the wall (as discussed by Speziale et al 1990) it does not seem reasonable to expect the one mode estimate of the triple moment and dissipation to have the right behavior. Only the results obtained using the experimental data of Herzog are discussed since the results from the simulations of Zang exhibit similar behavior as discussed above.

Glauser, Mark N.↗

Some observations concerning "laminarization" in heated vertical tubes

For low "turbulent" Reynolds numbers in strongly-heated vertical gas flow, fundamental results of existing direct numerical simulation (DNS) databases have been examined to deduce differences between cases which "revert" to turbulent and those that yield integral parameters which correspond to laminar flow (henceforth laminarizing or laminarized). The objectives of this paper are (1) to examine the streamwise evolution of near-wall flow structures and (2) to determine which, if any, proposed laminarization parameters could be used to discriminate between turbulent and laminarizing flows resulting. Views of streamlines in r-θ cross sections showed evidence of a ring of irregular-shaped streamwise vortices near the wall at all locations. The transverse advection of these vortices is hypothesized to provide a path of least resistance to the transport of streamwise momentum in the wall-normal direction. It is demonstrated that a steady streamwise laminar vortex can provide apparent turbulent momentum transport (as a consequence of the Reynolds decomposition) such that a reasonable mean velocity profile is predicted near the wall. For the present cases, the values of the non-dimensional radius (y c,w + = Re τ,w ) and of the modified wall Reynolds number did discriminate whether turbulent or laminarizing flow resulted.

99 GENERAL AND MISCELLANEOUS↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary-layer transition↗

Direct Numerical Simulation of hydrogen combustion at auto-ignitive conditions: Ignition, stability and turbulent reaction-front velocity

Direct Numerical Simulations (DNS) are performed to investigate the process of spontaneous ignition of hydrogen flames at laminar, turbulent, adiabatic and non-adiabatic conditions. Mixtures of hydrogen and vitiated air at temperatures representing gas-turbine reheat combustion are considered. Adiabatic spontaneous ignition processes are investigated first, providing a quantitative characterization of stable and unstable flames. Results indicate that, in hydrogen reheat combustion, compressibility effects play a key role in flame stability and that unstable ignition and combustion are consistently encountered for reactant temperatures close to the mixture’s characteristic crossover temperature. Furthermore, it is also found that the characterization of the adiabatic processes is also valid in the presence of non-adiabaticity due to wall heat-loss. Finally, a quantitative characterization of the instantaneous fuel consumption rate within the reaction front is obtained and of its ability, at auto-ignitive conditions, to advance against the approaching turbulent flow of the reactants, for a range of different turbulence intensities, temperatures and pressure levels.

42 ENGINEERING↗

EDQNM closure: A homogeneous simulation to support it. A quasi-homogeneous simulation to disprove it

It is known that two-point closures are useful tools for understanding and predicting turbulence. Among the various closures, the Eddy Damped Quasi-Normal Markovian (EDQNM) approach is one of the simplest and, at the same time, most useful. Direct numerical simulations (DNS) can provide information that can be used to test the validity of two-point theories. It is the purpose of the present work to use DNS to validate, or improve upon, EDQNM. A case was selected for which EDQNM is known to give satisfactory results: homogeneous isotropic turbulence. Quantities were then evaluated which may be used to test the assumptions of two-point closure approximations: spectral Lagrangian time scales. The goal was to make a careful and refined study to validate the EDQNM theory. A reference case was built for which EDQNM is likely to give poor results. An attempt to generate a quasi-homogeneous turbulent field containing organized structures, was built by artifically injecting them in the initial conditions. The results of direct simulations using such initial conditions are expected to provide a challenge for EDQNM since this kind of field is simple enough to allow comparisons with two-point theories, but at the same time contains coherent structures which cannot be expected to be accurately accounted for by closures based on expansions about Gaussianity.

Bertoglio, J. P.↗

The production of premixed flame surface area in turbulent shear flow

In the present work, we use three-dimensional Direct Numerical Simulation (DNS) of premixed flames in turbulent shear flow to characterize the effect of a mean shear motion on flame surface production. The shear is uniform in the unburnt gas, and simulations are performed for different values of the mean shear rate, S. The data base is then used to estimate and compare the different terms appearing in the Sigma-equation as a function of S. The analysis gives in particular the relative weights f the turbulent flow and mean flow components, a(sub T) and A(sub T), of the flame surface production term. This comparison indicates whether the dominant effects of a mean flow velocity gradient on flame surface area are implicit and scale with the modified turbulent flow parameters, kappa and epsilon, or explicit and scale directly with the rate of deformation.

Trouve, A.↗

Understanding the Role of CFD, Including RANS, DES, LES and DNS, in Hypersonic Propulsion Design and Analysis

Computational fluid dynamics (CFD) has matured to the point where it is becoming ubiquitous in the design and analysis of aerospace vehicles and their propulsion systems. The current state-of-the-art of CFD in hypersonic propulsion design and analysis as well as a historical perspective of its development over the last 50 years is presented. The current state-of-the-art for the design of finite rate chemically reacting turbulent hypersonic flow is shown to currently be Reynolds averaged simulations (RAS) where the Reynolds averaged Navier-Stokes (RAS) equations are solved to obtain steady state solutions. An overview of the physical modeling techniques employed by RAS for the simulation of hypersonic chemically reacting flow are presented. However, scale resolving simulation techniques, such as detached eddy simulation (DES), large eddy simulations (LES), hybrid RAS/LES and direct numerical simulations (DNS), are becoming becoming powerful tools for point design analysis. One example of such a scale resolving technique, hybrid RAS/LES, will be briefly described as well as how it can be used to help address model form uncertainties in RAS.

Hypersonic↗