DOCSLIB.ORG

Continuum Mechanics Lecture 4 Fluid Dynamics

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Continuum Mechanics Lecture 4 Fluid Dynamics Continuum Mechanics Lecture 4 Fluid dynamics Prof. Romain Teyssier http://www.itp.uzh.ch/~teyssier Continuum Mechanics 30/04/2013 Romain Teyssier Outline - Fluid kinematics - Mass and momentum conservation laws - The energy equation - Real fluids - Ideal fluids - Incompressible fluids Continuum Mechanics 30/04/2013 Romain Teyssier What is a fluid ? In thermodynamics or kinetic theory, a fluid is a collection of atoms or molecules, in liquid or gaseous form. In continuum mechanics, a fluid is a system that flows. The central property is the fluid velocity. In solid mechanics, we have studied various equilibrium solutions, for which the stress was related to the strain (static deformation): the elastic regime. Above a given threshold (the yield strength), the solid enters the plastic regime and eventually breaks. For a fluid, the deformation can be arbitrarily large. The stress is related to the rate of strain, not the strain. Note that the boundary between solid and fluid is fuzzy: waves, tooth paste, glaciers... It depends on the exact time scale and length scale we are interested in. Example: on short time scale, a glacier is a solid, but on longer time scale (months), it behaves like a fluid. A fluid will be described as a continuum if length scales are much larger than the atomic or molecular mean free path. Continuum Mechanics 30/04/2013 Romain Teyssier Fluid kinematics For each particle in the system, the velocity is defined as x(t +dt) x(t) v(t)= lim − dt 0 dt where the trajectory of the particle→ is x ( t ) . We consider the mapping between the initial positions y = x(t = 0) and the final position, defined as x ( t )= −→φ ( y, t ). The velocity field is the vector field defined by dx 1 v ( x, t )= ( y, t ) where y = −→φ − (x, t) dt We can now obtain the trajectory directory from the velocity field as t x(t)=x(0) + v(x(u),u)du 0 Streamlines are also used to visualize the flow. dx dy dz They are defined as = = vx(x, t) vy(x, t) vz(x, t) Only for a stationary flow are streamlines and trajectories equal. Continuum Mechanics 30/04/2013 Romain Teyssier Kinematics of a volume element Consider a scalar field Ψ ( x, t ) . Along a trajectory, we have Ψ x 0 ( t )= Ψ ( x ( t ) ,t ) . dΨ ∂Ψ DΨ Using the chain rule, we have x 0 ( t )= + v −→ Ψ = dt ∂t · ∇ Dt This is the Lagrangian derivative of the scalar function Ψ. We consider the scalar «color» function Ψ ( x )=1 for x V and Ψ(x)=0 x/V ∈ t ∈ t Using this color function, we can show that for any scalar field α we have: d ∂α d3xα(x, t)= d3x + −→ (αv) dt ∂t ∇ · Vt Vt This is Reynold’s transport theorem. d d ∂Ψ ∂α Proof: d3xα(x, t)= d3xΨ(x, t)α(x, t)= d3x α + Ψ dt dt 3 3 ∂t ∂t Vt R R DΨ ∂Ψ d ∂α We have = + v −→ Ψ =0 so d3xα(x, t)= d3x Ψ αv −→ Ψ Dt ∂t · ∇ dt 3 ∂t − · ∇ d Vt R ∂α Integrating by parts, we have d3xα(x, t)= d3x Ψ + Ψ−→ (αv) dt 3 ∂t ∇ · Vt R Continuum Mechanics 30/04/2013 Romain Teyssier Mass conservation The total mass in the volume element is obtained using α(x, t)=ρ(x, t) Since the total mass is conserved within the Lagrangian volume Vt, we have dM d ∂ρ = d3xρ(x, t)= d3x + −→ (ρv) =0 dt dt ∂t ∇ · Vt Vt Since this equation is true for any Lagrangian volume, we have proven the continuity equation (or mass conservation in continuous form) ∂ρ + −→ (ρv)=0 ∂t ∇ · We now consider the mass variation in a fixed Eulerian volume V0 dM d ∂ρ(x, t) 0 = d3xρ(x, t)= d3x = d3x−→ (ρv) dt dt ∂t − ∇ · V0 V0 V0 Using the divergence theorem, we have finally the equation for mass conservation in integral form d d3xρ(x, t)= ρv ndS dt − · V0 S0 Continuum Mechanics 30/04/2013 Romain Teyssier Lagrangian derivative We now compute the variation of a scalar field per unit mass β(x, t) d ∂β ∂ρ α = ρβ dx3ρβ = dx3 ρ + β + −→ (ρβv) −→ dt ∂t ∂t ∇ · Vt Vt Differentiating the product, we have −→ (ρβv)=β−→ (ρv)+ρv −→ β ∇ · ∇ · · ∇ and using the continuity equation, we have d ∂β Dβ dx3ρβ = dx3ρ + v −→ β = dx3ρ dt ∂t · ∇ Dt Vt Vt Vt We have used the Lagrangian derivative of the scalar field defined by Dβ ∂β = + v −→ β Dt ∂t · ∇ Scalar field per unit mass are also called specific quantities. For example, the internal energy per unit volume e and the internal energy per unit mass (or the specific energy) ε are related by e = ρ The fluid velocity is equivalent to a specific momentum. Continuum Mechanics 30/04/2013 Romain Teyssier Variation of the volume 3 We now used the scalar field α ( x, t )=1 . The total volume is trivially Vt = dx V dV t From Reynold’s transport theorem, we have: t = dx3−→ v dt ∇ · Vt 3 We define the specific volume V =1 / ρ so that Vt = dx ρV Vt DV 1 DV Using the Lagrangian derivative, we have ρ = = −→ v Dt V Dt ∇ · The velocity divergence is equal to the rate of change of the specific volume. Equivalently, the Lagrangian derivative of the density writes 1 Dρ = −→ v ρ Dt −∇ · 1 Dρ 1 DV An easy check: we have ρ V =1 so we obtain = = −→ v ρ Dt −V Dt −∇ · Continuum Mechanics 30/04/2013 Romain Teyssier Momentum conservation We use the dynamical equilibrium equation we derived in the previous lectures d ρvdV = ρ−→F dV + σndS dt Vt Vt St d Dv From the definition of the Lagrangian derivative, we get ρv dV = ρ i dV dt i Dt Vt Vt Using the divergence theorem, we get σndS = −→ σdV ∇ · St Vt So finally, we get the Euler equations in Lagrangian form Dv ρ = ρ−→F + −→ σ Dt ∇ · ∂v It writes in Eulerian form: ρ i + v −→ v = ρF + −→ σ ∂t · ∇ i i ∇ · −→ i Using the continuity equation, we derive the conservative form for the momentum ∂ρv i + −→ (ρv v σ )=ρF ∂t ∇ · i − −→ i i Momentum conservation in integral form in a fixed Eulerian volume: d dx3ρv + ρv v ndS σ ndS = dx3ρF dt i i · − −→ i · i V0 S0 S0 V0 Continuum Mechanics 30/04/2013 Romain Teyssier The internal energy equation From the first principle of thermodynamics, we have derived in the previous lecture dE = δQ + Tr σδ dV V For the small displacement field δ u = v d t , we define the rate of strain tensor 1 ∂v ∂v ˙ = i + j ij 2 ∂x ∂x j i The Lagrangian variation of the total internal energy is dE D Dq = dx3ρ = dx3ρ + dx3Tr(σ ˙) dt Dt Dt Vt Vt Vt We derive the internal energy equation in continuous form D Dq ρ = ρ +Tr(σ ˙) Dt Dt In case there is no external heat source or sink (adiabatic system), we have D ρ =Tr(σ ˙) Dt Warning: do not get confused between the specific energy and the strain tensor ! Continuum Mechanics 30/04/2013 Romain Teyssier The total energy equation dE dK δW δQ From the kinetic energy theorem, we have + = ext + dt dt dt dt δW The work per unit time for external forces is ext = dx3ρ−→F v + −→T vdS dt · · Vt St Expressing the work of the stress field as a volume integral (see previous lectures) −→T vdS = σn vdS = σv ndS = −→ σv dV · · · ∇ · St St St Vt we deduce the Lagrangian derivative of the specific total energy D v2 ρ + = −→ σv + ρv −→F Dt 2 ∇ · · v2 Define the total energy per unit volume as E = ρ + 2 We have the total energy conservation law in continuous form ∂E + −→ Ev σv = ρv −→F ∂t ∇ · − · Exercise: derive the total energy conservation in integral form in a fixed Eulerian volume. Continuum Mechanics 30/04/2013 Romain Teyssier Fluid dynamics equations Equations in Eulerian (conservative) form Equations in Lagrangian form ∂ρ 1 Dρ + −→ (ρv)=0 = −→ v ∂t ∇ · ρ Dt −∇ · ∂ρvi Dv + −→ (ρviv −→σi )=ρFi ρ = ρ−→F + −→ σ ∂t ∇ · − Dt ∇ · ∂E D + −→ Ev σv = ρv −→F ρ =Tr(σ ˙) ∂t ∇ · − · Dt To close the previous systems, we need the material law of the fluid σ = fonction (ρ, , ∂vi/∂xj,...) In general, we define the pressure and the viscous tensor as σ = p1+τ with in general Tr(τ)=0 − The pressure is given by the Equation-of-State (EoS) and depends on 2 thermodynamical quantities p = p ( ρ , ) or p = p(ρ,T) Continuum Mechanics 30/04/2013 Romain Teyssier Equation-of-State for Aluminum The SESAME library (LANL 1992) Continuum Mechanics 30/04/2013 Romain Teyssier Mie-Grüneisen EoS for real gases and fluids A general form for real fluid EoS: P P (ρ)=Γρ ( (ρ)) − c − c The main constituents are the cold pressure P c ( ρ ) and the cold energy c ( ρ ) . Both cold curves should satisfy the thermodynamical consistency condition: Pc(ρ) (ρ)= c ρ2 P Pc(ρ) The Mie-Grüneisen sound speed is cs = Pc(ρ)+(Γ + 1) − ρ Simple examples: 1- Ideal gas EoS: P =(γ 1)ρ − ρ γ 2- Polytropic EoS: P = P 0 ρ 0 2 3- Isothermal EoS: P = ρc0 Continuum Mechanics 30/04/2013 Romain Teyssier Real fluids and gases: linear viscosity Following the methodology of linear, isotropic and thermoelastic material, we consider a model where the stress tensor depends linearly on the rate of strain. ∂vi ∂vj τij = λ −→ v δij + µ + ∇ · ∂xj ∂xi where λ and µ are the equivalent of the Lamé coefficients (not the same units !).
Load more

Recommended publications
  • 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]
  • 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]
  • 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]
  • CONTINUITY EQUATION  Another Principle on Which We Can Derive a New Equation Is the Conservation of Mass
    ESCI 342 – Atmospheric Dynamics I Lesson 7 – The Continuity and Additional Equations Suggested Reading: Martin, Chapter 3 THE SYSTEM OF EQUATIONS IS INCOMPLETE The momentum equations in component form comprise a system of three equations with 4 unknown quantities (u, v, p, and ). Du1 p fv (1) Dt x Dv1 p fu (2) Dt y p g (3) z They are not a closed set, because there are four dependent variables (u, v, p, and ), but only three equations. We need to come up with some more equations in order to close the set. DERIVATION OF THE CONTINUITY EQUATION Another principle on which we can derive a new equation is the conservation of mass. The equation derived from this principle is called the mass continuity equation, or simply the continuity equation. Imagine a cube at a fixed point in space. The net change in mass contained within the cube is found by adding up the mass fluxes entering and leaving through each face of the cube.1 The mass flux across a face of the cube normal to the x-axis is given by u. Referring to the picture below, these fluxes will lead to a rate of change in mass within the cube given by m u yz u yz (4) t x xx 1 A flux is a quantity per unit area per unit time. Mass flux is therefore the rate at which mass moves across a unit area, and would have units of kg s1 m2. The mass in the cube can be written in terms of the density as m xyz so that m x y z .
    [Show full text]
  • 1 the Continuity Equation 2 the Heat Equation
    1 The Continuity Equation Imagine a fluid flowing in a region R of the plane in a time dependent fashion. At each point (x; y) R2 it has a velocity v = v (x; y; t) at time t. Let ρ = ρ(x; y; t) be the density of the fluid 2 −! −! at (x; y) at time t. Let P be any point in the interior of R and let Dr be the closed disk of radius r > 0 and center P . The mass of fluid inside Dr at any time t is ρ dxdy: ZZDr If matter is neither created nor destroyed inside Dr, the rate of decrease of this quantity is equal to the flux of the vector field −!J = ρ−!v across Cr, the positively oriented boundary of Dr. We therefore have d ρ dxdy = −!J −!N ds; dt − · ZZDr ZCr where N~ is the outer normal and ds is the element of arc length. Notice that the minus sign is needed since positive flux at time t represents loss of total mass at that time. Also observe that the amount of fluid transported across a small piece ds of the boundary of Dr at time t is ρ−!v −!N ds. Differentiating under the integral sign on the left-hand side and using the flux form of Green’s· Theorem on the right-hand side, we get @ρ dxdy = −!J dxdy: @t − r · ZZDr ZZDr Gathering terms on the left-hand side, we get @ρ ( + −!J ) dxdy = 0: @t r · ZZDr If the integrand was not zero at P it would be different from zero on Dr for some sufficiently small r and hence the integral would not be zero which is not the case.
    [Show full text]
  • Chapter 2 Wave Propagation in Viscous Fluid
    Chapter 2 Wave Propagation in Viscous Fluid This chapter summarizes with the derivation of the mathematical form of the acoustic wave propagation in the fluid. Before we derive the final form of the wave propagation equation in viscous fluid, we first look at two conservation (mass and momentum) of equations and state equation in the fluid. Detailed derivations can be found in the literatures [4, 6, 27, 28]. We limit our discussion only on the lossy 1-dimensional plane wave. 2.1 One-dimensional Viscous Wave Equation 2.1.1 One-dimensional Continuity Equation (Conversation of Mass) Consider the flow of a compressible fluid through a duct of arbitrary cross section (area S) in one dimension. The control volume (CV), is the segment between x and x + ∆x . 5 We want to know the rate at which the inside mass changes. First, we made two assumptions: 1. The CV is fixed in space 2. The flow is one-dimensional, so the mass flow only depends on t and x. Mass flow in Mass flow out = ρ uS |x = ρ uS |x+∆x x x+∆x , where ρ, u are the average mass flow density and average mass flow speed, respectively. For ∆x →0 , ρ and u will become a true point function. The time rate of increase of mass inside the CV is equal to net mass flow into the CV through the CV surfaces or in mathematical terms, ∂ (Sρ ∆x)= ρ uS | −ρ uS | ( 1 ) ∂t x x+∆x Since S is a constant and ∆x is not a function of time, Eq.
    [Show full text]
  • ESCI 342 – Atmospheric Dynamics I Lesson 10 – Vertical Motion, Pressure Coordinates
    ESCI 342 – Atmospheric Dynamics I Lesson 10 – Vertical Motion, Pressure Coordinates Reading: Martin, Section 4.1 PRESSURE COORDINATES Pressure is often a convenient vertical coordinate to use in place of altitude. If the hydrostatic approximation is used, the relationship between pressure and altitude is given by the hydrostatic equation, p g . (1) z In height coordinates the vertical velocity is defined as w Dz Dt . In pressure coordinates the vertical velocity is defined as Dp , (2) Dt and is commonly called simply omega. The units of are Pa/s (often microbars per second, b/s, is also used). Since pressure decreases upward, a negative omega means rising motion, while a positive omega means subsiding motion. w and are related as follows: Dp p p p p u v w . (3) Dt t x y z On the synoptic scale we can assume hydrostatic vertical balance, so that (3) becomes Dp p p p u v gw gw , (4) Dt t x y since the local pressure tendencies and horizontal pressure advection terms are much, much smaller in magnitude than the vertical pressure advection terms. The total derivative in pressure coordinates is D uv .1 (5) Dt t x y p The conversion of a height derivative to a pressure derivative is accomplished using the chain rule as follows z . (6) p p z g z In pressure coordinates, the directions of the unit vectors ( iˆ , ˆj , and kˆ ) are the same as in height coordinates. The x and y axes are still horizontal, and not oriented along the constant pressure surface.2 The vertical axis is still vertical (perpendicular to x and y.) 1 If you take the total derivative of pressure you end up with the seeming absurdity that p p p uv .
    [Show full text]
  • Fluid Mechanics – Continuity Equation Bernoulli’S Equation
    PART II Continuity and Conservation of Mass • Some applications of fluid mechanics • Flow through a pipe: • . Conservation. of mass for steady state (no storage) says m in = m out . m m _1A1Vm1 = _ 2A2Vm2 • For incompressible fluids, density does not changes (_ 1 = _ 2) so A1Vm1 = A2Vm2 = Q Fluid Mechanics – Continuity Equation Bernoulli’s equation • The equation of continuity states that for an Y incompressible fluid flowing in a tube of varying 2 A V cross-sectional area (A), the mass flow rate is the 2 2 same everywhere in the tube: . m1= m2 _ 1A1V1 = _ 2A2V2 _1A1V1= _2 A2V2 • Generally, the density stays constant and then it's Y For incompressible flow simply the flow rate (Av) that is constant. 1 A1 V1 A1V1= A2V2 Assume steady flow, V parallel to streamlines & no viscosity 1 Bernoulli Equation – energy Bernoulli Equation – work • Consider energy terms for steady flow: • Consider work done on the system is Force x distance • We write terms for KE and PE at each point • We write terms for force in terms of Pressure and area Y Y 2 2 Wi = FiVi dt =PiViAi dt A2 V2 A2 V2 E = KE + PE i i i Note ViAi dt = mi/_i 1 2 E1 = 2 m&1V1 + gm&1 y1 W1 = P1m&1 / ρ1 2 Y 1 Y 1 E2 = 2 m&2V2 + gm&2 y2 1 W2 = −P2m&2 / ρ2 A1 V1 A1 V1 As the fluid moves, work is being done by the external Now we set up an energy balance on the system.
    [Show full text]
  • Finite Element Simulation of Heat Transfer in Ferrofluid
    28 Finite Element Simulation of Heat Transfer in Ferrofluid Tomasz Strek Poznan University of Technology, Institute of Applied Mechanics Poland 1. Introduction During the last decades, an extensive research work has been done on the fluids dynamics in the presence of magnetic field (magnetorheological fluids and ferrofluids). The effect of magnetic field on fluids is worth investigating due to its innumerable applications in wide spectrum of fields. The study of interaction of the magnetic field or the electromagnetic field with fluids have been documented e.g. among nuclear fusion, chemical engineering, medicine, high speed noiseless printing and transformer cooling. Ferrofluids are industrially prepared magnetic fluids which consist of stable colloidal suspensions of small single-domain ferromagnetic particles in suitable carrier liquids. Usually these fluids do not conduct electric current and exhibit a nonlinear paramagnetic behavior. An external magnetic field imposed on a ferrofluid with varying susceptibility, e.g., due to a temperature gradient, results in a non-uniform magnetic body force, which leads to a form of heat transfer called thermomagnetic convection. This form of heat transfer can be useful when conventional convection heat transfer is inadequate, e.g., in miniature microscale devices or under reduced gravity conditions. This chapter contains numerical results of computer simulation of heat transfer through a ferrofluid in channel flow under the influence of magnetic dipole, ferrofluid cooling of heat- generating device and heat transfer through a ferrofluid in channel with porous walls. All considered problems describe two-dimensional time dependent heat transfer through a ferrofluid in channel in response to magnetic force arising from the magnetic field of the magnetic dipole.
    [Show full text]
  • Performance Analysis of Continuity Equation and Its Applications
    EJERS, European Journal of Engineering Research and Science Vol. 1, No. 5, November 2016 Performance Analysis of Continuity Equation and Its Applications Md. Towhiduzzaman continuity equation is applicable over non-homogeneous Abstract— This study has been focused on the Continuity transferred conditions, one can carry out a research program Equation and its applications. A continuity equation illustrates concerning data collection of density, flow, and speed which the conveyances of some quantity. Principally, it is is essential for diverged three mid-blocks. This process of uncomplicated as well as powerful when applied to a conserved data compilation took place when some un-congested quantity. Contrarily, it can be generalized to apply to any extensive capacity. In this respect, firstly, it has been analyzed situations are prevailed with the moderation to slightly high on the Basic Hydrodynamics and the Continuity Equations. It volumes. Relating the regular density, obtained from has also been demonstrated on the applications of Continuity experiential densities in the field, reveals whether the Equation in which one can find out the different continuities continuity equation exactly calculates the average density for different areas of taps at different heights. Besides, there is under non-homogeneous traffic conditions or not [3]-[4]. a comparison also on the Theoretical Flow Rate and A variety of scholars explain this continuity equation in Experimental Flow Rate of fluid under the sluice gate by using Continuity Equation and Bernoulli Equation. So, it can be several ways. In 1903, Resc elucidated the convergence experimented here on the Theoretical Flow Rate and explanations of the channel flow as well as he calculated the Experimental Flow Rate to make a good agreement on the figure of the total errors in the continuity equation for the subject matters.
    [Show full text]
  • Introduction to Fluid Flow
    Stanford University Department of Aeronautics and Astronautics AA200A Applied Aerodynamics Chapter 1 - Introduction to fluid flow Stanford University Department of Aeronautics and Astronautics 1.1 Introduction Compressible flows play a crucial role in a vast variety of man-made and natural phenomena. Propulsion and power systems High speed flight Star formation, evolution and death Geysers and geothermal vents Earth meteor and comet impacts Gas processing and pipeline transfer Sound formation and propagation Stanford University Department of Aeronautics and Astronautics 1.2 Conservation of mass Stanford University Department of Aeronautics and Astronautics Divide through by the volume of the control volume. 1.2.1 Conservation of mass - Incompressible flow If the density is constant the continuity equation reduces to Note that this equation applies to both steady and unsteady incompressible flow Stanford University Department of Aeronautics and Astronautics 1.2.2 Index notation and the Einstein convention Make the following replacements Using index notation the continuity equation is Einstein recognized that such sums from vector calculus always involve a repeated index. For convenience he dropped the summation symbol. Coordinate independent form Stanford University Department of Aeronautics and Astronautics 1.3 Particle paths and streamlines in 2-D steady flow The figure below shows the streamlines over a 2-D airfoil. The flow is irrotational and incompressible Stanford University Department of Aeronautics and Astronautics Streamlines Streaklines Stanford University Department of Aeronautics and Astronautics A vector field that satisfies can always be represented as the gradient of a scalar potential or If the vector potential is substituted into the continuity equation the result is Laplaces equation.
    [Show full text]
  • Chapter 4 Continuity, Energy, and Momentum Equations
    Ch 4. Continuity, Energy, and Momentum Equation Chapter 4 Continuity, Energy, and Momentum Equations 4.1 Conservation of Matter in Homogeneous Fluids • Conservation of matter in homogeneous (single species) fluid → continuity equation positive dA1 - perpendicular to the control surface negative 4.1.1 Finite control volume method-arbitrary control volume - Although control volume remains fixed, mass of fluid originally enclosed (regions A+B) occupies the volume within the dashed line (regions B+C). Since mass m is conserved: mmm m (4.1) ABBtttdt Ctdt m m mmBBtdt t A t C tdt (4.2) dt dt •LHS of Eq. (4.2) = time rate of change of mass in the original control volume in the limit mmBB m tdt t B dV (4.3) dt t t CV 4−1 Ch 4. Continuity, Energy, and Momentum Equation where dV = volume element •RHS of Eq. (4.2) = net flux of matter through the control surface = flux in – flux out = qdA qdA nn12 where q = component of velocity vector normal to the surface of CV q cos n dV qdAnn12 qdA (4.4) CV CS CS t ※ Flux (= mass/time) is due to velocity of the flow. • Vector form dV qdA (4.5) t CV CS where dA = vector differential area pointing in the outward direction over an enclosed control surface qdAq dA cos positive for an outflow from cv, 90 negative for inflow into cv, 90 180 If fluid continues to occupy the entire control volume at subsequent times → time independent 4−2 Ch 4. Continuity, Energy, and Momentum Equation LHS: dV dV (4.5a) ttCV CV Eq.
    [Show full text]