1 / 35

Numerical Hydraulics

Numerical Hydraulics. Lecture 1: The equations. Wolfgang Kinzelbach with Marc Wolf and Cornel Beffa. Contents of course. The equations Compressible flow in pipes Numerical treatment of the pressure surge Flow in open channels Numerical solution of the St. Venant equations Waves.

noura
Download Presentation

Numerical Hydraulics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Numerical Hydraulics Lecture 1: The equations Wolfgang Kinzelbach with Marc Wolf and Cornel Beffa

  2. Contents of course • The equations • Compressible flow in pipes • Numerical treatment of the pressure surge • Flow in open channels • Numerical solution of the St. Venant equations • Waves

  3. Basic equations of hydromechanics • The basic equations are transport equations for • Mass, momentum, energy … • General treatment • Transported extensive quantity m • Corresponding intensive quantity f (m/Volume) • Flux j of quantity m • Volume-sources/sinks s of quantity m

  4. Extensive/intensive quantities • Extensive quantities are additive • e.g. volume, mass, energy • Intensive quantities are specific quantities, they are not additive • e.g. temperature, density • Integration of an intensive quantity over a volume yields the extensive quantity

  5. Balance over a control volume unit normal to surface boundary G flux volume W Balance of quantity m: minus sign, as orientation of normal to surface and flux are in opposite direction

  6. Differential form • Using the Gauss integral theorem we obtain: The basic equations of hydromechanics follow from this equation for special choices of m, f, s and j

  7. Continuity equation • m = M (Mass), f = r (Density), j = ur (Mass flux) yields the continuity equation for the mass: For incompressible fluids (r = const.) we get: For compressible fluids an equation of state is required:

  8. Other approach: General principle: in 1D Time interval [t, t+Dt] Storage is change in extensive quantity Cross-sectional area A Volume V = ADx Fluxin Fluxout Gain/loss from volume sources/sinks Dx x+Dx x Conservation law in words:

  9. General principle in 1D Division by DtDxA yields: In the limit Dt, Dx to 0:

  10. General principle in 3D or

  11. Mass balance: in 1D Density assumed constant! Storage can be seen as change in intensive quantity Time interval [t, t+Dt] V=ADx Dx x+Dx x Conservation equation for water volume

  12. Mass balance: in 1D continued In the limit

  13. Generalization to 3D or

  14. Essential derivative The total or essential derivative of a time-varying field quantity is defined by The total derivative is the derivative along the trajectory given by the velocity vector field Using the total derivative the continuity equation can be written in a different way

  15. Momentum equation (equation of motion) • Example: momentum in x-direction • m = Mux (x-momentum), f = rux (density), (momentum flux), sx force density (volume- and surface forces) in x-direction inserted into the balance equation yields the x-component of the Navier-Stokes equations: pressure force gravity force friction force per unit volume In a rotating coordinate system the Coriolis-force has to be taken into account

  16. Momentum equation (equation of motion) • Using the essential derivative and the continuity equation we obtain: • The x-component of the pressure force per unit volume is • The x-component of gravity per unit volume is • The friction force per unit volume will be derived later Newton: Ma = F

  17. Momentum equation (equation of motion) • In analogy to the x-component the equations for the y- and z-component can be derived. Together they yield a vector equation:

  18. Momentum equation (equation of motion) • Writing out the essential derivative we get: • The friction term fR depends on the rate of deformation. The relation between the two is given by a material law.

  19. Friction force

  20. Friction force • The strain forms a tensor of 2nd rank The normal strain only concerns the deviations from the mean pressure p due to friction: deviatoric stress tensor. The tensor is symmetric. • The friction force per unit volume is

  21. The material law • Water is in a very good approximation a Newtonian fluid: strain tensor a tensor of deformation • Deformations comprise shear, rotation and compression

  22. y y y x x x rotation Deformation shearing compression

  23. Compression Relative volume change per time

  24. Shearing and rotation Dy Dx

  25. Shearing and rotation The shear rate is The angular velocity of rotation is

  26. General tensor of deformation rotation and shear components Symmetric part (shear velocity) contains the friction Anti-symmetric part (angular velocity of rotation) frictionless x,y,z represented by xi with i=1,2,3

  27. Material law according to NewtonMost general version with Three assumptions: Stress tensor is a linear function of the strain rates The fluid is isotropic For a fluid at rest must be zero so that hydrostatic pressure results h is the usual (first) viscosity, l is called second viscosity

  28. Resulting friction term for momentum equation Friction force on volume element Compression force due to friction It can be shown that If one assumes that during pure compression the entropy of a fluid does not increase (no dissipation).

  29. Navier-Stokes equations Under isothermal conditions (T = const.) one has thus together with the continuity equation 4 equations for the 4 unknown functions ux, uy, uz, and p in space and time. They are completed by the equation of state for r(p) as well as initial and boundary conditions.

  30. Vorticity • The vorticity is defined as the rotation of the velocity field

  31. Vorticity equation • Applying the operator to the Navier-Stokes equation and using various vector algebraic identities one obtains in the case of the incompressible fluid: • The Navier-Stokes equation is therefore also a transport equation (advection-diffusion equation) for vorticity. • Other approach: transport equation for angular momentum

  32. Vorticity equation • Pressure and gravity do not influence the vorticity as they act through the center of mass of the mass particles. • Under varying density a source term for vorticity has to be added which acts if the gravitational acceleration is not perpendicular to the surfaces of equal pressure (isobars). • In a rotating reference system another source term for the vorticity has to be added.

  33. Energy equation • m = E, f = r(e+u²/2) inner+kinetic energy per unit volume, j = fu=r(e+u²/2)u, s work done on the control volume by volume and surface forces, dissipation by heat conduction

  34. Energy equation • The new variable e requires a new material equation. It follows from the equation of state: e = e(T,p) • In the energy equation, additional terms can appear, representing adsorption of heat radiation

  35. Solute transport equation • m = Msolute, f = c concentration, (advection and diffusion), s solute sources and sinks Advection-diffusion equation for passive scalar transport in microscopic view.

More Related