MODELING AND SIMULATION OF FREE FLUID OVER A POROUS MEDIUM BY MESHLESS METHOD

Size: px
Start display at page:

Download "MODELING AND SIMULATION OF FREE FLUID OVER A POROUS MEDIUM BY MESHLESS METHOD"

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

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 information

A Permeability Function for Brinkman s Equation

A 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 information

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium

Effect 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 information

REE Internal Fluid Flow Sheet 2 - Solution Fundamentals of Fluid Mechanics

REE 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 information

On Temporal Instability of Electrically Forced Axisymmetric Jets with Variable Applied Field and Nonzero Basic State Velocity

On 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 information

CHAPTER 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 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 information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (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 information

2, where dp is the constant, R is the radius of

2, 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 information

UNSTEADY POISEUILLE FLOW OF SECOND GRADE FLUID IN A TUBE OF ELLIPTICAL CROSS SECTION

UNSTEADY 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 information

Shell Balances in Fluid Mechanics

Shell 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 information

Number 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 information

6.2 Governing Equations for Natural Convection

6.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 information

CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer

CENG 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 information

Flow through a Variable Permeability Brinkman Porous Core

Flow 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 information

Total energy in volume

Total 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 information

Fluid Mechanics II Viscosity and shear stresses

Fluid 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 information

Mobility of Power-law and Carreau Fluids through Fibrous Media

Mobility 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 information

Laminar Flow. Chapter ZERO PRESSURE GRADIENT

Laminar 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 information

Chapter 9: Differential Analysis

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 information

ENGR Heat Transfer II

ENGR 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 information

Problem 4.3. Problem 4.4

Problem 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 information

Chapter 9: Differential Analysis of Fluid Flow

Chapter 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 information

Fluid Dynamics for Ocean and Environmental Engineering Homework #2 Viscous Flow

Fluid 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 information

Hydrodynamic permeability of aggregates of porous particles with an impermeable core

Hydrodynamic 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 information

MHD Free convection flow of couple stress fluid in a vertical porous layer

MHD 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 information

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition

Fluid 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 information

Studies 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 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 information

FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER

FLUID 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 information

Candidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.

Candidates 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 information

EFFECTS OF STRONG TEMPERATURE GRADIENT ON A COMPRESSIBLE TURBULENT CHANNEL FLOW

EFFECTS 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 information

P = ρ{ g a } + µ 2 V II. FLUID STATICS

P = ρ{ 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 information

Fluid Dynamics Exercises and questions for the course

Fluid 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 information

6.1 Steady, One-Dimensional Rectilinear Flows Steady, Spherically Symmetric Radial Flows 42

6.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 information

The Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)

The 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 information

MHD FORCED FLOW OF A CONDUCTING VISCOUS FLUID THROUGH A POROUS MEDIUM INDUCED BY AN IMPERVIOUS ROTATING DISK

MHD 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 information

Multi-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 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 information

MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow

MOMENTUM 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 information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals 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 information

THE 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 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 information

Exercise 5: Exact Solutions to the Navier-Stokes Equations I

Exercise 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 information

Parallel Plate Heat Exchanger

Parallel 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 information

Joule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate

Joule 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 information

Numerical Analysis of MHD Flow of Fluid with One Porous Bounding Wall

Numerical 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 information

Design and Modeling of Fluid Power Systems ME 597/ABE Lecture 7

Design 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 information

Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions

Simulating 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 information

Turbulent Natural Convection in an Enclosure with Colliding Boundary Layers

Turbulent 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 information

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies

FLUID 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 information

Masters 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. 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 information

Hydrodynamic Lubrication

Hydrodynamic 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 information

Using LBM to Investigate the Effects of Solid-Porous Block in Channel

Using 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 information

2 Navier-Stokes Equations

2 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 information

10.52 Mechanics of Fluids Spring 2006 Problem Set 3

10.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 information

Homework #4 Solution. μ 1. μ 2

Homework #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 information

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads

Buckling 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 information

Rheometry. II.1 Introduction

Rheometry. 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 information

Numerical 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 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 information

Unsteady Flow of a Newtonian Fluid in a Contracting and Expanding Pipe

Unsteady 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 information

Fluid Mechanics. Chapter 9 Surface Resistance. Dr. Amer Khalil Ababneh

Fluid 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 information

UNIVERSITY OF EAST ANGLIA

UNIVERSITY 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 information

This section develops numerically and analytically the geometric optimisation of

This 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 information

CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP

CESSATION 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 information

2.25 Advanced Fluid Mechanics

2.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 information

CONVECTIVE HEAT TRANSFER

CONVECTIVE 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 information

Supporting Information. Technique for real-time measurements of endothelial permeability in a

Supporting 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 information

OE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004

OE4625 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 information

Turbulence Laboratory

Turbulence 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 information

Chapter 6: Incompressible Inviscid Flow

Chapter 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 information

150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces

150A 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 information

Principles of Convection

Principles 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 information

Basic Fluid Mechanics

Basic 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 information

Viscosity and Polymer Melt Flow. Rheology-Processing / Chapter 2 1

Viscosity 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 information

UNIT II CONVECTION HEAT TRANSFER

UNIT 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 information

Transport processes. 7. Semester Chemical Engineering Civil Engineering

Transport 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 information

4.2 Concepts of the Boundary Layer Theory

4.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 information

Analytical solutions of heat transfer for laminar flow in rectangular channels

Analytical 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 information

INTRODUCTION TO FLUID MECHANICS June 27, 2013

INTRODUCTION 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 information

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 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 information

Oscillatory MHD Mixed Convection Boundary Layer Flow of Finite Dimension with Induced Pressure Gradient

Oscillatory 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 information

Differential relations for fluid flow

Differential 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 information

Chapter 6: Momentum Analysis of Flow Systems

Chapter 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 information

The most common methods to identify velocity of flow are pathlines, streaklines and streamlines.

The 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 information

1. Tasks of designing

1. 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 information

The dependence of the cross-sectional shape on the hydraulic resistance of microchannels

The 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 information

Unsteady Couette Flow with Transpiration in a Rotating System

Unsteady 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 information

Unsteady Mhd Flow of a Non-Newtonian Fluid Down and Open Inclined Channel with Naturally Permeable Bed

Unsteady 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 information

UNIFORM FLOW CRITICAL FLOW GRADUALLY VARIED FLOW

UNIFORM 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 information

BOUNDARY LAYER ANALYSIS WITH NAVIER-STOKES EQUATION IN 2D CHANNEL FLOW

BOUNDARY 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 information

Interfacial dynamics

Interfacial dynamics Interfacial dynamics Interfacial dynamics = dynamic processes at fluid interfaces upon their deformation Interfacial rheological properties: elasticity, viscosity, yield stress, Relation between macroscopic

More information

Candidates must show on each answer book the type of calculator used. Only calculators permitted under UEA Regulations may be used.

Candidates 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 information

Number 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 information

Prediction of Coating Thickness

Prediction 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 information

g (z) = 1 (1 + z/a) = 1

g (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 information

2 Law of conservation of energy

2 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 information

Ring-shaped crack propagation in a cylinder under nonsteady cooling

Ring-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 information

Numerical 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 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 information

VERTICAL 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 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 information

PHYSICAL MECHANISM OF CONVECTION

PHYSICAL 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 information

Microfluidics 1 Basics, Laminar flow, shear and flow profiles

Microfluidics 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 information

Introduction to Marine Hydrodynamics

Introduction 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 information

NATURAL CONVECTION IN INCLINED RECTANGULAR POROUS CAVITY SUBJECT TO HEAT FLUXES ON THE LONG SIDE WALLS

NATURAL 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