MODELING AND SIMULATION OF FREE FLUID OVER A POROUS MEDIUM BY MESHLESS METHOD
|
|
- Earl Price
- 5 years ago
- Views:
Transcription
1 6th uropean Conference on Computational echanics (CC 6) 7th uropean Conference on Computational Fluid Dynamics (CFD 7) 5 June 8, Glasgow, UK ODLING AND SIULATION OF FR FLUID OVR A POROUS DIU BY SLSS TOD AGDALNA IRZWICZAK¹ ¹Poznan University of Technology 5. Słodowsa-Curie Square magdalena.mierzwicza@wp.pl, Key words: the Darcy equation, Brinman equation, effective viscosity, meshless method Astract. The effective viscosity in the Brinman filtration equation is determined y a numerical simulation. The porous medium is modelled y a undle of parallel fiers arranged in a regular square array. In the first step, in order to determine the permeaility, the flow driven y the pressure gradient in an unounded porous medium is considered. The flow is longitudinal with respect to fiers. In the second step, to determine the effective viscosity, the Poisseulle flow in a flat layer of porous medium is considered. Comparison of the tangential stress allows to determine the effective viscosity in filtration equation. In the third step, to determine the length of the area where the Brinman equation is applicale, the flow through a porous medium connected with a pure fluid region is considered The Trefftz method with the special purpose Trefftz functions is applied in numerical simulations. INTRODUCTION In the surrounding of a oundary of a porous region with a high porosity adjacent to the free flow area, two modeling methods are nown. In the first method, the Darcy filtration equation is used to incompressile fluid flow in porous region µ P q () where q is the macroscopic velocity, P is the pressure, µ is the dynamic viscosity of the fluid, is the permeaility of the porous medium. In the free flow area the Navier - Stoes equations are used. The Beavers - Joseph oundary condition is applied for the tangent component of the velocity vector etween this regions []. In such modeling a slip constant in the oundary condition appears. In paper [] this constant was determined in a physical experiment, however, in some papers this constant is determined in the numerical experiment. In the second modeling method, the Brinman filtration equation [] is used for a porous medium with very high porosity in the presence of a free fluid region, or for a wall-ounded porous medium µ q ~ µ P q ()
2 agdalena ierzwicza where ~ µ is the effective viscosity of the porous medium. As a matter of fact, this is the Stoes equation for creeping flow with the Darcy resistant term. In contrast to Darcy s equation, q. () includes the term related to the viscous transfer of momentum. The continuity of the filtration and the free flow velocity is assumed at the contact oundary of the considered areas []. In this modeling case, an effective viscosity in the Brinman equation appears. Almost all authors using the Brinman equation have supposed that the effective viscosity ~ µ is equal to the viscosity μ of a pure fluid. owever, there are some papers in which this viscosity is determined y numerical simulation [3]. The Beavers-Joseph oundary condition appiers in the firs method can e deduced as a consequence of what is now commonly called the Brinman equation. In the prolem of flow in a channel ounded y a thic porous wall one gets the same solution with the Brinman equation as with the Darcy equation together with the Beavers-Joseph condition, provided that one identifies the constant α in this condition as ~ µ µ. In this way, experimental or computational data for the constant α can e used for determination of the effective viscosity. immovale wall high porosity area free flow area immovale wall l l Figure. Porous medium with very high porosity region in the presence of free fluid region etween two impermeale walls. In this paper, a comination of these two modeling approaches is proposed. For this purpose, it is assumed that in the porous area with high porosity there is a layer near the oundary in which the flow is governed y the Brinman filtration equation (), while in the rest of the area the flow is governed y the Darcy filtration equation (). Then we consider three suareas. At the oundaries the continuity of the velocity vector is assumed. The prolem in such a modeling is connected with the width of this layer l and with the unnown effective viscosity µ ~. To demonstrating of the proposed new model the following numerical simulation is carried out. The longitudinal laminar flow (the Poisseulle flow) in a parallel-plate conduit is considered (Fig. ). The first half of the considered region is a porous medium and the second one pure fluid region. The porous medium is modeled as a undle of parallel fires arranged in a square array. The purpose of the present consideration is determination of the effective viscosity µ ~ in the Brinman equation and the width of the layer l in which this equation is to e used. But for do this; the permeaility of the porous medium is required. Then to determine the permeaility the flow with the pressure gradient in unounded porous medium was
3 agdalena ierzwicza considered. Numerical simulations are conducted using the meshless methods: the special purpose Trefftz functions (SPTF) []. Let us consider a plane infinite channel adjoining with a plane infinite porous region in which flow is governed y the Darcy filtration equation dp µ q D for x l (3) dz the Brinman filtration equation dp µ ~ d qb qb µ for l x () dz dx and the Stoes equation dp d w µ for x (5) dz dx For the dimensionless variales w x qb( D) ~ µ W,, Q B( D), K, α (6) dp dp µ µ dz µ dz the governing q. (3-5) otain the form Q D for L (7a) K d QB α Q B for L (7) d K d W for d (7c) which is solved with following oundary conditions Q D Q B for L (8a) Q W dq dw for d d (8) W for (8c) The exact solutions of the prolem (7-8) taes the form Q B ( ) D exp D exp K αk αk (9a) W ( ) C C (9) L where D exp ( K ) ; D exp ( K ) ; B αk B αk ( K ) L L K C exp exp ; C exp exp ; B αk αk αk B αk αk αk L B exp exp. αk αk αk αk 3
4 agdalena ierzwicza DTRINATION OF T PRABILITY One of the method for determination of the permeaility of the porous medium ased on the mesurment filtration velocity at the nown pressure gradient. In the present paper the longitudinal flow prolem is consider y using the special purpose Trefftz functions. The method is semi-analytical. Application of the method gives analytical form of the dimensionless permeaility of the porous medium. Only the unnown coefficients of the solution are otained numerically. a) a W / Figure : Unounded porous medium (a), repeated element Y ) n Ω F (). Let us consider steady, fully developed, laminar, incompressile viscous fluid flow driven y a constant pressure. A porous medium is modelled y a regular array of parallel fiers. The flow is longitudinal with respect to fiers. which are arranged in a regular square (Fig. a) arrays. The equation of motion is given y dp w, in Ω F () µ dz where w is the axial velocity, p is the pressure, μ is the viscosity of the fluid, and Ω F is the fluid domain. For the dimensionless variales w r a W, R,. () dp µ dz the governing q. () otain the form W. () quation () is solvd with the oundary conditions (no slip and symmetry conditon) presented on Fig.. The exact solution of q. () can e expressed using the special purpose Trefftz functions N 8 (, ) ( ) R W R θ R B cos( ) ln, R θ B N R (3) The unnown coefficients B (,, N) are determined y solving the system of linear equations resulting from satisfying of the oundary condition for, using the oundary collocation technique [5].
5 agdalena ierzwicza 5 Using the Darcy law () the longitudinal component of the filtration velocity can e related to the average velocity through the repeated element of the fier system dz dp F q z ϕ µ, () where Ω Ω sec tan, θ θ β β ϕ d RdR W ddy Y W F T F (5) tan β is parameter of the porous medium. The dimensionless component of the permeaility tensor in the direction parallel to the fiers is a function of the numer of collocation points N and can e calculated from 3 cos ln tan 3 tan cos 8 B G B F N N where cos sin,. cos sin G The value of the permeaility parameter is calculated as F β σ. (6) Figure 3. The non-dimensional component of the permeaility tensor in the direction parallel to the fiers and the permeaility parameter. The numerical results (mars) compared with the analytical solution (line)[9] The value of the dimensionless permeaility as a function of the fiers volume fraction is presented in Fig. 3. The prolem was solved for N 7 collocation points on the oundary permeaility; F(φ) φ permeaility parameter; σ(φ) φ
6 agdalena ierzwicza The numerical results (mars) are compared with the analytical presented y Drummond and Tahir [6] (line). Compatiility etween the results is great. 3 DTRINATION OF T FFCTIV VISCOSITY In the rotational viscometer the fluid viscosity is determined y measurement of the torque of the rotating cylindrical surface. The effective viscosity in the Brinman equation can e determined in a similar way. In this case the imaginary physical experiment is realized y a numerical simulation. The porous region is placed etween two parallel plates. One plate is fixed ( q z ) whereas the other moves with a constant velocity ( q z u ). The determination of the effective viscosity is possile y calculation of the shear stress on the moveale plate. The Brinman equation for one-dimensional shear flow in the asence of pressure gradient through the porous layer has the following form ~ d qz µ µ q, z in Ω F. (7) dx where q z is the filtration velocity in direction of fier axis, ~ µ is the effective viscosity of the fluid. The q. (7) can e written in the non-dimensional form d Q Q in d αβf q x ~ z µ h where Q,, α, F,. u µ The exact solution taes the form Ω F. (8) exp exp sinh αβf αβf αβf Q exp exp sinh αβf αβf αβf (9) The tangential stress on the movale plate can e expressed as dq z u τ ~ µ ~ µ coth dx x αβf αβf () On the other hand, the microstructural shear flow prolem etween two plates W W R R R R θ () with no-slip oundary conditions: W on the immovale wall and fiers and W on the movale wall is solved y means of the Trefftz method. The velocity of the flow is approximated y a linear comination of special purpose Trefftz functions. 6
7 agdalena ierzwicza a) y immovale wall x a LR Figure.. Wall-ounded porous medium. movale wall Since the array of fiers is streaed and periodic in one direction it is sufficient to consider the prolem only in one repeated strip. The repeated strip is divided into smaller elements associated with each of the fiers which are called large finite elements. W W W Y R θ W W Y R Figure. 5. The symmetry lane of porous medium Θ W For each large finite element the approximate solution is express as N ( ) R ( ) W ( R, θ ) C ln CR cos[ θ ] R which satisfies exactly the governing q. () and some of the oundary conditions. The unnown coefficients C (,,,N) are determined using the oundary collocation technique y satisfying the remaining oundary conditions (in particular the splitting oundary conditions etween the large finite elements). The tangential stress on the movale wall can e determined from u τ µ dw d Comparing qs. () and (3) the following relationship can e written ~ µ c( ϕ, LR, N ) αβf α µ coth αβf dy () (3) () 7
8 agdalena ierzwicza where ( LR, N ) dw c ϕ, dy and LR denotes the numer of fiers rows. d The value of the constant c and the dimensionless effective viscosity as a function of the fier volume fraction are presented in Fig. 6. The dimensionless effective viscosity was calculated for 5 collocation points on each oundary and 8 rows of fiers. The effective viscosity is smaller than the pure fluid viscosity. The ratio of the effective viscosity to the viscosity of the pure fluid; α at the eginning is decreasing, otains minimum and then is growing. the effective viscosity ; α(φ) φ..6.8 Figure 6. Dimensionless effective viscosity α, and the constant c. T ICROSTRUCTURAL SAR FLOW PROBL Let us consider a layer of the porous medium located etween two fixed plates (Fig. ). The microstructural flow in the layer of porous medium is governed y the dimensionless Poisson equation W (5) with the no slip oundary conditions: W on the immovale wall. the constant c(φ).. φ..6.8 fixed plate a free flow area fixed plate LR Figure 7. The repeating part of the considered channel 8
9 agdalena ierzwicza W W Y W Y R θ W Y R Θ W W W Y Figure 8. Large finite element with oundaty conditions Since the array of fiers is streaed and periodic in one direction it is sufficient to consider the prolem only in one repeated strip (Fig. 7). In all cases the repeated strip is divided into smaller elements associated with each of the fiers (Fig. 8). The flow prolem is solved y means of the Trefftz method. For each large finite element the special purpose Trefftz functions are used to express the approximate solution N (, ) ( ) R W R θ R B cos ln, R θ B N R (6) q. (6) satisfies exactly the governing q. (5) and the oundary conditions W for R and θ for θ {,}. For the free flow area the approximate solution is express as Y W (, Y ) cnfn (, Y ) d ngn (, Y ) n n (7) where F n (,Y) and G n (,Y) are the trial functions F (, Y ) n n n Y ( n )! ( ) G (, Y ) n n n ( n )! ( )! ;!. (8) The unnown coefficients B (,,,N), c n, d (,,,) are determined using the oundary collocation technique y satisfying the oundary conditions (Fig. 8). After solution of the microstructural oundary value prolem the average value of the velocity derivative can e calculated ω ( ) (9) The average velocity derivative ω() for the different value of the fier volume fraction φ are presented in Fig. 9. The prolem was calculated for 5 collocation points on each oundary and 8 rows of fiers. Thus, the width of the area 6. The width of the layer with the Brinman equation can e determined on the asis of the otained results, ω(). It depends on the fier volume fraction. The Darcy filtration equation is valid as long as ω() is close to zero. For greater porosity, this area is wider and approximately equals to half of the width of the porous area. The width of the Brinman area decreases as the the fier volume fraction increases. dw d dy Y 9
10 agdalena ierzwicza ω() φ Figure 9. The average of velocity derivative for different value of the fier volume fraction. ω() CONCLUSIONS In this paper to avoid calculation of the slip constant in the Beavers - Joseph oundary condition the Brinman filtration layer etween porous medium and pure flow fluid area was introduced. In the Brinman filtration equation there are two parameters of the porous medium: the permeaility and the effective viscosity. These two parameters were determined y numerical simulation of an imaginary physical experiment. The porous medium was modeled y a parallel undle of straight fiers arranged in a regular square array. In order to determine the permeaility, the flow driven y the pressure gradient in an unounded porous medium was considered. To determine the effective viscosity, the shear flow in a flat layer of porous medium was considered. To determined the width of the Brinman filtration area the microstructural shear flow was considered. In numerical simulations, the Trefftz method with the special-purpose Trefftz functions was applied. The permeaility of the porous medium decreases as the the fier volume fraction increases.the effective viscosity is lower than the viscosity of the pure fluid. The ratio of the effective viscosity to the viscosity of the pure fluid decreases as the porosity decreases.the width of the Brinman area decreases as the the fier volume fraction increases. ACLOWLGNTS This wor was supported y the grant 5/9/D/ST8/753 from the National Science Center, Poland. RFRNCS [] Beavers G. S. and Joseph D. O. Boundary conditions at a naturally permeale wall, J. Fluid ech. (967) 3:97
11 agdalena ierzwicza [] Brinman. C. A calculation of the viscous force exerted y a flowing fluid on a dense swarm of particles, Appl. Sci. Res. (97) A:7 [3] Kolodziej J.A., ierzwicza. and Grasi J.Computer simulation of the effective viscosity in rinman filtration equation using the Trefftz method, J. ech. ater. Struct. (7) (): 93 [] ierzwicza. and Kolodziej J.A. Comparison of different methods for choosing the collocation points in the oundary collocation method for D-harmonic prolems with special purpose Trefftz functions, ngineering Analysis with Boundary lements () 36: 883 [5] J.A. Kołodziej, A.P. Zielińsi, Boundary Collocation Techniques and their Application in ngineering, WIT Press, Southampton (9). [6] Drummond, J.., Tahir,.I.,. Laminar viscous flow through arrays of parallel solid cylinders. Int. J. ultiphase Flow (98) : 55
Emden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability Porous Layers
4 Int. J. Open Problems Compt. Math., Vol. 0, No., June 07 ISSN 998-66; Copyright ICSRS Publication, 07 www.i-csrs.org Emden-Fowler Equation and Inverse Analysis of Simple Flows through Variable Permeability
More informationA Permeability Function for Brinkman s Equation
Proceedings of the th WSEAS Int. Conf. on MATHEMATICAL METHOS, COMPUTATIONAL TECHNIQUES AN INTELLIGENT SYSTEMS A Permeability Function for Brinman s Equation M.H. HAMAN, M.T. KAMEL, and H.I. SIYYAM 3,
More informationEffect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium
IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. III (Jan. - Fe. 06), 5-65 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility
More informationREE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics
REE 307 - Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics 1. Is the following flows physically possible, that is, satisfy the continuity equation? Substitute the expressions for
More informationOn Temporal Instability of Electrically Forced Axisymmetric Jets with Variable Applied Field and Nonzero Basic State Velocity
Availale at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 19-966 Vol. 6, Issue 1 (June 11) pp. 7 (Previously, Vol. 6, Issue 11, pp. 1767 178) Applications and Applied Mathematics: An International Journal
More informationCHAPTER 4 ANALYTICAL SOLUTIONS OF COUPLE STRESS FLUID FLOWS THROUGH POROUS MEDIUM BETWEEN PARALLEL PLATES WITH SLIP BOUNDARY CONDITIONS
CHAPTER 4 ANALYTICAL SOLUTIONS OF COUPLE STRESS FLUID FLOWS THROUGH POROUS MEDIUM BETWEEN PARALLEL PLATES WITH SLIP BOUNDARY CONDITIONS Introduction: The objective of this chapter is to establish analytical
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 information2, where dp is the constant, R is the radius of
Dynamics of Viscous Flows (Lectures 8 to ) Q. Choose the correct answer (i) The average velocity of a one-dimensional incompressible fully developed viscous flow between two fixed parallel plates is m/s.
More informationUNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume, Numer 4/0, pp. 9 95 UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION
More informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
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 informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationFlow through a Variable Permeability Brinkman Porous Core
Journal of Applied Mathematics and Phsics, 06, 4, 766-778 Published Online April 06 in SciRes http://wwwscirporg/journal/jamp http://dxdoiorg/0436/jamp0644087 Flow through a Variable Permeabilit Brinman
More informationTotal energy in volume
General Heat Transfer Equations (Set #3) ChE 1B Fundamental Energy Postulate time rate of change of internal +kinetic energy = rate of heat transfer + surface work transfer (viscous & other deformations)
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationMobility of Power-law and Carreau Fluids through Fibrous Media
Mobility of Power-law and Carreau Fluids through Fibrous Media Setareh Shahsavari, Gareth H. McKinley Department of Mechanical Engineering, Massachusetts Institute of Technology September 3, 05 Abstract
More informationLaminar Flow. Chapter ZERO PRESSURE GRADIENT
Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More 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 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 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 informationFluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow
OCEN 678-600 Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow Date distributed : 9.18.2005 Date due : 9.29.2005 at 5:00 pm Return your solution either in class or in my mail
More informationHydrodynamic permeability of aggregates of porous particles with an impermeable core
Loughborough University Institutional Repository Hydrodynamic permeability of aggregates of porous particles with an impermeable core This item was submitted to Loughborough University's Institutional
More informationMHD Free convection flow of couple stress fluid in a vertical porous layer
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research,, (6:5- ISSN: 976-86 CODEN (USA: AASRFC MHD Free convection flow of couple stress fluid in a vertical porous layer
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationFLUID MECHANICS PROF. DR. METİN GÜNER COMPILER
FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 5. FLOW IN PIPES 5.1.3. Pressure and Shear Stress
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationEFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW
th International Symposium on Turulence and Shear Flo Phenomena (TSFP, Chicago, USA, July, 7 EFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW Mitsuhiro Nagata Mechanical
More informationP = ρ{ g a } + µ 2 V II. FLUID STATICS
II. FLUID STATICS From a force analysis on a triangular fluid element at rest, the following three concepts are easily developed: For a continuous, hydrostatic, shear free fluid: 1. Pressure is constant
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 information6.1 Steady, One-Dimensional Rectilinear Flows Steady, Spherically Symmetric Radial Flows 42
Contents 6 UNIDIRECTIONAL FLOWS 1 6.1 Steady, One-Dimensional Rectilinear Flows 6. Steady, Axisymmetric Rectilinear Flows 19 6.3 Steady, Axisymmetric Torsional Flows 8 6.4 Steady, Axisymmetric Radial Flows
More informationThe Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)
Chapter 9 Viscous flow The Navier-Stokes equations for an incompressible fluid in indicial notation are i u i = 0 9.1 ρ t + u j j u i = i p + µ 2 j u i 9.2 When these are applied to practical problems,
More informationMHD FORCED FLOW OF A CONDUCTING VISCOUS FLUID THROUGH A POROUS MEDIUM INDUCED BY AN IMPERVIOUS ROTATING DISK
FLUID DYNAMICS MAGNETOHYDRODYNAMICS MHD FORCED FLOW OF A CONDUCTING VISCOUS FLUID THROUGH A POROUS MEDIUM INDUCED BY AN IMPERVIOUS ROTATING DISK BHUPENDRA KUMAR SHARMA, ABHAY KUMAR JHA, R. C. CHAUDHARY
More informationMulti-Point Special Boundary-Value Problems and Applications to Fluid Flow Through Porous Media
Multi-Point Special Boundary-Value Problems and Applications to Fluid Flow Through Porous Media Mohamed A. Hajji Abstract In this paper we propose a numerical scheme based on finite differences for the
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 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 informationTHE EFFECT OF SLIP CONDITION ON UNSTEADY MHD OSCILLATORY FLOW OF A VISCOUS FLUID IN A PLANER CHANNEL
THE EFFECT OF SLIP CONDITION ON UNSTEADY MHD OSCILLATORY FLOW OF A VISCOUS FLUID IN A PLANER CHANNEL A. MEHMOOD, A. ALI Department of Mathematics Quaid-i-Azam University 4530, Islamabad 44000 Pakistan
More informationExercise 5: Exact Solutions to the Navier-Stokes Equations I
Fluid Mechanics, SG4, HT009 September 5, 009 Exercise 5: Exact Solutions to the Navier-Stokes Equations I Example : Plane Couette Flow Consider the flow of a viscous Newtonian fluid between two parallel
More informationParallel Plate Heat Exchanger
Parallel Plate Heat Exchanger Parallel Plate Heat Exchangers are use in a number of thermal processing applications. The characteristics are that the fluids flow in the narrow gap, between two parallel
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 informationNumerical Analysis of MHD Flow of Fluid with One Porous Bounding Wall
Numerical Analysis of MHD Flow of Fluid with One Porous Bounding Wall Ramesh Yadav 1 & Vivek Joseph 2 1Assistant Professor, Department of Mathematics BBDNITM Lucknow U P 2Professor, Department of Mathematics
More informationDesign and Modeling of Fluid Power Systems ME 597/ABE Lecture 7
Systems ME 597/ABE 591 - Lecture 7 Dr. Monika Ivantysynova MAHA Professor Fluid Power Systems MAHA Fluid Power Research Center Purdue University Content of 6th lecture The lubricating gap as a basic design
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 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 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 informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationHydrodynamic Lubrication
ME 383S Bryant February 15, 2005 1 Hydrodynamic Lubrication Fluid Lubricant: liquid or gas (gas bearing) Mechanism: Pressures separate surfaces o Normal loads on bodies o Convergent profile between surfaces
More informationUsing LBM to Investigate the Effects of Solid-Porous Block in Channel
International Journal of Modern Physics and Applications Vol. 1, No. 3, 015, pp. 45-51 http://www.aiscience.org/journal/ijmpa Using LBM to Investigate the Effects of Solid-Porous Bloc in Channel Neda Janzadeh,
More information2 Navier-Stokes Equations
1 Integral analysis 1. Water enters a pipe bend horizontally with a uniform velocity, u 1 = 5 m/s. The pipe is bended at 90 so that the water leaves it vertically downwards. The input diameter d 1 = 0.1
More information10.52 Mechanics of Fluids Spring 2006 Problem Set 3
10.52 Mechanics of Fluids Spring 2006 Problem Set 3 Problem 1 Mass transfer studies involving the transport of a solute from a gas to a liquid often involve the use of a laminar jet of liquid. The situation
More informationHomework #4 Solution. μ 1. μ 2
Homework #4 Solution 4.20 in Middleman We have two viscous liquids that are immiscible (e.g. water and oil), layered between two solid surfaces, where the top boundary is translating: y = B y = kb y =
More informationBuckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads
NASA Technical Paper 3659 Buckling Behavior of Long Symmetrically Laminated Plates Sujected to Shear and Linearly Varying Axial Edge Loads Michael P. Nemeth Langley Research Center Hampton, Virginia National
More informationRheometry. II.1 Introduction
II Rheometry II.1 Introduction Structured materials are generally composed of microstructures dispersed in a homogeneous phase [30]. These materials usually have a yield stress, i.e. a threshold stress
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 informationUnsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe
Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe T S L Radhika**, M B Srinivas, T Raja Rani*, A. Karthik BITS Pilani- Hyderabad campus, Hyderabad, Telangana, India. *MTC, Muscat,
More informationFluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh
Fluid Mechanics Chapter 9 Surface Resistance Dr. Amer Khalil Ababneh Wind tunnel used for testing flow over models. Introduction Resistances exerted by surfaces are a result of viscous stresses which create
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationCESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,
More information2.25 Advanced Fluid Mechanics
MIT Department of Mechanical Engineering.5 Advanced Fluid Mechanics Problem 6.0a This problem is from Advanced Fluid Mechanics Problems by A.H. Shapiro and A.A. Sonin Consider a steady, fully developed
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More informationSupporting Information. Technique for real-time measurements of endothelial permeability in a
Supporting Information Technique for real-time measurements of endothelial permeability in a microfluidic membrane chip using laser-induced fluorescence detection Edmond W.K. Young a,b,, Michael W.L. Watson
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationTurbulence Laboratory
Objective: CE 319F Elementary Mechanics of Fluids Department of Civil, Architectural and Environmental Engineering The University of Texas at Austin Turbulence Laboratory The objective of this laboratory
More informationChapter 6: Incompressible Inviscid Flow
Chapter 6: Incompressible Inviscid Flow 6-1 Introduction 6-2 Nondimensionalization of the NSE 6-3 Creeping Flow 6-4 Inviscid Regions of Flow 6-5 Irrotational Flow Approximation 6-6 Elementary Planar Irrotational
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
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 informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationViscosity and Polymer Melt Flow. Rheology-Processing / Chapter 2 1
Viscosity and Polymer Melt Flow Rheology-Processing / Chapter 2 1 Viscosity: a fluid property resistance to flow (a more technical definition resistance to shearing) Remember that: τ μ du dy shear stress
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationTransport processes. 7. Semester Chemical Engineering Civil Engineering
Transport processes 7. Semester Chemical Engineering Civil Engineering 1. Elementary Fluid Dynamics 2. Fluid Kinematics 3. Finite Control Volume Analysis 4. Differential Analysis of Fluid Flow 5. Viscous
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 informationAnalytical solutions of heat transfer for laminar flow in rectangular channels
archives of thermodynamics Vol. 35(2014), No. 4, 29 42 DOI: 10.2478/aoter-2014-0031 Analytical solutions of heat transfer for laminar flow in rectangular channels WITOLD RYBIŃSKI 1 JAROSŁAW MIKIELEWICZ
More informationINTRODUCTION TO FLUID MECHANICS June 27, 2013
INTRODUCTION TO FLUID MECHANICS June 27, 2013 PROBLEM 3 (1 hour) A perfect liquid of constant density ρ and constant viscosity µ fills the space between two infinite parallel walls separated by a distance
More informationch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows
ch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows ch-01.qxd 8/4/04 2:33 PM Page 3 Introduction 1 Summary The introduction chapter reviews briefly the basic fluid properties
More informationOscillatory MHD Mixed Convection Boundary Layer Flow of Finite Dimension with Induced Pressure Gradient
Journal of Applied Fluid Mechanics, Vol. 9, No., pp. 75-75, 6. Available online at www.jafmonline.net, ISSN 75-57, EISSN 75-65. DOI:.8869/acadpub.jafm.68.5.876 Oscillatory MHD Mixed Convection Boundary
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 informationChapter 6: Momentum Analysis of Flow Systems
Chapter 6: Momentum Analysis of Flow Systems Introduction Fluid flow problems can be analyzed using one of three basic approaches: differential, experimental, and integral (or control volume). In Chap.
More informationThe most common methods to identify velocity of flow are pathlines, streaklines and streamlines.
4 FLUID FLOW 4.1 Introduction Many civil engineering problems in fluid mechanics are concerned with fluids in motion. The distribution of potable water, the collection of domestic sewage and storm water,
More information1. Tasks of designing
1 Lecture #18(14) Designing calculation of cross section of a highly aspect ratio wing Plan: 1 Tass of designing Distribution of shear force between wing spars Computation of the elastic center 4 Distribution
More informationThe dependence of the cross-sectional shape on the hydraulic resistance of microchannels
3-weeks course report, s0973 The dependence of the cross-sectional shape on the hydraulic resistance of microchannels Hatim Azzouz a Supervisor: Niels Asger Mortensen and Henrik Bruus MIC Department of
More informationUnsteady Couette Flow with Transpiration in a Rotating System
Global Journal of Science Frontier Research Mathematics and Decision Sciences Volume 3 Issue 5 Version.0 Year 03 Type : Double Blind Peer Reviewed International Research Journal Publisher: Global Journals
More informationUnsteady Mhd Flow of a Non-Newtonian Fluid Down and Open Inclined Channel with Naturally Permeable Bed
International Journal of Fluids Engineering. ISSN 974-8 Volume 5, Number (), pp. 57-76 International Research Publication House http://www.irphouse.com Unsteady Mhd Flow of a Non-Newtonian Fluid Down and
More informationUNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW
UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW Derivation of uniform flow equation Dimensional analysis Computation of normal depth UNIFORM FLOW 1. Uniform flow is the flow condition obtained from a
More informationBOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW
Proceedings of,, BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW Yunho Jang Department of Mechanical and Industrial Engineering University of Massachusetts Amherst, MA 01002 Email:
More informationInterfacial dynamics
Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic
More informationCandidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2011 2012 FLUID DYNAMICS MTH-3D41 Time allowed: 3 hours Attempt FIVE questions. Candidates must show on each answer book the type
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationPrediction of Coating Thickness
Prediction of Coating Tickness Jon D. Wind Surface Penomena CE 385M 4 May 1 Introduction Tis project involves te modeling of te coating of metal plates wit a viscous liquid by pulling te plate vertically
More informationg (z) = 1 (1 + z/a) = 1
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More information2 Law of conservation of energy
1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion
More informationRing-shaped crack propagation in a cylinder under nonsteady cooling
High Performance Structures and Materials III 5 Ring-shaped crack propagation in a cylinder under nonsteady cooling V. A. Zhornik, Yu. A. Prokopenko, A. A. Rybinskaya & P. A. Savochka Department of Theoretical
More informationNumerical Analysis of Laminar flow of Viscous Fluid Between Two Porous Bounding walls
Numerical Analysis of Laminar flow of Viscous Fluid Between Two Porous Bounding walls Ramesh Yadav Department of Mathematics Babu Banarasi Das National Institute of Technology & Management Lucknow Uttar
More informationVERTICAL SCALES AND SHEAR STRESSES IN WAVE BOUNDARY LAYERS OVER MOVABLE BEDS. Peter Nielsen & Paul A Guard
VERTICAL SCALES AND SHEAR STRESSES IN WAVE BOUNDARY LAYERS OVER MOVABLE BEDS. Peter Nielsen & Paul A Guard ABSTRACT Unified scaling rules are provided for smooth and rough wave boundary layers. It is shown
More informationPHYSICAL MECHANISM OF CONVECTION
Tue 8:54:24 AM Slide Nr. 0 of 33 Slides PHYSICAL MECHANISM OF CONVECTION Heat transfer through a fluid is by convection in the presence of bulk fluid motion and by conduction in the absence of it. Chapter
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationNATURAL CONVECTION IN INCLINED RECTANGULAR POROUS CAVITY SUBJECT TO HEAT FLUXES ON THE LONG SIDE WALLS
Proceedings of 4 th ICCHMT May 7 0, 005, Paris-Cachan, FRANCE ICCHMT 05-53 NATURAL CONVECTION IN INCLINED RECTANGULAR POROUS CAVITY SUBJECT TO HEAT FLUXES ON THE LONG SIDE WALLS L. Storesletten*, D.A.S.
More information