Anisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
|
|
- Hector Fletcher
- 5 years ago
- Views:
Transcription
1 Center for Turbulence Research Annual Research Briefs Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale models have long been recognized as a natural alternative to the isotropic Smagorinsky model. A priori tests, based either on experiments (Liu et al. 1994) or on DNS data have often shown that these models have much better correlations than the Smagorinsky model. However their ecacy in large eddy simulations has so far been limited by two factors. First, these models are in general viewed, through Taylor expansions, as a way to compute the subgridscale contribution of the Leonard term. They are used with the coecient values which result from these Taylor expansions and thus depend on the form of the lter function. Secondly, anisotropic models are based on tensor forms of the turbulent viscosity and therefore may produce backscatter as well as dissipation. The clipping techniques, which are proposed (Liu et al. 1994, Vreman et al. 1997) to overcome the destabilizing eects of backscatter, are based on energy balances but ignore that dissipation and backscatter in general coexist at every point along dierent strain directions of the ow. They thus lead to models which are not dissipative enough, which explains why they are often complemented by Smagorinsky terms in so-called mixed models. These formulations, however, are not able to retain the potential gain oered by the anisotropy of the original model, in particular in laminar or wall-bounded ows. Our goal here is twofold. We rst derive formulas based on the quantities computed on the grid that facilitate the implementation of the anistropic model in any code. By distinguishing between backscatter and dissipation directions in the ow, we also present strictly dissipative formulas which lead to stable and truly anisotropic schemes. 2. Accomplishments 2.1 The model Our starting point will be the following model for the residual shear stresses, sometimes referred to as the self-similarity model (Liu et al. 1994) or the gradient model (Vreman et al. 1997): ij ' C 2 D ik ud jk u (1) where D ik u i and is the lter width. Throughout the paper, we the convention of summation of repeated indices. 1 Permanent address: LMC-IMAG, Universite Joseph Fourier, BP 53 Grenoble Cedex 9, France
2 114 G.-H. Cottet & A. A. Wray Integral approximation and grid formulas Z Our derivation is based on a numerical lter function satisfying the following moment conditions: x k x l (x) dx = kl k l=1 2 3 (2) where x k denote the components of x and kl is the Kronecker symbol. Such a lter can be easily constructed by proper rescaling of any positive function satisfying symmetry properties (e.g. functions with spherical symmetry or functions constructed through tensor product of one-dimensional even functions). We rst write (here and in the sequel we drop the overbar notation for the resolved elds) ;3Z 2 D ik u(x)d jk u(x) = D ik u(x)d jl u(x)(y k ; x k )(y l ; x l )( y ; x ) dy (3) where we recall that the summation of repeated indices is implied. In the above formula, is the lter width, which for the time being is assumed to be constant. By Taylor expansions of u i and u j around x, this yields 2 D ik u(x)d jk u(x) = ;3Z [u j (y) ; u j (x)][u i (y) ; u i (x)]( y ; x ) dy + O(2 ) : (4) Since ultimately it is the divergence of ij that we need to model in order to solve the ltered Navier-Stokes equations, we take the divergence of (4) to obtain, after cancellation of the term involving the divergence of j [D ik u(x)d jk u(x)] ' A i + B i where we used the j j and ;3Z A i = ; [u j (y) ; u j (x)]@ j u i (x)( y ; x ) dy ;4Z B i = ; [u j (y) ; u j (x)][u i (y) ; u i (x)]@ j ( y ; x ) dy R It is readily seen that A i are convective terms: if one sets = (y) dy and Z ^u(x) = 1 3 u(y)( y ; x ) dy then A i can be rewritten as (^u ; u)ru i.itthus does not contribute to the energy balance. Since the goal of SGS models is to model the transfer of energy between
3 Anisotropic subgrid-scale models 115 large and small scales, we will only be interested in B i. Hence the SGS model can be expressed as j ij ';C [u j (y) ; u j (x)][u i (y) ; u i (x)]@ j ( y ; x ) dy (5) Viscous scales and grid renement The above derivation assumes an isotropic lter. In practice, it is natural to use a function which matches the computational grid. For a variable grid such as those generally used for wall-bounded ows, one can then expect that grid renement in specic directions will translate into variable viscous length scales in the subgridscale model. To be more specic let us assume that the lter is built from a unique one-dimensional shape in all 3 directions and that 1(x1) 2(x2) 3(x3) are the local grid sizes in the 3 directions. The natural generalization of formulas (4) and (5) to this case is given by ij (x) j ij (x) ' C!(x) Z C 1!(x) j (x j ) [u j (y) ; u j (x)][u i (y) ; u i (x)](x y) dy Z [u j (y) ; u j (x)][u i (y) ; u i (x)] j (x y) dy where!(x) =1(x1)2(x2)3(x3) is the volume of the cell centered at x and (x y) =( x 1 ; y1 1(x) )( x 2 ; y2 2(x) )( x 3 ; y3 3(x) ) 1(x y) = 0 ( x 1 ; y1 1(x) )( x 2 ; y2 2(x) )( x 3 ; y3 3(x) ): and similar formulas for 2 3. It can be shown (see Cottet 1997 for details) by using a coordinate mapping between the mesh and a uniform isotropic grid that these formulas correspond to the dierential subgrid-scale model ij ' k (x) 2 D ik ud jk u In other words, the eddy-viscosity length scale is given by the grid spacing in the corresponding direction. In near wall regions, when grid renement is used, this produces an additional damping in the SGS dissipation which we believe is desirable. Note that the traditional derivation of subgrid-scale models does not apply to varying size lters, making the use of (1) questionable Anisotropic clipping The total subgrid-scale dissipation D associated to the model (5) is obtained by multiplying this formula by u(x) and integrating over x. Writing u(x) = 1 2 [u(x)+ u(y)] [u(x) ; u(y)] and using the symmetry of, we are left with ;4Z D = C [u(y) ; u(x)] r( y ; x 2 )ju(y) ; u(x)j2 dx dy
4 116 G.-H. Cottet & A. A. Wray 2 T 1 3 α =u (x ) u (x) T 2 T 2 N α =u (x ) u (x) E 1 E 1 α = u (x )+u (x) S 3 S 3... W S B E Figure 1. Finite-dierence stencil of the subgrid-scale model in the case of a uniform grid with a radial hat lter. T = E = denote the weights of the corresponding grid points. This formula gives a way to measure the local dissipation at a given point x as well by restricting the integral over y only. It also enables us to formulate a strictly dissipative model as follows j ij (x) ' [u(x) ; u(y)] r( x ; y + )] [u(x) ; u(y)] dy (6) where a+ = max(0 a). Unlike the traditional clipping strategies, which would consist in j ij by zero whenever the global energy budget would have the wrong sign, this technique respects the anisotropy of the original formula in the sense that it allows to dissipate in one or more directions while controlling the backscatter which would arise in the other directions. A clipping technique can be written along the same lines to make sure that when combining the subgrid-scale model with the eect of the molecular viscosity the method is strictly dissipative Isotropic turbulence 2.2 Large eddy simulations We rst examine the validity of our model in simulations of decaying isotropic turbulence. All calculations were done with a spectral code in a periodic box with full dealiasing. The subgrid-scale models is therefore implemented on a uniform grid on which velocity values are classically obtained through FFT. We have considered two possible choices for the numerical lter : a radial function 1(x) = f(jxj) and a tensor product function 2(x) = f(x1)f(x2)f(x3). The function f is a piecewise quadratic spline with support in [;1:5 +1:5], and
5 Anisotropic subgrid-scale models Resolved energy time Figure 2. Energy decay for DNS and LES with Smagorinsky and anisotropic models. : DNS results : Smagorinsky model : Anisotropic model with 27 points lter : Anisotropic model with 7 points lter. Resolved enstrophy time Figure 3. Enstrophy decay for DNS and LES with Smagorinsky and anisotropic models. See Fig. 2 for symbol legend.
6 118 G.-H. Cottet & A. A. Wray -2 Energy k Figure 4. legend. Energy spectra at t =4: for DNS and LES. See Fig. 2 for symbol -2 Energy k Figure 5. Energy spectra at t = 7. See Fig. 2 for symbol legend.
7 0.0 Anisotropic subgrid-scale models SGS Dissipation Y Figure 6. Subgrid-scale dissipation for the dynamic model and anistropic models. : Dynamic model : Unclipped C=.25 : Clipped C=.075. are 2 normalization constants needed to enforce (2). As a result, the stencil corresponding to the implementation of (5) involves 7 points in the rst case and 27 in the second one. Figure 1 sketches the particular form of the stencil in the rst case. We have compared the results of LES on a 32 3 grid using the clipped anisotropic model (6), as well as the Smagorinsky model, with a DNS. The coecients of the various models have all been tuned to yield an energy decay which matches the DNS results (Fig. 2). The enstrophy decay curves in Fig. 3 show that the anisotropic models behave better than the Smagorinsky model up to time T 5. This is conrmed by the energy spectra showed in Fig. 4. Past this time, all models are too dissipative in the high modes, with slightly better results for the anisotropic model (Fig. 5). The particular implementation chosen for the anisotropic model (isotropic vs tensor-product lter form) does not seem to signicantly aect the results Channel calculations This case is more challenging as it is well known that the Smagorinsky model is unable in this geometry to give good results in the absence of ad hoc damping at the walls or dynamic coecient calculation. We show comparisons between the anisotropic model and the so-called global
8 120 G.-H. Cottet & A. A. Wray VRMS WRMS URMS Y+ Figure 7. Turbulence intensities in wall coordinates top curves: streamwise velocity, middle curves: normal velocity, bottom curves: spanwise velocity. Symbols same as Fig. 6. dynamic model (Germano et al. 1991), which can considered as the best available model for this geometry. Our tests are done in the context of a Tchebytchev-Fourier collocation code (Kim et al. 1987). In this numerical scheme, the Crank-Nicolson time-advancing scheme is used for diusion in the wall normal direction together with a third order Runge-Kutta method to advance the nonlinear convection and SGS terms. Periodic boundary conditions are assumed in the streamwise and spanwise directions and no-slip conditions at the walls located at y = 1. Dealiasing is performed in the periodic directions, and the pressure is updated at each iteration to maintain a constant momentum throughout the calculation. We have focused on the case of a Reynolds number of Re =1 030 based on the shear velocity, which corresponds to a Reynolds number of about based on the centerline velocity and the channel half-width. The anistropic model has been implemented with the tensor-product lter 2. Note that the grid renement given by the Tchebytchev collocation points in the wall normal direction implies a variable eddy-viscosity length scale in this direction as explained in Two cases have been considered: in the rst one the method
9 Anisotropic subgrid-scale models < UV > Y+ Figure 8. Large scale shear stress in wall coordinates. Symbls same as Fig. 6. (5) is used without any clipping and a coecient C =0:25, while in the second a method similar to (6) is used in order to prevent the backscatter from overwhelming the molecular dissipation. In this latter case, a smaller value of the coecient was chosen (C = 0:075). Figure 6 shows the subgrid-scale dissipation averaged on planes parallel to the walls. The similarity of the proles obtained by the anisotropic models and the dynamic model is striking and conrms that the anisotropic model has a much better behavior in near wall regions than the Smagorinsky model. This observation is conrmed by the results obtained for the velocity uctuations in Fig. 7 and the shear stresses in Fig. 8. We refer to (Cottet 1997) for more numerical results for this case. 3. Future plans Tests of the anisotropic model have shown the superiority of this model over the Smagorinsky model. This superiority, which can be expected from the a priori tests that can be found in the literature, has been conrmed in LES using the present method. The ecacy of the model has gained from its ability to incorporate a truly anistropic backscatter control. The performance of the model both in isotropic and wall-bounded ows encourage trying it for other ows such as shear layers or jets.
10 122 G.-H. Cottet & A. A. Wray Another direction of research is the implementation of a dynamic procedure to compute the model coecient, along the lines suggested by Cottet (1997). On the basis of the present results, one can expect that this technique should not produce highly oscillatory values for the coecient. In other words, the anisotropic model should be better conditioned than the Smagorinsky model for a dynamic procedure, in particular when no averaging over homogeneous directions can be used to stabilize it. REFERENCES Cottet, G. H Anisotropic subgrid-scale models for Large Eddy Simulations of turbulent ows. submitted. Germano, M., Piomelli, U., Moin, P. & Cabot, W. H A dynamic subgrid-scale eddy viscosity model. Phys. Fluids A. 3, Kim, J., Moin, P. &Moser, R Turbulence statistics in fully developed channel ow atlow Reynolds number. J. Fluid Mech. 177, Liu, S., Meneveau, C. & Katz, J The property of similarity subgrid-scale models. J. Fluid Mech. 275, Vreman, V., Geurts, B.& Kuerten, H Large-eddy simulations of the turbulent mixing layer. J. Fluid Mech. 339,
model and its application to channel ow By K. B. Shah AND J. H. Ferziger
Center for Turbulence Research Annual Research Briefs 1995 73 A new non-eddy viscosity subgrid-scale model and its application to channel ow 1. Motivation and objectives By K. B. Shah AND J. H. Ferziger
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
More informationEnsemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3
Center for Turbulence Research Proceedings of the Summer Program 1996 237 Ensemble averaged dynamic modeling By D. Carati 1,A.Wray 2 AND W. Cabot 3 The possibility of using the information from simultaneous
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationOn the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows
PHYSICS OF FLUIDS VOLUME 11, NUMBER 5 MAY 1999 On the relationship between the mean flow and subgrid stresses in large eddy simulation of turbulent shear flows L. Shao a) Laboratoire de Mécanique des Fluides
More informationA priori testing of subgrid-scale. models for the velocity-pressure. and vorticity-velocity formulations
Center for Turbulence Research Proceedings of the Summer Program 1996 309 A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations By G. S. Winckelmans 1,
More informationOn the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows
Center for Turbulence Research Annual Research Briefs 1998 267 On the feasibility of merging LES with RANS for the near-wall region of attached turbulent flows By Jeffrey S. Baggett 1. Motivation and objectives
More informationLarge eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions
June 30 - July 3, 2015 Melbourne, Australia 9 P-26 Large eddy simulation of turbulent flow over a backward-facing step: effect of inflow conditions Jungwoo Kim Department of Mechanical System Design Engineering
More informationTesting of a new mixed model for LES: the Leonard model supplemented by a dynamic Smagorinsky term
Center for Turbulence Research Proceedings of the Summer Program 1998 367 Testing of a new mixed model for LES: the Leonard model supplemented by a dynamic Smagorinsky term By G. S. Winckelmans 1,A.A.Wray
More informationCOMPARISON OF DIFFERENT SUBGRID TURBULENCE MODELS AND BOUNDARY CONDITIONS FOR LARGE-EDDY-SIMULATIONS OF ROOM AIR FLOWS.
7 TH INTRNATINAL CNFRNC N AIR DISTRIBTIN IN RMS, RMVNT 2 pp. 31-36 CMPARISN F DIFFRNT SBGRID TRBLNC MDLS AND BNDARY CNDITINS FR LARG-DDY-SIMLATINS F RM AIR FLWS. D. Müller 1, L. Davidson 2 1 Lehrstuhl
More informationOn modeling pressure diusion. in non-homogeneous shear ows. By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3
Center for Turbulence Research Proceedings of the Summer Program 1996 63 On modeling pressure diusion in non-homogeneous shear ows By A. O. Demuren, 1 M. M. Rogers, 2 P. Durbin 3 AND S. K. Lele 3 New models
More informationThree-dimensional wall filtering formulation for large-eddy simulation
Center for Turbulence Research Annual Research Briefs 6 55 Three-dimensional wall filtering formulation for large-eddy simulation By M. Shoeybi AND J. A. Templeton 1. Motivation and objectives Large-eddy
More informationAnalysis and modelling of subgrid-scale motions in near-wall turbulence
J. Fluid Mech. (1998), vol. 356, pp. 327 352. Printed in the United Kingdom c 1998 Cambridge University Press 327 Analysis and modelling of subgrid-scale motions in near-wall turbulence By CARLOS HÄRTEL
More informationAn improved velocity increment model based on Kolmogorov equation of filtered velocity
An improved velocity increment model based on Kolmogorov equation of filtered velocity Le Fang, Liang Shao, Jean-Pierre Bertoglio, Guixiang X. Cui, Chun-Xiao Xu, Zhaoshun Zhang To cite this version: Le
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW DATA SET Data provenance: J. Graham 1, M. Lee 2, N. Malaya 2, R.D. Moser 2, G. Eyink 1 & C. Meneveau 1 Database ingest and Web Services: K. Kanov
More informationThe JHU Turbulence Databases (JHTDB)
The JHU Turbulence Databases (JHTDB) TURBULENT CHANNEL FLOW AT Re τ = 5200 DATA SET Data provenance: M. Lee 1 & R. D. Moser 1 Database ingest and Web Services: Z. Wu 2, G. Lemson 2, R. Burns 2, A. Szalay
More informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More informationSharath S. Girimaji and Ye Zhou. NASA Langley Research Center, Hampton, VA Abstract
Analysis and modeling of subgrid scalar mixing using numerical data Sharath S. Girimaji and Ye hou Institute for Computer Applications in Science and Engineering NASA Langley Research Center, Hampton,
More informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationValidation of an Entropy-Viscosity Model for Large Eddy Simulation
Validation of an Entropy-Viscosity Model for Large Eddy Simulation J.-L. Guermond, A. Larios and T. Thompson 1 Introduction A primary mainstay of difficulty when working with problems of very high Reynolds
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More informationLES AND RANS STUDIES OF OSCILLATING FLOWS OVER FLAT PLATE
LES AND RANS STUDIES OF OSCILLATING FLOWS OVER FLAT PLATE By Chin-Tsau Hsu, 1 Xiyun Lu, and Man-Kim Kwan 3 ABSTRACT: Oscillatory flows over a flat plate are studied by using Large Eddy Simulation (LES)
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationOn the Euler rotation angle and axis of a subgrid-scale stress model
Center for Turbulence Research Proceedings of the Summer Program 212 117 On the Euler rotation angle and axis of a subgrid-scale stress model By M. Saeedi, B.-C. Wang AND S. T. Bose In this research, geometrical
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More information43rd AIAA Aerospace Sciences Meeting and Exhibit, Jan 2005, Reno, Nevada
43rd AIAA Aerospace Sciences Meeting and Exhibit, 10-13 Jan 2005, Reno, Nevada A Dynamic Procedure for the Lagrangian Averaged Navier-Stokes-α Model of Turbulent Flows Kamran Mohseni and Hongwu Zhao Aerospace
More informationDatabase analysis of errors in large-eddy simulation
PHYSICS OF FLUIDS VOLUME 15, NUMBER 9 SEPTEMBER 2003 Johan Meyers a) Department of Mechanical Engineering, Katholieke Universiteit Leuven, Celestijnenlaan 300, B3001 Leuven, Belgium Bernard J. Geurts b)
More informationBy B. Perot 1 AND P. Moin 2. A new approach to Reynolds averaged turbulence modeling is proposed which has
Center for Turbulence Research Proceedings of the Summer Program 1996 35 A new approach to turbulence modeling By B. Perot 1 AND P. Moin 2 A new approach to Reynolds averaged turbulence modeling is proposed
More informationEffects of Eddy Viscosity on Time Correlations in Large Eddy Simulation
NASA/CR-2-286 ICASE Report No. 2- Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation Guowei He ICASE, Hampton, Virginia R. Rubinstein NASA Langley Research Center, Hampton, Virginia
More informationInvariant modeling in large-eddy. simulation of turbulence. By M. Oberlack
Center for Turbulence Research Annual Research Briefs 1997 3 Invariant modeling in large-eddy simulation of turbulence 1. Motivation and objectives By M. Oberlack Since the derivation of the Smagorinsky
More information2 GOVERNING EQUATIONS
2 GOVERNING EQUATIONS 9 2 GOVERNING EQUATIONS For completeness we will take a brief moment to review the governing equations for a turbulent uid. We will present them both in physical space coordinates
More informationApplication of turbulence. models to high-lift airfoils. By Georgi Kalitzin
Center for Turbulence Research Annual Research Briefs 997 65 Application of turbulence models to high-lift airfoils By Georgi Kalitzin. Motivation and objectives Accurate prediction of the ow over an airfoil
More informationGeneration of initial fields for channel flow investigation
Generation of initial fields for channel flow investigation Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (Mai 2) In the framework of the DFG-funded research
More informationShear-Improved Smagorinsky Model for Large-Eddy Simulation of Wall-Bounded Turbulent Flows
Under consideration for publication in J. Fluid Mech. Shear-Improved Smagorinsky Model for Large-Eddy Simulation of Wall-Bounded Turbulent Flows By E. Lévêque, F. Toschi,3, L. Shao 4 and J.-P. Bertoglio
More informationSome remarks on grad-div stabilization of incompressible flow simulations
Some remarks on grad-div stabilization of incompressible flow simulations Gert Lube Institute for Numerical and Applied Mathematics Georg-August-University Göttingen M. Stynes Workshop Numerical Analysis
More informationPredicting natural transition using large eddy simulation
Center for Turbulence Research Annual Research Briefs 2011 97 Predicting natural transition using large eddy simulation By T. Sayadi AND P. Moin 1. Motivation and objectives Transition has a big impact
More informationSubgrid-Scale Models for Compressible Large-Eddy Simulations
Theoret. Comput. Fluid Dynamics (000 13: 361 376 Theoretical and Computational Fluid Dynamics Springer-Verlag 000 Subgrid-Scale Models for Compressible Large-Eddy Simulations M. Pino Martín Department
More informationA Finite-Element based Navier-Stokes Solver for LES
A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,
More informationDefense Technical Information Center Compilation Part Notice
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP013647 TITLE: A Dynamic Procedure for Calculating the Turbulent Kinetic Energy DISTRIBUTION: Approved for public release, distribution
More informationTypsatt i AMS-L A TEX. Akademisk avhandling som med tillstand av Kungliga Tekniska Hogskolan i Stockholm framlagges till oentlig granskning for avlagg
Studies of turbulence and its modelling through large eddy- and direct numerical simulation by Krister Alvelius October 1999 Technical Reports from Royal Institute of Technology Department ofmechanics
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More information258 S. Ghosal & M. M. Rogers addition to the assumption of self-similarity, one often invokes additional restrictions inspired by the dictum \turbulen
Center for Turbulence Research Annual Research Briefs 1996 A numerical study of self-similarity in a turbulent plane wake using large-eddy simulation By Sandip Ghosal 1 AND Michael M. Rogers 2 257 1. Motivation
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationValidation of large-eddy simulation. in a plain asymmetric diuser. By M. Fatica, H.-J. Kaltenbach 1 AND R. Mittal 2
Center for Turbulence Research Annual Research Briefs 1997 23 Validation of large-eddy simulation in a plain asymmetric diuser By M. Fatica, H.-J. Kaltenbach 1 AND R. Mittal 2 Motivation The main motivation
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
More informationIMPACT OF LES TURBULENCE SUBGRID MODELS IN THE JET RELEASE SIMULATION
IMPACT OF LES TURBULENCE SUBGRID MODELS IN THE JET RELEASE SIMULATION E. S. FERREIRA JÚNIOR 1, S. S. V. VIANNA 1 1 State University of Campinas, Faculty of Chemical Engineering E-mail: elmo@feq.unicamp.br,
More informationLARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL
LARGE EDDY SIMULATION AND FLOW CONTROL OVER A 25 RAMP MODEL 09/11/2017 Paolo Casco Stephie Edwige Philippe Gilotte Iraj Mortazavi LES and flow control over a 25 ramp model : context 2 Context Validation
More informationCOMMUTATION ERRORS IN PITM SIMULATION
COMMUTATION ERRORS IN PITM SIMULATION B. Chaouat ONERA, 93 Châtillon, France Bruno.Chaouat@onera.fr Introduction Large eddy simulation is a promising route. This approach has been largely developed in
More informationJ OURNAL OF TURBULENCE. JoT 5 (2004) 020. Large-eddy simulation of heat transfer downstream of a backward-facing step
JOT J OURNAL OF TURBULENCE http://jot.iop.org/ Large-eddy simulation of heat transfer downstream of a backward-facing step A Keating 1,4, U Piomelli 2, K Bremhorst 1 andsnešić 3 1 Division of Mechanical
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationLarge Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena
29 Review Large Eddy Simulation as a Powerful Engineering Tool for Predicting Complex Turbulent Flows and Related Phenomena Masahide Inagaki Abstract Computational Fluid Dynamics (CFD) has been applied
More informationLES and unsteady RANS of boundary-layer transition induced by periodically passing wakes
Center for Turbulence Research Proceedings of the Summer Program 2 249 LES and unsteady RANS of boundary-layer transition induced by periodically passing wakes By F. E. Ham, F. S. Lien, X. Wu, M. Wang,
More informationA scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric boundary layer
J. Fluid Mech. (2000), vol. 415, pp. 261 284. Printed in the United Kingdom c 2000 Cambridge University Press 261 A scale-dependent dynamic model for large-eddy simulation: application to a neutral atmospheric
More informationTime / T p
Chapter 5 Single-Point QNS Measurements 5. Overview A total of 46 primary test cases were evaluated using the QNS technique. These cases duplicated the 8 studied by Ainley [] and extend the parameters
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationRECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD
RECONSTRUCTION OF TURBULENT FLUCTUATIONS FOR HYBRID RANS/LES SIMULATIONS USING A SYNTHETIC-EDDY METHOD N. Jarrin 1, A. Revell 1, R. Prosser 1 and D. Laurence 1,2 1 School of MACE, the University of Manchester,
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationDynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh
40th Fluid Dynamics Conference and Exhibit 8 June - July 00, Chicago, Illinois AIAA 00-506 Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows Xiaochuan Chai and Krishnan Mahesh University
More informationOn why dynamic subgrid-scale models work
/ / Center for Turbulence Research Annual Research Briefs 1995 25 On why dynamic subgrid-scale models work By J. Jim_nez 1 1. Motivation Dynamic subgrid models were introduced in (Germano et al. 1991)
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationarxiv:nlin/ v2 [nlin.cd] 1 Jun 2006
Under consideration for publication in J. Fluid Mech. arxiv:nlin/6553v [nlin.cd] Jun 6 Shear-Improved Smagorinsky Model for Large-Eddy Simulation of Wall-Bounded Turbulent Flows By E. Lévêque, F. Toschi,3,
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationEstimation of Turbulent Dissipation Rate Using 2D Data in Channel Flows
Proceedings of the 3 rd World Congress on Mechanical, Chemical, and Material Engineering (MCM'17) Rome, Italy June 8 10, 2017 Paper No. HTFF 140 ISSN: 2369-8136 DOI: 10.11159/htff17.140 Estimation of Turbulent
More informationHeating effects on the structure of noise sources of high-speed jets
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 2009, Orlando, Florida AIAA 2009-291 47th Aerospace Sciences Meeting and Exhibit 5 8 January 2009,
More informationwhere the viscosity coecients are functions of the invariants of, A = 1 + rv, and are
Appendix: Final form of the algebraic mass ux model The pertinent facts concerning the computational implementation of the mass ux representation are presented here. The general form of the model, valid
More informationMixing Models for Large-Eddy Simulation of Nonpremixed Turbulent Combustion
S. M. debruynkops Lecturer J. J. Riley Professor Department of Mechanical Engineering, University of Washington, Box 35600, Seattle, WA 98195-600 Mixing Models for Large-Eddy Simulation of Nonpremixed
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationA velocity-estimation subgrid model constrained by subgrid scale dissipation
Available online at www.sciencedirect.com Journal of Computational Physics 227 (28) 419 426 www.elsevier.com/locate/jcp A velocity-estimation subgrid model constrained by subgrid scale dissipation Noma
More informationA NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS
June 3 - July 3, 25 Melbourne, Australia 9 B-3 A NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS Ryan N. King Department of Mechanical Engineering University of Colorado Boulder, CO, 839, USA ryan.n.king@colorado.edu
More informationTurbulent boundary layer
Turbulent boundary layer 0. Are they so different from laminar flows? 1. Three main effects of a solid wall 2. Statistical description: equations & results 3. Mean velocity field: classical asymptotic
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationSubgrid-scale modeling of helicity and energy dissipation in helical turbulence
Subgrid-scale modeling of helicity and energy dissipation in helical turbulence Yi Li, 1 Charles Meneveau, 1 Shiyi Chen, 1,2 and Gregory L. Eyink 3 1 Department of Mechanical Engineering and Center for
More informationBoundary-Layer Transition. and. NASA Langley Research Center, Hampton, VA Abstract
Eect of Far-Field Boundary Conditions on Boundary-Layer Transition Fabio P. Bertolotti y Institut fur Stromungsmechanik, DLR, 37073 Gottingen, Germany and Ronald D. Joslin Fluid Mechanics and Acoustics
More informationNon-linear k;";v 2 modeling. with application to high-lift. By F. S. Lien 1 AND P. A. Durbin 2
Center for Turbulence Research Proceedings of the Summer Program 1996 5 Non-linear k;";v 2 modeling with application to high-lift By F. S. Lien 1 AND P. A. Durbin 2 The k;";v 2 model has been investigated
More informationProblem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600
Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation
More informationDrag reduction in turbulent MHD pipe ows. By P. Orlandi 1. This is a preliminary study devoted to verifying whether or not direct simulations
Center for Turbulence Research Proceedings of the Summer Program 1996 447 Drag reduction in turbulent MHD pipe ows By P. Orlandi 1 This is a preliminary study devoted to verifying whether or not direct
More informationLarge-Eddy Simulation of Laminar-Turbulent. Transition in a Swept-Wing Boundary Layer
AIAA 97-0750 Large-Eddy Simulation of Laminar-Turbulent Transition in a Swept-Wing Boundary Layer X. Huai University of Maryland College Park, MD R. D. Joslin NASA Langley Research Center Hampton, VA;
More informationSUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION
SUBGRID MODELS FOR LARGE EDDY SIMULATION: SCALAR FLUX, SCALAR DISSIPATION AND ENERGY DISSIPATION By Sergei G. Chumakov A dissertation submitted in partial fulfillment of the requirements for the degree
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More informationCHALMERS UNIVERSITY OF TECHNOLOGY
Report No. /3 CHALMERS UNIVERSITY OF TECHNOLOGY Department of Thermo and Fluid Dynamics Large-eddy simulations of turbulent flows in stationary and rotating channels and in a stationary square duct by
More informationActive Control of Instabilities in Laminar Boundary-Layer Flow { Part II: Use of Sensors and Spectral Controller. Ronald D. Joslin
Active Control of Instabilities in Laminar Boundary-Layer Flow { Part II: Use of Sensors and Spectral Controller Ronald D. Joslin Fluid Mechanics and Acoustics Division, NASA Langley Research Center R.
More informationCHAPTER 11: REYNOLDS-STRESS AND RELATED MODELS. Turbulent Flows. Stephen B. Pope Cambridge University Press, 2000 c Stephen B. Pope y + < 1.
1/3 η 1C 2C, axi 1/6 2C y + < 1 axi, ξ > 0 y + 7 axi, ξ < 0 log-law region iso ξ -1/6 0 1/6 1/3 Figure 11.1: The Lumley triangle on the plane of the invariants ξ and η of the Reynolds-stress anisotropy
More informationcfl Copyright by Fotini V. Katopodes 2000 All rights reserved.
A theory for the subfilter-scale model in large-eddy simulation Fotini V. Katopodes, Robert L. Street, Joel H. Ferziger March, 2000 Technical Report 2000-K1 Environmental Fluid Mechanics Laboratory Stanford,
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationn i,j+1/2 q i,j * qi+1,j * S i+1/2,j
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/TERMO-5-97 DATE: December 9,997 TITLE A comparison of complete vs. simplied viscous terms in boundary layer
More informationReliability of LES in complex applications
Reliability of LES in complex applications Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) DESIDER Symposium Corfu, June 7-8, 27 Sample of complex flow
More informationInfluence of high temperature gradient on turbulence spectra
NIA CFD Conference, Hampton, August 6-8, 2012 Future Directions in CFD Research, A Modeling and Simulation Conference Influence of high temperature gradient on turbulence spectra Sylvain Serra Françoise
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More informationOn the decay of two-dimensional homogeneous turbulence
On the decay of two-dimensional homogeneous turbulence J. R. Chasnov a) The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong Received 19 June 1996; accepted 3 September
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
More informationA DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW
A DEDICATED LES EXPERIMENTAL DATABASE FOR THE ASSESSMENT OF LES SGS MODELS: THE PULSATILE JET IMPINGEMENT IN TURBULENT CROSS FLOW Hubert Baya Toda Energy Applications Techniques IFP Energie Nouvelles Rueil
More informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
More informationAnalysis of the gradient-diffusion hypothesis in large-eddy simulations based on transport equations
PHYSICS OF FLUIDS 19, 035106 2007 Analysis of the gradient-diffusion hypothesis in large-eddy simulations based on transport equations Carlos B. da Silva a and José C. F. Pereira Instituto Superior Técnico,
More information(U c. t)/b (U t)/b
DYNAMICAL MODELING OF THE LARGE-SCALE MOTION OF A PLANAR TURBULENT JET USING POD MODES. S. Gordeyev 1 and F. O. Thomas 1 University of Notre Dame, Notre Dame, USA University of Notre Dame, Notre Dame,
More informationUncertainty quantification for RANS simulation of flow over a wavy wall
Uncertainty quantification for RANS simulation of flow over a wavy wall Catherine Gorlé 1,2,3, Riccardo Rossi 1,4, and Gianluca Iaccarino 1 1 Center for Turbulence Research, Stanford University, Stanford,
More information