State of the Art MHD Methods for Astrophysical Applications p.1/32
|
|
- April Lamb
- 6 years ago
- Views:
Transcription
1 State of the Art MHD Methods for Astrophysical Applications Scott C. Noble February 25, 2004 CTA, Physics Dept., UIUC State of the Art MHD Methods for Astrophysical Applications p.1/32
2 Plan of Attack Is NOT a survey of all possible methods! Highlight relatively new methods! Focus on tests and comparisons between the methods. Cover general idea of algorithms. Details are deferred to given references. State of the Art MHD Methods for Astrophysical Applications p.2/32
3 Outline Introduction to ideal Relativistic MagnetoHydroDynamics (RMHD). Two basic approaches: Non-conservative and Conservative Methods for preserving the Divergence Constraint Conclusions State of the Art MHD Methods for Astrophysical Applications p.3/32
4 Ideal Relativistic Magnetohydrodynamics Perfect fluid = Isotropic, inviscid, no thermal conductivity: T µν fluid = (P + ρ + ρɛ) uµ u ν + P g µν µ (ρu µ ) = 1 g µ ( gρu µ ) = 0 Electrically neutral, infinitely conducting fluid: T µν EM = F µα F ν α 1 4 gµν F αβ F αβ µ F µν = J ν µ F µν = 0, F µν = 1 2 ε µνκλf κλ Infinite conductance F µν u µ = 0 State of the Art MHD Methods for Astrophysical Applications p.4/32
5 Ideal RMHD (cont.) ν ν = = i t F µν u µ = 0 µ F µν = 0 µ F µν = J ν v j = u i /u t b µ F µν u ν = 1 2 εµνκλ F κλ u ν F µν = b µ u ν b ν u µ b µ u µ = 0 B i F it b t = B i u µ g iµ b i = 1 u t `Bi + b t u i t ` gb i = j ˆ g `Bj v i B i v j j ` gb j = 0 State of the Art MHD Methods for Astrophysical Applications p.5/32
6 Ideal RMHD (cont.) ν ν = = i t F µν u µ = 0 µ F µν = 0 µ F µν = J ν v j = u i /u t b µ F µν u ν = 1 2 εµνκλ F κλ u ν F µν = b µ u ν b ν u µ b µ u µ = 0 B i F it b t = B i u µ g iµ b i = 1 u t `Bi + b t u i t ` gb i = j ˆ g `Bj v i B i v j j ` gb j = 0 State of the Art MHD Methods for Astrophysical Applications p.5/32
7 Ideal RMHD (cont.) ν ν = i t F µν u µ = 0 µ F µν = 0 µ F µν = J ν v j = u i /u t b µ F µν u ν = 1 2 εµνκλ F κλ u ν F µν = b µ u ν b ν u µ b µ u µ = 0 B i F it b t = B i u µ g iµ b i = 1 u t `Bi + b t u i t ` gb i = j ˆ g `Bj v i B i v j j ` gb j = 0 State of the Art MHD Methods for Astrophysical Applications p.5/32
8 Ideal RMHD (cont.) ν ν = i t F µν u µ = 0 µ F µν = 0 µ F µν = J ν v j = u i /u t b µ F µν u ν = 1 2 εµνκλ F κλ u ν F µν = b µ u ν b ν u µ b µ u µ = 0 B i F it b t = B i u µ g iµ b i = 1 u t `Bi + b t u i t ` gb i = j ˆ g `Bj v i B i v j j ` gb j = 0 State of the Art MHD Methods for Astrophysical Applications p.5/32
9 Ideal RMHD (cont. II) T µν MHD = `ρ + ρɛ + P + b 2 u µ u ν + `P b2 δ µ ν b µ b ν µ T µ ν = 0 t g ρu t T t µ B i j g 6 4 ρu j T j µ B j v i B i v j = g T α βγ β αµ t ( gq) + j ` gf j = S(q, g µν ) j ` gb j = 0 State of the Art MHD Methods for Astrophysical Applications p.6/32
10 Non-conservative Methods t q + i f i = S(q, µ q) Experienced methods from hyrdro. simulations; Excellent for smooth flows; Computationally efficient; Use of upwind, monotonic methods for advection helps w/ shocks; Requires use of artificial viscosity, P P + Q; Still problems near shocks, especially as v 1; Norman and Winkler, in Astrophys. Radiation Hydro., 449 (1986) State of the Art MHD Methods for Astrophysical Applications p.7/32
11 ZEUS-like: FD/CT/MOC-CT Methods FD = Finite Difference, MOC = Method of Characteristics CT = Constrained Transport Wilson, in 7th Texas Symp. on Rel. Astrophys., [MHD, Schwarz. Accretion] Hawley, Smarr, Wilson, ApJ 277, 296 (1984). [Hydro. Eq.s] Hawley, Smarr, Wilson, ApJS 55, 211, (1984). [Hydro. Tests, Kerr Accretion] Evans and Hawley, ApJ 332, 659 (1988). [CT, 2D MHD, Schwarzschild Accretion] Stone and Norman, ApJS 80, 791 (1992). [ZEUS-2D, MOC-CT] ZEUS Code webpage, jstone/zeus.html Hawley, Stone, Comput. Phys. Commun. 89, 127 (1995). [MOC-CT] De Villiers and Hawley, ApJ 589, 458 (2003). [3D MHD, Kerr Accretion] State of the Art MHD Methods for Astrophysical Applications p.8/32
12 De Villiers and Hawley 2003 (DH) Non-conservative scheme: i t q + i hv i G(q, g µν ) = S(q, µ q, g µν, α g µν ) Evolution is broken down into several sub-steps: t q = L(q) MOC CT : Finds b µ q 1 = q 0 + tl 0 (q 0 ) : Transport L 0 i ( γv i H), H = {D, E, S j } q 2 = q 1 + tl 1 (q 1 ) : Source L 1 µ {W, P, g νλ, b ν }, q Staggered mesh: vectors at cell faces, scalars at cell centers; New version of MOC-CT for i B i = 0; Full 3D MHD in Boyer-Lindquist coordinates; L 0 found via monotonic upwind scheme (van Leer 1977); State of the Art MHD Methods for Astrophysical Applications p.9/32
13 Conservative Methods: Finite Volume t q + x f = S(q) q n+1 j q n j = t `F j+1/2 F j 1/2 n+1/2 + t S n+1/2 j n+1 q j n+1 F n+1/2 j 1/2 F n+1/2 j+1/2 n q j n State of the Art MHD Methods for Astrophysical Applications p.10/32
14 Conservative Methods: Godunov Schemes t q + i f i = S(q) Assume piecewise-constant data, q q State of the Art MHD Methods for Astrophysical Applications p.11/32
15 Conservative Methods: Godunov Schemes t q + i f i = S(q) Assume piecewise-constant data, q q For higher resolution, interpolate for q L and q R State of the Art MHD Methods for Astrophysical Applications p.11/32
16 Conservative Methods: Godunov Schemes t q + i f i = S(q) Assume piecewise-constant data, q q For higher resolution, interpolate for q L and q R Need to use monotonic (aka slope-limiting) interpolation; State of the Art MHD Methods for Astrophysical Applications p.11/32
17 Conservative Methods: Godunov Schemes t q + i f i = S(q) Assume piecewise-constant data, q q For higher resolution, interpolate for q L and q R Need to use monotonic (aka slope-limiting) interpolation; q L q R State of the Art MHD Methods for Astrophysical Applications p.12/32
18 Riemann Problem q L q R Find q(x, t) given t q + x f = 0 8 < q L for x < 0 q(x, t = 0) = : q R for x > 0 Riemann solution F n+1/2. State of the Art MHD Methods for Astrophysical Applications p.13/32
19 Hydrodynamic Riemann Solution Rarefaction Fan L * P *, v *, ρo Contact Discontinuity * R ρ * P, v *, o Shock v L L P L,, ρ o R R R P, v, ρo State of the Art MHD Methods for Astrophysical Applications p.14/32
20 MHD Riemann Solution F L A L S L C S R A R F R L q R q F L,R = Fast magnetosonic wave S L,R = Slow magnetosonic wave A L,R = Alfvenic wave C = Contact Discontinuity State of the Art MHD Methods for Astrophysical Applications p.15/32
21 Approximate Riemann Solvers: HLL t q + A i i q = 0 H L q* H R L q R q F n+1/2 j+1/2 = 8 >< >: 1 (H R H L ) f(q L ) ˆHR f L H L f R + H L H R `q R q L if 0 H L if H L 0 H R f(q R ) if H R 0 State of the Art MHD Methods for Astrophysical Applications p.16/32
22 Approximate Riemann Solvers: HLLC H* H L q L* R* q H R L q R q F n+1/2 j+1/2 = 8 >< >: f(q L ) if 0 H L f(q L ) + H `q L L q L if H L 0 H f(q R ) + H `q R R q R if H L 0 H f(q R ) if H R 0 State of the Art MHD Methods for Astrophysical Applications p.17/32
23 Magnetized Bondi Flow HARM DH 2003 State of the Art MHD Methods for Astrophysical Applications p.18/32
24 Gammie Inflow Steady-state, equatorial, magnetized inflow in Kerr inside marginally stable orbit; HARM B φ DH 2003 State of the Art MHD Methods for Astrophysical Applications p.19/32
25 Slow Shock Komissarov, MNRAS 303, 343 (1999). DH 2003 State of the Art MHD Methods for Astrophysical Applications p.20/32
26 Rarefaction Komissarov N = 150 DH 2003 N = 2048 State of the Art MHD Methods for Astrophysical Applications p.21/32
27 Qualities of Conservative GRMHD Schemes Improved behavior in shock tests Equivalent behavior in relativistic cases Better behavior for ultra-relativistic flows Still developmental Worse behavior in floored regions Worse behavior when P mag P Conservative Primitive variable calculation Slower code! State of the Art MHD Methods for Astrophysical Applications p.22/32
28 Other Schemes Del Zanna et al., AA 400, 397 (2003). 3rd-order CENO with HLL SRMHD Balsara, ApJS 132, 83 (2001). Roe-type Approx. Riemann Solver Interpolates characteristic variables SRMHD ATHENA, jstone/athena.html Roe-type Approx. Riemann Solver Not relativistic State of the Art MHD Methods for Astrophysical Applications p.23/32
29 Maintianing the Divergence Constraint Toth, J. Comp. Phys., 161, 605 (2000). TVD-Lax-Friedrich i B i 0 leads to: Non-perpendicular Lorentz forces to B i Inconsistency in MHD Instabilities and spurious effects E z = v x B y v y B x = f x CT : Evans and Hawley 1988 MOC-CT: Hawley and Stone 1995 MOC-CT2: De Villiers and Hawley 2003 State of the Art MHD Methods for Astrophysical Applications p.24/32
30 Flux/Field CT/CD for Conserv. Methods Field-CT: Interp. predicted B i and v i for E n+1/2 j+1 ; State of the Art MHD Methods for Astrophysical Applications p.25/32
31 Flux/Field CT/CD for Conserv. Methods Field-CT: Interp. predicted B i and v i for E n+1/2 j+1 ; Flux-CT: Interp. f i and f in for E n+1/2 j+1/2 ; State of the Art MHD Methods for Astrophysical Applications p.25/32
32 Flux/Field CT/CD for Conserv. Methods Field-CT: Interp. predicted B i and v i for E n+1/2 j+1 ; Flux-CT: Interp. f i and f in for E n+1/2 Field-CD: E n+1/2 j = 1 2 `E n j + Ej ; j+1/2 ; State of the Art MHD Methods for Astrophysical Applications p.25/32
33 Flux/Field CT/CD for Conserv. Methods Field-CT: Interp. predicted B i and v i for E n+1/2 j+1 ; Flux-CT: Interp. f i and f in for E n+1/2 Field-CD: E n+1/2 j = 1 2 `E n j Flux-CD: f i for E n+1/2 j ; + Ej ; j+1/2 ; State of the Art MHD Methods for Astrophysical Applications p.25/32
34 8-Wave Formulation K. G. Powell, P. L. Roe, et al., J. Comput. Phys. 154, 284 (1999). K. G. Powell, ICASE Report No , Langley, VA, MHD equations can be derived w/o assuming constraint source terms i B i ; With i B i source terms, conserves divergence constraint along flow lines: t ` i B i + `vj j i B i = 0 Supposed to preserve initially constrained data along flows; Any monopoles are transported away hopefully off the grid; Increases robustness, G. Tóth and D. Odstrěil, J. Comput. Phys. 128, 82 (1996). Formulation becomes non-conservative problems near strong shocks; State of the Art MHD Methods for Astrophysical Applications p.26/32
35 Divergence Cleaning Schemes J. U. Brackbill and D. C. Barnes, J. Comput. Phys. 35, 426 (1980). C.-D. Munz, J. Comput. Phys. 161, 484 (2000). A. Dedner, et al., J. Comput. Phys. 175, 645 (2002). S. S. Komissarov, Appendix C of astro-ph/ t B i j v i B j v j B i + g ij Ψ = 0 D(Ψ) + i B i = 0 t D(Ψ) = i i Ψ D(Ψ) = 0 i i Ψ = 0 Elliptic D(Ψ) = 1 p Ψ tψ = p i i Ψ Parabolic D(Ψ) = 1 h 2 t Ψ tt Ψ = h 2 i i Ψ Hyperbolic D(Ψ) = 1 h 2 t Ψ + 1 p 2 Ψ tt Ψ + h2 p 2 t Ψ = h 2 i i Ψ Mixed State of the Art MHD Methods for Astrophysical Applications p.27/32
36 Discretization Differences in Constraint State of the Art MHD Methods for Astrophysical Applications p.28/32
37 2D Shock Tube Test State of the Art MHD Methods for Astrophysical Applications p.29/32
38 2D Shock Tube Test State of the Art MHD Methods for Astrophysical Applications p.30/32
39 Orzag-Tang Vortex State of the Art MHD Methods for Astrophysical Applications p.31/32
40 State of the Art MHD Methods for Astrophysical Applications p.32/32
The RAMSES code and related techniques 2- MHD solvers
The RAMSES code and related techniques 2- MHD solvers Outline - The ideal MHD equations - Godunov method for 1D MHD equations - Ideal MHD in multiple dimensions - Cell-centered variables: divergence B
More informationApproximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD
Approximate Harten-Lax-Van Leer (HLL) Riemann Solvers for Relativistic hydrodynamics and MHD Andrea Mignone Collaborators: G. Bodo, M. Ugliano Dipartimento di Fisica Generale, Universita di Torino (Italy)
More informationMulti-D MHD and B = 0
CapSel DivB - 01 Multi-D MHD and B = 0 keppens@rijnh.nl multi-d MHD and MHD wave anisotropies dimensionality > 1 non-trivial B = 0 constraint even if satisfied exactly t = 0: can numerically generate B
More informationA Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes )
A Comparative Study of Divergence-Cleaning Techniques for Multi-Dimensional MHD Schemes ) Takahiro MIYOSHI and Kanya KUSANO 1) Hiroshima University, Higashi-Hiroshima 739-856, Japan 1) Nagoya University,
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationVarious Hydro Solvers in FLASH3
The Center for Astrophysical Thermonuclear Flashes Various Hydro Solvers in FLASH3 Dongwook Lee FLASH3 Tutorial June 22-23, 2009 An Advanced Simulation and Computing (ASC) Academic Strategic Alliances
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationThe PLUTO code for astrophysical gasdynamics
Mem. S.A.It. Suppl. Vol. 13, 67 c SAIt 009 Memorie della Supplementi The PLUTO code for astrophysical gasdynamics A. Mignone 1, 1 Dipartimento di Fisica Generale Amedeo Avogadro, Università degli Studi
More informationThe Center for Astrophysical Thermonuclear Flashes. FLASH Hydrodynamics
The Center for Astrophysical Thermonuclear Flashes FLASH Hydrodynamics Jonathan Dursi (CITA), Alan Calder (FLASH) B. Fryxell, T. Linde, A. Mignone, G. Wiers Many others! Mar 23, 2005 An Advanced Simulation
More informationPart 1: Numerical Modeling for Compressible Plasma Flows
Part 1: Numerical Modeling for Compressible Plasma Flows Dongwook Lee Applied Mathematics & Statistics University of California, Santa Cruz AMS 280C Seminar October 17, 2014 MIRA, BG/Q, Argonne National
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationComputational Astrophysics
16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing
More informationDivergence-free interpolation of vector fields from point values exact B = 0 in numerical simulations
Mon. Not. R. Astron. Soc. 413, L76 L80 (2011) doi:10.1111/j.1745-3933.2011.01037.x Divergence-free interpolation of vector fields from point values exact B = 0 in numerical simulations Colin P. McNally
More informationarxiv: v1 [astro-ph.im] 23 Feb 2011
Mon. Not. R. Astron. Soc. 000, 000 000 (0000) Printed 25 February 2011 (MN LATEX style file v2.2) Divergence-free Interpolation of Vector Fields From Point Values - Exact B = 0 in Numerical Simulations
More informationHARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics Charles F. Gammie 1,2,3 and Jonathan C. McKinney 2 Center for Theoretical Astrophysics, University of Illinois 1002 W. Green St.,
More informationAn efficient shock-capturing central-type scheme for multidimensional relativistic flows. II. Magnetohydrodynamics
A&A 400, 397 413 (2003) DOI: 10.1051/0004-6361:20021641 c ESO 2003 Astronomy & Astrophysics An efficient shock-capturing central-type scheme for multidimensional relativistic flows II. Magnetohydrodynamics
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationA Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows 4
draft of August 0, 998 A Divergence-Free Upwind Code for Multidimensional Magnetohydrodynamic Flows 4 Dongsu Ryu, Francesco Miniati,T.W.Jones, and Adam Frank 3 ABSTRACT A description is given for preserving
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationarxiv: v2 [math.na] 4 Mar 2014
Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics arxiv:1309.3344v [math.na] 4 Mar 014 Andrew J. Christlieb a, James A. Rossmanith
More informationLocally Divergence-Free Central Discontinuous Galerkin Methods for Ideal MHD Equations
Locally Divergence-Free Central Discontinuous Galerkin Methods for Ideal MHD Equations Sergey Yakovlev, Liwei Xu and Fengyan Li Abstract In this paper, we propose and numerically investigate a family of
More informationNumerical Analysis of the Ultrarelativistic and Magnetized Bondi Hoyle Problem
Numerical Analysis of the Ultrarelativistic and Magnetized Bondi Hoyle Problem by Andrew Jason Penner B.Sc., University of Manitoba, 2002 M.Sc., University of Manitoba, 2004 A THESIS SUBMITTED IN PARTIAL
More informationAn efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes
An efficient implementation of the divergence free constraint in a discontinuous Galerkin method for magnetohydrodynamics on unstructured meshes Christian Klingenberg, Frank Pörner, Yinhua Xia Abstract
More informationQuest for Osher-Type Riemann Solver
p./7 Quest for Osher-Type Riemann Solver for Ideal MHD Equations Project Proposal EMRAS 2 Hervé Guillard, arry Koren, oniface Nkonga IRM, July 9 August 27, 2 p.2/7 Wine guide to Riemann solvers Godunov:
More informationDivergence- and Curl-Preserving Prolongation and Restriction Formulas
Journal of Computational Physics 180, 736 750 00) doi:10.1006/jcph.00.710 Divergence- and Curl-Preserving Prolongation and Restriction Formulas G. Tóth and P. L. Roe University of Michigan, Ann Arbor,
More informationHLLC solver for ideal relativistic MHD
HLLC solver for ideal relativistic MHD V. Honkkila, P. Janhunen* University of Helsinki, Department of Physical Sciences, P.O.Box 64 FIN- 4, Helsinki, Finland *Also at Finnish Meteorological Institute,
More informationThe ECHO code: from classical MHD to GRMHD in dynamical spacetimes
The ECHO code: from classical MHD to GRMHD in dynamical spacetimes Luca Del Zanna Dipartimento di Fisica e Astronomia Università di Firenze Main collaborators: N. Bucciantini, O. Zanotti, S. Landi 09/09/2011
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationarxiv:astro-ph/ v1 28 Nov 2005
Relativistic MHD Simulations of Jets with Toroidal Magnetic Fields. arxiv:astro-ph/0511769v1 8 Nov 005 Andrea Mignone (mignone@to.astro.it), Silvano Massaglia and Gianluigi Bodo Dipartimento di Fisica
More informationABSTRACT. A new unsplit staggered mesh algorithm (USM) that solves multidimensional
ABSTRACT Title of dissertation: An Unsplit Staggered Mesh Scheme For Multidimensional Magnetohydrodynamics: A Staggered Dissipation-control Differencing Algorithm Dongwook Lee, Doctor of Philosophy, 2006
More informationHyperbolic Divergence Cleaning for the MHD Equations
Journal of Computational Physics 175, 645 673 (2002) doi:10.1006/jcph.2001.6961, available online at http://www.idealibrary.com on Hyperbolic Divergence Cleaning for the MHD Equations A. Dedner, F. Kemm,
More informationarxiv: v1 [astro-ph.ep] 30 Jun 2009
Astronomy & Astrophysics manuscript no. MS c ESO 2018 September 21, 2018 arxiv:0906.5516v1 [astro-ph.ep] 30 Jun 2009 A high-order Godunov scheme for global 3D MHD accretion disks simulations. I. The linear
More informationThe one-dimensional equations for the fluid dynamics of a gas can be written in conservation form as follows:
Topic 7 Fluid Dynamics Lecture The Riemann Problem and Shock Tube Problem A simple one dimensional model of a gas was introduced by G.A. Sod, J. Computational Physics 7, 1 (1978), to test various algorithms
More informationNUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
NUMERICAL SOLUTION OF HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS JOHN A. TRANGENSTEIN Department of Mathematics, Duke University Durham, NC 27708-0320 Ш CAMBRIDGE ЩР UNIVERSITY PRESS Contents 1 Introduction
More informationarxiv:astro-ph/ v1 17 May 2002
Astronomy & Astrophysics manuscript no. (will be inserted by hand later) An efficient shock-capturing central-type scheme for multidimensional relativistic flows arxiv:astro-ph/0205290v1 17 May 2002 I.
More informationNumerical Hydrodynamics in Special Relativity
Numerical Hydrodynamics in Special Relativity J. M. Martí Departamento de Astronomía y Astrofísica Universidad de Valencia 46100 Burjassot (Valencia), Spain martij@godunov.daa.uv.es E. Müller Max-Planck-Institut
More informationMATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS
K CHAPMAN & HALL/CRC Monographs and Surveys in Pure and Applied Mathematics I 18 MATHEMATICAL ASPECTS OF NUMERICAL SOLUTION OF HYPERBOLIC SYSTEMS ANDREI G. KULIKOVSKII NIKOLAI V. POGORELOV ANDREI YU. SEMENOV
More informationOutline of lectures. Adding more physics. Doing Better: source terms in CTU unsplit integrator
Grid-based methods for hydrodynamics, MHD, and radiation hydrodynamics. Outline of lectures (Four lectures) Lecture 1. Introduction to physics and numerics Lecture 2. Operator split (ZEUS-like) methods
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationExtremum-Preserving Limiters for MUSCL and PPM
arxiv:0903.400v [physics.comp-ph] 7 Mar 009 Extremum-Preserving Limiters for MUSCL and PPM Michael Sekora Program in Applied and Computational Mathematics, Princeton University Princeton, NJ 08540, USA
More informationFinite Volume for Fusion Simulations
Finite Volume for Fusion Simulations Elise Estibals, Hervé Guillard, Afeintou Sangam To cite this version: Elise Estibals, Hervé Guillard, Afeintou Sangam. Finite Volume for Fusion Simulations. Jorek Meeting
More informationarxiv: v3 [physics.comp-ph] 22 Mar 2018
Preprint 3rd March 08 Compiled using MNRAS LATEX style file v3.0 arxiv:504.0985v3 [physics.comp-ph] Mar 08 Fourth Order Accurate Finite Volume CWENO Scheme For Astrophysical MHD Problems Prabal Singh Verma,,4,
More informationNumerical solutions of magnetohydrodynamic equations
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 59, No., 11 DOI: 1.478/v1175-11-7-9 Numerical solutions of magnetohydrodynamic equations K. MURAWSKI Faculty of Mathematics, Physics
More informationRESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method
RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific
More informationPREPRINT November 5, 2007
AN UNSPLIT STAGGERED MESH SCHEME FOR MULTIDIMENSIONAL MAGNETOHYDRODYNAMICS WITH EFFICIENT DISSIPATION CONTROLS Dongwook Lee a and Anil Deane b, a ASC FLASH Center, University of Chicago, 5640 S. Ellis,
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationA constrained-transport magnetohydrodynamics algorithm with near-spectral resolution.
A constrained-transport magnetohydrodynamics algorithm with near-spectral resolution. Jason L. Maron 1 Mordecai-Mark Mac Low 2 arxiv:0709.1234v1 [astro-ph] 8 Sep 2007 Jeffrey S. Oishi 3 Department of Astrophysics,
More informationA Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws
A Central Rankine Hugoniot Solver for Hyperbolic Conservation Laws S. Jaisankar, S. V. Raghurama Rao CFD Centre, Department of Aerospace Engineering, Indian Institute of Science, Bangalore, India-56 Abstract
More informationMagnetized Binary Neutron Stars with Whisky
Magnetized Binary Neutron Stars with Whisky Bruno Giacomazzo (AEI, Potsdam, Germany) in collaboration with Luca Baiotti (Tokyo) and Luciano Rezzolla (AEI) Plan of the Talk Introduction The Whisky(MHD)
More informationWhat is ZEUS-3D? 1 Introduction. Contents
What is ZEUS-3D? David Clarke, December 2007 Institute for Computational Astrophysics Saint Mary s University Halifax, NS, B3H 3C3, Canada http://www.ica.smu.ca/zeus3d Contents 1. Introduction...............................
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationVector Potential Remap for 2d MHD
UCRL-ID-132123 Vector Potential Remap for 2d MHD P. W. Rambo October 13, 1998 Lawrence Livermore National Laboratory This is an informal report intended primarily for internal or limited external distribution.
More informationVariability of accreting black holes induced by shocks in low angular momentum flows
induced by shocks in low angular momentum flows Astronomical Institute of the CAS Prague, Czech Republic Cooperation with Agnieszka Janiuk, CFT PAN, Vladimír Karas ASU CAS 23.10.2017 Low angular momentum
More informationOn divergence-free reconstruction schemes for CED and MHD
On divergence-free reconstruction schemes for CED and MHD Praveen Chandrashekar praveen@math.tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore-560065, India
More informationAxisymmetric Shearing Box Models of Magnetized Disks
Axisymmetric Shearing Box Models of Magnetized Disks Xiaoyue Guan and Charles F. Gammie Astronomy Department, University of Illinois, Urbana, IL 61801 ABSTRACT The local model, or shearing box, has proven
More informationRecent Progress in FLASH: High-Energy- Density Physics Applications
The Center for Astrophysical Thermonuclear Flashes Recent Progress in FLASH: High-Energy- Density Physics Applications Dongwook Lee The Physics of Intracluster Medium: Theory & Computation University of
More informationMathematical foundations of relativistic hydrodynamics. José Antonio Font
Mathematical foundations of relativistic hydrodynamics José Antonio Font University of Valencia NewCompStar School 2017, University of Sofia Outline Lecture 1: Relativistic Hydrodynamics Lecture 2: Numerical
More informationGeneral Relativistic Magnetohydrodynamics Equations
General Relativistic Magnetohydrodynamics Equations Ioana Duţan IMPRS Students Workshop Berggasthof Lansegger, Austria, 28-31 Aug 2006 Ioana Duţan, GRMHD Equations p.1/22 Introductory remarks (General)
More informationArt Checklist. Journal Code: Article No: 6959
Art Checklist Journal Code: JCPH Article No: 6959 Disk Recd Disk Disk Usable Art # Y/N Format Y/N Remarks Fig. 1 Y PS Y Fig. 2 Y PS Y Fig. 3 Y PS Y Fig. 4 Y PS Y Fig. 5 Y PS Y Fig. 6 Y PS Y Fig. 7 Y PS
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationGlobal Simulations of Black Hole Accretion. John F. Hawley Department of Astronomy, University of Virginia
Global Simulations of Black Hole Accretion John F. Hawley Department of Astronomy, University of Virginia Collaborators and Acknowledgements Julian Krolik, Johns Hopkins University Scott Noble, JHU Jeremy
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationarxiv:astro-ph/ v1 27 Jan 2006
Mon. Not. R. Astron. Soc. 000, 7 005) Printed January 04 MN LATEX style file v.) An HLLC Solver for Relativistic Flows II. Magnetohydrodynamics A. Mignone and G. Bodo INAF Osservatorio Astronomico di Torino,
More informationGodunov methods in GANDALF
Godunov methods in GANDALF Stefan Heigl David Hubber Judith Ngoumou USM, LMU, München 28th October 2015 Why not just stick with SPH? SPH is perfectly adequate in many scenarios but can fail, or at least
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学,
RKK ( ) BHmag2012,, 2012.2.29 Outline Motivation and basis: Magnetic reconnection around astrophysical black holes Standard equations of resistive GRMHD Test calculations of resistive GRMHD A simulation
More informationarxiv:astro-ph/ v1 7 Oct 2003
1 On the Divergence-Free Condition in Godunov-Type Schemes for Ideal Magnetohydrodynamics: the Upwind Constrained Transport Method arxiv:astro-ph/0310183v1 7 Oct 2003 P. Londrillo INAF - Osservatorio di
More informationLES Simulations of Quiet Sun Magnetism
LES Simulations of Quiet Sun Magnetism Matthias Rempel HAO/NCAR Quiet sun magnetism Origin and spatial distribution of quiet sun field Small scale dynamo? Remnant field from large scale dynamo? Vögler,
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 informationHigh-order Godunov schemes for global 3D MHD simulations of accretion disks. I. Testing the linear growth of the magneto-rotational instability
A&A 516, A26 (2010) DOI: 10.1051/0004-6361/200912443 c ESO 2010 Astronomy & Astrophysics High-order Godunov schemes for global 3D MHD simulations of accretion disks I. Testing the linear growth of the
More informationAdvanced Methods for Numerical Fluid Dynamics and Heat Transfer (MVKN70)
2017-08-24 Advanced Methods for Numerical Fluid Dynamics and Heat Transfer (MVKN70) 1 Lectures and course plans Week 1, 2: RY: Dr. Rixin Yu: tel: 222 3814; e-mail: rixin.yu@energy.lth.se Week 3; BS: Prof.
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationSistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas
Sistemas Hiperbólicos no Lineales: Un Nuevo Método para Calcular Flujos Relativistas Pedro González-Casanova Henríquez Unidad de Investigación en Cómputo Aplicado DGSCA, UNAM Elvira Torondel y Ricard Garrido
More informationNumerical Hydrodynamics: A Primer. Wilhelm Kley Institut für Astronomie & Astrophysik & Kepler Center for Astro and Particle Physics Tübingen
A Primer Wilhelm Kley Institut für Astronomie & Astrophysik & Kepler Center for Astro and Particle Physics Tübingen Bad Honnef 2016 Hydrodynamics: Hydrodynamic Equations The Euler-Equations in conservative
More informationarxiv: v4 [physics.flu-dyn] 23 Mar 2019
High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows Man Long Wong a, Sanjiva K. Lele a,b arxiv:7.895v4 [physics.flu-dyn] 3 Mar 9
More informationLarge Time Step Scheme Behaviour with Different Entropy Fix
Proceedings of the Paistan Academy of Sciences: A. Physical and Computational Sciences 53 (): 13 4 (016) Copyright Paistan Academy of Sciences ISSN: 0377-969 (print), 306-1448 (online) Paistan Academy
More informationTHE RIEMANN FRAMEWORK AND HIGHER ORDER GODUNOV SCHEMES FOR PARALLEL, SELF-ADAPTIVE COMPUTATIONAL ASTROPHYSICS
RevMexAA (Serie de Conferencias), 9, 72 76 (2000) THE RIEMANN FRAMEWORK AND HIGHER ORDER GODUNOV SCHEMES FOR PARALLEL, SELF-ADAPTIVE COMPUTATIONAL ASTROPHYSICS Dinshaw Balsara N.C.S.A., University of Illinois
More informationFinite volume approximation of the relativistic Burgers equation on a Schwarzschild (anti-)de Sitter spacetime
Turkish Journal of Mathematics http:// journals. tubitak. gov. tr/ math/ Research Article Turk J Math 2017 41: 1027 1041 c TÜBİTAK doi:10.906/mat-1602-8 Finite volume approximation of the relativistic
More informationNUMERICAL METHODS IN ASTROPHYSICS An Introduction
-1 Series in Astronomy and Astrophysics NUMERICAL METHODS IN ASTROPHYSICS An Introduction Peter Bodenheimer University of California Santa Cruz, USA Gregory P. Laughlin University of California Santa Cruz,
More informationComparison of Approximate Riemann Solvers
Comparison of Approximate Riemann Solvers Charlotte Kong May 0 Department of Mathematics University of Reading Supervisor: Dr P Sweby A dissertation submitted in partial fulfilment of the requirement for
More informationDEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 26 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 26 DEVELOPMENT AND APPLICATION OF GENERALIZED MUSTA SCHEMES V.
More informationAMath 574 February 11, 2011
AMath 574 February 11, 2011 Today: Entropy conditions and functions Lax-Wendroff theorem Wednesday February 23: Nonlinear systems Reading: Chapter 13 R.J. LeVeque, University of Washington AMath 574, February
More informationChp 4: Non-linear Conservation Laws; the Scalar Case. By Prof. Dinshaw S. Balsara
Chp 4: Non-linear Conservation Laws; the Scalar Case By Prof. Dinshaw S. Balsara 1 4.1) Introduction We have seen that monotonicity preserving reconstruction and iemann solvers are essential building blocks
More informationPrimitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics
Primitive Variable Solvers for Conservative General Relativistic Magnetohydrodynamics Scott C. Noble Physics Department, University of Illinois, Urbana, IL 61801, U.S.A. scn@uiuc.edu and Charles F. Gammie
More informationGas Evolution Dynamics in Godunov-type Schemes and Analysis of Numerical Shock Instability
NASA/CR-1999-208985 ICASE Report No. 99-6 Gas Evolution Dynamics in Godunov-type Schemes and Analysis of Numerical Shock Instability Kun Xu ICASE, Hampton, Virginia Institute for Computer Applications
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationarxiv: v1 [astro-ph] 10 Nov 2008
Mon. Not. R. Astron. Soc. 000, 1 15 (2007) Printed 10 November 2008 (MN LATEX style file v2.2) A five-wave HLL Riemann solver for relativistic MHD arxiv:0811.1483v1 [astro-ph] 10 Nov 2008 A. Mignone 1,2,
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
More informationA constrained transport scheme for MHD on unstructured static and moving meshes
doi:10.1093/mnras/stu865 A constrained transport scheme for MHD on unstructured static and moving meshes Philip Mocz, 1 Mark Vogelsberger 2 and Lars Hernquist 1 1 Harvard Smithsonian Center for Astrophysics,
More informationA switch to reduce resistivity in smoothed particle magnetohydrodynamics
MNRAS 436, 2810 2817 (2013) Advance Access publication 2013 October 11 doi:10.1093/mnras/stt1776 A switch to reduce resistivity in smoothed particle magnetohydrodynamics Terrence S. Tricco and Daniel J.
More informationComputational Fluid Dynamics-1(CFDI)
بسمه تعالی درس دینامیک سیالات محاسباتی 1 دوره کارشناسی ارشد دانشکده مهندسی مکانیک دانشگاه صنعتی خواجه نصیر الدین طوسی Computational Fluid Dynamics-1(CFDI) Course outlines: Part I A brief introduction to
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationarxiv: v4 [astro-ph.im] 27 Apr 2011
A ROBUST NUMERICAL SCHEME FOR HIGHLY COMPRESSIBLE MAGNETOHYDRODYNAMICS: NONLINEAR STABILITY, IMPLEMENTATION AND TESTS K. WAAGAN, C. FEDERRATH, AND C. KLINGENBERG arxiv:1101.3007v4 [astro-ph.im] 27 Apr
More informationA simple robust and accurate a posteriori subcell finite
A simple robust and accurate a posteriori subcell finite volume limiter for the discontinuous Galerkin method Michael Dumbser ERC Grant STiMulUs FET-HPC Grant ExaHyPE 1/ 69 Séminaire LJLL, Paris, France,
More informationA Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1
A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2017 Randall J. LeVeque Applied Mathematics University of Washington January 4, 2017 http://faculty.washington.edu/rjl/classes/am574w2017
More information