CFD_RIT.ppt
- 1. 1 Computational Fluid Dynamics Prof. (Dr.) S. S. Mohite
- 2. Computational Fluid Dynamics • Outline: – A brief review of Governing equations of fluid dynamics – Classification of Flows – An introductory overview of CFD 2
- 3. Fluid dynamics • Fluid dynamics is the science of fluid motion. • Fluid flow is commonly studied in one of three ways: – Experimental fluid dynamics. – Theoretical fluid dynamics. – Numerically: computational fluid dynamics (CFD). 3
- 4. 4
- 5. 5
- 6. 6
- 7. 7
- 8. 8
- 9. 9
- 10. 10
- 11. 11
- 12. 12
- 13. 13 Newton’s second law • For a solid mass: F = m.a • For a continuum: • Expressed in terms of velocity field u(x,y,z,t). In this form the momentum equation is also called Cauchy’s law of motion. • For an incompressible Newtonian fluid, this becomes: • Here p is the pressure and m is the dynamic viscosity. In this form, the momentum balance is also called the Navier-Stokes equation. f t . . u u u acceleration Mass per volume (density) Force per area (stress tensor) Force per volume (body force) u u u u 2 . m p t
- 14. 14 The Navier-Stokes equations for incompressible flow involve four basic quantities: – Local (unsteady) acceleration. – Convective acceleration. – Pressure gradients. – Viscous forces. • The ease with which solutions can be obtained and the complexity of the resulting flows often depend on which quantities are important for a given flow. (steady laminar flow) (impulsively started) (boundary layer) (inviscid, impulsively started) (inviscid) (unsteady flow) (steady viscous flow)
- 15. 15 Steady laminar flow • Steady viscous laminar flow in a horizontal pipe involves a balance between the pressure forces along the pipe and viscous forces. • The local acceleration is zero because the flow is steady. • The convective acceleration is zero because the velocity profiles are identical at any section along the pipe.
- 16. 16 Flow past an impulsively started flat plate • Flow past an impulsively started flat plate of infinite length involves a balance between the local (unsteady) acceleration effects and viscous forces. Here, the development of the velocity profile is shown. • The pressure is constant throughout the flow. • The convective acceleration is zero because the velocity does not change in the direction of the flow, although it does change with time.
- 17. 17 Boundary layer flow along a flat plate • Boundary layer flow along a finite flat plate involves a balance between viscous forces in the region near the plate and convective acceleration effects. • The boundary layer thickness grows in the downstream direction. • The local acceleration is zero because the flow is steady.
- 18. 18 Inviscid flow past an airfoil • Inviscid flow past an airfoil involves a balance between pressure gradients and convective acceleration. • Since the flow is steady, the local (unsteady) acceleration is zero. • Since the fluid is inviscid (m=0) there are no viscous forces.
- 19. 19 Impulsively started flow of an inviscid fluid • Impulsively started flow of an inviscid fluid in a pipe involves a balance between local (unsteady) acceleration effects and pressure differences. • The absence of viscous forces allows the fluid to slip along the pipe wall, producing a uniform velocity profile. • The convective acceleration is zero because the velocity does not vary in the direction of the flow. • The local (unsteady) acceleration is not zero since the fluid velocity at any point is a function of time.
- 20. 20 Steady viscous flow past a cylinder • Steady viscous flow past a circular cylinder involves a balance among convective acceleration, pressure gradients, and viscous forces. • For the parameters of this flow (density, viscosity, size, and speed), the steady boundary conditions (i.e. the cylinder is stationary) give steady flow throughout. • For other values of these parameters the flow may be unsteady.
- 21. 21 Unsteady flow past an airfoil • Unsteady flow past an airfoil at a large angle of attack (stalled) is governed by a balance among local acceleration, convective acceleration, pressure gradients and viscous forces. • A wide variety of fluid mechanics phenomena often occurs in situations such as these where all of the factors in the Navier- Stokes equations are relevant.
- 22. 22 Flow classifications • Laminar vs. turbulent flow. – Laminar flow: fluid particles move in smooth, layered fashion (no substantial mixing of fluid occurs). – Turbulent flow: fluid particles move in a chaotic, “tangled” fashion (significant mixing of fluid occurs). • Steady vs. unsteady flow. – Steady flow: flow properties at any given point in space are constant in time, e.g. p = p(x,y,z). – Unsteady flow: flow properties at any given point in space change with time, e.g. p = p(x,y,z,t).
- 23. 23 Newtonian vs. non-Newtonian Newtonian (low μ) Newtonian (high μ) Bingham-plastic 0 c Casson fluid Pseudo-plastic (shear-thinning) Dilatant (shear-thickening) Strain rate (1/s) (Pa) • Newtonian fluids: water, air. • Pseudoplastic fluids: paint, printing ink. • Dilatant fluids: dense slurries, wet cement. • Bingham fluids: toothpaste, clay. • Casson fluids: blood, yogurt. • Visco-elastic fluids: polymers (not shown in graph because viscosity is not isotropic).
- 24. 24 Flow classifications • Incompressible vs. compressible flow. – Incompressible flow: volume of a given fluid particle does not change. • Implies that density is constant everywhere. • Essentially valid for all liquid flows. – Compressible flow: volume of a given fluid particle can change with position. • Implies that density will vary throughout the flow field. • Compressible flows are further classified according to the value of the Mach number (M), where. • M < 1 - Subsonic. • M > 1 - Supersonic. c V M
- 25. 25 Flow classifications • Single phase vs. multiphase flow. – Single phase flow: fluid flows without phase change (either liquid or gas). – Multiphase flow: multiple phases are present in the flow field (e.g. liquid-gas, liquid-solid, gas-solid). • Homogeneous vs. heterogeneous flow. – Homogeneous flow: only one fluid material exists in the flow field. – Heterogeneous flow: multiple fluid/solid materials are present in the flow field (multi-species flows).
- 26. 26 Flow configurations - external flow • Fluid flows over an object in an unconfined domain. • Viscous effects are important only in the vicinity of the object. • Away from the object, the flow is essentially inviscid. • Examples: flows over aircraft, projectiles, ground vehicles.
- 27. 27 Flow configurations - internal flow • Fluid flow is confined by walls, partitions, and other boundaries. • Viscous effects extend across the entire domain. • Examples: flows in pipes, ducts, diffusers, enclosures, nozzles. airflow temperature profile car interior
- 28. 28 Classification of partial differential equations • A general partial differential equation in coordinates x and y: • Characterization depends on the roots of the higher order (here second order) terms: – (b2-4ac) > 0 then the equation is called hyperbolic. – (b2-4ac) = 0 then the equation is called parabolic. – (b2-4ac) < 0 then the equation is called elliptic. • Note: if a, b, and c themselves depend on x and y, the equations may be of different type, depending on the location in x-y space. In that case the equations are of mixed type. 0 2 2 2 2 2 g f y e x d y c y x b x a
- 29. 29 Origin of the terms • The origin of the terms “elliptic,” “parabolic,” or “hyperbolic” used to label these equations is simply a direct analogy with the case for conic sections. • The general equation for a conic section from analytic geometry is: where if. – (b2-4ac) > 0 the conic is a hyperbola. – (b2-4ac) = 0 the conic is a parabola. – (b2-4ac) < 0 the conic is an ellipse. 0 2 2 f ey dx cy bxy ax
- 30. 30 Elliptic problems • Elliptic equations are characteristic of equilibrium problems, this includes many (but not all) steady state flows. • Examples are potential flow, the steady state temperature distribution in a rod of solid material, and equilibrium stress distributions in solid objects under applied loads. • For potential flows the velocity is expressed in terms of a velocity potential: u=. Because the flow is incompressible, .u=0, which results in 2=0. This is also known as Laplace’s equation: • The solution depends solely on the boundary conditions. This is also known as a boundary value problem. • A disturbance in the interior of the solution affects the solution everywhere else. The disturbance signals travel in all directions. • As a result, solutions are always smooth, even when boundary conditions are discontinuous. This makes numerical solution easier! 0 2 2 2 2 y x
- 31. 31 Parabolic problems • Parabolic equations describe marching problems. This includes time dependent problems which involve significant amounts of dissipation. Examples are unsteady viscous flows and unsteady heat conduction. Steady viscous boundary layer flow is also parabolic (march along streamline, not in time). • An example is the transient temperature distribution in a cooling down rod: • The temperature depends on both the initial and boundary conditions. This is also called an initial-boundary-value problem. • Disturbances can only affect solutions at a later time. • Dissipation ensures that the solution is always smooth. 2 2 x t T=T0 T=T0 t=0 t x=0 x=L
- 32. 32 Hyperbolic problems • Hyperbolic equations are typical of marching problems with negligible dissipation. • An example is the wave equation: • This describes the transverse displacement of a string during small amplitude vibrations. If y is the displacement, x the coordinate along the string, and a the initial amplitude, the solution is: • Note that the amplitude is independent of time, i.e. there is no dissipation. • Hyperbolic problems can have discontinuous solutions. • Disturbances may affect only a limited region in space. This is called the zone of influence. Disturbances propagate at the wave speed c. • Local solutions may only depend on initial conditions in the domain of dependence. • Examples of flows governed by hyperbolic equations are shockwaves in transonic and supersonic flows. 2 2 2 2 2 x c t L x L ct a t x y sin cos ) , (
- 33. 33 Classification of fluid flow equations Steady Flow Unsteady Flow Viscous flow Elliptic Parabolic Inviscid flow M < 1 (subsonic) Elliptic Hyperbolic M> 1 (supersonic) Hyperbolic Hyperbolic Thin shear layers Parabolic Parabolic • For inviscid flows at M<1, pressure disturbances travel faster than the flow speed. If M>1, pressure disturbances can not travel upstream. Limited zone of influence is a characteristic of hyperbolic problems. • In thin shear layer flows, velocity derivatives in flow direction are much smaller than those in the cross flow direction. Equations then effectively contain only one (second order) term and become parabolic. Also the case for other strongly directional flows such as fully developed duct flow and jets.
- 34. 34 Example: the blunt-nosed body Bow Shock M > 1 M > 1 M < 1 Elliptic region Hyperbolic region δ Sonic Line Blunt-nosed body • Blunt-nosed body designs are used for supersonic and hypersonic speeds (e.g. Apollo capsules and spaceshuttle) because they are less susceptible to aerodynamic heating than sharp nosed bodies. • There is a strong, curved bow shock wave, detached from the nose by the shock detachment distance δ. • Calculating this flow field was a major challenge during the 1950s and 1960s because of the difficulties involved in solving for a flow field that is elliptic in one region and hyperbolic in others. • Today’s CFD solvers can routinely handle such problems, provided that the flow is calculated as being transient.
- 35. 35 Initial and boundary conditions • Initial conditions for unsteady flows: – Everywhere in the solution region , u and T must be given at time t=0. • Typical boundary conditions for unsteady and steady flows: – On solid walls: • u = uwall (no-slip condition). • T=Twall (fixed temperature) or kT/n=-qwall (fixed heat flux). – On fluid boundaries. • For most flows, inlet: , u, and T must be known as a function of position. For external flows (flows around objects) and inviscid subsonic flows, inlet boundary conditions may be different. • For most flows, outlet: -p+m un/n=Fn and m ut/n=Fn (stress continuity). F is the given surface stress. For fully developed flow Fn=-p (un/n =0) and Ft=0. For inviscid supersonic flows, outlet conditions may be different. • A more detailed discussion of boundary conditions will be given later on.
- 36. 36 Summary • Fluid flows can be classified in a variety of ways: – Internal vs. external. – Laminar vs. turbulent. – Compressible vs. incompressible. – Steady vs. unsteady. – Supersonic vs. transonic vs. subsonic. – Single-phase vs. multiphase. – Elliptic vs. parabolic vs. hyperbolic.
- 37. What is computational fluid dynamics? • Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena by solving the mathematical equations which govern these processes using a numerical Methods. • The result of CFD analyses is relevant engineering data used in: – Conceptual studies of new designs. – Detailed product development. – Troubleshooting. – Redesign. • CFD analysis complements testing and experimentation. • Reduces the total effort required in the laboratory. 37
- 38. 38
- 39. 39
- 40. 40
- 41. 41
- 42. 42
- 43. 43
- 44. 44
- 45. 45
- 46. 46
- 47. 47
- 48. 48
- 49. 49
- 50. 50
- 51. 51
- 52. 52
- 53. 53
- 54. 54
- 55. 55
- 56. 56
- 57. 57
- 58. 58
- 59. 59
- 60. 60