Moment Models for kinetic chemotaxis/haptotaxis equations
|
|
- Imogen Stephens
- 5 years ago
- Views:
Transcription
1 Moment Models for kinetic chemotaxis/haptotaxis equations An overview Gregor Corbin March 20, 2017 Contents 1 Introduction The different model scales The kinetic cell model A moment model for the state variable 3 3 Moment models for the velocity The general procedure Some examples Discrete Ordinates The P 1 model in 1D The P N model in 3D The M 1 model in 1D The Kershaw closure Partial Moments models References 12 1 Introduction 1.1 The different model scales Microscale Equations of motion for individual cells/particles x i (t), v i (t), y i (t) many small particles Mesoscale Kinetic equation for the probability density function p(t, x, v, z) t ɛ 2 t, x ɛx Macroscale Advection-Diffusion equation for the density of particles ρ(t, x) = p Figure 1: The different model scales Here we consider the mesoscale only. The derivation of macroscopic equations from the kinetic levels is done in [2, 3] in the context of glioma tumor cell populations. 1
2 1.2 The kinetic cell model We consider the following kinetic model for glioma cell populations described in [2, 3]: t p + v x p }{{} + z ( }{{ Gp) } = L(p) }{{} straight movement change of state instant change of direction λ A p }{{} death or proliferation. (1) Here, p(t, x, v, z) is a distribution function for the number of cells at time t R +, position x D x R 3, with velocity v S 2 and internal state deviation z Z R. The internal state y in each cell follows an ODE ẏ = G(y) = k + (1 y)a(t, x) k y, (2) for some constant model parameters k +, k. The time and space dependent volume fraction of ECM (see [2]) is denoted A(t, x). Equation (2) has a steady state y, where G(y ) = 0. The variable z is then defined as the deviation from the steady state z = y y. With this substitution we obtain the G(z) G(z) = α(a(t, x))z + v x β(a(t, x)). (3) from equation (1) with α = k + A + k and β = k+ A α.(c.f [3]). The linear collision operator L corresponds to a local but instant change of velocity of the cells. It is given by L(p) = k z (v, v )p(v ) k z (v, v)p(v)dv, (4) S 2 = k z (v, v )p(v )dv λ(z)p(v). (5) S 2 The collision kernel k z (v, v ) has also to be modeled. It can be interpreted as the rate at which cell with velocity v change their velocity to v. In general it depends on the internal state of cells. It satisfies λ(z) := k z (v, v)dv, (6) S 2 for all v S 2. An important assumption on the collision kernel is, that it has an equilibrium distribution. Definition (Equilibrium distribution). A distribution F eq (v) that fulfills the detailed balance is an equilibrium distribution of the kernel k z (v, v ). k z (v, v )F eq (v ) = k z (v, v)f eq (v) v, v S 2 (7) The equilibrium could in general be different for every state z. Example (A simple kernel). where Q(v) is normalized, i.e. k z (v, v ) = λ(z)q(v), (8) S 2 Q(v)dv = 1. (9) The equilibrium distribution F eq (v) is obviously Q(v), independent of the state z. In the following we will only consider kernels of the above form (8). The model (8) comes from a haptotaxis setting, where cells orient themselves according to an underlying fibre distribution Q(v). With the notation := dv (10) the collision operator is V L = λ(z) [Q(v) p p]. (11) 2
3 2 A moment model for the state variable In this section we will derive the simplest possible moment model to discretize the state variable z of equation (1). This is done in the following steps: 1. define a basis b z = (b 0, b 1,...) of P (Z) 2. test (1) with all basis functions b i (z), i.e multiply by b i (z), and integrate ( t p + x (vp) + z (Gp)) b i (z)dz = L(p)b i (z)dz, i = 0, 1,... (12) Z 3. take only a finite number of basis functions b := (b 0,..., b N 1 ) to obtain a system of N equations for the moments f = (f 0, f 1,..., f N 1 ) := p(z)b z dz (13) As a basis of P (Z), we take the monomials b z = (1, z, z 2,...). We will use the simplest possible discrete moment model, with only a single basis function: How is this justified? Z Z b z = 1. (14) Assumption (Fast state dynamics). The change of state is fast, compared to the other effects. Therefore, most of the distribution p(z) is close to the steady state of G at z = 0. the solution is adequately described by its zeroth-order moment we neglect all second-order or higher moments f i, i 2 the distribution is negligibly small at the boundaries of Z: p(z min ) 0, p(z max ) 0 We expand the state dependent part of the kernel λ(z) in the basis Testing (1) by the finite basis b z = 1 yields t f 0 + x (vf 0 ) + Z λ(z) = λ i z i. (15) i=0 z (Gp)dz = i=0 λ i (Q(v) f i f i (v)) (16) where the integral over z (Gp) and the terms for i 2 vanish due to our assumptions: t f 0 + x (vf 0 ) = λ 0 (Q(v) f 0 f 0 (v)) + λ 1 (Q(v) f 1 f 1 (v)) (17) This equation for f 0 contains an additional unknown, namely the first moment f 1. We use a natural closure, meaning that we use the equation for f 1, which can be obtained by testing equation (1) by z: t f 1 + x (vf 1 ) Gpdz = λ 0 (Q(v) f 1 f 1 (v)). (18) Z Again using the assumption that most of the distribution is close to z = 0, we neglect the transport of the first moment f 1. Inserting the expression for G(z) yields: Solving this equation for f 1 gives αf 1 (v x β)f 0 = λ 0 (Q f 1 f 1 ) (19) f 1 = 1 α + λ 0 (λ 0 Q f 1 + (v x β)f 0 ). (20) Integrating equation (19) with respect to velocity gives an expression for f 1 f 1 = xβ α f 0v, (21) 3
4 where we used that Q = 1. Now we have an expression for the first-order moment f 1 that depends only on the zeroth-order moment f 0. Plugging equations (20) and (21) into the zeroth order equation (17) finally yields: t f 0 + x (vf 0 ) = λ H (Q f 0 f 0 ) + λ H (Q f 0 v f 0 v) (22) = L 0 (f 0 ) + L 1 (f 0 ) (23) where λ H := λ 0 and λ H := λ1 xβ λ 0+α. In contrast to the classical transport equation of radiation, there are now two collision terms. The first pushes the distribution f 0 towards the equilibrium Q, while the second does the same for f 0 v. In the following section we apply the method of moments again to discretize the velocity variable. 3 Moment models for the velocity 3.1 The general procedure All following considerations start from the simple moment model (23) from the previous section. We apply a similar procedure to discretize the velocity v. With the basis b for a finite-dimensional subspace of P (V ), the moments are defined as u = fb. (24) Then we take the moments of equation (23) w.r.t b t u + x v bf }{{} F = L0 (f) + L1 (f) } {{ } C (25) In general (but not always) F, C contain moments that are not components of u. We have to provide additional models for these moments. Usually, a reconstruction of the distribution ˆf(u) f is defined from the known! moments, with the condition that ˆfb = u. This is then used to compute approximations of F v b ˆf. The closed system then reads: t u + x v b ˆf = L0 ( ˆf) + L1 ( ˆf) (26) 3.2 Some examples Discrete Ordinates We choose b = (b 0,..., b N 1 ), b i = δ(v v i ) together with a linear ansatz for f: ˆf = i α iδ(v v i ). Then The coefficients α i follow from the moment constraints ˆfbi = α j b j b i The moment system becomes fb i = f(v i ) := f i = u i, (27) vb i f = v i f i, (28) j (29) = α i! = f i (30) t f i + x (v i f i ) = λ H (Q(v i ) f f i ) + λ H (Q(v i ) vf f i v i ). (31) The additional moments f, vf can be obtained from the ansatz: f ˆf = f i, (32) i vf v ˆf = v i f i. (33) i This is the discrete ordinates method(see also [7]). 4
5 3.2.2 The P 1 model in 1D For simplicity of notation, let us first consider a one-dimensional example, the so-called slab geometry, where x R and v V = [ 1, 1]. I. this setting the kinetic equation is t f + v x f = L 0 + L 1. (34) We choose the basis as b = (1, v), moments u = (u 0, u 1 ) = fb together with a modified linear ansatz ˆf = Q(v) α i b i = Q(v)(α 0 + α 1 v). (35) Remark. This ansatz incorporates the equilibrium distribution Q(v) of the kernel. It is thus able to represent the equilibrium exactly. Testing equation (23) by b gives t u 0 + x u 1 = λ H ( Q u 0 u 0 ) + λ H ( Q u 1 u 1 ), (36) t u 1 + x u 2 = λ H ( Qv u 0 u 1 ) + λ H ( Qv u 1 u 2 ), (37) There are two equations for the three unknowns u 0, u 1, u 2. We approximate the second moment u 2 from the ansatz u 2 (u 0, u 1 ) v 2 ˆf = Q(v)(α 0 + α 1 v)v 2, (38) where the multipliers α 0, α 1 have to be chosen, such that the ansatz has the coorect moments: condition leads to a linear system of equations:! ˆfb = u. This u 0 = α 0 Q + α 1 Qv, (39) u 1 = α 0 Qv + α 1 Qv 2. (40) In order to get the second moment u 2, we have to solve this system for α 0, α 1 and plug these values into the expression for u 2. In the special case of Q(v) = 3 2 v2, we have and therefore α 0 = u 0, (41) α 1 = 5 3 u 1, (42) u 2 (u 0, u 1 ) = 3 5 u 0 (43) The P N model in 3D Now we consider the general three-dimensional case, where x R 3 and v V = S 2. Also, we include higherorder functions in the basis. The following two choices are equivalent analytically (although not numerically): 1. b = (v k ) k for all multiindices k = (k 0, k 1, k 2 ) with k N. 2. b = (Yl m ) l,m, for l N, l m l, where Yl m are the spherical harmonics. Again, we choose a linear ansatz ˆf = Q α i b i = ˆαˆαˆα b. The multipliers α follow from the linear system of moment constraints:! ˆfb = u, (44) Qb b α! = u. (45) Note that in the special case of Q 1, and with the spherical harmonics, due to orthogonality, Qb b reduces to the identity matrix and therefore α = u. In the monomial basis, we have v b n f b n+1 f = u n+1, where the subscript n denotes all components with k = n. Only the highest-order moments v b N f have to be approximated by v b N ˆf = v b N (Qα b) = Qv b N b α. 5
6 3.2.4 The M 1 model in 1D Let us go back to the one-dimensional example x R, v V = [ 1, 1] with the first-order basis b = (1, v). But this time we choose a a nonlinear ansatz ˆf = Q exp( i α i b i ) = Q exp(α 0 + α 1 v) (46) Why should we do this? For Q 1, the ansatz minimizes the Maxwell-Boltzmann entropy under the moment constraints: Recall the H-Theorem for the following linear kinetic equation η(f) = f log f f, (47) ˆf = argmin{η(g) gb = u}. (48) g t f + v x f = λ 0 ( 1 f f). 4π (49) Definition (Entropy). For a convex function η, we call H(f) := η(f)dvdx (50) an entropy of the kinetic equation. X Theorem 1 (Boltzmann s H-Theorem). Any entropy for the linear kinetic equation (49) decreases monotonically with time d H(f) 0. (51) dt Thus, with this ansatz, the discretized system also dissipates the Maxwell-Boltzmann entropy. Additionally, it is guaranteed to be positive. For a detailed derivation of entropic moment closures, we refer to [1, 8]. Caution: In the context of haptotaxis, this ansatz is not rigourosly justified. There are several problems: V It is not known yet, if ˆf = Q exp(α b) even minimizes an entropy for (23). Also the H-Theorem does not hold for arbitrary entropies any more. The Maxwell-Boltzmann entropy is physically relevant for certain particle types. How sensible this choice is for biological cells we do not know. However: The M N model still guarantees positivity of the reconstruction, as opposed to the P N model. Intuitively, it makes sense to include Q into the ansatz to be able to reconstruct the equilibrium distribution correctly. Now, let us go back to compute the M 1 model. Writing out the moment constraints ˆfb = u leads to the following nonlinear system R u (α) = [ ] [ ] Q exp(α0 + α 1 v) u0! = 0 (52) Q exp(α 0 + α 1 v)v u 1 for α 0, α 1. Typically this is solved numerically, e.g. by Newton s method. Newton s method also needs the Jacobian of this expression: JR u = R [ ] u 1 v α = Q exp(α 0 + α 1 v) v v 2. (53) Note that R u (α) = 0 does not have a solution for all u. It turns out that it has a unique solution for the realizable moments. Definition (Realizable set). For a given basis b, the set of realizable moments is given by where P (V ) is the space of all non-negative distributions on V. R b = {u R n, f P (V ) : fb = u}, (54) The realizable set contains all those moments, for which there exists a non-negative distribution f on V. Conversely, if f is non-negative, its moments lie in the realizable set. Note that this definition only depends on the basis. In this case, the realizable set is simply ( ) 0 u0 R b = (55) u 1 u 0 a con as shown in figure 2. 6
7 Figure 2: The realizable set R b for b = (1, v) The Kershaw closure In the moment models that we have seen so far, the main desirable properties were representation of special distributions f (e.g. the equilibrium point) a closure that produces realizable higher-order moments The M N model naturally generates higher-order moments that are realizable. But this comes at the expense of having to solve a nonlinear system of equations for the multipliers α first. The Kershaw closure achieves these properties with a simpler ansatz that is therefore cheaper to compute. The reconstruction ˆf is done by interpolation between certain distributions f, that should be represented. To keep the higher-order moments realizable, explicit information on the higher-order realizable set is needed. For more theory on realizability and Kersahw closures, see [6, 9, 10] First-order Kershaw closure in 1D The first-order monomial basis in 1D is b = (1, v) with moments. u = (u 0, u 1 ). The Kershaw closure needs the realizable set R b + for the higher-order basis b+ = (1, v, v 2 ). The higher-order realizable set is 0 u 0 R b + u 1 u 0, (56) u 2 1 u 0 u 2 which is shown in figure 3. We get a realizable closure for u 2 by interpolation between the upper and lower bound on u 2 : u 2,K u 0 = ξ 1 + (1 ξ) ( u1 u 0 ) 2 (57) How to choose the interpolation coefficient ξ? At the equilibrium, we have u 0 = Q = 1 and u 1 = Qv = 0. Since we want to reproduce moments of the equilibrium, we need Qv 2! = u 2,K (58) Qv 2 ( ) 2 u1 = ξ + (1 ξ) (59) u 0 Qv 2 = ξ + (1 ξ) Qv 2 (60) Qv 2 Qv 2 ξ = 1 Qv 2 (61) If Q is symmetric around 0, then Qv = 0 and ξ = Qv 2. For example, using Q 1 yields ξ =
8 Figure 3: The realizable set R b for b = (1, v, v 2 ). The Kershaw model interpolates the second moment u ( ) 2 2 between the upper boundary u2 u 0 = 1 and the lower boundary u2 u u 0 = 1 u 0 of the realizable set First-order Kershaw closure in 2D The velocity space is V = S 2. The first-order monomial basis is b = (1, v) = (1, v 0, v 1 ), with moments u = (ρ, q 0, q 1 ). We need second-order moments P := v vf to close the system. Define the second order basis b + = (1, v, v v) = (1, v 0, v 1, v0, 2 v 0 v 1, v1) 2 and the corresponding moments u + = (ρ, q, P ) = (ρ, q 0, q 1, P 00, P 01, P 11 ). We want to construct P such that u + R b +. Additionally,the model should also incorporate our equilibrium Q. We calculate P for two extreme cases explicitly. 1. For the equilibrium solution f eq = ρq, we have q = 0. The second moments of this distribution are P eq = ρ v vq =: ρd F. 2. In the so-called free-streaming case the distribution f δ = ρδ(v q q q = ρ. The second moments are P δ = ρ v vδ(v q q ) = ρ q q q. 2 ) is a beam in direction q q. Therefore Then we interpolate between those two cases, s.t. realizability is preserved. Figure 4 shows the projection of R b + onto q0 ρ, q1 ρ and the interpolation procedure. The ansatz for P becomes ( P (ρ, q) = ρ α(ρ, q)d F + (1 α(ρ, q)) q q ) q 2 (62) Now, we have to choose α(ρ, q) s.t. u + is realizable. For example α = 1 q 2 ρ two extreme cases, this ansatz gives the correct second moments. achieves this. Note that for the Discussion The Kershaw closure lets us compute higher-order moments explicitly from the known moments, while preserving realizability and reconstructing moments of special distributions exactly. This means it is much faster to compute than the M N model. However, the generalization beyond a first-order model is hard, because the realizable set has to be known exactly. Already the second-order set is five-dimensional. For a second-order Kershaw closure, one needs to compute the third-order realizability domain in nine dimensions. Also, point-wise directional information on the reconstruction ˆf is needed for a kinetic upwind scheme Partial Moments models A simple first-order finite-volume method with forward Euler time discretization in 1D is [ ] u n+1 i = u n i + t vĝ i+ 1 b vĝ 2 i 1 b +... (63) 2 A kinetic scheme uses upwind point-wise for velocity in the numerical flux ĝ: ĝ i+ 1 2 = { ˆf(xi ) v > 0 ˆf(x i+1 ) v < 0. (64) 8
9 Figure 4: Kershaw closure on the realizable set R b for b = (1, v x, v y ): Shown are the normalized moment components qx ρ, qy ρ This cannot be done in the Kershaw model, since we do not have ˆf, but only the moments. Idea: Define a non-overlapping decomposition of V into patches. On each patch, use a low-order moment model. Advantage: Closure problem decouples. We can use a kinetic scheme, if patches do not cross quadrant boundaries. We can add more degrees of freedom relatively easy. [4, 5] 9
10 (a) flux at the right face (b) upwind patches at the right face (c) flux at the top face (d) upwind patches at the top face Figure 5: Upwind scheme for quarter-moments: For the right face of cell i, j, moments associated with positive velocity in x-direction(u ++, u + ) are taken from the left cell i, j. Moments associated with negative velocity in x-direction are taken from the right cell i + 1, j. Analogously in y-direction. 10
11 Figure 6: The unit sphere with quadrants 11
12 References [1] B. Dubroca and J. L. Feugeas. Entropic moment closure hierarchy for the radiative transfer equation. C. R. Acad. Sci. Paris Ser. I, pages , [2] C. Engwer, T. Hillen, M. Knappitsch, and C. Surulescu. Glioma follow white matter tracts: a multiscale dti-based model. Mathematical Biology, page , [3] Christian Engwer, Alexander Hunt, and Christina Surulescu. Effective equations for anisotropic glioma spread with proliferation: a multiscale approach and comparisons with previous settings. Mathematical Medicine and Biology, 33(4): , [4] M. Frank. Partial Moment Models for Radiative Transfer. Shaker-Verlag, Aachen, [5] Martin Frank, Bruno Dubroca, and Axel Klar. Partial moment entropy approximation to radiative heat transfer. Journal of Computational Physics, 218(1):1 18, [6] DS Kershaw. Flux limiting nature s own way. Lawrence Livermore National Laboratory, UCRL-78378, [7] E. W. Larsen and J.E.Morel. Advances in discrete-ordinates methodology. Nuclear computational Science, pages 1 84, [8] D. Levermore. Moment closure hierarchies for kinetic theories. Journal of Statistical Physics, [9] Philipp Monreal and Martin Frank. Moment realizability and kershaw closures in radiative transfer. Technical report, Lehr-und Forschungsgebiet Simulation in der Kerntechnik, [10] Florian Schneider. Kershaw closures for linear transport equations in slab geometry i: model derivation. Journal of Computational Physics, 322: ,
Fluid Dynamics from Kinetic Equations
Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationOptimal Prediction for Radiative Transfer: A New Perspective on Moment Closure
Optimal Prediction for Radiative Transfer: A New Perspective on Moment Closure Benjamin Seibold MIT Applied Mathematics Mar 02 nd, 2009 Collaborator Martin Frank (TU Kaiserslautern) Partial Support NSF
More informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More informationThe Boltzmann Equation and Its Applications
Carlo Cercignani The Boltzmann Equation and Its Applications With 42 Illustrations Springer-Verlag New York Berlin Heidelberg London Paris Tokyo CONTENTS PREFACE vii I. BASIC PRINCIPLES OF THE KINETIC
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationCORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire
CORBIS: Code Raréfié Bidimensionnel Implicite Stationnaire main ingredients: [LM (M3AS 00, JCP 00)] plane flow: D BGK Model conservative and entropic velocity discretization space discretization: finite
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
More informationAn electrokinetic LB based model for ion transport and macromolecular electrophoresis
An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained
More informationNumerical methods for kinetic equations
Numerical methods for kinetic equations Lecture 6: fluid-kinetic coupling and hybrid methods Lorenzo Pareschi Department of Mathematics and Computer Science University of Ferrara, Italy http://www.lorenzopareschi.com
More informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
More informationScattering of an α Particle by a Radioactive Nucleus
EJTP 3, No. 1 (6) 93 33 Electronic Journal of Theoretical Physics Scattering of an α Particle by a Radioactive Nucleus E. Majorana Written 198 published 6 Abstract: In the following we reproduce, translated
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationBERNOULLI EQUATION. The motion of a fluid is usually extremely complex.
BERNOULLI EQUATION The motion of a fluid is usually extremely complex. The study of a fluid at rest, or in relative equilibrium, was simplified by the absence of shear stress, but when a fluid flows over
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 7 (2/26/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects
.54 Neutron Interactions and Applications (Spring 004) Chapter 7 (/6/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects References -- J. R. Lamarsh, Introduction to Nuclear Reactor
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationAn asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations
An asymptotic preserving unified gas kinetic scheme for the grey radiative transfer equations Institute of Applied Physics and Computational Mathematics, Beijing NUS, Singapore, March 2-6, 2015 (joint
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationHierarchical Modeling of Complicated Systems
Hierarchical Modeling of Complicated Systems C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park, MD lvrmr@math.umd.edu presented
More informationLecture 6 Gas Kinetic Theory and Boltzmann Equation
GIAN Course on Rarefied & Microscale Gases and Viscoelastic Fluids: a Unified Framework Lecture 6 Gas Kinetic Theory and Boltzmann Equation Feb. 23 rd ~ March 2 nd, 2017 R. S. Myong Gyeongsang National
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationHydrodynamic Limits for the Boltzmann Equation
Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationSimplified Hyperbolic Moment Equations
Simplified Hyperbolic Moment Equations Julian Koellermeier and Manuel Torrilhon Abstract Hyperbolicity is a necessary property of model equations for the solution of the BGK equation to achieve stable
More informationQuantum Hydrodynamics models derived from the entropy principle
1 Quantum Hydrodynamics models derived from the entropy principle P. Degond (1), Ch. Ringhofer (2) (1) MIP, CNRS and Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse cedex, France degond@mip.ups-tlse.fr
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationGyrokinetic simulations of magnetic fusion plasmas
Gyrokinetic simulations of magnetic fusion plasmas Tutorial 2 Virginie Grandgirard CEA/DSM/IRFM, Association Euratom-CEA, Cadarache, 13108 St Paul-lez-Durance, France. email: virginie.grandgirard@cea.fr
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationSUPPLEMENTARY INFORMATION
doi: 1.138/nature5677 An experimental test of non-local realism Simon Gröblacher, 1, Tomasz Paterek, 3, 4 Rainer Kaltenbaek, 1 Časlav Brukner, 1, Marek Żukowski,3, 1 Markus Aspelmeyer, 1, and Anton Zeilinger
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationTHE GRADIENT PROJECTION ALGORITHM FOR ORTHOGONAL ROTATION. 2 The gradient projection algorithm
THE GRADIENT PROJECTION ALGORITHM FOR ORTHOGONAL ROTATION 1 The problem Let M be the manifold of all k by m column-wise orthonormal matrices and let f be a function defined on arbitrary k by m matrices.
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationDomain of definition of Levermore s five moment system
Domain of definition of Levermore s five moment system Michael Junk Institut für Techno und Wirtschaftsmathematik Kaiserslautern, Germany Abstract The simplest system in Levermore s moment hierarchy involving
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationPhysics GRE: Electromagnetism. G. J. Loges 1. University of Rochester Dept. of Physics & Astronomy. xkcd.com/567/
Physics GRE: Electromagnetism G. J. Loges University of Rochester Dept. of Physics & stronomy xkcd.com/567/ c Gregory Loges, 206 Contents Electrostatics 2 Magnetostatics 2 3 Method of Images 3 4 Lorentz
More informationKinetic models of Maxwell type. A brief history Part I
. A brief history Part I Department of Mathematics University of Pavia, Italy Porto Ercole, June 8-10 2008 Summer School METHODS AND MODELS OF KINETIC THEORY Outline 1 Introduction Wild result The central
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 7 (2/26/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects
.54 Neutron Interactions and Applications (Spring 004) Chapter 7 (/6/04) Neutron Elastic Scattering - Thermal Motion and Chemical Binding Effects References -- J. R. Lamarsh, Introduction to Nuclear Reactor
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationj=1 r 1 x 1 x n. r m r j (x) r j r j (x) r j (x). r j x k
Maria Cameron Nonlinear Least Squares Problem The nonlinear least squares problem arises when one needs to find optimal set of parameters for a nonlinear model given a large set of data The variables x,,
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationDerivation of the GENERIC form of nonequilibrium thermodynamics from a statistical optimization principle
Derivation of the GENERIC form of nonequilibrium thermodynamics from a statistical optimization principle Bruce Turkington Univ. of Massachusetts Amherst An optimization principle for deriving nonequilibrium
More informationParticle in Cell method
Particle in Cell method Birdsall and Langdon: Plasma Physics via Computer Simulation Dawson: Particle simulation of plasmas Hockney and Eastwood: Computer Simulations using Particles we start with an electrostatic
More informationPhysical models for plasmas II
Physical models for plasmas II Dr. L. Conde Dr. José M. Donoso Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Physical models,... Plasma Kinetic Theory
More informationKinetic theory of the ideal gas
Appendix H Kinetic theory of the ideal gas This Appendix contains sketchy notes, summarizing the main results of elementary kinetic theory. The students who are not familiar with these topics should refer
More informationTitle of communication, titles not fitting in one line will break automatically
Title of communication titles not fitting in one line will break automatically First Author Second Author 2 Department University City Country 2 Other Institute City Country Abstract If you want to add
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationBalázs Gyenis. SZTE TTIK Elméleti Fizika Tanszék,
Az egyensúly felé törekvés és a valószínűség interpretációja Balázs Gyenis Institute of Philosophy Hungarian Academy of Sciences, RCH SZTE TTIK Elméleti Fizika Tanszék, 2018.03.22. Introduction When two
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationLattice Boltzmann Method for Moving Boundaries
Lattice Boltzmann Method for Moving Boundaries Hans Groot March 18, 2009 Outline 1 Introduction 2 Moving Boundary Conditions 3 Cylinder in Transient Couette Flow 4 Collision-Advection Process for Moving
More informationPhysics Dec Time Independent Solutions of the Diffusion Equation
Physics 301 10-Dec-2004 33-1 Time Independent Solutions of the Diffusion Equation In some cases we ll be interested in the time independent solution of the diffusion equation Why would be interested in
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationChapter 18 Thermal Properties of Matter
Chapter 18 Thermal Properties of Matter In this section we define the thermodynamic state variables and their relationship to each other, called the equation of state. The system of interest (most of the
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. #12 Fundamentals of Discretization: Finite Volume Method
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationORDINARY DIFFERENTIAL EQUATION: Introduction and First-Order Equations. David Levermore Department of Mathematics University of Maryland
ORDINARY DIFFERENTIAL EQUATION: Introduction and First-Order Equations David Levermore Department of Mathematics University of Maryland 7 September 2009 Because the presentation of this material in class
More informationASTR 610 Theory of Galaxy Formation Lecture 4: Newtonian Perturbation Theory I. Linearized Fluid Equations
ASTR 610 Theory of Galaxy Formation Lecture 4: Newtonian Perturbation Theory I. Linearized Fluid Equations Frank van den Bosch Yale University, spring 2017 Structure Formation: The Linear Regime Thus far
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationORDINARY DIFFERENTIAL EQUATIONS: Introduction and First-Order Equations. David Levermore Department of Mathematics University of Maryland
ORDINARY DIFFERENTIAL EQUATIONS: Introduction and First-Order Equations David Levermore Department of Mathematics University of Maryland 1 February 2011 Because the presentation of this material in class
More informationA HYPERBOLIC RELAXATION MODEL FOR PRODUCT FLOW IN COMPLEX PRODUCTION NETWORKS. Ali Unver, Christian Ringhofer and Dieter Armbruster
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS Supplement 2009 pp. 790 799 A HYPERBOLIC RELAXATION MODEL FOR PRODUCT FLOW IN COMPLEX PRODUCTION NETWORKS Ali Unver, Christian Ringhofer
More informationTropical Cyclones: Steady State Physics
Tropical Cyclones: Steady State Physics Energy Production Carnot Theorem: Maximum efficiency results from a particular energy cycle: Isothermal expansion Adiabatic expansion Isothermal compression Adiabatic
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationReview of Electrostatics. Define the gradient operation on a field F = F(x, y, z) by;
Review of Electrostatics 1 Gradient Define the gradient operation on a field F = F(x, y, z) by; F = ˆx F x + ŷ F y + ẑ F z This operation forms a vector as may be shown by its transformation properties
More informationExponential methods for kinetic equations
Exponential methods for kinetic equations Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara, Italy http://utenti.unife.it/lorenzo.pareschi/ lorenzo.pareschi@unife.it Joint research
More informationLectures notes on Boltzmann s equation
Lectures notes on Boltzmann s equation Simone Calogero 1 Introduction Kinetic theory describes the statistical evolution in phase-space 1 of systems composed by a large number of particles (of order 1
More informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
More informationSpectral and Bulk Mass-Flux Convective Parameterizations
Spectral and Bulk Mass-Flux Convective Parameterizations Bob Plant Department of Meteorology, University of Reading COST ES0905 workshop on Concepts for Convective Parameterizations in Large-Scale Models:
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationInterpolation and Approximation
Interpolation and Approximation The Basic Problem: Approximate a continuous function f(x), by a polynomial p(x), over [a, b]. f(x) may only be known in tabular form. f(x) may be expensive to compute. Definition:
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationDirectional Field. Xiao-Ming Fu
Directional Field Xiao-Ming Fu Outlines Introduction Discretization Representation Objectives and Constraints Outlines Introduction Discretization Representation Objectives and Constraints Definition Spatially-varying
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationBrief Review of the R-Matrix Theory
Brief Review of the R-Matrix Theory L. C. Leal Introduction Resonance theory deals with the description of the nucleon-nucleus interaction and aims at the prediction of the experimental structure of cross
More informationReview of Electrostatics
Review of Electrostatics 1 Gradient Define the gradient operation on a field F = F(x, y, z) by; F = ˆx F x + ŷ F y + ẑ F z This operation forms a vector as may be shown by its transformation properties
More informationHamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization
NNP2017 11 th July 2017 Lawrence University Hamiltonian and Non-Hamiltonian Reductions of Charged Particle Dynamics: Diffusion and Self-Organization N. Sato and Z. Yoshida Graduate School of Frontier Sciences
More information