Flow and Transport. c(s, t)s ds,

Transcription

1 Flow and Transport 1. The Transport Equation We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section S. (The very same discussion applies to the description of the transport of gas by air moving through a pipe.) We identify G with the open interval (a, b), and the velocity ν > 0 is in the (rightward) positive direction of the -ais. Since the tube is very thin, we can assume that the concentration of the chemical is constant across the cross section S at each point G. Let c(, t) denote this concentration within the tube at a point G and at time t > 0. The Conservation Law. The amount of chemical stored in the tube within a section [, + h] of length h > 0 is given by +h c(s, t)s ds, The flu q(, t) at the point is the mass flow rate of the chemical to the right per unit area, Equating the rate at which the chemical is stored within the section [, + h] to the rate at which it flows into the section plus the rate at which the chemical is generated within this section, we arrive at d dt +h c(s, t)s ds = S ( q(, t) q( + h, t) ) + +h F (s, t)s ds where F (, t) represents the rate at which chemical is generated per unit volume. This source term is assumed to be a known function of space and time. We differentiate the integral on the left side and write the difference on the right side as an integral of a derivative to obtain +h (s, t) S ds = S +h 1 ( q(s, t) + F (s, t) ) ds.

2 2 FLOW AND TRANSPORT Dividing this by Sh and letting h 0 yield the conservation law (1) + q = F (, t). The transport equation. The flu q(, t) at the point is given by (2) q(, t) = νc(, t). This is the mass flow rate of the chemical due to convection, the direct transport by the moving water. Substituting this into the conservation law (1) yields the one dimensional transport equation (3) + ν = F (, t). This is also known as the first-order wave equation. It is just the differential form of the conservation equation (1) combined with the constitutive equation (2) of convection. It is assumed here that the velocity ν is sufficiently large that we can ignore the comparatively smaller effects of diffusion, i.e., the natural motion of the chemical from areas of high concentration to those of lower concentration. This will be a major topic later. Initial and Boundary conditions. Since the transport equation is first order in space and time, one may epect that in order to have a well posed problem, one boundary condition and one initial condition should be specified. We shall see that this is true here. We want to find a solution of (3) which satisfies an initial condition of the form c(, 0) = c 0 (), a < < b, where c 0 ( ) is given. That is, we need to specify the initial state of the model. Furthermore we shall specify the value of the concentration at the left end point where the substance is entering the tube: c(a, t) = c a (t), t > 0, where c a ( ) is given. This type of boundary condition arises when the value of the concentration at the end point is known, usually from a direct measurement. Such a condition arises when one sets the boundary concentration to a prescribed value, for eample, c a (t) = 0 when only pure water with no chemical is entering the left end of the tube.

3 1. THE TRANSPORT EQUATION 3 The Solution: introduction of characteristic curve. Consider the homogeneous transport equation c = c(, t) : c t + ν c = 0. If we make a change of variable to the new coordinates ξ, τ which are defined implicitly by = ξ + τ, ν t = ξ τ, and use the chain rule, ξ = ξ + ξ, we find that the transport equation is equivalent to c = c(ξ, τ) : ξ = 0. This shows that the solutions are given by c(ξ, τ) = f(τ), so the general solution of the transport equation is of the form c(, t) = f( ν t) for some function f( ). Now let s choose this function so that c(, ) satisfies the initial and boundary conditions. The initial condition requires c(, 0) = f() = c 0 (), a, and the boundary condition likewise requires c(a, t) = f(a ν t) = c a (t), t 0, so we must have f(s) = c a ( a s ν ), s a. Thus, the solution of the initialboundary-value problem is given by { c 0 ( ν t), ν t 0, c(, t) = c a ( a +ν t ν ), a ν t + a, Note that the important fact behind these calculations is that the solution was constant along the curves where τ is constant, i.e., along the curves where ν t is constant. It followed from this that the solution is a pure translation to the right with velocity ν, certainly no surprise in view of the origin of the transport equation. These special curves are the

4 4 FLOW AND TRANSPORT characteristic curves for the transport equation, and they will arise in our discussions of first order equations. Eercise 1. Suppose that c a ( ) and c 0 ( ) are continuous. Show that the solution of the initial-boundary-value problem is continuous if and only if c a (0 + ) = c 0 (a + ). Eercise 2. Find the solution of the initial-value problem for the nonhomogeneous transport equation (3) on the region {(, t) : t 0} with the initial condition c(, 0) = c 0 () on < < The Porous Medium Equation We shall describe the diffusion of fluid through the long thin tube G with uniform cross section S. As before, we identify G with the open interval (a, b), but here we assume that the tube contains a distribution of particles which impede the flow of the fluid. The fluid is constrained to flow in the complementary region of open channels and pores not occupied by the particles. This is the case for flow of fluids through soil or through any other such porous and permeable material. The conservation law. Let ρ(, t) denote the density of the fluid within the tube at a point G and at time t > 0. The mass of fluid stored in the section [, + h] of length h > 0 is given by +h φ(s)ρ(s, t)s ds, where φ() is the porosity of the porous medium at, i.e., the volume fraction of the medium occupied by the fluid. The fluid velocity ν(, t) at the point is the flow rate to the right per unit area measured in volume of fluid per unit area per time. Equating the rate at which fluid is stored within the section to the rate at which fluid flows into the section plus the fluid source rate within this section, we arrive at the fluid conservation equation for the section [, + h], +h φ(s)ρ(s, t)s ds = S ( ρ(, t)ν(, t) ρ( + h, t)ν( + h, t) ) + +h f(s, t)s ds,

5 2. THE POROUS MEDIUM EQUATION 5 where f(, t) represents the rate at which fluid is inserted per unit volume. This source term is assumed to be a known function of space and time. Dividing by Sh and letting h 0 yields the fluid conservation equation ρ(, t) (ρ(, t)ν(, t)) (4) φ() + = f(, t). Darcy s law. Let p(, t) denote the pressure of the fluid in the pores. This is measured in force per unit area. A fundamental eperimental observation in flow through porous media is that the fluid velocity ν(, t) at the point is proportional to the pressure gradient, (5) ν(, t) = k() p(, t) µ. The constant µ is the viscosity of the fluid, a measure of its resistance to shear, and this equation defines the permeability k() of the porous medium at the point (a, b). It is a measure of the conductivity of the medium, i.e., the inverse of resistance of the medium to internal flow. Since fluid flows in the direction of decreasing pressure, the minus sign is appropriate. In fact, if we write this in the form µ t) ν(, t) = p(, k(), it is a balance of forces on the fluid as it flows through the medium, and the coefficient µ/k() is the resistance to flow. By substituting Darcy s law (5) into the energy conservation law (1), we obtain the one dimensional porous medium equation (6) φ() ρ ( ρ k() µ ) p = f(, t). It remains to specify the state equation, the relation between density ρ and pressure p for the particular fluid. If the fluid is slightly compressible, then we write ρ = ρ 0 e c(p p0) where the constant c is the compressibility of the fluid. In this case the chain rule shows that ρp = 1 c ρ, so we obtain the linear diffusion equation for fluid density (7) cφ() ρ ( k() µ ρ ) = c f(, t).

6 6 FLOW AND TRANSPORT If the fluid is incompressible, i.e., if c = 0, then the density is constant and we obtain the equation (8) ( ) k() p ρ 0 = f(, t). µ Of course, either of these must be supplemented with appropriate initial and boundary conditions to get a well posed problem which determines the density and pressure along the length of the tube, and then these determine the Darcy velocity (5). This may be then used as data for the transport equation (3). 3. Dynamics of Chemical Adsorption We consider net the effect of an adsorbing medium that is distributed throughout the length of the tube. This material occurs as a distribution of small particles or fibres, such as a porous medium. In equilibrium conditions, the concentration a() of the chemical that is adsorbed onto the material at a point is a specific function of the concentration d() within the pores of the material. The quantities a and d are related by a given function a = a(d) called the adsorption isotherm. That is, a(d) is the concentration of chemical adsorbed in the medium which is at the concentration d. In dynamic conditions, the concentration d in the medium can be different from the concentration c in the surrounding water. The relation of the adsorbed concentration a(, t) to the eternal water concentration c(, t) is given by a(d) (9) = β(c d), where d(, t) is the internal concentration of the chemical in the medium. The relation (9) is the equation of kinetic adsorption, and it describes the dynamic response of the adsorbed concentration to a changing concentration in the surrounding pores. The constant β is called the kinetic coefficient of the process. The simplest form for the adsorption isotherm is the linear relationship a(d) = α d called the Henry isotherm in which the coefficient α is called the Henry coefficient. It is a reasonable assumption in regions for which

7 3. DYNAMICS OF CHEMICAL ADSORPTION 7 the variations of concentrations remain small. adsorption equation then takes the form α d + β(d c) = 0, for which the solution is given eplicitly by In this case, the kinetic d(t) = e β α t d(0) + t 0 β α e β α (t τ) c(τ) dτ. Thus, the internal medium concentration d follows the true pore concentration eponentially with time constant β α. In particular, in the limiting case β, we have instantaneous equilibrium and d c. Eercise 3. Note that the kernel appearing above is of the form k(t τ), where the function k(t) = β α e β α (t) satisfies (1) k(t) 0, (2) 0 k(t) dt = 1. Sketch k( ) for very large values of β. Transport with Adsorption. Consider again the transport of a chemical by water that is traveling with uniform velocity ν through a long thin tube G with uniform cross section S as before, but now the tube is filled with the adsorbing material. We let c(, t) be the concentration of chemical in the pores surounding the adsorbing material and a(, t) the concentration of chemical that is adsorbed onto the material. The chemical in the pores is carried by the water, but the adsorbed chemical is held stationary at the adsorption site. Thus, the flu is given eactly as before, but the total concentration of chemical at G is given by c(, t) + a(, t), so the conservation law leads to the transport equation with adsorption, (10) ( ) c + a(d) + ν = F (, t). In the case of instantaneous adsorption rates, β, we have c = d in the limit, and the adsorption equation becomes (11) ( ) c + a(c) + ν = F.

8 8 FLOW AND TRANSPORT This is a quasi-linear partial differential equation of first-order. More generally, the model of transport with kinetic adsorption leads to the system (12a) + ν + β(c d) = F, a(d) (12b) + β(d c) = 0, for the two unknown functions c(, t), d(, t). The appropriate initial and boundary conditions for this system are c(, 0) = c 0 (), d(, 0) = d 0, a < < b, c(a, t) = c a (t), t > 0. Note that the condition on d(, 0) is equivalent to specifying a corresponding condition on the adsorbed concentration, a(, 0). In the linear case of the Henry isotherm, we obtain the linear system of first order partial differential equations (13a) + ν + β(c d) = F, (13b) α d + β(d c) = 0. A Second order Wave Equation. We continue with an additional assumption, namely, that the concentration of chemical in the pores is much less than the amount adsorbed in the material. This amounts to assuming that most of the available volume in the tube is occupied by the adsorbing material, and it means that c << a, and so we approimate c+a a. Thus, we drop c from the sum c+a in the conservation equation to obtain in the linear case (14a) ν + β(c d) = F, α d (14b) + β(d c) = 0. Eercise 4. Eliminate d above to obtain a single equation of second order for the unknown c: 1 (15) β c t + 1 ν c t + 1 α c = 1 να F + 1 νβ F t. What are the appropriate initial and boundary conditions for this equation?

9 4. LONGITUDINAL VIBRATIONS 9 Eercise 5. Eliminate d from the system (13) to obtain a single equation of second order for the unknown c. What are the appropriate initial and boundary conditions for this equation? 4. Longitudinal Vibrations We describe the longitudinal vibrations in a long narrow cylindrical rod of cross section area S. The rod is located along the ais, and we identify it with the interval (a, b) in R. The rod is assumed to stretch or contract in the horizontal direction, and we assume that the vertical plane cross-sections of the rod move only horizontally. Denote by u(, t) the displacement in the positive direction from the point [a, b] at the time t > 0. The corresponding displacement rate or velocity is denoted by v(, t) u t (, t). Let σ(, t) denote the local stress, the force per unit area with which the part of the rod to the right of the point acts on the part to the left of. Since force is positive to the right, the stress is positive in conditions of tension. For a section of the rod, 1 < < 2, the total (rightward) force acting on that section due to the remainder of the rod is given by ( σ(2, t) σ( 1, t) ) S. If the density of the rod at is given by ρ 0 > 0, the momentum of this section is just 2 ρ 0 u t (, t) S d. 1 If we let F (, t) denote any eternal applied force per unit of volume in the positive direction, then we obtain from Newton s second law that d 2 ρ 0 u t (, t) S d = ( σ( 2, t) σ( 1, t) ) 2 S + F (, t) S d dt 1 for any such 1 < 2. For a sufficiently smooth displacement u(, t), we obtain the conservation of momentum equation (16) ρ 0 u tt (, t) σ (, t) = F (, t), a < < b, t > 0. The stress σ(, t) is determined by the type of material of which the rod is composed and the amount by which the neighboring region is stretched or compressed, i.e., on the elongation or strain, ε(, t). In order to define this, first note that a section [, + h] of the rod is deformed by the 1

10 10 FLOW AND TRANSPORT displacement to the new position [ + u(), ( + h) + u( + h)]. The elongation is the limiting increment of the change in the length due to the deformation as given by [u( + h) + ( + h)] [u() + ] h lim h 0 h = d u() d, so the strain is given by ε(, t) u (, t). The relation between the stress and strain is a constitutive law, usually determined by eperiment, and it depends on the type of material. In the simplest case, with small displacements, we find by eperiment that σ(, t) is proportional to ε(, t), i.e., that there is a constant k called Young s modulus for which (17) σ(, t) = k ε(, t). The constant k is a property of the material, and in this case we say the material is purely elastic. The partial differential equation for the longitudinal vibrations of the rod is obtained by substituting (17) into (16) to obtain (18) ρ 0 u tt ku = F (, t), a < < b, t > 0. This is the second-order wave equation. A rate-dependent component of the stress-strain relationship arises when the force generated by the elongation depends not only on the magnitude of the strain but also on the speed at which it is changed, i.e., on the strain rate ε t (, t) = v (, t). The simplest such case is that of a visco-elastic material defined by the linear constitutive equation σ(, t) = k ε(, t) + µ ε t (, t), in which the material constant µ is the viscosity or internal friction of the material. In terms of displacement, this has the form (19) σ(, t) = ku (, t) + µu t (, t). The partial differential equation for the longitudinal vibrations of the visco-elastic rod is obtained by substituting (19) into (16) to obtain (20) ρ 0 u tt µu t ku = F (, t), a < < b, t > 0.

11 4. LONGITUDINAL VIBRATIONS 11 For µ > 0 this is the viscous wave equation. In order to see the character of this equation, note that the highest order terms are of the form ρ 0 v t µv for the velocity v = u t, so this equation can be epected to be similar to the diffusion equation. For this reason, (20) is also called the strongly-damped wave equation.

Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

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

7.4 The Elementary Beam Theory In this section, problems involving long and slender beams are addressed. s with pressure vessels, the geometry of the beam, and the specific type of loading which will be

skiladæmi 10 Due: 11:59pm on Wednesday, November 11, 015 You will receive no credit for items you complete after the assignment is due Grading Policy Alternative Exercise 1115 A bar with cross sectional

Applications of PDEs 25.2 Introduction In this Section we discuss briefly some of the most important PDEs that arise in various branches of science and engineering. We shall see that some equations can

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

The acoustic characterization of porous media and its standards Luc JAOUEN 1, François-Xavier BECOT, Fabien CHEVILLOTTE Matelys, France ABSTRACT While there is a growing number of methods for the acoustic

St Olave s Grammar School Mock Practical skills.. a Can you design experiments, including ones to solve problems set in a practical context?.. b Can you identify the variables that must be controlled in

1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds

3. BEAMS: STRAIN, STRESS, DEFLECTIONS The beam, or flexural member, is frequently encountered in structures and machines, and its elementary stress analysis constitutes one of the more interesting facets

Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform

Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects

Lecture 4 Spring problem and conservation of mechanical energy Hooke's Law The restoring force exerted by the spring is directly proportional to its displacement. The restoring force acts in a direction

UNIVERSITY OF ALABAMA Department of Physics and Astronomy PH 15 / LeClair Spring 009 Exam II: Solutions 1. A block of mass m is released from rest at a height d=40 cm and slides down a frictionless ramp

Flow through porous media 1 Flow through porous media Introduction This report shows the capabilities of Tdyn for modelling the fluid flow through porous media. Modelling approach Modeling of flows through

CHAPTER THE EFFECTS OF FORCES ON MATERIALS EXERCISE 1, Page 50 1. A rectangular bar having a cross-sectional area of 80 mm has a tensile force of 0 kn applied to it. Determine the stress in the bar. Stress

(ISO/IEC - 700-005 Certified) SUMMER 7 EXAMINATION 70 Model ject Code: Important Instructions to examiners: ) The answers should be examined by key words and not as word-to-word as given in the model answer

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

Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) 5.1 DEFINITION A construction member is subjected to centric (axial) tension or compression if in any cross section the single distinct stress

Chapter II: Reversible process and work 1- Process Defined by change in a system, a thermodynamic process is a passage of a thermodynamic system from an initial to a final state of thermodynamic equilibrium.

Physics Including Human Applications 280 Chapter 13 ELASTIC PROPERTIES OF MATERIALS GOALS When you have mastered the contents of this chapter, you will be able to achieve the following goals: Definitions

Mechanics and Thermodynamics GEIT5 Particle dynamics. Basic concepts The aim of dynamics is to find the cause of motion and knowing the cause to give description of motion. The cause is always some interaction.

Chapter 3 Mass quantity of matter composing a body represented by m Kinetic Concepts for Analyzing Human Motion units are kg Inertia tendency to resist change in state of motion proportional to mass has

Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial

100406 Quiz 1. Introduction to Polymers 1) Polymers are different than low-molecular weight oligomers. For example an oligomeric polyethylene is wax, oligomeric polystyrene is similar to naphthalene (moth

PHY 101 DR M. A. ELERUJA KINETIC ENERGY AND WORK POTENTIAL ENERGY AND CONSERVATION OF ENERGY CENTRE OF MASS AND LINEAR MOMENTUM Work is done by a force acting on an object when the point of application

An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.

urther pplications of Newton s Laws - riction Static and Kinetic riction The normal force is related to friction. When two surfaces slid over one another, they experience a force do to microscopic contact

Summary. Shear is the gradient of velocity in a direction normal to the velocity. In the presence of shear, collisions among molecules in random motion tend to transfer momentum down-shear (from faster

1. A particle of mass 0.5 kg is attached to one end of a light elastic spring of natural length 0.9 m and modulus of elasticity λ newtons. The other end of the spring is attached to a fixed point O 3 on

Thermal physics revision questions ONE SECTION OF QUESTIONS TO BE COMPLETED AND MARKED EVERY WEEK AFTER HALF TERM. Section 1: Energy 1. Define the law of conservation of energy. Energy is neither created

Chapter 2 CONTINUUM MECHANICS PROBLEMS The concept of treating solids and fluids as though they are continuous media, rather thancomposedofdiscretemolecules, is one that is widely used in most branches

Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass

1. Explicit algebraic models Two explicit algebraic models are here compared in order to assess their predictive capabilities in the simulation of boundary layer flow cases. The studied models are both

MECHANICS OF STRUCTURES SCI 1105 COURSE MATERIAL UNIT - I Engineering Mechanics Branch of science which deals with the behavior of a body with the state of rest or motion, subjected to the action of forces.

Lecture 8 Equilibrium and Elasticity July 19 EQUILIBRIUM AND ELASTICITY CHAPTER 12 Give a sharp blow one end of a stick on the table. Find center of percussion. Baseball bat center of percussion Equilibrium

KINEMATICS & DYNAMICS BY ADVANCED DIFFERENTIAL EQUATIONS Question (**+) In this question take g = 0 ms. A particle of mass M kg is released from rest from a height H m, and allowed to fall down through

Objectives In this lecture you will learn the following We shall summarise the principles used in fluid mechanics and heat transfer. It is assumed that the student has already been exposed to courses in

Determining the ormal Modes of Vibration Introduction at the end of last lecture you determined the symmetry and activity of the vibrational modes of ammonia Γ vib 3 ) = A 1 IR, pol) + EIR,depol) the vibrational

KOZENY-CARMAN EQUATION REVISITED Jack Dvorkin -- 009 Abstract The Kozeny-Carman equation is often presented as permeability versus porosity, grain size, and tortuosity. When it is used to estimate permeability

Chapter 10 TH GNRAL LASTICITY PROBLM IN SOLIDS In Chapters 3-5 and 8-9, we have developed equilibrium, kinematic and constitutive equations for a general three-dimensional elastic deformable solid bod.

9 trength Theories of Lamina 9- TRENGTH O ORTHOTROPIC LAMINA or isotropic materials the simplest method to predict failure is to compare the applied stresses to the strengths or some other allowable stresses.

Review before first test Physical Mechanics Fall 000 Newton's Laws You should be able to state these laws using both words and equations. The nd law most important for meteorology. Second law: net force