Conservation of mass momentum and energy in fluid

  • Slides: 46
Download presentation
Conservation of mass, momentum and energy in fluid dynamics

Conservation of mass, momentum and energy in fluid dynamics

Classical mechanics (Newtonian) The motion of a point mass is governed by Newton’s second

Classical mechanics (Newtonian) The motion of a point mass is governed by Newton’s second law This is a set of ordinary differential equations (ODEs):

Classical mechanics (Newtonian) The equations of motion are supplemented with the initial conditions (for

Classical mechanics (Newtonian) The equations of motion are supplemented with the initial conditions (for the position and velocity of the point mass at the initial time t = t 0) that is

Classical mechanics (Newtonian) Example: One-body problem: determine the motion of a point mass (m)

Classical mechanics (Newtonian) Example: One-body problem: determine the motion of a point mass (m) moving in the gravity force field produces by a stationary mass (M). Let the mass M is placed in the center of a coordinate system XYZ. Then the equation of motion for the point mass m is This is a set of the following equations in 3 D: In fact, this system can be reduced to 2 D (planar) problem, because the trajectory (x(t), y(t), z(t)) lies in the plane defined by the initial position and velocity.

Classical mechanics (Newtonian) Example (cond. ) Possible orbits, due to universal gravity, of a

Classical mechanics (Newtonian) Example (cond. ) Possible orbits, due to universal gravity, of a small mass around a single large mass.

Classical mechanics (Newtonian) Energy in classical mechanics appears as follows: when the force is

Classical mechanics (Newtonian) Energy in classical mechanics appears as follows: when the force is conservative (rot F = 0), then it can be expressed as the gradient of a scalar function V This function is called the potential of the force or simply the potential energy. It is also denoted by Ep. Then the mechanical energy E is defined as the sum of the kinetic and potential energy Important property of E is that this quantity is preserved: for any motion r(t) in the force field F it holds that We say that the mechanical energy is conserved.

Classical mechanics (Newtonian) Example: Motion of a point mass near the Earth surface hence

Classical mechanics (Newtonian) Example: Motion of a point mass near the Earth surface hence The notion of such conserved quantity helps solving many problems. For example, what is the maximum height that a stone thrown vertically upward will attain if the initial velocity is v 0?

Classical mechanics (Newtonian) Example: In the motion of a point mass (m) in the

Classical mechanics (Newtonian) Example: In the motion of a point mass (m) in the gravitational field of a big stationary planet/star (M) the potential energy of the mass m is hence In this particular case the mass m disappears as we can divide by it to obtain Even if we do not know the solution r(t) we know that the above quantity must be constant with respect to time durrng the motion.

Classical mechanics (Newtonian) If a constant force is applied to a particle and a

Classical mechanics (Newtonian) If a constant force is applied to a particle and a displacement is r , the work is defined as In general case when the force may vary in space the work is defined by the line integral where AB is a path along which the mass is moved from point A to B. Note that in general this definition implies that the work may depend on the path. This is not true if the force is conservative. In the case of conservative force we have

Fluid mechanics (concept of a continuum) Materials (solids, liquids and gases) are composed of

Fluid mechanics (concept of a continuum) Materials (solids, liquids and gases) are composed of molecules separated by empty space. But the continuum model as a mathematical concept assumes that material exists as a continuous entity. It means that the matter in the body is continuously distributed and fills the entire region of space it occupies. A continuum body can, for example, be infinitely sub-divided into smaller and smaller elements preserving properties of the bulk material. It also means that two points in such body may be close at arbitrary small distance. Due to these assumption we can consider physical quantities (density, pressure, velocity, forces etc. ) as defined at every point as continuous functions of space position and that space derivatives can be also defined. Example (density). A function (x, y, z, t) is defined at every point (x, y, z) and is continuous. The meaning of this function is as follows: (x, y, z) W is the mass contained in the region W.

Fluid mechanics Let be a region in two dimensional (2 D) or three dimensional

Fluid mechanics Let be a region in two dimensional (2 D) or three dimensional (3 D) space filled with a fluid (gas or liquid). Our object is to describe the motion of such a fluid. Let r ∈ be a point and consider the particle of fluid moving through r at time t. In a standard Euclidean coordinates system in space, we write r = (x, y, z). This particle traverses a well-defined trajectory r(t)=(x(t), y(t), z(t). Let v(r, t) denote the velocity of the particle of fluid that is moving through r at time t. Thus, for each fixed time, v is a vector field on , as in figure below. We call v the (spatial) velocity field of the fluid.

Fluid mechanics (mass conservation) We will use the notion of some region W which

Fluid mechanics (mass conservation) We will use the notion of some region W which will move according to the velocity field. The symbol Wt means the set of points of the initial set (region) W displaced to new location after time t. Thus we have W 0 = W. The obvious property of such defined set Wt is This formula is just a statement of the mass conservation law. It is also called the continuity equation or mass balance equation. Let us notice that it can also be written as

Fluid mechanics (mass conservation) The form of mass conservation law from the previous slide

Fluid mechanics (mass conservation) The form of mass conservation law from the previous slide may be converted to more useful differential form. It can be done in two ways: (1) from the above integral form by applying Reynold’s transport theorem; (2) By direct application of Gauss’ divergence theorem. Th final form is Some textbooks use the notation for the divergence operator: u = div u.

Fluid mechanics (material derivative) The material derivative of a scalar function f(t, r) in

Fluid mechanics (material derivative) The material derivative of a scalar function f(t, r) in the region filled with a moving fluid with the velocity field v(t, r) is defined as It takes into account the fact that the fluid is moving and that the positions of fluid particles change with time. Thus, if r(t) = [x(t), y(t), z(t)] describes the movement of a fluid particle (the path followed by the fluid particle), then the change of any quantity f along this path can be computed as follows Hence we can write for the change of function f(t, r) along the path r(t):

Fluid mechanics (ideal fluid) For any continuum, forces acting on a piece of material

Fluid mechanics (ideal fluid) For any continuum, forces acting on a piece of material are of two types: • External, or body, forces such as gravity, intertial (pseudoforce) , a magnetic field, which exert a force per unit volume on the continuum • forces of stress, whereby the piece of material is acted on by forces across its surface by the rest of the continuum. An ideal fluid has the following property: For any motion of the fluid there is a scalar function p(t, r) called the pressure such that if S is a surface in the fluid with a normal vector n, the force of stress exerted across the surface S per unit area at r ∈ S at time t is p(t, r)n, i. e.

Fluid mechanics (ideal fluid – momentum equation) Intuitively, the absence of tangential forces implies

Fluid mechanics (ideal fluid – momentum equation) Intuitively, the absence of tangential forces implies that there is no way for rotation to start in a fluid, nor, if it is there at the beginning, to stop. This may be too much simplification and we can detect physical trouble for ideal fluids because of the abundance of rotation in real fluids (near the oars of a rowboat, in tornadoes, etc. ). If W is a region in the fluid at a particular instant of time t, the total force exerted on the fluid inside W by means of stress on its boundary is If e is any fixed vector in space, then the divergence theorem gives hence

Fluid mechanics (ideal fluid – momentum equation) If b(t, r) denotes the given body

Fluid mechanics (ideal fluid – momentum equation) If b(t, r) denotes the given body force per unit mass, then the total body force is Thus the total force acting on the fluid domain W in the case of ideal fluid is Hence

Fluid mechanics (ideal fluid) Now we want to apply Newton’s second law for the

Fluid mechanics (ideal fluid) Now we want to apply Newton’s second law for the fluid element W force = mass × acceleration Acceleration: let r(t) is the motion of the fluid particle. From the very definition of the velocity we get The acceleration a=d 2 r/dt 2 thus But, by the chain rule one can show

Fluid mechanics (ideal fluid – Euler’s equations) Let us summarize: From the above relations

Fluid mechanics (ideal fluid – Euler’s equations) Let us summarize: From the above relations we get The above equation is called Euler’s equation for ideal fluid. In fact it is the momentum balance for the ideal fluid.

Fluid mechanics (Navier Stokes equations) Now we consider a more general fluid than the

Fluid mechanics (Navier Stokes equations) Now we consider a more general fluid than the ideal one. Real fluids possess the phenomenon of viscosity. The microscopic basis of the effect of viscosity is connected with the transfer of momentum between layers of moving fluid. If the forces are only normal to S, there will be no transfer of momentum between the fluid volumes denoted by B 1 and B 2. However, in reality faster molecules from above (B 2) will diffuse across S and impart momentum to the fluid (B 1), and, likewise, slower molecules from below S will diffuse across S to slow down the fluid above S. This is how the viscosity arises.

Quantifying the dynamic viscosity In the previous slide the viscosity was explained in the

Quantifying the dynamic viscosity In the previous slide the viscosity was explained in the qualitative way. Here is how we can measure it. The dynamic viscosity (sometimes ) of a fluid describes its resistance to shearing flows. In the picture (the speed of the top plate is parallel to the x axis. We see that an external force F is required in order to keep the top plate moving at constant speed v.

Fluid mechanics (Navier Stokes equations) Instead of assuming that (ideal fluid) where n is

Fluid mechanics (Navier Stokes equations) Instead of assuming that (ideal fluid) where n is the normal vector to S, we now assume that (viscous fluid) where is the stress tensor which describes the shear forces (tensions) in the fluid. It is called the Cauchy stress tensor. In a fixed coordinate system the tensor stress is just a matrix 3 x 3. Hence the expression n should be understood as the matrix multiplication (n is written here as a column). From this formulation we see that shear force term is linear with respect to n. It is not obvious! It can be proved under a general assumption (continuous dependence on n) and using balance of momentum. This result is called Cauchy’s theorem.

Fluid mechanics (Navier Stokes equations) As in the case of the ideal fluid we

Fluid mechanics (Navier Stokes equations) As in the case of the ideal fluid we can formulate for viscous fluid Newton’s second law for the moving fluid particle Wt where d. V indicates the volume integral and d. A is the surface integral. The internal forces acting through surface are now composed of two terms: normal (represented by scalar function p – the pressure) and shear/tangential (represented by matrix function – the stress tensor). It can be also proved that the stress tensor is symmetric.

Fluid mechanics (Navier Stokes equations) To obtain a more specific equation for the motion

Fluid mechanics (Navier Stokes equations) To obtain a more specific equation for the motion of the fluid we have to analyze in more detail the stress tensor . The following assumptions are physically sound: 1. depends linearly on the velocity gradients v. 2. is invariant under rotations in space. 3. is symmetric (as was pointed out earlier this can be deduce as a consequence of balance of angular momentum. The velocity gradient is just a matrix of first order space derivatives of all componets of the velocity v=[vx, vy, vz]

Fluid mechanics (Navier Stokes equations) Under the assumptions 1, 2, 3 (previous slide) it

Fluid mechanics (Navier Stokes equations) Under the assumptions 1, 2, 3 (previous slide) it can be shown that the stress tensor has the following form Where and are non-negative parameters characterazing the viscous fluid and

Fluid mechanics (Navier Stokes equations) If we now employ the transport theorem and the

Fluid mechanics (Navier Stokes equations) If we now employ the transport theorem and the divergence theorem to the balance of momentum with the stress tensor = (div v)I + 2 D after some algebra we will obtain the famous Navier-Stokes equations (NS): This form is for a general fluid – at that stage we do not assume incompressibility, so the unknowns are: v(t, r), p(t, r). Thus we have five unknown functions, but the NS system gives only three equations. One more equation is obvious – it is the mass balance equation The fifth equation (for compressible fluid) comes usually in the form of a constitutive relation between pressure, density and temperature:

Fluid mechanics (incompressible Navier Stokes equations) Assumption: the viscous fluid is incompressible, i. e.

Fluid mechanics (incompressible Navier Stokes equations) Assumption: the viscous fluid is incompressible, i. e. = const. The mass balanace equation gives now and the term ( + ) (div v) vanishes from the NS system. Thus we obtain the Navier Stokes equations for incompressible fluid (usually written with the kinematic viscosity = / >0):

Fluid mechanics boundary conditions We will consider only two basic boundary conditions. Boundary is

Fluid mechanics boundary conditions We will consider only two basic boundary conditions. Boundary is denoted by . 1) The expression means that the normal component of the velocity on the boundary is zero. It is rather obvious statement that the wall is impenetrable. The tangential component may be non-zero (fluid can slip along the wall. This condition is consistent with the inviscid fluid and is used with Euler’sequation. 2) This condition says that the velocity on the boundary is completely zero. The tangential component must vanish, hence this is called non-slip boundary condition. It is consistent with viscous fluid ( > 0) and is used in the Navier. Stokes equations.

Fluid mechanics (energy conservation) For any moving fluid region of volume Wt we can

Fluid mechanics (energy conservation) For any moving fluid region of volume Wt we can define its energy as where is the density, v is the velocity field, and e is the density of internal energy. We identify here the kinetic energy of the moving volume Wt the internal energy of the moving volume Wt Now we are interested in how the quantity Wt changes with time when the fluid particle moves in the system. EWt

Fluid mechanics (energy) If we want to move the time derivative under the symbol

Fluid mechanics (energy) If we want to move the time derivative under the symbol of integral we must be cautious because the region of integration Wt also depend on time. It is not stationary! The rate of change with time can be rewritten Reynold’s transport theorem where we applied the transport theorem (Reynolds) valid for any quantity f(t, r)=f(t, x, y, z) that may depend on time and location. And moving region Wt. Here, Df/Dt is defined earlier the material derivative.

Fluid mechanics (energy conservation) The previous slide described a defintion of energy of a

Fluid mechanics (energy conservation) The previous slide described a defintion of energy of a moving element and how its change can be written with derivative under the integral. But these were merely mathematical manipulations of purely kinematical nature. Now we have to add some constitutive relations to bind the energy of the moving element Wt with other processes present in the system. This leads to the following general form of energy balance in the integral form: where: rate of work by (external) mass forces rate of work by (internal) pressure forces flux of the ”heat” energy gain/loss due the internal sources/sinks (e. g. reactions)

Fluid mechanics (energy conservation) The integral form: can be transformed to differential if we

Fluid mechanics (energy conservation) The integral form: can be transformed to differential if we apply some formulas of vector analysis. where the basic tool – Gauss’ divergence theorem was used: Here u stand for any vector field in space (3 D or 2 D).

Fluid mechanics (energy conservation) The previous slide contained the expression where a major assumption

Fluid mechanics (energy conservation) The previous slide contained the expression where a major assumption concerning the nature of the internal stress p was made. First, p=p(t, r, n) which means that it not only depend on time and location – but also on the normal to the surface! Second, it can be proved under some broad assumptions that it is a linear function of n, hence: where is called the stress tensor. It is symmetric tensor that can be expressed as a matrix

Fluid mechanics (energy conservation) Now we write all integrals on the right hand side

Fluid mechanics (energy conservation) Now we write all integrals on the right hand side as volume integrals and use the transport theorem Because this integral equality holds for any domain Wt we obtain following general form of energy conservation equation

Heat transport – elementary derivation The heat flux Jq in direction Ox. A domain

Heat transport – elementary derivation The heat flux Jq in direction Ox. A domain is in the form of a parallelepiped with edges parallel to the coordinate axes. cross-sectional are perpendicular to the flux (m 2) density of the body (kg/m 3) specific heat of the material (J/kg m 3) volumetric heat sources (J/m 3 s)] Basic definitions and properties: Change of internal energy (heat) with the change of temperature T Mass of a body with density and volume V The mount of energy which flows through the surface of unit area in unit time in the direction perpendicular to the surface The amount of energy which flows through area A in time t in a direction perpendicular to A The amount of energy produced/consumed in the volume of V in time t

Jq(x+ x, t) Qcałk = m cw T Jq(x, t) Qźr = S(x, t)

Jq(x+ x, t) Qcałk = m cw T Jq(x, t) Qźr = S(x, t) V t x x+ x Balance of energy in the domain in time t: Qtotal = (the flux through boundary) + (sources) m cw T = Jq(x, t) A t - Jq(x+x, t) A t + S(x, t) V t A x cw T = Jq(x, t) A t - Jq(x+x, t) A t + S(x, t) A x t Divide both sides by (A x t): Take the limit x 0, t 0

Heat transport – boundar conditions 1) Dirichlet boundary condition: temperature is given/controlled on the

Heat transport – boundar conditions 1) Dirichlet boundary condition: temperature is given/controlled on the boundary or part of it 2) System or part of it is closed to flux transfer (thermal isolation) 3) Newton’s law of cooling (the transfer of heat through boundary is proportional to the differnec between temperatures at both sides, of the boundary Here h is the heat transfer coefficient (W/m 2 K).

Chemical kinetics (introduction to formal descriptions of homogeneous reactions) Chemical reactions can be roughly

Chemical kinetics (introduction to formal descriptions of homogeneous reactions) Chemical reactions can be roughly divided in two broad categories: homogeneous and heterogeneous. Chemical reactions Homogeneous reactions Heterogeneous reactions Homogeneous reactions take place in the same phase in the whole volume. These reactions usually take place in a liquid or gaseous state. An example is the reaction between oxygen and nitrogen gases (without any catalyst) N 2 + O 2 2 NO Heterogeneous reactions usually take place at the boundaries of two different phases. In most cases one of the phase is in a solid state.

In a description of the chemical reaction we may be interested in how quickly

In a description of the chemical reaction we may be interested in how quickly it proceeds and how the concentration of reagents change with time. In an homogeneous reaction A+B C Is elementary, the we can assume the rate („speed”) of this reaction is given by v = k[A][B] Here we used a standard notation used in chemistry: [A] denotes the molar concentration per volume, i. e. [A] = n/V = (number of moles of A)/(volume) The SI unit of molar concentration is mole/m 3 but in laboratory practice the unit mole/dm 3 is used. The rate of reaction v means in this example that

The equation for [C] is an ordinary differential equation. But, we have to eliminate

The equation for [C] is an ordinary differential equation. But, we have to eliminate the functions [A] and [B ] to close the system. There are simple relations between [A], [B], and [C] coming from the fact that they are related by the transformation A + B C. If for simplicity we assume that at the beginning there was no product C in the system ([C]0 =0), then [A] = [A]0 – [C], [B] = [B]0 – [C]. where [A]0 = the initial (t=0) concentration for A, [B]0 = the initial (t=0) concentration for B. Thus we have ODE with initial condition problem where parameters [A]0, [B]0, and k are given.

Synthesis of hydrogen bromide It is a complex reaction with the balance formula The

Synthesis of hydrogen bromide It is a complex reaction with the balance formula The kinetics of this reaction (how the concentration of the product HBr changes with time) is given by the following ODE Sometimes written in equivalent form The kinetic constants k 1 and k 2 depend on temperature. In standard conditions k 2≈0, 1.

Synthesis of hydrogen bromide (contd. ) Let us denote y(t) = [HBr]. Taking the

Synthesis of hydrogen bromide (contd. ) Let us denote y(t) = [HBr]. Taking the overall reaction We have the molar balance equations Inserting these realtions int the equation fo d[HBr]/dt we obtain

Synthesis of hydrogen bromide (contd. ) The simulation of the reaction (solving the ODE

Synthesis of hydrogen bromide (contd. ) The simulation of the reaction (solving the ODE from the previous slide) True kinetics of HBr synthesis If the kinetics waw as for elementary reaction

Lotka-Volterra model (predator-prey system) In the lecture it was described in terms of the

Lotka-Volterra model (predator-prey system) In the lecture it was described in terms of the population dynamics context This is an overall reaction of converting A to B: A B. The presented mechanism of reaction leads to the following system of ODEs

The Lotka-Voltera system as a predator-prey model The are two interacting species in a

The Lotka-Voltera system as a predator-prey model The are two interacting species in a given region: the prey (for instance rabbits) and predators (for instance foxes). The prey have plenty of food and can easily multiply, but predators feed on the prey. Spatialy the system is homogeneous (the average density of both species is constant in the whole region). Neither species can leave nor enter the system from outside. Let us introduce: Thus we obtain the system of ODEs:

The Lotka-Voltera system – example of symulations In the calculations the following numbers were

The Lotka-Voltera system – example of symulations In the calculations the following numbers were assumed: