Kinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017

Similar documents
χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Mathematical Preparations

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

coordinates. Then, the position vectors are described by

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

In this section is given an overview of the common elasticity models.

Module 3: Element Properties Lecture 1: Natural Coordinates

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

Quantum Particle Motion in Physical Space

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Mechanics Physics 151

SYNChronised numerical methods

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Canonical transformations

PHYS 705: Classical Mechanics. Calculus of Variations II

MEV442 Introduction to Robotics Module 2. Dr. Santhakumar Mohan Assistant Professor Mechanical Engineering National Institute of Technology Calicut

Week 9 Chapter 10 Section 1-5

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

CONTROLLED FLOW SIMULATION USING SPH METHOD

Physics 181. Particle Systems

NUMERICAL DIFFERENTIATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Implicit Integration Henyey Method

Tensor Smooth Length for SPH Modelling of High Speed Impact

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Finite Element Modelling of truss/cable structures

Chapter 3. r r. Position, Velocity, and Acceleration Revisited

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

Affine and Riemannian Connections

Classical Field Theory

Rigid body simulation

PHYS 705: Classical Mechanics. Canonical Transformation II

11. Dynamics in Rotating Frames of Reference

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

Symmetric Lie Groups and Conservation Laws in Physics

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Solutions for Euler and Navier-Stokes Equations in Powers of Time

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD

Poisson brackets and canonical transformations

Arbitrary Lagrangian Eulerian Electromechanics in 3D

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

Lecture 20: Noether s Theorem

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Linear Approximation with Regularization and Moving Least Squares

Linear Momentum. Center of Mass.

Chapter 9: Statistical Inference and the Relationship between Two Variables

Solid Mechanics Z. Suo

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Constitutive Modelling of Superplastic AA-5083

CHAPTER 14 GENERAL PERTURBATION THEORY

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

Supplementary Notes for Chapter 9 Mixture Thermodynamics

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

6. Hamilton s Equations

Thermodynamics General

AP Physics 1 & 2 Summer Assignment

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Chapter 8. Potential Energy and Conservation of Energy

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

The Feynman path integral

MA 323 Geometric Modelling Course Notes: Day 13 Bezier Curves & Bernstein Polynomials

Lagrangian Field Theory

12. The Hamilton-Jacobi Equation Michael Fowler

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

Finite element method for structural dynamic and stability analyses

Introduction to elastic wave equation. Salam Alnabulsi University of Calgary Department of Mathematics and Statistics October 15,2012

THEOREMS OF QUANTUM MECHANICS

2 Finite difference basics

How Differential Equations Arise. Newton s Second Law of Motion

Inductance Calculation for Conductors of Arbitrary Shape

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Numerical Heat and Mass Transfer

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Phys304 Quantum Physics II (2005) Quantum Mechanics Summary. 2. This kind of behaviour can be described in the mathematical language of vectors:

6.3.7 Example with Runga Kutta 4 th order method

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Statistics Chapter 4

), it produces a response (output function g (x)

Linear Regression Analysis: Terminology and Notation

Geometry and Screw Theory for Robotics

Digital Signal Processing

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

Problem Points Score Total 100

= 1.23 m/s 2 [W] Required: t. Solution:!t = = 17 m/s [W]! m/s [W] (two extra digits carried) = 2.1 m/s [W]

Report on Image warping

6.3.4 Modified Euler s method of integration

Composite Hypotheses testing

PART 8. Partial Differential Equations PDEs

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

One Dimensional Axial Deformations

Transcription:

17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran way of analyzng flud moton. As the contnuum encompasses of several partcles, the Lagrangan analyss deals wth each partcles. To qucly descrbe certan quanttes, consder a contnuum n moton. Intally at tme t=to, the slope of the contnuum s as shown and s referred wth the orthogonal coordnates OX1XX3. Fg. 1: Representaton of the statement above (Source: Schaum s outlne of theory and problems of contnuum mechancs by George Mase) At a mathematcal pont P0, a partcle of the contnuum s assocated at tme t=t0. The poston vector of ths partcle s X X1I ˆ ˆ ˆ 1 X I X 3I3 The above poston vector s called materal coordnate. As the contnuum s movng, at a later tme t=t, the poston of the same partcle mght have changed. Even the contnuum also deforms shfts ts poston.

Let the poston at tme t be gven as n the fgure x x1eˆ ˆ ˆ 1 xe x3e3 Ths new poston descrpton s spatal coordnates. The same partcle that was present at P0 at tme t0 s now dsplaced and the dsplacement vector s gven as s. You can also see that the coordnates also shfted by a vector b. From vector algebra: s b x X If the coordnates OX1XX3 and ox1xx3 are merged, you get b 0 Hence, s x X (Ths means, x s the poston vector of the partcle at tme t, whose ntal poston s X ). In ndex notaton: s x X When the contnuum s n moton and deformaton, the partcles poston may be expressed n the form: x = x(x1,x,x3,t) or x x( X, t) You now, x present locaton of the partcle that occuped the pont (X1, X, X3) at tme t=to. (Ths s mappng the ntal confguraton wth the current confguraton). Such type of moton descrpton s Lagrangan formulaton. The dependent quantty s x and ndependent quantty s X. If the moton or deformaton s represented by the form: X = X(x1,x,x3,t) or X X ( x, t) where ndependent varable s x and t. Ths s Euleran formulaton. Ths descrpton provdes you the tracng of orgnal poston of the partcle that now occupes the spatal coordnate or locaton (x1, x, x3). The Lagrangan and Euleran mappngs are therefore nverse functons. For the nverse functons to exst, the necessary requrement s that the Jacoban must exst. Jacoban:

J x x x X X X 1 1 1 x x x x X X X X x x x X X X 3 3 3 If ths Jacoban (determnant) s zero, then unque nverse does not exst. From x = x(x1,x,x3,t), the Lagrangan form, you can form materal deformaton gradent by partally dfferentatng t wth X.e., x x x X X X x x x x X X X X x x x X X X F 1 1 1 3 3 3 Note: Ths s not a determnant. It s a tensor From X = X(x1,x,x3,t) The Euleran form, you get spatal deformaton gradent X X X x x x X X X X x x x x X X X x x x H 1 1 1 3 3 3 You can also form Materal Dsplacement Gradent. s x F X X (Recall s=x-x) As obvous, the materal dsplacement gradent s also a tensor. In smlar lnes, spatal dsplacement gradent tensor. can also be formulated as follows: s X H x x

The Deformaton Tensors To now about deformaton, the procedure s to see how much change s there between postons of two partcles from ther ntal confguraton (at t=t0) and later confguraton (at t=t). Consder the fgure below where the materal coordnates OX1XX3 and spatal coordnates ox1xx3 are merged. Fg. : The deformaton tensor representaton (Source: Schaum s outlne of theory and problems of contnuum mechancs by George Mase) There are two neghborng partcles that occupy postons P0 and Q0 ntally at tme t=t0. The dfferental elemental element length between two partcles s dx as per vector algebra. After a certan tme, at the nstant t=t, the contnuum has moved as well as deformed. The postons of those partcles are gven n spatal coordnates x and x dx The square of the dfferental element length between P0 and Q0 s: ( dx ). dx dx In ndex notaton, ( dx ). dx dx From X=X(x1,x,x3,t), you have seen:

X x dx ( ) F X x dx X dx dx dx x x X x C dx dx X X dx dx x Where C Cauchy s deformaton tensor. In the deformed confguratons, where the partcles are at postons P and Q, ( dx). dx dx ( dx) dxdx Also from Lagrangan expresson, x = x(x1,x,x3,t) x dx dx X ( ) dx dxdx dx dx X X ( dx) G dx dx x x Where G Green s deformaton tensor. The measure of deformaton s evaluated based on the dfference (dx) (dx) for the two neghborng partcles. ( ) ( ) dx dx dx dx dx dx X X x x dx dx X X L dx dx x x Here L Lagrangan or Green s fnte stran tensor L 1 x X x X

In a smlar way, you can form Euleran stran tensor E 1 X x X x Wth ths bref bacground nformaton on: 1. Materal coordnates OX1XX3. Spatal coordnates ox1xx3 x 3. Materal dervatve gradent X X 4. Spatal dervatve gradent x s x 5. Materal dsplacement gradent X X s X 6. Spatal dsplacement gradent x x X X 7. Cauchy s deformaton tensor x x x x 8. Green s deformaton tensor dxdx X X 9. Lagrangan s fnte stran tensor, etc. For fluds, we can descrbe propertes n Euleran or Lagrangan way..e., For example, the densty n the materal descrpton wll be: ρ = ρ(x1,x,x3,t).e. ρ = ρ(x, t) Ths wll be the densty of the flud partcle at the poston (X1,X,X3). In Euleran form : ρ= ρ(x(x1,x,x3,t), t) = ρ(x(x,t),t) = ρ * (x,t)