Multi-D MHD and B = 0
|
|
- Mitchell Dorsey
- 5 years ago
- Views:
Transcription
1 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 0 due to non-linearities of shock-capturing numerical methods B 0 build-up numerical instability (+ physical nonsense) need to control this somehow strategies for B = 0
2 CapSel DivB - 02 MHD wave signals locally (δ-function) perturb homogeneous magnetized plasma at rest take γ = 5/3, ρ = 1, p th = 0.6 and B = 0.9ê x (c = 1, b = 0.9) simulate on (x, y) [ 0.5, 0.5] 2 in 2.5D MHD (include v z, B z ) perturb at origin with δρ = 0.1, δv z = 0.01 and δp th = 0.06 MHD counterpart of throwing a stone in a puddle entropy, total pressure, B z at finite time
3 CapSel DivB waves: return in Friedrichs group diagram (theoretical one overplotted) group velocity for wave package (and energy flow direction) v gr = ω(k) k = ω ê x + ω ê y + ω ê z k x k y k z using wave frequency ω and wave vector k at origin: local entropy (density) disturbance remains (plasma at rest) slow magnetosonic wave signals: cusp-like features: anisotropic! fast magnetosonic waves: like sound waves (would be spherical signal) velocity c = γp/ρ along B [generally max(c, b = B/ ρ)] travel faster in perpendicular direction b 2 + c 2 Alfvén waves: extremely anisotropic: pure point-like along B travel at Alfvén speed b and sample magnetic field lines! note different polarizations: s only shows entropy wave, p shows slow and fast, Alfvén only visible in v z, B z
4 CapSel DivB - 04 B = 0 physically: only exact solenoidal field is allowed numerics: discretization error + machine precision unavoidable conservative form of momentum equation: uses divergence of Maxwell stress tensor equal to Lorentz force IF B = 0 since ) (I B2 2 B B = ( B ) B B ( B ) force orthogonal to B if solenoidal field would like discrete B = 0, orthogonal Lorentz force, conservative form difficult to satisfy all demands simultaneously but possible, cfr. Tóth JCP 182, 346 (2002)!!!
5 CapSel DivB - 05 Vector potential rewrite MHD equations in terms of vector potential A defined from B = A keeps B = 0 exactly analytically still need discrete (.) = 0 increases order of occuring spatial derivatives (accuracy loss) Boundary Conditions on A not always straightforward conflicts with current Roe-type solvers exploiting B
6 CapSel DivB - 06 Projection scheme should be sufficient to control numerical value of B enforce constraint in particular discretization to given accuracy Take corrective action by projection scheme [Brackbill & Barnes 1980] scheme yields B with B 0 correct to solenoidal B = B φ, solve 2 φ = B Poisson eqn., solve by (iterative) scheme up to desired accuracy projects B on subspace of zero divergence solutions no change in current density accuracy: does not need to be machine precision!
7 CapSel DivB - 07 does not violate conservation properties indirectly affects local thermal/magnetic balance in total energy same ideas in use for incompressible HD (where v = 0) uniform Cartesian grids: projection scheme makes smallest possible correction to remove divergence of B provided by base scheme given B, find closest field B from minimizing B B 2 under solenoidal constraint B = 0 does not affect order of accuracy of base scheme (remains same) proved by Tóth (2000): scheme is consistent for weak solutions Can use projection to eliminate finite divergence of discrete initial B
8 CapSel DivB - 08 Powell source terms solution is ok up to truncation error: same goes for B = 0 maintaining constraint to truncation error is sufficient system of 8 PDEs for ideal MHD equations in conservation form not Galilean invariant 8 8 system has eigenvalues zero, and v, v ± c s, v ± b, v ± c f spurious eigenvalue which conflicts with constraint B = 0 carries a jump in normal (to cell edge) component of B restore Galilean invariance by writing system as U t + ( F) = S sources S = (S ρ, S ρv, S e, S B ) proportional to B
9 CapSel DivB - 09 Powell suggests following sources: 0 ( B )B ( B )B v ( B )v restore Galilean invariance, replace zero with extra v eigenvalue introduces 8-wave approximate Riemann solver induction equation with source term B + (vb B v) = v( B ) t equivalent to evolution equation for B given by B + (v( B )) = 0 t passively convected scalar field B /ρ
10 CapSel DivB - 10 idea is that numerical divergence errors get advected away potential problems at internal flow stagnation points sources destroy conservation form potential violation of Rankine Hugoniot across shock case with incorrect jumps across discontinuities in Tóth (2000) usually does work correctly though! can be used with non-riemann solver based method (e.g TVDLF) restore conservation of momentum/energy by only taking S B along arguments given by Janhunen (2000) and Dellar (2001) Lorentz (instead of Galilean) invariance makes momentum and energy equations conservative
11 CapSel DivB - 11 Constrained Transport enforce B = 0 in one particular discretization kept zero to machine precision in one discretization must have initial field with zero B in chosen discretization must have BCs compatible with zero B in chosen discretization typical CT approaches employ staggered magnetic field representation with ρ, e at cell center, take B at cell edges fluxes in induction eqn at cell vertex (corners), B cell centered Tóth (2000): CT recast in FV sense: no staggering needed! VAC: 5 different variants of CT/central difference type schemes Tóth (2000): two new finite volume CD type schemes!
12 CapSel DivB - 12 Evans & Hawley CT idea place magnetic field components on cell interfaces 2D: take electric field E z Ω = (v B ) z at cell corners Ω b y Ω b x.b b x Ω b y Ω update cell interface B = (B x, B y ) from induction equation as B x,n+1 j+1/2,k Bx,n j+1/2,k t B y,n+1 j,k+1/2 By,n j,k+1/2 t = Ω j+1/2,k+1/2 Ω j+1/2,k 1/2 y = + Ω j+1/2,k+1/2 Ω j 1/2,k+1/2 x
13 CapSel DivB - 13 then B n+1 = 0 if B n = 0 with ( B ) j,k = Bx j+1/2,k Bx j 1/2,k x + By j,k+1/2 By j,k 1/2 y up to accuracy of round-off errors (machine precision) note: CT does not ensure B = 0 in any other discretization!
14 CapSel DivB - 14 variations of CT idea without staggering use spatio-temporal interpolations from cell-centered v and B B y v xb B y v xb Ω b y Ω v B b x v B b x B x B x B B x Ω b y Ω v B v B v xb B v xb y Field interpolated CT (spatial + temporal interpolation) versus Field interpolated Central Difference (latter uses only temporal interpolation) can both be written as FV conservative update on cell-centered B involves averaged fluxes for B update with certain stencil
15 CapSel DivB - 15 Diffusive treatment add diffusion type source terms to energy/induction equation diffuse B errors at maximal rate maximal rate = allowed by unchanged CFL condition for parabolic equation B t = η D 2 B get CFL constraint t < x 2 /η D take η D x 2 and use sources diffusion coefficient C d B C d x 2 ( B ) C d x 2 ( B )
16 CapSel DivB - 16 source terms destroy strict conservation (like Powell) could retain source term for induction equation alone parabolic approach: errors are damped at maximal rate hyperbolic divergence cleaning strategy: Dedner et al. (2002) B transported at maximal admissable speed AND damped simultaneously
17 CapSel DivB - 17 Numerical tests taken from Tóth (2000): 7 strategies on 9 2D MHD tests fair comparison: same (approximate Riemann solver) base scheme best strategy is projection/field-cd/flux-ct scheme need a strategy for B, usually does not matter too much which
18 CapSel DivB D rotated shock tube: 1D Riemann problem in 2D rotated over angle must maintain constant B (parallel to shock tube direction) RH violation by Powell source terms
19 CapSel DivB - 19a Orszag-Tang (1979) MHD vortex simulation (Picone & Dahlburg 1991) 2D domain [0, 2π] 2 with double periodic sides t = 0 uniform ρ = 25/9, p = 5/3, velocity vortex v = ( sin y, sin x) magnetic islands horizontal wavelength B = ( sin y, sin 2x) with γ = 5/3: Mach 1 flow conditions mimics evolution to compressible (supersonic) MHD turbulence
20 CapSel Orszag-Tang vortex problem, temperature at t = 3.14 DivB - 19b
21 CapSel DivB - 20 References G. Tóth, J. Comp. Phys. 161, 605 (2000) K.G. Powell et al., J. Comp. Phys. 154, 284 (1999) P.J. Dellar, J. Comp. Phys. 172, 392 (2001) P. Janhunen, J. Comp. Phys. 160, 649 (2000) G. Tóth, J. Comp. Phys. 182, 346 (2002) J.U. Brackbill, D.C. Barnes, J. Comp. Phys. 35, 426 (1980) A. Dedner et al., J. Comp. Phys. 175, 645 (2002) W. Dai, P.R. Woodward, Astrophysical Journal 494, 317 (1998) C.R. Evans, J.F. Hawley, Astrophysical Journal 332, 659 (1988)
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 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 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 informationPHYS 643 Week 4: Compressible fluids Sound waves and shocks
PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations
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 informationCapSel Euler The Euler equations. conservation laws for 1D dynamics of compressible gas. = 0 m t + (m v + p) x
CapSel Euler - 01 The Euler equations keppens@rijnh.nl conservation laws for 1D dynamics of compressible gas ρ t + (ρ v) x = 0 m t + (m v + p) x = 0 e t + (e v + p v) x = 0 vector of conserved quantities
More informationState of the Art MHD Methods for Astrophysical Applications p.1/32
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 Plan of Attack Is
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 informationRecapitulation: Questions on Chaps. 1 and 2 #A
Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?
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 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 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 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 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 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 informationMagnetohydrodynamic Waves
Magnetohydrodynamic Waves Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 17, 2016 These slides are largely based off of 4.5 and 4.8 of The Physics of
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 informationJet Stability: A computational survey
Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,
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 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 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 informationConstrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics
Project work at the Department of Mathematics, TUHH Constrained Transport Method for the Finite Volume Evolution Galerkin Schemes with Application in Astrophysics Katja Baumbach April 4, 005 Supervisor:
More informationNumerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence
Numerical Study of Compressible Isothermal Magnetohydrodynamic Turbulence Junseong Park, Dongsu Ryu Dept. of Physics, Ulsan National Institute of Science and Technology. Ulsan, Korea 2016 KNAG meeting
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 informationSOLAR MHD Lecture 2 Plan
SOLAR MHD Lecture Plan Magnetostatic Equilibrium ü Structure of Magnetic Flux Tubes ü Force-free fields Waves in a homogenous magnetized medium ü Linearized wave equation ü Alfvén wave ü Magnetoacoustic
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 informationProf. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit
Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n
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 informationRadiative & Magnetohydrodynamic Shocks
Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,
More informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationLinear stability of MHD configurations
Linear stability of MHD configurations Rony Keppens Centre for mathematical Plasma Astrophysics KU Leuven Rony Keppens (KU Leuven) Linear MHD stability CHARM@ROB 2017 1 / 18 Ideal MHD configurations Interested
More informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationMaxwell's Equations and Conservation Laws
Maxwell's Equations and Conservation Laws 1 Reading: Jackson 6.1 through 6.4, 6.7 Ampère's Law, since identically. Although for magnetostatics, generally Maxwell suggested: Use Gauss's Law to rewrite continuity
More informationIntroduction to Magnetohydrodynamics (MHD)
Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality
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 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 informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationIncompressible MHD simulations
Incompressible MHD simulations Felix Spanier 1 Lehrstuhl für Astronomie Universität Würzburg Simulation methods in astrophysics Felix Spanier (Uni Würzburg) Simulation methods in astrophysics 1 / 20 Outline
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationWebster s horn model on Bernoulli flow
Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular
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 informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationHow can jets survive MHD instabilities?
How can jets survive MHD instabilities? Hubert Baty Observatoire Astronomique, 11 Rue de l université 67000 Strasbourg, France Rony Keppens FOM-Institute for Plasma Physics Rijnhuizen, Association Euratom/FOM,
More informationMacroscopic plasma description
Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationThe 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 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 informationComputational Fluid Dynamics. PHY 688: Numerical Methods for (Astro)Physics
Computational Fluid Dynamics Hydrodynamics When we discussed PDEs, we focused so far on scalar PDEs Often we wish to study systems of PDEs. Here we'll look at the equations of hydrodynamics Nonlinear system
More informationIdeal Magnetohydrodynamics (MHD)
Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics
More informationWaves and characteristics: Overview 5-1
Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables
More informationConstrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds Terrence S. Tricco a,b,, Daniel J. Price c, Matthew R. Bate b a Canadian Institute for
More informationConstrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds Terrence S. Tricco a,b,, Daniel J. Price c, Matthew R. Bate b a Canadian Institute for
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 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 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 information13. ASTROPHYSICAL GAS DYNAMICS AND MHD Hydrodynamics
1 13. ASTROPHYSICAL GAS DYNAMICS AND MHD 13.1. Hydrodynamics Astrophysical fluids are complex, with a number of different components: neutral atoms and molecules, ions, dust grains (often charged), and
More informationTwo-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation
Mathematics and Computers in Simulation 76 (2007) 3 7 Two-scale numerical solution of the electromagnetic two-fluid plasma-maxwell equations: Shock and soliton simulation S. Baboolal a,, R. Bharuthram
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 information2 The incompressible Kelvin-Helmholtz instability
Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect
More information20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:
Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir
More informationarxiv: v1 [physics.comp-ph] 16 Nov 2017
Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations Dominik Derigs a,, Andrew R. Winters b, Gregor J. Gassner b, Stefanie Walch
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
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 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 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 informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More informationMHD WAVES AND GLOBAL ALFVÉN EIGENMODES
MHD WVES ND GLOBL LFVÉN EIGENMODES S.E. Sharapov Euratom/CCFE Fusion ssociation, Culham Science Centre, bingdon, Oxfordshire OX14 3DB, UK S.E.Sharapov, Lecture 3, ustralian National University, Canberra,
More informationLinear and non-linear evolution of the gyroresonance instability in Cosmic Rays
Linear and non-linear evolution of the gyroresonance instability in Cosmic Rays DESY Summer Student Programme, 2016 Olga Lebiga Taras Shevchenko National University of Kyiv, Ukraine Supervisors Reinaldo
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationThe Physics of Collisionless Accretion Flows. Eliot Quataert (UC Berkeley)
The Physics of Collisionless Accretion Flows Eliot Quataert (UC Berkeley) Accretion Disks: Physical Picture Simple Consequences of Mass, Momentum, & Energy Conservation Matter Inspirals on Approximately
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 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 informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationSpace Physics. ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen. Aalto University School of Electrical Engineering
Space Physics ELEC-E4520 (5 cr) Teacher: Esa Kallio Assistant: Markku Alho and Riku Järvinen Aalto University School of Electrical Engineering The 6 th week: topics Last week: Examples of waves MHD: Examples
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationarxiv: v1 [astro-ph.im] 17 Nov 2009
A Second-Order Unsplit Godunov Scheme for Cell-Centered MHD: the CTU-GLM scheme. Andrea Mignone,a, Petros Tzeferacos a a Dipartimento di Fisica Generale, Universitá degli studi di Torino arxiv:0911.3410v1
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationThe MRI in a Collisionless Plasma
The MRI in a Collisionless Plasma Eliot Quataert (UC Berkeley) Collaborators: Prateek Sharma, Greg Hammett, Jim Stone Modes of Accretion thin disk: energy radiated away (relevant to star & planet formation,
More informationDiffusion / Parabolic Equations. PHY 688: Numerical Methods for (Astro)Physics
Diffusion / Parabolic Equations Summary of PDEs (so far...) Hyperbolic Think: advection Real, finite speed(s) at which information propagates carries changes in the solution Second-order explicit methods
More informationComputational Astrophysics 7 Hydrodynamics with source terms
Computational Astrophysics 7 Hydrodynamics with source terms Oscar Agertz Outline - Optically thin radiative hydrodynamics - Relaxation towards the diffusion limit - Hydrodynamics with gravity source term
More informationProblem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600
Problem C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 6 Overview This problem is aimed at testing the accuracy and the performance of high-order methods on the direct numerical simulation
More informationEfficient Low Dissipative High Order Schemes for Multiscale MHD Flows, II: Minimization of B Numerical Error
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln NASA Publications National Aeronautics and Space Administration 10-2006 Efficient Low Dissipative High Order Schemes for
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationPiecewise Parabolic Method on a Local Stencil for Magnetized Supersonic Turbulence Simulation
Piecewise Parabolic Method on a Local Stencil for Magnetized Supersonic Turbulence Simulation Sergey D. Ustyugov, Mikhail V. Popov, Keldysh Institute of Applied Mathematics, Miusskaya Sq. 4, 125047, Moscow,
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 informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More information2D compressible vortex sheets. Paolo Secchi
2D compressible vortex sheets Paolo Secchi Department of Mathematics Brescia University Joint work with J.F. Coulombel EVEQ 2008, International Summer School on Evolution Equations, Prague, Czech Republic,
More informationComputation of Incompressible Flows: SIMPLE and related Algorithms
Computation of Incompressible Flows: SIMPLE and related Algorithms Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de SIMPLE-Algorithm I - - - Consider
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 informationShocks. Ellen Zweibel.
Shocks Ellen Zweibel zweibel@astro.wisc.edu Departments of Astronomy & Physics University of Wisconsin, Madison and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas Shocks
More informationarxiv:astro-ph/ v1 6 May 2003
A Free, Fast, Simple and Efficient TVD MHD Code Ue-Li Pen, Phil Arras, ShingKwong Wong,3 Canadian Institute for Theoretical Astrophysics, University of Toronto Kavli Institute for Theoretical Physics 3
More informationSolar Physics & Space Plasma Research Center (SP 2 RC) MHD Waves
MHD Waves Robertus vfs Robertus@sheffield.ac.uk SP RC, School of Mathematics & Statistics, The (UK) What are MHD waves? How do we communicate in MHD? MHD is kind! MHD waves are propagating perturbations
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More 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 informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationThe importance of including XMHD physics in HED codes
The importance of including XMHD physics in HED codes Charles E. Seyler, Laboratory of Plasma Studies, School of Electrical and Computer Engineering, Cornell University Collaborators: Nat Hamlin (Cornell)
More information