Quaternions and particle dynamics in the Euler fluid equations. J. D. Gibbon. Mathematics Department Imperial College London, London SW7 2AZ
|
|
- Mabel Anthony
- 5 years ago
- Views:
Transcription
1 Quaternions and particle dynamics in the Euler fluid equations J. D. Gibbon Mathematics Department Imperial College London, London SW7 2AZ Maryland CSCAMM Meeting; October 23-27, 2006 J. D. Gibbon, CSCAMM; October
2 Collaborators JDG & Darryl Holm 06: JDG, Holm, Kerr & Roulstone: Nonlinearity 19, , 2006 JDG, Physica D, 166, 17-28, Galanti, JDG & Heritage; Nonlinearity 10, 1675, Galanti, JDG & Kerr, in Turbulence structure & vortex dynamics, (pp 23-34, eds: Hunt & Vassilicos, CUP 2000). J. D. Gibbon, CSCAMM; October
3 Summary of this talk Question: Do the Euler equations possess some subtle geometric structure that guides the direction of vorticity see Peter Constantin, Geometric statistics in turbulence, SIAM Reviews, 36, Quaternions: what are they? 2. Lagrangian particle dynamics: We find explicit equations for the Lagrangian derivatives of an ortho-normal co-ordinate system at each point in space. (JDG/Holm 06) 3. For the 3D-Euler equations; Ertel s Theorem shows how Euler fits naturally into this framework (JDG, Holm, Kerr & Roulstone 2006). 4. Review of work on the direction of Euler vorticity, particularly that of Constantin, Fefferman & Majda 1996; Deng, Hou & Yu 2005/6 & Chae A different direction of vorticity result involving the pressure Hessian. J. D. Gibbon, CSCAMM; October
4 Lord Kelvin (William Thompson) once said: Quaternions came from Hamilton after his best work had been done, & though beautifully ingenious, they have been an unmixed evil to those who have touched them in any way. O Connor, J. J. & Robertson, E. F Sir William Rowan Hamilton, history/mathematicians/hamilton.html Kelvin was wrong because quaternions are now used in the computer animation, avionics & robotics industries to track objects undergoing sequences of tumbling rotations. Visualizing quaternions, by Andrew J. Hanson, MK-Elsevier, Quaternions & rotation Sequences: a Primer with Applications to Orbits, Aerospace & Virtual Reality, J. B. Kuipers, Princeton University Press, J. D. Gibbon, CSCAMM; October
5 What are quaternions? (Hamilton (1843)) Imperial College London Quaternions are constructed from a scalar p & a 3-vector q by forming the tetrad p = [p, q] = pi q σ, q σ = 3 q i σ i i=1 based on the Pauli spin matrices that obey the relations σ i σ j = δ ij ɛ ijk σ k ( ) ( ) ( ) 0 i 0 1 i 0 σ 1 = σ 2 = σ 3 =. i i Thus quaternions obey the multiplication rule p 1 p 2 = [p 1 p 2 q 1 q 2, p 1 q 2 + p 2 q 1 + q 1 q 2 ]. They are associative but obviously non-commutative. J. D. Gibbon, CSCAMM; October
6 Quaternions, Rotations and Cayley-Klein parameters Let ˆp = [p, q] be a unit quaternion with inverse ˆp = [p, q] with p 2 + q 2 = 1, which guarantees ˆp ˆp = [1, 0]. For a pure quaternion r = [0, r] there exists a transformation r r r = ˆp r ˆp = [0, (p 2 q 2 )r + 2p(q r) + 2q(r q)]. Now choose p = ± cos 1 2θ and q = ± ˆn sin 1 2θ, where ˆn is the unit normal to r r = ˆp r ˆp = [0, r cos θ + ( ˆn r) sin θ], where ˆp = ±[cos 1 2θ, ˆn sin 1 2θ]. This represents a rotation by an angle θ of the 3-vector r about its normal ˆn. The elements of the unit quaternion ˆp are the Cayley-Klein parameters from which the Euler angles can be calculated. All terms are quadratic in p and q, and thus allow a double covering (±) (see Whittaker 1945). J. D. Gibbon, CSCAMM; October
7 General Lagrangian evolution equations Imperial College London Consider the general Lagrangian evolution equation for a 3-vector w such that Dw D = a(x, t) = t + u transported by a velocity field u. Define the scalar α a the 3-vector χ a as α a = w 1 (ŵ a), χ a = w 1 (ŵ a). for w = 0. Via the decomposition a = α a w + χ a w, w & ŵ satisfy D w Dŵ = α a w, = χ a ŵ. α a is the growth rate (Constantin 1994) & χ a is the swing rate. The tetrads q a = [α a, χ a ], w = [0, w]. allow us to write this as Dw = q a w. J. D. Gibbon, CSCAMM; October
8 Theorem: (JDG/Holm 06) If a is differentiable in the Lagrangian sense s.t. Da = b(x, t), (i) For for w = 0, q a and q b satisfy the Ricatti equation Dq a + q a q a = q b ; (ii) At each point x there exists an ortho-normal frame (ŵ, ˆχ a, ŵ ˆχ a ) SO(3) whose Lagrangian time derivative is expressed as Dŵ = D ab ŵ, D(ŵ ˆχ a ) = D ab (ŵ ˆχ a ), D ˆχ a = D ab ˆχ a, where the Darboux angular velocity vector D ab is defined as D ab = χ a + c 1 ŵ, c 1 = ŵ ( ˆχ χ a χ b ). a J. D. Gibbon, CSCAMM; October
9 Lagrangian frame dynamics: tracking a particle Imperial College London ˆχ a ŵ ŵ ˆχ a ŵ ŵ ˆχ a ˆχ a t 1 t 2 The dotted line represents a particle ( ) trajectory moving from (x 1, t 1 ) to (x 2, t 2 ). The orientation of the orthonormal unit vectors {ŵ, ˆχ a, (ŵ ˆχ a )} is driven by the Darboux vector D ab = χ a + c 1 χ a ŵ where c 1 = ŵ ( ˆχ a χ b ). Thus we need the quartet of vectors to make this process work {u, w, a, b}. J. D. Gibbon, CSCAMM; October
10 Proof : (i) It is clear that (with q b = [α b, χ b ]) D 2 w = [0, b] = q 2 b w. Compatibility between this and the q-equation means that ( ) Dqa + q a q a q b w = 0, (ii) Now consider the ortho-normal frame (ŵ, ˆχ a, ŵ ˆχ a ) as in the Figure below. The evolution of χ a comes from ŵ. a ˆχa b...ŵ... ˆχ a Dq a + q a q a = q b, and gives Dχ a = 2α a χ a + χ b. J. D. Gibbon, CSCAMM; October
11 b can be expressed in this ortho-normal frame as the linear combination b = w [α b ŵ + c 1 ˆχ a + c 2 (ŵ ˆχ a )], χ b = c 1 (ŵ ˆχ a ) c 2 ˆχ a, where c 1 = ŵ ( ˆχ a χ b ) and c 2 = ( ˆχ a χ b ). From the Ricatti equation for the tetrad q a = [α a, χ a ] (where χ a = χ a ) Dχ a = 2α a χ a + χ b, Dχ a = 2α aχ a c 2, There follows D ˆχ a = c 1 χ 1 D(ŵ ˆχ a (ŵ ˆχ a ), a ) = χ a ŵ c 1 χ 1 a ˆχ a, which, together with Dŵ a = χ a ŵ, can be re-expressed in terms of the Darboux vector D a = χ a + c 1 χ a ŵ. J. D. Gibbon, CSCAMM; October
12 Ertel s Theorem & the 3D Euler equations Imperial College London Dω = ω u = Sω Euler in vorticity format Theorem: (Ertel 1942) If ω satisfies the 3D incompressible Euler equations then any arbitrary differentiable μ satisfies ( ) [ ] D Dμ D (ω μ) = ω =, ω = 0. Proof: Consider ω μ ω i μ,i D (ω i μ,i ) = Dω { ( ) } i Dμ μ,i + ω i u k,i μ, k x i ( ) Dμ = {ω j u i,j μ,i ω i u k,i μ, k } + ω }{{} i x i zero under summation In characteristic (Lie-derivative) form, ω x (t) = ω invariant (Cauchy 1859) and is frozen in. x (0) is a Lagrangian J. D. Gibbon, CSCAMM; October
13 Various references Ertel; Ein Neuer Hydrodynamischer Wirbelsatz, Met. Z. 59, , (1942). Truesdell & Toupin, Classical Field Theories, Encyclopaedia of Physics III/1, ed. S. Flugge, Springer (1960). Ohkitani; Phys. Fluids, A5, 2576, (1993). Kuznetsov & Zakharov; Hamiltonian formalism for nonlinear waves, Physics Uspekhi, 40 (11), (1997). Bauer s thesis 2000 (ETH-Berlin); Gradient entropy vorticity, potential vorticity and its history. Viudez; On the relation between Beltrami s material vorticity and Rossby- Ertel s Potential, J. Atmos. Sci. (2001). J. D. Gibbon, CSCAMM; October
14 Ohkitani s result & the pressure Hessian Imperial College London Define the Hessian matrix of the pressure { 2 } p P = {p,ij } = x i x j then Ohkitani took μ = u i (Phys. Fluids, A5, 2576, 1993). Result: The vortex stretching vector ω u = Sω obeys D(ω u) = D(Sω) ( ) Du = ω = P ω Thus for Euler, via Ertel s Theorem, we have the identification: w ω a ω u = Sω b P ω with a quartet (u, w, a, b) (u, ω, Sω, P ω). J. D. Gibbon, CSCAMM; October
15 Euler: the variables α(x, t) and χ(x, t) Imperial College London ˆω S ˆω ˆχ ˆω ˆχ S ˆω = α ˆω + χ ˆω See JDG, Holm, Kerr & Roulstone (α a ) α = ˆω S ˆω χ = ˆω S ˆω (χ a ) ( α b ) α p = ˆω P ˆω χ p = ˆω P ˆω ( χ b ) q = [α, χ] q b = q p = [α p, χ p ] Dq + q q + q p = 0, constrained by T rp = Δp = u i,j u j,i = 1 2ω 2 T rs 2. J. D. Gibbon, CSCAMM; October
16 Lagrangian frame dynamics: tracking an Euler fluid particle ˆχ ˆω ˆω ˆχ ˆω ˆω ˆχ ˆχ t 1 t 2 The dotted line represents the fluid packet ( ) trajectory moving from (x 1, t 1 ) to (x 2, t 2 ). The orientation of the orthonormal unit vectors {ˆω, ˆχ, (ˆω ˆχ)} is driven by the Darboux vector D = χ + c 1 χ ˆω, c 1 = ˆω ( ˆχ χ p ). Thus the pressure Hessian within c 1 drives the Darboux vector D. J. D. Gibbon, CSCAMM; October
17 The α and χ equations In terms of α and χ, the Ricatti equation for q Dq a + q a q a = q b ; becomes Dα = χ2 α 2 Dχ α p, = 2αχ χ p. (Galanti, JDG & Heritage; Nonlinearity 10, 1675, 1997). Stationary values are α = γ 0, χ = 0, α p = γ0 2 which correspond to Burgers -like vortices. When tubes & sheets bend & tangle then χ 0 and q becomes a full tetrad driven by q p which is coupled back through the elliptic pressure condition. Note: Off-diagonal elements of P change rapidly near intense vortical regions across which χ p and α p change rapidly. J. D. Gibbon, CSCAMM; October
18 Phase plane On Lagrangian trajectories, the α χ equations become α t = χ2 α 2 α p, where C p = ˆχ χ p. χ t = 2αχ + C p. In regions of the α χ phase plane where α p = const, C p = const there are 2 critical points: (α, χ) = (±α 0, χ 0 ) 2α0 2 = α p + [αp 2 + Cp] 2 1/2 The critical point in the LH-half-plane ( α 0, χ 0 ) is an unstable spiral; The critical point in the RH-half-plane is (α 0, χ 0 ) is a stable spiral. J. D. Gibbon, CSCAMM; October
19 The next few slides: remarks on the direction of vorticity in Euler 1. The BKM theorem 2. The work of Constantin, Fefferman & Majda 1996 and Constantin 1994 The work of Deng, Hou & Yu 2005/6 Can our quaternionic Ricatti equation give anything in terms of P? J. D. Gibbon, CSCAMM; October
20 The Beale-Kato-Majda Theorem (CMP 94, 61-6, 1984) Theorem: There exists a global solution of the Euler equations u C([0, ]; H s ) C 1 ([0, ]; H s 1 ) for s 3 if, for every t > 0, t 0 ω(τ) L (Ω) dτ <. The proof is based on u c ω [1 + log H 3 ]. Thus one needs to numerically monitor only t 0 ω(τ) dτ. Corollary: If a singularity is observed in a numerical experiment of the form ω (t t) β then β must lie in the range β 1 for the singularity to be genuine & not an artefact of the numerical calculation. J. D. Gibbon, CSCAMM; October
21 Constantin, Fefferman & Majda; Comm PDEs, 21, , 1996 The image W t of a set W 0 is given by W t = X(t, W 0 ). W 0 is said to be smoothly directed if there exists a length ρ > 0 and a ball 0 < r < 1 2ρ such that: 1. ˆω(, t) has a Lipschitz extn to the ball of radius 4ρ centred at X(q, t) & M = lim sup t T q W 0 t 0 ˆω(, t) 2 L (B 4ρ ) dt <. i.e. the direction of vorticity is well-behaved in the nbhd of a set of trajectories. 2. The condition sup B3r (W t ) ω(x, t) m sup Br (W t ) ω(x, t) holds for all t [0, T ) with m = const > 0 ; i.e. this nbhd captures large & growing vorticity but not so that it overlaps with another similar region & sup B4r (W t ) u(x, t) U(t) := sup x u(x, t) < (Cordoba & Fefferman 2001; for tubes). Theorem: (CFM 1996) Assume that W 0 is smoothly directed as in (i) (ii). Then a time τ > 0 & a constant Γ s.t. for any 0 t 0 < T and 0 t t 0 τ sup B r (W t ) ω(x, t) Γ sup B ρ (W t ) ω(x, t 0 ). J. D. Gibbon, CSCAMM; October
22 The work of Deng, Hou & Yu; Comm PDEs, 31, , 2006 Consider a family of vortex line segments L t in a region of max-vorticity. Denote by L(t) the arc length of L t, ˆn the unit normal & κ the curvature. DHY define Uˆω (t) max (u ˆω) (x, t) (u ˆω) (y, t), x,y L t U n (t) max Lt u ˆn, and M(t) max ( ˆω L (L t ), κ L (L t )). Theorem: (Deng, Hou & Yu 06) : Let A, B (0, 1) with B = 1 A, and C 0 be a positive constant. If 1. Uˆω (t) + U n (t) (T t) A, 2. M(t)L(t) C 0, 3. L(t) (T t) B, then there will be no blow-up up to time T. Also J. Deng, T. Y. Hou & X. Yu; Comm. PDEs, 30, , J. D. Gibbon, CSCAMM; October
23 Using the pressure Hessian Imperial College London (see also Chae: T 0 S ˆω P ˆω dτ < ; Comm. P&A-M., 109, 1 21, 2006). Theorem: (JDG, Holm, Kerr & Roulstone 06) : a global solution of the Euler equations, u C([0, ]; H s ) C 1 ([0, ]; H s 1 ) for s 3 if T 0 χ p L (D) dτ <, with the exception of when ˆω becomes collinear with an e-vec of P at t = T. Proof: With Sŵ 2 = α 2 + χ 2, D Sŵ α Sŵ + α α p + χ χ p. (α 2 + χ 2 ) 1/2 Because D ω / = α ω, our concern is with α 0 D Sŵ α p + χ p. Possible that P ˆω blows up simultaneously as the angle between ˆω and P ˆω 0 thus keeping χ p finite; i.e. t 0 χ p L (D)dτ < but t 0 α p L (D)dτ. J. D. Gibbon, CSCAMM; October
24 Frame dynamics & the Frenet-Serret equations Imperial College London With ŵ as the unit tangent vector, ˆχ as the unit bi-normal and ŵ ˆχ as the unit principal normal, the matrix N can be formed N = ( ŵ T, (ŵ ˆχ) T, ˆχ T), with 0 χ a 0 DN = GN, G = χ a 0 c 1 χ 1 a. 0 c 1 χ 1 a 0 The Frenet-Serret equations for a space-curve are 0 κ 0 dn ds = F N where F = κ 0 τ, 0 τ 0 where κ is the curvature and τ is the torsion. J. D. Gibbon, CSCAMM; October
25 The arc-length derivative d/ds is defined by d ds = ˆω. The evolution of the curvature κ and torsion τ may be obtained from Ertel s theorem expressed as the commutation of operators [ D, ω ] = 0 d [ D α a ds +, d ] = 0. ds This commutation relation immediately gives α a F + DF = dg ds + [G, F ]. Thus Ertel s Theorem gives explicit evolution equations for the curvature κ and torsion τ that lie within the matrix F and relates them to c 1, χ a and α a. J. D. Gibbon, CSCAMM; October
26 Mixing Consider a passive vector line-element δl in a flow transported by an independent velocity field u. For small δl we have the same equations as Euler for ω Dδl D = δl u = t + u Following the analogy with Euler, Ertel s Theorem holds so there is a b-field: ( ) D(δl u) Du = δl Du/ represents any dynamics one wishes to impose on the problem. Thus all the conditions hold for Theorem 1: 1. w = δl 2. a = δl u 3. b = δl ( ) Du a Ricatti equation plus an ortho-normal frame... J. D. Gibbon, CSCAMM; October
27 Ideal MHD Consider a magnetic field B coupled to a fluid (div u = 0 = div B) Du DB = B B p = B u Defining Elsasser variables with ±-material derivatives (two time-clocks) v ± = u ± B ; D ± = t + v± the magnetic field B and v ± satisfy with div v ± = 0 D ± v = p ; D ± B = B v ± Moffatt (1985) suggested that B takes the place of ω in ideal MHD. J. D. Gibbon, CSCAMM; October
28 Ertel s Theorem (proof omitted) for this system is D (B v ± ) = P B. With two time-clocks, we have the correspondence w B a ± B v ± b P B α pb = ˆB P ˆB Define tetrads q ± and q pb as follows χ pb = ˆB P ˆB q ± = [ α ±, χ ±] q pb = [ ] α pb, χ pb. The tetrads q ± satisfy the compatibility relation D q ± + q ± q + q pb = 0 J. D. Gibbon, CSCAMM; October
29 MHD-Lagrangian frame dynamics Imperial College London We have 2 sets of orthonormal vectors opposite Lagrangian time derivatives. ˆB, ( ˆB ˆχ ± ), ˆχ ± acted on by their D ˆB D ( ˆB ˆχ ± ) D ˆχ ± = D ˆB, = D ( ˆB ˆχ ± ), = D ˆχ ± where the pair of Darboux vectors D are defined as D = χ c 1 χ ˆB, c 1 = ˆB [ ˆχ ± (χ p b + α ± χ )]. J. D. Gibbon, CSCAMM; October
Vorticity alignment dynamics in fluids & MHD. NCAR - Geophysical Turbulence Program MHD Workshop Thursday June 29, 2006
Vorticity alignment dynamics in fluids & MHD DD Holm Mathematics Department Imperial College London, London SW7 2AZ d.holm@imperial.ac.uk Computer, Computational and Statistical Sciences Los Alamos National
More informationQuaternions and particle dynamics in the Euler fluid equations. J. D. Gibbon. Mathematics Department Imperial College London, London SW7 2AZ
Quaternions and particle dynamics in the Euler fluid equations J. D. Gibbon Mathematics Department Imperial College London, London SW7 2AZ j.d.gibbon@imperial.ac.uk D 2 H-Fest Lausanne, July 2007 In honour
More informationQuaternions and particle dynamics in the Euler fluid equations
23rd Feb 06: http://arxiv.org/abs/nlin.cd/0512034 Quaternions & particle dynamics... 1 Quaternions and particle dynamics in the Euler fluid equations J. D. Gibbon and D. D. Holm Department of Mathematics,
More informationThe three-dimensional Euler fluid equations : where do we stand? J. D. Gibbon. Mathematics Department Imperial College London, London SW7 2AZ
The three-dimensional Euler fluid equations : where do we stand? J. D. Gibbon Mathematics Department, London SW7 2AZ j.d.gibbon@imperial.ac.uk Euler Equations 250 years on 19th June 2007 J. D. Gibbon,
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationRemarks on the blow-up criterion of the 3D Euler equations
Remarks on the blow-up criterion of the 3D Euler equations Dongho Chae Department of Mathematics Sungkyunkwan University Suwon 44-746, Korea e-mail : chae@skku.edu Abstract In this note we prove that the
More informationThis note presents an infinite-dimensional family of exact solutions of the incompressible three-dimensional Euler equations
IMRN International Mathematics Research Notices 2000, No. 9 The Euler Equations and Nonlocal Conservative Riccati Equations Peter Constantin This note presents an infinite-dimensional family of exact solutions
More informationLevel Set Dynamics and the Non-blowup of the 2D Quasi-geostrophic Equation
Level Set Dynamics and the Non-blowup of the 2D Quasi-geostrophic Equation J. Deng, T. Y. Hou, R. Li, X. Yu May 31, 2006 Abstract In this article we apply the technique proposed in Deng-Hou-Yu [7] to study
More informationThe three-dimensional Euler equations: Where do we stand?
Physica D 237 (2008) 1894 1904 www.elsevier.com/locate/physd The three-dimensional Euler equations: Where do we stand? J.D. Gibbon Department of Mathematics, Imperial College London, London SW7 2AZ, United
More informationThe Euler Equations and Non-Local Conservative Riccati Equations
The Euler Equations and Non-Local Conservative Riccati Equations Peter Constantin Department of Mathematics The University of Chicago November 8, 999 The purpose of this brief note is to present an infinite
More informationTRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD
Proceedings of The Second International Conference on Vortex Methods, September 26-28, 21, Istanbul, Turkey TRANSIENT FLOW AROUND A VORTEX RING BY A VORTEX METHOD Teruhiko Kida* Department of Energy Systems
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationu t + u u = p (1.1) u = 0 (1.2)
METHODS AND APPLICATIONS OF ANALYSIS. c 2005 International Press Vol. 12, No. 4, pp. 427 440, December 2005 004 A LEVEL SET FORMULATION FOR THE 3D INCOMPRESSIBLE EULER EQUATIONS JIAN DENG, THOMAS Y. HOU,
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationWeek 3: Differential Geometry of Curves
Week 3: Differential Geometry of Curves Introduction We now know how to differentiate and integrate along curves. This week we explore some of the geometrical properties of curves that can be addressed
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationAn Eulerian-Lagrangian Approach for Incompressible Fluids: Local Theory
An Eulerian-Lagrangian Approach for Incompressible Fluids: Local Theory Peter Constantin Department of Mathematics The University of Chicago December 22, 2000 Abstract We study a formulation of the incompressible
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationSPLIT QUATERNIONS and CANAL SURFACES. in MINKOWSKI 3-SPACE
INTERNATIONAL JOURNAL OF GEOMETRY Vol. 5 (016, No., 51-61 SPLIT QUATERNIONS and CANAL SURFACES in MINKOWSKI 3-SPACE SELAHATTIN ASLAN and YUSUF YAYLI Abstract. A canal surface is the envelope of a one-parameter
More informationSolutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations
D. D. Holm Solutions to M3-4A16 Assessed Problems # 3 15 Dec 2010 1 Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations Vector notation for differential basis elements:
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationContents. 1. Introduction
FUNDAMENTAL THEOREM OF THE LOCAL THEORY OF CURVES KAIXIN WANG Abstract. In this expository paper, we present the fundamental theorem of the local theory of curves along with a detailed proof. We first
More informationBlow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation
Blow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation Dong Li a,1 a School of Mathematics, Institute for Advanced Study, Einstein Drive, Princeton, NJ 854,
More informationVortex knots dynamics and momenta of a tangle:
Lecture 2 Vortex knots dynamics and momenta of a tangle: - Localized Induction Approximation (LIA) and Non-Linear Schrödinger (NLS) equation - Integrable vortex dynamics and LIA hierarchy - Torus knot
More informationBlow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem
Blow-up or No Blow-up? Fluid Dynamic Perspective of the Clay Millennium Problem Thomas Y. Hou Applied and Comput. Mathematics, Caltech Research was supported by NSF Theodore Y. Wu Lecture in Aerospace,
More informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationOn Global Well-Posedness of the Lagrangian Averaged Euler Equations
On Global Well-Posedness of the Lagrangian Averaged Euler Equations Thomas Y. Hou Congming Li Abstract We study the global well-posedness of the Lagrangian averaged Euler equations in three dimensions.
More informationGoing with the flow: A study of Lagrangian derivatives
1 Going with the flow: A study of Lagrangian derivatives Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc/ 12 February
More informationTopics in Fluid Dynamics: Classical physics and recent mathematics
Topics in Fluid Dynamics: Classical physics and recent mathematics Toan T. Nguyen 1,2 Penn State University Graduate Student Seminar @ PSU Jan 18th, 2018 1 Homepage: http://math.psu.edu/nguyen 2 Math blog:
More informationPoint Vortex Dynamics in Two Dimensions
Spring School on Fluid Mechanics and Geophysics of Environmental Hazards 9 April to May, 9 Point Vortex Dynamics in Two Dimensions Ruth Musgrave, Mostafa Moghaddami, Victor Avsarkisov, Ruoqian Wang, Wei
More informationAPPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS
1 APPROXIMATE PERIODIC SOLUTIONS for the RAPIDLY ROTATING SHALLOW-WATER and RELATED EQUATIONS BIN CHENG Department of Mathematics University of Michigan Ann Arbor, MI 48109 E-mail:bincheng@umich.edu EITAN
More informationMHD dynamo generation via Riemannian soliton theory
arxiv:physics/0510057v1 [physics.plasm-ph] 7 Oct 2005 MHD dynamo generation via Riemannian soliton theory L.C. Garcia de Andrade 1 Abstract Heisenberg spin equation equivalence to nonlinear Schrödinger
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More information[#1] Exercises in exterior calculus operations
D. D. Holm M3-4A16 Assessed Problems # 3 Due when class starts 13 Dec 2012 1 M3-4A16 Assessed Problems # 3 Do all four problems [#1] Exercises in exterior calculus operations Vector notation for differential
More informationSurface x(u, v) and curve α(t) on it given by u(t) & v(t). Math 4140/5530: Differential Geometry
Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt Surface x(u, v) and curve α(t) on it given by u(t) & v(t). α du dv (t) x u dt + x v dt ( ds dt )2 Surface x(u, v)
More informationThe 3D Euler and 2D surface quasi-geostrophic equations
The 3D Euler and 2D surface quasi-geostrophic equations organized by Peter Constantin, Diego Cordoba, and Jiahong Wu Workshop Summary This workshop focused on the fundamental problem of whether classical
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationLecture Notes Introduction to Vector Analysis MATH 332
Lecture Notes Introduction to Vector Analysis MATH 332 Instructor: Ivan Avramidi Textbook: H. F. Davis and A. D. Snider, (WCB Publishers, 1995) New Mexico Institute of Mining and Technology Socorro, NM
More informationEvolution of sharp fronts for the surface quasi-geostrophic equation
Evolution of sharp fronts for the surface quasi-geostrophic equation José Luis Rodrigo University of Warwick supported by ERC Consolidator 616797 UIMP - Frontiers of Mathematics IV July 23rd 2015 José
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationGauge Fixing and Constrained Dynamics in Numerical Relativity
Gauge Fixing and Constrained Dynamics in Numerical Relativity Jon Allen The Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation is reviewed. Gauge freedom is discussed and
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationVector Algebra and Geometry. r = a + tb
Vector Algebra and Geometry Differentiation of Vectors Vector - valued functions of a real variable We have met the equation of a straight line in the form r = a + tb r therefore varies with the real variables
More informationA STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
A STOCHASTIC LAGRANGIAN REPRESENTATION OF THE 3-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS PETER CONSTANTIN AND GAUTAM IYER Abstract. In this paper we derive a probabilistic representation of the
More informatione 2 = e 1 = e 3 = v 1 (v 2 v 3 ) = det(v 1, v 2, v 3 ).
3. Frames In 3D space, a sequence of 3 linearly independent vectors v 1, v 2, v 3 is called a frame, since it gives a coordinate system (a frame of reference). Any vector v can be written as a linear combination
More informationStriated Regularity of Velocity for the Euler Equations
Striated Regularity of Velocity for the Euler Equations JIM KELLIHER 1 with HANTAEK BAE 2 1 University of California Riverside 2 Ulsan National Institute of Science and Technology, Korea 2015 SIAM Conference
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationThe Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations
The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan
More informationMATH 332: Vector Analysis Summer 2005 Homework
MATH 332, (Vector Analysis), Summer 2005: Homework 1 Instructor: Ivan Avramidi MATH 332: Vector Analysis Summer 2005 Homework Set 1. (Scalar Product, Equation of a Plane, Vector Product) Sections: 1.9,
More informationDIFFERENTIAL GEOMETRY AND THE QUATERNIONS. Nigel Hitchin (Oxford) The Chern Lectures Berkeley April 9th-18th 2013
DIFFERENTIAL GEOMETRY AND THE QUATERNIONS Nigel Hitchin (Oxford) The Chern Lectures Berkeley April 9th-18th 2013 16th October 1843 ON RIEMANNIAN MANIFOLDS OF FOUR DIMENSIONS 1 SHIING-SHEN CHERN Introduction.
More informationMathematical Hydrodynamics
Mathematical Hydrodynamics Ya G. Sinai 1. Introduction Mathematical hydrodynamics deals basically with Navier-Stokes and Euler systems. In the d-dimensional case and incompressible fluids these are the
More informationVortex stretching and anisotropic diffusion in the 3D NSE
University of Virginia SIAM APDE 2013 ; Lake Buena Vista, December 7 3D Navier-Stokes equations (NSE) describing a flow of 3D incompressible viscous fluid read u t + (u )u = p + ν u, supplemented with
More informationOn a variational principle for Beltrami flows
On a variational principle for Beltrami flows Rafael González, Andrea Costa, and E. Sergio Santini Citation: Physics of Fluids 22, 074102 (2010); doi: 10.1063/1.3460297 View online: http://dx.doi.org/10.1063/1.3460297
More informationNewman-Penrose formalism in higher dimensions
Newman-Penrose formalism in higher dimensions V. Pravda various parts in collaboration with: A. Coley, R. Milson, M. Ortaggio and A. Pravdová Introduction - algebraic classification in four dimensions
More informationUNIVERSITY OF DUBLIN
UNIVERSITY OF DUBLIN TRINITY COLLEGE JS & SS Mathematics SS Theoretical Physics SS TSM Mathematics Faculty of Engineering, Mathematics and Science school of mathematics Trinity Term 2015 Module MA3429
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationDynamics of the Coarse-Grained Vorticity
(B) Dynamics of the Coarse-Grained Vorticity See T & L, Section 3.3 Since our interest is mainly in inertial-range turbulence dynamics, we shall consider first the equations of the coarse-grained vorticity.
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationTurbulent spectra generated by singularities
Turbulent spectra generated by singularities p. Turbulent spectra generated by singularities E.A. Kuznetsov Landau Institute for Theoretical Physics, Moscow Russia In collaboration with: V. Naulin A.H.
More information1.1 Single Variable Calculus versus Multivariable Calculus Rectangular Coordinate Systems... 4
MATH2202 Notebook 1 Fall 2015/2016 prepared by Professor Jenny Baglivo Contents 1 MATH2202 Notebook 1 3 1.1 Single Variable Calculus versus Multivariable Calculus................... 3 1.2 Rectangular Coordinate
More informationInterpolated Rigid-Body Motions and Robotics
Interpolated Rigid-Body Motions and Robotics J.M. Selig London South Bank University and Yuanqing Wu Shanghai Jiaotong University. IROS Beijing 2006 p.1/22 Introduction Interpolation of rigid motions important
More informationVortex stretching in incompressible and compressible fluids
Vortex stretching in incompressible and compressible fluids Esteban G. Tabak, Fluid Dynamics II, Spring 00 1 Introduction The primitive form of the incompressible Euler equations is given by ( ) du P =
More informationMetrics and Curvature
Metrics and Curvature How to measure curvature? Metrics Euclidian/Minkowski Curved spaces General 4 dimensional space Cosmological principle Homogeneity and isotropy: evidence Robertson-Walker metrics
More informationBlow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations
Blow-up or No Blow-up? the Role of Convection in 3D Incompressible Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech Joint work with Zhen Lei; Congming Li, Ruo Li, and Guo
More informationGeometric Modelling Summer 2016
Geometric Modelling Summer 2016 Exercises Benjamin Karer M.Sc. http://gfx.uni-kl.de/~gm Benjamin Karer M.Sc. Geometric Modelling Summer 2016 1 Dierential Geometry Benjamin Karer M.Sc. Geometric Modelling
More informationGeometric approximation of curves and singularities of secant maps Ghosh, Sunayana
University of Groningen Geometric approximation of curves and singularities of secant maps Ghosh, Sunayana IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish
More informationTwisting versus bending in quantum waveguides
Twisting versus bending in quantum waveguides David KREJČIŘÍK Nuclear Physics Institute, Academy of Sciences, Řež, Czech Republic http://gemma.ujf.cas.cz/ david/ Based on: [Chenaud, Duclos, Freitas, D.K.]
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
More informationExercise 1 (Formula for connection 1-forms) Using the first structure equation, show that
1 Stokes s Theorem Let D R 2 be a connected compact smooth domain, so that D is a smooth embedded circle. Given a smooth function f : D R, define fdx dy fdxdy, D where the left-hand side is the integral
More informationSmoluchowski Navier-Stokes Systems
Smoluchowski Navier-Stokes Systems Peter Constantin Mathematics, U. of Chicago CSCAMM, April 18, 2007 Outline: 1. Navier-Stokes 2. Onsager and Smoluchowski 3. Coupled System Fluid: Navier Stokes Equation
More informationsin(α + θ) = sin α cos θ + cos α sin θ cos(α + θ) = cos α cos θ sin α sin θ
Rotations in the 2D Plane Trigonometric addition formulas: sin(α + θ) = sin α cos θ + cos α sin θ cos(α + θ) = cos α cos θ sin α sin θ Rotate coordinates by angle θ: 1. Start with x = r cos α y = r sin
More informationSpace curves, vector fields and strange surfaces. Simon Salamon
Space curves, vector fields and strange surfaces Simon Salamon Lezione Lagrangiana, 26 May 2016 Curves in space 1 A curve is the path of a particle moving continuously in space. Its position at time t
More informationOn explicit integration of two non-holonomic problems
On explicit integration of two non-holonomic problems Alexey V. Borisov 1 1 Institute of Computer Sciences, Izhevsk, Russia Generalized Chaplygin systems Equations of motion of the generalized Chaplygin
More informationFUNDAMENTAL CONCEPTS IN CONTINUUM MECHANICS
PART I FUNDAMENTAL CONCEPTS IN CONTINUUM MECHANICS CHAPTER ONE Describing the motion ofa system: geometry and kinematics 1.1. Deformations The purpose of mechanics is to study and describe the motion of
More informationOptimal Control Problems in Eye Movement. Bijoy K. Ghosh Department of Mathematics and Statistics Texas Tech University, Lubbock, TX
Optimal Control Problems in Eye Movement Bijoy K. Ghosh Department of Mathematics and Statistics Texas Tech University, Lubbock, TX We study the human oculomotor system as a simple mechanical control
More informationOHSx XM521 Multivariable Differential Calculus: Homework Solutions 14.1
OHSx XM5 Multivariable Differential Calculus: Homework Solutions 4. (8) Describe the graph of the equation. r = i + tj + (t )k. Solution: Let y(t) = t, so that z(t) = t = y. In the yz-plane, this is just
More informationUnderstanding Quaternions: Rotations, Reflections, and Perspective Projections. Ron Goldman Department of Computer Science Rice University
Understanding Quaternions: Rotations, Reflections, and Perspective Projections Ron Goldman Department of Computer Science Rice University The invention of the calculus of quaternions is a step towards
More informationDeformation and Symmetry in the Inviscid SQG and the 3D Euler Equations
J Nonlinear Sci 1 :665 688 DOI 1.17/s33-1-914-7 Deformation and Symmetry in the Inviscid SQG and the 3D Euler Equations Dongho Chae Peter Constantin Jiahong Wu Received: 1 July 11 / Accepted: February
More informationSuperintegrable 3D systems in a magnetic field and Cartesian separation of variables
Superintegrable 3D systems in a magnetic field and Cartesian separation of variables in collaboration with L. Šnobl Czech Technical University in Prague GSD 2017, June 5-10, S. Marinella (Roma), Italy
More information3 Rigid Spacecraft Attitude Control
1 3 Rigid Spacecraft Attitude Control Consider a rigid spacecraft with body-fixed frame F b with origin O at the mass centre. Let ω denote the angular velocity of F b with respect to an inertial frame
More informationEuler Characteristic of Two-Dimensional Manifolds
Euler Characteristic of Two-Dimensional Manifolds M. Hafiz Khusyairi August 2008 In this work we will discuss an important notion from topology, namely Euler Characteristic and we will discuss several
More informationQuasi-two Dimensional Hydrodynamics and Interaction of Vortex Tubes
Quasi-two Dimensional Hydrodynamics and Interaction of Vortex Tubes Vladimir Zakharov 1 Introduction This paper is long overdue. The most of the results presented here were obtained in 1986-87. Just a
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationTHE DARBOUX TRIHEDRONS OF REGULAR CURVES ON A REGULAR SURFACE
International lectronic Journal of eometry Volume 7 No 2 pp 61-71 (2014) c IJ TH DARBOUX TRIHDRONS OF RULAR CURVS ON A RULAR SURFAC MRAH TUNÇ AND MİN OZYILMAZ (Communicated by Levent KULA) Abstract In
More informationQuaternions and Arithmetic. Colloquium, UCSD, October 27, 2005
Quaternions and Arithmetic Colloquium, UCSD, October 27, 2005 This talk is available from www.math.mcgill.ca/goren Quaternions came from Hamilton after his really good work had been done; and, though beautifully
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationOn the leapfrogging phenomenon in fluid mechanics
On the leapfrogging phenomenon in fluid mechanics Didier Smets Université Pierre et Marie Curie - Paris Based on works with Robert L. Jerrard U. of Toronto) CIRM, Luminy, June 27th 2016 1 / 22 Single vortex
More informationON THE RULED SURFACES WHOSE FRAME IS THE BISHOP FRAME IN THE EUCLIDEAN 3 SPACE. 1. Introduction
International Electronic Journal of Geometry Volume 6 No.2 pp. 110 117 (2013) c IEJG ON THE RULED SURFACES WHOSE FRAME IS THE BISHOP FRAME IN THE EUCLIDEAN 3 SPACE ŞEYDA KILIÇOĞLU, H. HILMI HACISALIHOĞLU
More informationChapter 6: Vector Analysis
Chapter 6: Vector Analysis We use derivatives and various products of vectors in all areas of physics. For example, Newton s 2nd law is F = m d2 r. In electricity dt 2 and magnetism, we need surface and
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationPart II. Classical Dynamics. Year
Part II Year 28 27 26 25 24 23 22 21 20 2009 2008 2007 2006 2005 28 Paper 1, Section I 8B Derive Hamilton s equations from an action principle. 22 Consider a two-dimensional phase space with the Hamiltonian
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
More informationFINITE TIME BLOW-UP FOR A DYADIC MODEL OF THE EULER EQUATIONS
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 357, Number 2, Pages 695 708 S 0002-9947(04)03532-9 Article electronically published on March 12, 2004 FINITE TIME BLOW-UP FOR A DYADIC MODEL OF
More informationQuaternions. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Quaternions Semester 1, / 40
Quaternions Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-2017 B. Bona (DAUIN) Quaternions Semester 1, 2016-2017 1 / 40 Introduction Complex numbers with unit norm can be used as rotation operators
More informationBLACK HOLES (ADVANCED GENERAL RELATIV- ITY)
Imperial College London MSc EXAMINATION May 2015 BLACK HOLES (ADVANCED GENERAL RELATIV- ITY) For MSc students, including QFFF students Wednesday, 13th May 2015: 14:00 17:00 Answer Question 1 (40%) and
More informationPAPER 311 BLACK HOLES
MATHEMATICAL TRIPOS Part III Friday, 8 June, 018 9:00 am to 1:00 pm PAPER 311 BLACK HOLES Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationF1.9AB2 1. r 2 θ2 + sin 2 α. and. p θ = mr 2 θ. p2 θ. (d) In light of the information in part (c) above, we can express the Hamiltonian in the form
F1.9AB2 1 Question 1 (20 Marks) A cone of semi-angle α has its axis vertical and vertex downwards, as in Figure 1 (overleaf). A point mass m slides without friction on the inside of the cone under the
More information