DOCSLIB.ORG

Viscous Fluids

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Viscous Fluids Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass den- sity, velocity and pressure. An Eulerian description of uids focuses on the transport of a small xed volume. Real uids have the ability to resist stress and are known as viscous uids. Motion of these uids is governed by the Euler equation, where the Navier-Stokes equation is a simplica- tion of the Euler equation due to imcompressibilty. These equations are derived from conservation laws of mass, momentum and energy. 1 1 Introduction Viscosity is a measure of the resistance of a uid to a stress and diers among materials that may have similar properties such as density and weight. For example, water ows very dierently than oil. In order to analyze this, we can think of placing a material between two plates, where the bottom plate is xed and the top plate is allowed to move. If we have a uid between the plates and we apply a pressure to the upper plate, the plate will move continuously with a certain velocity. The uid in contact with the upper plate moves with a specic velocity, but the uid in contact with the lower plate has a zero velocity due to the solid boundary condition on the lower plate. This tells us that the uid between the plates developes a velocity gradient. The viscosity is the uid property that relates the shearing stress, or pressure per area, to the velocity gradient. If the shearing stress is linearly related to the velocity gradient, then the uid is known as Newtonian, whereas if the shearing stress is not linearly related to the velocity gradient, the uid is Non-Newtonian. In order to under- stand how the velocity varies, we need to be able to examine an innitesimal volume within the ow, which we do through dierential analysis. The equa- tions of motion for uids are governed by conservation of mass, momentum and energy. For incompressible, Newtonian uids, the Navier-Stokes equation is the dierential equation of motion. What follows is the derivation of these equations and their use in a simple example example describing uid ow in a pipe. 2 Derivations 2.1 Conservation Laws In order to look at what happens to a particular part of a uid as it moves through a system we need to look at the changes of an extensive property of the uid over a control volume, which is done using the Reynolds Transport theorem. The theorem states that these changes are equal to what is lost (or gained) through the boundary of the volume and what is created (or consumed) by the sources or sinks in the control volume. This is shown in the following integral where L is an extensive property of the uid, Ω is the control volume, v is the velocity and Q is the sink or source in the uid. d L dV = − Lv · n dA − Q dV (1) dt ˆΩ ˆ@Ω ˆΩ The divergence theorem can then be applied to change the surface integral,@Ω, to a volume integral,Ω. d L dV = − r · (Lv) dV − Q dV (2) dt ˆΩ ˆΩ ˆΩ Applying Leibniz's rule for dierentiation under an integral we get the following. 2 @L dV = − r · (Lv) dV − Q dV (3) ˆΩ @t ˆΩ ˆΩ And then combining the integrals we get, @L + r · (Lv) + Q dV = 0 (4) ˆΩ @t where the integrand is equally zero. @L + r · (Lv) + Q = 0: (5) @t This relationship gives us a general continuity equation in which we can then examine the conservation of mass, momentum and energy. 2.1.1 Conservation of Mass The continuity equation that describes liquids arises from looking at the con- servation of mass. Assuming Q = 0 or no sources or sinks of mass and looking at the mass density, ρ is the extensive property, the general continuity equation becomes @ρ + r · (ρv) = 0 (6) @t For an incompressible uid, ρ is constant and this equation simplies to r · v = 0 (7) which is a statement of conservation of volume. 2.1.2 Conservation of Momentum Conservation of momentum leads to discussion of a stress tensor and euler's equations of motion for the uid. Applying the continuity equation (Eq. 5) to momentum, ρv, we get the following equation. @ (ρv) + r · (ρvv) + Q = 0 (8) @t Expanding the derivatives and taking b to be a body force as a sink or source of momentum, @ρ @v v + ρ + r(ρv) · v + ρvr · v = b (9) @t @t we have expanded the divergence of a second rank tensor which gives a rst rank tensor @ρ @v v + ρ + r(ρ)v · v + ρr(v) · v + vρr · v = b (10) @t @t 3 if we rearrange this equation to be of the form @ρ @v v + ρ + vv · rρ + ρv · rv + ρvr · v = b (11) @t @t and pulling out v and ρ, @ρ @v v + v · rρ + ρr · v + ρ + v · rv = b (12) @t @t we can then recognize v · rρ + ρr · v = r · (ρv). @ρ @v v + r · (ρv) + ρ + v · rv = b (13) @t @t The left most term in Eq. 13 is the conservation of mass we saw earlier that is equal to zero (Eq. 6,7). Therefore, we are left with @v Dv ρ + v · rv = b ) ρ = b (14) @t Dt where D is the convective derivative. This equation is just an expression for Dt Newton's second law where we have body forces instead of point forces. And Euler's equations of motion follow directly from this, where b can be broken up into two terms, one describing forces resulting from stresses (σ is stress tensor) and the other describing forces such as gravity (f). Dv ρ = r · σ + f (15) Dt 2.1.3 Conservation of Energy Conservation of energy will lead to information pertaining to the thermody- namics of the system, which I am choosing not to focus on in this paper. For those interested, the extensive property energy density which has two parts, the kinetic energy density, 1 2, associated with the center-of-mass motion, and the 2 ρv internal energy density, ρ, where is the internal energy per unit mass will be used as L in Eq. 5. 2.2 Viscous Stress Tensor If we have a uid in a steady shear ow, where vx increases linearly with y, v = vx(y)^x (16) we have incompressible ow because the divergence is zero. The drag force acting on uid lying below area dA is then @v F = η x dA (17) x @y 4 where η is the viscosity. Refer to Fig. 1 for a schematic on viscous drag force. Figure 1: Viscous drag force exerted by uid above surface y=y0on uid below for linear shear ow. Since viscosity contributes to the surface stress, it is added in to the stress tensor. Newton's second law for a xed volume, subject to an external force f is as follows 3 @ (ρvk) 3 X 3 d x = − TkldA + ρf d x (18) ˆ @t ˆ ˆ k V l=1 A V And in dierential form this is exactly what we got out of applying the continuity equation to momentum (Eq. 15). 3 @ (ρvk) X @Tkl = − + ρfk (19) @t @xl l=1 The symmetric stress tensor for nonviscous or ideal uids has the form ideal Tkl = ρvkvl + pδkl (20) Where δkl is the Kronecker delta and we just need to add the part that describes the viscous forces. Since the uid is stationary, i.e. isotropic with no preferred spatial direction, the only terms that remain in the stress tensor are 5 @v @v k + l (21) @xl @xk 3 X @vj δkl = δkl(r:v) (22) @x j=1 j We assume the viscous stress tensor is a linear combination of these two. v @vk @vl 2δklr:v Tkl = −η + − − ζδkl(r:v) (23) @xl @xk 3 Where η is the viscosity and ζ is the bulk viscosity. This simplies for incom- pressible uids where the divergence of the velocity is zero meaning the bulk viscosity no longer contributes. v @vk @vl Tkl = −η + (24) @xl @xk So for our example shown in Figure 1, the viscous stress tensor becomes the following. @v T v = T v = −η x (25) xy yx @y The total stress tensor is then the sum of the ideal and viscous stress tensors. ideal v Tij = Tij + Tij (26) @vi @vj 2δijr:v Tij = ρvivj + pδij − η + − − ζδij(r:v) (27) @xj @xi 3 Using the expression for Newton's second law (Eq. 18) for a particular stress tensor. 3 3 2 3 2 X @ (ρvivj) @ (ρvi) η X @ vj X @ vi @p + = ζ + + η + ρf − (28) @x @t 3 @x @x @x 2 i @x j=1 j j=1 i j j=1 j i The left side simplies with the continuity equation (Eq. 5) 3 X @ (ρvj) @ρ + = 0 (29) @x @t j=1 j to give the following result. 3 3 2 3 2 X @vi @vi 1 η X @ vj η X @ vi 1 @p v + = ζ + + + f − (30) j @x @t ρ 3 @x @x ρ @x 2 i ρ @x j=1 j j=1 i j j=1 j i 6 And in vector form we get the generlized Euler equation for a vicous uid. @v 1 η η 1 (v:r)v + = ζ + r(r:v) + r2v + f − rp (31) @t ρ 3 ρ ρ With the simplication for incompressible ow, we get the Navier-Stokes equa- tion.
Load more

Recommended publications
  • 10. Collisions • Use Conservation of Momentum and Energy and The
    10. Collisions • Use conservation of momentum and energy and the center of mass to understand collisions between two objects. • During a collision, two or more objects exert a force on one another for a short time: -F(t) F(t) Before During After • It is not necessary for the objects to touch during a collision, e.g. an asteroid flied by the earth is considered a collision because its path is changed due to the gravitational attraction of the earth. One can still use conservation of momentum and energy to analyze the collision. Impulse: During a collision, the objects exert a force on one another. This force may be complicated and change with time. However, from Newton's 3rd Law, the two objects must exert an equal and opposite force on one another. F(t) t ti tf Dt From Newton'sr 2nd Law: dp r = F (t) dt r r dp = F (t)dt r r r r tf p f - pi = Dp = ò F (t)dt ti The change in the momentum is defined as the impulse of the collision. • Impulse is a vector quantity. Impulse-Linear Momentum Theorem: In a collision, the impulse on an object is equal to the change in momentum: r r J = Dp Conservation of Linear Momentum: In a system of two or more particles that are colliding, the forces that these objects exert on one another are internal forces. These internal forces cannot change the momentum of the system. Only an external force can change the momentum. The linear momentum of a closed isolated system is conserved during a collision of objects within the system.
    [Show full text]
  • AST242 LECTURE NOTES PART 3 Contents 1. Viscous Flows 2 1.1. the Velocity Gradient Tensor 2 1.2. Viscous Stress Tensor 3 1.3. Na
    AST242 LECTURE NOTES PART 3 Contents 1. Viscous flows 2 1.1. The velocity gradient tensor 2 1.2. Viscous stress tensor 3 1.3. Navier Stokes equation { diffusion 5 1.4. Viscosity to order of magnitude 6 1.5. Example using the viscous stress tensor: The force on a moving plate 6 1.6. Example using the Navier Stokes equation { Poiseuille flow or Flow in a capillary 7 1.7. Viscous Energy Dissipation 8 2. The Accretion disk 10 2.1. Accretion to order of magnitude 14 2.2. Hydrostatic equilibrium for an accretion disk 15 2.3. Shakura and Sunyaev's α-disk 17 2.4. Reynolds Number 17 2.5. Circumstellar disk heated by stellar radiation 18 2.6. Viscous Energy Dissipation 20 2.7. Accretion Luminosity 21 3. Vorticity and Rotation 24 3.1. Helmholtz Equation 25 3.2. Rate of Change of a vector element that is moving with the fluid 26 3.3. Kelvin Circulation Theorem 28 3.4. Vortex lines and vortex tubes 30 3.5. Vortex stretching and angular momentum 32 3.6. Bernoulli's constant in a wake 33 3.7. Diffusion of vorticity 34 3.8. Potential Flow and d'Alembert's paradox 34 3.9. Burger's vortex 36 4. Rotating Flows 38 4.1. Coriolis Force 38 4.2. Rossby and Ekman numbers 39 4.3. Geostrophic flows and the Taylor Proudman theorem 39 4.4. Two dimensional flows on the surface of a planet 41 1 2 AST242 LECTURE NOTES PART 3 4.5. Thermal winds? 42 5.
    [Show full text]
  • Derivation of Fluid Flow Equations
    TPG4150 Reservoir Recovery Techniques 2017 1 Fluid Flow Equations DERIVATION OF FLUID FLOW EQUATIONS Review of basic steps Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations, and constitutive equations for the fluids and the porous material involved. For simplicity, we will in the following assume isothermal conditions, so that we not have to involve an energy conservation equation. However, in cases of changing reservoir temperature, such as in the case of cold water injection into a warmer reservoir, this may be of importance. Below, equations are initially described for single phase flow in linear, one- dimensional, horizontal systems, but are later on extended to multi-phase flow in two and three dimensions, and to other coordinate systems. Conservation of mass Consider the following one dimensional rod of porous material: Mass conservation may be formulated across a control element of the slab, with one fluid of density ρ is flowing through it at a velocity u: u ρ Δx The mass balance for the control element is then written as: ⎧Mass into the⎫ ⎧Mass out of the ⎫ ⎧ Rate of change of mass⎫ ⎨ ⎬ − ⎨ ⎬ = ⎨ ⎬ , ⎩element at x ⎭ ⎩element at x + Δx⎭ ⎩ inside the element ⎭ or ∂ {uρA} − {uρA} = {φAΔxρ}. x x+ Δx ∂t Dividing by Δx, and taking the limit as Δx approaches zero, we get the conservation of mass, or continuity equation: ∂ ∂ − (Aρu) = (Aφρ). ∂x ∂t For constant cross sectional area, the continuity equation simplifies to: ∂ ∂ − (ρu) = (φρ) . ∂x ∂t Next, we need to replace the velocity term by an equation relating it to pressure gradient and fluid and rock properties, and the density and porosity terms by appropriate pressure dependent functions.
    [Show full text]
  • Law of Conversation of Energy
    Law of Conservation of Mass: "In any kind of physical or chemical process, mass is neither created nor destroyed - the mass before the process equals the mass after the process." - the total mass of the system does not change, the total mass of the products of a chemical reaction is always the same as the total mass of the original materials. "Physics for scientists and engineers," 4th edition, Vol.1, Raymond A. Serway, Saunders College Publishing, 1996. Ex. 1) When wood burns, mass seems to disappear because some of the products of reaction are gases; if the mass of the original wood is added to the mass of the oxygen that combined with it and if the mass of the resulting ash is added to the mass o the gaseous products, the two sums will turn out exactly equal. 2) Iron increases in weight on rusting because it combines with gases from the air, and the increase in weight is exactly equal to the weight of gas consumed. Out of thousands of reactions that have been tested with accurate chemical balances, no deviation from the law has ever been found. Law of Conversation of Energy: The total energy of a closed system is constant. Matter is neither created nor destroyed – total mass of reactants equals total mass of products You can calculate the change of temp by simply understanding that energy and the mass is conserved - it means that we added the two heat quantities together we can calculate the change of temperature by using the law or measure change of temp and show the conservation of energy E1 + E2 = E3 -> E(universe) = E(System) + E(Surroundings) M1 + M2 = M3 Is T1 + T2 = unknown (No, no law of conservation of temperature, so we have to use the concept of conservation of energy) Total amount of thermal energy in beaker of water in absolute terms as opposed to differential terms (reference point is 0 degrees Kelvin) Knowns: M1, M2, T1, T2 (Kelvin) When add the two together, want to know what T3 and M3 are going to be.
    [Show full text]
  • Continuity Equation in Pressure Coordinates
    Continuity Equation in Pressure Coordinates Here we will derive the continuity equation from the principle that mass is conserved for a parcel following the fluid motion (i.e., there is no flow across the boundaries of the parcel). This implies that δxδyδp δM = ρ δV = ρ δxδyδz = − g is conserved following the fluid motion: 1 d(δM ) = 0 δM dt 1 d()δM = 0 δM dt g d ⎛ δxδyδp ⎞ ⎜ ⎟ = 0 δxδyδp dt ⎝ g ⎠ 1 ⎛ d(δp) d(δy) d(δx)⎞ ⎜δxδy +δxδp +δyδp ⎟ = 0 δxδyδp ⎝ dt dt dt ⎠ 1 ⎛ dp ⎞ 1 ⎛ dy ⎞ 1 ⎛ dx ⎞ δ ⎜ ⎟ + δ ⎜ ⎟ + δ ⎜ ⎟ = 0 δp ⎝ dt ⎠ δy ⎝ dt ⎠ δx ⎝ dt ⎠ Taking the limit as δx, δy, δp → 0, ∂u ∂v ∂ω Continuity equation + + = 0 in pressure ∂x ∂y ∂p coordinates 1 Determining Vertical Velocities • Typical large-scale vertical motions in the atmosphere are of the order of 0. 01-01m/s0.1 m/s. • Such motions are very difficult, if not impossible, to measure directly. Typical observational errors for wind measurements are ~1 m/s. • Quantitative estimates of vertical velocity must be inferred from quantities that can be directly measured with sufficient accuracy. Vertical Velocity in P-Coordinates The equivalent of the vertical velocity in p-coordinates is: dp ∂p r ∂p ω = = +V ⋅∇p + w dt ∂t ∂z Based on a scaling of the three terms on the r.h.s., the last term is at least an order of magnitude larger than the other two. Making the hydrostatic approximation yields ∂p ω ≈ w = −ρgw ∂z Typical large-scale values: for w, 0.01 m/s = 1 cm/s for ω, 0.1 Pa/s = 1 μbar/s 2 The Kinematic Method By integrating the continuity equation in (x,y,p) coordinates, ω can be obtained from the mean divergence in a layer: ⎛ ∂u ∂v ⎞ ∂ω ⎜ + ⎟ + = 0 continuity equation in (x,y,p) coordinates ⎝ ∂x ∂y ⎠ p ∂p p2 p2 ⎛ ∂u ∂v ⎞ ∂ω = − ⎜ + ⎟ ∂p rearrange and integrate over the layer ∫p ∫ ⎜ ⎟ 1 ∂x ∂y p1⎝ ⎠ p ⎛ ∂u ∂v ⎞ ω(p )−ω(p ) = (p − p )⎜ + ⎟ overbar denotes pressure- 2 1 1 2 ⎜ ⎟ weighted vertical average ⎝ ∂x ∂y ⎠ p To determine vertical motion at a pressure level p2, assume that p1 = surface pressure and there is no vertical motion at the surface.
    [Show full text]
  • CONTINUUM MECHANICS for ENGINEERS Second Edition Second Edition
    CONTINUUM MECHANICS for ENGINEERS Second Edition Second Edition CONTINUUM MECHANICS for ENGINEERS G. Thomas Mase George E. Mase CRC Press Boca Raton London New York Washington, D.C. Library of Congress Cataloging-in-Publication Data Mase, George Thomas. Continuum mechanics for engineers / G. T. Mase and G. E. Mase. -- 2nd ed. p. cm. Includes bibliographical references (p. )and index. ISBN 0-8493-1855-6 (alk. paper) 1. Continuum mechanics. I. Mase, George E. QA808.2.M364 1999 531—dc21 99-14604 CIP This book contains information obtained from authentic and highly regarded sources. Reprinted material is quoted with permission, and sources are indicated. A wide variety of references are listed. Reasonable efforts have been made to publish reliable data and informa- tion, but the author and the publisher cannot assume responsibility for the validity of all materials or for the consequences of their use. Neither this book nor any part may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, microfilming, and recording, or by any information storage or retrieval system, without prior permission in writing from the publisher. The consent of CRC Press LLC does not extend to copying for general distribution, for promotion, for creating new works, or for resale. Specific permission must be obtained in writing from CRC Press LLC for such copying. Direct all inquiries to CRC Press LLC, 2000 N.W. Corporate Blvd., Boca Raton, Florida 33431. Trademark Notice: Product or corporate names may be trademarks or registered trademarks, and are only used for identification and explanation, without intent to infringe.
    [Show full text]
  • Heat and Energy Conservation
    1 Lecture notes in Fluid Dynamics (1.63J/2.01J) by Chiang C. Mei, MIT, Spring, 2007 CHAPTER 4. THERMAL EFFECTS IN FLUIDS 4-1-2energy.tex 4.1 Heat and energy conservation Recall the basic equations for a compressible fluid. Mass conservation requires that : ρt + ∇ · ρ~q = 0 (4.1.1) Momentum conservation requires that : = ρ (~qt + ~q∇ · ~q)= −∇p + ∇· τ +ρf~ (4.1.2) = where the viscous stress tensor τ has the components = ∂qi ∂qi ∂qk τ = τij = µ + + λ δij ij ∂xj ∂xi ! ∂xk There are 5 unknowns ρ, p, qi but only 4 equations. One more equation is needed. 4.1.1 Conservation of total energy Consider both mechanical ad thermal energy. Let e be the internal (thermal) energy per unit mass due to microscopic motion, and q2/2 be the kinetic energy per unit mass due to macroscopic motion. Conservation of energy requires D q2 ρ e + dV rate of incr. of energy in V (t) Dt ZZZV 2 ! = − Q~ · ~ndS rate of heat flux into V ZZS + ρf~ · ~qdV rate of work by body force ZZZV + Σ~ · ~qdS rate of work by surface force ZZX Use the kinematic transport theorm, the left hand side becomes D q2 ρ e + dV ZZZV Dt 2 ! 2 Using Gauss theorem the heat flux term becomes ∂Qi − QinidS = − dV ZZS ZZZV ∂xi The work done by surface stress becomes Σjqj dS = (σjini)qj dS ZZS ZZS ∂(σijqj) = (σijqj)ni dS = dV ZZS ZZZV ∂xi Now all terms are expressed as volume integrals over an arbitrary material volume, the following must be true at every point in space, 2 D q ∂Qi ∂(σijqi) ρ e + = − + ρfiqi + (4.1.3) Dt 2 ! ∂xi ∂xj As an alternative form, we differentiate the kinetic energy and get De
    [Show full text]
  • Chapter 3 Newtonian Fluids
    CM4650 Chapter 3 Newtonian Fluid 2/5/2018 Mechanics Chapter 3: Newtonian Fluids CM4650 Polymer Rheology Michigan Tech Navier-Stokes Equation v vv p 2 v g t 1 © Faith A. Morrison, Michigan Tech U. Chapter 3: Newtonian Fluid Mechanics TWO GOALS •Derive governing equations (mass and momentum balances •Solve governing equations for velocity and stress fields QUICK START V W x First, before we get deep into 2 v (x ) H derivation, let’s do a Navier-Stokes 1 2 x1 problem to get you started in the x3 mechanics of this type of problem solving. 2 © Faith A. Morrison, Michigan Tech U. 1 CM4650 Chapter 3 Newtonian Fluid 2/5/2018 Mechanics EXAMPLE: Drag flow between infinite parallel plates •Newtonian •steady state •incompressible fluid •very wide, long V •uniform pressure W x2 v1(x2) H x1 x3 3 EXAMPLE: Poiseuille flow between infinite parallel plates •Newtonian •steady state •Incompressible fluid •infinitely wide, long W x2 2H x1 x3 v (x ) x1=0 1 2 x1=L p=Po p=PL 4 2 CM4650 Chapter 3 Newtonian Fluid 2/5/2018 Mechanics Engineering Quantities of In more complex flows, we can use Interest general expressions that work in all cases. (any flow) volumetric ⋅ flow rate ∬ ⋅ | average 〈 〉 velocity ∬ Using the general formulas will Here, is the outwardly pointing unit normal help prevent errors. of ; it points in the direction “through” 5 © Faith A. Morrison, Michigan Tech U. The stress tensor was Total stress tensor, Π: invented to make the calculation of fluid stress easier. Π ≡ b (any flow, small surface) dS nˆ Force on the S ⋅ Π surface V (using the stress convention of Understanding Rheology) Here, is the outwardly pointing unit normal of ; it points in the direction “through” 6 © Faith A.
    [Show full text]
  • Lecture 5: Magnetic Mirroring
    !"#$%&'()%"*#%*+,-./-*+01.2(.*3+456789* !"#$%&"'()'*+,-".#'*/&&0&/-,' Dr. Peter T. Gallagher Astrophysics Research Group Trinity College Dublin :&2-;-)(*!"<-$2-"(=*%>*/-?"=)(*/%/="#* o Gyrating particle constitutes an electric current loop with a dipole moment: 1/2mv2 µ = " B o The dipole moment is conserved, i.e., is invariant. Called the first adiabatic invariant. ! o µ = constant even if B varies spatially or temporally. If B varies, then vperp varies to keep µ = constant => v|| also changes. o Gives rise to magnetic mirroring. Seen in planetary magnetospheres, magnetic bottles, coronal loops, etc. Bz o Right is geometry of mirror Br from Chen, Page 30. @-?"=)(*/2$$%$2"?* o Consider B-field pointed primarily in z-direction and whose magnitude varies in z- direction. If field is axisymmetric, B! = 0 and d/d! = 0. o This has cylindrical symmetry, so write B = Brrˆ + Bzzˆ o How does this configuration give rise to a force that can trap a charged particle? ! o Can obtain Br from " #B = 0 . In cylindrical polar coordinates: 1 " "Bz (rBr )+ = 0 r "r "z ! " "Bz => (rBr ) = #r "r "z o If " B z / " z is given at r = 0 and does not vary much with r, then r $Bz 1 2&$ Bz ) ! rBr = " #0 r dr % " r ( + $z 2 ' $z *r =0 ! 1 &$ Bz ) Br = " r( + (5.1) 2 ' $z *r =0 ! @-?"=)(*/2$$%$2"?* o Now have Br in terms of BZ, which we can use to find Lorentz force on particle. o The components of Lorentz force are: Fr = q(v" Bz # vzB" ) (1) F" = q(#vrBz + vzBr ) (2) (3) Fz = q(vrB" # v" Br ) (4) o As B! = 0, two terms vanish and terms (1) and (2) give rise to Larmor gyration.
    [Show full text]
  • Chapter 15 - Fluid Mechanics Thursday, March 24Th
    Chapter 15 - Fluid Mechanics Thursday, March 24th •Fluids – Static properties • Density and pressure • Hydrostatic equilibrium • Archimedes principle and buoyancy •Fluid Motion • The continuity equation • Bernoulli’s effect •Demonstration, iClicker and example problems Reading: pages 243 to 255 in text book (Chapter 15) Definitions: Density Pressure, ρ , is defined as force per unit area: Mass M ρ = = [Units – kg.m-3] Volume V Definition of mass – 1 kg is the mass of 1 liter (10-3 m3) of pure water. Therefore, density of water given by: Mass 1 kg 3 −3 ρH O = = 3 3 = 10 kg ⋅m 2 Volume 10− m Definitions: Pressure (p ) Pressure, p, is defined as force per unit area: Force F p = = [Units – N.m-2, or Pascal (Pa)] Area A Atmospheric pressure (1 atm.) is equal to 101325 N.m-2. 1 pound per square inch (1 psi) is equal to: 1 psi = 6944 Pa = 0.068 atm 1atm = 14.7 psi Definitions: Pressure (p ) Pressure, p, is defined as force per unit area: Force F p = = [Units – N.m-2, or Pascal (Pa)] Area A Pressure in Fluids Pressure, " p, is defined as force per unit area: # Force F p = = [Units – N.m-2, or Pascal (Pa)] " A8" rea A + $ In the presence of gravity, pressure in a static+ 8" fluid increases with depth. " – This allows an upward pressure force " to balance the downward gravitational force. + " $ – This condition is hydrostatic equilibrium. – Incompressible fluids like liquids have constant density; for them, pressure as a function of depth h is p p gh = 0+ρ p0 = pressure at surface " + Pressure in Fluids Pressure, p, is defined as force per unit area: Force F p = = [Units – N.m-2, or Pascal (Pa)] Area A In the presence of gravity, pressure in a static fluid increases with depth.
    [Show full text]
  • Fluid Dynamics and Mantle Convection
    Subduction II Fundamentals of Mantle Dynamics Thorsten W Becker University of Southern California Short course at Universita di Roma TRE April 18 – 20, 2011 Rheology Elasticity vs. viscous deformation η = O (1021) Pa s = viscosity µ= O (1011) Pa = shear modulus = rigidity τ = η / µ = O(1010) sec = O(103) years = Maxwell time Elastic deformation In general: σ ε ij = Cijkl kl (in 3-D 81 degrees of freedom, in general 21 independent) For isotropic body this reduces to Hooke’s law: σ λε δ µε ij = kk ij + 2 ij λ µ ε ε ε ε with and Lame’s parameters, kk = 11 + 22 + 33 Taking shear components ( i ≠ j ) gives definition of rigidity: σ µε 12 = 2 12 Adding the normal components ( i=j ) for all i=1,2,3 gives: σ λ µ ε κε kk = (3 + 2 ) kk = 3 kk with κ = λ + 2µ/3 = bulk modulus Linear viscous deformation (1) Total stress field = static + dynamic part: σ = − δ +τ ij p ij ij Analogous to elasticity … General case: τ = ε ij C'ijkl kl Isotropic case: τ = λ ε δ + ηε ij ' kk ij 2 ij Linear viscous deformation (2) Split in isotropic and deviatoric part (latter causes deformation): 1 σ ' = σ − σ δ = σ + pδ ij ij 3 kk ij ij ij 1 ε' = ε − ε δ ij ij 3 kk ij which gives the following stress: σ = − δ + ςε δ + ηε 'ij ( p p) ij kk ij 2 'ij ε = With compressibility term assumed 0 (Stokes condition kk 0 ) 2 ( ς = λ ' + η = bulk viscosity) 3 τ Now τ = 2 η ε or η = ij ij ij ε 2 ij η = In general f (T,d, p, H2O) Non-linear (or non-Newtonian) deformation General stress-strain relation: ε = Aτ n n = 1 : Newtonian n > 1 : non-Newtonian n → ∞ : pseudo-brittle 1 −1 1−n 1 −1/ n (1−n) / n Effective viscosity η = A τ = A ε eff 2 2 Application: different viscosities under oceans with different absolute plate motion, anisotropic viscosities by means of superposition (Schmeling, 1987) Microphysical theory and observations Maximum strength of materials (1) ‘Strength’ is maximum stress that material can resist In principle, viscous fluid has zero strength.
    [Show full text]
  • C:\Book\Booktex\Notind.Dvi 0
    362 INDEX A Cauchy stress law 216 Absolute differentiation 120 Cauchy-Riemann equations 293,321 Absolute scalar field 43 Charge density 323 Absolute tensor 45,46,47,48 Christoffel symbols 108,110,111 Acceleration 121, 190, 192 Circulation 293 Action integral 198 Codazzi equations 139 Addition of systems 6, 51 Coefficient of viscosity 285 Addition of tensors 6, 51 Cofactors 25, 26, 32 Adherence boundary condition 294 Compatibility equations 259, 260, 262 Aelotropic material 245 Completely skew symmetric system 31 Affine transformation 86, 107 Compound pendulum 195,209 Airy stress function 264 Compressible material 231 Almansi strain tensor 229 Conic sections 151 Alternating tensor 6,7 Conical coordinates 74 Ampere’s law 176,301,337,341 Conjugate dyad 49 Angle between vectors 80, 82 Conjugate metric tensor 36, 77 Angular momentum 218, 287 Conservation of angular momentum 218, 295 Angular velocity 86,87,201,203 Conservation of energy 295 Arc length 60, 67, 133 Conservation of linear momentum 217, 295 Associated tensors 79 Conservation of mass 233, 295 Auxiliary Magnetic field 338 Conservative system 191, 298 Axis of symmetry 247 Conservative electric field 323 B Constitutive equations 242, 251,281, 287 Continuity equation 106,234, 287, 335 Basic equations elasticity 236, 253, 270 Contraction 6, 52 Basic equations for a continuum 236 Contravariant components 36, 44 Basic equations of fluids 281, 287 Contravariant tensor 45 Basis vectors 1,2,37,48 Coordinate curves 37, 67 Beltrami 262 Coordinate surfaces 37, 67 Bernoulli’s Theorem 292 Coordinate transformations
    [Show full text]