Polar symmetric flow of a viscous compressible atmosphere; an application to Mars

NASA Technical Reports Server (NTRS)

Pirraglia, J. A.

1974-01-01

The atmosphere is assumed to be driven by a polar symmetric temperature field and the equations of motion in pressure ratio coordinates are linearized by considering the zero order in terms of a thermal Rossby number R delta I/(2a omega) sq where delta T is a measure of the latitudinal temperature gradient. When the eddy viscosity is greater than 1 million sq cm/sec, the boundary layer extends far up into the atmosphere, making the geostrophic approximation invalid for the bulk of the atmosphere. A temperature model for Mars was used which was based on Mariner 9 infrared spectral data with a 30% increase in the depth averaged temperature from the winter pole to the subsolar point. The results obtained for the increase in surface pressure from the subsolar point to the winter pole, as a function of eddy viscosity and with no-slip conditions imposed at the surface, are given.

Reference Solutions for Benchmark Turbulent Flows in Three Dimensions

NASA Technical Reports Server (NTRS)

Diskin, Boris; Thomas, James L.; Pandya, Mohagna J.; Rumsey, Christopher L.

2016-01-01

A grid convergence study is performed to establish benchmark solutions for turbulent flows in three dimensions (3D) in support of turbulence-model verification campaign at the Turbulence Modeling Resource (TMR) website. The three benchmark cases are subsonic flows around a 3D bump and a hemisphere-cylinder configuration and a supersonic internal flow through a square duct. Reference solutions are computed for Reynolds Averaged Navier Stokes equations with the Spalart-Allmaras turbulence model using a linear eddy-viscosity model for the external flows and a nonlinear eddy-viscosity model based on a quadratic constitutive relation for the internal flow. The study involves three widely-used practical computational fluid dynamics codes developed and supported at NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions computed with these three codes on families of consistently refined grids are presented. Grid-to-grid and code-to-code variations are described in detail.

Reynolds averaged turbulence modelling using deep neural networks with embedded invariance

DOE PAGES

Ling, Julia; Kurzawski, Andrew; Templeton, Jeremy

2016-10-18

There exists significant demand for improved Reynolds-averaged Navier–Stokes (RANS) turbulence models that are informed by and can represent a richer set of turbulence physics. This paper presents a method of using deep neural networks to learn a model for the Reynolds stress anisotropy tensor from high-fidelity simulation data. A novel neural network architecture is proposed which uses a multiplicative layer with an invariant tensor basis to embed Galilean invariance into the predicted anisotropy tensor. It is demonstrated that this neural network architecture provides improved prediction accuracy compared with a generic neural network architecture that does not embed this invariance property.more » Furthermore, the Reynolds stress anisotropy predictions of this invariant neural network are propagated through to the velocity field for two test cases. For both test cases, significant improvement versus baseline RANS linear eddy viscosity and nonlinear eddy viscosity models is demonstrated.« less

Reynolds averaged turbulence modelling using deep neural networks with embedded invariance

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ling, Julia; Kurzawski, Andrew; Templeton, Jeremy

There exists significant demand for improved Reynolds-averaged Navier–Stokes (RANS) turbulence models that are informed by and can represent a richer set of turbulence physics. This paper presents a method of using deep neural networks to learn a model for the Reynolds stress anisotropy tensor from high-fidelity simulation data. A novel neural network architecture is proposed which uses a multiplicative layer with an invariant tensor basis to embed Galilean invariance into the predicted anisotropy tensor. It is demonstrated that this neural network architecture provides improved prediction accuracy compared with a generic neural network architecture that does not embed this invariance property.more » Furthermore, the Reynolds stress anisotropy predictions of this invariant neural network are propagated through to the velocity field for two test cases. For both test cases, significant improvement versus baseline RANS linear eddy viscosity and nonlinear eddy viscosity models is demonstrated.« less

North Atlantic Coast Comprehensive Study (NACCS) Coastal Storm Model Simulations: Waves and Water Levels

DTIC Science & Technology

2015-08-01

published in the NGA’s DNCs, with distinct values assigned to areas of sand, gravel, clay , etc. ERDC/CHL TR-15-14 94 6.5.2 Lateral eddy viscosity As with...6.5.1 Manning’s n bottom friction coefficient ................................................................... 93 6.5.2 Lateral eddy viscosity ...this study include (1) Manning’s n bottom friction coefficient, (2) lateral eddy viscosity , (3) land cover effects on winds (also referred to as

Francis-99 turbine numerical flow simulation of steady state operation using RANS and RANS/LES turbulence model

NASA Astrophysics Data System (ADS)

Minakov, A.; Platonov, D.; Sentyabov, A.; Gavrilov, A.

2017-01-01

We performed numerical simulation of flow in a laboratory model of a Francis hydroturbine at three regimes, using two eddy-viscosity- (EVM) and a Reynolds stress (RSM) RANS models (realizable k-ɛ, k-ω SST, LRR) and detached-eddy-simulations (DES), as well as large-eddy simulations (LES). Comparison of calculation results with the experimental data was carried out. Unlike the linear EVMs, the RSM, DES, and LES reproduced well the mean velocity components, and pressure pulsations in the diffusor draft tube. Despite relatively coarse meshes and insufficient resolution of the near-wall region, LES, DES also reproduced well the intrinsic flow unsteadiness and the dominant flow structures and the associated pressure pulsations in the draft tube.

Eddy Flow during Magma Emplacement: The Basemelt Sill, Antarctica

NASA Astrophysics Data System (ADS)

Petford, N.; Mirhadizadeh, S.

2014-12-01

The McMurdo Dry Valleys magmatic system, Antarctica, forms part of the Ferrar dolerite Large Igneous Province. Comprising a vertical stack of interconnected sills, the complex provides a world-class example of pervasive lateral magma flow on a continental scale. The lowermost intrusion (Basement Sill) offers detailed sections through the now frozen particle macrostructure of a congested magma slurry1. Image-based numerical modelling where the intrusion geometry defines its own unique finite element mesh allows simulations of the flow regime to be made that incorporate realistic magma particle size and flow geometries obtained directly from field measurements. One testable outcome relates to the origin of rhythmic layering where analytical results imply the sheared suspension intersects the phase space for particle Reynolds and Peclet number flow characteristic of macroscopic structures formation2. Another relates to potentially novel crystal-liquid segregation due to the formation of eddies locally at undulating contacts at the floor and roof of the intrusion. The eddies are transient and mechanical in origin, unrelated to well-known fluid dynamical effects around obstacles where flow is turbulent. Numerical particle tracing reveals that these low Re number eddies can both trap (remove) and eject particles back into the magma at a later time according to their mass density. This trapping mechanism has potential to develop local variations in structure (layering) and magma chemistry that may otherwise not occur where the contact between magma and country rock is linear. Simulations indicate that eddy formation is best developed where magma viscosity is in the range 1-102 Pa s. Higher viscosities (> 103 Pa s) tend to dampen the effect implying eddy development is most likely a transient feature. However, it is nice to think that something as simple as a bumpy contact could impart physical and by implication chemical diversity in igneous rocks. 1Marsh, D.B. (2004), A magmatic mush column Rosetta stone: the McMurdo Dry Valleys of Antarcica. EOS, 85, 497-502. 2Petford, N. (2009), Which Effective Viscosity? Mineralogical Magazine, 73, 167-191. Fig. 1. Numerical simulation in the geometry showing magma flow field and eddy formation where circulating magma is trapped. Streamlines track particle orbits.

A priori study of subgrid-scale features in turbulent Rayleigh-Bénard convection

NASA Astrophysics Data System (ADS)

Dabbagh, F.; Trias, F. X.; Gorobets, A.; Oliva, A.

2017-10-01

At the crossroad between flow topology analysis and turbulence modeling, a priori studies are a reliable tool to understand the underlying physics of the subgrid-scale (SGS) motions in turbulent flows. In this paper, properties of the SGS features in the framework of a large-eddy simulation are studied for a turbulent Rayleigh-Bénard convection (RBC). To do so, data from direct numerical simulation (DNS) of a turbulent air-filled RBC in a rectangular cavity of aspect ratio unity and π spanwise open-ended distance are used at two Rayleigh numbers R a ∈{1 08,1 010 } [Dabbagh et al., "On the evolution of flow topology in turbulent Rayleigh-Bénard convection," Phys. Fluids 28, 115105 (2016)]. First, DNS at Ra = 108 is used to assess the performance of eddy-viscosity models such as QR, Wall-Adapting Local Eddy-viscosity (WALE), and the recent S3PQR-models proposed by Trias et al. ["Building proper invariants for eddy-viscosity subgrid-scale models," Phys. Fluids 27, 065103 (2015)]. The outcomes imply that the eddy-viscosity modeling smoothes the coarse-grained viscous straining and retrieves fairly well the effect of the kinetic unfiltered scales in order to reproduce the coherent large scales. However, these models fail to approach the exact evolution of the SGS heat flux and are incapable to reproduce well the further dominant rotational enstrophy pertaining to the buoyant production. Afterwards, the key ingredients of eddy-viscosity, νt, and eddy-diffusivity, κt, are calculated a priori and revealed positive prevalent values to maintain a turbulent wind essentially driven by the mean buoyant force at the sidewalls. The topological analysis suggests that the effective turbulent diffusion paradigm and the hypothesis of a constant turbulent Prandtl number are only applicable in the large-scale strain-dominated areas in the bulk. It is shown that the bulk-dominated rotational structures of vortex-stretching (and its synchronous viscous dissipative structures) hold the highest positive values of νt; however, the zones of backscatter energy and counter-gradient heat transport are related to the areas of compressed focal vorticity. More arguments have been attained through a priori investigation of the alignment trends imposed by existing parameterizations for the SGS heat flux, tested here inside RBC. It is shown that the parameterizations based linearly on the resolved thermal gradient are invalid in RBC. Alternatively, the tensor-diffusivity approach becomes a crucial choice of modeling the SGS heat flux, in particular, the tensorial diffusivity that includes the SGS stress tensor. This and other crucial scrutinies on a future modeling to the SGS heat flux in RBC are sought.

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

NASA Technical Reports Server (NTRS)

Winckelmans, G. S.; Lund, T. S.; Carati, D.; Wray, A. A.

1996-01-01

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

A kinematic eddy viscosity model including the influence of density variations and preturbulence

NASA Technical Reports Server (NTRS)

Cohen, L. S.

1973-01-01

A model for the kinematic eddy viscosity was developed which accounts for the turbulence produced as a result of jet interactions between adjacent streams as well as the turbulence initially present in the streams. In order to describe the turbulence contribution from jet interaction, the eddy viscosity suggested by Prandtl was adopted, and a modification was introduced to account for the effect of density variation through the mixing layer. The form of the modification was ascertained from a study of the compressible turbulent boundary layer on a flat plate. A kinematic eddy viscosity relation which corresponds to the initial turbulence contribution was derived by employing arguments used by Prandtl in his mixing length hypothesis. The resulting expression for self-preserving flow is similar to that which describes the mixing of a submerged jet. Application of the model has led to analytical predictions which are in good agreement with available turbulent mixing experimental data.

Renormalization-group theory for the eddy viscosity in subgrid modeling

NASA Technical Reports Server (NTRS)

Zhou, YE; Vahala, George; Hossain, Murshed

1988-01-01

Renormalization-group theory is applied to incompressible three-dimensional Navier-Stokes turbulence so as to eliminate unresolvable small scales. The renormalized Navier-Stokes equation now includes a triple nonlinearity with the eddy viscosity exhibiting a mild cusp behavior, in qualitative agreement with the test-field model results of Kraichnan. For the cusp behavior to arise, not only is the triple nonlinearity necessary but the effects of pressure must be incorporated in the triple term. The renormalized eddy viscosity will not exhibit a cusp behavior if it is assumed that a spectral gap exists between the large and small scales.

Simulations of Turbulent Flow Over Complex Terrain Using an Immersed-Boundary Method

NASA Astrophysics Data System (ADS)

DeLeon, Rey; Sandusky, Micah; Senocak, Inanc

2018-02-01

We present an immersed-boundary method to simulate high-Reynolds-number turbulent flow over the complex terrain of Askervein and Bolund Hills under neutrally-stratified conditions. We reconstruct both the velocity and the eddy-viscosity fields in the terrain-normal direction to produce turbulent stresses as would be expected from the application of a surface-parametrization scheme based on Monin-Obukhov similarity theory. We find that it is essential to be consistent in the underlying assumptions for the velocity reconstruction and the eddy-viscosity relation to produce good results. To this end, we reconstruct the tangential component of the velocity field using a logarithmic velocity profile and adopt the mixing-length model in the near-surface turbulence model. We use a linear interpolation to reconstruct the normal component of the velocity to enforce the impermeability condition. Our approach works well for both the Askervein and Bolund Hills when the flow is attached to the surface, but shows slight disagreement in regions of flow recirculation, despite capturing the flow reversal.

Simulations of Turbulent Flow Over Complex Terrain Using an Immersed-Boundary Method

NASA Astrophysics Data System (ADS)

DeLeon, Rey; Sandusky, Micah; Senocak, Inanc

2018-06-01

We present an immersed-boundary method to simulate high-Reynolds-number turbulent flow over the complex terrain of Askervein and Bolund Hills under neutrally-stratified conditions. We reconstruct both the velocity and the eddy-viscosity fields in the terrain-normal direction to produce turbulent stresses as would be expected from the application of a surface-parametrization scheme based on Monin-Obukhov similarity theory. We find that it is essential to be consistent in the underlying assumptions for the velocity reconstruction and the eddy-viscosity relation to produce good results. To this end, we reconstruct the tangential component of the velocity field using a logarithmic velocity profile and adopt the mixing-length model in the near-surface turbulence model. We use a linear interpolation to reconstruct the normal component of the velocity to enforce the impermeability condition. Our approach works well for both the Askervein and Bolund Hills when the flow is attached to the surface, but shows slight disagreement in regions of flow recirculation, despite capturing the flow reversal.

Large eddy simulation of spanwise rotating turbulent channel flow with dynamic variants of eddy viscosity model

NASA Astrophysics Data System (ADS)

Jiang, Zhou; Xia, Zhenhua; Shi, Yipeng; Chen, Shiyi

2018-04-01

A fully developed spanwise rotating turbulent channel flow has been numerically investigated utilizing large-eddy simulation. Our focus is to assess the performances of the dynamic variants of eddy viscosity models, including dynamic Vreman's model (DVM), dynamic wall adapting local eddy viscosity (DWALE) model, dynamic σ (Dσ ) model, and the dynamic volumetric strain-stretching (DVSS) model, in this canonical flow. The results with dynamic Smagorinsky model (DSM) and direct numerical simulations (DNS) are used as references. Our results show that the DVM has a wrong asymptotic behavior in the near wall region, while the other three models can correctly predict it. In the high rotation case, the DWALE can get reliable mean velocity profile, but the turbulence intensities in the wall-normal and spanwise directions show clear deviations from DNS data. DVSS exhibits poor predictions on both the mean velocity profile and turbulence intensities. In all three cases, Dσ performs the best.

Large eddy simulation of incompressible turbulent channel flow

NASA Technical Reports Server (NTRS)

Moin, P.; Reynolds, W. C.; Ferziger, J. H.

1978-01-01

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Hybrid LES of Detonations in Reacting Multi-Phase Mixtures

DTIC Science & Technology

2009-02-28

Distortion Theories and Linear Interaction Analyses in order to gain insight in the fundamental processes of compressible turbulence. This analytical...equilibrium. More insight into the development of supersonic mixing layers has been gained later from analytical results, Rapid Distortion Theory ...given by esgs s!t = 0.931i—=•—. Spectral closure theories (Kraichnan [1976]) can be used to evaluate the eddy viscosity formulation as ut = 0.441a -3/2

Numerical Studies into Flow Profiles in Confined Lubricant

NASA Astrophysics Data System (ADS)

di Mare, Luca; Ponjavic, Aleks; Wong, Janet

2013-03-01

This paper documents a computational study of flow profiles in confined fluids. The study is motivated by experimental evidence for deviation from Couette flow found by one of the authors (JSW). The computational study examines several possible stress-strain relations. Since a linear profile is the only possible solution for a constant stress layer even in presence of a power law, the study introduces a functional dependence of the fluid viscosity on the distance from the wall. Based on this dependence, a family of scaling laws for the velocity profile near the wall is derived which matches the measured profiles. The existence of this scaling law requires the viscosity of the fluid to increase at least linearly away from the wall. This behaviour is explained at a microscopic level by considerations on the mobility of long molecules near a wall. This behaviour is reminiscent of the variation of eddy length scales in near-wall turbulence.

Finite-element numerical modeling of atmospheric turbulent boundary layer

NASA Technical Reports Server (NTRS)

Lee, H. N.; Kao, S. K.

1979-01-01

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

On the Subgrid-Scale Modeling of Compressible Turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle; Zeman, Otto

1990-01-01

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

A Particle Representation Model for the Deformation of Homogeneous Turbulence

NASA Technical Reports Server (NTRS)

Kassinos, S. C.; Reynolds, W. C.

1996-01-01

In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.

The effect of artificial bulk viscosity in simulations of forced compressible turbulence

DOE Office of Scientific and Technical Information (OSTI.GOV)

Campos, A.; Morgan, B.

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performancemore » of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.« less

The effect of artificial bulk viscosity in simulations of forced compressible turbulence

DOE PAGES

Campos, A.; Morgan, B.

2018-05-17

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performancemore » of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.« less

Quasi-laminar stability and sensitivity analyses for turbulent flows: Prediction of low-frequency unsteadiness and passive control

NASA Astrophysics Data System (ADS)

Mettot, Clément; Sipp, Denis; Bézard, Hervé

2014-04-01

This article presents a quasi-laminar stability approach to identify in high-Reynolds number flows the dominant low-frequencies and to design passive control means to shift these frequencies. The approach is based on a global linear stability analysis of mean-flows, which correspond to the time-average of the unsteady flows. Contrary to the previous work by Meliga et al. ["Sensitivity of 2-D turbulent flow past a D-shaped cylinder using global stability," Phys. Fluids 24, 061701 (2012)], we use the linearized Navier-Stokes equations based solely on the molecular viscosity (leaving aside any turbulence model and any eddy viscosity) to extract the least stable direct and adjoint global modes of the flow. Then, we compute the frequency sensitivity maps of these modes, so as to predict before hand where a small control cylinder optimally shifts the frequency of the flow. In the case of the D-shaped cylinder studied by Parezanović and Cadot [J. Fluid Mech. 693, 115 (2012)], we show that the present approach well captures the frequency of the flow and recovers accurately the frequency control maps obtained experimentally. The results are close to those already obtained by Meliga et al., who used a more complex approach in which turbulence models played a central role. The present approach is simpler and may be applied to a broader range of flows since it is tractable as soon as mean-flows — which can be obtained either numerically from simulations (Direct Numerical Simulation (DNS), Large Eddy Simulation (LES), unsteady Reynolds-Averaged-Navier-Stokes (RANS), steady RANS) or from experimental measurements (Particle Image Velocimetry - PIV) — are available. We also discuss how the influence of the control cylinder on the mean-flow may be more accurately predicted by determining an eddy-viscosity from numerical simulations or experimental measurements. From a technical point of view, we finally show how an existing compressible numerical simulation code may be used in a black-box manner to extract the global modes and sensitivity maps.

Hybrid Reynolds-Averaged/Large Eddy Simulation of the Flow in a Model SCRamjet Cavity Flameholder

NASA Technical Reports Server (NTRS)

Baurle, R. A.

2016-01-01

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. Experimental data available for this configuration include velocity statistics obtained from particle image velocimetry. Several turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged/large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This e ort was undertaken to not only assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community, but to also begin to understand how this capability can best be used to augment standard Reynolds-averaged simulations. The numerical errors were quantified for the steady-state simulations, and at least qualitatively assessed for the scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results displayed a high degree of variability when comparing the flameholder fuel distributions obtained from each turbulence model. This prompted the consideration of applying the higher-fidelity scale-resolving simulations as a surrogate "truth" model to calibrate the Reynolds-averaged closures in a non-reacting setting prior to their use for the combusting simulations. In general, the Reynolds-averaged velocity profile predictions at the lowest fueling level matched the particle imaging measurements almost as well as was observed for the non-reacting condition. However, the velocity field predictions proved to be more sensitive to the flameholder fueling rate than was indicated in the measurements.

Large Eddy Simulation ... Where Do We Stand? International Workshop Held in St. Petersburg Beach, Florida on 19-21 December 1990.

DTIC Science & Technology

1990-01-01

S. Orszag, Chairman 1. P. Moin Some Issues in Computation of Turbulent Flows. 2. M. Lesieur, P. Comte, X. Normand, 0. Metais and A. Silveira Spectral...Richtmeyer’s computational experience with one-dimensional shock waves (1950) indicated the value of a non-linear artificial viscosity. Charney and... computer architecture and the advantages of semi-Lagrangian advective schemes may lure large-scale atmospheric modelers back to finite-difference

Vertical structure of mean cross-shore currents across a barred surf zone

USGS Publications Warehouse

Haines, John W.; Sallenger, Asbury H.

1994-01-01

Mean cross-shore currents observed across a barred surf zone are compared to model predictions. The model is based on a simplified momentum balance with a turbulent boundary layer at the bed. Turbulent exchange is parameterized by an eddy viscosity formulation, with the eddy viscosity Aυ independent of time and the vertical coordinate. Mean currents result from gradients due to wave breaking and shoaling, and the presence of a mean setup of the free surface. Descriptions of the wave field are provided by the wave transformation model of Thornton and Guza [1983]. The wave transformation model adequately reproduces the observed wave heights across the surf zone. The mean current model successfully reproduces the observed cross-shore flows. Both observations and predictions show predominantly offshore flow with onshore flow restricted to a relatively thin surface layer. Successful application of the mean flow model requires an eddy viscosity which varies horizontally across the surf zone. Attempts are made to parameterize this variation with some success. The data does not discriminate between alternative parameterizations proposed. The overall variability in eddy viscosity suggested by the model fitting should be resolvable by field measurements of the turbulent stresses. Consistent shortcomings of the parameterizations, and the overall modeling effort, suggest avenues for further development and data collection.

Thermo-fluid-dynamics of turbulent boundary layer over a moving continuous flat sheet in a parallel free stream

NASA Astrophysics Data System (ADS)

Afzal, Bushra; Noor Afzal Team; Bushra Afzal Team

2014-11-01

The momentum and thermal turbulent boundary layers over a continuous moving sheet subjected to a free stream have been analyzed in two layers (inner wall and outer wake) theory at large Reynolds number. The present work is based on open Reynolds equations of momentum and heat transfer without any closure model say, like eddy viscosity or mixing length etc. The matching of inner and outer layers has been carried out by Izakson-Millikan-Kolmogorov hypothesis. The matching for velocity and temperature profiles yields the logarithmic laws and power laws in overlap region of inner and outer layers, along with friction factor and heat transfer laws. The uniformly valid solution for velocity, Reynolds shear stress, temperature and thermal Reynolds heat flux have been proposed by introducing the outer wake functions due to momentum and thermal boundary layers. The comparison with experimental data for velocity profile, temperature profile, skin friction and heat transfer are presented. In outer non-linear layers, the lowest order momentum and thermal boundary layer equations have also been analyses by using eddy viscosity closure model, and results are compared with experimental data. Retired Professor, Embassy Hotel, Rasal Ganj, Aligarh 202001 India.

Effects of eddy viscosity and thermal conduction and Coriolis force in the dynamics of gravity wave driven fluctuations in the OH nightglow

NASA Technical Reports Server (NTRS)

Hickey, M. P.

1988-01-01

The chemical-dynamical model of Walterscheid et al. (1987), which describes wave-driven fluctuations in OH nightglow, was modified to include the effects of both eddy thermal conduction and viscosity, as well as the Coriolis force (with the shallow atmosphere approximation). Using the new model, calculations were performed for the same nominal case as used by Walterscheid et al. but with only wave periods considered. For this case, the Coriolis force was found to be unimportant at any wave period. For wave periods greater than 2 or 3 hours, the inclusion of thermal conduction alone greatly modified the results (in terms of a complex ratio 'eta' which expresses the relationship between the intensity oscillation about the time-averaged intensity and the temperature oscillation about the time-averaged temperature); this effect was reduced with the further inclusion of the eddy viscosity.

Numerical dissipation vs. subgrid-scale modelling for large eddy simulation

NASA Astrophysics Data System (ADS)

Dairay, Thibault; Lamballais, Eric; Laizet, Sylvain; Vassilicos, John Christos

2017-05-01

This study presents an alternative way to perform large eddy simulation based on a targeted numerical dissipation introduced by the discretization of the viscous term. It is shown that this regularisation technique is equivalent to the use of spectral vanishing viscosity. The flexibility of the method ensures high-order accuracy while controlling the level and spectral features of this purely numerical viscosity. A Pao-like spectral closure based on physical arguments is used to scale this numerical viscosity a priori. It is shown that this way of approaching large eddy simulation is more efficient and accurate than the use of the very popular Smagorinsky model in standard as well as in dynamic version. The main strength of being able to correctly calibrate numerical dissipation is the possibility to regularise the solution at the mesh scale. Thanks to this property, it is shown that the solution can be seen as numerically converged. Conversely, the two versions of the Smagorinsky model are found unable to ensure regularisation while showing a strong sensitivity to numerical errors. The originality of the present approach is that it can be viewed as implicit large eddy simulation, in the sense that the numerical error is the source of artificial dissipation, but also as explicit subgrid-scale modelling, because of the equivalence with spectral viscosity prescribed on a physical basis.

Computation of Turbulent Wake Flows in Variable Pressure Gradient

NASA Technical Reports Server (NTRS)

Duquesne, N.; Carlson, J. R.; Rumsey, C. L.; Gatski, T. B.

1999-01-01

Transport aircraft performance is strongly influenced by the effectiveness of high-lift systems. Developing wakes generated by the airfoil elements are subjected to strong pressure gradients and can thicken very rapidly, limiting maximum lift. This paper focuses on the effects of various pressure gradients on developing symmetric wakes and on the ability of a linear eddy viscosity model and a non-linear explicit algebraic stress model to accurately predict their downstream evolution. In order to reduce the uncertainties arising from numerical issues when assessing the performance of turbulence models, three different numerical codes with the same turbulence models are used. Results are compared to available experimental data to assess the accuracy of the computational results.

A Study of Grid Resolution, Transition and Turbulence Model Using the Transonic Simple Straked Delta Wing

NASA Technical Reports Server (NTRS)

Bartels, Robert E.

2001-01-01

Three-dimensional transonic flow over a delta wing is investigated using several turbulence models. The performance of linear eddy viscosity models and an explicit algebraic stress model is assessed at the start of vortex flow, and the results compared with experimental data. To assess the effect of transition location, computations that either fix transition aft of the leading edge or are fully turbulent are performed. These computations show that grid resolution, transition location and turbulence model significantly affect the 3D flowfield.

Local dynamic subgrid-scale models in channel flow

NASA Technical Reports Server (NTRS)

Cabot, William H.

1994-01-01

The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 1: Analysis description

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. The governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models are described in detail.

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 3D has been developed to solve the three dimensional, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort has been to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation have been emphasized. The governing equations are solved in generalized non-orthogonal body-fitted coordinates by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the Analysis Description, and presents the equations and solution procedure. It describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Convection without eddy viscosity: An attempt to model the interiors of giant planets

NASA Technical Reports Server (NTRS)

Ingersoll, A. P.

1986-01-01

In the theory of hydrostatic quasi-geostrophic flow in the Earth's atmosphere the principal results do not depend on the eddy viscosity. This contrasts with published theories of convection in deep rotating fluid spheres, where the wavelength of the fastest growing disturbance varies as E sup 1/3, where E, the Ekman number, is proportional to the eddy viscosity. A new theory of quasi-columnar motions in stably stratified fluid spheres attempts to capture the luck of the meteorologists. The theory allows one to investigate the stability of barotropic and baroclinic zonal flows that extend into the planetary interior. It is hypothesized that the internal heat Jupiter and Saturn comes out not radially but on sloping surfaces defined by the internal entropy distribution. To test the hypothesis one searches for basic states in which the wavelength of the fastest-growing disturbance remains finite as E tends to zero, and is which the heat flux vector is radially outward and poleward.

Influence of a thin compressible insoluble liquid film on the eddy currents generated by interacting surface waves

NASA Astrophysics Data System (ADS)

Parfenyev, Vladimir M.; Vergeles, Sergey S.

2018-06-01

Recently the generation of eddy currents by interacting surface waves was observed experimentally. The phenomenon provides the possibility for manipulation of particles which are immersed in the fluid. The analysis shows that the amplitude of the established eddy currents produced by stationary surface waves does not depend on the fluid viscosity in the free surface case. The currents become parametrically larger, being inversely proportional to the square root of the fluid viscosity in the case when the fluid surface is covered by an almost incompressible thin liquid (i.e., shear elasticity is zero) film formed by an insoluble agent with negligible internal viscous losses as compared to the dissipation in the fluid bulk. Here we extend the theory for a thin insoluble film with zero shear elasticity and small shear and dilational viscosities on the case of an arbitrary elastic compression modulus. We find both contributions into the Lagrangian motion of passive tracers, which are the advection by the Eulerian vertical vorticity and the Stokes drift. Whereas the Stokes drift contribution preserves its value for the free surface case outside a thin viscous sublayer, the Eulerian vertical vorticity strongly depends on the fluid viscosity at high values of the film compression modulus. The Stokes drift acquires a strong dependence on the fluid viscosity inside the viscous sublayer; however, the change is compensated by an opposite change in the Eulerian vertical vorticity. As a result, the vertical dependence of the intensity of eddy currents is given by a sum of two decaying exponents with both decrements being of the order of the wave number. The decrements are numerically different, so the Eulerian contribution becomes dominant at some depth for the surface film with any compression modulus.

An analytical model of capped turbulent oscillatory bottom boundary layers

NASA Astrophysics Data System (ADS)

Shimizu, Kenji

2010-03-01

An analytical model of capped turbulent oscillatory bottom boundary layers (BBLs) is proposed using eddy viscosity of a quadratic form. The common definition of friction velocity based on maximum bottom shear stress is found unsatisfactory for BBLs under rotating flows, and a possible extension based on turbulent kinetic energy balance is proposed. The model solutions show that the flow may slip at the top of the boundary layer due to capping by the water surface or stratification, reducing the bottom shear stress, and that the Earth's rotation induces current and bottom shear stress components perpendicular to the interior flow with a phase lag (or lead). Comparisons with field and numerical experiments indicate that the model predicts the essential characteristics of the velocity profiles, although the agreement is rather qualitative due to assumptions of quadratic eddy viscosity with time-independent friction velocity and a well-mixed boundary layer. On the other hand, the predicted linear friction coefficients, phase lead, and veering angle at the bottom agreed with available data with an error of 3%-10%, 5°-10°, and 5°-10°, respectively. As an application of the model, the friction coefficients are used to calculate e-folding decay distances of progressive internal waves with a semidiurnal frequency.

Water hammer prediction and control: the Green's function method

NASA Astrophysics Data System (ADS)

Xuan, Li-Jun; Mao, Feng; Wu, Jie-Zhi

2012-04-01

By Green's function method we show that the water hammer (WH) can be analytically predicted for both laminar and turbulent flows (for the latter, with an eddy viscosity depending solely on the space coordinates), and thus its hazardous effect can be rationally controlled and minimized. To this end, we generalize a laminar water hammer equation of Wang et al. (J. Hydrodynamics, B2, 51, 1995) to include arbitrary initial condition and variable viscosity, and obtain its solution by Green's function method. The predicted characteristic WH behaviors by the solutions are in excellent agreement with both direct numerical simulation of the original governing equations and, by adjusting the eddy viscosity coefficient, experimentally measured turbulent flow data. Optimal WH control principle is thereby constructed and demonstrated.

Moffatt eddies at an interface

NASA Astrophysics Data System (ADS)

Shtern, Vladimir

2014-12-01

It is shown that an infinite set of eddies can develop near the interface-wall intersection in a two-fluid flow. A striking feature is that the eddy occurrence depends on from what side of the interface the flow is driven. In air-water flows where the viscosity ratio is 0.018, the eddies develop if a driving source is located on (i) the air side for , (ii) any side for , and (iii) the water side for , where is the upper interface-wall angle.

Evaluation of the Momentum Closure Schemes in MPAS-Ocean

NASA Astrophysics Data System (ADS)

Zhao, Shimei; Liu, Yudi; Liu, Wei

2018-04-01

In order to compare and evaluate the performances of the Laplacian viscosity closure, the biharmonic viscosity closure, and the Leith closure momentum schemes in the MPAS-Ocean model, a variety of physical quantities, such as the relative reference potential energy (RPE) change, the RPE time change rate (RPETCR), the grid Reynolds number, the root mean square (RMS) of kinetic energy, and the spectra of kinetic energy and enstrophy, are calculated on the basis of results of a 3D baroclinic periodic channel. Results indicate that: 1) The RPETCR demonstrates a saturation phenomenon in baroclinic eddy tests. The critical grid Reynolds number corresponding to RPETCR saturation differs between the three closures: the largest value is in the biharmonic viscosity closure, followed by that in the Laplacian viscosity closure, and that in the Leith closure is the smallest. 2) All three closures can effectively suppress spurious dianeutral mixing by reducing the grid Reynolds number under sub-saturation conditions of the RPETCR, but they can also damage certain physical processes. Generally, the damage to the rotation process is greater than that to the advection process. 3) The dissipation in the biharmonic viscosity closure is strongly dependent on scales. Most dissipation concentrates on small scales, and the energy of small-scale eddies is often transferred to large-scale kinetic energy. The viscous dissipation in the Laplacian viscosity closure is the strongest on various scales, followed by that in the Leith closure. Note that part of the small-scale kinetic energy is also transferred to large-scale kinetic energy in the Leith closure. 4) The characteristic length scale L and the dimensionless parameter D in the Leith closure are inherently coupled. The RPETCR is inversely proportional to the product of D and L. When the product of D and L is constant, both the simulated RPETCR and the inhibition of spurious dianeutral mixing are the same in all tests using the Leith closure. The dissipative scale in the Leith closure depends on the parameter L, and the dissipative intensity depends on the parameter D. 5) Although optimal results may not be achieved by using the optimal parameters obtained from the 2D barotropic model in the 3D baroclinic simulation, the total energies are dissipative in all three closures. Dissipation is the strongest in the biharmonic viscosity closure, followed by that in the Leith closure, and that in the Laplacian viscosity closure is the weakest. Mesoscale eddies develop the fastest in the biharmonic viscosity closure after the baroclinic adjustment process finishes, and the kinetic energy reaches its maximum, which is attributed to the smallest dissipation of enstrophy in the biharmonic viscosity closure. Mesoscale eddies develop the slowest, and the kinetic energy peak value is the smallest in the Laplacian viscosity closure. Results in the Leith closure are between that in the biharmonic viscosity closure and the Laplacian viscosity closure.

Existence of k⁻¹ power-law scaling in the equilibrium regions of wall-bounded turbulence explained by Heisenberg's eddy viscosity.

PubMed

Katul, Gabriel G; Porporato, Amilcare; Nikora, Vladimir

2012-12-01

The existence of a "-1" power-law scaling at low wavenumbers in the longitudinal velocity spectrum of wall-bounded turbulence was explained by multiple mechanisms; however, experimental support has not been uniform across laboratory studies. This letter shows that Heisenberg's eddy viscosity approach can provide a theoretical framework that bridges these multiple mechanisms and explains the elusiveness of the "-1" power law in some experiments. Novel theoretical outcomes are conjectured about the role of intermittency and very-large scale motions in modifying the k⁻¹ scaling.

A Nonlinear Interactions Approximation Model for Large-Eddy Simulation

NASA Astrophysics Data System (ADS)

Haliloglu, Mehmet U.; Akhavan, Rayhaneh

2003-11-01

A new approach to LES modelling is proposed based on direct approximation of the nonlinear terms \\overlineu_iuj in the filtered Navier-Stokes equations, instead of the subgrid-scale stress, τ_ij. The proposed model, which we call the Nonlinear Interactions Approximation (NIA) model, uses graded filters and deconvolution to parameterize the local interactions across the LES cutoff, and a Smagorinsky eddy viscosity term to parameterize the distant interactions. A dynamic procedure is used to determine the unknown eddy viscosity coefficient, rendering the model free of adjustable parameters. The proposed NIA model has been applied to LES of turbulent channel flows at Re_τ ≈ 210 and Re_τ ≈ 570. The results show good agreement with DNS not only for the mean and resolved second-order turbulence statistics but also for the full (resolved plus subgrid) Reynolds stress and turbulence intensities.

Tests of dynamic Lagrangian eddy viscosity models in Large Eddy Simulations of flow over three-dimensional bluff bodies

NASA Astrophysics Data System (ADS)

Tseng, Yu-Heng; Meneveau, Charles; Parlange, Marc B.

2004-11-01

Large Eddy Simulations (LES) of atmospheric boundary-layer air movement in urban environments are especially challenging due to complex ground topography. Typically in such applications, fairly coarse grids must be used where the subgrid-scale (SGS) model is expected to play a crucial role. A LES code using pseudo-spectral discretization in horizontal planes and second-order differencing in the vertical is implemented in conjunction with the immersed boundary method to incorporate complex ground topography, with the classic equilibrium log-law boundary condition in the new-wall region, and with several versions of the eddy-viscosity model: (1) the constant-coefficient Smagorinsky model, (2) the dynamic, scale-invariant Lagrangian model, and (3) the dynamic, scale-dependent Lagrangian model. Other planar-averaged type dynamic models are not suitable because spatial averaging is not possible without directions of statistical homogeneity. These SGS models are tested in LES of flow around a square cylinder and of flow over surface-mounted cubes. Effects on the mean flow are documented and found not to be major. Dynamic Lagrangian models give a physically more realistic SGS viscosity field, and in general, the scale-dependent Lagrangian model produces larger Smagorinsky coefficient than the scale-invariant one, leading to reduced distributions of resolved rms velocities especially in the boundary layers near the bluff bodies.

Three-dimensional Navier-Stokes analysis of turbine passage heat transfer

NASA Technical Reports Server (NTRS)

Ameri, Ali A.; Arnone, Andrea

1991-01-01

The three-dimensional Reynolds-averaged Navier-Stokes equations are numerically solved to obtain the pressure distribution and heat transfer rates on the endwalls and the blades of two linear turbine cascades. The TRAF3D code which has recently been developed in a joint project between researchers from the University of Florence and NASA Lewis Research Center is used. The effect of turbulence is taken into account by using the eddy viscosity hypothesis and the two-layer mixing length model of Baldwin and Lomax. Predictions of surface heat transfer are made for Langston's cascade and compared with the data obtained for that cascade by Graziani. The comparison was found to be favorable. The code is also applied to a linear transonic rotor cascade to predict the pressure distributions and heat transfer rates.

A sensitivity analysis of low salinity habitats simulated by a hydrodynamic model in the Manatee River estuary in Florida, USA

NASA Astrophysics Data System (ADS)

Chen, XinJian

2012-06-01

This paper presents a sensitivity study of simulated availability of low salinity habitats by a hydrodynamic model for the Manatee River estuary located in the southwest portion of the Florida peninsula. The purpose of the modeling study was to establish a regulatory minimum freshwater flow rate required to prevent the estuarine ecosystem from significant harm. The model used in the study was a multi-block model that dynamically couples a three-dimensional (3D) hydrodynamic model with a laterally averaged (2DV) hydrodynamic model. The model was calibrated and verified against measured real-time data of surface elevation and salinity at five stations during March 2005-July 2006. The calibrated model was then used to conduct a series of scenario runs to investigate effects of the flow reduction on salinity distributions in the Manatee River estuary. Based on simulated salinity distribution in the estuary, water volumes, bottom areas and shoreline lengths for salinity less than certain predefined values were calculated and analyzed to help establish the minimum freshwater flow rate for the estuarine system. The sensitivity analysis conducted during the modeling study for the Manatee River estuary examined effects of the bottom roughness, ambient vertical eddy viscosity/diffusivity, horizontal eddy viscosity/diffusivity, and ungauged flow on the model results and identified the relative importance of these model parameters (input data) to the outcome of the availability of low salinity habitats. It is found that the ambient vertical eddy viscosity/diffusivity is the most influential factor controlling the model outcome, while the horizontal eddy viscosity/diffusivity is the least influential one.

A Realizable Reynolds Stress Algebraic Equation Model

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Zhu, Jiang; Lumley, John L.

1993-01-01

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Eddy damping, backscatter, and subgrid stresses in subgrid modeling of turbulence

NASA Technical Reports Server (NTRS)

Zhou, YE

1991-01-01

The Brief Report demonstrates the relationship of eddy-viscosity models to subgrid stresses. A formula that determines the relative importance of the cross stress and the Reynolds stress for the net eddy-damping and backscatter contributions is derived. The cross-stress term with sharp-cut filtering is identified as an important quantity to model. These concepts could prove useful as a basis for constructing specific models for the Reynolds and cross stresses.

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Benson, Thomas J.; Suresh, Ambady

1990-01-01

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Some effects of swirl on turbulent mixing and combustion

NASA Technical Reports Server (NTRS)

Rubel, A.

1972-01-01

A general formulation of some effects of swirl on turbulent mixing is given. The basis for the analysis is that momentum transport is enhanced by turbulence resulting from rotational instability of the fluid field. An appropriate form for the turbulent eddy viscosity is obtained by mixing length type arguments. The result takes the form of a corrective factor that is a function of the swirl and acts to increase the eddy viscosity. The factor is based upon the initial mixing conditions implying that the rotational turbulence decays in a manner similar to that of free shear turbulence. Existing experimental data for free jet combustion are adequately matched by using the modifying factor to relate the effects of swirl on eddy viscosity. The model is extended and applied to the supersonic combustion of a ring jet of hydrogen injected into a constant area annular air stream. The computations demonstrate that swirling the flow could: (1) reduce the burning length by one half, (2) result in more uniform burning across the annulus width, and (3) open the possibility of optimization of the combustion characteristics by locating the fuel jet between the inner wall and center of the annulus width.

Application of renormalization group theory to the large-eddy simulation of transitional boundary layers

NASA Technical Reports Server (NTRS)

Piomelli, Ugo; Zang, Thomas A.; Speziale, Charles G.; Lund, Thomas S.

1990-01-01

An eddy viscosity model based on the renormalization group theory of Yakhot and Orszag (1986) is applied to the large-eddy simulation of transition in a flat-plate boundary layer. The simulation predicts with satisfactory accuracy the mean velocity and Reynolds stress profiles, as well as the development of the important scales of motion. The evolution of the structures characteristic of the nonlinear stages of transition is also predicted reasonably well.

Improving a two-equation eddy-viscosity turbulence model to predict the aerodynamic performance of thick wind turbine airfoils

NASA Astrophysics Data System (ADS)

Bangga, Galih; Kusumadewi, Tri; Hutomo, Go; Sabila, Ahmad; Syawitri, Taurista; Setiadi, Herlambang; Faisal, Muhamad; Wiranegara, Raditya; Hendranata, Yongki; Lastomo, Dwi; Putra, Louis; Kristiadi, Stefanus

2018-03-01

Numerical simulations for relatively thick airfoils are carried out in the present studies. An attempt to improve the accuracy of the numerical predictions is done by adjusting the turbulent viscosity of the eddy-viscosity Menter Shear-Stress-Transport (SST) model. The modification involves the addition of a damping factor on the wall-bounded flows incorporating the ratio of the turbulent kinetic energy to its specific dissipation rate for separation detection. The results are compared with available experimental data and CFD simulations using the original Menter SST model. The present model improves the lift polar prediction even though the stall angle is still overestimated. The improvement is caused by the better prediction of separated flow under a strong adverse pressure gradient. The results show that the Reynolds stresses are damped near the wall causing variation of the logarithmic velocity profiles.

LARGE EDDY SIMULATION OF TURBULENT FLOW OVER MARGINALLY RESOLVED THREE BLUFF BODIES USING AN IMMERSED BOUNDARY METHOD AND LAGRANGIAN DYNAMIC EDDY-VISCOSITY MODELS. (R828771C004)

EPA Science Inventory

The perspectives, information and conclusions conveyed in research project abstracts, progress reports, final reports, journal abstracts and journal publications convey the viewpoints of the principal investigator and may not represent the views and policies of ORD and EPA. Concl...

Atmospheric flow over two-dimensional bluff surface obstructions

NASA Technical Reports Server (NTRS)

Bitte, J.; Frost, W.

1976-01-01

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

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

NASA Technical Reports Server (NTRS)

Baldwin, B. S.; Maccormack, R. W.

1976-01-01

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.

Multigrid calculation of three-dimensional viscous cascade flows

NASA Technical Reports Server (NTRS)

Arnone, A.; Liou, M.-S.; Povinelli, L. A.

1991-01-01

A 3-D code for viscous cascade flow prediction was developed. The space discretization uses a cell-centered scheme with eigenvalue scaling to weigh the artificial dissipation terms. Computational efficiency of a four stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full multigrid method. The Baldwin-Lomax eddy viscosity model is used for turbulence closure. A zonal, nonperiodic grid is used to minimize mesh distortion in and downstream of the throat region. Applications are presented for an annular vane with and without end wall contouring, and for a large scale linear cascade. The calculation is validated by comparing with experiments and by studying grid dependency.

Multigrid calculation of three-dimensional viscous cascade flows

NASA Technical Reports Server (NTRS)

Arnone, A.; Liou, M.-S.; Povinelli, L. A.

1991-01-01

A three-dimensional code for viscous cascade flow prediction has been developed. The space discretization uses a cell-centered scheme with eigenvalue scaling to weigh the artificial dissipation terms. Computational efficiency of a four-stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full-multigrid method. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. A zonal, nonperiodic grid is used to minimize mesh distortion in and downstream of the throat region. Applications are presented for an annular vane with and without end wall contouring, and for a large-scale linear cascade. The calculation is validated by comparing with experiments and by studying grid dependency.

Matrix exponential-based closures for the turbulent subgrid-scale stress tensor.

PubMed

Li, Yi; Chevillard, Laurent; Eyink, Gregory; Meneveau, Charles

2009-01-01

Two approaches for closing the turbulence subgrid-scale stress tensor in terms of matrix exponentials are introduced and compared. The first approach is based on a formal solution of the stress transport equation in which the production terms can be integrated exactly in terms of matrix exponentials. This formal solution of the subgrid-scale stress transport equation is shown to be useful to explore special cases, such as the response to constant velocity gradient, but neglecting pressure-strain correlations and diffusion effects. The second approach is based on an Eulerian-Lagrangian change of variables, combined with the assumption of isotropy for the conditionally averaged Lagrangian velocity gradient tensor and with the recent fluid deformation approximation. It is shown that both approaches lead to the same basic closure in which the stress tensor is expressed as the matrix exponential of the resolved velocity gradient tensor multiplied by its transpose. Short-time expansions of the matrix exponentials are shown to provide an eddy-viscosity term and particular quadratic terms, and thus allow a reinterpretation of traditional eddy-viscosity and nonlinear stress closures. The basic feasibility of the matrix-exponential closure is illustrated by implementing it successfully in large eddy simulation of forced isotropic turbulence. The matrix-exponential closure employs the drastic approximation of entirely omitting the pressure-strain correlation and other nonlinear scrambling terms. But unlike eddy-viscosity closures, the matrix exponential approach provides a simple and local closure that can be derived directly from the stress transport equation with the production term, and using physically motivated assumptions about Lagrangian decorrelation and upstream isotropy.

Numerical modeling of turbulent swirling flow in a multi-inlet vortex nanoprecipitation reactor using dynamic DDES

NASA Astrophysics Data System (ADS)

Hill, James C.; Liu, Zhenping; Fox, Rodney O.; Passalacqua, Alberto; Olsen, Michael G.

2015-11-01

The multi-inlet vortex reactor (MIVR) has been developed to provide a platform for rapid mixing in the application of flash nanoprecipitation (FNP) for manufacturing functional nanoparticles. Unfortunately, commonly used RANS methods are unable to accurately model this complex swirling flow. Large eddy simulations have also been problematic, as expensive fine grids to accurately model the flow are required. These dilemmas led to the strategy of applying a Delayed Detached Eddy Simulation (DDES) method to the vortex reactor. In the current work, the turbulent swirling flow inside a scaled-up MIVR has been investigated by using a dynamic DDES model. In the DDES model, the eddy viscosity has a form similar to the Smagorinsky sub-grid viscosity in LES and allows the implementation of a dynamic procedure to determine its coefficient. The complex recirculating back flow near the reactor center has been successfully captured by using this dynamic DDES model. Moreover, the simulation results are found to agree with experimental data for mean velocity and Reynolds stresses.

The tidally averaged momentum balance in a partially and periodically stratified estuary

USGS Publications Warehouse

Stacey, M.T.; Brennan, Matthew L.; Burau, J.R.; Monismith, Stephen G.

2010-01-01

Observations of turbulent stresses and mean velocities over an entire spring-neap cycle are used to evaluate the dynamics of tidally averaged flows in a partially stratified estuarine channel. In a depth-averaged sense, the net flow in this channel is up estuary due to interaction of tidal forcing with the geometry of the larger basin. The depth-variable tidally averaged flow has the form of an estuarine exchange flow (downstream at the surface, upstream at depth) and varies in response to the neap-spring transition. The weakening of the tidally averaged exchange during the spring tides appears to be a result of decreased stratification on the tidal time scale rather than changes in bed stress. The dynamics of the estuarine exchange flow are defined by a balance between the vertical divergence of the tidally averaged turbulent stress and the tidally averaged pressure gradient in the lower water column. In the upper water column, tidal stresses are important contributors, particularly during the neap tides. The usefulness of an effective eddy viscosity in the tidally averaged momentum equation is explored, and it is seen that the effective eddy viscosity on the subtidal time scale would need to be negative to close the momentum balance. This is due to the dominant contribution of tidally varying turbulent momentum fluxes, which have no specific relation to the subtidal circulation. Using a water column model, the validity of an effective eddy viscosity is explored; for periodically stratified water columns, a negative effective viscosity is required. ?? 2010 American Meteorological Society.

Calculation of eddy viscosity in a compressible turbulent boundary layer with mass injection and chemical reaction, volume 2. [computer programs

NASA Technical Reports Server (NTRS)

Omori, S.

1973-01-01

As described in Vol. 1, the eddy viscosity is calculated through the turbulent kinetic energy, in order to include the history of the flow and the effect of chemical reaction on boundary layer characteristics. Calculations can be performed for two different cooling concepts; that is, transpiration and regeneratively cooled wall cases. For the regenerative cooling option, coolant and gas side wall temperature and coolant bulk temperature in a rocket engine can be computed along the nozzle axis. Thus, this computer program is useful in designing coolant flow rate and cooling tube geometry, including the tube wall thickness as well as in predicting the effects of boundary layers along the gas side wall on thrust performances.

Calculation of eddy viscosity in a compressible turbulent boundary layer with mass injection and chemical reaction, volume 1. [theoretical analysis

NASA Technical Reports Server (NTRS)

Omori, S.

1973-01-01

The turbulent kinetic energy equation is coupled with boundary layer equations to solve the characteristics of compressible turbulent boundary layers with mass injection and combustion. The Reynolds stress is related to the turbulent kinetic energy using the Prandtl-Wieghardt formulation. When a lean mixture of hydrogen and nitrogen is injected through a porous plate into the subsonic turbulent boundary layer of air flow and ignited by external means, the turbulent kinetic energy increases twice as much as that of noncombusting flow with the same mass injection rate of nitrogen. The magnitudes of eddy viscosity between combusting and noncombusting flows with injection, however, are almost the same due to temperature effects, while the distributions are different. The velocity profiles are significantly affected by combustion; that is, combustion alters the velocity profile as if the mass injection rate is increased, reducing the skin-friction as a result of a smaller velocity gradient at the wall. If pure hydrogen as a transpiration coolant is injected into a rocket nozzle boundary layer flow of combustion products, the temperature drops significantly across the boundary layer due to the high heat capacity of hydrogen. At a certain distance from the wall, hydrogen reacts with the combustion products, liberating an extensive amount of heat. The resulting large increase in temperature reduces the eddy viscosity in this region.

Study of the far wake vortex field generated by a rectangular airfoil in a water tank

NASA Technical Reports Server (NTRS)

Lezius, D. K.

1973-01-01

Underwater towing experiments were carried out with a rectangular airfoil of aspect ratio 5.3 at 4 and 8 deg angles of attack and at chord-based Reynolds numbers between 2 x 100,000 and 7.5 x 100,000. Quantitative measurements by means of the hydrogen bubble technique indicated lower peak swirl velocities in the range of 100 to 1000 lenghts downstream than have been measured in wind tunnel of flight tests. The maximum circumferential velocity decayed whereas the turbulent eddy viscosity increased. This behavior and other known rates of vortex decay are explained in terms of an analytical solution for the vortex problem with time varying eddy viscosity. It is shown that this case corresponds to nonequilibrium turbulent vortex flow.

An Improved K-Epsilon Model for Near-Wall Turbulence and Comparison with Direct Numerical Simulation

NASA Technical Reports Server (NTRS)

Shih, T. H.

1990-01-01

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. The near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation is analyzed. Based on this analysis, a modified eddy viscosity model, having correct near-wall behavior, is suggested, and a model for the pressure transport term in the k-equation is proposed. In addition, a modeled dissipation rate equation is reformulated. Fully developed channel flows were used for model testing. The calculations using various k-epsilon models are compared with direct numerical simulations. The results show that the present k-epsilon model performs well in predicting the behavior of near-wall turbulence. Significant improvement over previous k-epsilon models is obtained.

On the physical basis of succussion.

PubMed

Torres, J L

2002-10-01

It is argued that succussion drives the homeopathic tincture undergoing potentisation to a turbulent regime, where vortices continually form and disappear, ranging in size from the linear extent of the container to a minimum scale determined by viscosity and the rate of energy dissipation. Input mechanical energy cascades down this population of eddies and becomes available at the microscopic level to perform work (chemical, electrical, etc). A structure generated in the tincture would be rupted by vortices smaller than it, and this sets definite limits on the strength of succussion, so the power input leads to larger vortices than the structures one is trying to create and preserve through potentisation. An experimental procedure to test this proposal is suggested, based on Rayleigh scattering.

PRATHAM: Parallel Thermal Hydraulics Simulations using Advanced Mesoscopic Methods

DOE Office of Scientific and Technical Information (OSTI.GOV)

Joshi, Abhijit S; Jain, Prashant K; Mudrich, Jaime A

2012-01-01

At the Oak Ridge National Laboratory, efforts are under way to develop a 3D, parallel LBM code called PRATHAM (PaRAllel Thermal Hydraulic simulations using Advanced Mesoscopic Methods) to demonstrate the accuracy and scalability of LBM for turbulent flow simulations in nuclear applications. The code has been developed using FORTRAN-90, and parallelized using the message passing interface MPI library. Silo library is used to compact and write the data files, and VisIt visualization software is used to post-process the simulation data in parallel. Both the single relaxation time (SRT) and multi relaxation time (MRT) LBM schemes have been implemented in PRATHAM.more » To capture turbulence without prohibitively increasing the grid resolution requirements, an LES approach [5] is adopted allowing large scale eddies to be numerically resolved while modeling the smaller (subgrid) eddies. In this work, a Smagorinsky model has been used, which modifies the fluid viscosity by an additional eddy viscosity depending on the magnitude of the rate-of-strain tensor. In LBM, this is achieved by locally varying the relaxation time of the fluid.« less

The role of bulk viscosity on the decay of compressible, homogeneous, isotropic turbulence

NASA Astrophysics Data System (ADS)

Johnsen, Eric; Pan, Shaowu

2016-11-01

The practice of neglecting bulk viscosity in studies of compressible turbulence is widespread. While exact for monatomic gases and unlikely to strongly affect the dynamics of fluids whose bulk-to-shear viscosity ratio is small and/or of weakly compressible turbulence, this assumption is not justifiable for compressible, turbulent flows of gases whose bulk viscosity is orders of magnitude larger than their shear viscosities (e.g., CO2). To understand the mechanisms by which bulk viscosity and the associated phenomena affect compressible turbulence, we conduct DNS of freely decaying compressible, homogeneous, isotropic turbulence for ratios of bulk-to-shear viscosity ranging from 0-1000. Our simulations demonstrate that bulk viscosity increases the decay rate of turbulent kinetic energy; while enstrophy exhibits little sensitivity to bulk viscosity, dilatation is reduced by an order of magnitude within the two eddy turnover time. Via a Helmholtz decomposition of the flow, we determined that bulk viscosity damps the dilatational velocity and reduces dilatational-solenoidal exchanges, as well as pressure-dilatation coupling. In short, bulk viscosity renders compressible turbulence incompressible by reducing energy transfer between translational and internal modes.

Simulation of Deep Convective Clouds with the Dynamic Reconstruction Turbulence Closure

NASA Astrophysics Data System (ADS)

Shi, X.; Chow, F. K.; Street, R. L.; Bryan, G. H.

2017-12-01

The terra incognita (TI), or gray zone, in simulations is a range of grid spacing comparable to the most energetic eddy diameter. Spacing in mesoscale and simulations is much larger than the eddies, and turbulence is parameterized with one-dimensional vertical-mixing. Large eddy simulations (LES) have grid spacing much smaller than the energetic eddies, and use three-dimensional models of turbulence. Studies of convective weather use convection-permitting resolutions, which are in the TI. Neither mesoscale-turbulence nor LES models are designed for the TI, so TI turbulence parameterization needs to be discussed. Here, the effects of sub-filter scale (SFS) closure schemes on the simulation of deep tropical convection are evaluated by comparing three closures, i.e. Smagorinsky model, Deardorff-type TKE model and the dynamic reconstruction model (DRM), which partitions SFS turbulence into resolvable sub-filter scales (RSFS) and unresolved sub-grid scales (SGS). The RSFS are reconstructed, and the SGS are modeled with a dynamic eddy viscosity/diffusivity model. The RSFS stresses/fluxes allow backscatter of energy/variance via counter-gradient stresses/fluxes. In high-resolution (100m) simulations of tropical convection use of these turbulence models did not lead to significant differences in cloud water/ice distribution, precipitation flux, or vertical fluxes of momentum and heat. When model resolutions are coarsened, the Smagorinsky and TKE models overestimate cloud ice and produces large-amplitude downward heat flux in the middle troposphere (not found in the high-resolution simulations). This error is a result of unrealistically large eddy diffusivities, i.e., the eddy diffusivity of the DRM is on the order of 1 for the coarse resolution simulations, the eddy diffusivity of the Smagorinsky and TKE model is on the order of 100. Splitting the eddy viscosity/diffusivity scalars into vertical and horizontal components by using different length scales and strain rate components helps to reduce the errors, but does not completely remedy the problem. In contrast, the coarse resolution simulations using the DRM produce results that are more consistent with the high-resolution results, suggesting that the DRM is a more appropriate turbulence model for simulating convection in the TI.

right-sized dimple evaluator

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rodriguez, Sal

2017-08-24

The code (aka computer program written as a Matlab script) uses a unique set of n independent equations to solve for n turbulence variables. The code requires the input of a characteristic dimension, a characteristic fluid velocity, the fluid dynamic viscosity, and the fluid density. Most importantly, the code estimates the size of three key turbulent eddies: Kolmogorov, Taylor, and integral. Based on the eddy sizes, dimples dimensions are prescribed such that the key eddies (principally Taylor, and sometimes Kolmogorov), can be generated by the dimple rim and flow unimpeded through the dimple’s concave cavity. It is hypothesized that turbulentmore » eddies are generated by the dimple rim at the dimple-surface interface. The newly-generated eddies in turn entrain the movement of surrounding regions of fluid, creating more mixing. The eddies also generate lift near the wall surrounding the dimple, as they accelerate and reduce pressure in the regions near and at the dimple cavity, thereby minimizing the fluid drag.« less

Very-high-Reynolds-number vortex dynamics via Coherent-vorticity-Preserving (CvP) Large-eddy simulations

NASA Astrophysics Data System (ADS)

Chapelier, Jean-Baptiste; Wasistho, Bono; Scalo, Carlo

2017-11-01

A new approach to Large-Eddy Simulation (LES) is introduced, where subgrid-scale (SGS) dissipation is applied proportionally to the degree of local spectral broadening, hence mitigated in regions dominated by large-scale vortical motion. The proposed CvP-LES methodology is based on the evaluation of the ratio of the test-filtered to resolved (or grid-filtered) enstrophy: σ = ξ ∧ / ξ . Values of σ = 1 indicate low sub-test-filter turbulent activity, justifying local deactivation of any subgrid-scale model. Values of σ < 1 span conditions ranging from incipient spectral broadening σ <= 1 , to equilibrium turbulence σ =σeq < 1 , where σeq is solely as a function of the test-to-grid filter-width ratio Δ ∧ / Δ , derived assuming a Kolmogorov's spectrum. Eddy viscosity is fully restored for σ <=σeq . The proposed approach removes unnecessary SGS dissipation, can be applied to any eddy-viscosity model, is algorithmically simple and computationally inexpensive. A CvP-LES of a pair of unstable helical vortices, representative of rotor-blade wake dynamics, show the ability of the method to sort the coherent motion from the small-scale dynamics. This work is funded by subcontract KSC-17-001 between Purdue University and Kord Technologies, Inc (Huntsville), under the US Navy Contract N68335-17-C-0159 STTR-Phase II, Purdue Proposal No. 00065007, Topic N15A-T002.

Turbulent Flow Structure Inside a Canopy with Complex Multi-Scale Elements

NASA Astrophysics Data System (ADS)

Bai, Kunlun; Katz, Joseph; Meneveau, Charles

2015-06-01

Particle image velocimetry laboratory measurements are carried out to study mean flow distributions and turbulent statistics inside a canopy with complex geometry and multiple scales consisting of fractal, tree-like objects. Matching the optical refractive indices of the tree elements with those of the working fluid provides unobstructed optical paths for both illuminations and image acquisition. As a result, the flow fields between tree branches can be resolved in great detail, without optical interference. Statistical distributions of mean velocity, turbulence stresses, and components of dispersive fluxes are documented and discussed. The results show that the trees leave their signatures in the flow by imprinting wake structures with shapes similar to the trees. The velocities in both wake and non-wake regions significantly deviate from the spatially-averaged values. These local deviations result in strong dispersive fluxes, which are important to account for in canopy-flow modelling. In fact, we find that the streamwise normal dispersive flux inside the canopy has a larger magnitude (by up to four times) than the corresponding Reynolds normal stress. Turbulent transport in horizontal planes is studied in the framework of the eddy viscosity model. Scatter plots comparing the Reynolds shear stress and mean velocity gradient are indicative of a linear trend, from which one can calculate the eddy viscosity and mixing length. Similar to earlier results from the wake of a single tree, here we find that inside the canopy the mean mixing length decreases with increasing elevation. This trend cannot be scaled based on a single length scale, but can be described well by a model, which considers the coexistence of multi-scale branches. This agreement indicates that the multi-scale information and the clustering properties of the fractal objects should be taken into consideration in flows inside multi-scale canopies.

Artificial fluid properties for large-eddy simulation of compressible turbulent mixing

NASA Astrophysics Data System (ADS)

Cook, Andrew W.

2007-05-01

An alternative methodology is described for large-eddy simulation (LES) of flows involving shocks, turbulence, and mixing. In lieu of filtering the governing equations, it is postulated that the large-scale behavior of a LES fluid, i.e., a fluid with artificial properties, will be similar to that of a real fluid, provided the artificial properties obey certain constraints. The artificial properties consist of modifications to the shear viscosity, bulk viscosity, thermal conductivity, and species diffusivity of a fluid. The modified transport coefficients are designed to damp out high wavenumber modes, close to the resolution limit, without corrupting lower modes. Requisite behavior of the artificial properties is discussed and results are shown for a variety of test problems, each designed to exercise different aspects of the models. When combined with a tenth-order compact scheme, the overall method exhibits excellent resolution characteristics for turbulent mixing, while capturing shocks and material interfaces in a crisp fashion.

Wave attenuation in the marginal ice zone during LIMEX

NASA Technical Reports Server (NTRS)

Liu, Antony K.; Vachon, Paris W.; Peng, Chih Y.; Bhogal, A. S.

1992-01-01

The effect of ice cover on ocean-wave attenuation is investigated for waves under flexure in the marginal ice zone (MIZ) with SAR image spectra and the results of models. Directional wavenumber spectra are taken from the SAR image data, and the wave-attenuation rate is evaluated with SAR image spectra and by means of the model by Liu and Mollo-Christensen (1988). Eddy viscosity is described by means of dimensional analysis as a function of ice roughness and wave-induced velocity, and comparisons are made with the remotely sensed data. The model corrects the open-water model by introducing the effects of a continuous ice sheet, and turbulent eddy viscosity is shown to depend on ice thickness, floe sizes, significant wave height, and wave period. SAR and wave-buoy data support the trends described in the model results, and a characteristic rollover is noted in the model and experimental wave-attenuation rates at high wavenumbers.

Computation of turbulent rotating channel flow with an algebraic Reynolds stress model

NASA Technical Reports Server (NTRS)

Warfield, M. J.; Lakshminarayana, B.

1986-01-01

An Algebraic Reynolds Stress Model has been implemented to modify the Kolmogorov-Prandtl eddy viscosity relation to produce an anisotropic turbulence model. The eddy viscosity relation becomes a function of the local turbulent production to dissipation ratio and local turbulence/rotation parameters. The model is used to predict fully-developed rotating channel flow over a diverse range of rotation numbers. In addition, predictions are obtained for a developing channel flow with high rotation. The predictions are compared with the experimental data available. Good predictions are achieved for mean velocity and wall shear stress over most of the rotation speeds tested. There is some prediction breakdown at high rotation (rotation number greater than .10) where the effects of the rotation on turbulence become quite complex. At high rotation and low Reynolds number, the laminarization on the trailing side represents a complex effect of rotation which is difficult to predict with the described models.

Design optimization using adjoint of Long-time LES for the trailing edge of a transonic turbine vane

NASA Astrophysics Data System (ADS)

Talnikar, Chaitanya; Wang, Qiqi

2017-11-01

Adjoint-based design optimization methods have been applied to low-fidelity simulation methods like Reynolds Averaged Navier-Stokes (RANS) and are useful for designing fluid machinery components. But to reliably capture the complex flow phenomena involved in turbomachinery, high fidelity simulations like large eddy simulation (LES) are required. Unfortunately due to the chaotic dynamics of turbulence, the unsteady adjoint method for LES diverges and produces incorrect gradients. Using a viscosity stabilized unsteady adjoint method developed for LES, the gradient can be obtained with reasonable accuracy. In this paper, design of the trailing edge of a gas turbine inlet guide vane is performed with the objective to reduce stagnation pressure loss and heat transfer over the surface of the vane. Slight changes in the shape of trailing edge can significantly impact these quantities by altering the boundary layer development process and separation points. The trailing edge is parameterized using a linear combination of 5 convex designs. Bayesian optimization is used as a global optimizer with the objective function evaluated from the LES and gradients obtained using the viscosity adjoint method. Results from the optimization, performed on the supercomputer Mira, are presented.

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 3: Programmer's reference

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This User's Guide describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

The Influence of Slope Breaks on Lava Flow Surface Disruption

NASA Technical Reports Server (NTRS)

Glaze, Lori S.; Baloga, Stephen M.; Fagents, Sarah A.; Wright, Robert

2014-01-01

Changes in the underlying slope of a lava flow impart a significant fraction of rotational energy beyond the slope break. The eddies, circulation and vortices caused by this rotational energy can disrupt the flow surface, having a significant impact on heat loss and thus the distance the flow can travel. A basic mechanics model is used to compute the rotational energy caused by a slope change. The gain in rotational energy is deposited into an eddy of radius R whose energy is dissipated as it travels downstream. A model of eddy friction with the ambient lava is used to compute the time-rate of energy dissipation. The key parameter of the dissipation rate is shown to be rho R(sup 2/)mu, where ? is the lava density and mu is the viscosity, which can vary by orders of magnitude for different flows. The potential spatial disruption of the lava flow surface is investigated by introducing steady-state models for the main flow beyond the steepening slope break. One model applies to slow-moving flows with both gravity and pressure as the driving forces. The other model applies to fast-moving, low-viscosity, turbulent flows. These models provide the flow velocity that establishes the downstream transport distance of disrupting eddies before they dissipate. The potential influence of slope breaks is discussed in connection with field studies of lava flows from the 1801 Hualalai and 1823 Keaiwa Kilauea, Hawaii, and 2004 Etna eruptions.

Comparative Study Of Four Models Of Turbulence

NASA Technical Reports Server (NTRS)

Menter, Florian R.

1996-01-01

Report presents comparative study of four popular eddy-viscosity models of turbulence. Computations reported for three different adverse pressure-gradient flowfields. Detailed comparison of numerical results and experimental data given. Following models tested: Baldwin-Lomax, Johnson-King, Baldwin-Barth, and Wilcox.

Wall function treatment for bubbly boundary layers at low void fractions.

PubMed

Soares, Daniel V; Bitencourt, Marcelo C; Loureiro, Juliana B R; Silva Freire, Atila P

2018-01-01

The present work investigates the role of different treatments of the lower boundary condition on the numerical prediction of bubbly flows. Two different wall function formulations are tested against experimental data obtained for bubbly boundary layers: (i) a new analytical solution derived through asymptotic techniques and (ii) the previous formulation of Troshko and Hassan (IJHMT, 44, 871-875, 2001a). A modified k-e model is used to close the averaged Navier-Stokes equations together with the hypothesis that turbulence can be modelled by a linear superposition of bubble and shear induced eddy viscosities. The work shows, in particular, how four corrections must the implemented in the standard single-phase k-e model to account for the effects of bubbles. The numerical implementation of the near wall functions is made through a finite elements code.

Ocean tides for satellite geodesy

NASA Technical Reports Server (NTRS)

Dickman, S. R.

1990-01-01

Spherical harmonic tidal solutions have been obtained at the frequencies of the 32 largest luni-solar tides using prior theory of the author. That theory was developed for turbulent, nonglobal, self-gravitating, and loading oceans possessing realistic bathymetry and linearized bottom friction; the oceans satisfy no-flow boundary conditions at coastlines. In this theory the eddy viscosity and bottom drag coefficients are treated as spatially uniform. Comparison of the predicted degree-2 components of the Mf, P1, and M2 tides with those from numerical and satellite-based tide models allows the ocean friction parameters to be estimated at long and short periods. Using the 32 tide solutions, the frequency dependence of tidal admittance is investigated, and the validity of sideband tide models used in satellite orbit analysis is examined. The implications of admittance variability for oceanic resonances are also explored.

Study of Near-Surface Models in Large-Eddy Simulations of a Neutrally Stratified Atmospheric Boundary Layer

NASA Technical Reports Server (NTRS)

Senocak, I.; Ackerman, A. S.; Kirkpatrick, M. P.; Stevens, D. E.; Mansour, N. N.

2004-01-01

Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.

Nonlinear closures for scale separation in supersonic magnetohydrodynamic turbulence

NASA Astrophysics Data System (ADS)

Grete, Philipp; Vlaykov, Dimitar G.; Schmidt, Wolfram; Schleicher, Dominik R. G.; Federrath, Christoph

2015-02-01

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment—magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures—models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The SGS energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ({{β }p}=0.25,1,2.5,5,25) and sonic Mach numbers ({{M}s}=2,2.5,4). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and α -β -γ closures. We find only mediocre performance of the kinetic eddy-viscosity and α -β -γ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlations and constant coefficient values over time and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

VALIDITY OF A TWO-DIMENSIONAL MODEL FOR VARIABLE-DENSITY HYDRODYNAMIC CIRCULATION

EPA Science Inventory

A three-dimensional model of temperatures and currents has been formulated to assist in the analysis and interpretation of the dynamics of stratified lakes. In this model, nonlinear eddy coefficients for viscosity and conductivities are included. A two-dimensional model (one vert...

A novel VLES model accounting for near-wall turbulence: physical rationale and applications

NASA Astrophysics Data System (ADS)

Jakirlic, Suad; Chang, Chi-Yao; Kutej, Lukas; Tropea, Cameron

2014-11-01

A novel VLES (Very Large Eddy Simulation) model whose non-resolved residual turbulence is modelled by using an advanced near-wall eddy-viscosity model accounting for the near-wall Reynolds stress anisotropy influence on the turbulence viscosity by modelling appropriately the velocity scale in the relevant formulation (Hanjalic et al., 2004) is proposed. It represents a variable resolution Hybrid LES/RANS (Reynolds-Averaged Navier-Stokes) computational scheme enabling a seamless transition from RANS to LES depending on the ratio of the turbulent viscosities associated with the unresolved scales corresponding to the LES cut-off and the `unsteady' scales pertinent to the turbulent properties of the VLES residual motion, which varies within the flow domain. The VLES method is validated interactively in the process of the model derivation by computing fully-developed flow in a plane channel (important representative of wall-bounded flows, underlying the log-law for the velocity field, for studying near-wall Reynolds stress anisotropy) and a separating flow over a periodic arrangement of smoothly-contoured 2-D hills. The model performances are also assessed in capturing the natural decay of the homogeneous isotropic turbulence. The model is finally applied to swirling flow in a vortex tube, flow in an IC-engine configuration and flow past a realistic car model.

A new time scale based k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1992-01-01

A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

Near-wall k-epsilon turbulence modeling

NASA Technical Reports Server (NTRS)

Mansour, N. N.; Kim, J.; Moin, P.

1987-01-01

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

New time scale based k-epsilon model for near-wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1993-01-01

A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

A defect stream function, law of the wall/wake method for compressible turbulent boundary layers

NASA Technical Reports Server (NTRS)

Barnwell, Richard W.; Dejarnette, Fred R.; Wahls, Richard A.

1989-01-01

The application of the defect stream function to the solution of the two-dimensional, compressible boundary layer is examined. A law of the wall/law of the wake formulation for the inner part of the boundary layer is presented which greatly simplifies the computational task near the wall and eliminates the need for an eddy viscosity model in this region. The eddy viscosity model in the outer region is arbitrary. The modified Crocco temperature-velocity relationship is used as a simplification of the differential energy equation. Formulations for both equilibrium and nonequilibrium boundary layers are presented including a constrained zero-order form which significantly reduces the computational workload while retaining the significant physics of the flow. A formulation for primitive variables is also presented. Results are given for the constrained zero-order and second-order equilibrium formulations and are compared with experimental data. A compressible wake function valid near the wall has been developed from the present results.

Some effects of horizontal discretization on linear baroclinic and symmetric instabilities

NASA Astrophysics Data System (ADS)

Barham, William; Bachman, Scott; Grooms, Ian

2018-05-01

The effects of horizontal discretization on linear baroclinic and symmetric instabilities are investigated by analyzing the behavior of the hydrostatic Eady problem in ocean models on the B and C grids. On the C grid a spurious baroclinic instability appears at small wavelengths. This instability does not disappear as the grid scale decreases; instead, it simply moves to smaller horizontal scales. The peak growth rate of the spurious instability is independent of the grid scale as the latter decreases. It is equal to cf /√{Ri} where Ri is the balanced Richardson number, f is the Coriolis parameter, and c is a nondimensional constant that depends on the Richardson number. As the Richardson number increases c increases towards an upper bound of approximately 1/2; for large Richardson numbers the spurious instability is faster than the Eady instability. To suppress the spurious instability it is recommended to use fourth-order centered tracer advection along with biharmonic viscosity and diffusion with coefficients (Δx) 4 f /(32√{Ri}) or larger where Δx is the grid scale. On the B grid, the growth rates of baroclinic and symmetric instabilities are too small, and converge upwards towards the correct values as the grid scale decreases; no spurious instabilities are observed. In B grid models at eddy-permitting resolution, the reduced growth rate of baroclinic instability may contribute to partially-resolved eddies being too weak. On the C grid the growth rate of symmetric instability is better (larger) than on the B grid, and converges upwards towards the correct value as the grid scale decreases.

Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES

NASA Astrophysics Data System (ADS)

Marras, Simone; Nazarov, Murtazo; Giraldo, Francis X.

2015-11-01

The high order spectral element approximation of the Euler equations is stabilized via a dynamic sub-grid scale model (Dyn-SGS). This model was originally designed for linear finite elements to solve compressible flows at large Mach numbers. We extend its application to high-order spectral elements to solve the Euler equations of low Mach number stratified flows. The major justification of this work is twofold: stabilization and large eddy simulation are achieved via one scheme only. Because the diffusion coefficients of the regularization stresses obtained via Dyn-SGS are residual-based, the effect of the artificial diffusion is minimal in the regions where the solution is smooth. The direct consequence is that the nominal convergence rate of the high-order solution of smooth problems is not degraded. To our knowledge, this is the first application in atmospheric modeling of a spectral element model stabilized by an eddy viscosity scheme that, by construction, may fulfill stabilization requirements, can model turbulence via LES, and is completely free of a user-tunable parameter. From its derivation, it will be immediately clear that Dyn-SGS is independent of the numerical method; it could be implemented in a discontinuous Galerkin, finite volume, or other environments alike. Preliminary discontinuous Galerkin results are reported as well. The straightforward extension to non-linear scalar problems is also described. A suite of 1D, 2D, and 3D test cases is used to assess the method, with some comparison against the results obtained with the most known Lilly-Smagorinsky SGS model.

A Lagrangian subgrid-scale model with dynamic estimation of Lagrangian time scale for large eddy simulation of complex flows

NASA Astrophysics Data System (ADS)

Verma, Aman; Mahesh, Krishnan

2012-08-01

The dynamic Lagrangian averaging approach for the dynamic Smagorinsky model for large eddy simulation is extended to an unstructured grid framework and applied to complex flows. The Lagrangian time scale is dynamically computed from the solution and does not need any adjustable parameter. The time scale used in the standard Lagrangian model contains an adjustable parameter θ. The dynamic time scale is computed based on a "surrogate-correlation" of the Germano-identity error (GIE). Also, a simple material derivative relation is used to approximate GIE at different events along a pathline instead of Lagrangian tracking or multi-linear interpolation. Previously, the time scale for homogeneous flows was computed by averaging along directions of homogeneity. The present work proposes modifications for inhomogeneous flows. This development allows the Lagrangian averaged dynamic model to be applied to inhomogeneous flows without any adjustable parameter. The proposed model is applied to LES of turbulent channel flow on unstructured zonal grids at various Reynolds numbers. Improvement is observed when compared to other averaging procedures for the dynamic Smagorinsky model, especially at coarse resolutions. The model is also applied to flow over a cylinder at two Reynolds numbers and good agreement with previous computations and experiments is obtained. Noticeable improvement is obtained using the proposed model over the standard Lagrangian model. The improvement is attributed to a physically consistent Lagrangian time scale. The model also shows good performance when applied to flow past a marine propeller in an off-design condition; it regularizes the eddy viscosity and adjusts locally to the dominant flow features.

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

NASA Technical Reports Server (NTRS)

Horiuti, K.; Mansour, N. N.; Kim, John J.

1993-01-01

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Analysis of turbulent heat transfer, mass transfer, and friction in smooth tubes at high Prandtl and Schmidt numbers

NASA Technical Reports Server (NTRS)

Deissler, Robert G

1955-01-01

The expression for eddy diffusivity from a previous analysis was modified in order to account for the effect of kinematic viscosity on the turbulence in the region close to a wall. By using the modified expression, good agreement was obtained between predicted and experimental results for heat and mass transfer at Prandtl and Schmidt numbers between 0.5 and 3000. The effects of length-to-diameter ratio and of variable viscosity were also investigated for a wide range of Prandtl numbers.

Interactions between a fractal tree-like object and hydrodynamic turbulence: flow structure and characteristic mixing length

NASA Astrophysics Data System (ADS)

Meneveau, C. V.; Bai, K.; Katz, J.

2011-12-01

The vegetation canopy has a significant impact on various physical and biological processes such as forest microclimate, rainfall evaporation distribution and climate change. Most scaled laboratory experimental studies have used canopy element models that consist of rigid vertical strips or cylindrical rods that can be typically represented through only one or a few characteristic length scales, for example the diameter and height for cylindrical rods. However, most natural canopies and vegetation are highly multi-scale with branches and sub-branches, covering a wide range of length scales. Fractals provide a convenient idealization of multi-scale objects, since their multi-scale properties can be described in simple ways (Mandelbrot 1982). While fractal aspects of turbulence have been studied in several works in the past decades, research on turbulence generated by fractal objects started more recently. We present an experimental study of boundary layer flow over fractal tree-like objects. Detailed Particle-Image-Velocimetry (PIV) measurements are carried out in the near-wake of a fractal-like tree. The tree is a pre-fractal with five generations, with three branches and a scale reduction factor 1/2 at each generation. Its similarity fractal dimension (Mandelbrot 1982) is D ~ 1.58. Detailed mean velocity and turbulence stress profiles are documented, as well as their downstream development. We then turn attention to the turbulence mixing properties of the flow, specifically to the question whether a mixing length-scale can be identified in this flow, and if so, how it relates to the geometric length-scales in the pre-fractal object. Scatter plots of mean velocity gradient (shear) and Reynolds shear stress exhibit good linear relation at all locations in the flow. Therefore, in the transverse direction of the wake evolution, the Boussinesq eddy viscosity concept is appropriate to describe the mixing. We find that the measured mixing length increases with increasing streamwise locations. Conversely, the measured eddy viscosity and mixing length decrease with increasing elevation, which differs from eddy viscosity and mixing length behaviors of traditional boundary layers or canopies studied before. In order to find an appropriate length for the flow, several models based on the notion of superposition of scales are proposed and examined. One approach is based on spectral distributions. Another more practical approach is based on length-scale distributions evaluated using fractal geometry tools. These proposed models agree well with the measured mixing length. The results indicate that information about multi-scale clustering of branches as it occurs in fractals has to be incorporated into models of the mixing length for flows through canopies with multiple scales. The research is supported by National Science Foundation grant ATM-0621396 and AGS-1047550.

An interior penalty stabilised incompressible discontinuous Galerkin-Fourier solver for implicit large eddy simulations

NASA Astrophysics Data System (ADS)

Ferrer, Esteban

2017-11-01

We present an implicit Large Eddy Simulation (iLES) h / p high order (≥2) unstructured Discontinuous Galerkin-Fourier solver with sliding meshes. The solver extends the laminar version of Ferrer and Willden, 2012 [34], to enable the simulation of turbulent flows at moderately high Reynolds numbers in the incompressible regime. This solver allows accurate flow solutions of the laminar and turbulent 3D incompressible Navier-Stokes equations on moving and static regions coupled through a high order sliding interface. The spatial discretisation is provided by the Symmetric Interior Penalty Discontinuous Galerkin (IP-DG) method in the x-y plane coupled with a purely spectral method that uses Fourier series and allows efficient computation of spanwise periodic three-dimensional flows. Since high order methods (e.g. discontinuous Galerkin and Fourier) are unable to provide enough numerical dissipation to enable under-resolved high Reynolds computations (i.e. as necessary in the iLES approach), we adapt the laminar version of the solver to increase (controllably) the dissipation and enhance the stability in under-resolved simulations. The novel stabilisation relies on increasing the penalty parameter included in the DG interior penalty (IP) formulation. The latter penalty term is included when discretising the linear viscous terms in the incompressible Navier-Stokes equations. These viscous penalty fluxes substitute the stabilising effect of non-linear fluxes, which has been the main trend in implicit LES discontinuous Galerkin approaches. The IP-DG penalty term provides energy dissipation, which is controlled by the numerical jumps at element interfaces (e.g. large in under-resolved regions) such as to stabilise under-resolved high Reynolds number flows. This dissipative term has minimal impact in well resolved regions and its implicit treatment does not restrict the use of large time steps, thus providing an efficient stabilization mechanism for iLES. The IP-DG stabilisation is complemented with a Spectral Vanishing Viscosity (SVV) method, in the z-direction, to enhance stability in the continuous Fourier space. The coupling between the numerical viscosity in the DG plane and the SVV damping, provides an efficient approach to stabilise high order methods at moderately high Reynolds numbers. We validate the formulation for three turbulent flow cases: a circular cylinder at Re = 3900, a static and pitch oscillating NACA 0012 airfoil at Re = 10000 and finally a rotating vertical-axis turbine at Re = 40000, with Reynolds based on the circular diameter, airfoil chord and turbine diameter, respectively. All our results compare favourably with published direct numerical simulations, large eddy simulations or experimental data. We conclude that the DG-Fourier high order solver, with IP-SVV stabilisation, proves to be a valuable tool to predict turbulent flows and associated statistics for both static and rotating machinery.

Evaluation of scale-aware subgrid mesoscale eddy models in a global eddy-rich model

NASA Astrophysics Data System (ADS)

Pearson, Brodie; Fox-Kemper, Baylor; Bachman, Scott; Bryan, Frank

2017-07-01

Two parameterizations for horizontal mixing of momentum and tracers by subgrid mesoscale eddies are implemented in a high-resolution global ocean model. These parameterizations follow on the techniques of large eddy simulation (LES). The theory underlying one parameterization (2D Leith due to Leith, 1996) is that of enstrophy cascades in two-dimensional turbulence, while the other (QG Leith) is designed for potential enstrophy cascades in quasi-geostrophic turbulence. Simulations using each of these parameterizations are compared with a control simulation using standard biharmonic horizontal mixing.Simulations using the 2D Leith and QG Leith parameterizations are more realistic than those using biharmonic mixing. In particular, the 2D Leith and QG Leith simulations have more energy in resolved mesoscale eddies, have a spectral slope more consistent with turbulence theory (an inertial enstrophy or potential enstrophy cascade), have bottom drag and vertical viscosity as the primary sinks of energy instead of lateral friction, and have isoneutral parameterized mesoscale tracer transport. The parameterization choice also affects mass transports, but the impact varies regionally in magnitude and sign.

Studies of the relationship between environmental forcing and the structure and dynamics of tornado-like vortices

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nolan, David S.; Almgren, Ann S.; Bell, John B.

Axisymmetric numerical simulations continue to provide insight into how the structure, dynamics, and maximum wind speeds of tornadoes, and other convectively-maintained vortices, are influenced by the surrounding environment. This work is continued with a new numerical model of axisymmetric incompressible flow that incorporates adaptive mesh refinement. The model dynamically increases or decreases the resolution in regions of interest as determined by a specified refinement criterion. Here, the criterion used is based on the cell Reynolds number dx dv / nu, so that the flow is guaranteed to be laminar on the scale of the local grid spacing. The model ismore » used to investigate how the altitude and shape of the convective forcing, the size of the domain, and the effective Reynolds number (based on the choice of the eddy viscosity nu) influence the structure and dynamics of the vortex. Over a wide variety of domain and forcing geometries,the vortex Reynolds number Gamma / nu (the ratio of the far-field circulation to the eddy viscosity) is shown to be the most important parameter for determining vortex structure and behavior. Furthermore,it is found that the vertical scale of the convective forcing only affects the vortex inasmuch as this vertical scale contributes to the total strength of the convective forcing. The horizontal scale of the convective forcing, however, is found to be the fundamental length scale in the problem, in that it can determine both the circulation of the fluid that is drawn into the vortex core, and also influences the depth of the swirling boundary layer. Higher mean wind speeds are sustained as the eddy viscosity is decreased; however, it is observed that the highest wind speeds are found in the high-swirl, two-celled vortex regime rather than in the low-swirl, one-celled regime, which is in contrast with some previous results. The conclusions drawn from these results are applied to dimensional simulations with scales similar to the mesocyclone/thunderstorm environment. Tornado-like vortices are reproduced, using a constant eddy viscosity with such values as 40 m2s-1, which have maximum wind speeds, radii of maximum winds, and boundary layer depths which are quite similar to those recently observed with portable Doppler radar. Based on the results of both nondimensional and tornado-scale simulations, scaling laws are empirically derived for the internal length scales in tornado-like vortices, such as the depth of the boundary layer and the radius of maximum winds.« less

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

NASA Technical Reports Server (NTRS)

Murman, Scott M.; Childs, Robert E.; Garcia, Joseph A.

2015-01-01

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Subgrid or Reynolds stress-modeling for three-dimensional turbulence computations

NASA Technical Reports Server (NTRS)

Rubesin, M. W.

1975-01-01

A review is given of recent advances in two distinct computational methods for evaluating turbulence fields, namely, statistical Reynolds stress modeling and turbulence simulation, where large eddies are followed in time. It is shown that evaluation of the mean Reynolds stresses, rather than use of a scalar eddy viscosity, permits an explanation of streamline curvature effects found in several experiments. Turbulence simulation, with a new volume averaging technique and third-order accurate finite-difference computing is shown to predict the decay of isotropic turbulence in incompressible flow with rather modest computer storage requirements, even at Reynolds numbers of aerodynamic interest.

Proteus two-dimensional Navier-Stokes computer code, version 2.0. Volume 2: User's guide

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Bui, Trong T.

1993-01-01

A computer code called Proteus 2D was developed to solve the two-dimensional planar or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body-fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. This is the User's Guide, and describes the program's features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

On Eddy Viscosity, Energy Cascades, and the Horizontal Resolution of Gridded Satellite Altimeter Products

DTIC Science & Technology

2013-02-01

Manuscript received 21 December 2011, in final form 23 August 2012) ABSTRACT Motivated by the recent interest in ocean energetics, the widespread use...Inhomogeneous two-dimensional turbu- lence in the atmosphere. Advances in Turbulence, G. Comte - Bellot and J. Mathieu, Eds., Springer-Verlag, 269-278

SGS Dynamics and Modeling near a Rough Wall.

NASA Astrophysics Data System (ADS)

Juneja, Anurag; Brasseur, James G.

1998-11-01

Large-eddy simulation (LES) of the atmospheric boundary layer (ABL) using classical subgrid-scale (SGS) models is known to poorly predict mean shear at the first few grid cells near the rough surface, creating error which can propogate vertically to infect the entire ABL. Our goal was to determine the first-order errors in predicted SGS terms that arise as a consequence of necessary under-resolution of integral scales and anisotropy which exist at the first few grid levels in LES of rough wall turbulence. Analyzing the terms predicted from eddy-viscosity and similarity closures with DNS anisotropic datasets of buoyancy- and shear-driven turbulence, we uncover three important issues which should be addressed in the design of SGS closures for rough walls and we provide a priori tests for the SGS model. Firstly, we identify a strong spurious coupling between the anisotropic structure of the resolved velocity field and predicted SGS dynamics which can create a feedback loop to incorrectly enhance certain components of the predicted resolved velocity. Secondly, we find that eddy viscosity and similarity SGS models do not contain enough degrees of freedom to capture, at a sufficient level of accuracy, both RS-SGS energy flux and SGS-RS dynamics. Thirdly, to correctly capture pressure transport near a wall, closures must be made more flexible to accommodate proper partitioning between SGS stress divergence and SGS pressure gradient.

Scale-Similar Models for Large-Eddy Simulations

NASA Technical Reports Server (NTRS)

Sarghini, F.

1999-01-01

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers

NASA Technical Reports Server (NTRS)

Kinney, R. B.; Taslim, M. E.; Hung, S. C.

1985-01-01

The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long.

3-D residual eddy current field characterisation: applied to diffusion weighted magnetic resonance imaging.

PubMed

O'Brien, Kieran; Daducci, Alessandro; Kickler, Nils; Lazeyras, Francois; Gruetter, Rolf; Feiweier, Thorsten; Krueger, Gunnar

2013-08-01

Clinical use of the Stejskal-Tanner diffusion weighted images is hampered by the geometric distortions that result from the large residual 3-D eddy current field induced. In this work, we aimed to predict, using linear response theory, the residual 3-D eddy current field required for geometric distortion correction based on phantom eddy current field measurements. The predicted 3-D eddy current field induced by the diffusion-weighting gradients was able to reduce the root mean square error of the residual eddy current field to ~1 Hz. The model's performance was tested on diffusion weighted images of four normal volunteers, following distortion correction, the quality of the Stejskal-Tanner diffusion-weighted images was found to have comparable quality to image registration based corrections (FSL) at low b-values. Unlike registration techniques the correction was not hindered by low SNR at high b-values, and results in improved image quality relative to FSL. Characterization of the 3-D eddy current field with linear response theory enables the prediction of the 3-D eddy current field required to correct eddy current induced geometric distortions for a wide range of clinical and high b-value protocols.

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.; Elmiligui, Alaa; Giriamaji, Sharath S.

2008-01-01

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

A Coherent vorticity preserving eddy-viscosity correction for Large-Eddy Simulation

NASA Astrophysics Data System (ADS)

Chapelier, J.-B.; Wasistho, B.; Scalo, C.

2018-04-01

This paper introduces a new approach to Large-Eddy Simulation (LES) where subgrid-scale (SGS) dissipation is applied proportionally to the degree of local spectral broadening, hence mitigated or deactivated in regions dominated by large-scale and/or laminar vortical motion. The proposed coherent-vorticity preserving (CvP) LES methodology is based on the evaluation of the ratio of the test-filtered to resolved (or grid-filtered) enstrophy, σ. Values of σ close to 1 indicate low sub-test-filter turbulent activity, justifying local deactivation of the SGS dissipation. The intensity of the SGS dissipation is progressively increased for σ < 1 which corresponds to a small-scale spectral broadening. The SGS dissipation is then fully activated in developed turbulence characterized by σ ≤σeq, where the value σeq is derived assuming a Kolmogorov spectrum. The proposed approach can be applied to any eddy-viscosity model, is algorithmically simple and computationally inexpensive. LES of Taylor-Green vortex breakdown demonstrates that the CvP methodology improves the performance of traditional, non-dynamic dissipative SGS models, capturing the peak of total turbulent kinetic energy dissipation during transition. Similar accuracy is obtained by adopting Germano's dynamic procedure albeit at more than twice the computational overhead. A CvP-LES of a pair of unstable periodic helical vortices is shown to predict accurately the experimentally observed growth rate using coarse resolutions. The ability of the CvP methodology to dynamically sort the coherent, large-scale motion from the smaller, broadband scales during transition is demonstrated via flow visualizations. LES of compressible channel are carried out and show a good match with a reference DNS.

Large-Eddy Simulations of Radiatively Driven Convection: Sensitivities to the Representation of Small Scales.

NASA Astrophysics Data System (ADS)

Stevens, Bjorn; Moeng, Chin-Hoh; Sullivan, Peter P.

1999-12-01

Large-eddy simulations of a smoke cloud are examined with respect to their sensitivity to small scales as manifest in either the grid spacing or the subgrid-scale (SGS) model. Calculations based on a Smagorinsky SGS model are found to be more sensitive to the effective resolution of the simulation than are calculations based on the prognostic turbulent kinetic energy (TKE) SGS model. The difference between calculations based on the two SGS models is attributed to the advective transport, diffusive transport, and/or time-rate-of-change terms in the TKE equation. These terms are found to be leading order in the entrainment zone and allow the SGS TKE to behave in a way that tends to compensate for changes that result in larger or smaller resolved scale entrainment fluxes. This compensating behavior of the SGS TKE model is attributed to the fact that changes that reduce the resolved entrainment flux (viz., values of the eddy viscosity in the upper part of the PBL) simultaneously tend to increase the buoyant production of SGS TKE in the radiatively destabilized portion of the smoke cloud. Increased production of SGS TKE in this region then leads to increased amounts of transported, or fossil, SGS TKE in the entrainment zone itself, which in turn leads to compensating increases in the SGS entrainment fluxes. In the Smagorinsky model, the absence of a direct connection between SGS TKE in the entrainment and radiatively destabilized zones prevents this compensating mechanism from being active, and thus leads to calculations whose entrainment rate sensitivities as a whole reflect the sensitivities of the resolved-scale fluxes to values of upper PBL eddy viscosities.

PSF mapping-based correction of eddy-current-induced distortions in diffusion-weighted echo-planar imaging.

PubMed

In, Myung-Ho; Posnansky, Oleg; Speck, Oliver

2016-05-01

To accurately correct diffusion-encoding direction-dependent eddy-current-induced geometric distortions in diffusion-weighted echo-planar imaging (DW-EPI) and to minimize the calibration time at 7 Tesla (T). A point spread function (PSF) mapping based eddy-current calibration method is newly presented to determine eddy-current-induced geometric distortions even including nonlinear eddy-current effects within the readout acquisition window. To evaluate the temporal stability of eddy-current maps, calibration was performed four times within 3 months. Furthermore, spatial variations of measured eddy-current maps versus their linear superposition were investigated to enable correction in DW-EPIs with arbitrary diffusion directions without direct calibration. For comparison, an image-based eddy-current correction method was additionally applied. Finally, this method was combined with a PSF-based susceptibility-induced distortion correction approach proposed previously to correct both susceptibility and eddy-current-induced distortions in DW-EPIs. Very fast eddy-current calibration in a three-dimensional volume is possible with the proposed method. The measured eddy-current maps are very stable over time and very similar maps can be obtained by linear superposition of principal-axes eddy-current maps. High resolution in vivo brain results demonstrate that the proposed method allows more efficient eddy-current correction than the image-based method. The combination of both PSF-based approaches allows distortion-free images, which permit reliable analysis in diffusion tensor imaging applications at 7T. © 2015 Wiley Periodicals, Inc.

Renormalization group estimates of transport coefficients in the advection of a passive scalar by incompressible turbulence

NASA Technical Reports Server (NTRS)

Zhou, YE; Vahala, George

1993-01-01

The advection of a passive scalar by incompressible turbulence is considered using recursive renormalization group procedures in the differential sub grid shell thickness limit. It is shown explicitly that the higher order nonlinearities induced by the recursive renormalization group procedure preserve Galilean invariance. Differential equations, valid for the entire resolvable wave number k range, are determined for the eddy viscosity and eddy diffusivity coefficients, and it is shown that higher order nonlinearities do not contribute as k goes to 0, but have an essential role as k goes to k(sub c) the cutoff wave number separating the resolvable scales from the sub grid scales. The recursive renormalization transport coefficients and the associated eddy Prandtl number are in good agreement with the k-dependent transport coefficients derived from closure theories and experiments.

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

NASA Technical Reports Server (NTRS)

Gurta, R. N.; Trimpi, R. L.

1974-01-01

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.

Computation of turbulent boundary layer flows with an algebraic stress turbulence model

NASA Technical Reports Server (NTRS)

Kim, Sang-Wook; Chen, Yen-Sen

1986-01-01

An algebraic stress turbulence model is presented, characterized by the following: (1) the eddy viscosity expression is derived from the Reynolds stress turbulence model; (2) the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale; and (3) the diffusion coefficients for turbulence equations are adjusted so that the kinetic energy profile extends further into the free stream region found in most experimental data. The turbulent flow equations were solved using a finite element method. Examples include: fully developed channel flow, fully developed pipe flow, flat plate boundary layer flow, plane jet exhausting into a moving stream, circular jet exhausting into a moving stream, and wall jet flow. Computational results compare favorably with experimental data for most of the examples considered. Significantly improved results were obtained for the plane jet flow, the circular jet flow, and the wall jet flow; whereas the remainder are comparable to those obtained by finite difference methods using the standard kappa-epsilon turbulence model. The latter seems to be promising with further improvement of the expression for the eddy viscosity coefficient.

An experimental investigation of turbulent boundary layers along curved surfaces

NASA Technical Reports Server (NTRS)

So, R. M. C.; Mellor, G. L.

1972-01-01

A curved wall tunnel was designed, and an equilibrium turbulent boundary layer was set up on the straight section preceding the curved test section. Turbulent boundary layer flows with uniform and adverse pressure distributions along convex and concave walls were investigated. Hot-wire measurements along the convex surface indicated that turbulent mixing between fluid layers was very much reduced. However, the law of the wall held and the skin friction, thus determined, correlated well with other measurements. Hot-wire measurements along the concave test wall revealed a system of longitudinal vortices inside the boundary layer and confirmed that concave curvature enhances mixing. A self-consistent set of turbulent boundary layer equations for flows along curved surfaces was derived together with a modified eddy viscosity. Solution of these equations together with the modified eddy viscosity gave results that correlated well with the present data on flows along the convex surface with arbitrary pressure distribution. However, it could only be used to predict the mean characteristics of the flow along concave walls because of the existence of the system of longitudinal vortices inside the boundary layer.

A low dimensional dynamical system for the wall layer

NASA Technical Reports Server (NTRS)

Aubry, N.; Keefe, L. R.

1987-01-01

Low dimensional dynamical systems which model a fully developed turbulent wall layer were derived.The model is based on the optimally fast convergent proper orthogonal decomposition, or Karhunen-Loeve expansion. This decomposition provides a set of eigenfunctions which are derived from the autocorrelation tensor at zero time lag. Via Galerkin projection, low dimensional sets of ordinary differential equations in time, for the coefficients of the expansion, were derived from the Navier-Stokes equations. The energy loss to the unresolved modes was modeled by an eddy viscosity representation, analogous to Heisenberg's spectral model. A set of eigenfunctions and eigenvalues were obtained from direct numerical simulation of a plane channel at a Reynolds number of 6600, based on the mean centerline velocity and the channel width flow and compared with previous work done by Herzog. Using the new eigenvalues and eigenfunctions, a new ten dimensional set of ordinary differential equations were derived using five non-zero cross-stream Fourier modes with a periodic length of 377 wall units. The dynamical system was integrated for a range of the eddy viscosity prameter alpha. This work is encouraging.

Subgrid-scale models for large-eddy simulation of rotating turbulent channel flows

NASA Astrophysics Data System (ADS)

Silvis, Maurits H.; Bae, Hyunji Jane; Trias, F. Xavier; Abkar, Mahdi; Moin, Parviz; Verstappen, Roel

2017-11-01

We aim to design subgrid-scale models for large-eddy simulation of rotating turbulent flows. Rotating turbulent flows form a challenging test case for large-eddy simulation due to the presence of the Coriolis force. The Coriolis force conserves the total kinetic energy while transporting it from small to large scales of motion, leading to the formation of large-scale anisotropic flow structures. The Coriolis force may also cause partial flow laminarization and the occurrence of turbulent bursts. Many subgrid-scale models for large-eddy simulation are, however, primarily designed to parametrize the dissipative nature of turbulent flows, ignoring the specific characteristics of transport processes. We, therefore, propose a new subgrid-scale model that, in addition to the usual dissipative eddy viscosity term, contains a nondissipative nonlinear model term designed to capture transport processes, such as those due to rotation. We show that the addition of this nonlinear model term leads to improved predictions of the energy spectra of rotating homogeneous isotropic turbulence as well as of the Reynolds stress anisotropy in spanwise-rotating plane-channel flows. This work is financed by the Netherlands Organisation for Scientific Research (NWO) under Project Number 613.001.212.

Eddy current testing probe with dual half-cylindrical coils

NASA Astrophysics Data System (ADS)

Bae, Byung-Hoon; Choi, Jung-Mi; Kim, Soo-Yong

2000-02-01

We have developed a new eddy current probe composed of a dual half-cylindrical (2HC) coil as an exciting coil and a sensing coil that is placed in the small gap of the 2HC coil. The 2HC coil induces a linear eddy current on the narrow region within the target medium. The magnitude of eddy current has a maximum peak with the narrow width, underneath the center of the exciting 2HC coil. Because of the linear eddy current, the probe can be used to detect not only the existence of a crack but also its direction in conducting materials. Using specimen with a machined crack, and varying the exciting frequency from 0.5 to 100 kHz, we investigated the relationships between the direction of crack and the output voltage of the sensing coil.

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

NASA Technical Reports Server (NTRS)

Wang, Chi R.

2005-01-01

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

A two-equation model for heat transport in wall turbulent shear flows

NASA Astrophysics Data System (ADS)

Nagano, Y.; Kim, C.

1988-08-01

A new proposal for closing the energy equation is presented at the two-equation level of turbulence modeling. The eddy diffusivity concept is used in modeling. However, just as the eddy viscosity is determined from solutions of the k and epsilon equations, so the eddy diffusivity for heat is given as functions of temperature variance, and the dissipation rate of temperature fluctuations, together with k and epsilon. Thus, the proposed model does not require any questionable assumptions for the 'turbulent Prandtl number'. Modeled forms of the equations are developed to account for the physical effects of molecular Prandtl number and near-wall turbulence. The model is tested by application to a flat-plate boundary layer, the thermal entrance region of a pipe, and the turbulent heat transfer in fluids over a wide range of the Prandtl number. Agreement with the experiment is generally very satisfactory.

LES FOR SIMULATING THE GAS EXCHANGE PROCESS IN A SPARK IGNITION ENGINE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ameen, Muhsin M; yang, xiaofeng; kuo, tang-wei

2015-01-01

The gas exchange process is known to be a significant source of cyclic variability in Internal Combustion Engines (ICE). Traditionally, Large Eddy Simulations (LES) are expected to capture these cycle-to-cycle variations. This paper reports a numerical effort to establish best practices for capturing cyclic variability with LES tools in a Transparent Combustion Chamber (TCC) spark ignition engine. The main intention is to examine the sensitivity of cycle averaged mean and Root Mean Square (RMS) flow fields and Proper Orthogonal Decomposition (POD) modes to different computational hardware, adaptive mesh refinement (AMR) and LES sub-grid scale (SGS) models, since these aspects havemore » received little attention in the past couple of decades. This study also examines the effect of near-wall resolution on the predicted wall shear stresses. LES is pursued with commercially available CONVERGE code. Two different SGS models are tested, a one-equation eddy viscosity model and dynamic structure model. The results seem to indicate that both mean and RMS fields without any SGS model are not much different than those with LES models, either one-equation eddy viscosity or dynamic structure model. Computational hardware results in subtle quantitative differences, especially in RMS distributions. The influence of AMR on both mean and RMS fields is negligible. The predicted shear stresses near the liner walls is also found to be relatively insensitive to near-wall resolution except in the valve curtain region.« less

Drag reduction by a linear viscosity profile.

PubMed

De Angelis, Elisabetta; Casciola, Carlo M; L'vov, Victor S; Pomyalov, Anna; Procaccia, Itamar; Tiberkevich, Vasil

2004-11-01

Drag reduction by polymers in turbulent flows raises an apparent contradiction: the stretching of the polymers must increase the viscosity, so why is the drag reduced? A recent theory proposed that drag reduction, in agreement with experiments, is consistent with the effective viscosity growing linearly with the distance from the wall. With this self-consistent solution the reduction in the Reynolds stress overwhelms the increase in viscous drag. In this Rapid Communication we show, using direct numerical simulations, that a linear viscosity profile indeed reduces the drag in agreement with the theory and in close correspondence with direct simulations of the FENE-P model at the same flow conditions.

Numerical solutions of the complete Navier-Strokes equations. no. 27

NASA Technical Reports Server (NTRS)

Hassan, H. A.

1996-01-01

This report describes the development of an enstrophy model capable of predicting turbulence separation and its application to two airfoils at various angles of attack and Mach numbers. In addition, a two equation kappa-xi model with a tensor eddy viscosity was developed. Plans call for this model to be used in calculating three dimensional turbulent flows.

Application of a range of turbulence energy models to the determination of M4 tidal current profiles

NASA Astrophysics Data System (ADS)

Xing, Jiuxing; Davies, Alan M.

1996-04-01

A fully nonlinear, three-dimensional hydrodynamic model of the Irish Sea, using a range of turbulence energy sub-models, is used to examine the influence of the turbulence closure method upon the vertical variation of the current profile of the fundamental and higher harmonics of the tide in the region. Computed tidal current profiles are compared with previous calculations using a spectral model with eddy viscosity related to the flow field. The model has a sufficiently fine grid to resolve the advection terms, in particular the advection of turbulence and momentum. Calculations show that the advection of turbulence energy does not have a significant influence upon the current profile of either the fundamental or higher harmonic of the tide, although the advection of momentum is important in the region of headlands. The simplification of the advective terms by only including them in their vertically integrated form does not appear to make a significant difference to current profiles, but does reduce the computational effort by a significant amount. Computed current profiles both for the fundamental and the higher harmonic determined with a prognostic equation for turbulence and an algebraic mixing length formula, are as accurate as those determined with a two prognostic equation model (the so called q2- q2l model), provided the mixing length is specified correctly. A simple, flow-dependent eddy viscosity with a parabolic variation of viscosity also performs equally well.

Viscous properties of isotropic fluids composed of linear molecules: departure from the classical Navier-Stokes theory in nano-confined geometries.

PubMed

Hansen, J S; Daivis, Peter J; Todd, B D

2009-10-01

In this paper we present equilibrium molecular-dynamics results for the shear, rotational, and spin viscosities for fluids composed of linear molecules. The density dependence of the shear viscosity follows a stretched exponential function, whereas the rotational viscosity and the spin viscosities show approximately power-law dependencies. The frequency-dependent shear and spin viscosities are also studied. It is found that viscoelastic behavior is first manifested in the shear viscosity and that the real part of the spin viscosities features a maximum for nonzero frequency. The calculated transport coefficients are used together with the extended Navier-Stokes equations to investigate the effect of the coupling between the intrinsic angular momentum and linear momentum for highly confined fluids. Both steady and oscillatory flows are studied. It is shown, for example, that the fluid flow rate for Poiseuille flow is reduced by up to 10% in a 2 nm channel for a buta-triene fluid at density 236 kg m(-3) and temperature 306 K. The coupling effect may, therefore, become very important for nanofluidic applications.

Computational and Experimental Flow Field Analyses of Separate Flow Chevron Nozzles and Pylon Interaction

NASA Technical Reports Server (NTRS)

Massey, Steven J.; Thomas, Russell H.; AbdolHamid, Khaled S.; Elmiligui, Alaa A.

2003-01-01

A computational and experimental flow field analyses of separate flow chevron nozzles is presented. The goal of this study is to identify important flow physics and modeling issues required to provide highly accurate flow field data which will later serve as input to the Jet3D acoustic prediction code. Four configurations are considered: a baseline round nozzle with and without a pylon, and a chevron core nozzle with and without a pylon. The flow is simulated by solving the asymptotically steady, compressible, Reynolds-averaged Navier-Stokes equations using an implicit, up-wind, flux-difference splitting finite volume scheme and standard two-equation kappa-epsilon turbulence model with a linear stress representation and the addition of a eddy viscosity dependence on total temperature gradient normalized by local turbulence length scale. The current CFD results are seen to be in excellent agreement with Jet Noise Lab data and show great improvement over previous computations which did not compensate for enhanced mixing due to high temperature gradients.

Predicting NonInertial Effects with Algebraic Stress Models which Account for Dissipation Rate Anisotropies

NASA Technical Reports Server (NTRS)

Jongen, T.; Machiels, L.; Gatski, T. B.

1997-01-01

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coriolis-modified eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the turbulent model validation. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.

On the calculation of the complex wavenumber of plane waves in rigid-walled low-Mach-number turbulent pipe flows

NASA Astrophysics Data System (ADS)

Weng, Chenyang; Boij, Susann; Hanifi, Ardeshir

2015-10-01

A numerical method for calculating the wavenumbers of axisymmetric plane waves in rigid-walled low-Mach-number turbulent flows is proposed, which is based on solving the linearized Navier-Stokes equations with an eddy-viscosity model. In addition, theoretical models for the wavenumbers are reviewed, and the main effects (the viscothermal effects, the mean flow convection and refraction effects, the turbulent absorption, and the moderate compressibility effects) which may influence the sound propagation are discussed. Compared to the theoretical models, the proposed numerical method has the advantage of potentially including more effects in the computed wavenumbers. The numerical results of the wavenumbers are compared with the reviewed theoretical models, as well as experimental data from the literature. It shows that the proposed numerical method can give satisfactory prediction of both the real part (phase shift) and the imaginary part (attenuation) of the measured wavenumbers, especially when the refraction effects or the turbulent absorption effects become important.

One-equation near-wall turbulence modeling with the aid of direct simulation data

NASA Technical Reports Server (NTRS)

Rodi, W.; Mansour, N. N.

1990-01-01

The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is also examined in which bar-(v'(exp 2))(exp 1/2) is used as velocity scale instead of k(exp 1/2). With this velocity scale, the length scales now appearing in the model follow very closely a linear relationship near the wall so that no damping is necessary. For the determination of bar-v'(exp 2) in the context of a one-equation model, a correlation is provided between bar-(v'(exp 2))/k and bar-(u'v')/k.

Modelling the transitional boundary layer

NASA Technical Reports Server (NTRS)

Narasimha, R.

1990-01-01

Recent developments in the modelling of the transition zone in the boundary layer are reviewed (the zone being defined as extending from the station where intermittency begins to depart from zero to that where it is nearly unity). The value of using a new non-dimensional spot formation rate parameter, and the importance of allowing for so-called subtransitions within the transition zone, are both stressed. Models do reasonably well in constant pressure 2-dimensional flows, but in the presence of strong pressure gradients further improvements are needed. The linear combination approach works surprisingly well in most cases, but would not be so successful in situations where a purely laminar boundary layer would separate but a transitional one would not. Intermittency-weighted eddy viscosity methods do not predict peak surface parameters well without the introduction of an overshooting transition function whose connection with the spot theory of transition is obscure. Suggestions are made for further work that now appears necessary for developing improved models of the transition zone.

Energy loss due to eddy current in linear transformer driver cores

NASA Astrophysics Data System (ADS)

Kim, A. A.; Mazarakis, M. G.; Manylov, V. I.; Vizir, V. A.; Stygar, W. A.

2010-07-01

In linear transformer drivers [Phys. Rev. ST Accel. Beams 12, 050402 (2009)PRABFM1098-440210.1103/PhysRevSTAB.12.050402; Phys. Rev. ST Accel. Beams 12, 050401 (2009)PRABFM1098-440210.1103/PhysRevSTAB.12.050401] as well as any other linear induction accelerator cavities, ferromagnetic cores are used to prevent the current from flowing along the induction cavity walls which are in parallel with the load. But if the core is made of conductive material, the applied voltage pulse generates the eddy current in the core itself which heats the core and therefore also reduces the overall linear transformer driver (LTD) efficiency. The energy loss due to generation of the eddy current in the cores depends on the specific resistivity of the core material, the design of the core, as well as on the distribution of the eddy current in the core tape during the remagnetizing process. In this paper we investigate how the eddy current is distributed in a core tape with an arbitrary shape hysteresis loop. Our model is based on the textbook knowledge related to the eddy current generation in ferromagnetics with rectangular hysteresis loop, and in usual conductors. For the reader’s convenience, we reproduce some most important details of this knowledge in our paper. The model predicts that the same core would behave differently depending on how fast the applied voltage pulse is: in the high frequency limit, the equivalent resistance of the core reduces during the pulse whereas in the low frequency limit it is constant. An important inference is that the energy loss due to the eddy current generation can be reduced by increasing the cross section of the core over the minimum value which is required to avoid its saturation. The conclusions of the model are confirmed with experimental observations presented at the end of the paper.

Effect of Turbulence Modeling on Hovering Rotor Flows

NASA Technical Reports Server (NTRS)

Yoon, Seokkwan; Chaderjian, Neal M.; Pulliam, Thomas H.; Holst, Terry L.

2015-01-01

The effect of turbulence models in the off-body grids on the accuracy of solutions for rotor flows in hover has been investigated. Results from the Reynolds-Averaged Navier-Stokes and Laminar Off-Body models are compared. Advection of turbulent eddy viscosity has been studied to find the mechanism leading to inaccurate solutions. A coaxial rotor result is also included.

A computer program for the calculation of laminar and turbulent boundary layer flows

NASA Technical Reports Server (NTRS)

Dwyer, H. A.; Doss, E. D.; Goldman, A. L.

1972-01-01

The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.

A Galilean and tensorial invariant k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1993-01-01

A k-epsilon model is proposed for wall bounded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation rate equation is reformulated using this time scale and no singularity exists at the wall. A new parameter R = k/S(nu) is introduced to characterize the damping function in the eddy viscosity. This parameter is determined by local properties of both the mean and the turbulent flow fields and is free from any geometry parameter. The proposed model is then Galilean and tensorial invariant. The model constants used are the same as in the high Reynolds number Standard k-epsilon Model. Thus, the proposed model will also be suitable for flows far from the wall. Turbulent channel flows and turbulent boundary layer flows with and without pressure gradients are calculated. Comparisons with the data from direct numerical simulations and experiments show that the model predictions are excellent for turbulent channel flows and turbulent boundary layers with favorable pressure gradients, good for turbulent boundary layers with zero pressure gradients, and fair for turbulent boundary layer with adverse pressure gradients.

Influence of magnet eddy current on magnetization characteristics of variable flux memory machine

NASA Astrophysics Data System (ADS)

Yang, Hui; Lin, Heyun; Zhu, Z. Q.; Lyu, Shukang

2018-05-01

In this paper, the magnet eddy current characteristics of a newly developed variable flux memory machine (VFMM) is investigated. Firstly, the machine structure, non-linear hysteresis characteristics and eddy current modeling of low coercive force magnet are described, respectively. Besides, the PM eddy current behaviors when applying the demagnetizing current pulses are unveiled and investigated. The mismatch of the required demagnetization currents between the cases with or without considering the magnet eddy current is identified. In addition, the influences of the magnet eddy current on the demagnetization effect of VFMM are analyzed. Finally, a prototype is manufactured and tested to verify the theoretical analyses.

Early stages of transition in viscosity-stratified channel flow

NASA Astrophysics Data System (ADS)

Govindarajan, Rama; Jose, Sharath; Brandt, Luca

2013-11-01

In parallel shear flows, it is well known that transition to turbulence usually occurs through a subcritical process. In this work we consider a flow through a channel across which there is a linear temperature variation. The temperature gradient leads to a viscosity variation across the channel. A large body of work has been done in the linear regime for this problem, and it has been seen that viscosity stratification can lead to considerable changes in stability and transient growth characteristics. Moreover contradictory effects of introducing a non uniform viscosity in the system have been reported. We conduct a linear stability analysis and direct numerical simulations (DNS) for this system. We show that the optimal initial structures in the viscosity-stratified case, unlike in unstratified flow, do not span the width of the channel, but are focussed near one wall. The nonlinear consequences of the localisation of the structures will be discussed.

Computer program for solving laminar, transitional, or turbulent compressible boundary-layer equations for two-dimensional and axisymmetric flow

NASA Technical Reports Server (NTRS)

Harris, J. E.; Blanchard, D. K.

1982-01-01

A numerical algorithm and computer program are presented for solving the laminar, transitional, or turbulent two dimensional or axisymmetric compressible boundary-layer equations for perfect-gas flows. The governing equations are solved by an iterative three-point implicit finite-difference procedure. The software, program VGBLP, is a modification of the approach presented in NASA TR R-368 and NASA TM X-2458, respectively. The major modifications are: (1) replacement of the fourth-order Runge-Kutta integration technique with a finite-difference procedure for numerically solving the equations required to initiate the parabolic marching procedure; (2) introduction of the Blottner variable-grid scheme; (3) implementation of an iteration scheme allowing the coupled system of equations to be converged to a specified accuracy level; and (4) inclusion of an iteration scheme for variable-entropy calculations. These modifications to the approach presented in NASA TR R-368 and NASA TM X-2458 yield a software package with high computational efficiency and flexibility. Turbulence-closure options include either two-layer eddy-viscosity or mixing-length models. Eddy conductivity is modeled as a function of eddy viscosity through a static turbulent Prandtl number formulation. Several options are provided for specifying the static turbulent Prandtl number. The transitional boundary layer is treated through a streamwise intermittency function which modifies the turbulence-closure model. This model is based on the probability distribution of turbulent spots and ranges from zero to unity for laminar and turbulent flow, respectively. Several test cases are presented as guides for potential users of the software.

A High Resolution Study of Black Sea Circulation and Hypothetical Oil Spills

NASA Astrophysics Data System (ADS)

Dietrich, D. E.; Bowman, M. J.; Korotenko, K. A.

2008-12-01

A 1/24 deg resolution adaptation of the DieCAST ocean model simulates a realistically intense Rim Current and ubiquitous mesoscale coastal anticyclonic eddies that result from anticyclonic vorticity generation by laterally differential bottom drag forces that are amplified near Black Sea coastal headlands. Climatological and synoptic surface forcings are compared. The effects of vertical momentum transfer by known (by Synop region fishermen, as reported by Ballard National Geographic article) big amplitude internal waves are parameterized by big vertical viscosity. Sensitivity to vertical viscosity is shown. Results of simulated hypothetical oil spills are shown. A simple method to nowcast/forecast the Black Sea currents is described and early results are shown.

A Dynamic Eddy Viscosity Model for the Shallow Water Equations Solved by Spectral Element and Discontinuous Galerkin Methods

NASA Astrophysics Data System (ADS)

Marras, Simone; Suckale, Jenny; Giraldo, Francis X.; Constantinescu, Emil

2016-04-01

We present the solution of the viscous shallow water equations where viscosity is built as a residual-based subgrid scale model originally designed for large eddy simulation of compressible [1] and stratified flows [2]. The necessity of viscosity for a shallow water model not only finds motivation from mathematical analysis [3], but is supported by physical reasoning as can be seen by an analysis of the energetics of the solution. We simulated the flow of an idealized wave as it hits a set of obstacles. The kinetic energy spectrum of this flow shows that, although the inviscid Galerkin solutions -by spectral elements and discontinuous Galerkin [4]- preserve numerical stability in spite of the spurious oscillations in the proximity of the wave fronts, the slope of the energy cascade deviates from the theoretically expected values. We show that only a sufficiently small amount of dynamically adaptive viscosity removes the unwanted high-frequency modes while preserving the overall sharpness of the solution. In addition, it yields a physically plausible energy decay. This work is motivated by a larger interest in the application of a shallow water model to the solution of tsunami triggered coastal flows. In particular, coastal flows in regions around the world where coastal parks made of mitigation hills of different sizes and configurations are considered as a means to deviate the power of the incoming wave. References [1] M. Nazarov and J. Hoffman (2013) "Residual-based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods" Int. J. Numer. Methods Fluids, 71:339-357 [2] S. Marras, M. Nazarov, F. X. Giraldo (2015) "Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES" J. Comput. Phys. 301:77-101 [3] J. F. Gerbeau and B. Perthame (2001) "Derivation of the viscous Saint-Venant system for laminar shallow water; numerical validation" Discrete Contin. Dyn. Syst. Ser. B, 1:89?102 [4] F. X. Giraldo and M. Restelli (2010) "High-order semi-implicit time-integrators for a triangular discontinuous Galerkin oceanic shallow water model. Int. J. Numer. Methods Fluids, 63:1077-1102

Downscaling modelling system for multi-scale air quality forecasting

NASA Astrophysics Data System (ADS)

Nuterman, R.; Baklanov, A.; Mahura, A.; Amstrup, B.; Weismann, J.

2010-09-01

Urban modelling for real meteorological situations, in general, considers only a small part of the urban area in a micro-meteorological model, and urban heterogeneities outside a modelling domain affect micro-scale processes. Therefore, it is important to build a chain of models of different scales with nesting of higher resolution models into larger scale lower resolution models. Usually, the up-scaled city- or meso-scale models consider parameterisations of urban effects or statistical descriptions of the urban morphology, whereas the micro-scale (street canyon) models are obstacle-resolved and they consider a detailed geometry of the buildings and the urban canopy. The developed system consists of the meso-, urban- and street-scale models. First, it is the Numerical Weather Prediction (HIgh Resolution Limited Area Model) model combined with Atmospheric Chemistry Transport (the Comprehensive Air quality Model with extensions) model. Several levels of urban parameterisation are considered. They are chosen depending on selected scales and resolutions. For regional scale, the urban parameterisation is based on the roughness and flux corrections approach; for urban scale - building effects parameterisation. Modern methods of computational fluid dynamics allow solving environmental problems connected with atmospheric transport of pollutants within urban canopy in a presence of penetrable (vegetation) and impenetrable (buildings) obstacles. For local- and micro-scales nesting the Micro-scale Model for Urban Environment is applied. This is a comprehensive obstacle-resolved urban wind-flow and dispersion model based on the Reynolds averaged Navier-Stokes approach and several turbulent closures, i.e. k -ɛ linear eddy-viscosity model, k - ɛ non-linear eddy-viscosity model and Reynolds stress model. Boundary and initial conditions for the micro-scale model are used from the up-scaled models with corresponding interpolation conserving the mass. For the boundaries a kind of Dirichlet condition is chosen to provide the values based on interpolation from the coarse to the fine grid. When the roughness approach is changed to the obstacle-resolved one in the nested model, the interpolation procedure will increase the computational time (due to additional iterations) for meteorological/ chemical fields inside the urban sub-layer. In such situations, as a possible alternative, the perturbation approach can be applied. Here, the effects of main meteorological variables and chemical species are considered as a sum of two components: background (large-scale) values, described by the coarse-resolution model, and perturbations (micro-scale) features, obtained from the nested fine resolution model.

Linearized simulation of flow over wind farms and complex terrains.

PubMed

Segalini, Antonio

2017-04-13

The flow over complex terrains and wind farms is estimated here by numerically solving the linearized Navier-Stokes equations. The equations are linearized around the unperturbed incoming wind profile, here assumed logarithmic. The Boussinesq approximation is used to model the Reynolds stress with a prescribed turbulent eddy viscosity profile. Without requiring the boundary-layer approximation, two new linear equations are obtained for the vertical velocity and the wall-normal vorticity, with a reduction in the computational cost by a factor of 8 when compared with a primitive-variables formulation. The presence of terrain elevation is introduced as a vertical coordinate shift, while forestry or wind turbines are included as body forces, without any assumption about the wake structure for the turbines. The model is first validated against some available experiments and simulations, and then a simulation of a wind farm over a Gaussian hill is performed. The speed-up effect of the hill is clearly beneficial in terms of the available momentum upstream of the crest, while downstream of it the opposite can be said as the turbines face a decreased wind speed. Also, the presence of the hill introduces an additional spanwise velocity component that may also affect the turbines' operations. The linear superposition of the flow over the hill and the flow over the farm alone provided a first estimation of the wind speed along the farm, with discrepancies of the same order of magnitude for the spanwise velocity. Finally, the possibility of using a parabolic set of equations to obtain the turbulent kinetic energy after the linearized model is investigated with promising results.This article is part of the themed issue 'Wind energy in complex terrains'. © 2017 The Author(s).

Linearized simulation of flow over wind farms and complex terrains

NASA Astrophysics Data System (ADS)

Segalini, Antonio

2017-03-01

The flow over complex terrains and wind farms is estimated here by numerically solving the linearized Navier-Stokes equations. The equations are linearized around the unperturbed incoming wind profile, here assumed logarithmic. The Boussinesq approximation is used to model the Reynolds stress with a prescribed turbulent eddy viscosity profile. Without requiring the boundary-layer approximation, two new linear equations are obtained for the vertical velocity and the wall-normal vorticity, with a reduction in the computational cost by a factor of 8 when compared with a primitive-variables formulation. The presence of terrain elevation is introduced as a vertical coordinate shift, while forestry or wind turbines are included as body forces, without any assumption about the wake structure for the turbines. The model is first validated against some available experiments and simulations, and then a simulation of a wind farm over a Gaussian hill is performed. The speed-up effect of the hill is clearly beneficial in terms of the available momentum upstream of the crest, while downstream of it the opposite can be said as the turbines face a decreased wind speed. Also, the presence of the hill introduces an additional spanwise velocity component that may also affect the turbines' operations. The linear superposition of the flow over the hill and the flow over the farm alone provided a first estimation of the wind speed along the farm, with discrepancies of the same order of magnitude for the spanwise velocity. Finally, the possibility of using a parabolic set of equations to obtain the turbulent kinetic energy after the linearized model is investigated with promising results. This article is part of the themed issue 'Wind energy in complex terrains'.

A new class of finite element variational multiscale turbulence models for incompressible magnetohydrodynamics

DOE PAGES

Sondak, D.; Shadid, J. N.; Oberai, A. A.; ...

2015-04-29

New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and highmore » Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.« less

Large-Eddy Simulation of Internal Flow through Human Vocal Folds

NASA Astrophysics Data System (ADS)

Lasota, Martin; Šidlof, Petr

2018-06-01

The phonatory process occurs when air is expelled from the lungs through the glottis and the pressure drop causes flow-induced oscillations of the vocal folds. The flow fields created in phonation are highly unsteady and the coherent vortex structures are also generated. For accuracy it is essential to compute on humanlike computational domain and appropriate mathematical model. The work deals with numerical simulation of air flow within the space between plicae vocales and plicae vestibulares. In addition to the dynamic width of the rima glottidis, where the sound is generated, there are lateral ventriculus laryngis and sacculus laryngis included in the computational domain as well. The paper presents the results from OpenFOAM which are obtained with a large-eddy simulation using second-order finite volume discretization of incompressible Navier-Stokes equations. Large-eddy simulations with different subgrid scale models are executed on structured mesh. In these cases are used only the subgrid scale models which model turbulence via turbulent viscosity and Boussinesq approximation in subglottal and supraglottal area in larynx.

A family of dynamic models for large-eddy simulation

NASA Technical Reports Server (NTRS)

Carati, D.; Jansen, K.; Lund, T.

1995-01-01

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Finite-mode spectral model of homogeneous and isotropic Navier-stokes turbulence: a rapidly depleted energy cascade.

PubMed

Lévêque, E; Koudella, C R

2001-04-30

An eddy-viscous term is added to Navier-Stokes dynamics at wave numbers k greater than the inflection point kc of the energy flux F(log(k)). The eddy viscosity is fixed so that the energy spectrum satisfies E(k) = E(kc) (k/kc)(-3) for k>kc. This resulting forcing induces a rapid depletion of the energy cascade at k>kc. It is observed numerically that the model reproduces turbulence energetics at k< or =kc and statistics of two-point velocity correlations at scales r>lambda (Taylor microscale). Compared to a direct numerical simulation of R(lambda) = 130 an equivalent run with the present model results in a gain of a factor 20 in CPU time.

A comparison of three approaches to compute the effective Reynolds number of the implicit large-eddy simulations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhou, Ye; Thornber, Ben

2016-04-12

Here, the implicit large-eddy simulation (ILES) has been utilized as an effective approach for calculating many complex flows at high Reynolds number flows. Richtmyer–Meshkov instability (RMI) induced flow can be viewed as a homogeneous decaying turbulence (HDT) after the passage of the shock. In this article, a critical evaluation of three methods for estimating the effective Reynolds number and the effective kinematic viscosity is undertaken utilizing high-resolution ILES data. Effective Reynolds numbers based on the vorticity and dissipation rate, or the integral and inner-viscous length scales, are found to be the most self-consistent when compared to the expected phenomenology andmore » wind tunnel experiments.« less

Image-based modelling of lateral magma flow: the Basement Sill, Antarctica

PubMed Central

Mirhadizadeh, Seyed

2017-01-01

The McMurdo Dry Valleys magmatic system, Antarctica, provides a world-class example of pervasive lateral magma flow on a continental scale. The lowermost intrusion (Basement Sill) offers detailed sections through the now frozen particle microstructure of a congested magma slurry. We simulated the flow regime in two and three dimensions using numerical models built on a finite-element mesh derived from field data. The model captures the flow behaviour of the Basement Sill magma over a viscosity range of 1–104 Pa s where the higher end (greater than or equal to 102 Pa s) corresponds to a magmatic slurry with crystal fractions varying between 30 and 70%. A novel feature of the model is the discovery of transient, low viscosity (less than or equal to 50 Pa s) high Reynolds number eddies formed along undulating contacts at the floor and roof of the intrusion. Numerical tracing of particle orbits implies crystals trapped in eddies segregate according to their mass density. Recovered shear strain rates (10−3–10−5 s−1) at viscosities equating to high particle concentrations (around more than 40%) in the Sill interior point to shear-thinning as an explanation for some types of magmatic layering there. Model transport rates for the Sill magmas imply a maximum emplacement time of ca 105 years, consistent with geochemical evidence for long-range lateral flow. It is a theoretically possibility that fast-flowing magma on a continental scale will be susceptible to planetary-scale rotational forces. PMID:28573002

Image-based modelling of lateral magma flow: the Basement Sill, Antarctica.

PubMed

Petford, Nick; Mirhadizadeh, Seyed

2017-05-01

The McMurdo Dry Valleys magmatic system, Antarctica, provides a world-class example of pervasive lateral magma flow on a continental scale. The lowermost intrusion (Basement Sill) offers detailed sections through the now frozen particle microstructure of a congested magma slurry. We simulated the flow regime in two and three dimensions using numerical models built on a finite-element mesh derived from field data. The model captures the flow behaviour of the Basement Sill magma over a viscosity range of 1-10 4  Pa s where the higher end (greater than or equal to 10 2  Pa s) corresponds to a magmatic slurry with crystal fractions varying between 30 and 70%. A novel feature of the model is the discovery of transient, low viscosity (less than or equal to 50 Pa s) high Reynolds number eddies formed along undulating contacts at the floor and roof of the intrusion. Numerical tracing of particle orbits implies crystals trapped in eddies segregate according to their mass density. Recovered shear strain rates (10 -3 -10 -5  s -1 ) at viscosities equating to high particle concentrations (around more than 40%) in the Sill interior point to shear-thinning as an explanation for some types of magmatic layering there. Model transport rates for the Sill magmas imply a maximum emplacement time of ca 10 5 years, consistent with geochemical evidence for long-range lateral flow. It is a theoretically possibility that fast-flowing magma on a continental scale will be susceptible to planetary-scale rotational forces.

Classical closure theory and Lam's interpretation of epsilon-RNG

NASA Technical Reports Server (NTRS)

Zhou, YE

1995-01-01

Lam's phenomenological epsilon-renormalization group (RNG) model is quite different from the other members of that group. It does not make use of the correspondence principle and the epsilon-expansion procedure. We demonstrate that Lam's epsilon-RNG model is essentially the physical space version of the classical closure theory in spectral space and consider the corresponding treatment of the eddy viscosity and energy backscatter.

Comparison of Laminar and Linear Eddy Model Closures for Combustion Instability Simulations

DTIC Science & Technology

2015-07-01

14. ABSTRACT Unstable liquid rocket engines can produce highly complex dynamic flowfields with features such as rapid changes in temperature and...applicability. In the present study, the linear eddy model (LEM) is applied to an unstable single element liquid rocket engine to assess its performance and to...Sankaran‡ Air Force Research Laboratory, Edwards AFB, CA, 93524 Unstable liquid rocket engines can produce highly complex dynamic flowfields with features

Stationary Waves of the Ice Age Climate.

NASA Astrophysics Data System (ADS)

Cook, Kerry H.; Held, Isaac M.

1988-08-01

A linearized, steady state, primitive equation model is used to simulate the climatological zonal asymmetries (stationary eddies) in the wind and temperature fields of the 18 000 YBP climate during winter. We compare these results with the eddies simulated in the ice age experiments of Broccoli and Manabe, who used CLIMAP boundary conditions and reduced atmospheric CO2 in an atmospheric general circulation model (GCM) coupled with a static mixed layer ocean model. The agreement between the models is good, indicating that the linear model can be used to evaluate the relative influences of orography, diabatic heating, and transient eddy heat and momentum transports in generating stationary waves. We find that orographic forcing dominates in the ice age climate. The mechanical influence of the continental ice sheets on the atmosphere is responsible for most of the changes between the present day and ice age stationary eddies. This concept of the ice age climate is complicated by the sensitivity of the stationary eddies to the large increase in the magnitude of the zonal mean meridional temperature gradient simulated in the ice age GCM.

Development of a One-Equation Eddy Viscosity Turbulence Model for Application to Complex Turbulent Flows

NASA Astrophysics Data System (ADS)

Wray, Timothy J.

Computational fluid dynamics (CFD) is routinely used in performance prediction and design of aircraft, turbomachinery, automobiles, and in many other industrial applications. Despite its wide range of use, deficiencies in its prediction accuracy still exist. One critical weakness is the accurate simulation of complex turbulent flows using the Reynolds-Averaged Navier-Stokes equations in conjunction with a turbulence model. The goal of this research has been to develop an eddy viscosity type turbulence model to increase the accuracy of flow simulations for mildly separated flows, flows with rotation and curvature effects, and flows with surface roughness. It is accomplished by developing a new zonal one-equation turbulence model which relies heavily on the flow physics; it is now known in the literature as the Wray-Agarwal one-equation turbulence model. The effectiveness of the new model is demonstrated by comparing its results with those obtained by the industry standard one-equation Spalart-Allmaras model and two-equation Shear-Stress-Transport k - o model and experimental data. Results for subsonic, transonic, and supersonic flows in and about complex geometries are presented. It is demonstrated that the Wray-Agarwal model can provide the industry and CFD researchers an accurate, efficient, and reliable turbulence model for the computation of a large class of complex turbulent flows.

An improved k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Shih, T. H.; Hsu, Andrew T.

1991-01-01

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. In the first part of this work, the near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation are analyzed. Based on these analyses, a modified eddy viscosity model with the correct near-wall behavior is suggested, and a model for the pressure transport term in the k-equation is proposed. In addition, a modeled dissipation rate equation is reformulated, and a boundary condition for the dissipation rate is suggested. In the second part of the work, one of the deficiencies of the existing k-epsilon models, namely, the wall distance dependency of the equations and the damping functions, is examined. An improved model that does not depend on any wall distance is introduced. Fully developed turbulent channel flows and turbulent boundary layers over a flat plate are studied as validations for the proposed new models. Numerical results obtained from the present and other previous k-epsilon models are compared with data from direct numerical simulation. The results show that the present k-epsilon model, with added robustness, performs as well as or better than other existing models in predicting the behavior of near-wall turbulence.

A comparison of capillary and rotational viscometry of aqueous solutions of hypromellose.

PubMed

Sklubalová, Z; Zatloukal, Z

2007-10-01

A comparison of capillary and rotational viscometry of gentle pseudoplastic solutions of hypromellose (HPMC 4000) by using only single-point value of viscosity is difficult. Single-point comparison becomes topical in consequence to the pharmacopoeial requirement that the apparent viscosity of 2% hypromellose solution should be read at the shear rate of approximately 10 s(-1). This communication is focused on the estimation of the suitable shear rate, D eta, at which the apparent viscosity read using the rotational viscometer is numerically equal to the dynamic viscosity read using a capillary viscometer. For the solutions of HPMC in concentrations up to 2% w/v, the non-linear regression equations generated showed the influencing of the D eta value by the dynamic viscosity and/or by the originally derived linear velocity of the solution flowing through the capillary viscometer tube. To compare the apparent viscosity read using the rotational viscometer with the dynamic viscosity read using capillary viscometer, the exact estimation of the shear rate D eta at which both viscosities are numerically equal is essential since it is markedly affected by the concentration of HPMC solution.

Notes on rotating turbulence

NASA Technical Reports Server (NTRS)

Zeman, Otto

1994-01-01

This work investigates the turbulent constitutive relation when turbulence is subjected to solid body rotation. Laws regarding spectra and asymptotic decay of rotating homogeneous turbulence were confirmed through large-eddy simulation (LES) computations. Rotating turbulent flows exist in many industrial, geophysical, and astrophysical applications. From Lagrangian analysis a relation between turbulent stress and strain in rotating homogeneous turbulence was inferred. This relation was used to derive the spectral energy flux and, ultimately, the energy spectrum form. If the rotation wavenumber k(sub Omega) lies in the inertial subrange, then for wavenumbers less than k(sub Omega) the turbulence motions are affected by rotation and the energy spectrum slope is modified. Energy decay laws inferred in other reports and the present results suggest a modification of the epsilon model equation and eddy viscosity in k-epsilon models.

Shock wave-droplet interaction

NASA Astrophysics Data System (ADS)

Habibi Khoshmehr, Hamed; Krechetnikov, Rouslan

2016-11-01

Disintegration of a liquid droplet under the action of a shock wave is experimentally investigated. The shock wave-pulse is electromagnetically generated by discharging a high voltage capacitor into a flat spiral coil, above which an isolated circular metal membrane is placed in a close proximity. The Lorentz force arising due to the eddy current induced in the membrane abruptly accelerates it away from the spiral coil thus generating a shock wave. The liquid droplet placed at the center of the membrane, where the maximum deflection occurs, is disintegrated in the process of interaction with the shock wave. The effects of droplet viscosity and surface tension on the droplet destruction are studied with high-speed photography. Water-glycerol solution at different concentrations is used for investigating the effect of viscosity and various concentrations of water-sugar and water-ethanol solution are used for studying the effect of surface tension. Here we report on how the metamorphoses, which a liquid drop undergoes in the process of interaction with a shock wave, are affected by varied viscosity and surface tension.

Linear simulation of the stationary eddies in a GCM. II - The 'Mountain' model

NASA Technical Reports Server (NTRS)

Nigam, Sumant; Held, Isaac M.; Lyons, Steven W.

1988-01-01

Linear stationary wave theory is used to account for zonal asymmetries of the winter-averaged tropospheric circulation obtained in a GCM. The eddy zonal velocity field in the upper troposphere indicates that the orographic and thermal plus transient contributions are nearly equal in amplitude, while the eddy meridional velocity field (which is dominated by shorter zonal scales) shows the orographic contribution to be dominant. The two contributions are found to be roughly in phase over the east Asian coast, and they contribute roughly equal amounts to the low level Siberian high. Results indicate that the 300 mb extratropical response to tropical forcing reaches 50 gpm over Alaska, and that the responses to sensible heating and lower tropospheric transients are strongly anticorrelated.

Impacts of mesoscale eddies in the South China Sea on biogeochemical cycles

NASA Astrophysics Data System (ADS)

Guo, Mingxian; Chai, Fei; Xiu, Peng; Li, Shiyu; Rao, Shivanesh

2015-09-01

Biogeochemical cycles associated with mesoscale eddies in the South China Sea (SCS) were investigated. The study was based on a coupled physical-biogeochemical Pacific Ocean model (Regional Ocean Model System-Carbon, Silicate, and Nitrogen Ecosystem, ROMS-CoSiNE) simulation for the period from 1991 to 2008. A total of 568 mesoscale eddies with lifetime longer than 30 days were used in the analysis. Composite analysis revealed that the cyclonic eddies were associated with abundance of nutrients, phytoplankton, and zooplankton while the anticyclonic eddies depressed biogeochemical cycles, which are generally controlled by the eddy pumping mechanism. In addition, diatoms were dominant in phytoplankton species due to the abundance of silicate. Dipole structures of vertical fluxes with net upward motion in cyclonic eddies and net downward motion in anticyclonic eddies were revealed. During the lifetime of an eddy, the evolutions of physical, biological, and chemical structures were not linearly coupled at the eddy core where plankton grew, and composition of the community depended not only on the physical and chemical processes but also on the adjustments by the predator-prey relationship.

Retrieval of Urban Boundary Layer Structures from Doppler Lidar Data. Part I: Accuracy Assessment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xia, Quanxin; Lin, Ching Long; Calhoun, Ron

2008-01-01

Two coherent Doppler lidars from the US Army Research Laboratory (ARL) and Arizona State University (ASU) were deployed in the Joint Urban 2003 atmospheric dispersion field experiment (JU2003) held in Oklahoma City. The dual lidar data are used to evaluate the accuracy of the four-dimensional variational data assimilation (4DVAR) method and identify the coherent flow structures in the urban boundary layer. The objectives of the study are three-fold. The first objective is to examine the effect of eddy viscosity models on the quality of retrieved velocity data. The second objective is to determine the fidelity of single-lidar 4DVAR and evaluatemore » the difference between single- and dual-lidar retrievals. The third objective is to correlate the retrieved flow structures with the ground building data. It is found that the approach of treating eddy viscosity as part of control variables yields better results than the approach of prescribing viscosity. The ARL single-lidar 4DVAR is able to retrieve radial velocity fields with an accuracy of 98% in the along-beam direction and 80-90% in the cross-beam direction. For the dual-lidar 4DVAR, the accuracy of retrieved radial velocity in the ARL cross-beam direction improves to 90-94%. By using the dual-lidar retrieved data as a reference, the single-lidar 4DVAR is able to recover fluctuating velocity fields with 70-80% accuracy in the along-beam direction and 60-70% accuracy in the cross-beam direction. Large-scale convective roll structures are found in the vicinity of downtown airpark and parks. Vortical structures are identified near the business district. Strong updrafts and downdrafts are also found above a cluster of restaurants.« less

Effective Inflow Conditions for Turbulence Models in Aerodynamic Calculations

NASA Technical Reports Server (NTRS)

Spalart, Philippe R.; Rumsey, Christopher L.

2007-01-01

The selection of inflow values at boundaries far upstream of an aircraft is considered, for one- and two-equation turbulence models. Inflow values are distinguished from the ambient values near the aircraft, which may be much smaller. Ambient values should be selected first, and inflow values that will lead to them after the decay second; this is not always possible, especially for the time scale. The two-equation decay during the approach to the aircraft is shown; often, the time scale has been set too short for this decay to be calculated accurately on typical grids. A simple remedy for both issues is to impose floor values for the turbulence variables, outside the viscous sublayer, and it is argued that overriding the equations in this manner is physically justified. Selecting laminar ambient values is easy, if the boundary layers are to be tripped, but a more common practice is to seek ambient values that will cause immediate transition in boundary layers. This opens up a wide range of values, and selection criteria are discussed. The turbulent Reynolds number, or ratio of eddy viscosity to laminar viscosity has a huge dynamic range that makes it unwieldy; it has been widely mis-used, particularly by codes that set upper limits on it. The value of turbulent kinetic energy in a wind tunnel or the atmosphere is also of dubious value as an input to the model. Concretely, the ambient eddy viscosity must be small enough to preserve potential cores in small geometry features, such as flap gaps. The ambient frequency scale should also be small enough, compared with shear rates in the boundary layer. Specific values are recommended and demonstrated for airfoil flows

Quantifying residual, eddy, and mean flow effects on mixing in an idealized circumpolar current

DOE PAGES

Wolfram, Phillip J.; Ringler, Todd D.

2017-07-13

Meridional diffusivity is assessed in this paper for a baroclinically unstable jet in a high-latitudeIdealized Circumpolar Current (ICC) using the Model for Prediction Across Scales-Ocean (MPAS-O) and the online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) diagnostic via space-time dispersion of particle clusters over 120 monthly realizations of O(10 6) particles on 11 potential density surfaces. Diffusivity in the jet reaches values of O(6000 m 2 s -1) and is largest near the critical layer supporting mixing suppression and critical layer theory. Values in the vicinity of the shelf break are suppressed to O(100 m 2 s -1) due tomore » the presence of westward slope front currents. Diffusivity attenuates less rapidly with depth in the jet than both eddy velocity and kinetic energy scalings would suggest. Removal of the mean flow via high-pass filtering shifts the nonlinear parameter (ratio of the eddy velocity to eddy phase speed) into the linear wave regime by increasing the eddy phase speed via the depth-mean flow. Low-pass filtering, in contrast, quantifies the effect of mean shear. Diffusivity is decomposed into mean flow shear, linear waves, and the residual nonhomogeneous turbulence components, where turbulence dominates and eddy-produced filamentation strained by background mean shear enhances mixing, accounting for ≥ 80% of the total diffusivity relative to mean shear [O(100 m 2 s -1)], linear waves [O(1000 m 2 s -1)], and undecomposed full diffusivity [O(6000 m 2 s -1)]. Finally, diffusivity parameterizations accounting for both the nonhomogeneous turbulence residual and depth variability are needed.« less

Eddy, drift wave and zonal flow dynamics in a linear magnetized plasma

PubMed Central

Arakawa, H.; Inagaki, S.; Sasaki, M.; Kosuga, Y.; Kobayashi, T.; Kasuya, N.; Nagashima, Y.; Yamada, T.; Lesur, M.; Fujisawa, A.; Itoh, K.; Itoh, S.-I.

2016-01-01

Turbulence and its structure formation are universal in neutral fluids and in plasmas. Turbulence annihilates global structures but can organize flows and eddies. The mutual-interactions between flow and the eddy give basic insights into the understanding of non-equilibrium and nonlinear interaction by turbulence. In fusion plasma, clarifying structure formation by Drift-wave turbulence, driven by density gradients in magnetized plasma, is an important issue. Here, a new mutual-interaction among eddy, drift wave and flow in magnetized plasma is discovered. A two-dimensional solitary eddy, which is a perturbation with circumnavigating motion localized radially and azimuthally, is transiently organized in a drift wave – zonal flow (azimuthally symmetric band-like shear flows) system. The excitation of the eddy is synchronized with zonal perturbation. The organization of the eddy has substantial impact on the acceleration of zonal flow. PMID:27628894

Boundary Layers for the Navier-Stokes Equations Linearized Around a Stationary Euler Flow

NASA Astrophysics Data System (ADS)

Gie, Gung-Min; Kelliher, James P.; Mazzucato, Anna L.

2018-03-01

We study the viscous boundary layer that forms at small viscosity near a rigid wall for the solution to the Navier-Stokes equations linearized around a smooth and stationary Euler flow (LNSE for short) in a smooth bounded domain Ω \\subset R^3 under no-slip boundary conditions. LNSE is supplemented with smooth initial data and smooth external forcing, assumed ill-prepared, that is, not compatible with the no-slip boundary condition. We construct an approximate solution to LNSE on the time interval [0, T], 0

Development of Interior Permanent Magnet Motors with Concentrated Windings for Reducing Magnet Eddy Current Loss

NASA Astrophysics Data System (ADS)

Yamazaki, Katsumi; Kanou, Yuji; Fukushima, Yu; Ohki, Shunji; Nezu, Akira; Ikemi, Takeshi; Mizokami, Ryoichi

In this paper, we present the development of interior magnet motors with concentrated windings, which reduce the eddy current loss of the magnets. First, the mechanism of the magnet eddy current loss generation is investigated by a simple linear magnetic circuit. Due to the consideration, an automatic optimization method using an adaptive finite element method is carried out to determine the stator and rotor shapes, which decrease the eddy current loss of the magnet. The determined stator and rotor are manufactured in order to proof the effectiveness by the measurement.

Options for refractive index and viscosity matching to study variable density flows

NASA Astrophysics Data System (ADS)

Clément, Simon A.; Guillemain, Anaïs; McCleney, Amy B.; Bardet, Philippe M.

2018-02-01

Variable density flows are often studied by mixing two miscible aqueous solutions of different densities. To perform optical diagnostics in such environments, the refractive index of the fluids must be matched, which can be achieved by carefully choosing the two solutes and the concentration of the solutions. To separate the effects of buoyancy forces and viscosity variations, it is desirable to match the viscosity of the two solutions in addition to their refractive index. In this manuscript, several pairs of index matched fluids are compared in terms of viscosity matching, monetary cost, and practical use. Two fluid pairs are studied in detail, with two aqueous solutions (binary solutions of water and a salt or alcohol) mixed into a ternary solution. In each case: an aqueous solution of isopropanol mixed with an aqueous solution of sodium chloride (NaCl) and an aqueous solution of glycerol mixed with an aqueous solution of sodium sulfate (Na_2SO_4). The first fluid pair allows reaching high-density differences at low cost, but brings a large difference in dynamic viscosity. The second allows matching dynamic viscosity and refractive index simultaneously, at reasonable cost. For each of these four solutes, the density, kinematic viscosity, and refractive index are measured versus concentration and temperature, as well as wavelength for the refractive index. To investigate non-linear effects when two index-matched, binary solutions are mixed, the ternary solutions formed are also analyzed. Results show that density and refractive index follow a linear variation with concentration. However, the viscosity of the isopropanol and NaCl pair deviates from the linear law and has to be considered. Empirical correlations and their coefficients are given to create index-matched fluids at a chosen temperature and wavelength. Finally, the effectiveness of the refractive index matching is illustrated with particle image velocimetry measurements performed for a buoyant jet in a linearly stratified environment. The creation of the index-matched solutions and linear stratification in a large-scale experimental facility are detailed, as well as the practical challenges to obtain precise refractive index matching.

Film cooling from inclined cylindrical holes using large eddy simulations

NASA Astrophysics Data System (ADS)

Peet, Yulia V.

2006-12-01

The goal of the present study is to investigate numerically the physics of the flow, which occurs during the film cooling from inclined cylindrical holes, Film cooling is a technique used in gas turbine industry to reduce heat fluxes to the turbine blade surface. Large Eddy Simulation (LES) is performed modeling a realistic film cooling configuration, which consists of a large stagnation-type reservoir, feeding an array of discrete cooling holes (film holes) flowing into a flat plate turbulent boundary layer. Special computational methodology is developed for this problem, involving coupled simulations using multiple computational codes. A fully compressible LES code is used in the area above the flat plate, while a low Mach number LES code is employed in the plenum and film holes. The motivation for using different codes comes from the essential difference in the nature of the flow in these different regions. Flowfield is analyzed inside the plenum, film hole and a crossflow region. Flow inside the plenum is stagnating, except for the region close to the exit, where it accelerates rapidly to turn into the hole. The sharp radius of turning at the trailing edge of the plenum pipe connection causes the flow to separate from the downstream wall of the film hole. After coolant injection occurs, a complex flowfield is formed consisting of coherent vortical structures responsible for bringing hot crossflow fluid in contact with the walls of either the film hole or the blade, thus reducing cooling protection. Mean velocity and turbulent statistics are compared to experimental measurements, yielding good agreement for the mean flowfield and satisfactory agreement for the turbulence quantities. LES results are used to assess the applicability of basic assumptions of conventional eddy viscosity turbulence models used with Reynolds-averaged (RANS) approach, namely the isotropy of an eddy viscosity and thermal diffusivity. It is shown here that these assumptions do not hold for the film cooling flows. Comparison of film cooling effectiveness with experiments shows fair agreement for the centerline and laterally-averaged effectiveness. Lateral growth of the jet as judged from the lateral distribution of effectiveness is predicted correctly.

Progress in simulating industrial flows using two-equation models: Can more be achieved with further research?

NASA Technical Reports Server (NTRS)

Haroutunian, Vahe

1995-01-01

This viewgraph presentation provides a brief review of two-equation eddy-viscosity models (TEM's) from the perspective of applied CFD. It provides objective assessment of both well-known and newer models, compares model predictions from various TEM's with experiments, identifies sources of modeling error and gives historical perspective of their effects on model performance and assessment, and recommends directions for future research on TEM's.

Calculation of Suspended Sediment Transport by Combined Wave-Current Flows.

DTIC Science & Technology

1994-11-01

Anderson, and Silberg (1985) presented a model that had an eddy viscosity and boundary layer thickness that varied in time. The reference concentration was...sediment model. This model, along with that of Fredsoe, Anderson, and Silberg (1985), are the only two models that account for both the sediment and the...ignores any correlation between the periodic components of the velocity and the concentration. Even in the model of Fredsoe, Anderson, and Silberg (1985

Application of a new K-tau model to near wall turbulent flows

NASA Technical Reports Server (NTRS)

Thangam, S.; Abid, R.; Speziale, Charles G.

1991-01-01

A recently developed K-tau model for near wall turbulent flows is applied to two severe test cases. The turbulent flows considered include the incompressible flat plate boundary layer with the adverse pressure gradients and incompressible flow past a backward facing step. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

Modeling of structural uncertainties in Reynolds-averaged Navier-Stokes closures

NASA Astrophysics Data System (ADS)

Emory, Michael; Larsson, Johan; Iaccarino, Gianluca

2013-11-01

Estimation of the uncertainty in numerical predictions by Reynolds-averaged Navier-Stokes closures is a vital step in building confidence in such predictions. An approach to model-form uncertainty quantification that does not assume the eddy-viscosity hypothesis to be exact is proposed. The methodology for estimation of uncertainty is demonstrated for plane channel flow, for a duct with secondary flows, and for the shock/boundary-layer interaction over a transonic bump.

Anisotropic diffusion in mesh-free numerical magnetohydrodynamics

NASA Astrophysics Data System (ADS)

Hopkins, Philip F.

2017-04-01

We extend recently developed mesh-free Lagrangian methods for numerical magnetohydrodynamics (MHD) to arbitrary anisotropic diffusion equations, including: passive scalar diffusion, Spitzer-Braginskii conduction and viscosity, cosmic ray diffusion/streaming, anisotropic radiation transport, non-ideal MHD (Ohmic resistivity, ambipolar diffusion, the Hall effect) and turbulent 'eddy diffusion'. We study these as implemented in the code GIZMO for both new meshless finite-volume Godunov schemes (MFM/MFV). We show that the MFM/MFV methods are accurate and stable even with noisy fields and irregular particle arrangements, and recover the correct behaviour even in arbitrarily anisotropic cases. They are competitive with state-of-the-art AMR/moving-mesh methods, and can correctly treat anisotropic diffusion-driven instabilities (e.g. the MTI and HBI, Hall MRI). We also develop a new scheme for stabilizing anisotropic tensor-valued fluxes with high-order gradient estimators and non-linear flux limiters, which is trivially generalized to AMR/moving-mesh codes. We also present applications of some of these improvements for SPH, in the form of a new integral-Godunov SPH formulation that adopts a moving-least squares gradient estimator and introduces a flux-limited Riemann problem between particles.

Computational Analysis of the Effect of Porosity on Shock Cell Strength at Cruise

NASA Technical Reports Server (NTRS)

Massey, Steven J.; Elmiligui, Alaa A.; Pao, S. Paul; Abdol-Hamid, Khaled S.; Hunter, Craig A.

2006-01-01

A computational flow field analysis is presented of the effect of core cowl porosity on shock cell strength for a modern separate flow nozzle at cruise conditions. The goal of this study was to identify the primary physical mechanisms by which the application of porosity can reduce shock cell strength and hence the broadband shock associated noise. The flow is simulated by solving the asymptotically steady, compressible, Reynoldsaveraged Navier-Stokes equations on a structured grid using an implicit, up-wind, flux-difference splitting finite volume scheme. The standard two-equation k - epsilon turbulence model with a linear stress representation is used with the addition of a eddy viscosity dependence on total temperature gradient normalized by local turbulence length scale. Specific issues addressed in this study were the optimal area required to weaken a shock impinging on the core cowl surface and the optimal level of porosity and placement of porous areas for reduction of the overall shock cell strength downstream. Two configurations of porosity were found to reduce downstream shock strength by approximately 50%.

DPW-VI Results Using FUN3D with Focus on k-kL-MEAH2015 (k-kL) Turbulence Model

NASA Technical Reports Server (NTRS)

Abdol-Hamid, K. S.; Carlson, Jan-Renee; Rumsey, Christopher L.; Lee-Rausch, Elizabeth M.; Park, Michael A.

2017-01-01

The Common Research Model wing-body configuration is investigated with the k-kL-MEAH2015 turbulence model implemented in FUN3D. This includes results presented at the Sixth Drag Prediction Workshop and additional results generated after the workshop with a nonlinear Quadratic Constitutive Relation (QCR) variant of the same turbulence model. The workshop provided grids are used, and a uniform grid refinement study is performed at the design condition. A large variation between results with and without a reconstruction limiter is exhibited on "medium" grid sizes, indicating that the medium grid size is too coarse for drawing conclusions in comparison with experiment. This variation is reduced with grid refinement. At a fixed angle of attack near design conditions, the QCR variant yielded decreased lift and drag compared with the linear eddy-viscosity model by an amount that was approximately constant with grid refinement. The k-kL-MEAH2015 turbulence model produced wing root junction flow behavior consistent with wind tunnel observations.

Grid Convergence for Turbulent Flows(Invited)

NASA Technical Reports Server (NTRS)

Diskin, Boris; Thomas, James L.; Rumsey, Christopher L.; Schwoppe, Axel

2015-01-01

A detailed grid convergence study has been conducted to establish accurate reference solutions corresponding to the one-equation linear eddy-viscosity Spalart-Allmaras turbulence model for two dimensional turbulent flows around the NACA 0012 airfoil and a flat plate. The study involved three widely used codes, CFL3D (NASA), FUN3D (NASA), and TAU (DLR), and families of uniformly refined structured grids that differ in the grid density patterns. Solutions computed by different codes on different grid families appear to converge to the same continuous limit, but exhibit different convergence characteristics. The grid resolution in the vicinity of geometric singularities, such as a sharp trailing edge, is found to be the major factor affecting accuracy and convergence of discrete solutions, more prominent than differences in discretization schemes and/or grid elements. The results reported for these relatively simple turbulent flows demonstrate that CFL3D, FUN3D, and TAU solutions are very accurate on the finest grids used in the study, but even those grids are not sufficient to conclusively establish an asymptotic convergence order.

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

NASA Technical Reports Server (NTRS)

Steinthorsson, E.; Liou, M. S.; Povinelli, L. A.

1993-01-01

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping, and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results shown are a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

NASA Technical Reports Server (NTRS)

Steinthorsson, E.; Liou, M.-S.; Povinelli, L. A.

1993-01-01

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results are shown a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

NASA Technical Reports Server (NTRS)

Towne, Charles E.; Schwab, John R.; Benson, Thomas J.; Suresh, Ambady

1990-01-01

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Turbulence Model Predictions of Strongly Curved Flow in a U-Duct

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Gatski, Thomas B.; Morrison, Joseph H.

2000-01-01

The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied. An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation. In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does. Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model. Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature. By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall.

On the eddy-resolving capability of high-order discontinuous Galerkin approaches to implicit LES / under-resolved DNS of Euler turbulence

NASA Astrophysics Data System (ADS)

Moura, R. C.; Mengaldo, G.; Peiró, J.; Sherwin, S. J.

2017-02-01

We present estimates of spectral resolution power for under-resolved turbulent Euler flows obtained with high-order discontinuous Galerkin (DG) methods. The '1% rule' based on linear dispersion-diffusion analysis introduced by Moura et al. (2015) [10] is here adapted for 3D energy spectra and validated through the inviscid Taylor-Green vortex problem. The 1% rule estimates the wavenumber beyond which numerical diffusion induces an artificial dissipation range on measured energy spectra. As the original rule relies on standard upwinding, different Riemann solvers are tested. Very good agreement is found for solvers which treat the different physical waves in a consistent manner. Relatively good agreement is still found for simpler solvers. The latter however displayed spurious features attributed to the inconsistent treatment of different physical waves. It is argued that, in the limit of vanishing viscosity, such features might have a significant impact on robustness and solution quality. The estimates proposed are regarded as useful guidelines for no-model DG-based simulations of free turbulence at very high Reynolds numbers.

Structure-Preserving Variational Multiscale Modeling of Turbulent Incompressible Flow with Subgrid Vortices

NASA Astrophysics Data System (ADS)

Evans, John; Coley, Christopher; Aronson, Ryan; Nelson, Corey

2017-11-01

In this talk, a large eddy simulation methodology for turbulent incompressible flow will be presented which combines the best features of divergence-conforming discretizations and the residual-based variational multiscale approach to large eddy simulation. In this method, the resolved motion is represented using a divergence-conforming discretization, that is, a discretization that preserves the incompressibility constraint in a pointwise manner, and the unresolved fluid motion is explicitly modeled by subgrid vortices that lie within individual grid cells. The evolution of the subgrid vortices is governed by dynamical model equations driven by the residual of the resolved motion. Consequently, the subgrid vortices appropriately vanish for laminar flow and fully resolved turbulent flow. As the resolved velocity field and subgrid vortices are both divergence-free, the methodology conserves mass in a pointwise sense and admits discrete balance laws for energy, enstrophy, and helicity. Numerical results demonstrate the methodology yields improved results versus state-of-the-art eddy viscosity models in the context of transitional, wall-bounded, and rotational flow when a divergence-conforming B-spline discretization is utilized to represent the resolved motion.

Subgrid-scale models for large-eddy simulation of rotating turbulent flows

NASA Astrophysics Data System (ADS)

Silvis, Maurits; Trias, Xavier; Abkar, Mahdi; Bae, Hyunji Jane; Lozano-Duran, Adrian; Verstappen, Roel

2016-11-01

This paper discusses subgrid models for large-eddy simulation of anisotropic flows using anisotropic grids. In particular, we are looking into ways to model not only the subgrid dissipation, but also transport processes, since these are expected to play an important role in rotating turbulent flows. We therefore consider subgrid-scale models of the form τ = - 2νt S +μt (SΩ - ΩS) , where the eddy-viscosity νt is given by the minimum-dissipation model, μt represents a transport coefficient; S is the symmetric part of the velocity gradient and Ω the skew-symmetric part. To incorporate the effect of mesh anisotropy the filter length is taken in such a way that it minimizes the difference between the turbulent stress in physical and computational space, where the physical space is covered by an anisotropic mesh and the computational space is isotropic. The resulting model is successfully tested for rotating homogeneous isotropic turbulence and rotating plane-channel flows. The research was largely carried out during the CTR SP 2016. M.S, and R.V. acknowledge the financial support to attend this Summer Program.

Viscosity-dependent diffusion of fluorescent particles using fluorescence correlation spectroscopy.

PubMed

Jung, Chanbae; Lee, Jaeran; Kang, Manil; Kim, Sok Won

2014-11-01

Fluorescent particles show the variety characteristics by the interaction with other particles and solvent. In order to investigate the relationship between the dynamic properties of fluorescent particles and solvent viscosity, particle diffusion in various solvents was evaluated using a fluorescence correlation spectroscopy. Upon analyzing the correlation functions of AF-647, Q-dot, and beads with different viscosity values, the diffusion time of all particles was observed to increase with increasing solvent viscosity, and the ratio of diffusion time to solvent viscosity, τ D /η, showed a linear dependence on particle size. The particle diffusion coefficients calculated from the diffusion time decreased with increasing solvent viscosity. Further, the hydrodynamic radii of AF-647, Q-dot, and beads were 0.98 ± 0.1 nm, 64.8 ± 3.23 nm, and 89.8 ± 4.91 nm, respectively, revealing a linear dependence on τ D /η, which suggests that the hydrodynamic radius of a particle strongly depends on both the physical size of the particle and solvent viscosity.

Anomalous intrinsic viscosity of octadecylamine-functionalised carbon nanotubes in suspension.

PubMed

Donovan, K J; Scott, K

2013-06-28

Single walled carbon nanotubes, SWCNTs, are used as a model cylinder of nanoscopic dimensions for testing rheological theories of how addition of cylindrical particles affects the viscosity of a suspension of such particles. Using the rate of growth of the accompanying induced linear dichroism following application of an applied electric field, the dynamics of carbon nanotube alignment is studied in suspensions of octadecylamine functionalised single walled carbon nanotubes, ODA-SWCNTs, in 1,2 dichloroethane. From such measurements the viscosity of the suspension is measured as the concentration of the suspension is varied. While working within the dilute limit the viscosity is found to increase linearly with concentration and the intrinsic viscosity of the suspension is found to be 8000. This anomalously high intrinsic viscosity is compared with the predictions of various models for a rigid cylinder and found to be incompatible with any of the current models. Some suggestions are made as to the way this ODA-SWCNT result may be eventually accommodated within other models.

Testing Munk's hypothesis for submesoscale eddy generation using observations in the North Atlantic

NASA Astrophysics Data System (ADS)

Buckingham, Christian E.; Khaleel, Zammath; Lazar, Ayah; Martin, Adrian P.; Allen, John T.; Naveira Garabato, Alberto C.; Thompson, Andrew F.; Vic, Clément

2017-08-01

A high-resolution satellite image that reveals a train of coherent, submesoscale (6 km) vortices along the edge of an ocean front is examined in concert with hydrographic measurements in an effort to understand formation mechanisms of the submesoscale eddies. The infrared satellite image consists of ocean surface temperatures at ˜390 m resolution over the midlatitude North Atlantic (48.69°N, 16.19°W). Concomitant altimetric observations coupled with regular spacing of the eddies suggest the eddies result from mesoscale stirring, filamentation, and subsequent frontal instability. While horizontal shear or barotropic instability (BTI) is one mechanism for generating such eddies (Munk's hypothesis), we conclude from linear theory coupled with the in situ data that mixed layer or submesoscale baroclinic instability (BCI) is a more plausible explanation for the observed submesoscale vortices. Here we assume that the frontal disturbance remains in its linear growth stage and is accurately described by linear dynamics. This result likely has greater applicability to the open ocean, i.e., regions where the gradient Rossby number is reduced relative to its value along coasts and within strong current systems. Given that such waters comprise an appreciable percentage of the ocean surface and that energy and buoyancy fluxes differ under BTI and BCI, this result has wider implications for open-ocean energy/buoyancy budgets and parameterizations within ocean general circulation models. In summary, this work provides rare observational evidence of submesoscale eddy generation by BCI in the open ocean.