Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Turbulence simulations”

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

A Numerical Study of Self-Similarity in a Turbulent Plane Wake Using Large-Eddy Simulation

Turbulent wakes are known to develop self-similarly sufficiently far downstream from obstacles that generate them. It has long been assumed that the spreading rate of the wake in the self-similar regime is independent of the details of the body generating the wake, being dependent only on the total drag (or momentum deficit). This assumption seems to be in contradiction with some recent experiments. In this study we attempt to complement these experimental investigations through a numerical study of a time-developing wake. A numerical study has the advantage of eliminating many of the uncontrolled factors present in experiments and allowing precise control of initial conditions. Large-eddy simulations employing the recently developed dynamic localization model are used to extend previous results from direct numerical simulations. The large-eddy simulation results are compared to the direct numerical simulation database, wherever such comparisons are feasible, as a check of the method. Like the experiments, the large-eddy simulations suggest that non-unique self-similar states, characterized by different spreading rates and turbulent statistics, are possible and that they can be maintained for significant time periods. The study also demonstrates the predictive capability of the dynamic localization subgrid model.

Ghosal, Sandip

Study of Lagrangian characteristic times using direct numerical simulation of turbulence

Direct numerical simulations and large eddy simulations of homogeneous isotropic turbulence are used to compute Lagrangian statistics of turbulence and, in particular, its time scales. The computed time scales are compared with the spectral times scales that are frequently used in Eddy Damped Quasi-Normal Markovian calculations of the spectrum. The time scale models are rather good at high wavenumber and the results point to directions for improvement of the time scales at low wavenumber.

Lee, C. H.

Advances in turbulent transport modeling based on direct simulations of turbulence

Recent work using numerical simulations to study the modeling of scalar transport is outlined. A model developed by reference to simulations for homogeneous shear flow is presented. This model predicts cross-gradient transports that can give rise to the illusion of counter-gradient transport. The model is found to work remarkably well when compared to simulations for channel flow, but fails to predict important features of the scalar fluxes in mixing layers. Reasons for this failure are discussed and suggestions are made for future work.

Reynolds, W. C.

Development and application of an atmospheric turbulence model for use in flight simulators in flight simulators

The influence of simulated turbulence on aircraft handling qualities was investigated. Pilot opinion of the handling qualities of a light general aviation aircraft were evaluated in a motion-base simulator using a simulated turbulence environment. A realistic representation of turbulence disturbances is described in terms of rms intensity and scale length and their random variations with time. The time histories generated by the proposed turbulence models showed characteristics which appear to be more similar to real turbulence than the frequently-used Gaussian turbulence model. In addition, the proposed turbulence models can flexibly accommodate changes in atmospheric conditions and be easily implemented in flight simulator studies. Six turbulence time histories, including the conventional Gaussian model, were used in an IFR-tracking task. The realism of each of the turbulence models and the handling qualities of the simulated airplane were evaluated. Analysis of pilot opinions shows that at approximately the same rms intensities of turbulence, the handling quality ratings transit from the satisfactory level, for the simple Gaussian model, to an unacceptable level for more realistic and compositely structured turbulence models.

Jacobson, I. D.

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble averaged results for a widely used 1-D turbulence approximation.

J. Riaz

A new method for simulating atmospheric turbulence for rotorcraft applications

Simulation of atmospheric turbulence as seen by a rotating blade element involves treatment of cyclostationary processes. Conventional filtering techniques do not lend themselves well to the generation of such turbulence sample functions as are required in rotorcraft flight dynamics simulation codes. A method to generate sample functions containing second-order statistics of mean and covariance is presented. Compared to ensemble averaging involving excessive computer time, the novelty is to exploit cycloergodicity and thereby, replace ensemble averaging by averaging over a single-path sample function of long duration. The method is validated by comparing its covariance results with the analytical and ensemble-averaged results for a widely used one-dimensional turbulence approximation.

Prasad, J. V. R.

Numerical simulations of turbulent premixed H2/O2/N2 flames with complex chemistry

Premixed stoichiometric H2/O2/N2 flames propagating in two-dimensional turbulence were studied using direct numerical simulation (simulations in which all fluid and thermochemical scales are fully resolved) including realistic chemical kinetics and molecular transport. Results are compared with earlier zero-chemistry (flame sheet) and one-step chemistry simulations. Consistent with the simpler models, the turbulent flame with realistic chemistry aligns preferentially with extensive strain rates in the tangent plane and flame curvature probability density functions are close to symmetric with near-zero means. By contrast to simple-chemistry results with non-unity Lewis numbers (ratio of thermal to species diffusivity), local flame structure does not correlate with curvature but rather with tangential strain rate. Turbulent straining results in substantial thinning of the flame relative to the steady unstrained laminar case. Heat release and H2O2 contours remain thin and connected ('flamelet-like') while species including H-atom and OH are more diffuse. Peak OH concentration occurs well behind the peak heat-release zone. The feasibility of incorporating realistic chemistry into full turbulence simulations to address issues such as pollutant formation in hydrocarbon-air flames is suggested.

Baum, M.

Simulations of turbulent mixing and reacting flows and their applications to turbulence modeling

The method of full simulation is applied to reacting turbulent flows. Full simulation has proven of great value as a complement to experiments for the study of nonreacting turbulent flows. It provides insight into the physics of turbulent flows and their modeling. It is natural to try to extend these methods to the simulation of reacting turbulent flows. Because this is one of the first attempts at this type of simulation, a subsidiary goal of this work is to demonstrate the feasibility of using simulation to study turbulent reacting flows. In addition, it is shown that such simulations can be used to provide physical insight into the nature of turbulent combustion and to provide data that will help to construct models that can be used in engineering simulations of turbulent reacting flows.

Ferziger, J. H.

Simulations of transition and turbulence on the Navier-Stokes computer

The Navier-Stokes Computer (NSC) consists of multiple local memory parallel processors interconnected in a hypercube network. Efficient implementation of algorithms on the NSC thus requires the effective utilization of both the coarse and fine grain paralelism inherent in the architectural design. The basic approach to implementing an algorithm on the NSC is presented herein. The particular finite-difference algorithm considered was developed for performing transition and turbulence simulations by direct solution of the time-dependent incompressible Navier-Stokes equations. The suitability of this algorithm for performing simulations of the isotropic turbulence problem is verified from computations performed on a Cray 2. Projected timing results for the algorithm on the NSC itself are presented for both the isotropic turbulence and laminar turbulent transition problems.

Krist, S. E.

Beam stresses induced by a turbulent boundary layer and simulated with a reverberant acoustic field

Unsteady pressure fluctuations are a feature of the flow field of an airplane or a space vehicle in atmospheric flight. Undesirable effects of these fluctuations range from internal (cabin) noise to structural fatigue and damage of avionic and ordnance systems. The integrity to withstand these fluctuating loads is established through reverberant room acoustical qualification testing. The purpose of this paper is to develop a methodology of simulating turbulent boundary layer fluctuating pressure loading for a simply supported beam with a reverberant acoustic field. This goal was accomplished using normal mode vibration analysis. The main results were developed in terms of dimensionless quantities such as turbulent boundary layer spectrums, reverberant acoustical field simulation spectrums, and stress response spectrums with parameters: dimensionless fundamental frequency, Mach number, and relative boundary layer displacement thickness.

Schutzenhofer, L. A.

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan

Turbulence-Resolving Simulations of Hurricane Laura (2020): Insights Into Extreme Winds and Eyewall Turbulence

Extreme weather events pose significant risks to coastal and offshore energy infrastructure. In this work, we investigate the structure of mean winds and turbulence near the surface ( z < 300 m) that occurred during Hurricane Laura in 2020 on the US Gulf Coast. To this end, we perform turbulence-resolving simulations (..delta..x = 33.33 m) that encompass the entire eyewall of Hurricane Laura by introducing a novel method that we will refer to as Moving-Mesoscale to Static-LES Integrated Coupling (MOSAIC). The simulation results are validated against surface, boundary-layer, and mid-tropospheric observations. Our analysis of the extreme wind conditions near the surface indicates that the mean wind and turbulence profiles vary greatly at and near the eyewall, resulting in extreme values of mean wind speed (U > 50 m * s-1) and turbulence kinetic energy (k ~ 50 m2 * s-2) at altitudes above 50 m . To investigate in detail the nature of turbulence within the eyewall of the storm, we conduct higher-resolution simulations (..delta..x = 11.11 m) of a portion of the hurricane. We provide a comprehensive statistical description of turbulence in the eyewall region, highlighting the need to employ ..delta..x ~ 10 m grid spacing to capture the integral length-scales near the surface, second- and higher-order central moments (i.e., skewness and kurtosis), and spectral coherence in the flow accurately. These numerical simulations provide the most detailed description of mean winds and turbulent conditions within the eyewall of a historical tropical cyclone to date, illustrating how high-resolution simulations can enhance limited turbulence measurements in hurricanes and giving insight into the extreme conditions threatening offshore and coastal infrastructure.

17 WIND ENERGY

Simulation of turbulent wall pressure

A Monte Carlo procedure was developed to simulate turbulent boundary layer wall pressure fluctuations. The approach utilizes much of the newly available conditional sampling information to construct the required distribution functions. Various disturbance wave forms were examined, as well as the effect of frequency-dependent decay. Good agreement between the simulation and experimental data was achieved for root mean square pressure level, power spectrum, and space time correlation.

Ash, R. L.

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.

Large eddy numerical simulations of turbulent flows

Large eddy simulations are a numerical technique in which large scale turbulent structures are explicitly computed and the small structures are modelled. Arguments for believing this method to be superior to more conventional approaches are given, the basis of the method is given, and some typical results displayed. The results show that the method does have enormous promise, but much further development is required.

Ferziger, J. H.

Simulation of turbulent shear flows at Stanford and NASA-Ames - What can we do and what have we learned?

The capabilities and limitations of large eddy simulation (LES) and full turbulence simulation (FTS) are outlined. It is pointed out that LES, although limited at the present time by the need for periodic boundary conditions, produces large-scale flow behavior in general agreement with experiments. What is more, FTS computations produce small-scale behavior that is consistent with available experiments. The importance of the development work being done on the National Aerodynamic Simulator is emphasized. Studies at present are limited to situations in which periodic boundary conditions can be applied on boundaries of the computational domain where the flow is turbulent.

Reynolds, W. C.

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.