Centrum voor Wiskunde en Informatica
|
|
- Pearl Gardner
- 5 years ago
- Views:
Transcription
1 Centrum voor Wiskune en Informatica Moelling, Analysis an Simulation Moelling, Analysis an Simulation Conservation properties of smoothe particle hyroynamics applie to the shallow water equations J.E. Frank, S. Reich REPORT MAS-R11 AUGUST 20
2 CWI is the National Research Institute for Mathematics an Computer Science. It is sponsore by the Netherlans Organization for Scientific Research (NWO). CWI is a founing member of ERCIM, the European Research Consortium for Informatics an Mathematics. CWI's research has a theme-oriente structure an is groupe into four clusters. Liste below are the names of the clusters an in parentheses their acronyms. Probability, Networks an Algorithms (PNA) Software Engineering (SEN) Moelling, Analysis an Simulation (MAS) Information Systems (INS) Copyright 20, Stichting Centrum voor Wiskune en Informatica P.O. Box 94079, 1090 GB Amsteram (NL) Kruislaan 413, 1098 SJ Amsteram (NL) Telephone Telefax ISSN
3 Conservation Properties of Smoothe Particle Hyroynamics Applie to the Shallow Water Equations Jason Frank CWI P.O. Box 94079, 1090 GB Amsteram, The Netherlans Sebastian Reich epartment of Mathematics, Imperial College 180 Queen s Gate, Lonon, SW7 2BZ, UK s.reich@ic.ac.uk ABSTRACT Conservation of potential vorticity (PV) an Kelvin s circulation theorem express two of the funamental concepts in ieal flui ynamics. In this note, we iscuss these two concepts in the context of the Smoothe Particle Hyroynamics (SPH) metho. We show that, if interprete in an appropriate way, one can make a statement about conservation of circulation an PV in a particle flui moel such as SPH. We also inicate some limitations where the analogy with ieal flui ynamics breaks own Mathematics Subject Classification: 76M28 Keywors an Phrases: geophysical flui ynamics, potential vorticity conserving methos, geometric methos, smoothe particle hyroynamics Note: Work of J. Frank carrie out uner project MAS Atmospheric Flow an Transport Problems. Partial support by GM is gratefully acknowlege. Note: Partial financial support of S. Reich by EPSRC Grant GR/R09565/ an by European Commission funing for the Research Training Network Mechanics an Symmetry in Europe is gratefully acknowlege. 1. Introuction Smoothe Particle Hyroynamics (SPH) [7, 10] is a popular particle base metho for computational flui ynamics. It is well-known that the SPH metho can be erive from a variational principle (see, e.g., [10]) or, in other wors, can be given a Hamiltonian structure. In this note, we aress the important question of conservation of potential vorticity an circulation in SPH simulations. These conserve quantities follow from the unerlying relabeling symmetry of ieal flui ynamics an are funamental to the long-time solution behavior [8]. For simplicity of exposition, we consier the two-imensional shallow water equations (SWEs): t u = c 0 x h, (1.1) t h = h x u, (1.2) where u =(u,v) T is the horizontal velocity fiel, h is the layer epth, c 0 > 0 is an appropriate constant, an t = t + u x is the material time erivative. 2. Conserve Quantities in Ieal Shallow Water Flows The SWEs (1.1)-(1.2) possess a number of conserve quantities which are important for the longtime ynamics. Let us start with the vorticity ζ = x u. Using x t u = x u t +( x u)( x u)+u x ( x u)=0,
4 2. Conserve Quantities in Ieal Shallow Water Flows 2 it is easy to conclue that vorticity satisfies the continuity equation t ζ = ζ x u. (2.1) The ratio of ζ to h, i.e., q = ζ/h, is calle the potential vorticity (PV) [8]. The PV fiel q is materially conserve since { t q = h 1 t ζ q } t h =0. Using this result, it follows that (1.2) an (2.1) are special cases of an infinite family of continuity equations t [hf(q)] = x [hf(q)u], (2.2) where f is an arbitrary (smooth) function of q. Next we introuce flui particle labels a =(a,b) R 2 an efine particle positions x(a,t) R 2 via a iffeomorphism. The labels are fixe for each particle [8], i.e. t a = 0. The eterminant of the 2 2 Jacobian matrix satisfies x a = (x, y) (a, b) h x a = h o, where h o (a) is a time-inepenent function [8]. Let us now iscuss the concept of circulation. Take a close loop S = {a(s)} s S 1 in label space an consier the particle locations x(s) =x(a(s)) parameterize by s S 1. By efinition, the loop {x(s)} s S 1 in configuration space is avecte along the velocity fiel, i.e. t x(s)=u(x(s)). Kelvin s circulation theorem [8] states that u x s s =0. (2.3) t Inee, we obtain ( ) ( ) u x s s = t t u x s s + u s t x s = c 0 x h x s s + u u s s ( ) 1 = 2 (u u) s c 0 h s s = 0. Let A enote the area enclose by S in label space an R its image in x-space. Then Stokes theorem applie to (2.3) yiels x u) x a a b = t A( ( x u)x y =0. (2.4) t R
5 3. Raial Basis Functions an Smoothe Particle Hyroynamics 3 Since A is arbitrary, the left sie of this equation yiels another statement of PV conservation since ( x u) x a = h o ζ h = h oq, an h o /t = 0. Similarly, after applying the transport theorem [4] to the right equality in (2.4), a secon appeal to the arbitaryness of A yiels the vorticity equation (2.1). The SWEs also conserve energy E = 1 h o {u u + c 0 h} a b = 1 h {u u + c 0 h} x y 2 2 an the symplectic two-form ω := (u x)a b [8]. 3. Raial Basis Functions an Smoothe Particle Hyroynamics Assume that a set of Lagrangian particles with positions {X k (t); X k R 2 } is given as a function of time an that t X k = U k, where U k is the velocity of the particle. Then the time evolution of a quantity g, satisfying a continuity equation t g = x [gu], can be approximate by (3.1) g(x,t)= k γ k ψ( x X k (t) 2 ). Here {γ k } are constants etermine by the initial g(x)fielanψ(r) is an appropriate raial basis function [3], for example, ψ(r 2 )= ( ( r ) ) 2 1/2 + c 2, R c, R > 0 two parameters. Let us apply this iea to the layer-epth h, i.e., we introuce the approximation h(x,t)= k m k ψ( x X k (t) 2 ) (3.2) an assume that h(x,t) > 0. Then each particle contributes the fraction ρ k (x,t):= m kψ( x X k (t) 2 ) h(x,t) (3.3) to the total layer-epth. These fractions form a partition of unity, i.e. ρ k (x,t)=1. k Hence they can be use to interpolate ata from the particle locations to any x R 2.Inparticular, we efine an interpolate Eulerian velocity fiel u(x,t):= k ρ k (x,t)u k (t) (3.4)
6 4. Conservation of Enstrophy an Circulation 4 an a layer epth flux ensity h(x,t)u(x,t)= k m k ψ( x X k (t) 2 )U k (t). Using (3.1), it is now easily verifie that t h(x,t)+ x [h(x,t)u(x,t)] = 0. (3.5) It follows that the layer epth approximation (3.2) exactly satisfies the continuity equation (1.2) uner the flow of the formally efine velocity fiel (3.4). In general the particle avection velocity is ifferent from the interpolate velocity, i.e U k u(x k ). We note that the moification suggeste in [6] to avoi penetration in SPH correspons 1 to avecting the particles in the velocity fiel (3.4). Hence, (3.1) an (3.2) provie an approximation to the continuity equation (1.2). To get a close system of iscretize equations, we still have to approximate the momentum equation (1.1). For example, one can use t U k = c 0 x h(x = X k,t)= c 0 m j Xk ψ( X k X j 2 ). (3.6) j The equations (3.1), (3.2), an (3.6) are commonly known as a smoothe particle hyroynamics (SPH) approximation [7] to the SWEs (1.1)-(1.2). A variant of these equations is obtaine by replacing (3.6) by t U k = c 0 2 (m k + m j ) Xk ψ( X k X j 2 ). (3.7) j The equations (3.1), (3.2), an (3.7) are now canonical with Hamiltonian (energy) H = 1 2 k U k 2 + c 0 2 m k ψ( X l X k 2 ) l,k an symplectic structure ω = k U k X k. If we assume that the particle locations {X k } have been chosen such that m k = m = const., then (3.6) an (3.7) coincie. A numerical time-stepping scheme is obtaine by noting that t X k = U k, t U k = 0 can be solve exactly an that the associate time evolution of h(x,t) exactly satisfies (3.5). Similarly, equation (3.7) an t X k = 0 can also be integrate exactly since h t = 0. A composition of these exact propagators leas to a symplectic time-stepping scheme [9] implying goo long-time energy conservation [2]. The same approach using (3.6) can be mae time-symmetric but is not symplectic, in general. 4. Conservation of Enstrophy an Circulation The PV fiel q is materially conserve. This suggests assigning a constant PV value q k to each particle location X k, thereby trivially enforcing the PV conservation law q/t =0. Usingthe fractions ρ k, we obtain the interpolate fiels f(q(x,t)) := k f(q k )ρ k (x,t), (4.1) 1 More precisely, the formulation Eqn. (2.6) in [6] avocates particle avection in an interpolate velocity fiel base on a generic kernel. Taking this kernel to be ψ yiels the avection fiel u(x k ).
7 4. Conservation of Enstrophy an Circulation 5 where f is again an arbitrary function of q. It is easily verifie that the approximation (4.1) satisfies the continuity equation { } t [hf(q)](x,t)= x m k f(q k )ψ( x X k (t) 2 )U k k corresponing to (2.2). Uner the given perioic bounary conitions, this iscrete conservation law implies the exact conservation of the generalize enstrophies Q f = hf(q)x y. Since this is also true for the split equations of motion use for the time-stepping, the overall space-time iscretization conserves enstrophy inepenently of whether (3.6) or (3.7) is use. One shoul, however, be aware that the prouct hq is not equivalent to the vorticity of the interpolate velocity fiel u as is true for the SWEs. See below. The circulation is also conserve, but for a ifferent velocity fiel U(x,t) efine as follows: Let U(x,0) be any initial velocity fiel satisfying U(X k,0) = U k (0) at the particle locations. (For example, suppose the particle velocities at t = 0 are given as a continuous function.) An let U(x,t) evolve uner the solution of the SPH flow ue to (3.6) accoring to t U(x)= c 0 X h(x,t)= c 0 m k x ψ( x X k (t) 2 ). k Note that for this velocity fiel it oes hol that U(X k (t)) = U k (t). Figure 1 illustrates the relationship between the velocity fiels u an U. X k U k u(x(s)) u(x k ) S X(s) U(s) X l u(x l) U l Figure 1: A close curve avecte with the flow, illustrating the velocity fiels u an U. Now, we efine a curve of Lagrangian points X(s)=X(a(s)) with s S 1 an S = {a(s)} s S 1 being a close loop in label space. The associate loop {X(s)} s S 1 in configuration space is propagate in the velocity fiel U(x,t) accoring to t X(s)=U(X(s)).
8 4. Conservation of Enstrophy an Circulation 6 We assume that X(s)anU(s) are sufficiently ifferentiable. Then Kelvin s circulation theorem (2.3) becomes U X s s =0. (4.2) t Inee, we obtain ( ) ( ) U X s s = t t U X s s + U s t X s = c 0 X h X s s + U U s s ( ) 1 = 2 (U U) s c 0 h s s = 0. This reasoning is no longer vali if the momentum equation (3.7) is use. Conservation of circulation is thus in conflict with conservation of symplectic structure. If we now give eachparticle insie S alabela an let A enote the area enclose by S an let R enote the image of A in x-space, then applying Stokes theorem to (4.2) yiels t A( x U) X a a b = t for which the left sie implies R ( x U)x y =0, (4.3) t {( x U) X a } =0, (4.4) since A is arbitrary. Note that, in principle, the eterminant X a can be compute along a given particle path X k (t) by integrating the linearize SPH equations along X k (t). However, the prouct of X a an h is not equivalent to a time-inepenent function h o (a) an, hence, equation (4.4) is not equivalent to conservation of PV in the stanar sense, i.e. q/t =0. Applying the transport theorem to the right equality of (4.3), an again noting that A is arbitrary, yiels a continuity equation for the relative vorticity of the velocity fiel U: ζ t = x (ζu), ζ = x U, cf. (2.1). However, the vorticity ζ(x,t) will be ifferent from the quantity h(x,t)q(x,t)= k q k m k ψ( x X k (t) 2 ), even if the two are initially equal, because the latter is avecte by the interpolate velocity fiel u. Hence the agreement between u an U is a funamental measure of the quality of an SPH simulation. We woul a that (4.2) an (4.3) are preserve uner time iscretization via a splitting as escribe in the previous section. For a numerical verification of (4.2) one has to represent the loop {X(s)} by a sufficient number of particles { X l } with associate velocities {Ū l }. The integral (4.2) is then approximate by U X s s Ū l ( X l+1 X l ). l Note that the particles { X l } are not part of the SPH iscretization but are instea obtaine via post-processing of the SPH solutions.
9 5. Concluing Remarks 7 5. Concluing Remarks In this note, we have shown that the SPH metho with (3.6) satisfies a Kelvin circulation law. However, conservation of circulation is in conflict with conservation of symplecticness unless the weights m k are all equal. We have also shown that conservation of circulation implies a form of PV conservation an that the generalize enstrophies Q f are also preserve. Note that the numerical approximations o not, in general, satisfy the funamental relations ζ = hq an h x a = h o.the consequences of this limitation require further investigations. The results of this paper easily generalize to the rotating SWEs t u = f 0u c 0 x h, t h = h x u, where u =( v,u) T an f 0 /2 is the angular velocity of the reference plane. Potential vorticity is now efine by q = x u + f 0 h an Kelvin s circulation theorem becomes t ( u + f 0 2 x ) x s s = t A ( x u + f 0 ) x a a b = t R ( x u + f 0 ) x y =0. We wish to mention the Balance Particle-Mesh (BPM) metho of [5] which uses raial basis functions to approximate the absolute vorticity ω = x u + f 0. See [5] for the geometric properties of the BPM metho. Kelvin s circulation theorem also applies to three-imensional ieal fluis while conservation of PV takes a more complicate form (see [8]). Again, conservation of circulation can be shown for the SPH metho in the same manner as outline in this note for two-imensional fluis. In fact, the concept of circulation even applies to molecular simulations of a mono-atomic liqui [1] with Hamiltonian H = 1 P k 2 + φ( X k X l ), 2m l>k k where m is the atomic mass an φ(r) an interaction potential. We introuce the function ρ(x,t)= l φ( x X l (t) ), x R 3, an note that Newton s law is equivalent to t P k = Xk H = x ρ(x = X k,t). We also have t X k = 1 m P k. Applying the notations of 4 an replacing U by P, we obtain the circulation theorem t P X s s =0 an, in two imensions, conservation of vorticity per control area, i.e., ( x P)x y. t R One shoul keep in min that φ(r) is often singular at r = 0 an, hence, ρ(x,t) is not efine for x = X k. However, one can replace φ(r) by a smooth truncation φ(r) such that φ(r) =φ(r) for r r o an ρ (0) = 0, ρ(0) <. Herer o is chosen such that X i (t) X j (t) >r o for all t 0 an all i j.
10 References 8 References 1. M. P. Allen &. J. Tilesley, Computer Simulation of Liquis, Oxfor University Press, Oxfor, G. Benettin & A. Giorgilli, On the Hamiltonian interpolation of near to the ientity symplectic mappings with application to symplectic integration algorithms, J. Stat. Phys. 74 (1994), M.. Buhmann, Raial basis functions, Acta Numerica (2000), A.J.Chorin&J.E.Marsen, A Mathematical Introuction to Flui Mechanics, Springer- Verlag, New York, J. Frank & S. Reich, A particle-mesh metho for the shallow water equations near geostrophic balance, submitte. 6. J. J. Monaghan, On the problem of penetration in particle methos, J. Comput. Phys. 82 (1989), J. J. Monaghan, Smoothe particle hyroynamics, Ann. Rev. Astron. Astrophys. 30 (1992), R. Salmon, Lectures on Geophysical Flui ynamics, Oxfor University Press,Oxfor, J. M. Sanz-Serna & M. P. Calvo, Numerical Hamiltonian Problems, Chapman an Hall, Lonon, H. Yserentant, A new class of particle methos., Numer. Math. 76, 1997,
CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationThe Hamiltonian particle-mesh method for the spherical shallow water equations
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 5: 89 95 (004) Publishe online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.100/asl.70 The Hamiltonian particle-mesh metho for the spherical
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationPhysics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1
Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems
More informationEuler Equations: derivation, basic invariants and formulae
Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationIntegration Review. May 11, 2013
Integration Review May 11, 2013 Goals: Review the funamental theorem of calculus. Review u-substitution. Review integration by parts. Do lots of integration eamples. 1 Funamental Theorem of Calculus In
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationGeometric Algebra Approach to Fluid Dynamics
Geometric Algebra Approach to Flui Dynamics Carsten Cibura an Dietmar Hilenbran Abstract In this work we will use geometric algebra to prove a number of well known theorems central to the fiel of flui
More informationWUCHEN LI AND STANLEY OSHER
CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationSolution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010
NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationPDE Notes, Lecture #11
PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationDarboux s theorem and symplectic geometry
Darboux s theorem an symplectic geometry Liang, Feng May 9, 2014 Abstract Symplectic geometry is a very important branch of ifferential geometry, it is a special case of poisson geometry, an coul also
More informationDiscrete Operators in Canonical Domains
Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:
More informationThe proper definition of the added mass for the water entry problem
The proper efinition of the ae mass for the water entry problem Leonaro Casetta lecasetta@ig.com.br Celso P. Pesce ceppesce@usp.br LIE&MO lui-structure Interaction an Offshore Mechanics Laboratory Mechanical
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS
ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics
More informationA new proof of the sharpness of the phase transition for Bernoulli percolation on Z d
A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition
More informationON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS
ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ADRIAN CONSTANTIN AND BORIS KOLEV Abstract. Accoring to the principle of least action, the spatially perioic motions of one-imensional
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationC e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a Modelling, Analysis and Simulation Modelling, Analysis and Simulation Bilinear forms for the recovery-based discontinuous Galerkin method
More information2.20 Marine Hydrodynamics Lecture 3
2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationThe canonical controllers and regular interconnection
Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,
More informationAbstract A nonlinear partial differential equation of the following form is considered:
M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent
More informationCurvature, Conformal Mapping, and 2D Stationary Fluid Flows. Michael Taylor
Curvature, Conformal Mapping, an 2D Stationary Flui Flows Michael Taylor 1. Introuction Let Ω be a 2D Riemannian manifol possibly with bounary). Assume Ω is oriente, with J enoting counterclockwise rotation
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationA geometric theory of selective decay in fluids with advected quantities
A geometric theory of selective ecay in fluis with avecte quantities arxiv:1310.4543v2 [physics.plasm-ph] 20 Oct 2018 François Gay-Balmaz 1 an arryl. Holm 2 PACS Numbers: Geometric mechanics 02.40.Yy;
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationSemicompressible Ocean Thermodynamics and Boussinesq Energy Conservation
Article Semicompressible Ocean Thermoynamics an Boussinesq Energy Conservation William K. Dewar 1, *, Joseph Schoonover 1, Trevor McDougall 2 an Rupert Klein 3 1 Department of Earth, Ocean, an Atmospheric
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationWELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES
Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationCentrum voor Wiskunde en Informatica
Centrum voor Wiskunde en Informatica Modelling, Analysis and Simulation Modelling, Analysis and Simulation Two-point Taylor expansions of analytic functions J.L. Lópe, N.M. Temme REPORT MAS-R0 APRIL 30,
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa
Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik
More informationBEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS. Mauro Boccadoro Magnus Egerstedt Paolo Valigi Yorai Wardi
BEYOND THE CONSTRUCTION OF OPTIMAL SWITCHING SURFACES FOR AUTONOMOUS HYBRID SYSTEMS Mauro Boccaoro Magnus Egerstet Paolo Valigi Yorai Wari {boccaoro,valigi}@iei.unipg.it Dipartimento i Ingegneria Elettronica
More informationarxiv:nlin/ v1 [nlin.cd] 21 Mar 2002
Entropy prouction of iffusion in spatially perioic eterministic systems arxiv:nlin/0203046v [nlin.cd] 2 Mar 2002 J. R. Dorfman, P. Gaspar, 2 an T. Gilbert 3 Department of Physics an Institute for Physical
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationThe Generalized Incompressible Navier-Stokes Equations in Besov Spaces
Dynamics of PDE, Vol1, No4, 381-400, 2004 The Generalize Incompressible Navier-Stokes Equations in Besov Spaces Jiahong Wu Communicate by Charles Li, receive July 21, 2004 Abstract This paper is concerne
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationLECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES
LECTURE 1: BASIC CONCEPTS, PROBLEMS, AND EXAMPLES WEIMIN CHEN, UMASS, SPRING 07 In this lecture we give a general introuction to the basic concepts an some of the funamental problems in symplectic geometry/topology,
More informationSpectral Flow, the Magnus Force, and the. Josephson-Anderson Relation
Spectral Flow, the Magnus Force, an the arxiv:con-mat/9602094v1 16 Feb 1996 Josephson-Anerson Relation P. Ao Department of Theoretical Physics Umeå University, S-901 87, Umeå, SWEDEN October 18, 2018 Abstract
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson
JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises
More informationTransmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency
Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationProblem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs
Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationDelocalization of boundary states in disordered topological insulators
Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators
More informationLeft-invariant extended Kalman filter and attitude estimation
Left-invariant extene Kalman filter an attitue estimation Silvere Bonnabel Abstract We consier a left-invariant ynamics on a Lie group. One way to efine riving an observation noises is to make them preserve
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationRemarks on Singularities, Dimension and Energy Dissipation for Ideal Hydrodynamics and MHD
Commun. Math. Phys. 184, 443 455 (1997) Communications in Mathematical Physics c Springer-Verlag 1997 Remarks on Singularities, imension an Energy issipation for Ieal Hyroynamics an MH Russel E. Caflisch,
More informationCalculus in the AP Physics C Course The Derivative
Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationDiscrete Hamilton Jacobi Theory and Discrete Optimal Control
49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL. Qiang Du, Manlin Li and Chun Liu
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 8, Number 3, October 27 pp. 539 556 ANALYSIS OF A PHASE FIELD NAVIER-STOKES VESICLE-FLUID INTERACTION MODEL Qiang
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationTractability results for weighted Banach spaces of smooth functions
Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March
More informationNonlinear Schrödinger equation with a white-noise potential: Phase-space approach to spread and singularity
Physica D 212 (2005) 195 204 www.elsevier.com/locate/phys Nonlinear Schröinger equation with a white-noise potential: Phase-space approach to sprea an singularity Albert C. Fannjiang Department of Mathematics,
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More information