Integration of Vlasov-type equations
|
|
- Mervyn Elliott
- 6 years ago
- Views:
Transcription
1 Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017
2 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma Source: Hubble telescope Source: Luc Viatour /
3 Vlasov equation Evolution of (collisionless) plasma is described by t f (t, x, v) + v x f (t, x, v) + F v f (t, x, v) = 0, where f... particle density function v... velocity x... position F... force field Coupled to the self-consistent electric (and magnetic) field via F = e (E + v B) m yields the Vlasov Poisson (Vlasov Maxwell) equations. Applications: Tokamaks (fusion), plasma-laser interactions,...
4 Vlasov equation, cont. Anatoly Alexandrovich Vlasov ( ), Moscow State U The vibrational properties of an electron gas (1938, 1968) e t f (t, x, v )+v x f (t, x, v )+ E +v B v f (t, x, v ) = 0 m
5 Filamentation of phase space Simulation of the filamentation of phase space as is studied, for example, in the context of Landau damping. The horizontal axis represents space and the vertical one velocity.
6 Bump-on-tail instability The bump-on-tail instability leads to a traveling vortex in phase space. The horizontal axis represents space and the vertical one velocity.
7 Outline Strang splitting for Vlasov-type equations dg space discretization Numerical examples References
8 Abstract formulation Abstract initial value problem where t f (t) = (A + B)f (t), f (0) = f 0, A is a linear differential operator; the nonlinear operator B has the form Bf = B(f )f, where B(f ) is an (unbounded) linear operator. Formulation comprises the Vlasov Poisson and the Vlasov Maxwell equations as special cases: take A = v x and choose B appropriately
9 Example: Vlasov Poisson equation in 1+1 d Recall the abstract formulation t f (t) = (A + B)f (t), f (0) = f 0. Example: Vlasov Poisson equations in 1+1 dimensions t f (t, x, v) = v x f (t, x, v) E(f (t,, ), x) v f (t, x, v) x E(f (t,, ), x) = f (t, x, v) dv 1, E = x Φ R f (0, x, v) = f 0 (x, v) with periodic boundary conditions in space f (t, x, v) = f (t, x + L, v).
10 Strang splitting Problem: t f (t) = (A + B)f (t), Bf = B(f )f. Denote f k f (t k ) at time t k = kτ with step size τ. Solve t f = Af, t g = B k+1/2 g, where B k+1/2 B ( f (t k + τ 2 )) and let f k+1 = S k f k, where the (nonlinear) splitting operator S k is given by S k = e τ 2 A e τb k+1/2 e τ 2 A. For grid-based Vlasov solvers: Cheng and Knorr 1976; Mangeney et al. 2002
11 Computational efficiency and stability Splitting is very efficient in our application: t f = v x f f (τ, x, v) = e τv x f (0, x, v) = f (0, x τv, v) and t g = E k+1/2 (x) v g g(τ, x, v) = e τe k+1/2(x) v g(0, x, v) = g(0, x, v τe k+1/2 (x)) imply S k = 1. Moreover, we have ( ) e τb k+1/2 e τb(f (t k+τ/2)) e τa f (t k ) Cτ f k+1/2 f (t k +τ/2) (requires some smoothness of f ).
12 Consistency technical tools Lemma [Gröbner Alekseev formula] Consider Then f (t) = G(f (t)) + R(f (t)), f (0) = f 0 g (t) = G(g(t)) with solution E G (t, g 0 ) f (t) = E G (t, f 0 ) + t 0 2 E G (t s, f (s)) R (f (s)) ds. Proof. Fundamental theorem of calculus + trick. Lemma. Let E generate a (semi)group. Then e τe = m 1 k=0 τ k k! E k + τ m E m 1 e (1 s)τe 0 s m 1 (m 1)! ds Proof. Integration by parts.
13 Proof of the Gröbner Alekseev formula Fix t and let u(t) be a solution of u (t) = G (u(t)). Next differentiate the identity with respect to s to get E G (t s, u(s)) = u(t) 1 E G (t s, u(s)) + 2 E G (t s, u(s)) G (u(s)) = 0. The initial value of u is now chosen such that u(s) = f (s) which implies 1 E G (t s, f (s)) + 2 E G (t s, f (s)) G (f (s)) = 0.
14 Proof of the Gröbner Alekseev formula, cont. Let ϕ(s) = E G (t s, f (s)). By the fundamental theorem of calculus f (t) E G (t, f 0 ) = = = t 0 t 0 t 0 ϕ (s) ds ( ) 1 E G (t s, f (s)) + 2 E G (t s, f (s))f (s) ds 2 E G (t s, f (s)) R (f (s)) ds, where we have used f (s) = G(f (s)) + R(f (s)) and 1 E G (t s, f (s)) + 2 E G (t s, f (s)) G (f (s)) = 0.
15 Consistency, cont. Expansion of the exact solution f (τ) = E B (τ, f 0 ) τ τ1 0 0 τ τ1 τ2 0 0 τ 0 2 E B (τ τ 1, f (τ 1 ))AE B (τ 1, f 0 ) dτ 1 2 E B (τ τ 1, f (τ 1 ))A 2 E B (τ 1 τ 2, f (τ 2 ))AE B (τ 2, f 0 )dτ 2 dτ 1 ( 2 ) 2 E B (τ k τ k+1, f (τ k+1 ))A f (τ 3 ) dτ 3 dτ 2 dτ 1, 0 where τ 0 = τ. k=0 Expansion of the numerical solution S k f 0 = e τ 2 A e τb k+1/2 e τ 2 A f 0 S 0 f 0 = e τb 1/2 f 0 + τ { } A, e τb 1/2 f0 + τ 2 { { }} A, A, e τb 1/2 f0 +R 3 f 0 2 8
16 Consistency, cont. Compare the terms by employing quadrature rules (Jahnke, Lubich 2000). Treat τ term with trapezoidal rule τ 0 2 E B (τ τ 1, f (τ 1 ))AE B (τ 1, f 0 ) dτ 1 = and compare with τ 2 2E B (τ, f 0 )Af 0 + τ 2 AE B(τ, f 0 )f 0 + R 1 τ { } A, e τb 1/2 f0 = τ 2 2 A e τb 1/2 f 0 + τ 2 e τb 1/2 Af 0
17 Consistency, cont. Estimate the difference ( e τb 1/2 2 E B (τ, f 0 ) ) Af 0 = 2 ( e τb 1/2 E B (τ, f 0 ) ) Af 0. Apply again the variation of constants formula E B (τ, f 0 ) e τb 1/2 f 0 = τ 0 e (τ s)b 1/2 (B B 1/2 )E B (s, f 0 )ds. Gives condition on B 1/2. Treat all terms in this way and estimate remainders.
18 Regularity of exact solution Consider Vlasov Poisson problem as before. Theorem (Glassey 1996) Assume that f 0 C m per,c is non-negative, then f C m (0, T ; C m per,c), E(f (t,, ), x) C m (0, T ; C m per). Moreover, there exists Q(T ) such that for all t [0, T ] and x R it holds that supp f (t, x, ) [ Q(T ), Q(T )]. Important result: Regularity of initial data is preserved.
19 Abstract convergence result Theorem (L. Einkemmer, A.O., SIAM J. Numer. Anal., 2014a) Let f 0 C 3 per,c be non-negative f k+1/2 f (t k + τ ) Cτ 2 Then Strang splitting for the Vlasov Poisson equations is second-order convergent. Possible choices for f k+1/2 : first-order Lie Trotter splitting f k+1/2 = e τ 2 B(fk) e τ 2 A f k due to the special structure of the E field f k+1/2 = e τ 2 A f k linear extrapolation using f k 1 and f k
20 Outline Strang splitting for Vlasov-type equations dg space discretization Numerical examples References
21 dg approximation in space and velocity Choose grid in (x, v) plane with cells R ij. Project onto space of Legendre polynomials of order l. The projection is given by Pg Rij = l l b ij km P(1) k k=0 m=0 (x)p(2) m (v) with coefficients b ij (2k + 1)(2m + 1) km = h x h v P (1) k (x)p(2) m (v)g(x, v) d(x, v), R ij where P (1) k and P m (2) are translated and scaled Legendre polynomials.
22 Space discretization Pe τ 2 A g(x, v) = Pg ( x τ v, v) 2 cell 1 cell 2 cell 1 cell 2 Polynomial approximation of functions with a small jump discontinuity; Jump heights ε (k) = g (k) (x 0 +) g (k) (x 0 ) satisfy ε (k) ch l k+1 for all k {0,..., l}; Then g (k) (Pg) (k) Ch l k+1, 0 k l.
23 Space and velocity discretization Pe τ 2 A g(x, v) = Pg ( x τ v, v) 2 Pe τb(pf k+1/2) g(x, v) = Pg ( x, v τe k+1/2 (x) ) evaluate integrals exactly with Gauss Legendre quadrature Spatially discretized Strang splitting: S k = Pe τ 2 A Pe τb(pf k+1/2) Pe τ 2 A
24 Convergence of full discretization Theorem (L. Einkemmer, A.O., SIAM J. Numer. Anal., 2014b) Suppose that f 0 C max{l+1,3}, nonnegative and compactly supported. Then ( n 1 ) ) S k Pf 0 f (nτ) (τ C 2 + hl+1 + h l+1. τ k=0 Proof. Combine time and space error
25 Stability, revisited In advection dominated problems instabilities can occur if the numerical integration is not performed exactly. Example (Molenkamp Crowley test problem) t f (t, x, y) = 2π(y x x y )f (t, x, y) f (0, x, v) = f 0 (x, y), with off-centered Gaussian as initial value. Initial value turns around origin with period 1. Finite element schemes turn out to be unstable for most quadrature rules [Morton, Priestley, Süli (1988)].
26 Molenkamp Crowley test problem, cont. Morton, Priestley, Süli: Stability of the Lagrange Galerkin method with non-exact integration, RAIRO (1988)
27 Molenkamp Crowley test problem, cont. 1 1 max: 1.0 min: max: min: y y x x Projected initial value (left); numerical solution after 60 revolutions with τ = 0.02, N x = N v = 40, and l = 2 (right). Our dg method is stable.
28 Outline Strang splitting for Vlasov-type equations dg space discretization Numerical examples References
29 Example: nonlinear Landau damping Vlasov Poisson equations in 1+1 dimensions on [0, 4π] R together with the initial value f 0 (x, v) = 1 ( e v 2 /2 1 + α cos x ). 2π 2 Landau damping is a popular test problem. weak or linear Landau damping α = 0.01 strong or nonlinear Landau damping α = 0.5 Consider Strang and Lie Trotter splitting S L,k = e τb k e τa.
30 Nonlinear Landau damping; time error error (discrete L1 norm) e-005 Strang splitting Lie-Trotter splitting order 1 order 2 1e step size Error (in discrete L 1 norm) of the particle density function f (1,, ) for Strang and Lie Trotter splitting, respectively. h
31 Nonlinear Landau damping; space error error (infinity norm) e-05 1e-06 l=0 l=1 l=2 order 1 order 2 order 3 1e Spatial error (in maximum norm) of the particle density function for dg discretizations of Strang splitting. N
32 Recurrence phenomenon Most easily understood for a transport equation. Consider t f = v x f, f 0 (x, v) = e v 2 /2 2π ( cos(0.5x) ) with exact solution f (t, x, v) = e v 2 /2 2π ( cos(0.5x 0.5vt) ). The electric energy E(t) = 4π decays exponentially in time. 0 E(t, x) 2 dx = π 1250 e 0.25t2
33 Recurrence, cont. For the numerical approximation of f (t, x, v) = e v 2 /2 2π ( cos(0.5x 0.5vt) ). consider a piecewise constant approximation in v with grid size h v. Let the time t p be given by h v t p = 4π, v j t p = jh v t p = 4jπ. Obviously, the numerical solution is periodic with period t p. The numerical approximation to the electric energy E is thus periodic in time as well.
34 Recurrence for the advection equation Nv=32, l=0 Nv=64, l=0 Nv=32, l=2 exact solution electric energy e-06 1e-08 1e-10 1e-12 1e-14 1e t There is no periodicity in the approximation of degree 2. Nevertheless, a recurrence-like effect is still visible.
35 Recurrence, weak Landau damping The electric energy decays asymptotically as e 2γt, γ We compare our numerical solutions with this function Nv=64, l=1 Nv=64, l=2 Nv=64, l=4 analytic decay rate Nv=128, l=1 Nv=128, l=2 analytic decay rate electric energy 1e-06 1e-08 1e-10 electric energy 1e-06 1e-08 1e-10 1e-12 1e-12 1e-14 1e-14 1e e t t The decay of the electric energy is shown for N x = N v = 64 (left) and N x = N v = 128 (right). In all cases a relatively large time step of τ = 0.2 is employed.
36 A strategy to suppress the recurrence effect Recurrence is not a consequence of the lack of accuracy of the underlying space discretization, but a result of aliasing. higher and higher frequencies in phase space are created; aliasing of the high frequencies introduces an error in the macroscopic quantities (such as the electric energy). Strategy to avoid recurrence (Einkemmer, A.O., EJPD 2014): damp highest frequencies; in our dg implementation: FFTs of coefficients of Legendre polynomials; after a cutoff time t c, cut off the highest M frequencies cutoff (t c,m)
37 Recurrence, weak Landau damping piecewise constant approximation, N=64 electric energy e-06 1e-08 1e-10 1e-12 1e-14 1e time no cutoff, N=128 cutoff (1.2,5) piecewise linear aproximation, N=64 electric energy e-06 1e-08 1e-10 1e-12 1e-14 1e time no cutoff cutoff (1.2,5) piecewise quadratic Integration approximation, of Vlasov-type N=64 equations
38 Example: plasma echo phenomenon Initial value f 0 (x, v) = 1 2π e v 2 /2 (1 + α cos(k 1 x)) corresponds to an excitation with wavenumber k 1 at time t = 0. Furthermore, at time t = t 2 we excite a second perturbation with wavenumber k 2 ; that is, we superimpose α 2π e v 2 /2 cos(k 2 x) on the (numerical) solution at time t 2. Data: α = 10 3, k 1 = 12π/100, k 2 = 25π/100, and t 2 = 200
39 Echo, low resolution no cutoff, 5000 points in v-space (piecewise-constant) electric energy e-06 1e-08 1e-10 1e-12 1e-14 echo time no cutoff, 500 points in v-space (piecewise-constant) electric energy e-06 1e-08 1e-10 1e-12 1e-14 echo time electric energy e-06 1e-08 1e-10 1e-12 1e-14 cutoff (1,20), 500 points in v-space (piecewise-constant) echo time
40 electric e 1e-08 1e-10 1e-12 Echo, higher 1e-14resolution time electric energy electric energy cutoff no cutoff, (1,20), 5000 points in in v-space (piecewise-constant) e-06 echo 1e-08 1e-10 1e-12 1e time no cutoff, points in v-space (piecewise-constant) e-06 echo 1e-08 1e-10 1e-12 1e time electric energy d) cutoff (1,20), points in v-space (piecewise-constant) e-06 echo 1e-08 1e-10 1e-12 1e time
41 References L. Einkemmer, A. Ostermann. Convergence analysis of Strang splitting for Vlasov-type equations. SIAM J. Numer. Anal. 52, (2014) L. Einkemmer, A. Ostermann. Convergence analysis of a discontinuous Galerkin / Strang splitting approximation for the Vlasov Poisson equations. SIAM J. Numer. Anal. 52, (2014) L. Einkemmer, A. Ostermann. A strategy to suppress recurrence in grid-based Vlasov solvers. Eur. Phys. J. D 68, 197 (2014)
Splitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationarxiv: v2 [math.na] 5 Jul 2016
AN ASYMPTOTIC PRESERVING SCHEME FOR THE RELATIVISTIC VLASOV MAXWELL EQUATIONS IN THE CLASSICAL LIMIT by arxiv:1602.09062v2 [math.na] 5 Jul 2016 Nicolas Crouseilles, Lukas Einkemmer & Erwan Faou Abstract.
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
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 informationequations in the semiclassical regime
Adaptive time splitting for nonlinear Schrödinger equations in the semiclassical regime Local error representation and a posteriori error estimator Workshop, Numerical Analysis of Evolution Equations Vill,
More informationKinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence
2013 US-Japan JIFT workshop on New Aspects of Plasmas Kinetic Simulation NIFS, November 22-23, 2013 Kinetic damping in gyro-kinetic simulation and the role in multi-scale turbulence cf. Revisit for Landau
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 informationA particle-in-cell method with adaptive phase-space remapping for kinetic plasmas
A particle-in-cell method with adaptive phase-space remapping for kinetic plasmas Bei Wang 1 Greg Miller 2 Phil Colella 3 1 Princeton Institute of Computational Science and Engineering Princeton University
More informationDIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS. Eskil Hansen and Alexander Ostermann
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume X, Number 0X, XX 00X pp. 1 11 DIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS Eskil Hansen and Alexander Ostermann Institut
More informationNhung Pham 1, Philippe Helluy 2 and Laurent Navoret 3
ESAIM: PROCEEDINGS AND SURVEYS, September 014, Vol. 45, p. 379-389 J.-S. Dhersin, Editor HYPERBOLIC APPROXIMATION OF THE FOURIER TRANSFORMED VLASOV EQUATION Nhung Pham 1, Philippe Helluy and Laurent Navoret
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 informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationApplying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models
0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,
More informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
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 informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationarxiv: v1 [math.na] 6 Nov 2017
Efficient boundary corrected Strang splitting Lukas Einkemmer Martina Moccaldi Alexander Ostermann arxiv:1711.02193v1 [math.na] 6 Nov 2017 Version of November 6, 2017 Abstract Strang splitting is a well
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
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 informationMeshfree Exponential Integrators
Meshfree joint work with A. Ostermann (Innsbruck) M. Caliari (Verona) Leopold Franzens Universität Innsbruck Innovative Integrators 3 October 2 Meshfree Problem class: Goal: Time-dependent PDE s with dominating
More informationEnergy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media
Energy Stable Discontinuous Galerkin Methods for Maxwell s Equations in Nonlinear Optical Media Yingda Cheng Michigan State University Computational Aspects of Time Dependent Electromagnetic Wave Problems
More informationExponential integrators
Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1 Outline Time dependent Schrödinger equations
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationA Discontinuous Galerkin Method for Vlasov Systems
A Discontinuous Galerkin Method for Vlasov Systems P. J. Morrison Department of Physics and Institute for Fusion Studies The University of Texas at Austin morrison@physics.utexas.edu http://www.ph.utexas.edu/
More informationLARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION
Differential Integral Equations Volume 3, Numbers 1- (1), 61 77 LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION Robert Glassey Department of Mathematics, Indiana University
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationVolterra Integral Equations of the First Kind with Jump Discontinuous Kernels
Volterra Integral Equations of the First Kind with Jump Discontinuous Kernels Denis Sidorov Energy Systems Institute, Russian Academy of Sciences e-mail: contact.dns@gmail.com INV Follow-up Meeting Isaac
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationNumerical Methods I Orthogonal Polynomials
Numerical Methods I Orthogonal Polynomials Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course G63.2010.001 / G22.2420-001, Fall 2010 Nov. 4th and 11th, 2010 A. Donev (Courant Institute)
More informationWhat place for mathematicians in plasma physics
What place for mathematicians in plasma physics Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 15-19 September 2008 Eric Sonnendrücker (U. Strasbourg) Math
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationMath 2C03 - Differential Equations. Slides shown in class - Winter Laplace Transforms. March 4, 5, 9, 11, 12, 16,
Math 2C03 - Differential Equations Slides shown in class - Winter 2015 Laplace Transforms March 4, 5, 9, 11, 12, 16, 18... 2015 Laplace Transform used to solve linear ODEs and systems of linear ODEs with
More informationSmall BGK waves and nonlinear Landau damping (higher dimensions)
Small BGK waves and nonlinear Landau damping higher dimensions Zhiwu Lin and Chongchun Zeng School of Mathematics Georgia Institute of Technology Atlanta, GA, USA Abstract Consider Vlasov-Poisson system
More informationBenchmarks in Computational Plasma Physics
Benchmarks in Computational Plasma Physics P. Londrillo INAF, Bologna, Italie S. Landi Università di Firenze, Italie What you compute when you do computations of the Vlasov equation? Overview A short review
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More information1 Introduction and Notations he simplest model to describe a hot dilute electrons moving in a fixed ion background is the following one dimensional Vl
Numerical Study on Landau Damping ie Zhou 1 School of Mathematical Sciences Peking University, Beijing 100871, China Yan Guo 2 Lefschetz Center for Dynamical Systems Division of Applied Mathematics Brown
More informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationSpatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations
Spatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations N. Chalmers and L. Krivodonova March 5, 014 Abstract We apply the discontinuous Galerkin finite element method
More informationLANDAU DAMPING A. GRIFFIN, J.J. KOSKELA, F.Q. NAZAR, T. RANNER
LANDAU DAMPING A. GRIFFIN, J.J. KOSKELA, F.Q. NAZAR, T. RANNER 1. Introduction Landau damping is a phenomena observed in plasma wherein there is an exponential decay in the oscillations of the number density
More informationu = (A + F )u, u(0) = η, (1.1)
A CONVERGENCE ANALYSIS OF THE PEACEMAN RACHFORD SCHEME FOR SEMILINEAR EVOLUTION EQUATIONS ESKIL HANSEN AND ERIK HENNINGSSON Abstract. The Peaceman Rachford scheme is a commonly used splitting method for
More informationNumerical Methods I: Numerical Integration/Quadrature
1/20 Numerical Methods I: Numerical Integration/Quadrature Georg Stadler Courant Institute, NYU stadler@cims.nyu.edu November 30, 2017 umerical integration 2/20 We want to approximate the definite integral
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationMATHia Unit MATHia Workspace Overview TEKS
1 Function Overview Searching for Patterns Exploring and Analyzing Patterns Comparing Familiar Function Representations Students watch a video about a well-known mathematician creating an expression for
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More information1.5 Approximate Identities
38 1 The Fourier Transform on L 1 (R) which are dense subspaces of L p (R). On these domains, P : D P L p (R) and M : D M L p (R). Show, however, that P and M are unbounded even when restricted to these
More informationMassively parallel semi-lagrangian solution of the 6d Vlasov-Poisson problem
Massively parallel semi-lagrangian solution of the 6d Vlasov-Poisson problem Katharina Kormann 1 Klaus Reuter 2 Markus Rampp 2 Eric Sonnendrücker 1 1 Max Planck Institut für Plasmaphysik 2 Max Planck Computing
More informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
More informationNumerical Integration
Numerical Integration Sanzheng Qiao Department of Computing and Software McMaster University February, 2014 Outline 1 Introduction 2 Rectangle Rule 3 Trapezoid Rule 4 Error Estimates 5 Simpson s Rule 6
More informationElectromagnetic Relaxation Time Distribution Inverse Problems in the Time-domain
Electromagnetic Relaxation Time Distribution Inverse Problems in the Time-domain Prof Nathan L Gibson Department of Mathematics Joint Math Meeting Jan 9, 2011 Prof Gibson (OSU) Inverse Problems for Distributions
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationHybrid Discontinuous/Continuous Galerkin Algorithms for (Gyro) Kinetic Simulations of Plasmas
Hybrid Discontinuous/Continuous Galerkin Algorithms for (Gyro) Kinetic Simulations of Plasmas A. H. Hakim G. W. Hammett Eric Shi Princeton Plasma Physics Laboratory, Princeton, NJ ammar@princeton.edu http://www.ammar-hakim.org/sj
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationAdaptive semi-lagrangian schemes for transport
for transport (how to predict accurate grids?) Martin Campos Pinto CNRS & University of Strasbourg, France joint work Albert Cohen (Paris 6), Michel Mehrenberger and Eric Sonnendrücker (Strasbourg) MAMCDP
More informationNEW FUNCTIONAL INEQUALITIES
1 / 29 NEW FUNCTIONAL INEQUALITIES VIA STEIN S METHOD Giovanni Peccati (Luxembourg University) IMA, Minneapolis: April 28, 2015 2 / 29 INTRODUCTION Based on two joint works: (1) Nourdin, Peccati and Swan
More informationHarbor Creek School District
Unit 1 Days 1-9 Evaluate one-sided two-sided limits, given the graph of a function. Limits, Evaluate limits using tables calculators. Continuity Evaluate limits using direct substitution. Differentiability
More informationEnergy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas
Energy-Conserving Numerical Simulations of Electron Holes in Two-Species Plasmas Yingda Cheng Andrew J. Christlieb Xinghui Zhong March 18, 2014 Abstract In this paper, we apply our recently developed energy-conserving
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More information2.6. Graphs of Rational Functions. Copyright 2011 Pearson, Inc.
2.6 Graphs of Rational Functions Copyright 2011 Pearson, Inc. Rational Functions What you ll learn about Transformations of the Reciprocal Function Limits and Asymptotes Analyzing Graphs of Rational Functions
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationHeat kernels of some Schrödinger operators
Heat kernels of some Schrödinger operators Alexander Grigor yan Tsinghua University 28 September 2016 Consider an elliptic Schrödinger operator H = Δ + Φ, where Δ = n 2 i=1 is the Laplace operator in R
More informationCompensation. June 29, Arizona State University. Edge Detection from Nonuniform Data via. Convolutional Gridding with Density.
Arizona State University June 29, 2012 Overview gridding with density compensation is a common method for function reconstruction nonuniform Fourier data. It is computationally inexpensive compared with
More information5 Numerical Integration & Dierentiation
5 Numerical Integration & Dierentiation Department of Mathematics & Statistics ASU Outline of Chapter 5 1 The Trapezoidal and Simpson Rules 2 Error Formulas 3 Gaussian Numerical Integration 4 Numerical
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationDipl.-Ing. Lukas Einkemmer, PhD MSc BSc
Dipl.-Ing., PhD MSc BSc Thurnfelsstraße 39 Völs 6176, Austria lukas.einkemmer@uibk.ac.at +43664/2254179 Education 2011 2014 PhD in applied mathematics University of Innsbruck, University of California
More informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
More informationApplied Math 205. Full office hour schedule:
Applied Math 205 Full office hour schedule: Rui: Monday 3pm 4:30pm in the IACS lounge Martin: Monday 4:30pm 6pm in the IACS lounge Chris: Tuesday 1pm 3pm in Pierce Hall, Room 305 Nao: Tuesday 3pm 4:30pm
More informationAdaptive simulation of Vlasov equations in arbitrary dimension using interpolatory hierarchical bases
Adaptive simulation of Vlasov equations in arbitrary dimension using interpolatory hierarchical bases Erwan Deriaz Fusion Plasma Team Institut Jean Lamour (Nancy), CNRS / Université de Lorraine VLASIX
More informationON THE VLASOV-MAXWELL SYSTEM WITH A STRONG EXTERNAL MAGNETIC FIELD
ON THE VLASOV-MAXWELL SYSTEM WITH A STONG EXTENAL MAGNETIC FIELD Francis Filbet, Tao Xiong, Eric Sonnendr ucker To cite this version: Francis Filbet, Tao Xiong, Eric Sonnendr ucker. ON THE VLASOV-MAXWELL
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationPolynomial Preserving Jump-Diffusions on the Unit Interval
Polynomial Preserving Jump-Diffusions on the Unit Interval Martin Larsson Department of Mathematics, ETH Zürich joint work with Christa Cuchiero and Sara Svaluto-Ferro 7th General AMaMeF and Swissquote
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
More informationHUDSONVILLE HIGH SCHOOL COURSE FRAMEWORK
HUDSONVILLE HIGH SCHOOL COURSE FRAMEWORK COURSE / SUBJECT A P C a l c u l u s ( B C ) KEY COURSE OBJECTIVES/ENDURING UNDERSTANDINGS OVERARCHING/ESSENTIAL SKILLS OR QUESTIONS Limits and Continuity Derivatives
More informationLocally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials
Locally Linearized Euler Equations in Discontinuous Galerkin with Legendre Polynomials Harald Klimach, Michael Gaida, Sabine Roller harald.klimach@uni-siegen.de 26th WSSP 2017 Motivation Fluid-Dynamic
More informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
More informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationExponential Runge-Kutta methods for parabolic problems
Exponential Runge-Kutta methods for parabolic problems Marlis Hochbruck a,, Alexander Ostermann b a Mathematisches Institut, Heinrich-Heine Universität Düsseldorf, Universitätsstraße 1, D-4225 Düsseldorf,
More informationEvaluating Singular and Nearly Singular Integrals Numerically
Evaluating Singular and early Singular Integrals umerically Matthew J. Atwood Duke University Mentor: J. Thomas Beale This work is supported by SF VIGRE grant number DMS-998332. Abstract We look for a
More informationk=1 = e λ λ k λk 1 k! (k + 1)! = λ. k=0
11. Poisson process. Poissonian White Boise. Telegraphic signal 11.1. Poisson process. First, recall a definition of Poisson random variable π. A random random variable π, valued in the set {, 1,, }, is
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationComposite Membranes. 8th August 2006
Composite Membranes Aldo Lopes UFMG, aldoelopes@hotmail.com Russell E. Howes, Brigham Young University rhowes@byu.edu Cynthia Shepherd, California State University - Northridge Cynthia.Shepherd.87@csun.edu
More information