n i,j+1/2 q i,j * qi+1,j * S i+1/2,j
|
|
- Cody Copeland
- 5 years ago
- Views:
Transcription
1 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 ow and lid-driven cavity ow. AUTHOR() Petri Majander ABTRACT Viscous terms are derived in curved co-ordinates using nite volume formulation. An approximation in common use called Thin hear Layer approximation is derived. Turbulent ow over a at plate with a k ; -model as well as a laminar lid-driven ow are computed and the results are compared between the full viscous terms and the TL approximation. MAIN REULT An implementation of full viscous terms in FINFLO solver. PAGE 4 KEY WORD Thin hear Layer, viscous ow, boundary layer ow, lid-driven cavity ow APPROVED BY Timo iikonen December 9,997
2 Introduction At high Reynolds numbers wall shear layers, wakes or free shear layers will be of limited size, and if the viscous region remains limited then the dominating inuence of the shear layers comes from the gradients transverse to the main ow direction. If we consider an arbitrary curvilinear system of co-ordinates with and along the surface and directed towards the normal of the surface, then the thin shear layer (TL) approximation of the Navier-tokes equations consists of neglecting all and -derivatives occurring in the turbulent and viscous shear stress terms. Also, at high Reynolds numbers the mesh is made dense in the direction normal to the shear layer, and therefore the neglected terms are computed with a lower accuracy []. In the following applications derivatives are computed in a simplied manner in all curved co-ordinade directions which is referred to as TL approximation. Finite Volume Formulation. Derivation of the viscous terms Transport equation in conservation form can be written as d dt V UdV + which in discretized form is written as X f aces d dt (V i iu i )+ The ux at the surface is written as k ~F (U) d~ V QdV () ^F (U) k V i Q i : () ^F (F ; F v )~n x +(G ; G v )~n y +(H ; H v )~n z () where the components of the unit normal vector are written as ~n x jj~i, x ~n y ~j y and ~n jj z jj~ z k. Each surface of the cell has a unit normal ~n (see Fig. ). Viscous uxes F v, G v and H v are " # T F v 0 xx xy xz u xx + v xy + w xz ; q x x " # T G v 0 yx yy yz u yx + v yy + w yz ; q y y " # T H v 0 zx zy zz u xz + v yz + w zz ; q z (4) z
3 n i,j+/ q i,j * qi+,j * i+/,j i,j / i+/,j Fig. : Two neighbouring gridcells in the -direction are shown to clarify the notations. where the last term is viscous part of generic scalar equation, including turbulent kinetic energy and dissipation of kinetic energy. The diusion coecient is approximated as + T, where T is the turbulent viscosity and is appropriate chmidt's number. Heat conduction ux is ~q ;(k + k T )rt ; c p Pr + c p T rt (5) Pr T where k, k T are laminar and turbulent thermal conductivities respectivily, Pr,Pr T are laminar and turbulent Prandl numbers respectivily and c p is specic heat at constant pressure. A stress tensor for newtonian viscous uid is written as ij u j + u i ; x i x j (r ~V ) ij ; u 0 iu 0 j (6) where for the Reynolds stresses Boussinesq's approximation is used u ;u 0 iu 0 j j T + u i ; x i x j (r ~V ) ij The stress tensor is symmetric ij ji. Viscous uxes are obtained from Eqs. () and (4) ^F v 0 ^F v ^F v ^F 4 v ^F 5 v ( xx n x + yx n y + zx n z )d ( xy n x + yy n y + zy n z )d ( xz n x + yz n y + zz n z )d [(u xx + v xy + w xz )n x +(u xy + v yy + w yz )n y (7)
4 + (u xz + v yz + w zz )n z ; q x n x ; q y n y ; q z n z ]d ^F 6 v [u( xx n x + xy n y + xz n z )+v( xy n x + yy n y + yz n z ) + w( zx n x + zy n y + zz n z ) ; q x n x ; q y n y ; q z n z ]d x n x + y n y + z n z d: (8) After inserting the stress tensors dened by Eqs. (6) and (7) the uxes in the direction and on the surface i+ become " ^F ( + v i+ T ) u v w x ; y ; n z x i+ + u y + v n x y i+ w + x + u # n z z i+ i+ ^F v i+ ^F 4 v i+ " u ( + T ) y + v x v + z + w n y z i+ ( + T ) " u z + w x + w u v z ; x ; y n x i+ + v u w y ; x ; z # i+ n x i+ + v n z i+ # i+ z + w y n y i+ n y i+ ^F 5 v i+ ^F 6 v i+ u i+ ; qx i+ ^F v i+ + v i+ n x i+ ^F v i+ + w i+ + q y i+ n y i+ ^F 4 v i+ + q z i+ x n x i+ + y n y i+ + z n z i+ n z i+ i+ i+ (9) These equations are valid on every surface as corresbonding components of the surface areas l and normals n l are used. A detail concerning a computation of incompressible ow is mentioned here. The stress tensor of Eq. (6) can be divided into incompressible and compressible parts u j ij ( + T ) + u i ; x i x j (r ~V ) ij ij ; ( + T ) (r ~V ) ij : (0) Inserting this form to the uxes of Eq. (8) gives for ^F v ^F v xx ; ( + T ) (r ~V ) ^F 5 v u xx ; ( + T ) (r ~V ) + v xy n x + yy ; ( + T ) (r ~V ) and ^F 5 v n x + yx n y + zx n z d n x + xy n y + xz n z n y + yz n z
5 + w zx n x + zy n y + zz ; ( + T ) (r ~V ) n z ; q x n x ; q y n y ; q z n z ] d () For an incompressible ow r ~ V 0 mathematically. In conservative equations the mass ux is in balance when R ~Vd 0. In a pseudo-compressibility method density is not necessarily constant. Even with a constant density, the mass ux balance does not indicate that r~v 0, since the uxes are numerically calculated from dierent formulae from r~v 0. Hence compressible viscous eects may be present to some extent. For truly incompressible i.e. constant density ows this eect is of the same order as the truncation order. In order to compute the cartesian derivatives a somewhat less known version of the generalized divergence theorem of Gauss is applied (see e.g. []), V rdv d~ () Integration is performed over the volume V ijk on the left side, V ijk rdv V ijk x ~ i + y ~ j + z ~ k and over the corresponding faces on the right side, d ~ x i+ j k i+ j k ; x i; j k i; j k ~ i + x i j+ x i j k+ y i+ y i j+ y i j k+ z i+ z i j+ ijk k i j+ ; k x i j; k i j; ~ i + k ij k+ ; x i j k; i j k; j k i+ j k ; y i; j k i; 4 () ~ i + j k ~ j + k i j+ ; k y i j; k i j; ~ j + k i j k+ ; y i j k; i j k; j k i+ j k ; z i; j k i; z i j k+ i j k+ k i j+ k ; z i j; k i j; ; z i j k; i j k; Hence the derivative of in x -direction can be solved V ijk x ijk ~ j + j k ~k + ~k + k ~k: (4) x i+ j k i+ j k ; x i; j k i; j k + x i j+ k i j+ k ; x i j; k i j; + k x i j k+ ij k+ ; x i j k; i j k; (5) If the variables are dened in the middle of volumes, then the values at surfaces must be interpolated as averages. These integrations can also be performed over a dierent volume e.g. V i+ j k which is equivalent to increasing all i -indeces by in Eqs. (4) and (5).
6 5. Thin hear Layer Approximation Computing derivatives accurately is computationally heavy and the required memory size is increased, which of course also depends on the structure of the code. At high Reynolds numbers only the terms in the direction of curved co-ordinates are preserved. This is obtained by dropping the cross derivatives from equation (5) []. The computational grid should be as orthogonal and dense as possible in the directions against the surfaces, where the dominating viscous eects are present. Also, in the following, the integration is performed over the volume V i+ j k, V i+ j k x i+ j k x i+ j k i+ j k ; x i j k i j k : (6) In the thin-layer approximation the surface components are replaced by the interpolated surface component between them. This is usually readily available and the volume V i+ j k can be approximated as an average (V i + V i+ ). The Eq.( 6) can be further modied as x j k x i+ ( i+ j k ; i j k ) V i+ j k i+ j k i+ j k V i+ j k n x i+ j k ( i+ j k ; i j k ) n x i+ j k ( i+ j k ; i j k ) j k : (7) As the approximate derivatives of Eq. (7) are inserted into the uxes given by Eqs.( 9) ^F v i+ ; 8 < : n z i+ d i+ + 4 n y i+ + 4 n x i+ 4 n x i+ 4 4 (w i+ ; w i ) (u i+ ; u i ) ; 5 nx i+ (u i+ ; u i )+ n x i+ n y i+ (v i+ ; v i ) (v i+ ; v i ) (w i+ ; w i )+ n z i+ (u i+ ; u i ) n x i+ + n y i+ + n z i+ + n x i+ n y i+ (v i+ ; v i ) 5 ny i+ 5 nz i+ (u i+ ; u i ) 9 i+
7 + n x i+ n z i+ (w i+ ; w i ) 5 i+ Other momentum ux components are derived in a similar manner, ^F v i+ ^F 4 v i+ + n y i+ 4 n x i+ (u i+ ; u i ) 4 n + n + n y i+ x i+ z i+ + n z i+ n y i+ (w i+ ; w i ) + 4 n x i+ n z i+ 5 i+ (v i+ ; v i ) (u i+ ; u i )+ n y i+ n z i+ (v i+ ; v i ) 4 n z i+ + n x i+ + n z i+ (w i+ ; w i ) 5 i+ 6 (8) : (9) The TL approximation for the energy and the generic scalar equation becomes ^F 5 v i+ u i+ ^F + v v i+ i+ ^F + w v i+ i+ ^F 4 v i+ ; (k + k T ) T x n x i+ + T y n y i+ + T z n z i+ i+ ^F 6 v i+ u i+ ^F + v v i+ i+ c p ; + Pr T Prdi+ Pr T ^F + w v i+ i+ ^F 4 v i+ i+ (T i+ ; T i ) i+ x n x i+ + y n y i+ + z n z i+ ( i+ ; i ) i+ i+ i+ i+ The momentum uxes can also be written in a compact form ^F i+ v i+ i+ (u i+ ; u i )+ n x i+ (u i+ ; u i ) ^F (v v i+ i+ ; v i )+ n y i+ (u i+ ; u i ) ^F 4 v i+ i+ i+ (w i+ ; w i )+ n z i+ (u i+ ; u i ) (0) () where u i and u i+ are the scaled contravariant velocities in the direction of the surface ~ i+ u i n x i+ u i + n y i+ v i + n z i+ w i ()
8 7 Fig. : Grid of the at plate boundary ow. In the implicit phase the equations are linearized in a simplied manner. The eect of the viscous terms T is added on the Jacobian matrices A FU, B GU and C HU as A R A ( A + TI)R ; A () B R B ( B + TI)R ; B C R C ( C + TI)R ; C : Here R is the transformation matrix from conservative variables U to characteristic variables W, W R ; U. are the diagonal matrices containing the positive and the negative eigenvalues of the Jacobians. The diagonal weight factor T is computed as T ( + T ) (4) d where d is the lenght of the cell in the direction of the computational sweep. Test cases. Flat Plate The case is a ow over at plate with a high free-stream turbulent intensity.inlet velocity was 9.4 m/s and the pressure gradient was zero. Measurements have been made down to x.495 m that corresponds to Re x Upstream turbulence intensity is Tu 6.0%.Dissipation is set so that decay of free stream is in balance with production. The calculation is started 6 cm before the plate. The lenght of the at plate is.6 m. The height is0cmand the height of the rst row of cells was :5 0 ;5 [4] that is equal to y + 0:7 at most of the domain (at the leading edge y + :). Grid is clustered near the wall though the nearest three rows are kept constant. The ratio between neighbouring cells is y n+ y n :5. The grid size is and the grid is seen in gure.
9 8 u w Fig. : Grid of the lid driven cavity ow.. Lid-driven cavity ow The second test case is an isothermal cavity ow -problem. The grid size for the problem is and the grid is seen in gure. The top wall is moving with velocity m/s and the resulting Reynolds number is Results 4. Flow over the Flat Plate The second-order upwind discretization with Roe's splitting was used for inviscid uxes. The case was scaled so that the free stream velocity U e was increased from 9:4ms to 94ms to ensure convergence. Computation with full uxes was done with four multigrid levels whereas the TL computation converged only with three levels maximum. These cases were iterated excessive 0000 iterations at the nest grid level, which wasduetopoorconvergence of some properties in the TL computation shown in Fig. (4). The computation of the full viscous uxes took about 5% more CPU -time/cycle and the memory need was.09 MB for the full terms and 9.08 MB for the TL approximation (9.47 MB with the same four levels). The results were nearly identical from both calculations. kin friction c f w U e, where w is a shear stress at the wall, is presented in Fig. 5. The velocity proles are compared in Fig. 6 in terms q of u +, which is a dimensionless velocity dened as u + uu, where u w is a friction velocity. Dimensionless distance
10 9 a) b) Fig. 4: Convergence histories of drag coecients, total mass and turbulent kinetic energy with full viscous terms (a) and with TL approximation (b). Fig. 5: kin friction coecient.
11 0 from the q wall is approximated as y + y n u. Turbulence level is dened as Tu ku e, where k is the turbulent kinetic energy. The turbulence level proles are presented in Fig. 7 In this case a pseudo-compressibility method was used with the second-order upwind discretization iterations were computed. The u/v-velocity distributions were plotted at y/x -positions 0:065L, 0:5L, 0:5L, 0:75L and 0:975L respectively, where L stands for the lenght and the height of the cube. The plots are represented in Fig. 8. If u-velocity results are compared with the reference solution [5] at y 0:5L, a suspection arised that numerical dissipation has spoilt the solution. The computatation with the full viscous terms was continued with 5000 iterations with the third-order upwind biased uxes and these results are closer to the reference solution. In this case the full viscous terms give results closer to the reference solution than the TL approximation. 5 Discussion In the boundary layer ow there was practically no dierence between the results computed with the TL approximation and the full viscous terms. In this case the dominating viscous forces might be those in the direction of the plate surface and hence in direction of the grid line. In the cavity ow there was a small dierence between the two cases and the full viscous terms produced slightly better results. This might result from the nature of the case, at small Reynolds numbers the TL approximation is not valid. If the last 'contravariant' terms are neglected in equations (), then correct terms are obtained for any orthogonal grid with incompressible ow. It was rst thought that this 'compressible' term might not vanish as well as rv -term in the full viscous terms, but both terms were of the same order of magnitude. Also, the pseudo-compressibility method was suspected to cause numerical dissipation in the solution.
12 Fig. 6: Dimensionless velocity proles at various downstream positions.
13 Fig. 7: Turbulence intensities at various downstream positions.
14 Fig. 8: caled u/v -velocity proles at dierent x/y -positions.
15 4 References [] Hirsch C., Numerical Computation of Internal and External Flows, Volume : Fundamentals of Numerical Discretization, John Wiley & ons Ltd. [] Arfken G., Mathematical Methods for Physicists, Academic Press Inc., 985, an Diego. [] iikonen T., Computation of Viscous Fluxes in Curved Co-ordinates., memo M-0-88, Laboratory of Aerodynamics, Helsinki University of Technology (in Finnish). [4] Rautaheimo P., Test Calculation of Various Two-equation Low Reynolds Number Turbulence Models,Memo CFD/TERMO-96, Laboratory of Applied Thermodynamics, Helsinki University of Technology. [5] Ghia, U. and Ghia, K.N. and hin, C.T.,High-Re olutions for Incompressible Flow Using the Navier-tokes Equations and a Multigrid Method, Journal of Computational Physics, 48, 98.
Chapter 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 informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
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 informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
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 informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. 12/21/01 topic7_ns_equations 1
AE/ME 339 Professor of Aerospace Engineering 12/21/01 topic7_ns_equations 1 Continuity equation Governing equation summary Non-conservation form D Dt. V 0.(2.29) Conservation form ( V ) 0...(2.33) t 12/21/01
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationk-" Turbulence Model to Attached and
NASA Technical Paper 348 Application of Navier-Stokes Code PAB3D With k-" Turbulence Model to Attached and Separated Flows Khaled S. Abdol-Hamid, B. Lakshmanan, and John R. Carlson January 1995 NASA Technical
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 informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationSpatial discretization scheme for incompressible viscous flows
Spatial discretization scheme for incompressible viscous flows N. Kumar Supervisors: J.H.M. ten Thije Boonkkamp and B. Koren CASA-day 2015 1/29 Challenges in CFD Accuracy a primary concern with all CFD
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationAE/ME 339. Computational Fluid Dynamics (CFD) K. M. Isaac. Momentum equation. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 9//005 Topic7_NS_ F0 1 Momentum equation 9//005 Topic7_NS_ F0 1 Consider the moving fluid element model shown in Figure.b Basis is Newton s nd Law which says
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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information12. Stresses and Strains
12. Stresses and Strains Finite Element Method Differential Equation Weak Formulation Approximating Functions Weighted Residuals FEM - Formulation Classification of Problems Scalar Vector 1-D T(x) u(x)
More informationGoverning Equations of Fluid Dynamics
Chapter 3 Governing Equations of Fluid Dynamics The starting point of any numerical simulation are the governing equations of the physics of the problem to be solved. In this chapter, we first present
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 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 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 informationIntroduction For a long time, so-called wall functions has been used for calculating wall-bounded turbulent ows. It means that boundar conditions can
Helsini Universit of Technolog CFD-group/ Laborator of Applied Thermodnamics MEMO No CFD/TERMO-9-96 DATE: Ma 24,996 TITLE TEST CALCULATIONS OF VARIOUS TWO-EQUATION LOW REYNOLDS NUMBER TURBULENCE MODELS
More information1 1 Introduction A centrifugal process pump made by Ahlstrom is calculated for industrial evaluation purposes of a FINFLO [1] ow solver. Two turbulenc
Helsinki University of Technology CFD-group/ The Laboratory of Applied Thermodynamics MEMO No CFD/THERMO-23-97 DATE: November 17, 1997 TITLE Flow Simulation of a Centrifugal Process Pump AUTHOR(S) Juha
More informationTensor Visualization. CSC 7443: Scientific Information Visualization
Tensor Visualization Tensor data A tensor is a multivariate quantity Scalar is a tensor of rank zero s = s(x,y,z) Vector is a tensor of rank one v = (v x,v y,v z ) For a symmetric tensor of rank 2, its
More informationHelsinki University of Technology Laboratory of Applied Thermodynamics
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
More informationA DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any
Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»
More informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationChapter 2: Basic Governing Equations
-1 Reynolds Transport Theorem (RTT) - Continuity Equation -3 The Linear Momentum Equation -4 The First Law of Thermodynamics -5 General Equation in Conservative Form -6 General Equation in Non-Conservative
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationSOE3211/2 Fluid Mechanics lecture 2
SOE311/ Fluid Mechanics lecture example example Fluid ows governed by conservation of mass, momentum. We can use this to solve ow problems. Draw box (control volume) around region of interest, then equate
More information1 PART1: Bratu problem
D9: Advanced Numerical Analysis: 3 Computer Assignment 3: Recursive Projection Method(RPM) Joakim Möller PART: Bratu problem. Background The Bratu problem in D is given by xx u + yy u + λe u =, u Γ = ()
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium
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 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 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 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 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 informationViscous flow along a wall
Chapter 8 Viscous flow along a wall 8. The no-slip condition All liquids and gases are viscous and, as a consequence, a fluid near a solid boundary sticks to the boundary. The tendency for a liquid or
More informationLaminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition
Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 18-05 Efficient and robust Schur complement approximations in the augmented Lagrangian preconditioner for high Reynolds number laminar flows X. He and C. Vuik ISSN
More informationIntroduction to Fluid Mechanics
Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the
More informationTHERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW
THERMAL PERFORMANCE EVALUATION OF AN INNOVATIVE DOUBLE GLAZING WINDOW Luigi De Giorgi, Carlo Cima, Emilio Cafaro Dipartimento di Energetica, Politecnico di Torino, Torino, Italy Volfango Bertola School
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 informationNumerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders
Numerical Investigation of Thermal Performance in Cross Flow Around Square Array of Circular Cylinders A. Jugal M. Panchal, B. A M Lakdawala 2 A. M. Tech student, Mechanical Engineering Department, Institute
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationJoule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate
Nonlinear Analysis: Modelling and Control, 27, Vol. 12, No. 3, 37 316 Joule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate M. A. Alim
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 informationCorrected Copy. Flow and Turbulence Modeling and Computation of Shock Buffet Onset for Conventional and Supercritical Airfoils NASA/TP
NASA/TP-1998-2698 Flow and Turbulence Modeling and Computation of Shock Buffet Onset for Conventional and Supercritical Airfoils Robert E. Bartels Langley Research Center, Hampton, Virginia Corrected Copy
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 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 informationModelling of dispersed, multicomponent, multiphase flows in resource industries. Section 3: Examples of analyses conducted for Newtonian fluids
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 3: Examples of analyses conducted for Newtonian fluids Globex Julmester 017 Lecture # 04 July 017 Agenda Lecture
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
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 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 informationM E 320 Professor John M. Cimbala Lecture 10
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Finish our example problem rates of motion and deformation of fluid particles Discuss the Reynolds Transport Theorem (RTT) Show how the RTT
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 informationMathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts 1 and Stephan Matthai 2 3rd Febr
HIGH RESOLUTION POTENTIAL FLOW METHODS IN OIL EXPLORATION Stephen Roberts and Stephan Matthai Mathematics Research Report No. MRR 003{96, Mathematics Research Report No. MRR 003{96, HIGH RESOLUTION POTENTIAL
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationModule 2: Governing Equations and Hypersonic Relations
Module 2: Governing Equations and Hypersonic Relations Lecture -2: Mass Conservation Equation 2.1 The Differential Equation for mass conservation: Let consider an infinitely small elemental control volume
More informationDynamics of Ice Sheets and Glaciers
Dynamics of Ice Sheets and Glaciers Ralf Greve Institute of Low Temperature Science Hokkaido University Lecture Notes Sapporo 2004/2005 Literature Ice dynamics Paterson, W. S. B. 1994. The Physics of
More informationENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES
Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,
More informationEKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009)
EKC314: TRANSPORT PHENOMENA Core Course for B.Eng.(Chemical Engineering) Semester II (2008/2009) Dr. Mohamad Hekarl Uzir-chhekarl@eng.usm.my School of Chemical Engineering Engineering Campus, Universiti
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 informationAnumerical and analytical study of the free convection thermal boundary layer or wall jet at
REVERSED FLOW CALCULATIONS OF HIGH PRANDTL NUMBER THERMAL BOUNDARY LAYER SEPARATION 1 J. T. Ratnanather and P. G. Daniels Department of Mathematics City University London, England, EC1V HB ABSTRACT Anumerical
More informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
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 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 informationFREE SURFACE HYDRODYNAMICS (module 2)
Introduction to FREE SURFACE HYDRODYNAMICS (module 2) Prof Arthur E Mynett Professor of Hydraulic Engineering in addition to the lecture notes by Maskey et al. Introduction to FREE SURFACE HYDRODYNAMICS
More informationAerodynamics. Lecture 1: Introduction - Equations of Motion G. Dimitriadis
Aerodynamics Lecture 1: Introduction - Equations of Motion G. Dimitriadis Definition Aerodynamics is the science that analyses the flow of air around solid bodies The basis of aerodynamics is fluid dynamics
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 informationBoundary layer flow on a long thin cylinder
PHYSICS OF FLUIDS VOLUME 14, NUMBER 2 FEBRUARY 2002 O. R. Tutty a) and W. G. Price School of Engineering Sciences, University of Southampton, Southampton SO17 1BJ, United Kingdom A. T. Parsons QinetiQ,
More informationImplicit numerical scheme based on SMAC method for unsteady incompressible Navier-Stokes equations
172 Pet.Sci.(28)5:172-178 DOI 1.17/s12182-8-27-z Implicit numerical scheme based on SMAC method for unsteady incompressible Navier-Stokes equations Li Zhenlin and Zhang Yongxue School of Mechanical and
More informationTurbulent 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 informationOCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction
OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres
More informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
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 information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationApplication of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations. Xiaoyue Shen
Application of High-Order Summation-by-Parts Operators to the Steady Reynolds-Averaged Navier-Stokes Equations Xiaoyue Shen Supervisor: Prof. David W. Zingg A thesis submitted in conformity with the requirements
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationGG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS
GG303 Lecture 6 8/27/09 1 SCALARS, VECTORS, AND TENSORS I Main Topics A Why deal with tensors? B Order of scalars, vectors, and tensors C Linear transformation of scalars and vectors (and tensors) II Why
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationChapter 2 CONTINUUM MECHANICS PROBLEMS
Chapter 2 CONTINUUM MECHANICS PROBLEMS The concept of treating solids and fluids as though they are continuous media, rather thancomposedofdiscretemolecules, is one that is widely used in most branches
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationContinuum mechanism: Stress and strain
Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the
More informationFINITE-DIFFERENCE IMPLEMENTATION OF LATTICE BOLTZMANN METHOD FOR USE WITH NON-UNIFORM GRIDS
7. ANKARA INTERNATIONAL AEROSPACE CONFERENCE AIAC-2013-143 11-13 September 2013 - METU, Ankara TURKEY FINITE-DIFFERENCE IMPLEMENTATION OF LATTICE BOLTZMANN METHOD FOR USE WITH NON-UNIFORM GRIDS Fatih ÇEVİK
More informationNOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7 Number Pages 3 39 c Institute for Scientific Computing and Information NOVEL FINITE DIFFERENCE SCHEME FOR THE NUMERICAL SOLUTION OF TWO-DIMENSIONAL
More informationOutlines. simple relations of fluid dynamics Boundary layer analysis. Important for basic understanding of convection heat transfer
Forced Convection Outlines To examine the methods of calculating convection heat transfer (particularly, the ways of predicting the value of convection heat transfer coefficient, h) Convection heat transfer
More informationLecture 3: The Navier-Stokes Equations: Topological aspects
Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without
More informationAn Unstructured Grid Incompressible Navier-Stokes Algorithm for Convective Heat. Transfer Based on Artificial Compressibility. Jessica Elaine Kress
An Unstructured Grid Incompressible Navier-Stokes Algorithm for Convective Heat Transfer Based on Artificial Compressibility By Jessica Elaine Kress Approved: Lafayette Taylor Professor of Computational
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 informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationThe conservation equations
Chapter 5 The conservation equations 5.1 Leibniz rule for di erentiation of integrals 5.1.1 Di erentiation under the integral sign According to the fundamental theorem of calculus if f is a smooth function
More informationModications and Clarications for the Implementation of the Spalart-Allmaras Turbulence Model
Modications and Clarications for the Implementation of the Spalart-Allmaras Turbulence Model Steven R. Allmaras, Forrester T. Johnson and Philippe R. Spalart corresponding author: steven.r.allmaras@gmail.com
More informationin order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
G25.2651: Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i = @H @p i _p i =? @H @q i We noted early in the course that an ensemble
More informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More information