On two fractional step finite volume and finite element schemes for reactive low Mach number flows
|
|
- Cameron Shepherd
- 5 years ago
- Views:
Transcription
1 Fourth International Symposium on Finite Volumes for Complex Applications - Problems and Perspectives - July 4-8, 2005 / Marrakech, Morocco On two fractional step finite volume and finite element schemes for reactive low Mach number flows Fabrice Babik *, Thierry Gallouët **, Jean-Claude Latché *, Sylvain Suard *, and Didier Vola * * Institut de Radioprotection et de Sûreté Nucléaire, France ** Laboratoire d Analyse, Topologie, Probabilité, UMR 6632, CMI, France 1/21
2 Outline 1. Introduction 2. Finite Volumes formulation 3. Two numerical schemes 4. Numerical tests Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 2/21
3 Introduction 1. Introduction 1.1 Presentation 1.2 Toward a numerical scheme 2. Finite Volumes formulation 3. Two numerical schemes 4. Numerical tests Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 3/21
4 1.1 Presentation Governing equations: ϱ t + ϱu = 0 ϱu + (ϱu u) = τ p + (ρ ρ )g t ϱz + ϱzu = (κ z) t ϱ = G(z) + boundary and initial conditions Several physical problems enter this abstract framework: Natural convection flows in the low Mach number limit: z stands for the temperature, G(z) is the equation of state of perfect gaz (at fixed pressure in case of open boundaries; otherwise, the preceding system must be supplemented by an evolution equation for the thermodynamic pressure). Non-premixed laminar flame: z is a Zeldovitch s variable called mixture-fraction, G(z) derived from T (z), Y F (z), Y Ox (z) and Y P (z). Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 4/21
5 1.2 Toward a numerical scheme (1/2) Particular features and requirements: the unknown field z must have convergent and stability properties construction of a FV numerical scheme which reproduces these features at the discrete level ϱ = G(z): no independant evolution equation can be stated for the fluid density the mass balance equation acts as a constraint on the velocity - use of fractional step schemes, namely projection methods FEs for an accurate velocity and pressure computation - nonconforming velocity approximations (dof located at the center of the faces) for a coupling with a FV method for the advection-diffusion of z Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 5/21
6 1.2 Toward a numerical scheme (2/2) Overshoots or undershoots (even small) have a dramatic effect in chemically reacting flows From both a physical and numerical viewpoint, all the species concentrations remain positive A maximum principle stands for the unknown field z As the velocity is not divergence free, the conservation equation for z does not imply any maximum principle... But it can be combined with the mass balance equation as follows: ϱz t + ϱzu = ϱ( z t + u z) + z ( ϱ t + ϱu) }{{} =0 = (κ z) For this latter equation, if z(x, 0) [0, 1] and z takes values in [0, 1] on Dirichlet boundaries, then z(x, t) [0, 1], x, t. Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 6/21
7 Mixed FV / FE formulation 1. Introduction 2. Finite Volumes formulation 2.1 A finite volumes scheme for the transport of z 2.2 High order interpolation FV schemes 3. Two numerical schemes 4. Numerical tests Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 7/21
8 2.1 A finite volumes scheme for the transport of z (1/4) Definition: M is a M-matrix if M ii > 0, 1 i N, M ij 0, for i j M ii M ij 0, i j i M ii M ij > 0, for one node j i If M is a M-matrix: if MX 0, = X 0 If M t is an M-matrix, MX 0 = X 0 Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 8/21
9 2.1 A finite volumes scheme for the transport of z (2/4) Scheme: solve in sequence a prediction for the density, then the (now linear) convection-diffusion equation for z, and, finally, update the density using the constitutive equation. m K ϱ K ϱ n K k ϱ K z n+1 K m K k ϱ n+1 K = G(zn+1 K ) + σ E K m σ v σ ϱ σ = 0 ϱn K zn K + m σ v σ ϱ σ zσ n+1 σ E K + σ E K m σ κ (z n+1 K d σ zn+1 L ) = 0 where : v σ stands for the normal velocity to the edge σ ϱ σ and z n+1 σ stand for the standard upwind approximation for ϱ and z n+1. Of course, the prediction step makes sense only if the mass balance holds. Exemple: see the full scheme hereafter. Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 9/21
10 2.1 A finite volumes scheme for the transport of z (3/4) The discrete operator associated to the first equation is the transposed of an M-matrix. The discrete operator associated to the second equation is an M-matrix: Σ= L A K,L = m K ϱ K k + m σ v σ,k ϱ σ = m ϱn K K k > 0 σ E K Let z n and z n be defined by z n = max K zn K, zn = min K zn K. Then: Idea of the proof: K, z n+1 K [zn, z n ] ϱ K (z n+1 K m zn ) K k + σ E K m σ + m σ v σ ϱ σ (zσ n+1 z n ) σ E K κ [(z n+1 K d zn ) (z n+1 L σ z n )] = m K ϱ n K zn K k Σ z n = m K ϱ n K (zn K zn ) k }{{} 0 Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 10/21
11 2.1 A finite volumes scheme for the transport of z (4/4) A L 2 stability property. Multiplying by z the advection-diffusion equation for z and integrating over Ω yields: Ω [ ] ϱz t + (ϱzu) = 1 2 Ω = 1 2 Ω = 1 ϱ z2 2 t Ω z = Ω ϱ z2 t + ϱu (z2 ) ϱ z2 t z2 (ϱu) + z2 ϱ t = 1 2 t [ ϱ t + (ϱu)] z2 + ϱz z Ω t Ω ϱz 2 = Ω + ϱzu z (mass balance equation) (integration by parts) (mass balance equation) κ z 2 The same computation can be followed "line by line" in the discrete case, which provides: a stability result, a tool for proving an error estimate? Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 11/21
12 2.2 High order interpolation FV schemes A second order FV scheme for diffusive cases: Numerical experiments show that the spatial accuracy is limited to first order as soon as an upwind discretization is used for the prediction of the density. Use of a centered deferred correction scheme (2 nd order): ρ ρ n + U (ρv) n+1 + β ( C (ρv) n + U (ρv) n ) = 0 The β parameter is chosen to ensure that ρ remains positive. Use of an hybrid FV interpolation scheme for the transport of z ϱz t + H (ϱzu) = (κ z) with H (φ) = α C (φ)+(1 α) U (φ), α [0, 1] Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 12/21
13 Two numerical schemes 1. Introduction 2. Finite Volumes formulation 3. Two numerical schemes 3.1 A non-conservative scheme 3.2 A conservative scheme 4. Numerical tests Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 13/21
14 3.1 A non-conservative scheme solve for ϱ ϱ ϱ n + U ϱu n = 0 solve for z n+1 ϱz n+1 ϱ n z n + U ϱz n+1 u n C,κ z n+1 = 0 ϱ n+1 given by ϱ n+1 = G(z n+1 ) solve for ũ solve for p n+1, u n+1 ϱũ ϱ n u n + ( ϱu n ũ) τ(ũ) + p n = g n+1 ϱ n+1 u n+1 ϱũ + (p n+1 p n ) = 0 ϱ n+1 u n+1 = ϱn+1 ϱ n the density is updated after the computation of z n+1 this scheme does not conserve the quantity ϱz Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 14/21
15 3.2 A conservative scheme solve for ϱ solve for z ϱ n+1 given by solve for ũ solve for p n+1, u n+1 solve for z n+1 ϱ ϱ n + U ϱu n = 0 ϱ z ϱ n z n + U ϱ zu n C,κ z = 0 ϱ n+1 = G( z) ϱũ ϱ n u n + ( ϱu n ũ) τ(ũ) + p n = g n+1 ϱ n+1 u n+1 ϱũ + (p n+1 p n ) = 0 ϱ n+1 u n+1 = ϱn+1 ϱ n ϱ n+1 z n+1 ϱ n z n + CU ϱ n+1 z n+1 u n+1 C,κ z n+1 = 0 stability of the last step requires: - a centered approximation for the density - an upwind approximation for the variable z Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 15/21
16 Numerical tests 1. Introduction 2. Finite Volumes formulation 3. Two numerical schemes 4. Numerical tests 4.1 Analytical premixed flame propagation 4.2 Planar non-premixed flame propagation 4.3 Natural convection flows Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 16/21
17 4.1 Analytical premixed flame propagation (1/2) Governing equations: ϱ t + ϱu = 0 ϱu + (ϱu u) = τ p t ϱz + ϱzu = (κ z) + f(z) t ϱ = 1/(1 γ + γz) Analytical solution: z(x, t) = 1 2 ( [ ]) x ut 1 tanh δ propagation burnt mixture fresh mixture reaction zone Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 17/21
18 4.1 Analytical premixed flame propagation (2/2) z, N-C S z, C S rho, N-C S, order 2 z, N-C S, order 2 slope order 2 L 2 error norm 0.1 slope order 1 N = 125, N-C S C S N = 500, N-C S C S N = 2000, N-C S C S Time-step L 2 error norm Dx L 2 error norm obtained with the non-conservative (N-C S) and conservative scheme (C S), as a function of: left: time-step for differents grid size; right: grid size for a time-step of Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 18/21
19 4.2 Planar non-premixed flame propagation Infinitely fast chemistry assumption: Burke-Schumann approximation The variable z stands for the mixture fraction A very stiff problem: Large gradients of density and temperature appear near z S 0.05 Mixtures fraction and density mixture fraction, dt = s dt = s dt = s density, dt = s dt = s dt = s Mixtures fraction and density mixture fraction, N = 200 N = 400 N = 800 density, N = 200 N = 400 N = Reduced abscissa Reduced abscissa Distributions of mixture fraction and density along the flow direction obtained with the non-conservative scheme for three time steps (left) and three mesh sizes (right). Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 19/21
20 4.3 Natural convection flows (1/2) A buoyancy-driven flow in a two-dimensional square enclosure [0,L] [0,L] with large temperature differences Two insulated walls and two vertical walls heated to T h respectively cooled down to T c. zero heat flux Flow specification: Dimension 2D T h hot fluid g T c Geometry Temporal Cartesian Steady laminar & Thermal cold fluid Gas law Perfect gas zero heat flux Nusselt number: Nu(y) = L k 0 (T h T c ) k T x Average Nusselt number: y=l Nu = 1 L y=0 Nu h + Nu c = 0 Nu(y) dy, w Governing equations: ϱ t + ϱu = 0 ϱu + (ϱu u) = τ p + (ρ ρ )g t ϱh + ϱhu = (κ h) + dp t dt P = ϱrt Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 20/21
21 4.3 Natural convection flows (2/2) Test case T2 (H. Paillère and P. Le Quéré): Ra = P r g ρ2 0 (T h T c ) L 3 T 0 µ 2 0 = 10 6, Sutherland s Law, ɛ = T h T c T h + T c = 0.6 Initial conditions: The problem is completely defined by the Rayleigh number, the value of ɛ and the following coefficient: P 0 = Pa, T 0 = 600 K, ρ 0 = P 0 /(RT 0 ), P r = 0.71, γ = 1.4, g = 9.81 m/s 2 Boundary conditions: hot wall: T h = T 0 (1 + ɛ), cold wall: T c = T 0 (1 ɛ) horizontal walls: adiabatic conditions, all walls: no-slip condition for the velocity Ref relative error Nu h e := 1 + Nu c / Nu h O(h) for scheme 1 Fourth International Symposium on Finite Volumes for Complex Applications, July 4-8, 2005 / Marrakech 21/21
Convergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationTowards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below
Copyright 2008 Tech Science Press FDMP, vol.4, no.4, pp.263-269, 2008 Towards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below G. Accary 1, S. Meradji
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 informationEffect of an adiabatic fin on natural convection heat transfer in a triangular enclosure
American Journal of Applied Mathematics 2013; 1(4): 78-83 Published online November 10, 2013 (http://www.sciencepublishinggroup.com/j/ajam) doi: 10.11648/j.ajam.20130104.16 Effect of an adiabatic fin on
More informationLecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory
Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where
More information10. Buoyancy-driven flow
10. Buoyancy-driven flow For such flows to occur, need: Gravity field Variation of density (note: not the same as variable density!) Simplest case: Viscous flow, incompressible fluid, density-variation
More informationarxiv: v1 [math.na] 16 Oct 2007
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher arxiv:070.2987v [math.na] 6 Oct 2007 AN UNCONDITIONNALLY STABLE PRESSURE CORRECTION
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 5 NATURAL CONVECTION HEAT TRANSFER BASIC CONCEPTS MECHANISM OF NATURAL
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
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 informationNatural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationVector and scalar variables laminar natural convection in 2D geometry arbitrary angle of inclination
Vector and scalar variables laminar natural convection in 2D geometry arbitrary angle of inclination Miomir Raos and Ljubiša Nešić Abstract The aim of this study 1 is to analyze vector and scalar depending
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationBenchmark solutions for the natural convective heat transfer problem in a square cavity
Benchmark solutions for the natural convective heat transfer problem in a square cavity J. Vierendeels', B.Merci' &L E. Dick' 'Department of Flow, Heat and Combustion Mechanics, Ghent University, Belgium.
More informationAbstract. Introduction
Combined forced and natural convection in a square cavity - numerical solution and scale analysis A.T. Franco/ M.M. Ganzarolli'' "DAMEC, CEFET, PR 80230-901, Curitiba, PR Brasil >>DE, FEM, UNICAMP 13081-970,
More informationNUMERICAL SIMULATION OF THERMAL CONVECTION IN A CLOSED CAVITY IN THE PRESENCE OF A THIN HORIZONTAL HEATED PLATE
Proceedings of CHT-1 ICHMT International Symposium on Advances in Computational Heat Transfer July 1-6, 1, Bath, England CHT-1-NC17 NUMERICAL SIMULATION OF THERMAL CONVECTION IN A CLOSED CAVITY IN THE
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
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 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 informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
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 informationDimensionless Numbers
1 06.10.2017, 09:49 Dimensionless Numbers A. Salih Dept. of Aerospace Engineering IIST, Thiruvananthapuram The nondimensionalization of the governing equations of fluid flow is important for both theoretical
More informationFREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER
Free Convective Heat Transfer From an Object at Low Rayleigh Number FREE CONVECTIVE HEAT TRANSFER FROM AN OBJECT AT LOW RAYLEIGH NUMBER Md. Golam Kader and Khandkar Aftab Hossain * Department of Mechanical
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 informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More information5. Coupling of Chemical Kinetics & Thermodynamics
5. Coupling of Chemical Kinetics & Thermodynamics Objectives of this section: Thermodynamics: Initial and final states are considered: - Adiabatic flame temperature - Equilibrium composition of products
More informationConvection Workshop. Academic Resource Center
Convection Workshop Academic Resource Center Presentation Outline Understanding the concepts Correlations External Convection (Chapter 7) Internal Convection (Chapter 8) Free Convection (Chapter 9) Solving
More informationParabolic Flow in Parallel Plate Channel ME 412 Project 4
Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationAbstract. Nomenclature. A^ = dimensionless half amplitude ' = half amplitude of heat flux variation * - thermal conductivity
Unsteady natural connective flow in an enclosure with a periodically varying side wall heat flux PH. Oosthuizen, J.T. Paul Heat Transfer Laboratory, Department ofmechanical Engineering, Queen's University,
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationModelling Turbulent Heat Transfer in a Natural Convection Flow
Journal of Applied Mathematics and Physics, 2014, 2, 662-670 Published Online June 2014 in SciRes. http://www.scirp.org/journal/jamp http://dx.doi.org/10.4236/jamp.2014.27073 Modelling Turbulent Heat Transfer
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationFundamental Concepts of Convection : Flow and Thermal Considerations. Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.
Fundamental Concepts of Convection : Flow and Thermal Considerations Chapter Six and Appendix D Sections 6.1 through 6.8 and D.1 through D.3 6.1 Boundary Layers: Physical Features Velocity Boundary Layer
More informationOn Clean Cooling Systems for Wind Turbine Nacelle operating in Hot Climate
International Conférence en Clean Cooling Technologiesin the ME NA Regions ICT3_MENA 201 Bou Smail, W. Tipaza, 5-6 October 2015 On Clean Cooling Systems for Wind Turbine Nacelle operating in Hot Climate
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 informationENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD
HEFAT 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 6 8 July Malta ENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD Nawaf
More informationBounds on Surface Stress Driven Flows
Bounds on Surface Stress Driven Flows George Hagstrom Charlie Doering January 1, 9 1 Introduction In the past fifteen years the background method of Constantin, Doering, and Hopf has been used to derive
More informationNATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION
NATURAL CONVECTION WITHIN TRAPEZOIDAL ENCLOSURE WITH TWO BAFFLES: EFFECT OF VARIOUS ANGLES OF INCLINATION Éliton Fontana, eliton_fontana@hotmail.com Universidade Federal de Santa Catarina - UFSC Adriano
More informationLAMINAR NATURAL CONVECTION IN VERTICAL 2D GLAZING CAVITIES
Mechanical and Industrial Engineering University of Massachusetts, Amherst AMINAR NATURA CONVECTION IN VERTICA 2D GAZING CAVITIES Bhaskar Adusumalli ABSTRACT Finite element predictions of natural convection
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationInfluence of Heat Transfer Process in Porous Media with Air Cavity- A CFD Analysis
Proceedings of the 4 th International Conference of Fluid Flow, Heat and Mass Transfer (FFHMT'17) Toronto, Canada August 21 23, 2017 Paper No. 161 DOI: 10.11159/ffhmt17.161 Influence of Heat Transfer Process
More informationNATURAL CONVECTION OF AIR IN TILTED SQUARE CAVITIES WITH DIFFERENTIALLY HEATED OPPOSITE WALLS
Proceedings of the International onference on Mechanical Engineering 0 (IME0 8-0 December 0, Dhaka, Bangladesh IME- NATURAL ONVETION OF AIR IN TILTED SQUARE AVITIES WIT DIFFERENTIALLY EATED OPPOSITE WALLS
More informationPHYSICAL MECHANISM OF NATURAL CONVECTION
1 NATURAL CONVECTION In this chapter, we consider natural convection, where any fluid motion occurs by natural means such as buoyancy. The fluid motion in forced convection is quite noticeable, since a
More informationINSTRUCTOR: PM DR MAZLAN ABDUL WAHID
SMJ 4463: HEAT TRANSFER INSTRUCTOR: PM ABDUL WAHID http://www.fkm.utm.my/~mazlan TEXT: Introduction to Heat Transfer by Incropera, DeWitt, Bergman, Lavine 5 th Edition, John Wiley and Sons Chapter 9 Natural
More informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationFlow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers ( )
Advances in Fluid Mechanics VII 391 Flow patterns and heat transfer in square cavities with perfectly conducting horizontal walls: the case of high Rayleigh numbers (10 6 10 9 ) R. L. Frederick & S. Courtin
More informationCombined Natural Convection and Thermal Radiation in an Inclined Cubical Cavity with a Rectangular Pins Attached to Its Active Wall
Periodicals of Engineering and Natural Sciences ISSN 2303-4521 Vol.5, No.3, November 2017, pp. 347~354 Available online at:http://pen.ius.edu.ba Combined Natural Convection and Thermal Radiation in an
More informationHEFAT th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics September 2005, Cairo, Egypt AA10
HEFAT5 4 th International Conference on Heat Transfer, Fluid Mechanics, and Thermodynamics 9- September 5, Cairo, Egypt AA Numerical Study of Natural Convection Heat Transfer in Enclosures with Conducting
More informationNUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS
NUMERICAL STUDY OF MELTING OF TIN WITHIN A RECTANGULAR CAVITY INCLUDING CONVECTIVE EFFECTS Christiano Garcia da Silva Santim, chrisoff22@yahoo.com.br Luiz Fernando Milanez, milanez@fem.unicamp.br Universidade
More informationVALIDATION PROCESS OF THE ISIS CFD SOFTWARE FOR FIRE SIMULATION. C. Lapuerta, F. Babik, S. Suard, L. Rigollet
VALIDATION PROCESS OF THE ISIS CFD SOFTWARE FOR FIRE SIMULATION C. Lapuerta, F. Babik, S. Suard, L. Rigollet Institut de Radioprotection et de Sûreté Nucléaire (IRSN) DPAM, SEMIC, LIMSI BP3 13115 Saint
More informationInternational Journal of Thermal Sciences
International Journal of Thermal Sciences 105 (2016) 137e158 Contents lists available at ScienceDirect International Journal of Thermal Sciences journal homepage: www. elsevier. com/ locate/ ijts Natural
More informationMixing Efficiency of Inertia-Gravity Wave Breaking: an oceanographic perspective
Mixing Efficiency of Inertia-Gravity Wave Breaking: an oceanographic perspective Pascale Lelong, NWRA, Seattle USA Pascale Bouruet-Aubertot, U. of Paris VI, France Sylvain Beaufils, ENS-Cachan, France
More informationCombustion Behind Shock Waves
Paper 3F-29 Fall 23 Western States Section/Combustion Institute 1 Abstract Combustion Behind Shock Waves Sandeep Singh, Daniel Lieberman, and Joseph E. Shepherd 1 Graduate Aeronautical Laboratories, California
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationCoupling of Lattice Boltzmann Equation and Finite Volume Method to Simulate Heat Transfer in a Square Cavity
Copyright 2009 Tech Science Press FDMP, vol.5, no.3, pp.283-295, 2009 Coupling of Lattice Boltzmann Equation and Finite Volume Method to Simulate Heat Transfer in a Square Cavity Ahmed Mezrhab 1 and Hassan
More informationIntroduction to Heat and Mass Transfer. Week 9
Introduction to Heat and Mass Transfer Week 9 補充! Multidimensional Effects Transient problems with heat transfer in two or three dimensions can be considered using the solutions obtained for one dimensional
More informationTransport processes. 7. Semester Chemical Engineering Civil Engineering
Transport processes 7. Semester Chemical Engineering Civil Engineering 1 Course plan 1. Elementary Fluid Dynamics 2. Fluid Kinematics 3. Finite Control Volume nalysis 4. Differential nalysis of Fluid Flow
More informationComputational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor
Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor SD Ravi 1, NKS Rajan 2 and PS Kulkarni 3 1 Dept. of Aerospace Engg., IISc, Bangalore, India. ravi@cgpl.iisc.ernet.in
More informationTURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows
TURBINE BURNERS: Engine Performance Improvements; Mixing, Ignition, and Flame-Holding in High Acceleration Flows Presented by William A. Sirignano Mechanical and Aerospace Engineering University of California
More informationNumerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 541 552 Numerical Investigation of Combined Buoyancy and Surface Tension Driven Convection in an Axi-Symmetric Cylindrical Annulus M. Sankar
More informationGeneralizing the Boussinesq Approximation to Stratified Compressible Flow
Generalizing the Boussinesq Approximation to Stratified Compressible Flow Dale R. Durran a Akio Arakawa b a University of Washington, Seattle, USA b University of California, Los Angeles, USA Abstract
More informationINDIAN INSTITUTE OF TECHNOOGY, KHARAGPUR Date: -- AN No. of Students: 5 Sub. No.: ME64/ME64 Time: Hours Full Marks: 6 Mid Autumn Semester Examination Sub. Name: Convective Heat and Mass Transfer Instructions:
More informationNumerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall
Heat Transfer Research, 2011, Vol. 42, No. 3 Numerical Study of Free Convection Heat Transfer in a Square Cavity with a Fin Attached to Its Cold Wall SAEID JANI, 1* MEYSAM AMINI, 2 and MOSTAFA MAHMOODI
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationTHE INFLUENCE OF HEIGHT RATIO ON RAYLEIGH-NUMBER SCALING AND STABILITY OF HORIZONTAL CONVECTION
Seventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 9-11 December 2009 THE INFLUENCE OF HEIGHT RATIO ON RAYLEIGH-NUMBER SCALING AND STABILITY OF HORIZONTAL
More informationNumerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles
Numerical Investigation of Fluid and Thermal Flow in a Differentially Heated Side Enclosure walls at Various Inclination Angles O. A. SHAHRUL Faculty of Mechanical & Manufacturing Engineering Universiti
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationIntroduction to Heat and Mass Transfer. Week 14
Introduction to Heat and Mass Transfer Week 14 HW # 7 prob. 2 Hot water at 50C flows through a steel pipe (thermal conductivity 14 W/m-K) of 100 mm outside diameter and 8 mm wall thickness. During winter,
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
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 informationDIRECT NUMERICAL SIMULATION OF NATURAL CONVECTION OVER HORIZONTAL PLATES
XXIII Congresso Nazionale UIT sulla Trasmissione del Calore Parma, - Giugno DIRECT NUMERICAL SIMULATION OF NATURAL CONVECTION OVER HORIZONTAL PLATES Daniele Venturi, Beatrice Pulvirenti, Sandro Salvigni
More informationME 331 Homework Assignment #6
ME 33 Homework Assignment #6 Problem Statement: ater at 30 o C flows through a long.85 cm diameter tube at a mass flow rate of 0.020 kg/s. Find: The mean velocity (u m ), maximum velocity (u MAX ), and
More informationA Stochastic Projection Method for Fluid Flow
Journal of Computational Physics 8, 9 44 (22) doi:.6/jcph.22.74 A Stochastic Projection Method for Fluid Flow II. Random Process Olivier P. Le Maître, Matthew T. Reagan, Habib N. Najm, Roger G. Ghanem,
More informationTurbulent Natural Convection in an Enclosure with Colliding Boundary Layers
Turbulent Natural Convection in an Enclosure with Colliding Boundary Layers Abstract Mutuguta John Wanau 1* 1. School of Pure and Applied Sciences, Murang a University of Technology, P.O box 75-10200,
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 informationHeat transfer increase with thin fins in three dimensional enclosures
157 Heat transfer increase with thin fins in three dimensional enclosures R. L. Frederick & S. Samper Universidad de Chile, Departamento de Ingeniería Mecánica, Santiago, Chile Abstract Heat transfer enhancement
More informationSolving PDEs with freefem++
Solving PDEs with freefem++ Tutorials at Basque Center BCA Olivier Pironneau 1 with Frederic Hecht, LJLL-University of Paris VI 1 March 13, 2011 Do not forget That everything about freefem++ is at www.freefem.org
More informationHeat processes. Heat exchange
Heat processes Heat exchange Heat energy transported across a surface from higher temperature side to lower temperature side; it is a macroscopic measure of transported energies of molecular motions Temperature
More informationTwo Phase Thermal Boundary Layer Flow
Two Phase Thermal Boundary Layer Flow P. K. Tripathy Deptt. Of Mathematics & Science, U. C. P. Engg. School, Berhampur (Orissa), INDIA S. K. Mishra Deptt. Of Mathematics, R. C. M. Science College, Khallokote
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 informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
More informationHeat Transfer Modeling using ANSYS FLUENT
Lecture 1 - Introduction 14.5 Release Heat Transfer Modeling using ANSYS FLUENT 2013 ANSYS, Inc. March 28, 2013 1 Release 14.5 Outline Modes of Heat Transfer Basic Heat Transfer Phenomena Conduction Convection
More informationHybrid Discontinuous Galerkin methods for incompressible flow problems
Hybrid Discontinuous Galerkin methods incompressible flow problems Christoph Lehrenfeld, Joachim Schöberl MathCCES Computational Mathematics in Engineering Workshop Linz, May 31 - June 1, 2010 Contents
More informationFlame / wall interaction and maximum wall heat fluxes in diffusion burners
Flame / wall interaction and maximum wall heat fluxes in diffusion burners de Lataillade A. 1, Dabireau F. 1, Cuenot B. 1 and Poinsot T. 1 2 June 5, 2002 1 CERFACS 42 Avenue Coriolis 31057 TOULOUSE CEDEX
More informationTOPICAL PROBLEMS OF FLUID MECHANICS 97
TOPICAL PROBLEMS OF FLUID MECHANICS 97 DOI: http://dx.doi.org/10.14311/tpfm.2016.014 DESIGN OF COMBUSTION CHAMBER FOR FLAME FRONT VISUALISATION AND FIRST NUMERICAL SIMULATION J. Kouba, J. Novotný, J. Nožička
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationNatural Convection Inside a Rectangular Enclosure with Two Discrete Heat Sources Placed on Its Bottom Wall
Proceedings of the 2 nd World Congress on Momentum, Heat and Mass Transfer (MHMT 17) Barcelona, Spain April 6 8, 2017 Paper No. ENFHT 108 ISSN: 2371-5316 DOI: 10.11159/enfht17.108 Natural Convection Inside
More informationBy ablation we mean the recession of a surface due to heating, usually by a hot gas. It is the key process for
16.50 Lecture 15 Subject: Ablative cooling By ablation we mean the recession of a surface due to heating, usually by a hot gas. It is the key process for a) Re-entry heat shields b) Solid propellant nozzles
More informationOn Computing Reactive Flows Using Slow Manifolds
On Computing Reactive Flows Using Slow Manifolds by S. Paolucci 1, J. M. Powers 2, and S. Singh 3 Department of Aerospace and Mechanical Engineering University of Notre Dame Notre Dame, Indiana 46556-5637,
More informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More information