SEARCH · Engineering Papers

Results for “turbulence 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

The Influence of Adaptive Mesh Refinement on the Prediction of Vortex Interactions about a Generic Missile Airframe

The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.

slender body aerodynamics↗

DOE OSTI2026

Modification of ion-temperature-gradient turbulence by impurities in stellarator plasmas

Recent nonlinear gyrokinetic simulations have shown that impurities can strongly modify the turbulent heat flux in stellarator plasmas. Here, the ion-temperature-gradient (ITG) dispersion relation in a plasma containing impurities is analytically solved in certain limits and an expression for the modification of the ITG growth rate by impurities is derived. The analytical expression is the sum of three terms corresponding to three different physical causes (impurity density gradient, impurity temperature gradient and dilution) of the change in the growth rate. The scalings predicted analytically for the modification of the growth rate are shown to be reproduced by linear gyrokinetic simulations. The conditions for reduction or increase of the ITG growth by impurities are also correctly predicted by the analytical solution to the dispersion relation. Finally, a remarkable correlation is found between the analytical expression for the modification of the growth rate and the modification of the turbulent heat flux obtained from nonlinear gyrokinetic simulations.

Calvo, Iván [Research Centre for Energy, Environme↗

NASA NTRS1999

Streamwise Vorticity Generation in Laminar and Turbulent Jets

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Demuren, Aodeji O.↗

DOE OSTI2020

Simulating Real Atmospheric Boundary Layers at Gray-Zone Resolutions: How Do Currently Available Turbulence Parameterizations Perform?

Recent computational and modeling advances have led a diverse modeling community to experiment with atmospheric boundary layer (ABL) simulations at subkilometer horizontal scales. Accurately parameterizing turbulence at these scales is a complex problem. The modeling solutions proposed to date are still in the development phase and remain largely unvalidated. This work assesses the performance of methods currently available in the Weather Research and Forecasting (WRF) model to represent ABL turbulence at a gray-zone grid spacing of 333 m. We consider three one-dimensional boundary layer parameterizations (MYNN, YSU and Shin-Hong) and coarse large-eddy simulations (LES). The reference dataset consists of five real-case simulations performed with WRF-LES nested down to 25 m. Results reveal that users should refrain from coarse LES and favor the scale-aware, Shin-Hong parameterization over traditional one-dimensional schemes. Overall, the spread in model performance is large for the cellular convection regime corresponding to the majority of our cases, with coarse LES overestimating turbulent energy across scales and YSU underestimating it and failing to reproduce its horizontal structure. Despite yielding the best results, the Shin-Hong scheme overestimates the effect of grid dependence on turbulent transport, highlighting the outstanding need for improved solutions to seamlessly parameterize turbulence across scales.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2007

Tandem Cylinder Noise Predictions

In an effort to better understand landing-gear noise sources, we have been examining a simplified configuration that still maintains some of the salient features of landing-gear flow fields. In particular, tandem cylinders have been studied because they model a variety of component level interactions. The present effort is directed at the case of two identical cylinders spatially separated in the streamwise direction by 3.7 diameters. Experimental measurements from the Basic Aerodynamic Research Tunnel (BART) and Quiet Flow Facility (QFF) at NASA Langley Research Center (LaRC) have provided steady surface pressures, detailed off-surface measurements of the flow field using Particle Image Velocimetry (PIV), hot-wire measurements in the wake of the rear cylinder, unsteady surface pressure data, and the radiated noise. The experiments were conducted at a Reynolds number of 166 105 based on the cylinder diameter. A trip was used on the upstream cylinder to insure a fully turbulent shedding process and simulate the effects of a high Reynolds number flow. The parallel computational effort uses the three-dimensional Navier-Stokes solver CFL3D with a hybrid, zonal turbulence model that turns off the turbulence production term everywhere except in a narrow ring surrounding solid surfaces. The current calculations further explore the influence of the grid resolution and spanwise extent on the flow and associated radiated noise. Extensive comparisons with the experimental data are used to assess the ability of the computations to simulate the details of the flow. The results show that the pressure fluctuations on the upstream cylinder, caused by vortex shedding, are smaller than those generated on the downstream cylinder by wake interaction. Consequently, the downstream cylinder dominates the noise radiation, producing an overall directivity pattern that is similar to that of an isolated cylinder. Only calculations based on the full length of the model span were able to capture the complete decay in the spanwise correlation, thereby producing reasonable noise radiation levels.

Lockard, David P.↗

DOE OSTI2024

Planck Dust Polarization Power Spectra Are Consistent with Strongly Supersonic Turbulence

Abstract The polarization of the cosmic microwave background is rich in information but obscured by foreground emission from the Milky Way’s interstellar medium (ISM). To uncover relationships between the underlying turbulent ISM and the foreground power spectra, we simulated a suite of driven, magnetized, turbulent models of the ISM, varying the fluid properties via the sonic Mach number, M S , and magnetic (Alfvén) Mach number, M A . We measure the power spectra of density ( ρ ), velocity ( v ), magnetic field ( H ), total projected intensity ( T ), parity-even polarization ( E ), and parity-odd polarization ( B ). We find that the slopes of all six quantities increase with M S . Most increase with M A , while the magnetic field spectrum steepens with M A . By comparing spectral slopes of E and B to those measured by Planck, we infer typical values of M S and M A for the ISM. As the fluid velocity increases, M S > 4 , the ratio of BB power to EE power increases to approach a constant value near the Planck-observed value of ∼0.5, regardless of the magnetic field strength. We also examine correlation coefficients between projected quantities, and find that r TE ≈ 0.3, in agreement with Planck, for appropriate combinations of M S and M A . Finally, we consider parity-violating correlations r TB and r EB .

Astronomy & Astrophysics↗

NASA NTRS1991

The compressible mixing layer

The capability of turbulence modeling correctly to handle natural unsteadiness appearing in compressible turbulent flows is investigated. Physical aspects linked to the unsteadiness problem and the role of various flow parameters are analyzed. It is found that unsteady turbulent flows can be simulated by dividing these motions into an 'organized' part for which equations of motion are solved and a remaining 'incoherent' part represented by a turbulence model. Two-equation turbulence models and second-order turbulence models can yield reasonable results. For specific compressible unsteady turbulent flow, graphic presentations of different quantities may reveal complementary physical features. Strong compression zones are observed in rapid flow parts but shocklets do not yet occur.

Vandromme, Dany↗

NASA NTRS2011

Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow

The proposed paper reports advances in developing a method for high Reynolds number compressible viscous flow simulations using a Cartesian cut-cell method with embedded boundaries. This preliminary work focuses on accuracy of the discretization near solid wall boundaries. A model problem is used to investigate the accuracy of various difference stencils for second derivatives and to guide development of the discretization of the viscous terms in the Navier-Stokes equations. Near walls, quadratic reconstruction in the wall-normal direction is used to mitigate mesh irregularity and yields smooth skin friction distributions along the body. Multigrid performance is demonstrated using second-order coarse grid operators combined with second-order restriction and prolongation operators. Preliminary verification and validation for the method is demonstrated using flat-plate and airfoil examples at compressible Mach numbers. Simulations of flow on laminar and turbulent flat plates show skin friction and velocity profiles compared with those from boundary-layer theory. Airfoil simulations are performed at laminar and turbulent Reynolds numbers with results compared to both other simulations and experimental data

Berger, Marsha↗

DOE OSTI2025

Assessment and Improvement of the SST-Gamma Transition Model in Nalu-Wind

We conduct laminar–turbulent boundary-layer transition simulations using a local correlation-based transition model for two-dimensional incompressible flow and present enhancements to improve the accuracy of transition predictions. Menter’s Galilean-invariant 𝛾 transition model is implemented in the incompressible, unstructured-grid flow solver Nalu-Wind and is validated against experimental data and results from NASA’s flow solvers. The test cases of the AIAA Transition Prediction and Modeling Workshop are investigated, namely, the T3A/T3B flat plates and the NLF(1)-0416 and S809 airfoils. Based on the results, best practices for transition simulations, particularly for an unstructured-grid flow solver, are identified. Additional airfoil simulations are conducted for two wind turbine airfoils, S822 at Reynolds numbers of 𝑂⁡(10 5 ) and DU00-W-212 at Reynolds numbers of 𝑂⁡(10 7 ), to assess the model at low and high Reynolds numbers. Furthermore, through this work, we propose several approaches to enhance transition simulations, including 1) enforcing positivity of the implicit operator for the source terms of the transition model, 2) employing a constant turbulence intensity in stationary external flow simulations, and 3) recommending meshing for unstructured-grid flow solvers. Finally, we provide detailed documentation of the validation and data for the canonical cases to the transition modeling community.

17 WIND ENERGY↗

DOE OSTI2022

The role of the turbulence driving mode for the initial mass function

Turbulence is a critical ingredient for star formation, yet its role for the initial mass function (IMF) is not fully understood. Here we perform magnetohydrodynamical (MHD) simulations of star cluster formation including gravity, turbulence, magnetic fields, stellar heating, and outflow feedback to study the influence of the mode of turbulence driving on IMF. We find that simulations that employ purely compressive turbulence driving (COMP) produce a higher fraction of low-mass stars as compared to simulations that use purely solenoidal driving (SOL). The characteristic (median) mass of the sink particle (protostellar) distribution for COMP is shifted to lower masses by a factor of ∼1.5 compared to SOL. Our simulation IMFs capture the important features of the observed IMF form. We find that turbulence-regulated theories of the IMF match our simulation IMFs reasonably well in the high-mass and low-mass range, but underestimate the number of very low-mass stars, which form towards the later stages of our simulations and stop accreting due to dynamical interactions. Our simulations show that for both COMP and SOL, the multiplicity fraction is an increasing function of the primary mass, although the multiplicity fraction in COMP is higher than that of SOL for any primary mass range. We find that binary mass ratio distribution is independent of the turbulence driving mode. The average specific angular momentum of the sink particles in SOL is a factor of 2 higher than that for COMP. Overall, we conclude that the turbulence driving mode plays a significant role in shaping the IMF.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1998

Effect of Swirl on Turbulent Structures in Supersonic Jets

Direct numerical simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial.helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moser are observed. The jet core breakdown stage exhibits increased acoustic radiations.

Rao, Ram Mohan↗

NASA NTRS1998

Effect of Swirl on Turbulent Structures in Supersonic Jets

Direct Numerical Simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moserl are observed. The jet core breakdown stage exhibits increased acoustic radiations.

Rao, Ram Mohan↗

NASA NTRS1998

Effect of Swirl on Turbulent Structures in Supersonic Jets

Rao, Ram Mohan↗

DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: CCN Distribution

Abstract Following our previous investigation of the turbulence impact on cloud-base single-size CCN activation, this study considers a similar problem assuming CCN size distribution obtained from field measurements. The total CNN concentration is taken as either 200 cm −3 to represent clean conditions, or as 2000 cm −3 to represent polluted conditions. CCN is assumed to be sodium chloride. The CCN activation in the rising nonturbulent adiabatic parcel is contrasted with the activation within a rising adiabatic parcel filled with inertial-range homogeneous isotropic turbulence. The turbulent parcel of 64 3 m 3 and the turbulent kinetic energy dissipation rate of 10 −3 m −2 s −3 are used in most of the simulations. Results for a range of mean parcel ascent rates, between 0.125 and 8 m s −1 , are discussed. Overall, the adiabatic turbulent parcel simulations show results consistent with the adiabatic nonturbulent parcel, with higher activated CCN concentrations for stronger parcel ascent rates. The key difference is a blurriness of the separation between dry CCN size bins featuring activated and nonactivated (haze) CCN, especially for weak mean ascent rates. The blurriness comes from CCN getting activated and subsequently deactivated in the fluctuating supersaturation field, instead of all becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the nonturbulent parcel when activation is completed. Modeling results are discussed in the context of the impact of turbulent fluctuations on CCN activation documented in laboratory experiments using the Pi chamber.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗

DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Ion and Electron Acceleration in Fully Kinetic Plasma Turbulence

Turbulence is often invoked to explain the origin of nonthermal particles in space and astrophysical plasmas. By means of 3D fully kinetic particle-in-cell simulations, we demonstrate that turbulence in low-β plasmas (β is the ratio of plasma pressure to magnetic pressure) accelerates ions and electrons into a nonthermal energy distribution with a power-law energy range. The ion spectrum is harder than the electron one, and both distributions get steeper for higher β. We show that the energization of electrons is accompanied by a significant energy-dependent pitch-angle anisotropy, with most electrons moving parallel to the local magnetic field, while ions stay roughly isotropic. We demonstrate that particle injection from the thermal pool occurs in regions of high current density. Parallel electric fields associated with magnetic reconnection are responsible for the initial energy gain of electrons, whereas perpendicular electric fields control the overall energization of ions. Our findings have important implications for the origin of nonthermal particles in space and astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2013

Solutions of the Taylor-Green Vortex Problem Using High-Resolution Explicit Finite Difference Methods

A computational fluid dynamics code that solves the compressible Navier-Stokes equations was applied to the Taylor-Green vortex problem to examine the code s ability to accurately simulate the vortex decay and subsequent turbulence. The code, WRLES (Wave Resolving Large-Eddy Simulation), uses explicit central-differencing to compute the spatial derivatives and explicit Low Dispersion Runge-Kutta methods for the temporal discretization. The flow was first studied and characterized using Bogey & Bailley s 13-point dispersion relation preserving (DRP) scheme. The kinetic energy dissipation rate, computed both directly and from the enstrophy field, vorticity contours, and the energy spectra are examined. Results are in excellent agreement with a reference solution obtained using a spectral method and provide insight into computations of turbulent flows. In addition the following studies were performed: a comparison of 4th-, 8th-, 12th- and DRP spatial differencing schemes, the effect of the solution filtering on the results, the effect of large-eddy simulation sub-grid scale models, and the effect of high-order discretization of the viscous terms.

DeBonis, James R.↗