Coupling Solar Simulations with Space Weather: Code Comparison in - - PowerPoint PPT Presentation

▶
coupling solar simulations with space weather code
SMART_READER_LITE
LIVE PREVIEW

Coupling Solar Simulations with Space Weather: Code Comparison in - - PowerPoint PPT Presentation

Dec 20, 2022 506 likes •723 views

Coupling Solar Simulations with Space Weather: Code Comparison in SWIFF Consortium Vyacheslav Olshevsky 1 Center for Plasma Astrophysics, Mathematics Department, KU Leuven 2 Main Astronomical Observatory, Kyiv, Ukraine

slide-1
SLIDE 1

Coupling Solar Simulations with Space Weather: Code Comparison in SWIFF Consortium

Vyacheslav Olshevsky

1Center for Plasma Astrophysics, Mathematics Department, KU Leuven 2Main Astronomical Observatory, Kyiv, Ukraine

Vyacheslav.Olshevsky@wis.kuleuven.be Co-authors: Anna Lisa Restante1, Giovanni Lapenta1 and Rony Keppens1

slide-2
SLIDE 2

Introduction

20.06.2012 First results of SWIFF 2

SWIFF – Space Weather Integrated Forecasting Framework.

The principal goal of SWIFF is development of mathematical models and computational methods to handle multiple physics and multiple scales of space weather.

slide-3
SLIDE 3

Multi-physics in space weather

20.06.2012 First results of SWIFF 3

slide-4
SLIDE 4

Multi-scale in space weather

20.06.2012 First results of SWIFF 4

slide-5
SLIDE 5

Participants

20.06.2012 First results of SWIFF 5

5 scientific and 1 management work packages (WPs):

  • WP1. Coordination: Giovanni Lapenta, Katholieke

Universiteit Leuven (KU Leuven, Belgium)

  • WP2. Multiscale-Multiphysics Modelling: Rony Keppens,

Katholieke Universiteit Leuven (KU Leuven, Belgium)

  • WP3. Coupling at the Sun: Åke Nordlund, University of

Copenhagen (UCPH, Denmark)

  • WP4. Coupling in space: Francesco Califano, Università di

Pisa (UNIPISA, Italy)

  • WP5. Coupling at the Earth: Viviene Pierrard, Belgian

Institute for Space Aeronomy (BISA, Belgium)

slide-6
SLIDE 6

20.06.2012 First results of SWIFF 6

slide-7
SLIDE 7

Numerical Codes

20.06.2012 First results of SWIFF 7

  • FlipMHD – viscous resistive MHD code

(Brackbill, 1991)

  • MPI-AMRVAC – MPI implementation of the

Versatile Advection Code with adaptive Mesh Refinement (Keppens et al 2012)

  • iPic3D (KU Leuven)
  • Stagger (UCPH, Denmark)
  • UNIPI two-fluid code
  • ASI hybrid PIC code
slide-8
SLIDE 8

Benchmarks

20.06.2012 First results of SWIFF 8

  • Transition to turbulent reconnection (2D)
  • Longcope-Strauss problem (2D)
  • Magnetopause challenge (2D)
  • CME initiation challenge (3D)
slide-9
SLIDE 9

Flux Rope Emergence Setup

20.06.2012 First results of SWIFF 9

𝑪 = 𝜶 × 𝐵 𝑠, 𝜄 𝑠 sin 𝜄 𝜒 + 𝐶𝜒 𝑠, 𝜄 𝜒 , A 𝑠, 𝜄 =

1 2 𝑟𝑏2𝐶𝑢exp − 𝜕 2 𝑠,𝜄 𝑏2

, 𝐶𝜒 𝑠, 𝜄 = 𝑏𝐶𝑢 𝑠 sin 𝜄 exp − 𝜕 2 𝑠, 𝜄 𝑏2 .

Semi-circular rope (Fan & Gibson 2004) Embedded to the background sheared arcade (Aschwanden 2004)

𝐶𝑦 = 𝐶𝑦0 sin 𝑙𝑙 exp −𝑚𝑚 , 𝐶𝑧 = 𝐶𝑧0 sin 𝑙𝑙 exp −𝑚𝑚 , 𝐶𝑨 = 𝐶𝑨0 cos 𝑙𝑙 −𝑚𝑚 , 𝐶𝑦0 = 𝑚 𝑙 𝐶0, 𝐶𝑧0 = 𝛽 𝑙 𝐶0, 𝑙2 − 𝑚2 − 𝛽2 = 0.

slide-10
SLIDE 10

Setup Parameters

20.06.2012 First results of SWIFF 10

𝐶0 = 𝑞0 = 𝜍0 = 1, 𝑣0 = 0, 𝐶𝑢 = 9𝐶0, 𝑆 = 0.375, 𝑟 = −1, 𝑏 = 0.1, 𝑀𝑦 = 1.5, 𝑀𝑧 = 1, 𝑀𝑨 = 1.25, 𝑜𝑦 = 30, 𝑜𝑧 = 20, 𝑜𝑨 = 25. Boundary conditions: X – open Y – periodic Top – open Bottom – constant Bz

slide-11
SLIDE 11

Results: FlipMHD

20.06.2012 First results of SWIFF 11

Viscous No viscosity

Blue represent magnetic field lines; red contour is at constant density 0.8; The grayscale on the bottom is Bz. Note the destruction

  • f density “tube” in the second case!
slide-12
SLIDE 12

FlipMHD vs MPI-AMRVAC

20.06.2012 First results of SWIFF 12

FlipMHD (no viscosity) MPI-AMRVAC

Blue represent magnetic field lines; red contour is at constant density 0.8; The grayscale on the bottom is Bz. In both cases, the tube is destroyed. The behavior of field lines is very similar.

slide-13
SLIDE 13

Density Variation in XZ Plane

20.06.2012 First results of SWIFF 13

slide-14
SLIDE 14

Results: rise of the rope

20.06.2012 First results of SWIFF 14

slide-15
SLIDE 15

Conclusions

20.06.2012 First results of SWIFF 15

The initial configuration leads to immediate emergence of the flux rope, and is handy for comparison of numerical codes The overall behavior is qualitatively the same in FlipMHD and MPI-AMRVAC simulations Rise and rotation of the flux rope is faster in the simulations of MPI-AMRVAC

slide-16
SLIDE 16

Future? Particle-in-cell (PIC)

20.06.2012 First results of SWIFF 16

slide-17
SLIDE 17

Equations of implicit moment PIC

20.06.2012 First results of SWIFF 17

Equations of motion (moments of Vlasov) Interpolation of fields to particles Coupling of fields and particles is nonlinear:

Markidis, Lapenta, Rizwan-uddin (2010), Lapenta (2012)

slide-18
SLIDE 18

iPic3D code

20.06.2012 First results of SWIFF 18

  • Implicit time-stepping (“implicit moment method”)
  • Generalized Minimal Residual (GMRes) solver for the coupled

equations of motion – Maxwell system

  • Divergence cleaning (Conjugate Gradient) – charge density

continuation assurance

Applications

  • Reconnection on the Sun; magnetosphere; anywhere
  • Experiments in plasma devices
  • Ionization of spacecrafts
  • PIC – hybrid and PIC – fluid coupling
  • Mass loss in “hot jupiters”?
slide-19
SLIDE 19

http://www.swiff.eu

20.06.2012 First results of SWIFF 19

slide-20
SLIDE 20

Annex 1. Timeline of the SWIFF project

20.06.2012 First results of SWIFF 20