CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
|
|
- Deborah Allen
- 6 years ago
- Views:
Transcription
1 BIT /00/ $ ?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN REICH 2 1 CWI, P.O. Box 94079, 1090 GB Amsteram. jason@cwi.nl 2 epartment of Mathematics, Imperial College 180 Queen s Gate, Lonon, SW7 2BZ. s.reich@ic.ac.uk Abstract. Kelvin s circulation theorem an its implications for potential vorticity (PV) conservation are among the most funamental concepts in ieal flui ynamics. In this note, we iscuss the numerical treatment of these concepts with the Smoothe Particle Hyroynamics (SPH) an relate methos. We show that SPH satisfies an exact circulation theorem in an interpolate velocity fiel, an that, when appropriately interprete, this leas to statements of conservation of PV an generalize enstrophies. We also inicate some limitations where the analogy with ieal flui ynamics breaks own. AMS subject classification: 76M28. Key wors: geophysical flui ynamics, potential vorticity conserving methos, geometric methos, smoothe particle hyroynamics. 1 Introuction an Lagrangian Equations of Motion Large scale geophysical flows in the atmosphere an ocean are incompressible an nearly two-imensional. As a result, vorticity plays a central role in geophysical flui ynamics. In the Lagrangian flui escription, conservation of vorticity an circulation follow from the fact that flui particles are ientical an thus there is a great egree of freeom in labeling particles; that is, vorticity conservation follows from the particle nature of the flui ( 2). One might therefore expect that a Lagrangian, particle-base approach woul lea to goo vorticity conservation also in a computational setting. In this article, we consier the conservation properties of the popular Smoothe Particle Hyroynamics metho (SPH) [9, 14], as briefly outline in 3. We will show that, inee, the continuous velocity fiel that interpolates the SPH particle velocities exactly satisfies Kelvin s circulation theorem ( 4). In turn, Stoke s theorem implies that absolute vorticity is also exactly conserve. Recently, Monaghan [10] has suggeste that circulation is conserve approximately by SPH ue to a iscrete relabeling symmetry. We stress that in this article we prove exact conservation of circulation for a continuous interpolate velocity fiel, which implies the convergence of the iscrete integral introuce in [10]. The Hamiltonian Particle-Mesh metho evelope by the authors in [7] inherits conservation
2 2 J. FRANK AN S. REICH of circulation from SPH, an in 5 we provie a numerical illustration of this result using the HPM metho. Aitional conservation properties of importance in geophysical flui ynamics are mass an energy conservation. Mass conservation is intrinsic to the SPH formulation. It is also well-known that SPH can be erive from a variational principle (see, e.g., [14, 3, 10]); i.e. it can be given a Hamiltonian structure, for which symplectic integrators may be use to obtain excellent energy conservation [13]. For simplicity of exposition, we consier the two-imensional shallow water equations (SWEs): (1.1) (1.2) t u = c 0 x h, t h = h x u, 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. The results of this paper immeiately generalize to the rotating SWEs. This is briefly iscusse in 6. In the Lagrangian escription [12], the positions of all flui particles are given as a time epenent iffeomorphism from the flui label space A R 2 to R 2 : X = X(a, t), a = (a, b) A, X = (X, Y ) R 2. The labels are fixe for each particle, i.e. ta = 0, an the flui layer epth h is efine, as a function of the eterminant of the 2 2 Jacobian matrix through the relation (1.3) X a = (X, Y ) (a, b), h(x(a, t), t) X a = h o (a), where h o (a) is a time-inepenent function [12]. ifferentiation of (1.3) with respect to time yiels an expression that is equivalent to the continuity equation (1.2). Hence (1.3) an (1.2) are essentially equivalent statements. A natural choice for the labels a is given by a = X(0) = x which we assume from now on. Then the matrix X a is the ientity an h o the layer-epth at t = 0. Consier the integral ientity efining the layer-epth h at time t an Eulerian position x h(x, t) = h(x, t) δ(x X) X, where δ is the irac elta function. Using (1.3) we can pull this integral back to label space, arriving at the relation (1.4) h(x, t) = h o (a) δ(x X(a, t)) a,
3 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS 3 which can be taken as the efinition of the layer epth in the Eulerian reference frame. The motion of the SWEs is given as the stationary state δl = 0 of the action integral { 1 L = h o 2Ẋ Ẋ c } 0 (1.5) 2 h a t with respect to amissible variations δx(a, t), where the ot notation refers to time ifferentiation for fixe label, i.e. the material time erivative t. The equations of motion become (1.6) t X = u (1.7) t u = c 0 X h where h is efine by (1.3) or (1.4). Since the SWEs (1.6) (1.7) are erive from a Lagrangian variational form, they are Hamiltonian, conserving the energy E = 1 h o {u u + c 0 h} a = 1 p p a + c 0 h 2 x, 2 2 h o 2 p = h o u, an the symplectic two-form ω := (p X) a [11]. 2 Vorticity Conservation in Ieal Shallow Water Flows The action integral (1.5) epens on the labels a explicitly only through the layer epth h = h o X a 1. As a result, the action is invariant to any transformation of label space (i.e. relabeling ) that leaves the eterminant X a unchange (e.g. any ivergence-free relabeling (a, b )/ (a, b) = 1 will suffice.) By Noether s theorem, this particle relabeling symmetry implies a conserve quantity of the ynamics, which turns out to be the conservation of vorticity in its many forms as outline below. Note, however, that this symmetry with implie conservation law is very much tie into the particle nature of the flow in the Lagrangian escription; it follows essentially because particles are inistinguishable from one other. It is this observation which motivates us to consier the vorticity conservation properties of Lagrangian methos such as SPH in this article. For further iscussion of the particle relabeling symmetry an its relation to vorticity, we highly recommen Salmon s monograph [12]. Consier the vorticity ζ = x u. Using x t u = x u t + ( x u)( x u) + u x ( x u) = 0,
4 4 J. FRANK AN S. REICH it is easy to conclue that vorticity satisfies the continuity equation (2.1) t ζ = ζ x u. The ratio of ζ to h, i.e., q = ζ/h, is calle the potential vorticity (PV) [12]. The PV fiel q is materially conserve since, using (1.2) an (2.1), { t q = h 1 t ζ q } t h = 0. As a result, it follows that (1.2) an (2.1) are special cases of an infinite family of continuity equations (2.2) t [hf(q)] = x [hf(q)u], where f is an arbitrary (smooth) function of q. 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. X(s) = u(x(s)). t Kelvin s circulation theorem [12] states that (2.3) u X s s = 0. 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 V 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 = t V( (2.4) ( x u) x = 0. t R Because V 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 [5] to the right equality in (2.4), a secon appeal to the arbitrariness of V yiels the vorticity equation (2.1).
5 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS 5 3 Review of Smoothe Particle Hyroynamics The Smoothe Particle Hyroynamics (SPH) metho [9] is a Lagrangian metho for flui ynamics, in which the flui mass is istribute over a number of smooth, compactly supporte, particle-centere basis functions. A similar metho is the Finite Mass Metho [14]. 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 (3.1) 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 can be approximate by t g = x [gu], ḡ(x, t) = k γ k ψ(x X k (t)). Here {γ k } are constants etermine by the initial g(x) fiel an ψ is an appropriate basis function. Typically, SPH is implemente with raially symmetric basis functions, i.e. ψ(x) = Ψ( x ) an Ψ(r) is either a Gaussian, a compactly supporte raial basis function [4], or a raial spline [9]. Let us apply this iea to the layer-epth h, i.e., we introuce the approximation (3.2) h(x, t) = k m k ψ(x X k (t)) an assume that h(x, t) > 0. Then each particle contributes the fraction (3.3) ρ k (x, t) := m kψ(x X k (t)) h(x, t) 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 approximate ata from the particle locations to any x R 2. In particular, we efine an approximate Eulerian velocity fiel (3.4) ū(x, t) := k ρ k (x, t) u k (t) with layer epth flux ensity (inserting (3.3)) h(x, t) ū(x, t) = k m k ψ(x X k (t)) u k (t).
6 6 J. FRANK AN S. REICH Using (3.1), it is now easily verifie that (3.5) t h(x, t) + x [ h(x, t)ū(x, t)] = 0. 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 approximate velocity, i.e u k ū(x k ). We note that the moification suggeste in [8] 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 (3.6) t u k = c 0 x h(x, t) x=xk = c 0 m j Xk ψ(x k X j ). The equations (3.1), (3.2), an (3.6) comprise the stanar SPH approximation to the SWEs (1.1) (1.2). We introuce the canonical momenta p k = m k u k. The equations (3.1), (3.2), an (3.6) are now canonical with Hamiltonian (energy) (3.7) H = 1 2 k 1 p k 2 + c 0 m k 2 j m k m l ψ(x l X k ) an symplectic structure ω = k p k X k. A numerical time-stepping scheme is obtaine by noting that t X k = u k, l,k t u k = 0 can be solve exactly an that the implie time evolution of h(x, t) exactly satisfies (3.5). Similarly, equation (3.6) an t X k = 0 can also be integrate exactly since in this case h t = 0. A composition of these exact propagators leas to a symplectic time-stepping scheme [13] implying goo long-time energy conservation [2]. Bonet & Lok [3] have also iscusse conservation properties of SPH in a variational formulation. In particular they show conservation of linear an angular momentum, provie the basis function ψ is raially symmetric. These follow from the fact that the Hamiltonian (3.7) is invariant uner translations an rotations in the Lagrangian particle positions {X k } an subsequent application of Noether s theorem [11]. 1 More precisely, the formulation Eqn. (2.6) in [8] avocates particle avection in an approximate velocity fiel base on a generic kernel. Taking this kernel to be ψ yiels the avection fiel ū( ).
7 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS 7 4 Vorticity Conservation Properties of SPH Circulation is conserve, using SPH, for an interpolate velocity fiel u(x) efine as follows: at time t = 0, let u(x) be any initial velocity fiel satisfying u(x k ) = u k (0) at the particle locations X k. (For example, suppose the particle velocities at t = 0 are given as a continuous function.) efine natural labels a = X(0) = x an let the fiel of particle locations X(a) an particle velocities u(a) = u(x(a)) evolve uner the solution of the SPH flow ue to (3.6) accoring to the orinary ifferential equations (i.e. ecouple in label space) (4.1) t u = c 0 X h(x, t) = c0 m k X ψ( X X k (t) ), k t X = u. Note that for this velocity fiel it oes hol that u(x k ) = u k for all time t, in contrast to the approximate velocity fiel ū(x) of the continuity equation (3.4). Figure 4.1 illustrates the relationship between the velocity fiels u an ū. S u(x k ) = uk X(s) u(x(s)) X k u(x ) k Figure 4.1: A close curve X(s) avecte with the flow, illustrating the velocity fiels u an ū. The SPH particle at X k (support inicate by the otte line) with velocity u k remains on the curve throughout the integration. 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(s) = u(x(s)) accoring to X(s) = u(s). t
8 8 J. FRANK AN S. REICH We assume that X(s) an u(s) are sufficiently ifferentiable. circulation theorem (2.3) becomes (4.2) u X s s = 0. t Inee, we obtain ( ) ( ) u X s s = t t u X s s + u s t X s = c 0 X h Xs s + u u s s ( ) 1 = 2 (u u) s c 0 hs s = 0. Then Kelvin s It is important to notice that the circulation theorem above inuces a true constraint on the numerical solution, since any SPH particle that is initially locate on the loop S will remain on the loop as the integration procees, an furthermore, the particle velocity u k will be exactly interpolate by u(x k ). If we now give each particle insie S a label a an let V enote the area enclose by S an let R enote the image of V in x-space, then applying Stokes theorem to (4.2) yiels (4.3) t for which the left sie implies (4.4) V( X u) X a a = t R t {( X u) X a } = 0, ( x u) x = 0, since V is arbitrary. Let h(x(a)) enote the layer-epth approximation obtaine as the solution of the continuity equation (1.2) along the interpolate velocity fiel u(x(a)). Then h X a = h o an equation (4.4) implies conservation of the PV fiel q = ( X u)/h, i.e. q/t = 0. Furthermore, applying the transport theorem to the right equality of (4.3), an again noting that V is arbitrary, yiels a continuity equation for the absolute vorticity of the velocity fiel u: t ζ = x (ζu), ζ = x u, cf. (2.1). 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 can represent the loop {X(s)} by a sufficient number of particles { ˆX l } with associate velocities {û l } satisfying
9 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS 9 (4.1). Note that the particles { ˆX l } may be compute along with the SPH simulation as passive tracers with zero mass. The integral (4.2) is approximate by (4.5) u X s s û l ( ˆX l+1 ˆX l ). l One shoul observe that the variation of this integral in time converges to zero as the number of points iscretizing the loop is increase. Such a verification is inclue in 5 for the Hamiltonian Particle-Mesh metho. Monaghan [10] motivates the approximate conservation of (4.5), appealing to Noether s theorem for a iscrete relabeling (change of inex) of particles of equal mass. This reasoning is limite because Noether s theorem applies only to continuous symmetries. However, our result shows that the approximate conservation of (4.5) is an implication of a stronger result, namely the exact conservation of (4.2). We have seen that the PV fiel q = ( X u)/h is exactly conserve in SPH. It is also easy to verify that the prouct hf(q), where f is an arbitrary function of q, exactly satisfies the continuity equation (1.2) in the interpolate velocity fiel u(x). For iagnostic purposes, one can obtain a computable continuous approximation of hf(q) by again taking a sufficient number of passive tracer particles { ˆX l } with associate velocities {û l } satisfying (4.1). Following the iscussion for the layer-epth h in 3, we efine a continuous approximation hf(q)(x, t) = l α l ψ(x ˆX l (t). Given a particle mass m l an PV value q l for each tracer particle ˆX l, as etermine by the initial ata, the weights {α l } can be efine by α l = m l f(q l ). It is straightforwar to show that hf(q) satisfies a moifie continuity equation of type (3.5). Uner perioic bounary conitions, this continuity equation implies the exact conservation of the generalize enstrophies Q f = hf(q) x. Since this is also true for the split equations of motion use for the time-stepping, the overall space-time approximation conserves enstrophy. The number of tracer particles can be chosen to be quite large if, for example, the statistical/spectral properties of the SPH approximation to hq 2 are of interest. 5 Hamiltonian Particle-Mesh Metho an Numerical Verification The SPH equations of motion are equivalent to the simulation of a molecular flui with a softene repulsive pair potential given by the SPH basis function ψ(x) = Ψ( x ) [1]. In general, such flows ten to a statistically uniform state of isorer. In practice, therefore, SPH is use with some form of artificial viscosity, an this results in a monotone ecrease in circulation along any loop, an a monotone loss of energy.
10 10 J. FRANK AN S. REICH To improve the stability of SPH, we suggeste in [7] working with an average SWE (5.1) t X = u, (5.2) t u = c 0 X (A h), where A is some smoothing operator. These equations are still canonical with Hamiltonian H = 1 p p a + c 0 (5.3) h(a h) x 2 h o 2 an circulation preserving. Let {X k } be an initially equi-istribute set of points with an associate area A. Numerically the layer-epth h is now approximate by the singular measure h(x, t) = k m k δ(x X k (t)), m k = h o (a k ) A, a k = X k (0), which approximates the integral (1.4) in a weak sense. If, for example, A is chosen to be convolution with the SPH basis function ψ such that (A h) (x, t) = h(x, t) = k m k ψ(x X k (t)), then the layer-epth epenent part of the Hamiltonian (5.3) becomes c 0 h(a 2 h) x = c 0 h h x = c 0 m k m l ψ(x k X l ) 2 2 an the stanar SPH metho is recovere. However, typical SPH basis functions ψ o not provie enough smoothing, which results in the above mentione tenency to a state of isorer. The layer-epth approximation h satisfies the continuity equation (1.2) in a weak sense. Furthermore, the following integral version of (1.3) is easily shown: h X a a = h x = V R l,k k:x k R m k an m k h o (a) a. k:x k R V The set of particles over which the sum is performe is constant provie no particle enters or leaves the omain R which woul correspon to a singular X a an, hence, to a non-physical state.
11 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS 11 The Hamiltonian Particle-Mesh metho (HPM), introuce in [7], iffers from SPH primarily in the construction of the smoothe layer epth A h. Specifically, we efine a uniform gri with gri points x ij an gri spacing x. Let h SPH (x, t) be the SPH approximation to the layer epth (3.2) with the SPH basis function replace by a tensor prouct basis function ψ(x X k ) := φ( x X k )φ( y Y k ), where φ(r) is given by the cubic spline 2 3 r r3, r R 1 φ(r) = 6 (2 r)3, R < r 2R 0, r > 2R for R = x. efine the grie layer epth values h ij SPH (t) = h SPH (x ij, t) an let the matrix S = {Sij mn } enote the representation of a spatial averaging operator S over the given gri {x ij }. Since the cubic splines form a partition of unity on the gri, we can efine a continuous approximation of a smoothe/average layer-epth in space (A h) (x, t) = h HPM (x, t) = ij,mn h ij SPH (t) Smn ij ψ(x mn x). This approximation, use in the HPM metho, can be viewe as a spatial averaging over short wave-length isturbances in stanar SPH. One can also think of the HPM metho as an efficient implementation of the SPH metho for a globally supporte basis function ψ efine by A δ = ψ := S ψ, ψ a stanar SPH basis function. For a more etaile escription of the HPM metho, incluing its Hamiltonian structure, see [7]. The HPM metho conserves circulation using the same proof as for SPH in 4. The essential observation is that particles are avecte in a velocity fiel that exactly evolves in some continuous approximate layer epth. We have performe an experiment with HPM to verify the conservation of circulation. The flow moels the interaction of two positively oriente vortices in a rotating reference plane. We iscretize this flow using HPM on a gri with particles. We intentionally chose a fairly coarse iscretization for this problem to illustrate that the circulation theorem hols inepenent of the precision of the iscretization. Initially a circular loop of M evenly space particles of zero mass was place in the flow. The solution incluing the loop particles was evolve over time intervals of T = 3, T = 6 an T = 15 revolutions of the plane. The experiments
12 12 J. FRANK AN S. REICH were repeate, each time refining the iscretization, for M = 100, 200, 400, 800 an 1600 particles. The circulation integral was approximate using (cf. (4.5)) C M (t) = m û m ( ˆX m+1 ˆX m 1 ) Figure 5.1 shows the eformation of the loop at time T = 15, compute using M = 3200 particles. The loop is superimpose over a contour plot of potential vorticity, an its interior is shae. In Figure 5.2 we see secon orer convergence of (C M (t) C M (0))/C M (0) to zero as M increases. Convergence of this sum as M is implie by (4.2) Figure 5.1: Final eformation of the circular loop at time T = 15, using M = 3200 particles. 6 Concluing Remarks In this article, we have shown that the SPH metho with (3.6) satisfies a Kelvin circulation law. The results are base on the introuction of a continuous velocity fiel u(x) which interpolates the particle velocities u k for all time
13 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS T=15 (C M (T) C M (0))/C M (0) T=3 T= M Figure 5.2: Convergence of the iscretization of the circulation integral. an is avecte in the flow of the continuous SPH approximation (4.1). This velocity fiel conserves circulation (4.2) an, by Stokes theorem, absolute vorticity (4.3). Furthermore, we can formally efine a layer epth approximation h such that h satisfies the continuity equation (1.2) in the interpolate velocity fiel. efining the potential vorticity with respect to this layer epth yiels exact PV conservation. One can expect a particle metho to converge to the solutions of the average SWEs (5.1)-(5.2) for an appropriate smoothing operator A an in the limit of large particle numbers. The necessary amount of regularity an the impact of the smoothing operator A on the long term ynamics of the SWEs are not yet clear. 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
14 14 J. FRANK AN S. REICH an Kelvin s circulation theorem becomes 0 = ( u + f ) 0 t 2 X X s s = ( X u + f 0 ) X a a t V = ( x u + f 0 ) x. t R We wish to mention the Balance Particle-Mesh (BPM) metho of [6] which uses raial basis functions to approximate the absolute vorticity ω = x u + f 0. See [6] 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 [12]). 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 2m p k 2 + φ( X k X l ), l>k k where m is the atomic mass an φ(r) an interaction potential. We introuce the function ρ(x, t) = l an note that Newton s law is equivalent to φ( x X l (t) ), x R 3, 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 efining u = p/m, we obtain the circulation theorem u X s s = 0 t an, in two imensions, conservation of vorticity per control area, i.e., ( x u) x. 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) <. Here r o is chosen such that X i (t) X j (t) > r o for all t 0 an all i j. One shoul also note that, contrary to flui ynamics, the prouct ρ X a nee not to be approximately conserve an that X a can become singular.
15 CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS 15 Acknowlegements. J. Frank gratefully acknowleges partial funing by GM, Bonn. Partial financial support of S. Reich by EPSRC Grant GR/R09565/01 an by European Commission funing for the Research Training Network Mechanics an Symmetry in Europe is gratefully acknowlege. 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) J. Bonet & T.-S.L. Lok, Variational an momentum preservation aspects of smooth particle hyroynamic formulations, Comput. Methos Appl. Mech. Engrg., 180 (1999) 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, J. Sci. Comput., to appear. 7. J. Frank, G. Gottwal & S. Reich, The Hamiltonian particle-mesh metho, Meshfree Methos for Partial ifferential Equations, Lecture Notes in Computational Science an Engineering, Vol. 26, Springer, 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) J. J. Monaghan, SPH an the alpha turbulence moel, preprint, P. J. Morrison, Hamiltonian escription of the ieal flui, Rev. Moern Phys., 70 (1998) no. 2, 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)
Centrum voor Wiskunde en Informatica
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.
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
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 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 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 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 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 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 informationG4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.
G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether
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 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 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 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 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 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 informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
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 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 information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
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 informationFree rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012
Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate
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 informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationApplication of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate
Freun Publishing House Lt., International Journal of Nonlinear Sciences & Numerical Simulation, (9), -, 9 Application of the homotopy perturbation metho to a magneto-elastico-viscous flui along a semi-infinite
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 informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
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 informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
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 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 information1. At time t = 0, the wave function of a free particle moving in a one-dimension is given by, ψ(x,0) = N
Physics 15 Solution Set Winter 018 1. At time t = 0, the wave function of a free particle moving in a one-imension is given by, ψ(x,0) = N where N an k 0 are real positive constants. + e k /k 0 e ikx k,
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 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 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 informationTRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS
TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email
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 informationAn energy and potential enstrophy conserving numerical scheme for the multi-layer shallow water equations with complete Coriolis force
An energy an potential enstrophy conserving numerical scheme for the multi-layer shallow water equations with complete Coriolis force Anrew L. Stewart a, Paul J. Dellar b a Department of Atmospheric an
More informationFINITE POINTSET METHOD (FPM): A MESHFREE APPROACH FOR INCOMPRESSIBLE FLOW SIMULATIONS APPLIED TO COMPOSITE MATERIALS
FINITE POINTSET METHOD (FPM): A MESHFREE APPROACH FOR INCOMPRESSIBLE FLOW SIMULATIONS APPLIED TO COMPOSITE MATERIALS Alain Trameçon 1, Patrick e Luca 1, Christophe Binetruy 2 an Jörg Kuhnert 3 1 ESI-GROUP
More informationFluid Mechanics EBS 189a. Winter quarter, 4 units, CRN Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 Veihmeyer Hall.
Flui Mechanics EBS 189a. Winter quarter, 4 units, CRN 52984. Lecture TWRF 12:10-1:00, Chemistry 166; Office hours TH 2-3, WF 4-5; 221 eihmeyer Hall. Course Description: xioms of flui mechanics, flui statics,
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 information1.2 - Stress Tensor Marine Hydrodynamics Lecture 3
13.021 Marine Hyroynamics, Fall 2004 Lecture 3 Copyright c 2004 MIT - Department of Ocean Engineering, All rights reserve. 1.2 - Stress Tensor 13.021 Marine Hyroynamics Lecture 3 Stress Tensor τ ij:. The
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 information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationThe Ehrenfest Theorems
The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent
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 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 informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationNoether s theorem applied to classical electrodynamics
Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance
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 informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationCHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold
CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
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 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 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 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 informationMomentum and Energy. Chapter Conservation Principles
Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationThe influence of the equivalent hydraulic diameter on the pressure drop prediction of annular test section
IOP Conference Series: Materials Science an Engineering PAPER OPEN ACCESS The influence of the equivalent hyraulic iameter on the pressure rop preiction of annular test section To cite this article: A
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 informationRelation between the propagator matrix of geodesic deviation and the second-order derivatives of the characteristic function
Journal of Electromagnetic Waves an Applications 203 Vol. 27 No. 3 589 60 http://x.oi.org/0.080/0920507.203.808595 Relation between the propagator matrix of geoesic eviation an the secon-orer erivatives
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 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 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 informationMATHEMATICS BONUS FILES for faculty and students
MATHMATI BONU FIL for faculty an stuents http://www.onu.eu/~mcaragiu1/bonus_files.html RIVD: May 15, 9 PUBLIHD: May 5, 9 toffel 1 Maxwell s quations through the Major Vector Theorems Joshua toffel Department
More informationBasic Thermoelasticity
Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................
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 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 informationOptimal Variable-Structure Control Tracking of Spacecraft Maneuvers
Optimal Variable-Structure Control racking of Spacecraft Maneuvers John L. Crassiis 1 Srinivas R. Vaali F. Lanis Markley 3 Introuction In recent years, much effort has been evote to the close-loop esign
More informationDomain Decomposition and Model Reduction of Systems with Local Nonlinearities
Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
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 Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More information