Review of last lecture I
|
|
- Ethelbert Strickland
- 6 years ago
- Views:
Transcription
1 Review of last lecture I t f + v x f = Q γ (f,f), (1) where f (x,v,t) is a nonegative function that represents the density of particles in position x at time t with velocity v. In (1) Q γ is the so-called granular collision operator, which is a quadratic integral operator descrbing the change in the density function due to creation and annihilation of particles in binary collisions: { v Q γ (f,f) = β (θ) w 1 } R R + J f(v )f(w ) v w f(v)f(w) dw Where J is the jacobian of the transformation (v,w ) (v,w). Here we consider (v,w ) (v,w), and v = 1 2 ( v + w ) ( v w ) h
2 Review of last lecture II w = 1 2 ( v + w ) 1 2 ( v w ) h h = 0 means the two particles become one after the collision and share the same velocity. And h = 1 corresponds the elasticity case. And we give the expression of h In this case J = h. h = θ
3 Splitting algorithm I The main part of our analysis will be devoted to the space homogeneous equation t f = Q γ (f,f), (2) It is well known, in fact, that by the standard splitting algorithm we may consider seperately the transport t f + v x f = 0, (3) and the relaxation given by (2). The overall accuracy of this simple splitting is first order in time. A second order generalization of this method is given by Strang s splitting [54]. Further discretization techniques can be found in [44].
4 Splitting algorithm II Since the main difficulties are represented by the discretization of the collision operator, it is clear that after the splitting the most difficult part relies in the approximation of the collision step (2). For the operator Q γ, there is conservation both of the mass ρ(x,t) = f(x,v,t)dv, and of the momentum u(x,t) = while the energy is dissipated. E (x,t) = R R R vf(x,v,t)dv v 2 f(x,v,t)dv
5 A spectral method I In this section we derive the Fourier spectral method for the homogeneous Boltzmann equation (7.2)-(7.9) following [43,46]. To this aim we define with, the inner product in L 1 (R) and consider the weak form of the equation t f,ϕ = Q γ (f,f),ϕ = β (θ) v w [ ϕ ( v ) ϕ(v) ] R + R 2 f(v,t)f(w,t)dwdθ, for t > 0 and all test functions ϕ. Here (v,w) (v,w ), that is, we have v = 1 2 (v + w) + 1 (v w) h (4) 2 w = 1 2 (v + w) 1 (v w) h (5) 2
6 A spectral method II A simple change of variables permits to write Q γ (f,f),ϕ = β (θ) q [ ϕ(v + q + ) ϕ(v) ] R + R 2 f(v,t)f(v + q,t)dqdθdv, where q = w v is the relative velocity, and the vectors q + and q that parametrize the post-collisional velocities are given by q + = q ( ), q = q ( 1 1 ) θ θ We point out that the possibility to integrate the collision operator over the relative velocity is essential in the derivation of the method. We consider now an initial density function f 0 (v,t) with compact support Supp (f 0 (v,t)) [ R,R]. The solution to (2)
7 A spectral method III has compact support for later times. In fact by (4,5), if v, w R then v 1 2 v (1 + h) + 1 w (1 h) R 2 and similary we get w R. In addition q, q +, q 2R, thus we have the following Lemma 1 If the function f(v,t) is such that Supp(f(v,t)) [ R,R] then i) Supp(Q γ ( f,f)(v,t)) [ R,R], ii) Q γ (f,f),ϕ = β (θ) q [ ϕ(v + q + ) ϕ(v) ] R + v R q 2R f(v,t)f(v + q,t)dqdθdv,
8 A spectral method IV with v + q +, v + q [ 3R,3R]. Remark: The previous result shows that for compactly supported functions f in order to evaluate Q γ (f,f) by a psectral method without aliasing error we can consider the density function f restricted on the interval [ 2R,2R], and extend it by periodicity to a periodic function on [ 2R,2R]. Spectral projection of the equation To simplify the notation let us take 2R = π. The approximate function f N is represented as the truncated Fourier series f N (v) = N k= N ˆf k = 1 f(v)e ikv dv. 2π [ π,π] ˆf k e ikv, (6)
9 A spectral method V A Fourier-Galerkin method is obtained by considering the projection of the homogeneous Boltzmann equation on the space of trigonometric polynomials of degree N. Hence, taking f = f N and ϕ = e ikv for k = N,,N we have [ t f N Q γ (f N,f N )]e ikv dv = 0. (7) [ π,π] By substituting expression (6) into (7) we get a set of ordinary differential equations for the Fourier coefficients t ˆfk = N l,m= N l+m=k ˆf l ˆfm ˆβ (l,m) δk 2 ˆfk, k = N,,N,
10 A spectral method VI where the Boltzmann kernel modes ˆβ (l,m) are given by ˆβ (l,m) = dθβ (θ) q [ cos ( lq + mq ) cos (lq) ] dq. R + q π In fact, by evaluating (7.14) for ϕ = e ikv and f = f N, one obtains [ ] ˆβ (l,m) = dθβ (θ) q e ikq+ 1 e ilq dq, (8) R + q π and (8) follows by using the parities of the trignometric functions. Note that (8) is a real quantity completely independent of the argument v, dpending on just the particular kernel structure. This property is strictly related to the use of a Fourier spectral method. Other spectral methods may be developed, however they do not lead to this simplification.
11 A spectral method VII In practice all the information characterizing the kinetic eqation is now contained in the kernel modes. Clearly, these quantities can be computed in advance and then stored in a two-dimensional matrix of size 2N. Thanks to symmetry considerations the effective number of kernel modes that need to be computed and stored for the implementation of the method is reduced in practice since ˆβ (l,m) = ˆβ ( l, m). Note that an analytic expression for ˆβ(l,m) can be readily comuted as [ ˆβ (l,m) = β (θ)π 2 2Sinc (p) Sinc ( p) 2 2 R + 2Sinc (l) +Sinc ( ) l 2 ]dθ 2
12 A spectral method VIII where p = ((l m) + (l + m)h) /2 and Sinc(x) = sin(πx) /(πx). Finally we can rewrite scheme (7.18) as t ˆfk = N m= N ˆf k m ˆfm ˆβ (k m,m), k = N,,N, In the previous expression we assume that thefourier coefficents are extended to zero fo k > N. The eavaluation of (7.22) requires exactly O ( N 2) opeartions which is smaller than the cost of a standard method based on N parameters for f in the velocity space since we gain the integration over the variable Thus the straightforwrad evaluation of (7.22) is slightly less expensive than a usual discrete-velocity algorithm.
13 Monte Carlo Method I Prof. Pareschi s ppt PartIII.pdf.
14 Homework I where t f = Q(f,f) Q(f,f) = R R + β (θ) { v w 1 } J f(v )f(w ) v w f(v)f(w) dwd β (θ) = exp( θ), h = 1 θ 1 + θ, J = h v = 1 2 (v + w) + 1 (v w) h 2 w = 1 2 (v + w) + 1 (v w) h 2
15 Homework II ( f (v,0) = exp (2v 2) 2) ( + exp (2v + 2) 2). t = 0.5, t = 8 f (v) v [ π,π]. Monte Carlo
16 Monte Carlo Method for this hw I tf = Q(f,f) can be written as t f = Q(f,f) = Q + (f,f) µf where Q + (f,f) = β (θ) v w 1 R R + J f(v )f(w )dw = β (θ)dθ v w f(v )f(w )dw R + R π = β (θ)dθ v w f(v )f(w )dw R + π π µ = β (θ)dθ v w f(w)dw (9) R + π
17 Monte Carlo Method for this hw II R v w f(w)dw O ( N 2) Using Euler s method to discretize the Boltzmann equation yields f n+1 = (1 µ t)f n + µ t Q+ (f n,f n ) µ where µ is given by (9). But µ is a integral and directly computing will result into O ( N 2) operation for all the v. v [ π,π]. v w 2π = κ. And we can give an estimation of µ π µ = β (θ)dθ v w f(w)dw R + π π β (θ)dθκ f(w)dw R + π
18 Monte Carlo Method for this hw III = κρ β (θ)dθ R + And we use the notation = κρ β (θ)dθ R + Now we can modify the Algorithm [Nanbu-Babovsky for VHS molecules] on page 21 of PartIII.pdf (Prof. Pareschi) to solve our problem. 1. Compute the initial velocity of the particles, { v 0 i,i = 1,,N } by sampling them from the initial density f 0 (v) 2. set time counter t c = 0 3. for n=1 to n tot - Set N c = Iround ( N 2 t)
19 Monte Carlo Method for this hw IV - select N c random pairs (i,j) uniformly within all possible N (N 1) /2 pairs for each pair, find a θ by sampling from the density β (θ)/ R β (θ)dθ. + - Compute the relative velocity B ij = v i w j - get a random number r by sampling a uniformly distributed density function on the interval [0,κ]. - if r < B ij perform the collision between i and j, and compute v i and according to the collision law v j v = 1 2 (v + w) + 1 (v w)h 2 w = 1 2 (v + w) 1 (v w)h 2
20 Monte Carlo Method for this hw V set vi n+1 = v i, and vn+1 j = v j - set vi n+1 = vi n for all the particles that did not collide. end for
Monte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 2: Monte Carlo simulation methods Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
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 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 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 informationA fast spectral method for the inelastic Boltzmann collision operator and application to heated granular gases
A fast spectral method for the inelastic Boltzmann collision operator and application to heated granular gases Jingwei Hu a, Zheng Ma a a Department of Mathematics, Purdue University 150 N. University
More informationON THE SELF-SIMILAR ASYMPTOTICS FOR GENERALIZED NON-LINEAR KINETIC MAXWELL MODELS
ON THE SELF-SIMILAR ASYMPTOTICS FOR GENERALIZED NON-LINEAR KINETIC MAXWELL MODELS A.V. BOBYLEV( ), C. CERCIGNANI( ), I.M. GAMBA( ) Abstract. Maxwell models for nonlinear kinetic equations have many applications
More informationExponential tail behavior for solutions to the homogeneous Boltzmann equation
Exponential tail behavior for solutions to the homogeneous Boltzmann equation Maja Tasković The University of Texas at Austin Young Researchers Workshop: Kinetic theory with applications in physical sciences
More informationHybrid and Moment Guided Monte Carlo Methods for Kinetic Equations
Hybrid and Moment Guided Monte Carlo Methods for Kinetic Equations Giacomo Dimarco Institut des Mathématiques de Toulouse Université de Toulouse France http://perso.math.univ-toulouse.fr/dimarco giacomo.dimarco@math.univ-toulouse.fr
More informationA conservative spectral method for the Boltzmann equation with anisotropic scattering and the grazing collisions limit
Work supported by NSF grants DMS-636586 and DMS-1217254 Jeff Haack (UT Austin) Conservative Spectral Boltzmann 7 Jun 213 1 / 21 A conservative spectral method for the Boltzmann equation with anisotropic
More informationc 2012 International Press FAST COMMUNICATION
COMMU. MATH. SCI. Vol. 1, o. 3, pp. 989 999 c 1 International Press FAST COMMUICATIO A FAST SPECTRAL ALGORITHM FOR THE QUATUM BOLTZMA COLLISIO OPERATOR JIGWEI HU AD LEXIG YIG Abstract. This paper introduces
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 informationExponential moments for the homogeneous Kac equation
Exponential moments for the homogeneous Kac equation Maja Tasković University of Pennsylvania Young Researchers Workshop: Stochastic and deterministic methods in kinetic theory, Duke University, November
More informationMonte Carlo method with negative particles
Monte Carlo method with negative particles Bokai Yan Joint work with Russel Caflisch Department of Mathematics, UCLA Bokai Yan (UCLA) Monte Carlo method with negative particles 1/ 2 The long range Coulomb
More informationYAN GUO, JUHI JANG, AND NING JIANG
LOCAL HILBERT EXPANSION FOR THE BOLTZMANN EQUATION YAN GUO, JUHI JANG, AND NING JIANG Abstract. We revisit the classical ork of Caflisch [C] for compressible Euler limit of the Boltzmann equation. By using
More informationPhysics 212: Statistical mechanics II Lecture IV
Physics 22: Statistical mechanics II Lecture IV Our program for kinetic theory in the last lecture and this lecture can be expressed in the following series of steps, from most exact and general to most
More informationA Conservative Discontinuous Galerkin Scheme with O(N 2 ) Operations in Computing Boltzmann Collision Weight Matrix
A Conservative Discontinuous Galerkin Scheme with O(N 2 ) Operations in Computing Boltzmann Collision Weight Matrix Irene M. Gamba, and Chenglong Zhang ICES, The University of Texas at Austin, 201 E. 24th
More informationEXPONENTIAL AND ALGEBRAIC RELAXATION IN KINETIC MODELS FOR WEALTH DISTRIBUTION
1 EXPONENTIAL AND ALGEBAIC ELAXATION IN KINETIC MODELS FO WEALTH DISTIBUTION B. DÜING Institut für Analysis und Scientific Computing, Technische Universität Wien, 1040 Wien, Austria. D. MATTHES and G.
More informationHypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th
Hypocoercivity and Sensitivity Analysis in Kinetic Equations and Uncertainty Quantification October 2 nd 5 th Department of Mathematics, University of Wisconsin Madison Venue: van Vleck Hall 911 Monday,
More informationOn the Optimal Choice of Coefficients in a Truncated Wild Sum and Approximate Solutions for the Kac Equation
Journal of Statistical Physics, Vol. 109, Nos. 1/2, October 2002 ( 2002) On the Optimal Choice of Coefficients in a Truncated Wild Sum and Approximate Solutions for the Kac Equation Eric A. Carlen 1 and
More informationEntropic structure of the Landau equation. Coulomb interaction
with Coulomb interaction Laurent Desvillettes IMJ-PRG, Université Paris Diderot May 15, 2017 Use of the entropy principle for specific equations Spatially Homogeneous Kinetic equations: 1 Fokker-Planck:
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 informationA Selective Modern History of the Boltzmann and Related Eq
A Selective Modern History of the Boltzmann and Related Equations October 2014, Fields Institute Synopsis 1. I have an ambivalent relation to surveys! 2. Key Words, Tools, People 3. Powerful Tools, I:
More informationlaplace s method for ordinary differential equations
Physics 24 Spring 217 laplace s method for ordinary differential equations lecture notes, spring semester 217 http://www.phys.uconn.edu/ rozman/ourses/p24_17s/ Last modified: May 19, 217 It is possible
More informationOn high energy tails in inelastic gases arxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Oct 2005
On high energy tails in inelastic gases arxiv:cond-mat/0510108v1 [cond-mat.stat-mech] 5 Oct 2005 R. Lambiotte a,b, L. Brenig a J.M. Salazar c a Physique Statistique, Plasmas et Optique Non-linéaire, Université
More informationOn the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions
On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions Vladislav A. Panferov Department of Mathematics, Chalmers University of Technology and Göteborg
More informationON A NEW CLASS OF WEAK SOLUTIONS TO THE SPATIALLY HOMOGENEOUS BOLTZMANN AND LANDAU EQUATIONS
ON A NEW CLASS OF WEAK SOLUTIONS TO THE SPATIALLY HOMOGENEOUS BOLTZMANN AND LANDAU EQUATIONS C. VILLANI Abstract. This paper deals with the spatially homogeneous Boltzmann equation when grazing collisions
More informationFast methods for the Boltzmann collision integral
Fast methods for the Boltzmann collision integral Clément Mouhot Lorenzo Pareschi Abstract In this note we present methods for the development of fast numerical schemes for the Boltzmann collision integral.
More informationCalculus II Practice Test 1 Problems: , 6.5, Page 1 of 10
Calculus II Practice Test Problems: 6.-6.3, 6.5, 7.-7.3 Page of This is in no way an inclusive set of problems there can be other types of problems on the actual test. To prepare for the test: review homework,
More informationarxiv: v2 [math.na] 12 Mar 2019
Approximation of the Boltzmann Collision Operator Based on Hermite Spectral Method Yanli Wang, Zhenning Cai March 3, 9 arxiv:83.9v [math.na] Mar 9 Abstract Based on the Hermite expansion of the distribution
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 informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationarxiv: v1 [math.na] 25 Oct 2018
Multi-scale control variate methods for uncertainty quantification in kinetic equations arxiv:80.0844v [math.na] 25 Oct 208 Giacomo Dimarco and Lorenzo Pareschi October 26, 208 Abstract Kinetic equations
More informationANSWER KEY. Page 1 Page 2 cake key pie boat glue cat sled pig fox sun dog fish zebra. Page 3. Page 7. Page 6
ANWE EY P 1 P 2 d fx d fh z P 3 P 4 m, m, m, m, m, m P 6 d d P 7 m P 5 m m P 10 d P 8 P 9 h,,, d d,, h v P 11 P 12,, m, m, m,, dd P 13 m f m P 18 h m P 22 f fx f fh fm fd P 14 m h P 19 m v x h d m P 15
More information08257 ANSWER KEY. Page 1 Page 2 cake key pie boat glue cat sled pig fox sun dog fish zebra. Page 3. Page 7. Page 6
P 1 P 2 d fx d fh z P 3 P 4 m, m, m, m, m, m P 6 d d P 7 m P 5 m m P 10 P 8 P 9 h,,, d d,, d v P 11 P 12,, m, m, m,, dd P 13 m f m P 18 h m P 22 f fx f fh fm fd P 14 m h P 19 m v x h d m P 15 m m P 20
More informationA QUICK INTRODUCTION TO KINETIC THEORY
A QUICK INTRODUCTION TO KINETIC THEORY GOVIND MENON Abstract. The purpose of these notes is to advertise kinetic theory to beginning graduate students with an interest in analysis and differential equations.
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 informationAlignment processes on the sphere
Alignment processes on the sphere Amic Frouvelle CEREMADE Université Paris Dauphine Joint works with : Pierre Degond (Imperial College London) and Gaël Raoul (École Polytechnique) Jian-Guo Liu (Duke University)
More informationCONVERGENCE OF THE FOURIER SERIES
CONVERGENCE OF THE FOURIER SERIES SHAW HAGIWARA Abstract. The Fourier series is a expression of a periodic, integrable function as a sum of a basis of trigonometric polynomials. In the following, we first
More informationLecture 22. m n c (k) i,j x i x j = c (k) k=1
Notes on Complexity Theory Last updated: June, 2014 Jonathan Katz Lecture 22 1 N P PCP(poly, 1) We show here a probabilistically checkable proof for N P in which the verifier reads only a constant number
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 informationSlowing down the neutrons
Slowing down the neutrons Clearly, an obvious way to make a reactor work, and to make use of this characteristic of the 3 U(n,f) cross-section, is to slow down the fast, fission neutrons. This can be accomplished,
More information18.01 Single Variable Calculus Fall 2006
MIT OpenCourseWare http://ocw.mit.edu 18.01 Single Variable Calculus Fall 2006 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Exam 4 Review 1. Trig substitution
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationFrom Hamiltonian particle systems to Kinetic equations
From Hamiltonian particle systems to Kinetic equations Università di Roma, La Sapienza WIAS, Berlin, February 2012 Plan of the lectures 1 Particle systems and BBKGY hierarchy (the paradigm of the Kinetic
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 informationhal , version 1-18 Sep 2013
Author manuscript, published in "Multiscale Modeling & Simulation 0 (202 792-87" ON DETERMINISTIC APPROXIMATION OF THE BOLTZMANN EQUATION IN A BOUNDED DOMAIN FRANCIS FILBET Abstract. In this paper we present
More informationA comparative study of no-time-counter and majorant collision frequency numerical schemes in DSMC
Purdue University Purdue e-pubs School of Aeronautics and Astronautics Faculty Publications School of Aeronautics and Astronautics 2012 A comparative study of no-time-counter and majorant collision frequency
More informationCalculus II Practice Test Problems for Chapter 7 Page 1 of 6
Calculus II Practice Test Problems for Chapter 7 Page of 6 This is a set of practice test problems for Chapter 7. This is in no way an inclusive set of problems there can be other types of problems on
More informationNumerical methods for plasma physics in collisional regimes
J. Plasma Physics (15), vol. 81, 358116 c Cambridge University Press 1 doi:1.117/s3778176 1 Numerical methods for plasma physics in collisional regimes G. Dimarco 1,Q.Li,L.Pareschi 1 and B. Yan 3 1 Department
More informationChapter 2 A Brief Introduction to the Mathematical Kinetic Theory of Classical Particles
Chapter 2 A Brief Introduction to the Mathematical Kinetic Theory of Classical Particles 2.1 Plan of the Chapter The contents of this chapter are motivated by the second key question posed in Section 1.3,
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 informationAn Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation
An Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation Wei Ren a, Hong Liu a,, Shi Jin b,c a J C Wu Center for Aerodynamics, School of Aeronautics and Aerospace, Shanghai Jiao Tong University,
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 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 informationMath 21B - Homework Set 8
Math B - Homework Set 8 Section 8.:. t cos t dt Let u t, du t dt and v sin t, dv cos t dt Let u t, du dt and v cos t, dv sin t dt t cos t dt u v v du t sin t t sin t dt [ t sin t u v ] v du [ ] t sin t
More informationAccurate Deterministic Solutions for the Classic Boltzmann Shock Profile
City University of New York (CUNY) CUNY Academic Works Dissertations, Theses, and Capstone Projects Graduate Center 9-2015 Accurate Deterministic Solutions for the Classic Boltzmann Shock Profile Yubei
More informationBOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES
BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES Vicente Garzó Departamento de Física, Universidad de Extremadura Badajoz, SPAIN Collaborations Antonio Astillero, Universidad de Extremadura José
More informationMath 102 Spring 2008: Solutions: HW #3 Instructor: Fei Xu
Math Spring 8: Solutions: HW #3 Instructor: Fei Xu. section 7., #8 Evaluate + 3 d. + We ll solve using partial fractions. If we assume 3 A + B + C, clearing denominators gives us A A + B B + C +. Then
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. The 8th International Conference on Computational
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationMath 1B Final Exam, Solution. Prof. Mina Aganagic Lecture 2, Spring (6 points) Use substitution and integration by parts to find:
Math B Final Eam, Solution Prof. Mina Aganagic Lecture 2, Spring 20 The eam is closed book, apart from a sheet of notes 8. Calculators are not allowed. It is your responsibility to write your answers clearly..
More informationQuadratic reciprocity (after Weil) 1. Standard set-up and Poisson summation
(September 17, 010) Quadratic reciprocity (after Weil) Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ I show that over global fields (characteristic not ) the quadratic norm residue
More informationRecent advances in kinetic theory for mixtures of polyatomic gases
Recent advances in kinetic theory for mixtures of polyatomic gases Marzia Bisi Parma University, Italy Conference Problems on Kinetic Theory and PDE s Novi Sad (Serbia), September 25 27, 2014 M. Bisi,
More informationMATH 220 solution to homework 4
MATH 22 solution to homework 4 Problem. Define v(t, x) = u(t, x + bt), then v t (t, x) a(x + u bt) 2 (t, x) =, t >, x R, x2 v(, x) = f(x). It suffices to show that v(t, x) F = max y R f(y). We consider
More informationTwo Simple Nonlinear Edge Detectors
Two Simple Nonlinear Edge Detectors Mladen Victor Wickerhauser Washington University in St. Louis, Missouri victor@math.wustl.edu http://www.math.wustl.edu/~victor International Conference on Wavelet Analysis
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationA Time Splitting for the Semiclassical Schrödinger Equation
A Time Splitting for the Semiclassical Schrödinger Equation V. Gradinaru Seminar for Applied Mathematics, ETH Zürich, CH 809, Zürich, Switzerland vasile.gradinaru@sam.math.ethz.ch and G. A. Hagedorn Department
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 informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationFluid 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 informationMath 162: Calculus IIA
Math 62: Calculus IIA Final Exam ANSWERS December 9, 26 Part A. (5 points) Evaluate the integral x 4 x 2 dx Substitute x 2 cos θ: x 8 cos dx θ ( 2 sin θ) dθ 4 x 2 2 sin θ 8 cos θ dθ 8 cos 2 θ cos θ dθ
More informationSpectral - Lagrangian methods for Collisional Models of Non - Equilibrium Statistical States
Spectral - Lagrangian methods for Collisional Models of Non - Equilibrium Statistical States Irene M. Gamba Dept. of Mathematics & Institute of Computational Engineering and Sciences, University of Texas
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 informationTopics in Fourier analysis - Lecture 2.
Topics in Fourier analysis - Lecture 2. Akos Magyar 1 Infinite Fourier series. In this section we develop the basic theory of Fourier series of periodic functions of one variable, but only to the extent
More informationAccepted Manuscript. Deterministic numerical solutions of the Boltzmann equation using the fast spectral method
Accepted Manuscript Deterministic numerical solutions of the Boltzmann equation using the fast spectral method Lei Wu, Craig White, Thomas J. Scanlon, Jason M. Reese, Yonghao Zhang PII: S1-9991(13)37-6
More informationHOMEWORK 3 MA1132: ADVANCED CALCULUS, HILARY 2017
HOMEWORK MA112: ADVANCED CALCULUS, HILARY 2017 (1) A particle moves along a curve in R with position function given by r(t) = (e t, t 2 + 1, t). Find the velocity v(t), the acceleration a(t), the speed
More informationLast lecture: Recurrence relations and differential equations. The solution to the differential equation dx
Last lecture: Recurrence relations and differential equations The solution to the differential equation dx = ax is x(t) = ce ax, where c = x() is determined by the initial conditions x(t) Let X(t) = and
More informationDevelopment of a 3D Spectral Boltzmann Solver for Plasma Modeling
Development of a 3D Spectral Boltzmann Solver for Plasma Modeling IEPC-9-38 Presented at the 3 st International Electric Propulsion Conference, University of Michigan, Ann Arbor, Michigan, USA Jacques
More informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
More informationPhysics Dec The Maxwell Velocity Distribution
Physics 301 7-Dec-2005 29-1 The Maxwell Velocity Distribution The beginning of chapter 14 covers some things we ve already discussed. Way back in lecture 6, we calculated the pressure for an ideal gas
More informationAn asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling
An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3 Saint-Malo 13 December 2016 1 Université
More informationNoether s Theorem. 4.1 Ignorable Coordinates
4 Noether s Theorem 4.1 Ignorable Coordinates A central recurring theme in mathematical physics is the connection between symmetries and conservation laws, in particular the connection between the symmetries
More informationHölder regularity estimation by Hart Smith and Curvelet transforms
Hölder regularity estimation by Hart Smith and Curvelet transforms Jouni Sampo Lappeenranta University Of Technology Department of Mathematics and Physics Finland 18th September 2007 This research is done
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. Workshop Asymptotic-Preserving schemes, Porquerolles.
More informationA discontinuous Galerkin fast spectral method for the multi-species Boltzmann equation
A discontinuous Galerkin fast spectral method for the multi-species Boltzmann equation Shashank Jaiswal a, Alina A. Alexeenko a, Jingwei Hu b, a School of Aeronautics and Astronautics, Purdue University,
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationDIFFERENTIAL GEOMETRY HW 4. Show that a catenoid and helicoid are locally isometric.
DIFFERENTIAL GEOMETRY HW 4 CLAY SHONKWILER Show that a catenoid and helicoid are locally isometric. 3 Proof. Let X(u, v) = (a cosh v cos u, a cosh v sin u, av) be the parametrization of the catenoid and
More informationRandom Averaging. Eli Ben-Naim Los Alamos National Laboratory. Paul Krapivsky (Boston University) John Machta (University of Massachusetts)
Random Averaging Eli Ben-Naim Los Alamos National Laboratory Paul Krapivsky (Boston University) John Machta (University of Massachusetts) Talk, papers available from: http://cnls.lanl.gov/~ebn Plan I.
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 7 Beyond Atoms: Simulating Molecules David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu Review Fundamentals units, properties, statistical
More information11.D.2. Collision Operators
11.D.. Collision Operators (11.94) can be written as + p t m h+ r +p h+ p = C + h + (11.96) where C + is the Boltzmann collision operator defined by [see (11.86a) for convention of notations] C + g(p)
More informationVector Fields. It is standard to define F µν = µ ϕ ν ν ϕ µ, so that the action may be written compactly as
Vector Fields The most general Poincaré-invariant local quadratic action for a vector field with no more than first derivatives on the fields (ensuring that classical evolution is determined based on the
More informationSpring 2017 Econ 574 Roger Koenker. Lecture 14 GEE-GMM
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 14 GEE-GMM Throughout the course we have emphasized methods of estimation and inference based on the principle
More informationMonte Carlo simulation inspired by computational optimization. Colin Fox Al Parker, John Bardsley MCQMC Feb 2012, Sydney
Monte Carlo simulation inspired by computational optimization Colin Fox fox@physics.otago.ac.nz Al Parker, John Bardsley MCQMC Feb 2012, Sydney Sampling from π(x) Consider : x is high-dimensional (10 4
More informationFix(g). Orb(x) i=1. O i G. i=1. O i. i=1 x O i. = n G
Math 761 Fall 2015 Homework 4 Drew Armstrong Problem 1 Burnside s Lemma Let X be a G-set and for all g G define the set Fix(g : {x X : g(x x} X (a If G and X are finite, prove that Fix(g Stab(x g G x X
More informationDYNAMICS OF SERIAL ROBOTIC MANIPULATORS
DYNAMICS OF SERIAL ROBOTIC MANIPULATORS NOMENCLATURE AND BASIC DEFINITION We consider here a mechanical system composed of r rigid bodies and denote: M i 6x6 inertia dyads of the ith body. Wi 6 x 6 angular-velocity
More informationAnalisis para la ecuacion de Boltzmann Soluciones y Approximaciones
Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Irene M. Gamba Department of Mathematics and ICES The University of Texas at Austin Buenos Aires, June 2012 Collaborators: R. Alonso,
More informationKey to Complex Analysis Homework 1 Spring 2012 (Thanks to Da Zheng for providing the tex-file)
Key to Complex Analysis Homework 1 Spring 212 (Thanks to Da Zheng for providing the tex-file) February 9, 212 16. Prove: If u is a complex-valued harmonic function, then the real and the imaginary parts
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5
More information