Engineering Papers⌕ Search

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.

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

At least 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS2023

Advances in Aircraft System Identification at NASA Langley Research Center

Advances in aircraft system identification at NASA Langley Research Center are discussed. The relevant time period includes the years since the last summary paper of this kind, which was published in the Journal of Aircraft in 2005. Research advances were achieved in flight test experiment design, frequency-domain modeling, real-time autonomous global modeling, rapid simulation development and updating, dynamic modeling in turbulence, flight data corrections, model uncertainty characterization, and aeroelastic modeling using distributed sensing. Possible future developments in the field are identified.

Aircraft system identification↗
NASA NTRS2024

Computational Techniques to Generate Space Launch System Aerodynamic Databases

This document describes the reasoning and trade studies used to evaluate tools for constructing the aerodynamic lineload databases for the liftoff and transition phases of flight for the Space Launch System. Three computational fluid dynamics codes (USM3D, FUN3D, Kestrel) were investigated with various turbulence models, as well as detached eddy simulation variants for the launch vehicle in free air and in proximity to the tower. Decisions were made mostly based on results from brief developmental studies performed in response to specific, unforeseen challenges that were encountered in the analysis of a given configuration. The need to develop databases in a timely manner, as well as accurately capture the expected leeward-wake flowfield characteristics, led to the selection of the Kestrel flow solver with its delayed detached eddy simulation method, the Spalart-Allmaras turbulence model, and the adaptive mesh refinement capability in the off-body Cartesian grid region.

Karen A. Deere↗
NASA NTRS2015

Higher-Order Methods for Compressible Turbulent Flows Using Entropy Variables

Turbulent flows have a large range of spatial and temporal scales which need to be resolved in order to obtain accurate predictions. Higher-order methods can provide greater efficiency for simulations requiring high spatial and temporal resolution, allowing for solutions with fewer degrees of freedom and lower computational cost than traditional second-order computational fluid dynamics (CFD) methods.1 Higher-order methods have been widely used for turbulent flows. However, the reduced numerical stabilization present in higher-order schemes implies that special care needs to be taken in the development of numerical methods to suppress nonlinear instabilities.2–6 In this work we present the development of a higher-order space-time discontinuous Galerkin method with a focus on the aspects of our numerical scheme required for ensuring nonlinear stability for turbulent simulations at high Reynolds numbers.

Diosady, Laslo T.↗
NASA NTRS1979

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional multi-element airfoils

The development of a numerical simulation of time dependent, turbulent, compressible flow about two dimensional multi-element airfoils of arbitrary shape is described. The basis of this simulation is a technique of automatic numerical generation of coordinate systems fitted to the multiple bodies regardless of their number or shape. Procedures developed whereby the coordinate lines are automatically concentrated in the boundary layer at any Reynolds number are discussed. The compressible turbulent solution involves an algebraic eddy viscosity turbulence model. The laminar version was run for transonic flow at free stream Mach numbers up to 0.9.

Thompson, J. F.↗
NASA NTRSDate not supplied

Split-Step Simulations of Sonic Boom Propagation Beyond the Lateral Cutoff in a Turbulent Atmosphere

Recent flight tests during the Quiet Supersonic Flights 2018 (QSF18) study reported sonic booms heard outside of the primary carpet region. In the absence of turbulence, the lateral cutoff region separates the primary sonic boom carpet from the shadow zone, where the sonic boom signal experiences significant attenuation. However, when turbulence is present in the atmospheric boundary layer (ABL), additional scattering of the sonic boom to the shadow zone region occurs. A method is presented for simulating sonic boom propagation in a turbulent atmospheric boundary layer beyond the lateral cutoff region into the shadow zone. A split-step method is used to integrate a partially one-way equation for the acoustic pressure. Inhomogeneous turbulence, representative of the ABL, is generated in the computational domain with a Fourier synthesis approach. Distributions of several loudness metrics in the shadow zone region for a sonic boom N-wave and a shaped boom are examined. Increasing both turbulence root-mean-square velocity and integral length scale are found to increase the average loudness of booms in the shadow zone. (This research is supported by the Commercial Supersonic Technology Project of the National Aeronautics and Space Administration under Grant No. 80NSSC19K1685.)

sonic boom↗
NASA NTRS1979

Computational simulation of corotational vortex merger using 0, 1, and 2 equation turbulence models

A computational study has been conducted to gain a better understanding of the role of turbulence in the merger of corotating vortices and in the decay of a single vortex. The flow field is simulated by using an inviscid flow model and zero-, one-, and two-equation turbulent-flow models. An implicit finite-difference procedure is employed to integrate the unsteady, two-dimensional, incompressible flow equations in a cross-plane. In the zero- and one-equation models, a spatially varying mixing-length model that simulates the streamline curvature effect is used to define an eddy viscosity. The computational merging distance criterion is in good agreement with the experimental one. In the case of the two-equation model, the turbulent kinetic energy and a modified rate of dissipation equation incorporating a curvature correction term are solved. A comparative study of the computational results corresponding to the various flow models and the available experimental data is presented.

Raj, P.↗
NASA NTRS1987

Numerical simulation of nonlinear interactions in channel and boundary-layer transition

The basic numerical techniques used for the direct simulation of transition to turbulence are described. Applications are presented for the computation of finite amplitude, two-dimensional steady states, the simulation of the secondary instability of channels and boundary-layers, and the simulation of the breakdown of detached shear layers.

Zang, T. A.↗
NASA NTRS1988

Helical magnetohydrodynamic turbulence and the coronal heating problem

Numerical simulations are used to investigate the relaxation of an unconfined, helically turbulent, fully three-dimensional magnetofluid, with conditions similar to those which are thought to result in the heating of the solar corona. In these simulations, the system evolves through a succession of force free states. After a relatively quiescent period of Ohmic decay, a phase of accelerated magnetic energy dissipation occurs. Some magnetic energy is transformed into kinetic energy, and the magnitude of entrophy created is a nontrival fraction of the mean square electric current. Concentrated vorticity structures are seen to play almost as important a role as electric current sheets in the heating process. Coincident with this accelerated dissipation process, a reorganization of the magnetic fields occurs, with transfer of magnetic energy to both shorter and longer wavelength modes than are initially present. The ratio of the magnetic field to the electric current density, alpha does not in general tend to assume a constant value in the force free regions during the evolution of the magnetofluid.

Dahlburg, R. B.↗
NASA NTRS1994

A modified restricted Euler equation for turbulent flows with mean velocity gradients

The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.

Girimaji, Sharath S.↗
NASA NTRS2012

Reynolds-Averaged Turbulence Model Assessment for a Highly Back-Pressured Isolator Flowfield

The use of computational fluid dynamics in scramjet engine component development is widespread in the existing literature. Unfortunately, the quantification of model-form uncertainties is rarely addressed with anything other than sensitivity studies, requiring that the computational results be intimately tied to and calibrated against existing test data. This practice must be replaced with a formal uncertainty quantification process for computational fluid dynamics to play an expanded role in the system design, development, and flight certification process. Due to ground test facility limitations, this expanded role is believed to be a requirement by some in the test and evaluation community if scramjet engines are to be given serious consideration as a viable propulsion device. An effort has been initiated at the NASA Langley Research Center to validate several turbulence closure models used for Reynolds-averaged simulations of scramjet isolator flows. The turbulence models considered were the Menter BSL, Menter SST, Wilcox 1998, Wilcox 2006, and the Gatski-Speziale explicit algebraic Reynolds stress models. The simulations were carried out using the VULCAN computational fluid dynamics package developed at the NASA Langley Research Center. A procedure to quantify the numerical errors was developed to account for discretization errors in the validation process. This procedure utilized the grid convergence index defined by Roache as a bounding estimate for the numerical error. The validation data was collected from a mechanically back-pressured constant area (1 2 inch) isolator model with an isolator entrance Mach number of 2.5. As expected, the model-form uncertainty was substantial for the shock-dominated, massively separated flowfield within the isolator as evidenced by a 6 duct height variation in shock train length depending on the turbulence model employed. Generally speaking, the turbulence models that did not include an explicit stress limiter more closely matched the measured surface pressures. This observation is somewhat surprising, given that stress-limiting models have generally been developed to better predict shock-separated flows. All of the models considered also failed to properly predict the shape and extent of the separated flow region caused by the shock boundary layer interactions. However, the best performing models were able to predict the isolator shock train length (an important metric for isolator operability margin) to within 1 isolator duct height.

Baurle, Robert A.↗
NASA NTRS2003

Characterizing the Severe Turbulence Environments Associated with Commercial Aviation Accidents: Hydrostatic Mesobeta Scale Numerical Simulations of Supergradient Wind Flow and Streamwise Ageostrophic Frontogenesis - Part 2

Simulation experiments reveal key processes that organize a hydrostatic environment conducive to severe turbulence. The paradigm requires juxtaposition of the entrance region of a curved jet stream, which is highly subgeostrophic, with the entrance region of a straight jet stream, which is highly supergeostrophic. The wind and mass fields become misphased as the entrance regions converge resulting in the significant spatial variation of inertial forcing, centripetal forcing, and along- and cross-stream pressure gradient forcing over a mesobeta scale region. This results in frontogenesis and the along-stream divergence of cyclonic and convergence of cyclonic ageostrophic vertical vorticity. The centripetally forced mesoscale front becomes the locus of large gradients of ageostrophic vertical vorticity along an overturning isentrope. This region becomes favorable for streamwise vorticity gradient formation enhancing the environment for organization of horizontal vortex tubes in the presence of buoyant forcing.

Kaplan, Michael L.↗
NASA NTRS1993

Study Of Compressibility Corrections To Turbulence Models

Effects on shear layers in simulated confined and unconfined flows studied. Report presents comparative study of some terms that correct for effects of compressibility in standard k-epsilon mathematical model of turbulence where k denotes turbulence kinetic energy and epsilon denotes rate of dissipation of turbulence kenetic energy. Involved simulation of flows by numerical solution of Reynolds-averaged Navier-Stokes equations.

Viegas, J. R.↗
NASA NTRS1990

Modeling the turbulent kinetic energy equation for compressible, homogeneous turbulence

The turbulent kinetic energy transport equation, which is the basis of turbulence models, is investigated for homogeneous, compressible turbulence using direct numerical simulations performed at CTR. It is shown that the partition between dilatational and solenoidal modes is very sensitive to initial conditions for isotropic decaying turbulence but not for sheared flows. The importance of the dilatational dissipation and of the pressure-dilatation term is evidenced from simulations and a transport equation is proposed to evaluate the pressure-dilatation term evolution. This transport equation seems to work well for sheared flows but does not account for initial condition sensitivity in isotropic decay. An improved model is proposed.

Aupoix, B.↗
NASA NTRS1997

Computational Fluid Dynamics Analysis Method Developed for Rocket-Based Combined Cycle Engine Inlet

Renewed interest in hypersonic propulsion systems has led to research programs investigating combined cycle engines that are designed to operate efficiently across the flight regime. The Rocket-Based Combined Cycle Engine is a propulsion system under development at the NASA Lewis Research Center. This engine integrates a high specific impulse, low thrust-to-weight, airbreathing engine with a low-impulse, high thrust-to-weight rocket. From takeoff to Mach 2.5, the engine operates as an air-augmented rocket. At Mach 2.5, the engine becomes a dual-mode ramjet; and beyond Mach 8, the rocket is turned back on. One Rocket-Based Combined Cycle Engine variation known as the "Strut-Jet" concept is being investigated jointly by NASA Lewis, the U.S. Air Force, Gencorp Aerojet, General Applied Science Labs (GASL), and Lockheed Martin Corporation. Work thus far has included wind tunnel experiments and computational fluid dynamics (CFD) investigations with the NPARC code. The CFD method was initiated by modeling the geometry of the Strut-Jet with the GRIDGEN structured grid generator. Grids representing a subscale inlet model and the full-scale demonstrator geometry were constructed. These grids modeled one-half of the symmetric inlet flow path, including the precompression plate, diverter, center duct, side duct, and combustor. After the grid generation, full Navier-Stokes flow simulations were conducted with the NPARC Navier-Stokes code. The Chien low-Reynolds-number k-e turbulence model was employed to simulate the high-speed turbulent flow. Finally, the CFD solutions were postprocessed with a Fortran code. This code provided wall static pressure distributions, pitot pressure distributions, mass flow rates, and internal drag. These results were compared with experimental data from a subscale inlet test for code validation; then they were used to help evaluate the demonstrator engine net thrust.

Source record↗
NASA NTRS2001

Design and Calibration of Autonomous Coherent Doppler Lidar for Space Missions

Developed a new algorithm for the simulation of three dimensional homogeneous turbulent velocity fields. For typical atmospheric conditions it is impossible to produce a simulated velocity field that simultaneously satisfy a given spatial correlation and the corresponding spatial spectrum because of spectral aliasing. The new algorithms produce a turbulent velocity field which has accurate spatial correlations which is required for performance predictions from space-based systems. Developed a new algorithm for extracting the spatial statistics of the atmospheric velocity field using coherent Doppler lidar. The performance of the algorithm was compared with past methods and the new algorithm produces useful results for space-based data, which was not possible before. Developed new methods for verification of the errors in ground-based and space-based Doppler lidar wind measurements. These new methods do not require independent in situ data. This is an important issue for the verification of space-based Doppler lidar measurements of the global wind field. The performance of the new algorithm was compared with past results for both space-based and ground-based operation. The new algorithm has the best performance and is the only algorithm that performed satisfactory for spacebased operation. The performance of coherent Doppler lidar for a space missions with various scanning geometries was determined using computer simulation which contained the effects of random instrumental velocity errors, wind shear, wind variability along the range-gate and from shot-to-shot, and random variations in atmospheric aerosol backscatter over the measurement volume. The bias in the velocity estimates was small and the accuracy in the is typically less than 0.5 m/s for high signal conditions. For a large number of shot per velocity estimate, the threshold signal level for acceptable estimates is proportional to the number of shots to the minus one half power. This agrees with previous results determined for ground-based measurements. The use of multi-element optical detectors for autonomous operation of coherent Doppler lidar was shown to be a very promising technique. Optimal detector geometries were determined by computer simulation of performance: for ground-based testing with a fixed calibration target and for space-based operation using the random surface returns. The effects of refractive turbulence on ground-based calibration of coherent Doppler lidar was determined by computer simulations and compared with theoretical predictions. New techniques were required to correctly predict performance for the focused beam geometry commonly used for verification of space-based operation. An improved velocity estimator was evaluated for space-based applications were signal shot measurements are used to produce vector wind measurements. This permits more accurate measurements when the signal level is not known a priori or not available from multiple shot measurements. The average Doppler lidar signal spectrum including the effects of velocity turbulence was derived and calculated. This permits new estimation algorithms for turbulence based on spectral estimates. In situ atmospheric measurements were conducted and analyzed using an instrumented kite-platform. This work helps provide the required in situ data for verification of Doppler lidar velocity statistics.

Frehlich, Rod G.↗
NASA NTRS1975

Computation of separated transonic turbulent flows

The two-dimensional Reynolds-averaged compressible Navier-Stokes equations are solved using MacCormack's second-order-accurate explicit finite difference method to simulate the separated transonic turbulent flow field over an airfoil. Four different algebraic eddy viscosity models are tested for viability to achieve turbulence closure for the class of flows considered. These models range from an unmodified boundary-layer mixing-length model to a relaxation model incorporating special considerations for the separation bubble region. Results of this study indicate the necessity for special attention to the separated flow region and suggest limits of applicability of algebraic turbulence models to these separated flow fields.

Deiwert, G. S.↗
NASA NTRS1994

Langmuir turbulence in the auroral ionosphere 2: Nonlinear theory and simulations

A theoretical interpretation of sounding rocket measurements of intense Langmuir wave fields (less than or equal to 500 mV/m) driven by a stream of 20 eV to 4 keV electrons in the lower auroral zone is developed. This interpretation is based on the ability of the 10 microseconds sampling rate of the wave detector to temporally resolve the structure of the Langmuir wave field envelope. A modified form of the Zakharov equations is used to numerically study beam-driven Langmuir turbulence in the presence of a moderate magnetic field (OMEGA (sub e) approximately equals Omega (sub pe). Strong Landau damping on observed nonthermal scattered electrons, which is treated in a companion paper (Newman et al., this issue), plays an important role by inhibiting backscatter cascade and the development of strong turbulence. A parameterized model of the linear electron stream-driven wave instability is introduced, which incorporates limited quasilinear plateau formation. A reasonable set of parameters is found that yields semiquantitative agreement between observed properties of the Langmuir fields and the results of Zakharov equation simulations, including the amplitude and characteristic frequency of the electric field envelope modulations.

Newman, D. L.↗