1 / 48

HE 316 Term Project Presentation Symmetry Analysis in Fluid Dynamics

HE 316 Term Project Presentation Symmetry Analysis in Fluid Dynamics. Saikishan Suryanarayanan Engineering Mechanics Unit JNCASR. Outline. Introduction to Symmetry Analysis Lie series, Group operator and infinitesimal invariance condition for functions.

manjit
Download Presentation

HE 316 Term Project Presentation Symmetry Analysis in Fluid Dynamics

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. HE 316 Term Project PresentationSymmetry Analysis in Fluid Dynamics Saikishan Suryanarayanan Engineering Mechanics Unit JNCASR

  2. Outline • Introduction to Symmetry Analysis • Lie series, Group operator and infinitesimal invariance condition for functions. • Applications to ODEs and PDEs – general remarks • Boundary Layer • Using invariance group to identify self-similarity and solve Blasius boundary layer • Incompressible Navier Stokes Equation. • Group of invariance • Application : Viscous flow in a diverging channel (Jeffery-Hamel fl0w) • Compressible Euler Equation • Group of invariance

  3. Introduction to Symmetry Analysis

  4. Introduction to Symmetry analysis • Symmetry analysis involves identifying and studying the group of symmetries leaves a function/ equation invariant. • Most physical laws have associated symmetry groups, for eg. the Schrödinger equation is invariant under rotation. • When a differential equation and Boundary conditions are invariant under a certain group, one expects the solution also be invariant under the same symmetry group. • Identifying the symmetry groups can provide physical insight and help solve practical problems, • ODE - Order reduction • PDE - Reduction in # of dimensions, some cases reduction to a ODE. • Better numerical schemes can be developed for solving ODE/PDE.

  5. Lie Groups, Lie Series and Condition for Invariance Symmetry groups of differential equations are continuous groups. Illustration of 1 parameter Lie group – The transformation x -> x’ is defined as withSince the group is continuous wrt. s, the transformation can be expressed as Infinitesimal of the group is defined asConsider a function y [x]. The function transforms as where

  6. Lie Groups, Lie Series and Condition for Invariance The function is invariant , ie. y [x’] = y [x] if and only if Therefore a group can be expressed in form of its infinitesimals. Another common way of expressing a Lie Group is by Transformation of the function can be expressed asThe number of operators (Xs) is equal to the number of parameters of the Lie group. The condition for invariance of a differential equation is similar to the condition for invariance of functions, except that it is applied to the extended tangent space of differential functions.

  7. Laminar Boundary Layer Derivation of self – similar solutions (generally obtained by physical insight / trial and error) using Symmetry Group Analysis.

  8. 2D Navier Stokes Equations • The Continuity (conservation of mass) and Navier – Stokes Equations (conservation of linear momentum) in two dimensions are

  9. The Boundary Layer • In a seminal paper in 1904, Prandtl proposed that the changes in velocity from the free stream to the velocity at the surface take place in a thin layer near the wall termed as the ‘boundary layer’. . Figure from : http://www-scf.usc.edu/~tchklovs/Proposal.htm

  10. Boundary Layer Equations • In the Boundary layer, the changes in the stream-wise (x) direction are small compared to changes in the wall normal (y) direction. Consequently (from continuity) v << u. • Under this ansatz and under steady state and constant free stream velocity (zero pressure gradient) the governing equations are simplified to

  11. Stream Function Formulation • Introducing the stream function y , defined as • The Boundary Layer equation and boundary conditions in terms of stream function are : • This is a third order PDE with 2 independent variables.

  12. Symmetry group of the governing equation Consider the dilatation group, The governing equation transforms as If a= b – c, the original form is recovered, therefore the governing eqn is invariant under the two parameter dilation group :

  13. Symmetry group of the Boundary condition • The inhomogeneous boundary conditions transform as • In order to leave the b.c. invariant, there is an additional constraint on the group parameters, b = 2c • The equation, the boundary condition are invariant under the following transformation • So it is expected that the solution will also be invariant under this lie group.

  14. Dimensional Analysis and Reduction to ODE • The stream function can be expressed as a function of the co-ordinates and problem parameters as : • A naïve dimensional analysis (using Buckingham Pi theorem) leads to : • On the transformation by the symmetry group of the governing PDE and BC

  15. Dimensional Analysis and Reduction to ODE • y/x is not invariant to the group transformation, therefore the solution, should be of the form • On substitution into the original PDE, it leads to the following ODEand Boundary conditions

  16. Comments, implications and solution • Solution is Self- Similar – the velocity profile only gets stretched from one stream-wise location to another. • Without using group theory, this problem is solved by physical insight/ experimental observation of self similar profiles/ guess work. • Numerical Solution to the ODE. • u/Ue = f [y/d] d ~ x / (Rex)1/2

  17. Symmetry group of the incompressible 3D Navier-Stokes Equations

  18. Incompressible Navier Stokes Equations • For incompressible flow, the conservation equations of mass and momentum are given by • Note that the above form applies to both two and three dimensions. Einstein convention is used, dummy indices are summed over. • The pressure satisfies a Possion equation obtained by taking divergence of the momentum equation.

  19. Time Translation • The Equations are invariant under the following transformation • This can be shown as • The infinitesimal of this group is given by • Also holds for compressible flows (shown later)

  20. Arbitrary function of time added to pressure • The Equations are invariant under the following transformation • This can be shown as • In operator form, • Doesn’t apply for compressible flows.

  21. Rotation/ Reflection • The Navier stokes equations are invariant under uniform rotation/reflection of the position and velocity vectors where In Cartesian indicial notation, since

  22. Rotation • The invariance is demonstrated as follows :

  23. Rotation/Reflection In operator form, Rotation about z- axis Rotation about x- axis Rotation about y-axisApplies for compressible flows also.

  24. Non-uniform Translation • The Incompressible Navier Stokes equations are invariant not only under constant translation, but also translations of the form: • This is demonstrated as follows :

  25. Non – uniform translation • In operator form, the non-uniform translation about x, y and z axes are: • Physically, invariance under non-uniform translation implies that the incompressible NS equations are valid even in accelerating frames of reference, provided the acceleration is irrotational. • This invariance does not hold for compressible flows, unless d2a/dt2 = 0

  26. Dilatation • The full NS equations are invariant under the following one –parameter dilatation group (note: the boundary layer equations were invariant under a two parameter dilatation group) • The partial derivatives transform as

  27. Dilatation • The continuity equation is invariant under the dilatation group transformation as shown below • The invariance of the momentum eqn is shown as follows:

  28. Dilatation • In operator form, • Aside : Euler Equations (when n = 0) admit a two parameter dilatation group.

  29. Commutator Table

  30. Jeffery-Hamel Flow

  31. Jeffery-Hamel Flow • 2D viscous flow through converging/ diverging channel.

  32. Jeffery-Hamel Flow The governing equation for the stream function (2D INS) is which is invariant under the dilatation group as shown in the last section Dimensional analysis of and choosing dimensionless parameters that are invariant to the group transformations leads to

  33. In polar co-ordinates, • On substitution into the original PDE, the following ODE is obtained

  34. Solution • The ODE can be solved numerically, to obtain the following self-similar solution. http://www.cfd-online.com/Wiki/Jeffery-Hamel_flow

  35. Symmetry Group of the Compressible Euler Equations

  36. Compressible Euler Equations • For compressible flow of an inviscid fluid, the conservation equations of mass, momentum and energy are given by

  37. Time Translation

  38. Rotations

  39. Translation

  40. Dilatation • Since all terms in all equations have either a single time derivative or a single space derivative, it can be seen that they are invariant under the above dilatation group – not valid if viscosity is present. • Equations admit additional symmetry groups based on functional form of F[p,r]

  41. References • Brian J. Cantwell. 2002. Introduction to Symmetry Analysis . Cambridge University Press .

More Related