Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulent Simulation”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Asymptotic Dynamics of High Dynamic Range Stratified Turbulence

Direct numerical simulations of homogeneous sheared and stably stratified turbulence are considered to probe the asymptotic high-dynamic range regime suggested by Gargett et al. [1] and Shih et al. [2]. We consider statistically stationary configurations of the flow that span three decades in dynamic range defined by the separation between the Ozmidov length scale, L O = √ ε/N 3 , and the Kolmogorov length scale, L K = (ν 3 /ε) 1/4 , up to Re b ≡ (L O /L K ) 4/3 = ε/(νN 2 ) ~ O(1000), where ε is the mean turbulent kinetic energy dissipation rate, ν is the kinematic viscosity, and N is the buoyancy frequency. We isolate the effects of Reb, particularly on irreversible mixing, from the effects of other flow parameters of stratified and sheared turbulence [3]. Specifically, we evaluate the influence of dynamic range independent of initial conditions. We present evidence that the flow approaches an asymptotic state for Re b &#10886 300, characterized both by an asymptotic partitioning between the potential and kinetic energies and by the approach of components of the dissipation rate to their expected values under the assumption of isotropy. As Re b increases above 100, there is a slight decrease in the turbulent flux coefficient Γ = χ/ε, where χ is the dissipation rate of buoyancy variance, but, for this flow, there is no evidence of the commonly suggested Γ ∝ Re b -1/2 dependence when 100 ≤ Re b ≤ 1000.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Subgrid-scale backscatter in transitional and turbulent flows

Most subgrid-scale (SGS) models for large-eddy simulations are absolutely dissipative (that is, they remove energy from the large scales at each point in the physical space). The actual SGS stresses, however, may transfer energy to the large scales (backscatter) at a given location. Direct numerical simulations of turbulent channel flow and compressible isotropic turbulence are used to study the backscatter phenomena. In all flows considered roughly 50 percent of the grid points were experiencing backscatter when a Fourier cutoff filter was used. The backscatter fraction was less with a Gaussian filter, and intermediate with a box filter in physical space. Moreover, the backscatter and forward scatter contributions to the SGS dissipation were comparable, and each was often much larger than the total SGS dissipation. The SGS dissipation (normalized by total dissipation) increased with filter width almost independently of filter type and Reynolds number. The amount of backscatter showed an increasing trend with Reynolds numbers. In the near-wall region of the channel, events characterized by strong Reynolds shear stress correlated fairly well with areas of high SGS dissipation (both forward and backward). In compressible isotropic turbulence similar results were obtained, independent of fluctuation Mach number.

Piomelli, Ugo↗
NASA NTRS1990

On the Subgrid-Scale Modeling of Compressible Turbulence

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Squires, Kyle↗
DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗
NASA NTRSDate not supplied

Turbulence Effects on Shaped Booms: Propagation Simulations Using KZKFourier

Upcoming X-59 aircraft flight tests as part of NASA’s Quesst Mission are expected to occur in a range of atmospheric conditions. Sensitivity of ground waveform acoustic metrics to turbulent perturbations during propagation complicates the determination of noise levels. A series of simulations through turbulence was executed using the KZKFourier model of Stout et al. to develop a database of results for expanding the functionality of NASA tools for estimating turbulence effects on shaped-boom ground waveforms. For 45 cases covering a seven-factor design space, multiple realizations of turbulence were generated to characterize statistical turbulence effects on levels of six acoustic metrics. An individual simulation in the study produced over one thousand waveforms across a virtual microphone array. Different measures were used to evaluate refinement of results with increasing number of simulations. Data suggest that mean effects and variability were most strongly influenced by propagation distance and velocity fluctuation intensity. Although local increases in acoustic metrics were common, the overall average result in all cases was a reduction in metric levels. Sensitivity of acoustic metrics varied, with mean reductions of up to 2 dB in Perceived Level occurring across the 45 cases.

turbulence↗
NASA NTRS1991

Simulation and modeling of homogeneous compressible turbulence under isotropic mean compression

Compressible homogeneous turbulence subjected to isotropic mean compression is simulated numerically for high and low turbulent Mach numbers, at various compression rates. It is found that at low Mach numbers, the effects of viscosity variations on the development of the turbulent kinetic energy can be significant. A model consistent with the invariance of the Navier-Stokes equations to spherical compression, that takes into account variable viscosity effects is proposed. At high Mach numbers, the distribution of energy between the acoustic and solenoid fields is found to depend upon initial conditions. The simulation results are also used to evaluate two-equation compressible turbulence models.

Coleman, G. N.↗
DOE OSTI2021

Investigation of grid-based vorticity-velocity large eddy simulation off-body solvers for application to overset CFD

Accurately predicting unsteady wakes and vortex-dominated flows is essential to a wide range of engineering applications, including aircraft, rotorcraft, shipboard operations, bio-inspired unsteady flight and propulsion, wind turbines, and urban flows. While current CFD software can model the complete flow field and wake system, the computational costs incurred in high Reynolds number unsteady turbulent flow simulations often remain prohibitive for routine engineering use, particularly for applications involving moving components. Prior work has demonstrated that by adopting a vorticity-velocity formulation in a grid-based off-body flow solver (VorTran-M and VorTran-M2) one can lower these costs by several orders of magnitude when compared to conventional approaches. This paper describes the extensions made to VorTran-M2 to support turbulent flows, and associated benchmarking activity to assess its performance for problems involving strong stretching and diffusion processes, whose competing contributions to the vorticity field are core drivers of turbulent flow evolution. Predictions are presented for: (i) the Kida-Pelz problem whose inviscid form is of mathematical interest due to its apparent formation of singular flow in finite time; and (ii) the Taylor Green vortex arrangement, which has been extensively studied as a fundamental simulation challenge in the turbulent modeling community. Here, the results are used to evaluate the overall predictive ability and performance of two sub-grid scale models incorporated into VorTran-M2. Results indicate that the computational cost savings seen previously for inviscid and convection dominated problems extend to turbulent flow simulations supporting the viability of VorTran-M2 as a low cost means for accurately modeling the far-field and background flow, particularly when long duration vorticity evolution is of interest.

42 ENGINEERING↗
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
DOE OSTI2020

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

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Development and application of a non-Gaussian atmospheric turbulence model for use in flight simulators

A method is described for generating time histories which model the frequency content and certain non-Gaussian probability characteristics of atmospheric turbulence including the large gusts and patchy nature of turbulence. Methods for time histories using either analog or digital computation are described. A STOL airplane was programmed into a 6-degree-of-freedom flight simulator, and turbulence time histories from several atmospheric turbulence models were introduced. The pilots' reactions are described.

Reeves, P. M.↗
DOE OSTI2013

VLASOV SIMULATIONS OF MULTI-ION PLASMA TURBULENCE IN THE SOLAR WIND

Hybrid Vlasov-Maxwell simulations are employed to investigate the role of kinetic effects in a two-dimensional turbulent multi-ion plasma, composed of protons, alpha particles, and fluid electrons. In the typical conditions of the solar-wind environment, and in situations of decaying turbulence, the numerical results show that the velocity distribution functions of both ion species depart from the typical configuration of thermal equilibrium. These non-Maxwellian features are quantified through the statistical analysis of the temperature anisotropy, for both protons and alpha particles, in the reference frame given by the local magnetic field. Anisotropy is found to be higher in regions of high magnetic stress. Both ion species manifest a preferentially perpendicular heating, although the anisotropy is more pronounced for the alpha particles, according to solar wind observations. The anisotropy of the alpha particle, moreover, is correlated to the proton anisotropy and also depends on the local differential flow between the two species. Evident distortions of the particle distribution functions are present, with the production of bumps along the direction of the local magnetic field. The physical phenomenology recovered in these numerical simulations reproduces very common measurements in the turbulent solar wind, suggesting that the multi-ion Vlasov model constitutes a valid approach to understanding the nature of complex kinetic effects in astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Beta dependence of electron heating in decaying whistler turbulence: Particle-in-cell simulations

Two-dimensional particle-in-cell simulations have been carried out to study electron beta dependence of decaying whistler turbulence and electron heating in a homogeneous, collisionless magnetized plasma. Initially, applied whistler fluctuations at relatively long wavelengths cascade their energy into shorter wavelengths. This cascade leads to whistler turbulence with anisotropic wavenumber spectra which are broader in directions perpendicular to the background magnetic field than in the parallel direction. Comparing the development of whistler turbulence at different electron beta values, it is found that both the wavenumber spectrum anisotropy and electron heating anisotropy decrease with increasing electron beta. This indicates that higher electron beta reduces the perpendicular energy cascade of whistler turbulence. Fluctuation energy dissipation by electron Landau damping responsible for the electron parallel heating becomes weaker at higher electron beta, which leads to more isotropic heating. It suggests that electron kinetic processes are important in determining the properties of whistler turbulence. This kinetic property is applied to discuss the generation of suprathermal strahl electron distributions in the solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Numerical investigation of supersonic turbulent boundary layers with high wall temperature

A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Guo, Y.↗
DOE OSTI2012

FY12 Report: Meteorological Observations for Renewable Energy Applications at Site 300

Accurate, high-resolution meteorological measurements of wind flow in the planetary boundary layer (PBL) and surface-atmosphere energy exchange are required for understanding the properties and quality of available wind power at Site 300. Wind speeds at heights found in a typical wind turbine rotor disk (~ 40- 120 m) are driven by the synergistic impacts of atmospheric stability, orography, and land-surface characteristics on the mean wind flow and related turbulence structures in the PBL. In FY11, two Laser and Detection Ranging (LIDAR) units, a 3 m tall flux tower, and a 3 m tall meteorological tower were installed in the northern section of Site 300 as a first step in development of a renewable energy testbed facility. The LIDAR units were programmed to measure wind speed and direction up to 200 m. The flux tower measured land surface-atmospheric energy exchanges and provided information on the mechanisms for turbulence production in the atmosphere and atmospheric stability. Atmospheric stability largely determines wind speed and turbulence profiles in the lower part of the boundary layer where wind turbines reside. In FY12 the SMS project was further expanded by the earlier addition of solar power to the LIDAR units. This enabled the units to be moved to different areas within Site 300 in order to gather a larger spatial range of measurements. Additional LIDAR locations were chosen based on (1) feasibility and logistics (e.g., away from the blast zone and special habitat zones) and (2) areas forecasted to have high wind speeds and/or high turbulence from the Weather Research and Forecasting (WRF) model simulations. High turbulence is of interest because very high amounts of turbulent energy can be damaging to turbine components. Between December 2011-April 2012, the two LIDAR units were rotated amongst three locations: the tall (52 m) meteorological tower, the radio tower, and the Matthew area (named after Matthew Simpson) (Figure 1). The radio tower location was in the center of Site 300 and represented the site with highest elevation and potentially highest wind speeds as forecasted by WRF. The Matthew area in the southeast represented a site with high wind speeds and potentially high levels of turbulence. In FY11 we presented a quality control report which illustrated the accuracy of the LIDAR wind speed and direction measurements as compared to the 52 m tall meteorological tower and provided annual statistics for wind speed and direction up to 200 m for the meteorological tower location in the northern section of Site 300. In the FY12 report we present: 1. a second year of LIDAR and flux tower statistics for 2012, 2. an enhanced spatial dataset of LIDAR observations with the addition of two new sites, and 3. a look at off-site LIDAR field campaigns.

17 WIND ENERGY↗
NASA NTRS1992

Computed Turbulent Free Shear Flow Of Air

Standard k-epsilon model of turbulence yields fairly accurate results. Symposium paper discusses numerical simulation of turbulent free shear flow of nonreacting compressible fluid. Ability to compute such flows essential to advances in design.

Viegas, J. R.↗
DOE OSTI2013

MODELING STATISTICAL PROPERTIES OF SOLAR ACTIVE REGIONS THROUGH DIRECT NUMERICAL SIMULATIONS OF 3D-MHD TURBULENCE

Statistical properties of the Sun's photospheric turbulent magnetic field, especially those of the active regions (ARs), have been studied using the line-of-sight data from magnetograms taken by the Solar and Heliospheric Observatory and several other instruments. This includes structure functions and their exponents, flatness curves, and correlation functions. In these works, the dependence of structure function exponents ({zeta}{sub p}) of the order of the structure functions (p) was modeled using a non-intermittent K41 model. It is now well known that the ARs are highly turbulent and are associated with strong intermittent events. In this paper, we compare some of the observations from Abramenko et al. with the log-Poisson model used for modeling intermittent MHD turbulent flows. Next, we analyze the structure function data obtained from the direct numerical simulations (DNS) of homogeneous, incompressible 3D-MHD turbulence in three cases: sustained by forcing, freely decaying, and a flow initially driven and later allowed to decay (case 3). The respective DNS replicate the properties seen in the plots of {zeta}{sub p} against p of ARs. We also reproduce the trends and changes observed in intermittency in flatness and correlation functions of ARs. It is suggested from this analysis that an AR in the onset phase of a flare can be treated as a forced 3D-MHD turbulent system in its simplest form and that the flaring stage is representative of decaying 3D-MHD turbulence. It is also inferred that significant changes in intermittency from the initial onset phase of a flare to its final peak flaring phase are related to the time taken by the system to reach the initial onset phase.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Comment on "Heterodyne Lidar Returns in the Turbulent Atmosphere: Performance Evaluation of Simulated Systems"

The explanation for the difference between simulation and the zero-order theory for heterodyne lidar returns in a turbulent atmosphere proposed by Belmonte and Rye is incorrect. The theoretical expansion is not developed under a square- law-structure function approximation (random wedge atmosphere). Agreement between the simulations and the zero-order term of the theoretical expansion is produced for the limit of statistically independent paths (bi-static operation with large transmitter-receiver separation) when the simulations correctly include the large-scale gradients of the turbulent atmosphere.

Frehlich, Rod↗
NASA NTRS1990

Physical considerations governing computation of turbulent flows over bodies at large incidence

The physical aspects governing accurate numerical simulation of turbulent flows having large regions of crossflow separation are re-examined. Time-accurate, three-dimensional fine-grid Navier-Stokes solutions were obtained for turbulent subsonic and supersonic flows over a slender ogive-cylinder body of revolution at large angles of attack. These flowfields are complex and contain regions of crossflow separation and an organized leeward-side vortex structure. An algebraic eddy-viscosity turbulence model has been modified to correctly account for the effects of the vortices on the underlying viscous layers. The numerical results show that the vortical flow structure of both subsonic and supersonic high incidence flows are similar, and can be accurately treated by the turbulence model. As the angle of attack is increased, the effectiveness of the model increases since the influence of the vortices on the underlying boundary layer decreases.

Degani, David↗