Helsinki University of Technology Laboratory of Applied Thermodynamics
|
|
- Tyler Douglas
- 6 years ago
- Views:
Transcription
1 Helsinki University of Technology Laboratory of pplied Thermodynamics umerical Methods for Coupling the Reynolds-veraged avier Stokes Equations with the Reynolds-Stress Turbulence Model Patrik Rautaheimo 1 and Timo Siikonen 2 Helsinki University of Technology, Espoo, Finland Report o taniemi ISB ISS Research Scientist, Laboratory of pplied Thermodynamics 2 ssociate Professor, Laboratory of pplied Thermodynamics
2 1 Contents omenclature 2 1 Introduction 4 2 Governing Equations 6 21 Flow Equations 6 22 Reynolds-Stress Turbulence Model 7 23 Dissipation Transport Equation 10 3 umerical Method Spatial Discretization Diagonalization of the Flow Equations Rotation perator Boundary Conditions Time Integration Method 20 4 Test Calculations Channel Flow 22 5 Summary 27
3 3 4 ; >? B C KL K 2 omenclature Jacobian matrix Jacobian matrix Jacobian matrix Courant number difference between and ; also Jacobian of the source term total internal energy per unit volume "$% flux vectors in -, - and ' -directions total entalphy *, matrix of Reynolds stresses flux in a given direction in space right eigenvector matrix with the primitive variables transformation matrix from the conservative to the primitive variables 4: Mach number production of the kinetic energy of turbulence; also boundary condition matrix ;=< Prandtl number source term right eigenvector matrix with the conservative variables; also residual?=@ Reynolds number strain-rate; also cell-face area temperature; also rotation matrix; turbulent transport rotation operator vector of the conservative variables 8 cell volume; also vector of the primitive variables D vector of the characteristic variables EGF specific heat in a constant pressure EIH specific heat in a constant volume diffusion of the Reynolds-stresses; also cell thickness K internal energy per unit mass, ie specific internal energy M unit vectors in Cartesian coordinate system heat conductivity; also kinetic energy of turbulence, ; stability factor K unit vector normal to a cell face + static pressure
4 B L L < M L 3 ' + effective + derivative of pressure in constant density + + derivative of pressure in constant internal energy +1 heat flux time velocity components in -, - and -directions ' Cartesian coordinates non-dimensional normal distance from the surface diagonal eigenvalue matrix the velocity pressure-gradient correlation characteristic variable ratio of specific heats E%H Kronecker s delta dissipation of kinetic energy of turbulence eigenvalue of dynamic viscosity kinematic viscosity density Schmidt s number normal or shear stress specific dissipation rate L M Subscripts I ' Superscripts turbulent conditions L -index; summation index M -component of a matrix grid coordinate directions viscous coordinate directions dissipation left-hand side of the cell face right-hand side of the cell face local coordinates; also Roe averaging; corrected dissipation Favre time-averaging operator Time-averaging operator Fluctuating component Fluctuating component
5 4 1 Introduction In recent years, computational fluid dynamics CFD has begun to be an engineering design tool for practical flow problems Some complex aerodynamical configurations have been calculated, eg a delta wing [1] Unfortunately, many systems involve turbulent or transitional flow that cannot be fully understood, and thus, the case is not realistically simulated Direct numerical simulation DS gives an accurate solution also in complex flow situations, but cannot be applied for practical calculations because of an enormous computing cost Large eddy simulation LES has produced promising results, but like DS, LES requires a lot of computing resources and also the theory is not yet fully evaluated Thus, turbulence modeling for Reynolds-averaged avier Stokes equations still plays a very important role in CFD In a simple case, eg, in the case of a 2-D boundary layer flow at a constant pressure, there is no need to use a complicated turbulence model like the Reynoldsstress model RSM For example, the algebraic turbulence model of Baldwin and Lomax [2] has been successfully used for various airfoil sections [3] lgebraic models utilize universal near-wall functions that are not valid in more complex 3- D flows The limitations of the algebraic models drive us to use more general and, unfortunately, also more complex turbulence models Two-equation models, as and models are more general because they take the history of turbulence into account The weakness of these Boussinesq-approximation-based models is that they do not include the effect of different Reynolds stresses Furthermore, the Reynolds-stress tensor is aligned with the mean strain tensor If curvature, pressure change or rotation is introduced, the standard Boussinesq eddy-viscosity approximation fails if no ad hoc modifications are used Unfortunately, these ad hoc modifications are not general The main advantages of the RSM in comparison with the Boussinesq approximation can be found in the history effects of different stress components and in anisotropic assumption RSM also provides more data on turbulence itself than the lower order models From a theoretical point of view, RSM has the potential to solve many complicated engineering flow problems However, RSM also has disadvantages because of complexity, terms that need modelling and many closure coefficients, and an increased number of differential equations In spite of the apparent difficulties, RSM has been used to calculate internal flows for a long time [4] Recently some calculations have been made for external flow as flows over airfoils [5] RSM has mainly been used with incompressible flow In this case, the Reynolds stresses act only in the momentum equations However, in the present case with a
6 compressible flow assumption, the Reynolds stresses act not only in the momentum equation, but also in the energy equation This makes the implementation more difficult Furthermore, with the Reynolds-stress equations the coupling with the flow equations is anisotropic In general, an isotropic assumption is made [6] also with compressible flow nisotropic coupling makes the implicit time integration and upwinding difficult to apply, but in some cases this approach may accelerate the convergence In this study, Shima s low-reynolds number RSM [7] is coupled with a compressible flow solver [8] based on Roe s method [9] new anisotropic coupling method that utilizes the exact equation for the production of turbulence is introduced The method is applied and the results are presented for a flow in a plane channel 5
7 H H H > 6 2 Governing Equations 21 Flow Equations The Reynolds-averaged avier Stokes equations RS, and the equations for the Reynolds stresses and dissipation of turbulence can be written in the following form H H 1' H 21 where is vector of conservative variables,, and are the inviscid fluxes and, and are the viscous fluxes The source-term > has non-zero components for the Reynolds-stress equations Vectors and are % where is density, + is pressure,, and are velocity components in the Cartesian coordinate system,,,,, and are the Reynolds stresses, and the total internal energy defined L is the specific internal energy and is the summation index Here double prime is the Favre-averaged fluctuation component in the direction, tilde denotes Favre averaging, and the bar means ordinary time averaging Vectors and can be written similarly
8 H B B 7 The viscous flux in the L direction is J 24 1% where the viscous stress tensor is defined as 25 and J is the diffusion flux of the kinetic energy of turbulence Diffusion of the Reynolds stresses are included in the source vector > The diffusion coefficients of the dissipation as 26 where is the appropriate Schmidt s number, and is the turbulent viscosity of the fluid determined with the applied turbulence model Heat flux K is obtained by using Boussinesq s approximation and is written as K E F ;=< EGF ;=< 27 where EGF is specific heat The pressure is calculated from an equation of state +, which, for a perfect gas, is written + where is the ratio of specific heats E F E%H Reynolds-Stress Turbulence Model The Reynolds-stress transport can be written in the following Cartesian tensor form ; B J 29
9 8 where ;,, B, and J are the production, the velocity pressure-gradient correlation, the turbulent transport, the dissipation rate and the diffusion terms, respectively The exact form of the source-terms can be written as J is the mean strain-rate ; + + B s can be seen, the turbulent transport, the velocity pressure-gradient and the dissipation rate must be modelled The viscous diffusion takes a simpler form, if the flow is incompressible or weakly compressible [6] J 216 Shima s Model for the Velocity Pressure-Gradient and the Dissipation Rate Term In this study, the modelling is performed utilizing Shima s approach [7] However, the dissipation equation is adapted from Chien s model [10], and also some modifications were needed in the velocity pressure-gradient term in order to produce reasonable results In Shima s model the velocity pressure-gradient and the dissipation rate term are connected as where ; 218 ;
10 ; 9 ; ; ; 219 ; 220 where ; is the production term of Eq210, and ; are 0 0 varies with the wall function as The multiplyer in Eq224 was, after test calculation, modified to be It should also be noted that in this study the tangential components of the velocity pressure-gradient wall term were omitted in order to get a correct velocity profile and shear-stress distribution close to the wall Diffusion Term In the diffusion term, the pressure diffusion is ignored or it is assumed to be part of the turbulent transport term In Shima s model, the diffusion term of Hanjalic and Launder s [11] is applied E Convergence problems were experienced using Eq225 as a consequence of the destabilizing cross-derivatives Because of this, two diffusion methods were tested The first one is a scalar diffusion term [12] B J 226 Diffusion in the normal direction is usually much more important than in the streamwise direction This can be taken into account using the model of Daly and Harlow s [13] + 0 bove turbulent eddy viscosity + 0 E 0 is defined as E 228
11 E where the kinetic energy of turbulence is The equation for contains empirical coefficients These are given by E E 0 23 Dissipation Transport Equation The dissipation transport equation was taken from Chien s $ model [10] because it was experienced to be stable and well behaved Chien s dissipation transport equation can be written as = E ; E where is the normal distance from the wall, and is defined by K 5 8 The production ; and is is a trace of the production tensor ; The relationship between where is kinematic viscosity With this assumption, has a value of zero at the wall and has an exact value with respect to However, does not have the desired behaviour near the wall This results in a non-zero diffusion flux at the wall The equation for contains empirical coefficients These are given by 234 where the turbulence Reynolds number is defined as?=@ 235 Chien proposed slightly different forms for E and E Since the computations performed for the flat plate boundary layer [8] appeared to be insensitive to the mod- ifications, the formulas above were based on the most commonly used coefficients E and E
12 J J J C C 11 3 umerical Method The finite-volume CFD program for complex three-dimensional geometries [8] was used in the present calculations The program utilizes Cartesian velocity components in a cell-centred approach In the evaluation of the inviscid fluxes, Roe s method [9] is applied For a spatial discretization MUSCL-type TVD-scheme to approximate advective volume-face fluxes is applied The discretized equations are integrated in time by applying the 3 -factorization [14] The code utilizes a multigrid V-cycle for the acceleration of convergence Complicated geometries can be handled with multiblock grids 31 Spatial Discretization In the present solution, a finite-volume technique is applied The flow equations have an integral form J 8 K J K > J 8 31 for an arbitrary fixed region 8 a computational cell L yields 8 with a boundary Performing the integrations for J 8 > 32 where the sum is taken over the faces of the computational cell The flux for the face is defined 33 Here, and are the fluxes defined by Eqs22 to 24 in the, and ' directions, respectively In the evaluation of the inviscid fluxes, Roe s method [9] is applied rotation operator is used for the velocity components and also for the Reynolds stresses The flux is calculated as 34 where C is a rotation operator that transforms the dependent variables to a local coordinate system normal to the cell surface In this way, only the Cartesian form
13 of the inviscid flux is needed This is calculated from < 35 where and are the solution vectors evaluated on the left and right sides of the cell surface, < is a the right eigenvector of the Jacobian matrix 1, the corresponding eigenvalue is, and is the corresponding characteristic variable obtained from?, where MUSCL-type discretization is used for the evaluation of and In the evaluation of and, primary flow variables 1, and conservative turbulent variables are utilized Jacobian matrix can be split in the following way?? 4 36 where? and? are the right and left eigenvector matrices of, and are the corresponding matrices with respect to the primitive variables, is the diagonal eigenvalue matrix, and 4 and 4 are the transformation matrices between the conservative and the primitive variables 32 Diagonalization of the Flow Equations coupling between the avier Stokes and the Reynolds-stress equations is introduced, L since the Reynolds stresses may be connected with the pressure [6] In the -momentum equation, the resulting effective pressure can be defined as + 37 In order to utilize Roe s method, the Jacobian of the flux vectors must be diagonalized This requires that the Jacobian matrix of the flux vector has a complete set of eigenvectors Unfortunately, linearly independent eigenvectors cannot be found if the anisotropic pressure field of Eq37 is applied Since the anisotropic pressure field is difficult to handle, the turbulent pressure is usually approximated by the mean of three components Using this, the flux vector can be divided into the isotropic and anisotropic parts, and the Jacobian of the isotropic part can be diagonalized The second method of diagonalization utilizes the production of the turbulence ; in the vector The production term is exact in RSM and it can be included in vector This way, independent eigenvectors can be found In the following, both of the diagonalization approaches are described
14 @ > ; 13 Isotropic Diagonalization Vector can be divided into two parts where vector corresponds to the isotropic part of turbulence and anisotropic part contains the 310 The Jacobian of the vector can be diagonalized similarly to that of the model [8] The effect of on a solution is small, and, consequently, can be evaluated using central differences In this approach, there is no need to rotate Reynolds stresses into the local cell face coordinates nisotropic Diagonalization The second method of diagonalization utilizes the production term ; The production term is exact in RSM and it can be separated from the other source-terms > Production is included in vector This is not a conservative form of the vector but RSM is never in a strong conservation law form because of the source-terms In many cases the Jacobian takes the simplest form if the primary variables are used Here, the selected primary variables are ote that bars and tildes are dropped out of the flow variables for simplicity fter some mathematical manipulation, the Jacobian 1 can be found to be in the
15 0 F E 14 primitive variable form as F F I H 1 In this case the eigenvalues, ie the characteristic speeds, are E where E is the speed of sound For an arbitrary equation of state, the speed of sound is E ++,+ otations + and + are % It is seen that the Reynolds stresses have an effect on the characteristic speeds and on the definition of the speed of sound Using the primitive variables, the characteristic variables are D?
16 $ $ $ $ $ $, + ' $ $ $ $ $ $ ' ' ' 15 where % " % " " and matrix " * % $ * 0 $ $ $ " " " " " "- "- " is * " * " " " 0 0 ' 318
17 16 The right eigenvector matrix is " 0 0 * * 0 * * * * * 0 319
18 E E 17 where < Ë 5 E < Ë E E E * 320 In matrices and? there are terms that have to be limited to avoid unnatural behaviour between turbulent and laminar regions For example, the term or cannot get very large values This can be limited by using following inequality using this 2 can be limited as Proof for Eq321 can be obtained in the following way: In a 2-dimensional case fluctuating velocities and can be rotated into an other coordinate system and, which make an angle with the axes and The Reynolds stresses can be obtained in the new coordinate system as If we choose so that 0 0 so-called principal axes and after some calculation Eq321 is reduced to the following form 325 which proves the inequality The eigenvectors and the characteristic variables have a fairly complex form In a computational approach this form can be rearranged and simplified to some degree
19 9 3 B 9 B Rotation perator In the anisotropic case, the Reynolds stresses as well as velocity components, must be rotated into a local coordinate system For this purpose a rotation operator has an effect on the velocity components and also on the Reynolds stresses if the anisotropic diagonalization is applied The normal of the cell face is known and two tangent directions must be determined There are many ways to do this but all of them have singularity points In the present calculations the singular direction of the normal component is The grid must be checked to be such that no singularities exist, ie, the normal vectors of the surface are not parallel to the vector rotation matrix is built by taking a vector product between the normal vector and the vector The resulting vector is normalized The third directional vector is obtained by a vector product of the two first ones The resulting rotation matrix B can be written as where and Velocity components are rotated from the global to local coordinate system as 329 where the hats denote local coordinates The Reynolds stresses are rotated from the global to the local coordinate system by the following formula [15] 3 B:3 B 330 where 331
20 3 B B ; 19 Transformation from the local to the global coordinate system is given by B 3 B 332 The rotation operator C is built by using this rotation matrix B for the velocities and the Reynolds stresses 34 Boundary Conditions The boundary values are given in ghost cells so that the actual boundary conditions are satisfied on the cell faces t the free-stream boundary, the values of the dependent variables are kept as constants In the flow field, and are limited to their free-stream values t the wall, the velocity components are set to zero In the present cases, the wall is also assumed to be adiabatic ll turbulent quantities are set to zero, also Because of this, these variables within the ghost cells are set to be of the opposite sign to the values in the cell adjacent to the surface The ghost cell values are applied for the calculation of the flux adjacent to the surface For the calculation of the surface fluxes themselves a second-order extrapolation is applied for the evaluation of the wall pressure and one-sided formulas are used for the derivatives at the wall Velocities at the symmetry walls can be calculated by first rotating the velocity components into the local Cartesian coordinate system " that has a coordinate direction normal to the boundary face The sign of the velocity component in this direction is changed and the three velocity components are then rotated back into the original coordinate system This can be expressed in a mathematical form as B 5 B ; fter some manipulations the following equation is obtained: where vector is the normal vector of the symmetry plane The resulting velocity vector 1 is put into the ghost cell This formulation for the velocity components was also applied by Batina [16] In general, the symmetry conditions of the stress tensor are very complicated formulations In the local coordinate system where the L direction is normal to the symmetry plain, this condition can be written ; 335
21 3 3 3 B 3 > 20 3 B 3 B B ; 3 ; B B ; B 3 B ; B ; 3 By applying Eqs330, 332 and 335, the symmetry formulation in the global coordinate system is 336 It is regonized that in this equation is the same as in Eq333 Therefore, the final form of the symmetry formula for the stress tensor is given by Time Integration Method The discretized equations are integrated in time by applying the 3 -factorization [14] This is based on the approximate factorization and on the splitting of the Jacobians of the flux terms In the implicit stage the factorization is done isotropicly The implicit stage consists of a backward and forward sweep in every coordinate direction The sweeps are based on a first-order upwind differencing In addition, the linearization of the source-term is factored out of the spatial sweeps The boundary conditions are treated explicitly, and a spatially varying time step is utilized The implicit stage can be written after factorization as follows ? 338 where and are first-order spatial difference operators in the L, M and directions,, and are the corresponding Jacobian matrices,, and? is the right-hand side of Eq32 The Jacobians are calculated as? 3? 339 where are diagonal matrices containing the positive and negative eigenvalues, and is a factor to ensure the stability of the viscous term [17] $ J 340 where J is the height of the cell The idea of the diagonally dominant factorization is to put as much weight on the diagonal as possible In the L direction the tridiagonal
22 > > 21 equation set resulting from Eq338 is replaced by two bidiagonal sweeps and a matrix multiplication [8] The matrix inversion resulting from the source-term linearization is performed before the spatial sweeps In order to improve stability, only negative source-terms can be linearized lthough the form of the source-term indicates that equations may become stiff near the walls, the terms related to the walls are not linearized Beacause of the complexity of the source-terms, the matrix is approximated by using the following trick 341 In this way, the maximum change of caused by > is limited to value of is evaluated using the current values of as The 342 was set to after test calculations Calculation of the and is described in [8] multigrid method is used to accelerate the convergence The Jameson s method [18] with a simple V-cycle has been adopted The implementation for the multigrid cycling is described in [19] and [8] In order to enhance the stability of the multigrid cycling, the size of correction from a coarse to a finer grid level is recalculated using the current value of as 343 fter test calculations, and were assigned to be and, respectively Variable assures a possible change of sign in the tangential components of Reynolds stresses
23 22 4 Test Calculations 41 Channel Flow The model was checked by calculating a fully developed flow in a plane channel The results were compared with the DS data of Kim et al [20], and the Reynoldsstress budgets were compared with Mansour et al [21] data The DS data is at?=@ 5 where, and are the mean velocity, the channel halfwidth and molecular viscosity Because the flow solver uses compressible methods the Mach number was set to This introduced a change in density across the channel The mesh is rectangular The height of the first row of cells is or nly half of the channel is modelled The length of the computational mesh is The calculations were performed using cyclic boundary conditions fter having a converged result, the solution was taken from the downstream boundary and utilized as the upstream boundary condition of the next run Fully developed flow was obtained after computations, which corresponds to the length of channel widths The convergence of the results was checked by using a two-times denser grid The results obtained with the two grid densities are practically identical Several solution methods described earlier were compared The numbering of the test cases can be seen in Table 41 The calculated mean flow variables are presented in Table 42 s can be seen in Table 42, the difference between the isotropic and anisotropic flux-difference splittings, case and case, is small The displacement thickness and the momentum thickness are defined as J 41 Table 41 Description of the test cases Case Multigrid Flux splitting Diffusion 1 nisotropic Daly et al 2 nisotropic Scalar diffusion 3 Isotropic Scalar diffusion 4 Isotropic Daly et al
24 23 Table 42 Mean flow variables DS Case 1 Case 2 Case 3 E?=@ 5?=@?=@ =5 Fig 41 Convergence of the norm of the momentum and the J residuals 42 where is a half-width of the channel and velocity at the centre line The values of and are approximately smaller than in DS The calculations were performed with five grid levels at The multigrid corrections were omitted on the first two cells from the solid wall The secondorder upwind scheme was used First, iteration cycles were performed with the model at and the Reynolds stresses were uncoupled from the flow converged solution was obtained after to iteration cycles The convergence history of the norm of the momentum and the residuals are shown in Fig 41 In this case the difference in convergence rate between the isotropic and anisotropic methods is marginal The anisotropic method converged roughly in cycles faster than the isotropic one The computing time with RSM is increased by if the anisotropic coupling is applied The increased convergence rate was obtained by the multigrid acceleration, as can be seen in Fig41 Without the multigrid cycling it takes about iteration sweeps to get a converged solution The convergence rate is about the same with Chien s model and with RSM aturally, in the calculations with Chien s model there is no transient after cycles Chiens model runs two times faster
25 24 Fig 42 Mean velocity profiles in wall coordinates Fig 43 Comparision of the calculated Reynolds stresses and the DS-data in a plane channel per iteration sweep than RSM The velocity profiles are compared in Fig 42 in terms of, which is a universal dimensionless velocity defined as, where is a friction velocity The velocity profiles in a viscous sublayer agree well with DS and universal profiles The velocity profiles are not completely satisfactory in outer layers The Reynolds stresses can be seen in Fig43 where Turbulent intensities agree well with DS data except that the peak level is low and the near-wall values are not satisfactory The fluctuating component normal to the wall should be damped rapidly close to the wall, but the simulation entirely misses this effect This has also an effect on the source-term distributions The shear-stress distribution in Fig 43 agrees with the DS results The source-term distributions are presented in Fig 44 The source terms are
26 25 Fig 44 Budgets of Reynolds stresses Symbols are from the DS calculation non-dimensionalized using, where is kinematic viscosity Close to the wall, the source terms do not agree with the DS data, but closer to the centre of the channel the agreement is good The velocity pressure-gradient term does not behave correctly, and there are also problems in the dissipation term close to the wall The wall correction in the velocity pressure-gradient was omitted because it gave a totally unsatisfactory distribution of the shear stress Production of the shear stress ; is not correct close to the wall This is a consequence of the normal Reynolds-stress component perpendicular to the wall having too high values in the viscous sublayer in Fig 43 The dissipation term could have been modelled differently by using a Taylor series near the wall [22] That approach could have damped near the wall more efficiently lthough not shown, anisotropic and isotropic flux-difference splittings do not introduce changes in the distributions of the source-terms The calculation method for diffusion does not have a strong effect on the convergence rate The convergence rate of the scalar diffusion case case 2 in Fig 41 is similar to that obtained using the Daly and Harlow method case 1 However, the mean flow variables are changed slightly in Table 42 Especially, the choice of the diffusion model changed the friction coefficient Some difference can also be seen
27 26 Fig 45 Comparision of the turbulent diffusion Symbols are from the DS calculation in the velocity profiles in Fig 42, and in the shear stresses in Fig 43 lthough there are differences in the velocity profiles between these methods, the turbulent diffusion B does not exhibit large differences in Fig 45 It can be seen from Fig 45 that the turbulent diffusion rates are not satisfactorily modelled The Daly and Harlow turbulent diffusion term of Eq227 is only slightly better than the scalar diffusion term of Eq226
28 27 5 Summary The Reynolds-averaged avier Stokes equations with a low-reynolds number RSM have been solved using an implicit method with a multigrid acceleration for convergence In the evaluation of fluxes the turbulence equations are coupled with the inviscid part of the flow equations, and Roe s method is applied new anisotropic coupling of the avier Stokes and the Reynolds-stress equations is introduced lso a new method of treating symmetry boundaries is presented The solution methods have been tested using Shima s low-reynolds number model The developed numerical scheme appears to be stable and efficient Per iteration cycle the calculation with the Reynolds-stress model takes only about twice as long as the calculation with Chien s low-reynolds number model Using a multigrid only a few hundred iteration cycles are required for the solution of a flow in a plane channel This paper has focused attention on the problem of coupling the Reynolds stresses with avier Stokes equation with a compressible flow assumption The applied closure model is a relatively old one Hence, a very good agreement with the DS data was not even expected The new coupling method introduced only small improvements in the convergence rate in the present incompressible case The differences between the coupling methods may become larger in a case of super or hypersonic flow
29 28 Bibliography [1] Siikonen, T, Kaurinkoski, P, and Laine, S, Transonic Flow over a Delta Wing Using a Turbulence Model, in Proceedings of the 19th ICS Congress, naheim, pp , Sept 1994 ICS Paper [2] Baldwin, B and Lomax, H, Thin Layer pproximation and lgebraic Model for Separated Turbulent Flows, Jan 1978 I Paper [3] Holst, T, Viscous Transonic irfoil Workshop Compendium of Results, Journal of ircraft, Vol 25, Dec 1988, pp [4] Leschziner, M, Modelling Engineering Flows with Reynolds Stress Turbulence Closure, Journal of Wind Engineering and Industrial erodynamics, Vol 35, 1990, pp [5] Davidson, L, Prediction of the Flow round an irfoil Using a Reynolds Stress Transport Model, Journal of Fluids Engineering, Vol 117, March 1995, pp [6] Vandromme, D, Turbulence Modeling for Turbulent Flows and Implementation in avier Stokes Solvers, in Introduction to the Modeling of Turbulence, von Karman Institute for Fluid Dynamics Lecture Series , 1991 [7] Shima,, Reynolds-Stress Model for ear-wall and Low-Reynolds- umber Regions, Journal of Fluids Engineering, Vol 110, 1988, pp [8] Siikonen, T, n pplication of Roe s Flux Difference Splitting for the Turbulence Model, Helsinki University of Technology, Laboratory of erodynamics, 1994 ISB [9] Roe, P, pproximate Riemann Solvers, Parameter Vectors, and Difference Schemes, Journal of Computational Physics, Vol 43, 1981, pp [10] Chien, K-Y, Predictions of Channel and Boundary-Layer Flows with a Low-Reynods-umber Turbulence Model, I Journal, Vol 20, Jan 1982, pp [11] Hanjalic, K and Launder, B E, Reynolds Stress Model of Turbulence and its pplication to Thin Shear Flows, Journal of Fluid Mechanics, Vol 52, o 4, 1972, pp
30 29 [12] Davidson, L and Rizzi,, avier Stokes stall predictions using an lbebraic Reynolds-Stress model, Journal of Spacecraft and Rockets, Vol 29, o 6, 1992, pp [13] Daly, B and Harlow, F, Transport Equations of Turbulence, Physics of Fluids, Vol 13, 1970, pp [14] Lombard, C, Bardina, J, Venkatapathy, E, and liger, J, Multi- Dimensional Formulation of CSCM n Upwind Flux Difference Eigenvector Split Method for the Compressible avier Stokes Equations, in 6th I Computational Fluid Dynamics Conference, Danvers, Massachusetts, pp , Jul 1983 I Paper CP [15] Salonen, E-M, Dynamiikka Espoo: takustantamo, 1980 ISB [16] Batina, J, ccuracy of an Unstructured-Grid Upwind-Euler lgorithm for the ER M6 Wing, Journal of ircraft, Vol 28, June 1991, pp [17] bayashi, S and Kuwahara, K, n pproximate LU Factorization Method for the Compressible avier Stokes Equations, Journal of Computational Physics, Vol 63, March 1986, pp [18] Jameson, and Yoon, S, Multigrid Solution of the Euler Equations Using Implicit Schemes, I Journal, Vol 24, o 11, 1986 [19] Siikonen, T, Hoffren, J, and Laine, S, Multigrid 7 Factorization Scheme for the Thin-Layer avier Stokes Equations, in Proceedings of the 17th ICS Congress, Stockholm, pp , Sept 1990 ICS Paper [20] Kim, J, Moin, P, and Moser, R, Turbulence Statistics in Fully Developed Channel Flow at Low Reynolds umber, Journal of Fluid Mechanics, Vol 177, 1987, pp [21] Mansour,, Kim, J, and Moin, P, Reynolds-stress and Dissipation-rate Budgets in a Turbulent Channel Flow, Journal of Fluid Mechanics, Vol 194, 1988, pp [22] Launder, B E and Tselepidakis, D, Directions in Second-Moment Modelling of ear-wall Turbulence, in 29th erospace Sciences Meeting, evada, Reno, Jan 1991
Turbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
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 informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
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 informationAalto University School of Science and Technology CFD-group/ Department of Applied Mechanics. MEMO No CFD/MECHA DATE: March 15, 2012
Aalto University School of Science and Technology CFD-group/ Department of Applied Mechanics MEMO No CFD/MECHA-2-212 DATE: March 15, 212 TITLE The Effect of Free-Stream Turbulence Parameters on the SST
More informationFluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 49 Introduction to Turbulent Flow part -II Good morning I welcome you all to this
More informationWall treatments and wall functions
Wall treatments and wall functions A wall treatment is the set of near-wall modelling assumptions for each turbulence model. Three types of wall treatment are provided in FLUENT, although all three might
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationHelsinki University of Technology Laboratory of Applied Thermodynamics/CFD-Group. MEMO No CFD/TERMO DATE: November 21, 1997
Helsinki University of Technology Laboratory of Applied Thermodynamics/CFD-Group MEMO No CFD/TERMO-22-97 DATE: November 21, 1997 TITLE GRV-60 Impeller with airfoil-type and unvaned diffusor AUTHOR(S) Harri
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationCOMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL
COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the
More informationTurbulence - Theory and Modelling GROUP-STUDIES:
Lund Institute of Technology Department of Energy Sciences Division of Fluid Mechanics Robert Szasz, tel 046-0480 Johan Revstedt, tel 046-43 0 Turbulence - Theory and Modelling GROUP-STUDIES: Turbulence
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationNumerical Simulation of Flow Around An Elliptical Cylinder at High Reynolds Numbers
International Journal of Fluids Engineering. ISSN 0974-3138 Volume 5, Number 1 (2013), pp. 29-37 International Research Publication House http://www.irphouse.com Numerical Simulation of Flow Around An
More informationDAY 19: Boundary Layer
DAY 19: Boundary Layer flat plate : let us neglect the shape of the leading edge for now flat plate boundary layer: in blue we highlight the region of the flow where velocity is influenced by the presence
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationSIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS
ICAS 2000 CONGRESS SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS H Yao, R K Cooper, and S Raghunathan School of Aeronautical Engineering The Queen s University of Belfast, Belfast BT7 1NN,
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationWall turbulence with arbitrary mean velocity profiles
Center for Turbulence Research Annual Research Briefs 7 Wall turbulence with arbitrary mean velocity profiles By J. Jiménez. Motivation The original motivation for this work was an attempt to shorten the
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
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 informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More information1. Introduction, tensors, kinematics
1. Introduction, tensors, kinematics Content: Introduction to fluids, Cartesian tensors, vector algebra using tensor notation, operators in tensor form, Eulerian and Lagrangian description of scalar and
More informationNumerical investigation of swirl flow inside a supersonic nozzle
Advances in Fluid Mechanics IX 131 Numerical investigation of swirl flow inside a supersonic nozzle E. Eslamian, H. Shirvani & A. Shirvani Faculty of Science and Technology, Anglia Ruskin University, UK
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
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 informationThe mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):
8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation
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 informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More information4.2 Concepts of the Boundary Layer Theory
Advanced Heat by Amir Faghri, Yuwen Zhang, and John R. Howell 4.2 Concepts of the Boundary Layer Theory It is difficult to solve the complete viscous flow fluid around a body unless the geometry is very
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 informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
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 informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationCOMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE
COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE L. Velázquez-Araque 1 and J. Nožička 2 1 Division of Thermal fluids, Department of Mechanical Engineering, National University
More informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationTable of Contents. Foreword... xiii. Preface... xv
Table of Contents Foreword.... xiii Preface... xv Chapter 1. Fundamental Equations, Dimensionless Numbers... 1 1.1. Fundamental equations... 1 1.1.1. Local equations... 1 1.1.2. Integral conservation equations...
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationBOUNDARY LAYER FLOWS HINCHEY
BOUNDARY LAYER FLOWS HINCHEY BOUNDARY LAYER PHENOMENA When a body moves through a viscous fluid, the fluid at its surface moves with it. It does not slip over the surface. When a body moves at high speed,
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 informationcompression corner flows with high deflection angle, for example, the method cannot predict the location
4nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 4, Reno, Nevada Modeling the effect of shock unsteadiness in shock-wave/ turbulent boundary layer interactions AIAA 4-9 Krishnendu Sinha*, Krishnan
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationModeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique
Modeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique Torsten Palenschat 317220 30.04.2014 Master seminar CES WS 2013/2014 Center for Computational Engineering Sciences RWTH Aachen
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.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 informationProblem 4.3. Problem 4.4
Problem 4.3 Problem 4.4 Problem 4.5 Problem 4.6 Problem 4.7 This is forced convection flow over a streamlined body. Viscous (velocity) boundary layer approximations can be made if the Reynolds number Re
More informationHeat Transfer from An Impingement Jet onto A Heated Half-Prolate Spheroid Attached to A Heated Flat Plate
1 nd International Conference on Environment and Industrial Innovation IPCBEE vol.35 (1) (1) IACSIT Press, Singapore Heat Transfer from An Impingement Jet onto A Heated Half-Prolate Spheroid Attached to
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationEffect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process
, pp. 844 850 Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process Baokuan LI and Fumitaka TSUKIHASHI 1) Department of Thermal Engineering, The School
More informationAnisotropic grid-based formulas. for subgrid-scale models. By G.-H. Cottet 1 AND A. A. Wray
Center for Turbulence Research Annual Research Briefs 1997 113 Anisotropic grid-based formulas for subgrid-scale models By G.-H. Cottet 1 AND A. A. Wray 1. Motivations and objectives Anisotropic subgrid-scale
More informationAA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations
AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental
More informationFluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay. Lecture - 17 Laminar and Turbulent flows
Fluid Mechanics Prof. T.I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay Lecture - 17 Laminar and Turbulent flows Welcome back to the video course on fluid mechanics. In
More informationfluctuations based on the resolved mean flow
Temperature Fluctuation Scaling in Reacting Boundary Layers M. Pino Martín CTR/NASA Ames, Moffett Field, CA 94035 Graham V. Candler Aerospace Engineering and Mechanics University of Minnesota, Minneapolis,
More informationENGINEERING MECHANICS 2012 pp Svratka, Czech Republic, May 14 17, 2012 Paper #195
. 18 m 2012 th International Conference ENGINEERING MECHANICS 2012 pp. 309 315 Svratka, Czech Republic, May 14 17, 2012 Paper #195 NUMERICAL SIMULATION OF TRANSITIONAL FLOWS WITH LAMINAR KINETIC ENERGY
More informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
More informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
More informationNUMERICAL METHOD FOR THREE DIMENSIONAL STEADY-STATE TWO-PHASE FLOW CALCULATIONS
' ( '- /A NUMERCAL METHOD FOR THREE DMENSONAL STEADY-STATE TWO-PHASE FLOW CALCULATONS P. Raymond,. Toumi (CEA) CE-Saclay DMT/SERMA 91191 Gif sur Yvette, FRANCE Tel: (331) 69 08 26 21 / Fax: 69 08 23 81
More informationFundamentals of Aerodynamics
Fundamentals of Aerodynamics Fourth Edition John D. Anderson, Jr. Curator of Aerodynamics National Air and Space Museum Smithsonian Institution and Professor Emeritus University of Maryland Me Graw Hill
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II External Flows 1 Introduction In this chapter we will consider several fundamental flows, namely: the flat plate, the cylinder, the sphere, several other body shapes, and banks
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 informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationReview and Assessment of Turbulence Transition Models
International Journal of Engineering Research and Development e-issn: 2278-67X, p-issn: 2278-8X, www.ijerd.com Volume 13, Issue 43 (April 217), PP.32-57 Review and Assessment of Turbulence Transition Models
More informationComparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2
Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2 Department of Mechanical and Materials Engineering, Military Engineering Institute,
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationMODELLING OF INFLUENCE OF TURBULENT TRANSITION ON HEAT TRANSFER CONDITIONS KRZYSZTOF BOCHON, WŁODZIMIERZ WRÓBLEWSKI
TASK QUARTERLY 12 No 3, 173 184 MODELLING OF INFLUENCE OF TURBULENT TRANSITION ON HEAT TRANSFER CONDITIONS KRZYSZTOF BOCHON, WŁODZIMIERZ WRÓBLEWSKI AND SŁAWOMIR DYKAS Institute of Power Engineering and
More information41st Aerospace Sciences Meeting and Exhibit 6-9 January 2003, Reno, Nevada Modeling shock unsteadiness in shock/turbulence interaction
4st Aerospace Sciences Meeting and Exhibit 6-9 January 003, Reno, Nevada Modeling shock unsteadiness in shock/turbulence interaction AIAA 003-65 Krishnendu Sinha*, Krishnan Mahesh and Graham V. Candler
More informationAPPLICATION OF THE DEFECT FORMULATION TO THE INCOMPRESSIBLE TURBULENT BOUNDARY LAYER
APPLICATION OF THE DEFECT FORMULATION TO THE INCOMPRESSIBLE TURBULENT BOUNDARY LAYER O. ROUZAUD ONERA OAt1 29 avenue de la Division Leclerc - B.P. 72 92322 CHATILLON Cedex - France AND B. AUPOIX AND J.-PH.
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 information(Refer Slide Time: 2:14)
Fluid Dynamics And Turbo Machines. Professor Dr Shamit Bakshi. Department Of Mechanical Engineering. Indian Institute Of Technology Madras. Part A. Module-1. Lecture-3. Introduction To Fluid Flow. (Refer
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 informationA Solution Method for the Reynolds-Averaged Navier-Stokes Equation
A Solution Method for the Reynolds-Averaged Navier-Stokes Equation T.-W. Lee Mechanical and Aerospace Engineering, Arizona State University, Tempe, AZ, 8587 Abstract- Navier-Stokes equations are difficult
More informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationExplicit algebraic Reynolds stress models for boundary layer flows
1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both
More informationPEMP ACD2505. M.S. Ramaiah School of Advanced Studies, Bengaluru
Governing Equations of Fluid Flow Session delivered by: M. Sivapragasam 1 Session Objectives -- At the end of this session the delegate would have understood The principle of conservation laws Different
More information