Search NASA⌕ Search

SEARCH · Search NASA

Results for “eddying model”

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 307 records · Page 17

Modifications of the law of the wall and algebraic turbulence modelling for separated boundary layers

Various modifications of the conventional algebraic eddy viscosity turbulence model are investigated for application to separated flows. Friction velocity is defined in a way that avoids singular behavior at separation and reattachment but reverts to the conventional definition for flows with small pressure gradients. This leads to a modified law of the wall for separated flows. The effect on the calculated flow field of changes in the model that affect the eddy viscosity at various distances from the wall are determined by (1) switching from Prandtl's form to an inner layer formula due to Clauser at various distances from the wall, (2) varying the constant in the Van Driest damping factor, (3) using Clauser's inner layer formula all the way to the wall, and (4) applying a relaxation procedure in the evaluation of the constant in Clauser's inner layer formula. Numerical solutions of the compressible Navier-Stokes equations are used to determine the effects of the modifications. Experimental results from shock-induced separated flows at Mach numbers 2.93 and 8.45 are used for comparison. For these cases improved predictions of wall pressure distribution and positions of separation and reattachment are obtained from the relaxation version of the Clauser inner layer eddy viscosity formula.

Baldwin, B. S.↗

Morphology of the Venus ultraviolet night airglow

Images of the nightside of Venus in the (0,1)delta band of nitric oxide have been obtained by the Pioneer Venus orbiter ultraviolet spectrometer (OUVS). The emission, which is produced by radiative association of N and O, shows a bright spot reaching 5 kR and located at 2 a.m. local solar time just south of the equator. The emitting layer is at 111 + or - 7-km altitude. A one dimensional vertical transport model shows that the hemispheric average brightness of 0.8 kR is consistent with the orbiter neutral mass spectrometer (ONMS) measurements of N and O near 167 km, and that the altitude of the emitting layer is consistent with the eddy mixing model proposed to explain the dayside helium profile measured by the bus neutral mass spectrometer. In the model, N reaches a peak of 7 x 10 to the 8th per cu cm at 114 km, and O reaches a peak of 2.6 x 10 to the 11th per cu cm at 106 km. There is a fair degree of consistency between the ONMS, OUVS, and other airglow measurements, except as regards the local time dependence.

Stewart, A. I. F.↗

A theory for turbulent curved wall jets

An attempt is made to formulate the simplest possible model for the flow of a wall jet emanating from a two dimensional source into quiescent surroundings in the presence of a wall of arbitrary shape. The method used self similar profiles for the mean velocity together with a simple eddy viscosity model. The streamwise and radial momentum equations are integrated across the wall jet flow to give an expression for the momentum balance including the effect of the pressure gradient induced by the rate of change of surface curvature. The streamwise momentum equation is also evaluated at the point of maximum velocity to provide a second equation and thereby permit a solution for the two unknown quantities, the jet half width and the maximum velocity. This approach provides approximate closed form solutions for the flow of the wall jet over surfaces of various shapes, and in particular permits a direct comparison with the available experimental results for plane, cylindrical or logarithmic spiral surfaces.

Roberts, Leonard↗

Wall-modeled LES of the Three-dimensional Speed Bump Experiment

We evaluate the performance of wall-modeled large eddy simulation (WMLES) in predicting smooth-body turbulent flow separation for a three-dimensional Gaussian-shaped speed bump geometry. The Reynolds number based on the bump length is 2 million, and an unstructured compressible finite-volume solver is used with an equilibrium wall model and dynamic subgrid scale model. Spanwise periodic simulations of the centerline two-dimensional bump were used to assess grid resolution requirements, and resolving the thin internal layer in the accelerating region was found to be necessary to correctly capture the downstream separated flow region. Based on these insights, an optimized grid that was smaller by a factor of two provided comparable accuracy to a finer grid that has been used by us and other researchers in past studies. Using insights gained from the spanwise periodic simulation, an unstructured polyhderal grid with about 250 million cells was used for the three-dimensional (3D) configuration with inviscid tunnel side and top walls. Detailed comparisons of wall skin-friction coefficient, wall pressure, velocity and turbulent stresses with available experimental data indicated excellent agreement. While the low-Reynolds Number Spalart-Allmaras RANS model with rotation/curvature correction gave qualitatively good agreement with experiments, WMLES showed significantly more accurate quantitative predictions of the flowfield. The 3D WMLES showed good agreement with experiments in the separated region, and the centerplane results indicated a different separation topology compared to the spanwise periodic simulation.

Computational Fluid Dynamics↗

Relaxation of an unsteady turbulent boundary layer on a flat plate in an expansion tube

An analysis is presented for the relaxation of a turbulent boundary layer on a semi-infinite flat plate after passage of a shock wave and a trailing driver gas-driven gas interface. The problem has special application to expansion-tube flows. The flow-governing equations have been transformed into the Crocco variables, and a time-similar solution is presented in terms of the dimensionless distance-time variable alpha and the dimensionless velocity variable beta. An eddy-viscosity model, similar to that of time-steady boundary layers, is applied to the inner and outer regions of the boundary layer. A turbulent Prandtl number equal to the molecular Prandtl number is used to relate the turbulent heat flux to the eddy viscosity. The numerical results, obtained by using the Gauss-Seidel line-relaxation method, indicate that a fully turbulent boundary layer relaxes faster to the final steady-state values of heat transfer and skin friction than a laminar boundary layer. The results also give a fairly good estimate of the local skin friction and heat transfer for near steady-flow conditions.

Gurta, R. N.↗

Transonic Turbulent Flow Predictions With Two-Equation Turbulence Models

Solutions of the Favre-averaged Navier-Stokes equations for two well-documented transonic turbulent flows are compared in detail with existing experimental data. While the boundary layer in the first case remains attached, a region of extensive flow separation has been observed in the second case. Two recently developed k-epsilon, two-equation, eddy-viscosity models are used to model the turbulence field. These models satisfy the realizability constraints of the Reynolds stresses. Comparisons with the measurements are made for the wall pressure distribution, the mean streamwise velocity profiles, and turbulent quantities. Reasonably good agreement is obtained with the experimental data.

Liou, William W.↗

Sensitivity to Snapshot Frequency in the POD-Based Reduced-Order Modeling of Flow over a Gaussian Bump

The flow over a Gaussian bump at a Reynolds number of 10^6, based on the bump length and freestream conditions, and a Mach number of 0.2, was simulated using wall-modeled large eddy simulation (WMLES). The incoming turbulent flow exhibits relaminarization upstream of the bump apex, due to the favourable pressure gradient acceleration induced by the bump curvature, as well as separation, reattachment and a detached shear layer downstream of the apex. A Galerkin reduced-order model (ROM) was then constructed using a proper orthogonal decomposition (POD) modal basis. The complex, three-dimensional flow features are captured by spatial POD modes differently based on the frequency of the high-fidelity snapshots originally provided. POD modes constructed from differently sampled snapshot matrices were then truncated to retain 90% of the turbulence kinetic energy in the domain of interest. A parametric study is conducted to evaluate the effect of the frequency of snapshots used to construct the POD-ROM.

turbulence↗

Second-Moment RANS Model Verification and Validation Using the Turbulence Modeling Resource Website (Invited)

The implementation of the SSG/LRR-omega differential Reynolds stress model into the NASA flow solvers CFL3D and FUN3D and the DLR flow solver TAU is verified by studying the grid convergence of the solution of three different test cases from the Turbulence Modeling Resource Website. The model's predictive capabilities are assessed based on four basic and four extended validation cases also provided on this website, involving attached and separated boundary layer flows, effects of streamline curvature and secondary flow. Simulation results are compared against experimental data and predictions by the eddy-viscosity models of Spalart-Allmaras (SA) and Menter's Shear Stress Transport (SST).

Eisfeld, Bernhard↗

Plasma discharge in N2 + CH4 at low pressures - Experimental results and applications to Titan

Results are reported from laboratory continuous-flow plasma-discharge experiments designed to simulate the formation of hydrocarbons and nitriles from N2 and CH4 in the atmosphere of Titan. Gas-chromatography and mass-spectrometry data were obtained in experiments lasting up to 100 h at temperature 295 K and pressure 17 or 0.24 mbar, modeling (1) cosmic-ray-induced processes in the Titan troposphere and (2) processes related to stratospheric aurorae excited by energetic electrons and ions from the Saturn magnetosphere, respectively. The results are presented in extensive tables and graphs, and the 0.24-mbar yields are incorporated into an eddy-mixing model to give stratospheric column abundances and mole fractions in good agreement with Voyager IRIS observations.

Thompson, W. Reid↗

Experiments on near-wall structure of three-dimensional boundary layers

Investigations of three-dimensional turbulent boundary layers have shown basic differences between two- and three-dimensional flows. These differences can significantly impact the modeling of three-dimensional flows since many flow models are based on results from two-dimensional boundary layers. In many cases the shear stress vector direction has been shown to lag relative to the direction of the mean velocity gradient as the cross flow grows downstream. Coincidence of these vectors is necessary for a scalar eddy viscosity modeling assumption. A second effect is a reduction in magnitude of the shear stress and/or the shear stress to turbulence energy ratio, a(sub 1). This reduction has been observed in several experiments. Recent numerical simulations also indicate wall-layer structural differences between two- and three-dimensional boundary layers. The differences in structure between two- and three-dimensional boundary layers was also addressed in the experiment of Littell & Eaton. The experiment used two-point correlations to investigate the vortical structures in a three dimensional boundary layer on a spinning disk. It was found that each sign of longitudinal vortex is equally likely to exist, but one sign of vorticity is associated with a structure which is better at producing ejections. The goal of the current investigation is to study the structure of the inner layers. Among other questions, the differences between the effects deduced from the three-dimensional flow simulations and the effects seen in experiments can be examined. The research concentrates on the structure of the wall-layer through flow visualization and direct turbulence measurements down to y(+) = 5.

Flack, Karen A.↗

Numerical investigation of turbulence models for shock separated boundary-layer flows

Numerical solutions of the Navier-Stokes equations for shock separated turbulent boundary-layer flows are presented. Several turbulence models are investigated and assessed by their ability to predict the physical phenomena associated with two extensively documented experiments. The experimental flows consist of shock-wave boundary-layer interactions in axisymmetric internal and external geometries at Mach numbers of 1.5 and 7, respectively. Algebraic and one-equation eddy viscosity models are used to describe the Reynolds shear stress. Calculated values of skin friction, wall pressure distribution, kinetic energy of turbulence, and heat transfer are compared with measurements.

Viegas, J. R.↗

Effects of turbulence mixing, variable properties, and vaporization on spray droplet combustion

Combustion of liquid fuels in the form of spray droplets is simulated numerically. Various vaporization models are examined as to their performance in finite element calculations involving a turbulent flow field. The Eulerian coordinate for the gas and Lagrangian coordinate for the liquid spray droplets are coupled through source terms being updated in the equations of continuity, momentum, and energy. The k-epsilon and modified eddy breakup models are used for simulating turbulent spray combustion flow field. Numerical results for the droplet trajectories, droplet heating, recirculation characteristics, and effects of evaporation models are evaluated. It is also shown that the finite element method is advantageous in dealing with complex geometries, complex boundary conditions, adaptive unstructured grids.

Kim, Y. M.↗

Turbulence modeling for hypersonic flows

Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Marvin, J. G.↗

Turbulence modeling for hypersonic flows

Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Marvin, J. G.↗

Turbulence modeling for hypersonic flows

Turbulence modeling for high-speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models, and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary-layer flows, shock-wave boundary-layer interactions, and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Marvin, J. G.↗

Sensitivity of Middle Atmospheric Temperature and Circulation in the UIUC Mesosphere-Stratosphere-Troposphere GCM to the Treatment of Subgrid-Scale Gravity-Wave Breaking

The sensitivity of the middle atmospheric temperature and circulation to the treatment of mean- flow forcing due to breaking gravity waves was investigated using the University of Illinois at Urbana-Champaign 40-layer Mesosphere-Stratosphere-Troposphere General Circulation Model (MST-GCM). Three GCM experiments were performed. The gravity-wave forcing was represented first by Rayleigh friction, and then by the Alexander and Dunkerton (AD) parameterization with weak and strong breaking effects of gravity waves. In all experiments, the Palmer et al. parameterization was included to treat the breaking of topographic gravity waves in the troposphere and lower stratosphere. Overall, the experiment with the strong breaking effect simulates best the middle atmospheric temperature and circulation. With Rayleigh friction and the weak breaking effect, a large warm bias of up to 60 C was found in the summer upper mesosphere and lower thermosphere. This warm bias was linked to the inability of the GCM to simulate the reversal of the zonal winds from easterly to westerly crossing the mesopause in the summer hemisphere. With the strong breaking effect, the GCM was able to simulate this reversal, and essentially eliminated the warm bias. This improvement was the result of a much stronger meridional transport circulation that possesses a strong vertical ascending branch in the summer upper mesosphere, and hence large adiabatic cooling. Budget analysis indicates that 'in the middle atmosphere the forces that act to maintain a steady zonal-mean zonal wind are primarily those associated with the meridional transport circulation and breaking gravity waves. Contributions from the interaction of the model-resolved eddies with the mean flow are small. To obtain a transport circulation in the mesosphere of the UIUC MST-GCM that is strong enough to produce the observed cold summer mesopause, gravity-wave forcing larger than 100 m/s/day in magnitude is required near the summer mesopause. In the tropics, only with the AD parameterization can the model produce realistic semiannual oscillations.

Yang, Fanglin↗

Second order closure modeling of turbulent buoyant wall plumes

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Zhu, Gang↗

The Role of Gravity Waves in Generating Equatorial Oscillations in Modulating Atmospheric Tides

We discuss a Numerical Spectral Mode (NSM) that extends from the ground up into the thermosphere and incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations in the mean zonal circulation, the semi-annual and quasi-biennial oscillations (SAO and QBO), as well as the tides and planetary waves in the middle atmosphere. Initial results showed that this model can reproduce the salient features observed, including the QBO extending into the upper mesosphere inferred from UARS measurements. The model has now been extended to simulate also: (a) the zonal circulation of the lower stratosphere and upper troposphere, and (b) the upwelling at equatorial latitudes associated with the Brewer Dobsen circulation that affects the dynamics significantly as pointed out by Dunkerton. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase in the model the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. This development is conducive to extending the QBO and SAO to higher latitudes through global scale momentum redistribution. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. In a 3D version of the model, wave momentum is absorbed and dissipated by tides and planetary waves. A somewhat larger GW source (well within the DSP range) is then required to generate realistic QBO and SAO amplitudes. Since GW momentum is deposited in the altitude regime of increasing winds, the amplitude of the diurnal tide is amplified and its vertical wavelength is reduced at altitudes between 70 and 120 km. Wave filtering by the mean zonal circulation causes the GW flux to peak during equinox, and this produces a large semi-annual variation in the tide that has been observed on UARS. Without the diurnal tide, the semidiurnal tide would also be modulated in this way. But the diurnal tide filters out the GW preferentially during equinox, so that the semidiurnal tide tends to peak during solstice. Under the influence of GW, the tides are modulated significantly by planetary waves that are generated preferentially during solstice in part due to baroclinic instability.

Mayr, H. G.↗