CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics

Size: px
Start display at page:

Download "CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS. Multimedia Course on Continuum Mechanics"

Transcription

1 CH.9. CONSTITUTIVE EQUATIONS IN FLUIDS Multimedia Course on Continuum Mechanics

2 Overview Introduction Fluid Mechanics What is a Fluid? Pressure and Pascal s Law Constitutive Equations in Fluids Fluid Models Newtonian Fluids Constitutive Equations of Newtonian Fluids Relationship between Thermodynamic and Mean Pressures Components of the Constitutive Equation Stress, Dissipative and Recoverable Power Dissipative and Recoverable Powers Thermodynamic Considerations Limitations in the Viscosity Values Lecture 1 Lecture 3 Lecture 2 Lecture 4 Lecture 5 Lecture 6 2

3 9.1 Introduction Ch.9. Constitutive Equations in Fluids 3

4 What is a fluid? Fluids can be classified into: Ideal (inviscid) fluids: Also named perfect fluid. Only resists normal, compressive stresses (pressure). No resistance is encountered as the fluid moves. Real (viscous) fluids: Viscous in nature and can be subjected to low levels of shear stress. Certain amount of resistance is always offered by these fluids as they move. 5

5 9.2 Pressure and Pascal s Law Ch.9. Constitutive Equations in Fluids 6

6 Pascal s Law Pascal s Law: In a confined fluid at rest, pressure acts equally in all directions at a given point. 7

7 Consequences of Pascal s Law In fluid at rest: there are no shear stresses only normal forces due to pressure are present. The stress in a fluid at rest is isotropic and must be of the form: σ = p 1 σ ij 0 { } = p0 δ i, j 1, 2,3 ij Where p 0 is the hydrostatic pressure. 8

8 Pressure Concepts Hydrostatic pressure, fluid in equilibrium. p 0 : normal compressive stress exerted on a Mean pressure, p : minus the mean stress. 1 p = σm = Tr σ 3 ( ) REMARK Tr ( σ) is an invariant, thus, so are and p. σ m Thermodynamic pressure, p : Pressure variable used in the constitutive equations. It is related to density and temperature through the kinetic equation of state. REMARK F ( ρ, p, θ ) = 0 In a fluid at rest, p = p = p 0 9

9 Pressure Concepts Barotropic fluid: pressure depends only on density. ( ρ,p) = 0 = ( ρ) F p f Incompressible fluid: particular case of a barotropic fluid in which density is constant. ( ) ( ) F ρ,p, θ F ρ = ρ k = 0 ρ = k = const. 10

10 9.3 Constitutive Equations Ch.9. Constitutive Equations in Fluids 11

11 Reminder Governing Eqns. Governing equations of the thermo-mechanical problem: ρ+ ρ v = 0 Conservation of Mass. Continuity Equation. 1 eqn. σ + ρb= ρv Linear Momentum Balance. Cauchy s Motion Equation. 3 eqns. 8 PDE + 2 restrictions T σ = σ ρu = σ : d+ ρr q Angular Momentum Balance. Symmetry of Cauchy Stress Tensor. Energy Balance. First Law of Thermodynamics. 3 eqns. 1 eqn. ( u θ s ) σ :d 0 ρ + 1 q θ 0 2 ρθ Clausius-Planck Inequality. Heat flux Inequality. Second Law of Thermodynamics. 2 restrictions 19 scalar unknowns: ρ, v, σ, u, q, θ, s. 12

12 Reminder Constitutive Eqns. Constitutive equations of the thermo-mechanical problem: (, θ, ) σ = σ v ζ Thermo-Mechanical Constitutive Equations. 6 eqns. (19+p) PDE + (19+p) unknowns Entropy s = s( v, θ, ζ ) Constitutive Equation. 1 eqn. ( θ) K q= q = θ Thermal Constitutive Equation. Fourier s Law of Conduction. 3 eqns. i (,,, ) u = f ρ v θ ζ ( ρθ,, ζ ) = 0 { 1, 2,..., } F i p Caloric Kinetic State Equations. (1+p) eqns. set of new thermodynamic variables: ζ = ζ1, ζ 2,..., ζ p. { } The mechanical and thermal problem can be uncoupled if the temperature distribution is known a priori or does not intervene in the constitutive eqns. and if the constitutive eqns. involved do not introduce new thermodynamic variables. 13

13 Constitutive Equations Constitutive equations Together with the remaining governing equations, they are used to solve the thermo/mechanical problem. 14 In fluid mechanics, these are grouped into: Thermo-mechanical constitutive equations (, ρθ, ) ( d ) σ = p1 + f d σij = pδij + f ij, ρθ, i, j 1, 2, 3 Entropy constitutive equation (,, ) s = s d ρθ Fourier s Law q = k θ θ qi = k i, j 1, 2,3 xi Caloric equation of state ( ρθ) u = g, Kinetic equation of state F ( ρ, p, θ ) = 0 REMARK d v = ( ) s v

14 Viscous Fluid Models General form of the thermo-mechanical constitutive equations: σ = p1 + f d ij ij ij (, ρθ, ) f ( d,, ) i, j { 1, 2, 3} σ = pδ + ρθ Depending on the nature of f d, ρθ,, fluids are classified into : 1. Perfect fluid: ( ) ( ) f d, ρθ, = 0 σ = p1 2. Newtonian fluid: f is a linear function of the strain rate 3. Stokesian fluid: f is a non-linear function of its arguments 15

15 9.4. Newtonian Fluids Ch.9. Constitutive Equations in Fluids 16

16 Constitutive Equations of Newtonian Fluids Mechanic constitutive equations: σ = p1 + C :d σ where is the 4 th -order constant (viscous) constitutive tensor. Assuming: an isotropic medium the stress tensor is symmetrical Substitution of C into the constitutive equation gives: ( d) σ = p1+ λ Tr 1+ 2µ d { } σ = pδ + λd δ + 2 µ d i, j 1, 2,3 ij ij ll ij ij { } = pδ + C d i, j 1, 2,3 ij ij ijkl kl C C = λ µ I C ( ) = λδ δ + µ δ δ + δ δ ijkl ij kl ik jl il jk REMARK λ and µ { } i, jkl,, 1, 2, 3 are not necessarily constant. Both are a function of and. ρ θ 17

17 Relationship between Thermodynamic and Mean Pressures Taking the mechanic constitutive equation, Setting i=j, summing over the repeated index, and noting that δ ii = 3, we obtain 1 σ = 3p+ 3λ + 2µ d = 3p ( p = σ ) ii ( ) ll 3 p Tr( d) { } σ = pδ + λd δ + 2 µ d i, j 1, 2,3 ij ij ll ij ij 3 ii bulk viscosity κ 2 p = p + ( λ + µ ) Tr ( d) = p + κtr ( d 2 ) κ = λ + µ

18 Relationship between Thermodynamic and Mean Pressures Considering the continuity equation, dρ 1 dρ + ρ v= 0 v= dt ρ dt And the relationship Tr v x i ( d) = d = = v ii REMARK For a fluid at rest, v = 0 p= p = p 0 d For an incompressible fluid, ρ = 0 p = p dt For a fluid with, κ = 0 2 λ = µ p = p Stokes' 3 condition i p = p + κ Tr ( d) κ d ρ p = p+ κ v = p ρ dt 19

19 9.5 Components of the Constitutive Equations Ch.9. Constitutive Equations in Fluids 20

20 Components of the Constitutive Equation Given the Cauchy stress tensor, the following may be defined: σ = p1+ λ Tr ( d) 1+ 2µ d σ = σ sph + σ = p1 SPHERICAL PART mean pressure p = p κ v = p κtr ( d) 21 DEVIATORIC PART 1 λ ( d) 1 2µ d 1 σ σ= ( p p) 1+ Tr ( d) 1+ 2 p + Tr + = p + 2 σ= ( λ + µ ) Tr ( d) 1+ λ Tr ( d) 1+ 2µ d 3 1 σ=2 µ ( d ( ) ) 3 Tr d 1 = 2 µ d = d p = p κ Tr ( d) λ µ d 2 κ = λ + µ 3 deviatoric part of the rate of strain tensor

21 Components of the Constitutive Equation Given the Cauchy stress tensor, the following may be defined: SPHERICAL PART mean pressure p p = p κ v = p κ Tr ( d) p κ DEVIATORIC PART deviator stress tensor σ ij Tr ( d) 22 σ = 2µ d The stress tensor is then 1 σ = ( ) 3 Tr σ 1+ σ = p 1+ σ = 3p from the definition of mean pressure 2µ d ij REMARK κ Note that is not a function of d, while. µ = µ ( d)

22 9.6 Stress, Dissipative and Recoverable Powers Ch.9. Constitutive Equations in Fluids 23

23 Reminder Stress Power Mechanical Energy Balance: d 1 2 Pe ( t) = ρb vdv + ds ρ v dv dv V t v = + V dt 2 σ :d V V V t external mechanical power entering the medium d Pe ( t) = K ( t) + P dt kinetic energy σ stress power REMARK The stress power is the mechanical power entering the system which is not spent in changing the kinetic energy. It can be interpreted as the work per unit of time done by the stress in the deformation process of the medium. A rigid solid will have zero stress power. 24

24 Dissipative and Recoverable Powers Stress Power ( d) = σ :d dv 1 V d= ( ) 3 Tr d1+ d σ = p 1 + σ 1 σ: d= ( p1+ σ ) : Tr ( d) 1+ d = 3 = 3 = Tr ( d ) = = ptr ( d) 11 : + σ : d p1: d + Tr ( d) σ : 1= 3 3 = Tr σ = = ptr + σ : d σ = 2µ d 2 σ ptr ( ) Tr ( ) p = p κtr ( d) ( ) 0 : d= d + κ d + 2 µ d : d RECOVERABLE STRESS POWER,. W R DISSIPATIVE STRESS POWER,. 2W D 25

25 Dissipative and Recoverable Parts of the Cauchy Stress Tensor Associated to the concepts of recoverable and dissipative powers, the Cauchy stress tensor is split into: σ = p1+ λ Tr ( d) 1+ 2µ d RECOVERABLE PART,. σ R DISSIPATIVE PART,. σ D And the recoverable and dissipative powers are rewritten as: R ( ) 2 ( ) W = ptr d = p1 :d = σ :d 2 W = κ Tr d + 2µ d : d = σ :d D R D REMARK For an incompressible fluid, W = ptr d = 0 R ( ) 26

26 Thermodynamic considerations Specific recoverable stress power is an exact differential, 1 1 dg W R R :d ρ = σ ρ = dt (exact differential) Then, the recoverable stress work per unit mass in a closed cycle is zero: B A B A B A 1 1 W Rdt = σ R:d dt = dg = GB A GA = 0 ρ ρ A A A This justifies the denomination recoverable stress power. 28

27 Thermodynamic Considerations According to the 2 nd Law of Thermodynamics, the dissipative power is necessarily non-negative, ( ) = κ d + 2µ d d = d= 2 2WD 0 2W D Tr : 0 0 In a closed cycle, the work done by the dissipative stress per unit mass will, in general, be different to zero: B B 1 σ D :d dt ρ A 2WD > 0 > 0 This justifies the denomination dissipative power. 29

28 Limitations in the Viscosity Values The thermodynamic restriction, 2 κ Tr d 2µ d d ( ) 2W = + : 0 D introduces limitations in the values of the viscosity parameters κ, λ and µ : 1. For a purely spherical deformation rate tensor: Tr ( d) 0 d = 0 ( d) = κ 2 2WD Tr 0 2 κ = λ + µ For a purely deviatoric deformation rate tensor: Tr ( d) = 0 d 0 2WD = 2 µ d : d = 2µ dd ij ij 0 µ 0 > 0 30

29 Chapter 9 Constitutive Equations in Fluids 9.1 Concept of Pressure Several concepts of pressure are used in continuum mechanics (hydrostatic pressure, mean pressure and thermodynamic pressure) which, in general, do not coincide Hydrostatic Pressure Definition 9.1. Pascal s law In a confined fluid at rest, the stress state on any plane containing a given point is the same and is characterized by a compressive normal stress. In accordance with Pascal s law, the stress state of a fluid at rest is characterized by a stress tensor of the type Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar σ = p 0 1 σ ij = p 0 δ ij i, j {1,2,3}, (9.1) where p 0 is denoted as hydrostatic pressure (see Figure 9.1). Definition 9.2. The hydrostatic pressure is the compressive normal stress, constant on any plane, that acts on a fluid at rest. 439

30 440 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS Figure 9.1: Stress state of a fluid at rest. Figure 9.2: Mohr s circle of the stress tensor of a fluid at rest. Remark 9.1. The stress tensor of a fluid at rest is a spherical tensor and its representation in the Mohr s plane is a point (see Figure 9.2). Consequently, any direction is a principal stress direction and the stress state is constituted by the state defined in Section 4.8 of Chapter 4 as hydrostatic stress state Mean Pressure Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar Definition 9.3. The mean stress σ m is defined as σ m = 1 3 Tr (σ)=1 3 σ ii. The mean pressure p is defined as minus the mean stress, p def = mean pressure = σ m = 1 3 Tr (σ)= 1 3 σ ii. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

31 Concept of Pressure 441 Remark 9.2. In a fluid at rest, the mean pressure p coincides with the hydrostatic pressure p 0, σ = p 0 1 = σ m = 1 3 ( 3p 0)= p 0 = p = p 0. Generally, in a fluid in motion the mean pressure and the hydrostatic pressure do not coincide. Remark 9.3. The trace of the Cauchy stress tensor is a stress invariant. Consequently, the mean stress and the mean pressure are also stress invariants and, therefore, their values do not depend on the Cartesian coordinate system used Thermodynamic Pressure. Kinetic Equation of State A new thermodynamic pressure variable, named thermodynamic pressure and denoted as p, intervenes in the constitutive equations of fluids or gases. Definition 9.4. The thermodynamic pressure is the pressure variable that intervenes in the constitutive equations of fluids and gases, and is related to the density ρ and the absolute temperature θ by means of the kinetic equation of state, F (p,ρ,θ)=0. Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar Example 9.1 The ideal gas law is a typical example of kinetic equation of state: F (p,ρ,θ) p ρrθ = 0 = p = ρrθ, where p is the thermodynamic pressure and R is the universal gas constant. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

32 442 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS Remark 9.4. In a fluid at rest, the hydrostatic pressure p 0, the mean pressure p and the thermodynamic pressure p coincide. Fluid at rest : p 0 = p = p Generally, in a fluid in motion the hydrostatic pressure, the mean pressure and the thermodynamic pressure do not coincide. Remark 9.5. A barotropic fluid is defined by a kinetic equation of state in which the temperature does not intervene. Barotropic fluid : F (p,ρ)=0 = p = f (ρ) = ρ = g(p) Remark 9.6. An incompressible fluid is a particular case of barotropic fluid in which density is constant (ρ (x, t)= k = const.). In this case, the kinetic equation of state can be written as F (p,ρ,θ) ρ k = 0 and does not depend on the pressure or the temperature. Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar 9.2 Constitutive Equations in Fluid Mechanics Here, the set of equations, generically named constitutive equations, that must be added to the balance equations to formulate a problem in fluid mechanics (see Section 5.13 in Chapter 5) is considered. These equations can be grouped as follows: a) Thermo-mechanical constitutive equation This equation expresses the Cauchy stress tensor in terms of the other thermodynamic variables, typically the thermodynamic pressure p, the strain rate tensor d (which can be considered an implicit function of the velocity, d(v)= S v), the density ρ and the absolute temperature θ. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

33 Constitutive Equation in Viscous Fluids 443 Thermo-mechanical constitutive equation: σ = p1 + f(d,ρ,θ) 6 equations (9.2) b) Entropy constitutive equation An algebraic equation that provides the specific entropy s in terms of the strain rate tensor, the density and the absolute temperature. Entropy constitutive equation: s = s(d,ρ,θ) 1 equation (9.3) c) Thermodynamic constitutive equations or equations of state These are typically the caloric equation of state, which defines the specific internal energy u, and the kinetic equation of state, which provides an equation for the thermodynamic pressure. Caloric equation of state: Kinetic equation of state: d) Thermal constitutive equations u = g(ρ,θ) F (ρ, p,θ)=0 2 equations (9.4) The most common one is Fourier s law, which defines the heat flux by conduction q as q = k θ Fourier s θ 3 equations (9.5) law: q i = k ij i {1,2,3} x j Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar where k is the (symmetrical second-order) tensor of thermal conductivity, which is a property of the fluid. For the isotropic case, the thermal conductivity tensor is a spherical tensor k = k 1 and depends on the scalar parameter k, which is the thermal conductivity of the fluid. 9.3 Constitutive Equation in Viscous Fluids The general form of the thermo-mechanical constitutive equation (see (9.2)) for a viscous fluid is σ = p 1 + f(d,ρ,θ) σ ij = p δ ij + f ij (d,ρ,θ) i, j {1,2,3}, (9.6) X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

34 444 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS where f is a symmetrical tensor function. According to the character of the function f, the following models of fluids are defined: a) Stokesian or Stokes fluid: the function f is a non-linear function of its arguments. b) Newtonian fluid: the function f is a linear function of its arguments. c) Perfect fluid: the function f is null. In this case, the mechanical constitutive equation is σ = p1. In the rest of this chapter, only the cases of Newtonian and perfect fluids will be considered. Remark 9.7. The perfect fluid hypothesis is frequently used in hydraulic engineering, where the fluid under consideration is water. 9.4 Constitutive Equation in Newtonian Fluids The mechanical constitutive equation 1 for a Newtonian fluid is σ = p 1 + C : d σ ij = p δ ij + C ijkl d kl i, j {1,2,3}, (9.7) where C is a constant fourth-order (viscosity) constitutive tensor. A linear dependency of the stress tensor σ on the strain rate tensor d is obtained as a result of (9.7). For an isotropic Newtonian fluid, the constitutive tensor C is an isotropic fourth-order tensor. { C = λ μI C ijkl = λδ ij δ kl + μ ( ) (9.8) δ ik δ jl + δ il δ jk Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar i, j,k,l {1,2,3} Replacing (9.8) in the mechanical constitutive equation (9.7) yields σ = p 1 +(λ μI) : d = p 1 + λ Tr (d)1 + 2μ d, (9.9) which corresponds to the constitutive equation of an isotropic Newtonian fluid. 1 Note that the thermal dependencies of the constitutive equation are not considered here and, thus, the name mechanical constitutive equations. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

35 Constitutive Equation in Newtonian Fluids 445 { Constit. eqn. of σ = p 1 + λ Tr (d)1 + 2μ d an isotropic Newtonian fluid σ ij = p δ ij + λ d ll δ ij + 2μ d ij i, j {1,2,3} (9.10) Remark 9.8. Note the parallelism that can be established between the constitutive equation of a Newtonian fluid and that of a linear elastic solid (see Chapter 6): Newtonian fluid { σ = p 1 + C : d σ ij = p δ ij + C ijkl d kl Linear elastic solid { σ = C : ε σ ij = C ijkl ε kl Remark 9.9. The parameters λ and μ physically correspond to the viscosities, which are understood as material properties. In the most general case, they may not be constant and can depend on other thermodynamic variables, λ = λ (ρ,θ) and μ = μ (ρ,θ). A typical example is the dependency of the viscosity on the temperature in the form μ (θ) =μ 0 e α(θ θ 0), which establishes that the fluid s viscosity decreases as temperature increases (see Figure 9.3). Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar Figure 9.3: Possible dependency of the viscosity μ on the absolute temperature θ. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

36 446 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS Relation between the Thermodynamic and Mean Pressures In general, the thermodynamic pressure p and the mean pressure p in a Newtonian fluid in motion will be different but are related to each other. From the (mechanical) constitutive equation of a Newtonian fluid (9.10), σ = p 1 + λ Tr (d)1 + 2μ d = Tr (σ) = p Tr (1)+λ Tr (d)tr (1)+2μTr (d)= 3p +(3λ + 2μ)Tr (d) = }{{} 3 p p = p + (λ + 23 ) μ Tr (d)= p + K Tr (d) }{{} K (9.11) where K is denoted as bulk viscosity. Bulk viscosity : K = λ μ (9.12) Using the mass continuity equation (5.24), results in Then, considering the relation and replacing in (9.11), yields dρ dt + ρ v = 0 = v = 1 dρ ρ dt (9.13) Tr (d)=d ii = v i = v (9.14) x i Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar p = p + K v = p K ρ dρ dt (9.15) which relates the mean and thermodynamic pressures. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

37 Constitutive Equation in Newtonian Fluids 447 Remark In accordance with (9.15), the thermodynamic pressure and the mean pressure in a Newtonian fluid will coincide in the following cases: Fluid at rest: v = 0 = p = p = p 0 dρ Incompressible fluid: = 0 = p = p dt Fluid with null bulk viscosity K (Stokes condition 2 ): K = 0 = λ = 2 μ = p = p Constitutive Equation in Spherical and Deviatoric Components Spherical part From (9.15), the following relation is deduced. p = p K v = p K Tr (d) (9.16) Deviatoric part Using the decomposition of the stress tensor σ and the strain rate tensor d in its spherical and deviator components, and replacing in the constitutive equation (9.10), results in σ = 1 Tr (σ) 1 + σ = p1 + σ = p1 + λ Tr (d)1 + 2μ d 3 }{{} = σ = 3 p ( p p) 1 + λ Tr (d)1 + 2μd = ( ) λ K Tr (d)1 + 2μd }{{}}{{} K Tr (d) = Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar λ μ σ = 2 ( 3 μ Tr (d)1 + 2μd = 2μ d 1 ) 3 Tr (d)1 } {{ } d = (9.17) 2 Stokes condition is assumed in certain cases because the results it provides match the experimental observations. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

38 448 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS where (9.16) and (9.12) have been taken into account. σ = 2μd (9.18) Stress Power, Recoverable Power and Dissipative Power Using again the decomposition of the stress and strain rate tensors in their spherical and deviatoric components yields σ = p1 + σ and d = 1 3 Tr (d)1 + d, (9.19) and replacing in the expression of the stress power density (stress power per unit of volume) σ : d, results in 3 ( ) 1 σ : d =( p1 + σ ) : Tr (d)1 + d = 3 = 1 p Tr (d) 1 : 1 +σ : d p 1 : d }{{}}{{} 3 Tr ( d ) 3 Tr (d) σ : 1 = }{{} = 0 Tr ( σ ) (9.20) = 0 = p Tr (d)+σ : d. Replacing (9.16) and (9.17) in(9.20) produces ( ) σ : d = p K Tr (d) Tr (d)+2μd : d. (9.21) σ : d = p Tr (d) + KTr 2 (d)+2μd : d = W R + 2W D }{{}}{{} (9.22) recoverable power dissipative power W R 2W D Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar Recoverable power density: W R = p Tr (d) (9.23) Dissipative power density: 2W D = KTr 2 (d)+2μd : d 3 The property that the trace of a deviator tensor is null is used here. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

39 Constitutive Equation in Newtonian Fluids 449 Associated with the concepts of recoverable and dissipative powers, the recoverable and dissipative parts of the stress tensor, σ R and σ D, respectively, are defined as σ = p1 +λ Tr (d)1 + 2μ d = σ = σ R + σ D. (9.24) }{{}}{{} σ R σ D Using the aforementioned notation, the recoverable, dissipative and total power densities can be rewritten as W R = p Tr (d)= p1: d = σ R : d, 2W D = KTr 2 (d)+2μd : d = σ D : d, σ : d =(σ R + σ D ) : d = σ R : d + σ D : d = W R + 2W D. Remark In an incompressible fluid, the recoverable power is null. In effect, since the fluid is incompressible, dρ/dt = 0, and considering the mass continuity equation (5.24), v = 1 dρ ρ dt = 0 = Tr (d) = W R = p Tr (d)=0. Remark Introducing the decomposition of the stress power (9.25), the balance of mechanical energy (5.73) becomes P e = dk + σ : d dv = dk + σ R : d dv + σ D : d dv dt dt (9.25) Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar V V V P e = dk dt + V W R dv + V 2W D dv, which indicates that the mechanical power entering the fluid P e is invested in modifying the kinetic energy K and creating recoverable power and dissipative power. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

40 450 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS Thermodynamic Considerations 1) It can be proven that, under general conditions, the specific recoverable power (recoverable power per unit of mass) is an exact differential 1 ρ W R = 1 ρ σ R : d = dg. (9.26) dt In this case, the recoverable work per unit of mass performed in a closed cycle will be null (see Figure 9.4), B A A 1 ρ W R dt = B A A 1 ρ σ R : d dt = B A which justifies the denomination of W R as recoverable power. Figure 9.4: Closed cycle. A dg = G B A G A = 0, (9.27) 2) The second law of thermodynamics allows proving that the dissipative power 2W D in (9.25) is always non-negative, Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar 2W D 0 ; 2W D = 0 d = 0 (9.28) and, therefore, in a closed cycle the work performed per unit of mass by the dissipative stresses will, in general, not be null, B A 1 ρ σ D : d }{{} dt > 0. (9.29) 2W D > 0 This justifies the denomination of 2W D as (non-recoverable) dissipative power. The dissipative power is responsible for the dissipation (or loss of energy) phenomenon in fluids. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

41 Constitutive Equation in Newtonian Fluids 451 Example 9.2 Explain why an incompressible Newtonian fluid in motion that is not provided with external power (work per unit of time) tends to reduce its velocity to a complete stop. Solution The recoverable power in an incompressible fluid is null (see Remark 9.11). In addition, the dissipative power 2W D is known to be always non-negative (see (9.28)). Finally, applying the balance of mechanical energy (see Remark 9.12) results in 0 = P e = dk + W R dv + 2W D dv = dt }{{} V = 0 V dk = d 1 dt dt 2 ρv2 dv = 2 W D dv < 0 }{{} V V > 0 and, therefore, the fluid looses (dissipates) kinetic energy and the velocity of its particles decreases Limitations in the Viscosity Values Due to thermodynamic considerations, the dissipative power 2W D in (9.25) has been seen to always be non-negative, 2W D = KTr 2 (d)+2μ d : d 0. (9.30) This thermodynamic restriction introduces limitations in the admissible values of the viscosity parameters K, λ and μ of the fluid. In effect, given a certain fluid, the aforementioned restriction must be verified for all motions (that is, Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar for all velocity fields v) that the fluid may possibly have. Therefore, it must be verified for any arbitrary value of the strain rate tensor d = S (v). Consider, in particular, the following cases: a) The strain rate tensor d is a spherical tensor. In this case, from (9.30) results Tr (d) 0; d = 0 = 2W D = KTr 2 (d) 0 = K = λ μ 0 (9.31) such that only the non-negative values of the bulk viscosity K are feasible. X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

42 452 CHAPTER 9. CONSTITUTIVE EQUATIONS IN FLUIDS b) The strain rate tensor d is a deviatoric tensor. This type of flow is schematically represented in Figure 9.5. In this case, from (9.30) results Tr (d)=0; d 0 = 2W D = 2μ d : d = 2μ d ij : d ij 0 = }{{} > 0 (9.32) μ 0 1 v x v x (y) 0 0 v(x,y)= 0 ; d = 2 y 1 v x y = d Figure 9.5: Flow characterized by a deviatoric strain rate tensor. Continuum Mechanics for Engineers Theory and Problems X. Oliver and C. Agelet de Saracibar X. Oliver and C. Agelet de Saracibar Continuum Mechanics for Engineers.Theory and Problems doi: /rg

CH.5. BALANCE PRINCIPLES. Multimedia Course on Continuum Mechanics

CH.5. BALANCE PRINCIPLES. Multimedia Course on Continuum Mechanics CH.5. BALANCE PRINCIPLES Multimedia Course on Continuum Mechanics Overview Balance Principles Convective Flux or Flux by Mass Transport Local and Material Derivative of a olume Integral Conservation of

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

CH.1. THERMODYNAMIC FOUNDATIONS OF CONSTITUTIVE MODELLING. Computational Solid Mechanics- Xavier Oliver-UPC

CH.1. THERMODYNAMIC FOUNDATIONS OF CONSTITUTIVE MODELLING. Computational Solid Mechanics- Xavier Oliver-UPC CH.1. THERMODYNAMIC FOUNDATIONS OF CONSTITUTIVE MODELLING Computational Solid Mechanics- Xavier Oliver-UPC 1.1 Dissipation approach for constitutive modelling Ch.1. Thermodynamical foundations of constitutive

More information

Continuum Mechanics Fundamentals

Continuum Mechanics Fundamentals Continuum Mechanics Fundamentals James R. Rice, notes for ES 220, 12 November 2009; corrections 9 December 2009 Conserved Quantities Let a conseved quantity have amount F per unit volume. Examples are

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

Chapter 2: Fluid Dynamics Review

Chapter 2: Fluid Dynamics Review 7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading

More information

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature

Chapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte

More information

In this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case

In this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field

More information

Physical Conservation and Balance Laws & Thermodynamics

Physical Conservation and Balance Laws & Thermodynamics Chapter 4 Physical Conservation and Balance Laws & Thermodynamics The laws of classical physics are, for the most part, expressions of conservation or balances of certain quantities: e.g. mass, momentum,

More information

CH.6. LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics

CH.6. LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics CH.6. LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Hypothesis of the Linear Elasticity Theory Linear Elastic Constitutive Equation Generalized Hooke s Law Elastic Potential Isotropic

More information

Navier-Stokes Equation: Principle of Conservation of Momentum

Navier-Stokes Equation: Principle of Conservation of Momentum Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation

More information

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms

Continuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive

More information

Math 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations

Math 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The

More information

CH.3. COMPATIBILITY EQUATIONS. Multimedia Course on Continuum Mechanics

CH.3. COMPATIBILITY EQUATIONS. Multimedia Course on Continuum Mechanics CH.3. COMPATIBILITY EQUATIONS Multimedia Course on Continuum Mechanics Overview Introduction Lecture 1 Compatibility Conditions Lecture Compatibility Equations of a Potential Vector Field Lecture 3 Compatibility

More information

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)

- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum) 2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities

More information

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for

More information

Summary of the Equations of Fluid Dynamics

Summary of the Equations of Fluid Dynamics Reference: Summary of the Equations of Fluid Dynamics Fluid Mechanics, L.D. Landau & E.M. Lifshitz 1 Introduction Emission processes give us diagnostics with which to estimate important parameters, such

More information

CH.11. VARIATIONAL PRINCIPLES. Multimedia Course on Continuum Mechanics

CH.11. VARIATIONAL PRINCIPLES. Multimedia Course on Continuum Mechanics CH.11. ARIATIONAL PRINCIPLES Multimedia Course on Continuum Mechanics Overview Introduction Functionals Gâteaux Derivative Extreme of a Functional ariational Principle ariational Form of a Continuum Mechanics

More information

Course Syllabus: Continuum Mechanics - ME 212A

Course Syllabus: Continuum Mechanics - ME 212A Course Syllabus: Continuum Mechanics - ME 212A Division Course Number Course Title Academic Semester Physical Science and Engineering Division ME 212A Continuum Mechanics Fall Academic Year 2017/2018 Semester

More information

Differential relations for fluid flow

Differential relations for fluid flow Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow

More information

CH.4. STRESS. Continuum Mechanics Course (MMC)

CH.4. STRESS. Continuum Mechanics Course (MMC) CH.4. STRESS Continuum Mechanics Course (MMC) Overview Forces Acting on a Continuum Body Cauchy s Postulates Stress Tensor Stress Tensor Components Scientific Notation Engineering Notation Sign Criterion

More information

CH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics

CH.7. PLANE LINEAR ELASTICITY. Multimedia Course on Continuum Mechanics CH.7. PLANE LINEAR ELASTICITY Multimedia Course on Continuum Mechanics Overview Plane Linear Elasticit Theor Plane Stress Simplifing Hpothesis Strain Field Constitutive Equation Displacement Field The

More information

INDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226

INDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226 INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence

More information

MHA042 - Material mechanics: Duggafrågor

MHA042 - Material mechanics: Duggafrågor MHA042 - Material mechanics: Duggafrågor 1) For a static uniaxial bar problem at isothermal (Θ const.) conditions, state principle of energy conservation (first law of thermodynamics). On the basis of

More information

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

Contents. I Introduction 1. Preface. xiii

Contents. I Introduction 1. Preface. xiii Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................

More information

AA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations

AA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental

More information

Continuum mechanism: Stress and strain

Continuum mechanism: Stress and strain Continuum mechanics deals with the relation between forces (stress, σ) and deformation (strain, ε), or deformation rate (strain rate, ε). Solid materials, rigid, usually deform elastically, that is the

More information

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 5

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 5 .9 Numerical Fluid Mechanics Fall 011 Lecture 5 REVIEW Lecture 4 Roots of nonlinear equations: Open Methods Fixed-point Iteration (General method or Picard Iteration), with examples Iteration rule: x g(

More information

CHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:

CHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation: CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have

More information

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet

Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet Dynamics of the Mantle and Lithosphere ETH Zürich Continuum Mechanics in Geodynamics: Equation cheat sheet or all equations you will probably ever need Definitions 1. Coordinate system. x,y,z or x 1,x

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

Mechanics of solids and fluids -Introduction to continuum mechanics

Mechanics of solids and fluids -Introduction to continuum mechanics Mechanics of solids and fluids -Introduction to continuum mechanics by Magnus Ekh August 12, 2016 Introduction to continuum mechanics 1 Tensors............................. 3 1.1 Index notation 1.2 Vectors

More information

Elements of Rock Mechanics

Elements of Rock Mechanics Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider

More information

1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity

1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity 1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses

More information

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

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Review of fluid dynamics

Review of fluid dynamics Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence

More information

4 Constitutive Theory

4 Constitutive Theory ME338A CONTINUUM MECHANICS lecture notes 13 Tuesday, May 13, 2008 4.1 Closure Problem In the preceding chapter, we derived the fundamental balance equations: Balance of Spatial Material Mass ρ t + ρ t

More information

Engineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.

Engineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February. Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems

More information

Computer Fluid Dynamics E181107

Computer Fluid Dynamics E181107 Computer Fluid Dynamics E181107 2181106 Transport equations, Navier Stokes equations Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav

More information

Computational Astrophysics

Computational Astrophysics Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture

More information

Constitutive models: Incremental plasticity Drücker s postulate

Constitutive models: Incremental plasticity Drücker s postulate Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager

More information

Introduction to Fluid Mechanics

Introduction to Fluid Mechanics Introduction to Fluid Mechanics Tien-Tsan Shieh April 16, 2009 What is a Fluid? The key distinction between a fluid and a solid lies in the mode of resistance to change of shape. The fluid, unlike the

More information

Basic hydrodynamics. David Gurarie. 1 Newtonian fluids: Euler and Navier-Stokes equations

Basic hydrodynamics. David Gurarie. 1 Newtonian fluids: Euler and Navier-Stokes equations Basic hydrodynamics David Gurarie 1 Newtonian fluids: Euler and Navier-Stokes equations The basic hydrodynamic equations in the Eulerian form consist of conservation of mass, momentum and energy. We denote

More information

Other state variables include the temperature, θ, and the entropy, S, which are defined below.

Other state variables include the temperature, θ, and the entropy, S, which are defined below. Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive

More information

Vít Průša (joint work with K.R. Rajagopal) 30 October Mathematical Institute, Charles University

Vít Průša (joint work with K.R. Rajagopal) 30 October Mathematical Institute, Charles University On models for viscoelastic fluid-like materials that are mechanically incompressible and thermally compressible or expansible and their Oberbeck Boussinesq type approximations Vít Průša (joint work with

More information

1 Gauss integral theorem for tensors

1 Gauss integral theorem for tensors Non-Equilibrium Continuum Physics TA session #1 TA: Yohai Bar Sinai 16.3.216 Index Gymnastics: Gauss Theorem, Isotropic Tensors, NS Equations The purpose of today s TA session is to mess a bit with tensors

More information

Hydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010

Hydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal

More information

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations 6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =

More information

Quick Recapitulation of Fluid Mechanics

Quick Recapitulation of Fluid Mechanics Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change

More information

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.

By drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ. Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,

More information

CONTINUUM MECHANICS. lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern

CONTINUUM MECHANICS. lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern CONTINUUM MECHANICS lecture notes 2003 jp dr.-ing. habil. ellen kuhl technical university of kaiserslautern Contents Tensor calculus. Tensor algebra.................................... Vector algebra.................................

More information

A short review of continuum mechanics

A short review of continuum mechanics A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material

More information

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion

Conservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation

More information

Notes 4: Differential Form of the Conservation Equations

Notes 4: Differential Form of the Conservation Equations Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.

More information

CH.1. DESCRIPTION OF MOTION. Multimedia Course on Continuum Mechanics

CH.1. DESCRIPTION OF MOTION. Multimedia Course on Continuum Mechanics CH.1. DESCRIPTION OF MOTION Multimedia Course on Continuum Mechanics Overview 1.1. Definition of the Continuous Medium 1.1.1. Concept of Continuum 1.1.2. Continuous Medium or Continuum 1.2. Equations of

More information

7 The Navier-Stokes Equations

7 The Navier-Stokes Equations 18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and

More information

Chapter 3 Stress, Strain, Virtual Power and Conservation Principles

Chapter 3 Stress, Strain, Virtual Power and Conservation Principles Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress

More information

From the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)

From the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1) Lecture 9 9. Administration None. 9. Continuation of energy equation From the last time, we ended with an expression for the energy equation ρ D (e + ) u = ρg u + (τ u) q (9.) Where ρg u changes in potential

More information

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3

Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso

More information

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow

AA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit

More information

2 Law of conservation of energy

2 Law of conservation of energy 1 Newtonian viscous Fluid 1 Newtonian fluid For a Newtonian we already have shown that σ ij = pδ ij + λd k,k δ ij + 2µD ij where λ and µ are called viscosity coefficient. For a fluid under rigid body motion

More information

KINEMATICS OF CONTINUA

KINEMATICS OF CONTINUA KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation

More information

The Navier-Stokes Equations

The Navier-Stokes Equations s University of New Hampshire February 22, 202 and equations describe the non-relativistic time evolution of mass and momentum in fluid substances. mass density field: ρ = ρ(t, x, y, z) velocity field:

More information

Continuum Mechanics. Continuum Mechanics and Constitutive Equations

Continuum Mechanics. Continuum Mechanics and Constitutive Equations 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

More information

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace

Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and

More information

Dynamics of Ice Sheets and Glaciers

Dynamics of Ice Sheets and Glaciers Dynamics of Ice Sheets and Glaciers Ralf Greve Institute of Low Temperature Science Hokkaido University Lecture Notes Sapporo 2004/2005 Literature Ice dynamics Paterson, W. S. B. 1994. The Physics of

More information

CHAPTER 8 ENTROPY GENERATION AND TRANSPORT

CHAPTER 8 ENTROPY GENERATION AND TRANSPORT CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?

More information

Entropy generation and transport

Entropy generation and transport Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?

More information

Lecture 8: Tissue Mechanics

Lecture 8: Tissue Mechanics Computational Biology Group (CoBi), D-BSSE, ETHZ Lecture 8: Tissue Mechanics Prof Dagmar Iber, PhD DPhil MSc Computational Biology 2015/16 7. Mai 2016 2 / 57 Contents 1 Introduction to Elastic Materials

More information

Detailed Outline, M E 521: Foundations of Fluid Mechanics I

Detailed Outline, M E 521: Foundations of Fluid Mechanics I Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic

More information

( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #

( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint.  =  x,t,#, # Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue

More information

1. Background. is usually significantly lower than it is in uniaxial tension

1. Background. is usually significantly lower than it is in uniaxial tension NOTES ON QUANTIFYING MODES OF A SECOND- ORDER TENSOR. The mechanical behavior of rocks and rock-like materials (concrete, ceramics, etc.) strongly depends on the loading mode, defined by the values and

More information

Physical Law Eulerian Form N Lagrangian Form n = V, 3,

Physical Law Eulerian Form N Lagrangian Form n = V, 3, Chapter 5 Constitutive Modelling 5.1 Introduction Thus far we have established four conservation (or balance) equations, an entropy inequality and a number of kinematic relationships. The overall set of

More information

Chapter 16: Elastic Solids

Chapter 16: Elastic Solids Chapter 16: Elastic Solids Chapter 16: Elastic Solids... 366 16.1 Introduction... 367 16.2 The Elastic Strain... 368 16.2.1 The displacement vector... 368 16.2.2 The deformation gradient... 368 16.2.3

More information

Basic concepts to start Mechanics of Materials

Basic concepts to start Mechanics of Materials Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen

More information

Physics of Continuous media

Physics of Continuous media Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had

More information

Dynamics of Glaciers

Dynamics of Glaciers Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers

More information

Chapter 1 Fluid Characteristics

Chapter 1 Fluid Characteristics Chapter 1 Fluid Characteristics 1.1 Introduction 1.1.1 Phases Solid increasing increasing spacing and intermolecular liquid latitude of cohesive Fluid gas (vapor) molecular force plasma motion 1.1.2 Fluidity

More information

Mathematical Background

Mathematical Background CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous

More information

in this web service Cambridge University Press

in this web service Cambridge University Press CONTINUUM MECHANICS This is a modern textbook for courses in continuum mechanics. It provides both the theoretical framework and the numerical methods required to model the behavior of continuous materials.

More information

Chapter 5. The Differential Forms of the Fundamental Laws

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

Modeling viscous compressible barotropic multi-fluid flows

Modeling viscous compressible barotropic multi-fluid flows Modeling viscous compressible barotropic multi-fluid flows Alexander Mamontov, Dmitriy Prokudin 1 arxiv:1708.07319v2 [math.ap] 25 Aug 2017 August 25, 2017 Lavrentyev Institute of Hydrodynamics, Siberian

More information

Constitutive Equations

Constitutive Equations Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module

More information

Localization in Undrained Deformation

Localization in Undrained Deformation Localization in Undrained Deformation J. W. Rudnicki Dept. of Civil and Env. Engn. and Dept. of Mech. Engn. Northwestern University Evanston, IL 621-319 John.Rudnicki@gmail.com Fourth Biot Conference on

More information

2. Conservation Equations for Turbulent Flows

2. Conservation Equations for Turbulent Flows 2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging

More information

L8. Basic concepts of stress and equilibrium

L8. Basic concepts of stress and equilibrium L8. Basic concepts of stress and equilibrium Duggafrågor 1) Show that the stress (considered as a second order tensor) can be represented in terms of the eigenbases m i n i n i. Make the geometrical representation

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

Professor George C. Johnson. ME185 - Introduction to Continuum Mechanics. Midterm Exam II. ) (1) x

Professor George C. Johnson. ME185 - Introduction to Continuum Mechanics. Midterm Exam II. ) (1) x Spring, 997 ME85 - Introduction to Continuum Mechanics Midterm Exam II roblem. (+ points) (a) Let ρ be the mass density, v be the velocity vector, be the Cauchy stress tensor, and b be the body force per

More information

Continuum Mechanics Lecture 4 Fluid dynamics

Continuum Mechanics Lecture 4 Fluid dynamics Continuum Mechanics Lecture 4 Flui ynamics Prof. http://www.itp.uzh.ch/~teyssier Outline - Flui kinematics - Mass an momentum conservation laws - The energy equation - Real fluis - Ieal fluis - Incompressible

More information

A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations

A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm

More information

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

More information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,

More information

Chapter 9: Differential Analysis of Fluid Flow

Chapter 9: Differential Analysis of Fluid Flow of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known

More information

Continuum Mechanics Lecture 5 Ideal fluids

Continuum Mechanics Lecture 5 Ideal fluids Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation

More information

Level Set Tumor Growth Model

Level Set Tumor Growth Model Level Set Tumor Growth Model Andrew Nordquist and Rakesh Ranjan, PhD University of Texas, San Antonio July 29, 2013 Andrew Nordquist and Rakesh Ranjan, PhD (University Level Set of Texas, TumorSan Growth

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost 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

More information

Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount

Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic

More information

Subduction II Fundamentals of Mantle Dynamics

Subduction II Fundamentals of Mantle Dynamics 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 η

More information