The dependence of the cross-sectional shape on the hydraulic resistance of microchannels
|
|
- Philippa Arnold
- 6 years ago
- Views:
Transcription
1 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 Micro and Nanotechnology Technical University of Denmark 5 July 004
2 Contents 1 Introduction 3 Theory 4.1 Poiseuille Flow in a circular pipe Velocity field Flow rate Poiseuille Flow in an elliptic pipe Velocity field Flow rate Poiseuille Flow in a rectangular pipe - Fourier Transformation Velocity field Flow rate Hydraulic resistance Elliptic cross-section Rectangular cross-section FEMLAB simulations 13 4 Conclusion 16
3 Chapter 1 Introduction For a steady incompressible flow through a channel the hydraulic resistance of the channel depends on the cross-section shape of the channel. The objective of this project is to investigate the dependence between the hydraulic resistance and the cross-section of the channel in detail as well as to investigate the hydraulic resistance for various cross-section geometries. In the following, Chapter, a description of the Poiseuille flow through a channel with a series of cross-section geometries will be given and a definition of the hydraulic resistance will be introduced. The flow through channels will be simulated by use of FEMLAB. The results obtained in the FEMLAB simulation will be compared with those predicted by theory and discussed in Chapter 3. Finally, conclusions are given. 3
4 Chapter Theory.1 Poiseuille Flow in a circular pipe.1.1 Velocity field We consider an incompressible fluid flowing through a pipe (cylindrical tube) of length L in the x direction and a circular cross section with radius a in the yz plane; see Figure.1. The pressure drop over the length L is p, i.e p(0)=p 0 + p and p(l)=p 0. We introduce cylindrical polar coordinates in the yz-plane, y=rcosφ and z=rsinφ, where r is the radial distance from the cylinder axis and φ is the azimuthal angle (the angle around from the z axis). a Figure.1: A cylindrical tube with radius a. We study a flow along the axis of the tube, so that v = v x (r)e x = (v x (r), 0, 0). (.1) This field depends only on r and it is invariant under translations along the cylinder axis as well as rotations. 4
5 .1. POISEUILLE FLOW IN A CIRCULAR PIPE CHAPTER. THEORY In applying the Navier-Stokes equation 1 we assume that there is neither external forces nor instantaneous acceleration affecting the flow i.e f=0 and v t =0, respectively, giving p = η v (.3) where η is the viscosity of the fluid. Since v=e x v x (r) it follows from projecting of the Navier-Stokes equation p=η v on e r and e φ that p r = p = 0. (.4) φ This means that pressure only depends on the longitudinal coordinate p = p(x). Calculating the Laplacian term acting on a function of r in the Navier-Stokes equation we finally obtain ( ) dvx (r) v x (r) = ( v x (r)) = e x dr dp(x) dx = η ( d v x (r) dr + 1 r = d v x (r) dr + 1 r dv x dr (.5) ) dv x. (.6) dr The left hand side of this equation depends only on x whereas the right hand side depends only on r. Since the pressure difference over the length of the pipe is given by p L it follows that dp(x) dx = p (.7) L where the minus sign is due to the assumption that if the fluid is supposed to run in the positive x-direction the pressure must be higher at the beginning (x = 0) than at the end (x = L). Thus we obtain, with a little rewriting, the following Navier-Stokes equation equation ( 1 d r dv ) x(r) = η p r dr dr L. (.8) The solution to this equation with no-slip boundary conditions, i.e v x (0)=v x (a)=0, is ( ) v x (r) = v 0 1 r a (.9) where v 0 is the integration constant. This can easily be shown to be the solution by inserting dv x (r) = v 0r dr a (.10) 1 Navier-Stokes equation is a partial differential equation expressing the local balance of momentum in a fluid around any point of space at any time. It is given by ρ v t = p + η v + f (.) where the terms v is the instantaneous acceleration, p is the dynamical pressure gradient, t η v is the shear viscosity term and f is the external body force. 5
6 .. POISEUILLE FLOW IN AN ELLIPTIC PIPE CHAPTER. THEORY in the finally obtained Navier-Stokes. The integration constant is then determined to be v 0 = 1 p 4η L a. (.11) Thus, the solution to our Navier-Stokes equation is v x (r) = 1 ) p (1 4η L a r a. (.1).1. Flow rate The total flow rate Q (volume pr. time) carried through the pipe is given by Q = ρ dydz v x (r). (.13) circle In cylindrical polar coordinates it can be rewritten as to get Q = ρ a π 0 0 rdrdφ v x (r) (.14) Q = ρ p 8η L πa4. (.15). Poiseuille Flow in an elliptic pipe..1 Velocity field For an elliptic pipe, as shown in Figure., with an elliptic cross section having axes a and b along the y and the z axis, respectively, the same assumptions as for the circular pipe are valid. Instead of using the cylindrical polar coordinates we introduce the cartesian coordinates. The velocity field is again assumed longitudinal of the form v = v x (y, z)e x = (v x (y, z), 0, 0). (.16) The Navier-Stokes equation is given by the same equation as in a circular pipe. Choosing the cartesian coordinates the Navier-Stokes equation can be expressed as η v x (y, z) = p (.17) where the Laplacian term can be calculated by the same formalism of section.1 to give v x (y, z) = Thus we obtain the following Navier-Stokes equation 1 η ( ) v x (y, z) y + v x (y, z) z. (.18) ( p ) L = v x (y, z) y + v x (y, z) z 6 (.19)
7 .. POISEUILLE FLOW IN AN ELLIPTIC PIPE CHAPTER. THEORY a b Figure.: A pipe with an elliptic cross-section having axes a and b with the solution, satisfying the no-slip condition, given by Inserting v x(y,z) integration constant = v dy 0 a ( ) v x (y, z) = v 0 1 y a z b. (.0) and v x(y,z) dz = v 0 b v o = 1 η Thus, the longitudinal velocity profile in an elliptic pipe becomes v x (y, z) = 1 η p L in the Navier-Stokes equation we obtain the p L (ab) a + b (ab) a + b. (.1) ) (1 y a z b. (.) Notice that for a = b the integration constant for the elliptic pipe becomes the one for the circular pipe indicating that circular pipe is a special case of the elliptical... Flow rate The flow rate is given by Q = ρ dydz v x (y, z). (.3) ellipse Introducing the transformation variables Ỹ = 1 b z and Z= 1 ay, and inserting them in the integral for the flow rate we obtain Q = abv 0 ρ ellipse dỹ d Z 1 (1 ỹ z ) = πabv 0 rdr (1 r ) (.4) After the integration we finally obtain the flow rate for the elliptic pipe Q = ρ 4η p L 0 π(ab) 3 a + b (.5) 7
8 .3. POISEUILLE FLOW IN A RECTANGULAR PIPE - FOURIER TRANSFORMATION CHAPTER. THEORY.3 Poiseuille Flow in a rectangular pipe - Fourier Transformation.3.1 Velocity field The solution of the flow in the rectangular pipe is not simple to obtain. The no-slip boundary condition does not in this case allow for a separation of the D problem into two independent 1D problems. Instead we solve the problem by Fourier analysis, where only sine functions are involved since they guarantee fulfilment of the no-slip boundary condition. We Fourier-expand in the z-axis, so that v x (y, 0) = v x (y, h) = 0. The solution to the Navier-Stokes equation is of the form In applying the Navier-stokes equation v x (y, z) = f n (y) sin(nπ z h ) (.6) n=1 v x (y, z) = P ηl (.7) we get v x (y, z) = [ d f n (y) n=1 dy n π h ] f n(y) sin(nπ z h ) (.8) In the Navier-Stokes equation we also have to Fourier transform the constant P/ηL, or simpler just the constant 1 which can be multiplied by any constant. The Fourier transform of the constant P/ηL thus becomes The solution is only possible if P ηl = P ηl 1 = P 4 ηl π n,odd 1 n sin(nπ z ). (.9) h f n (y) = 0 for n even (.30) d f n (y) dy n π h f n(y) = P 4 for n odd (.31) ηl nπ The first step to solve this nonhomogeneous equation is to find a general solution of the homogeneous equation f n (y) h. The characteristic equation of the homogeneous equation is λ (nπ/h) λ = 0. (.3) It has the roots λ 1 = 0 and λ = nπ/h. Hence a real general solution of the homogenous equation is f n (y) h = Ae yλ 1 + Be yλ. (.33) 8
9 .3. POISEUILLE FLOW IN A RECTANGULAR PIPE - FOURIER TRANSFORMATION CHAPTER. THEORY W H Figure.3: A rectangular pipe of height h and width w. Using the definitions of the hyperbolic functions sinh x = ex e x cosh x = ex + e x the general solution of the homogeneous equation can be rewritten as ( f n (y) h = B cosh nπ y ) ( + B sinh nπ y ) h h (.34) (.35) (.36) The second step is to find a particular solution f n (y) p of the nonhomogeneous equation. Since the equation on the right has the derivative zero, we say f n (y) p = C. (.37) Substitution gives (n π h ) C = P 4 P ηl nπ. By comparison, C = ηl solution of the nonhomogeneous equation is 4h π 3 1 n 3. Hence a general f n (y) = f n (y) h + f n (y) p (.38) ( = B cosh nπ y ) ( + B sinh nπ y ) P 4h 1 h h η π 3 n 3. (.39) From the boundary condition f n ( ω/) = f n (ω/) = 0, we get B = P η resulting in [ ( ) ] cosh f n (y) = P η 4 h π 3 n 3 cosh nπ y z ( ) 1 nπ w h 4h π 3 1 n 3 1 cosh(nπ w h ). (.40) Thus the solution of the flow in the rectangular pipe is v x (y, z) = P 4 [ ( ) 1 cosh nπ y ] ( η h z π 3 n 3 ( ) 1 sin nπ y ). (.41) n,odd cosh nπ w z h 9
10 .4. HYDRAULIC RESISTANCE CHAPTER. THEORY.3. Flow rate The flow rate for a rectangular pipe is given by Using the following we get w 0 h 0 w h Q = ρ dy dz v x (y, z). (.4) 0 0 ( dy cosh nπ y ) h ( dz sinh nπ z ) h n,odd Q = P [ 1η ρwh3 1 = h ( nπ sinh nπ y ) (.43) h = h for n odd (.44) nπ 1 n 4 = π4 96. (.45) n,odd 19 π 5 h w 1 (nπ n 5 tanh w ) ]. (.46) h.4 Hydraulic resistance In analogy with Ohm s law for an electric resistor, it is natural to define the hydraulic resistance R hyd of a pipe carrying some liquid as p = R hyd Q (.47) Since we are interested in the relation between the hydraulic resistance and the crosssection geometry of the pipe it is appropriate to reexpress the hydraulic resistance as R = α η ρa (.48) where α is the correction coefficient. We introduce the compactness β which the measure of the square of the perimeter l over the area A of the cross-section geometry Ω, β = l A = ( Ω ) ( 1 dl dr) (.49) Ω.4.1 Elliptic cross-section Expressing the correction coefficient as a function of the compactness we get for the circular cross-section that α(β) = 8π = β. 10
11 .4. HYDRAULIC RESISTANCE CHAPTER. THEORY Dd Figure.4: A domain ω. subdomain is the area A. The red boundary is the perimeter l of the domain and the blue For the elliptic cross-section the compactness is given by ( ) π 0 1 dθ (1 b /a ) sin θ β = (.50) πab = 16 ( π a dθ 1 (1 b π b /a ) sin θ) (.51) 0 16 a, a b. (.5) π b Thus, the correction coefficient for the elliptic cross-section is given by α(β) = 4π a + b ab π β. (.53) 4.4. Rectangular cross-section ( ) Assuming that w > h then tanh nπ w h 1 and we get the flow rate to Q P ( 1η ρwh ) π 5 ζ(5)w. (.54) h The correction coefficient is then given by α 1π 5 γ π 5 γ 186ζ(5), (.55) where γ is the width-to-hight ratio w/h. Setting the compactness β = 8 + 4γ + 4/γ for a rectangle in the expression of α we get α = From the expression of α it can be shown that for 1π 5 (β/4 1/γ) π 5 (β/4 1/γ) 186ζ(5). (.56) α lim β β = 3 (.57) 11
12 .4. HYDRAULIC RESISTANCE CHAPTER. THEORY we have the asymptotic dependence α 3β + c. Since 3ζ(5) lim (α 3β) = 4 + β π 5 (.58) we get the correction coefficient for the rectangular class to α = 3β (.59) 1
13 Chapter 3 FEMLAB simulations FEMLAB is a partial differential equation solver that builds on MATLAB. FEMLAB - which stands for Finite Element Modelling LABoratory- is an advanced software package for modeling and simulation of any physical process that can be described with partial differential equations. The differential equation must be specified in some geometry, the so-called computational domain, together with the boundary conditions along the surface of the domain as well as the initial values at the starting time. The computational domain is divided into a huge number of finite-sized elements forming the basis computation. Having introduced the analytic description of the Poiseuille flow through a pipe with various cross-sections, we are now interested in simulating the flow in FEMLAB for the same cross-section geometries. The numerical calculations are then compared to those predicted by the theory. Before starting the simulation, a differential equation must be specified. Since the equation expressing the velocity field is of the form P = η v, the the simulation differential model was chosen to be a Piosson s equation in the D-space. To draw the geometry in which the partial differential equation must be specified, we click on Draw Mode on the tool bar. In this window the desired geometry can now be designed. Hereafter, the no-slip boundary is specified by clicking on Boundary Settings under Physics. Specifying the subdomain settings is done by clicking on Subdomain Settings under Physics. The geometry is then refined and the PDE solved. The solution plotted, displays the velocity field; see Figure 3.1. The velocity field of the flow has its maximum in the center of the pipe, but begins to decrease when moving radially from the center. When in contact with the cross-section boundary the velocity of the flow is at its minimum, fulfilling the no-slip condition. Since we also are interested in calculating the the hydraulic resistance, the correction coefficient 1 and the compactness is determined by making subdomain A, boundary l and velocity field integration. The results depicting the correction coefficient versus the com- 1 We saw that the hydraulic resistance was defined as P/L = R hyd Q. This can be rewritten as P = η α ρ dydzv L ρ A x (y, z). Since all the values are required dimensionless in FEMLAB we set the factor A P/Lη = 1 giving us the correction coefficient: α = dydzv x(y,z). 13
14 CHAPTER 3. FEMLAB SIMULATIONS Figure 3.1: A diagram showing the simulation results of the velocity field for a circular, elliptical, rectangular and triangular cross-section. pactness for various cross-sections is shown in Figure 3.. It emerges clearly from Figure 3. that the relation between the correction coefficient and the compactness is linear for the three cross-section classes: circular, elliptic and rectangular. The numeric proportionality factor for all cross-section groups is in good agreement with those predicted by theory. For the elliptic it is.55, close to the theoretical value of.47 The deviation is due to the assumption: a b, made in our calculation. For the rectangular class, the value of the proportionality factor is found to be 3. This value agree completely with the theory. Since the analytic calculations are not quite possible for triangle classes, only simulation was made. The results for this cross-section class is also shown in Figure 3.. Also here, the relation between α and β is linear. The proportionality factor is less than the other cross-section classes, namely For cross-section classes more complicated than the elliptical (circular), α increases for small values of β, leading to an increase of the hydraulic resistance. 14
15 CHAPTER 3. FEMLAB SIMULATIONS Correction coefficient vs. Compactness Theory for rectangular class. Theory for ellipse class. Linear fit to triangular class. 100 For circular l /A Figure 3.: Correction coefficient versus compactness for the elliptical, rectangular and triangular class. The result for the circle is also indicated. 15
16 Chapter 4 Conclusion We have demonstrated a close to linear relation between the hydraulic resistance and the compactness for the following cross-section geometries: circular, elliptic, triangular and rectangular. The results obtained were in good agreement with theory. For geometries more complicated than the circular the resistance increases, since a big part of the fluid is in contact with the inner surface of the channel. 16
Convection-Diffusion in Microchannels
33443 Microfluidics Theory and Simulation Convection-Diffusion in Microchannels Catharina Hartmann (s97176) Luigi Sasso (s06094) Institute of Micro and Nanotechnology, DTU, 800 Kongens Lyngby, DK June,
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 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 informationENGI Multiple Integration Page 8-01
ENGI 345 8. Multiple Integration Page 8-01 8. Multiple Integration This chapter provides only a very brief introduction to the major topic of multiple integration. Uses of multiple integration include
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationASTR 320: Solutions to Problem Set 2
ASTR 320: Solutions to Problem Set 2 Problem 1: Streamlines A streamline is defined as a curve that is instantaneously tangent to the velocity vector of a flow. Streamlines show the direction a massless
More informationAssignment 2 Continuous Matter 4335/8191(B)
Assignment 2 Continuous Matter 4335/8191(B Tyler Shendruk October 29, 2012 1 Problem 1 Given the stream function ψ = x x 3 y/2, determine if mass is conserved the flow is irrotational and the velocity
More informationFluid mechanics and living organisms
Physics of the Human Body 37 Chapter 4: In this chapter we discuss the basic laws of fluid flow as they apply to life processes at various size scales For example, fluid dynamics at low Reynolds number
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
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 informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
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 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 informationNote: Each problem is worth 14 points except numbers 5 and 6 which are 15 points. = 3 2
Math Prelim II Solutions Spring Note: Each problem is worth points except numbers 5 and 6 which are 5 points. x. Compute x da where is the region in the second quadrant between the + y circles x + y and
More informationInjection Molding. Figure 1: Principles of injection molding. Injection molding cycle: part solidifies. Open Mold Eject Part Close Mold
Injection Molding Figure 1: Principles of injection molding. Injection molding cycle: Extruder Pressure Extrude Mold Inject Pack Solidify Open Mold Eject Part Close Mold gate solidifies part solidifies
More informationCreated by T. Madas LINE INTEGRALS. Created by T. Madas
LINE INTEGRALS LINE INTEGRALS IN 2 DIMENSIONAL CARTESIAN COORDINATES Question 1 Evaluate the integral ( x + 2y) dx, C where C is the path along the curve with equation y 2 = x + 1, from ( ) 0,1 to ( )
More informationDUE: WEDS MARCH 26TH 2018
HOMEWORK # 2: FINITE DIFFERENCES MAPPING AND TWO DIMENSIONAL PROBLEMS DUE: WEDS MARCH 26TH 2018 NOTE: In this homework, you will choose ONE of the following three questions to perform and hand-in. Each
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 informationImplementation of TOPOPT in a Microfluidic Heat Transfer Model
Microfluidics Theory and Simulation, c33443 Implementation of TOPOPT in a Microfluidic Heat Transfer Model Kristian Krüger Grøndahl, (s011426) and Thomas Glasdam Jensen, (s011427) Supervisor: Henrik Bruus
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 informationMULTIVARIABLE INTEGRATION
MULTIVARIABLE INTEGRATION (PLANE & CYLINDRICAL POLAR COORDINATES) PLANE POLAR COORDINATES Question 1 The finite region on the x-y plane satisfies 1 x + y 4, y 0. Find, in terms of π, the value of I. I
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 informationMathematical Methods and its Applications (Solution of assignment-12) Solution 1 From the definition of Fourier transforms, we have.
For 2 weeks course only Mathematical Methods and its Applications (Solution of assignment-2 Solution From the definition of Fourier transforms, we have F e at2 e at2 e it dt e at2 +(it/a dt ( setting (
More informationKinematics (2) - Motion in Three Dimensions
Kinematics (2) - Motion in Three Dimensions 1. Introduction Kinematics is a branch of mechanics which describes the motion of objects without consideration of the circumstances leading to the motion. 2.
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More informationProblem Solving 1: Line Integrals and Surface Integrals
A. Line Integrals MASSACHUSETTS INSTITUTE OF TECHNOLOY Department of Physics Problem Solving 1: Line Integrals and Surface Integrals The line integral of a scalar function f ( xyz),, along a path C is
More informationMathematics (Course B) Lent Term 2005 Examples Sheet 2
N12d Natural Sciences, Part IA Dr M. G. Worster Mathematics (Course B) Lent Term 2005 Examples Sheet 2 Please communicate any errors in this sheet to Dr Worster at M.G.Worster@damtp.cam.ac.uk. Note that
More informationPHY321 Homework Set 10
PHY321 Homework Set 10 1. [5 pts] A small block of mass m slides without friction down a wedge-shaped block of mass M and of opening angle α. Thetriangular block itself slides along a horizontal floor,
More informationConsider a volume Ω enclosing a mass M and bounded by a surface δω. d dt. q n ds. The Work done by the body on the surroundings is.
The Energy Balance Consider a volume Ω enclosing a mass M and bounded by a surface δω. δω At a point x, the density is ρ, the local velocity is v, and the local Energy density is U. U v The rate of change
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 informationc2 2 x2. (1) t = c2 2 u, (2) 2 = 2 x x 2, (3)
ecture 13 The wave equation - final comments Sections 4.2-4.6 of text by Haberman u(x,t), In the previous lecture, we studied the so-called wave equation in one-dimension, i.e., for a function It was derived
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 information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More information!! +! 2!! +!"!! =!! +! 2!! +!"!! +!!"!"!"
Homework 4 Solutions 1. (15 points) Bernoulli s equation can be adapted for use in evaluating unsteady flow conditions, such as those encountered during start- up processes. For example, consider the large
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 informationSupplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer
Supplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer Thomas P. Forbes and Jason G. Kralj National Institute of Standards and Technology,
More information4. Analysis of heat conduction
4. Analysis of heat conduction John Richard Thome 11 mars 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Conduction 11 mars 2008 1 / 47 4.1 The well-posed problem Before we go further with
More information14. Rotational Kinematics and Moment of Inertia
14. Rotational Kinematics and Moment of nertia A) Overview n this unit we will introduce rotational motion. n particular, we will introduce the angular kinematic variables that are used to describe the
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 informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More information1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is
1. The value of the double integral (a) 15 26 (b) 15 8 (c) 75 (d) 105 26 5 4 0 1 1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is 2. What is the value of the double integral interchange the order
More informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationProblems in Magnetostatics
Problems in Magnetostatics 8th February 27 Some of the later problems are quite challenging. This is characteristic of problems in magnetism. There are trivial problems and there are tough problems. Very
More information14.1. Multiple Integration. Iterated Integrals and Area in the Plane. Iterated Integrals. Iterated Integrals. MAC2313 Calculus III - Chapter 14
14 Multiple Integration 14.1 Iterated Integrals and Area in the Plane Objectives Evaluate an iterated integral. Use an iterated integral to find the area of a plane region. Copyright Cengage Learning.
More informationWeek 7: Integration: Special Coordinates
Week 7: Integration: Special Coordinates Introduction Many problems naturally involve symmetry. One should exploit it where possible and this often means using coordinate systems other than Cartesian coordinates.
More informationLECTURE NOTES - III. Prof. Dr. Atıl BULU
LECTURE NOTES - III «FLUID MECHANICS» Istanbul Technical University College of Civil Engineering Civil Engineering Department Hydraulics Division CHAPTER KINEMATICS OF FLUIDS.. FLUID IN MOTION Fluid motion
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 informationMAC2313 Final A. (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative.
MAC2313 Final A (5 pts) 1. How many of the following are necessarily true? i. The vector field F = 2x + 3y, 3x 5y is conservative. ii. The vector field F = 5(x 2 + y 2 ) 3/2 x, y is radial. iii. All constant
More informationCurvilinear coordinates
C Curvilinear coordinates The distance between two points Euclidean space takes the simplest form (2-4) in Cartesian coordinates. The geometry of concrete physical problems may make non-cartesian coordinates
More informationSOLUTIONS TO HOMEWORK ASSIGNMENT #2, Math 253
SOLUTIONS TO HOMEWORK ASSIGNMENT #, Math 5. Find the equation of a sphere if one of its diameters has end points (, 0, 5) and (5, 4, 7). The length of the diameter is (5 ) + ( 4 0) + (7 5) = =, so the
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 4 HEAT TRANSFER IN CHANNEL FLOW BASIC CONCEPTS BASIC CONCEPTS Laminar
More informationf dr. (6.1) f(x i, y i, z i ) r i. (6.2) N i=1
hapter 6 Integrals In this chapter we will look at integrals in more detail. We will look at integrals along a curve, and multi-dimensional integrals in 2 or more dimensions. In physics we use these integrals
More informationStaple or bind all pages together. DO NOT dog ear pages as a method to bind.
Math 3337 Homework Instructions: Staple or bind all pages together. DO NOT dog ear pages as a method to bind. Hand-drawn sketches should be neat, clear, of reasonable size, with axis and tick marks appropriately
More informationME 509, Spring 2016, Final Exam, Solutions
ME 509, Spring 2016, Final Exam, Solutions 05/03/2016 DON T BEGIN UNTIL YOU RE TOLD TO! Instructions: This exam is to be done independently in 120 minutes. You may use 2 pieces of letter-sized (8.5 11
More informationEELE 3331 Electromagnetic I Chapter 3. Vector Calculus. Islamic University of Gaza Electrical Engineering Department Dr.
EELE 3331 Electromagnetic I Chapter 3 Vector Calculus Islamic University of Gaza Electrical Engineering Department Dr. Talal Skaik 2012 1 Differential Length, Area, and Volume This chapter deals with integration
More informationMA 441 Advanced Engineering Mathematics I Assignments - Spring 2014
MA 441 Advanced Engineering Mathematics I Assignments - Spring 2014 Dr. E. Jacobs The main texts for this course are Calculus by James Stewart and Fundamentals of Differential Equations by Nagle, Saff
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2013
Lecture 1 3/13/13 University of Washington Department of Chemistry Chemistry 53 Winter Quarter 013 A. Definition of Viscosity Viscosity refers to the resistance of fluids to flow. Consider a flowing liquid
More informationName: Answer Key David Arnold. Math 50B Integral Calculus May 13, Final Exam
Math 5B Integral Calculus May 3, 7 Final Exam Name: Answer Key David Arnold Instructions. (9 points) Follow the directions exactly! Whatever you are asked to do, you must do to receive full credit for
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationMATH 18.01, FALL PROBLEM SET # 6 SOLUTIONS
MATH 181, FALL 17 - PROBLEM SET # 6 SOLUTIONS Part II (5 points) 1 (Thurs, Oct 6; Second Fundamental Theorem; + + + + + = 16 points) Let sinc(x) denote the sinc function { 1 if x =, sinc(x) = sin x if
More informationLecture 10 - Moment of Inertia
Lecture 10 - oment of Inertia A Puzzle... Question For any object, there are typically many ways to calculate the moment of inertia I = r 2 dm, usually by doing the integration by considering different
More informationLecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018
Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent
More informationAP Calculus BC Chapter 4 AP Exam Problems A) 4 B) 2 C) 1 D) 0 E) 2 A) 9 B) 12 C) 14 D) 21 E) 40
Extreme Values in an Interval AP Calculus BC 1. The absolute maximum value of x = f ( x) x x 1 on the closed interval, 4 occurs at A) 4 B) C) 1 D) 0 E). The maximum acceleration attained on the interval
More informationHeat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection
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 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 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 informationPractice Problems for the Final Exam
Math 114 Spring 2017 Practice Problems for the Final Exam 1. The planes 3x + 2y + z = 6 and x + y = 2 intersect in a line l. Find the distance from the origin to l. (Answer: 24 3 ) 2. Find the area of
More informationPRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.
PRACTICE PROBLEMS Please let me know if you find any mistakes in the text so that i can fix them. 1.1. Let Show that f is C 1 and yet How is that possible? 1. Mixed partial derivatives f(x, y) = {xy x
More informationFurther Linear Elasticity
Torsion of cylindrical bodies Further Linear Elasticity Problem Sheet # 1. Consider a cylindrical body of length L, the ends of which are subjected to distributions of tractions that are statically equivalent
More informationxy 2 e 2z dx dy dz = 8 3 (1 e 4 ) = 2.62 mc. 12 x2 y 3 e 2z 2 m 2 m 2 m Figure P4.1: Cube of Problem 4.1.
Problem 4.1 A cube m on a side is located in the first octant in a Cartesian coordinate system, with one of its corners at the origin. Find the total charge contained in the cube if the charge density
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationFigure 21:The polar and Cartesian coordinate systems.
Figure 21:The polar and Cartesian coordinate systems. Coordinate systems in R There are three standard coordinate systems which are used to describe points in -dimensional space. These coordinate systems
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 informationChapter 1: Concept of similarity. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 1: Concept of similarity Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 1: Concept of similarity A few definitions Geometrical similarity Fractal similarity Physical
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationG G. G. x = u cos v, y = f(u), z = u sin v. H. x = u + v, y = v, z = u v. 1 + g 2 x + g 2 y du dv
1. Matching. Fill in the appropriate letter. 1. ds for a surface z = g(x, y) A. r u r v du dv 2. ds for a surface r(u, v) B. r u r v du dv 3. ds for any surface C. G x G z, G y G z, 1 4. Unit normal N
More informationy=1 1 J (a x ) dy dz dx dz 10 4 sin(2)e 2y dy dz sin(2)e 2y
Chapter 5 Odd-Numbered 5.. Given the current density J = 4 [sin(x)e y a x + cos(x)e y a y ]ka/m : a) Find the total current crossing the plane y = in the a y direction in the region
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationFinal Review Worksheet
Score: Name: Final Review Worksheet Math 2110Q Fall 2014 Professor Hohn Answers (in no particular order): f(x, y) = e y + xe xy + C; 2; 3; e y cos z, e z cos x, e x cos y, e x sin y e y sin z e z sin x;
More informationNo calculators, cell phones or any other electronic devices can be used on this exam. Clear your desk of everything excepts pens, pencils and erasers.
Name: Section: Recitation Instructor: READ THE FOLLOWING INSTRUCTIONS. Do not open your exam until told to do so. No calculators, cell phones or any other electronic devices can be used on this exam. Clear
More informationENGI Partial Differentiation Page y f x
ENGI 3424 4 Partial Differentiation Page 4-01 4. Partial Differentiation For functions of one variable, be found unambiguously by differentiation: y f x, the rate of change of the dependent variable can
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 informationFlow past a slippery cylinder
Faculty of Mathematics, University of Waterloo, Canada EFMC12, September 9-13, 2018, Vienna, Austria Problem description and background Conformal mapping Boundary conditions Rescaled equations Asymptotic
More informationMath Review for Exam Compute the second degree Taylor polynomials about (0, 0) of the following functions: (a) f(x, y) = e 2x 3y.
Math 35 - Review for Exam 1. Compute the second degree Taylor polynomial of f e x+3y about (, ). Solution. A computation shows that f x(, ), f y(, ) 3, f xx(, ) 4, f yy(, ) 9, f xy(, ) 6. The second degree
More informationFluid Mechanics Answer Key of Objective & Conventional Questions
019 MPROVEMENT Mechanical Engineering Fluid Mechanics Answer Key of Objective & Conventional Questions 1 Fluid Properties 1. (c). (b) 3. (c) 4. (576) 5. (3.61)(3.50 to 3.75) 6. (0.058)(0.05 to 0.06) 7.
More informationMath 10C - Fall Final Exam
Math 1C - Fall 217 - Final Exam Problem 1. Consider the function f(x, y) = 1 x 2 (y 1) 2. (i) Draw the level curve through the point P (1, 2). Find the gradient of f at the point P and draw the gradient
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 informationContinuum Mechanics Examination exercises
Continuum Mechanics Examination exercises Name:... Result: 1:1... 1:... 1:3... :1... :... :3... :4... 3:1... 3:... 3:3... 4:1... 4:... 5:1... 5:... LTH, Box 118, S-1 00 LUND Sweden Tel: +46 (0)46 3039
More informationReview Sheet for the Final
Review Sheet for the Final Math 6-4 4 These problems are provided to help you study. The presence of a problem on this handout does not imply that there will be a similar problem on the test. And the absence
More informationDavid Abrecht. February 17, 2012
Addendum to An analytical solution method for the unsteady, unbounded, incompressible three-dimensional Navier-Stokes equations in Cartesian coordinates using coordinate axis symmetry degeneracy David
More informationENGI Duffing s Equation Page 4.65
ENGI 940 4. - Duffing s Equation Page 4.65 4. Duffing s Equation Among the simplest models of damped non-linear forced oscillations of a mechanical or electrical system with a cubic stiffness term is Duffing
More informationAE301 Aerodynamics I UNIT B: Theory of Aerodynamics
AE301 Aerodynamics I UNIT B: Theory of Aerodynamics ROAD MAP... B-1: Mathematics for Aerodynamics B-: Flow Field Representations B-3: Potential Flow Analysis B-4: Applications of Potential Flow Analysis
More informationMcGill University December Intermediate Calculus. Tuesday December 17, 2014 Time: 14:00-17:00
McGill University December 214 Faculty of Science Final Examination Intermediate Calculus Math 262 Tuesday December 17, 214 Time: 14: - 17: Examiner: Dmitry Jakobson Associate Examiner: Neville Sancho
More informationAppendix A Vector Analysis
Appendix A Vector Analysis A.1 Orthogonal Coordinate Systems A.1.1 Cartesian (Rectangular Coordinate System The unit vectors are denoted by x, ŷ, ẑ in the Cartesian system. By convention, ( x, ŷ, ẑ triplet
More informationPressure Losses for Fluid Flow Through Abrupt Area. Contraction in Compact Heat Exchangers
Pressure Losses for Fluid Flow Through Abrupt Area Contraction in Compact Heat Exchangers Undergraduate Research Spring 004 By Bryan J. Johnson Under Direction of Rehnberg Professor of Ch.E. Bruce A. Finlayson
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 information