hydrodynamics code CMacIonize Summary and recent developments Bert - - PowerPoint PPT Presentation

β–Ά
hydrodynamics code cmacionize
SMART_READER_LITE
LIVE PREVIEW

hydrodynamics code CMacIonize Summary and recent developments Bert - - PowerPoint PPT Presentation

Jan 22, 2023 461 likes β€’642 views

The Monte Carlo radiation hydrodynamics code CMacIonize Summary and recent developments Bert Vandenbroucke (bv7@st-andrews.ac.uk) Kenneth Wood (St Andrews) Nina Sartorio (So Jos dos Campos) Slide 2 of 16 Photoionization WHAM survey (2003)

slide-1
SLIDE 1

The Monte Carlo radiation hydrodynamics code CMacIonize

Summary and recent developments

Bert Vandenbroucke (bv7@st-andrews.ac.uk) Kenneth Wood (St Andrews) Nina Sartorio (SΓ£o JosΓ© dos Campos)

slide-2
SLIDE 2

Photoionization

NGC 604

HST, 1995

WHAM survey (2003) Slide 2 of 16

slide-3
SLIDE 3

The problem

π‘œπ‘Œπ‘—π½π‘Œπ‘—β†’π‘Œπ‘—+1 = π‘œπ‘Œπ‘—+1π‘œπ‘“π›½π‘Œπ‘—+1β†’π‘Œπ‘— π‘ˆ

π½π‘Œπ‘—β†’π‘Œπ‘—+1 is the total ionizing intensity π›½π‘Œπ‘—+1β†’π‘Œπ‘— π‘ˆ is the total recombination rate We need to solve to obtain equilibrium temperatures and ion densities These equations are highly non-linear; the ionizing intensity is non-local

Slide 3 of 16

slide-4
SLIDE 4

The solution

Regular (or AMR) grid Voronoi grid

𝐽𝑗 = ΰΆ±

πœ‰π‘’ ∞ 4πœŒπΎπ‘— πœ‰

β„Žπœ‰ 𝜏 πœ‰ dπœ‰ β‰ˆ 𝑅 𝑋

π‘’π‘Š 𝑗

෍

π‘˜

π‘₯

π‘˜π‘šπ‘—π‘˜πœ πœ‰π‘˜

Vandenbroucke & Wood (2018)

Slide 4 of 16

slide-5
SLIDE 5

The solution (2)

Iteratively repeat until convergence (takes approx. 10 iterations)

Slide 5 of 16

slide-6
SLIDE 6

The design

Modularity:

  • different grid types
  • different source types
  • different spectra
  • different physics
  • ...

User-friendliness:

  • explicit use of units
  • class structure maps to

parameter file

  • extensive documentation
  • library functionality

Reproducibility:

  • (automated) unit tests
  • extensive logging
  • benchmarks
  • publicly available

Vandenbroucke & Wood (2018) https://github.com/bwvdnbro/CMacIonize

Slide 6 of 16

slide-7
SLIDE 7

Example: turbulent HII regions

Uniform HII region Turbulent HII region

based on Wood et al. (2013), Vandenbroucke & Sartorio (in prep.) Slide 7 of 16

slide-8
SLIDE 8

Example: spectrum tracking

No diffuse All radiation

based on Wood & Mathis (2004), Vandenbroucke & Sartorio (in prep.)

Spectral hardening at work

Slide 8 of 16

slide-9
SLIDE 9

Example: escape fraction

Vandenbroucke & Sartorio, in prep.

Radiation escaping from turbulent HII clouds is softer than radiation escaping from a uniform cloud

Slide 9 of 16

slide-10
SLIDE 10

Example: SILCC(2) HΞ± emission

Vandenbroucke et al. (2018) Slide 10 of 16

slide-11
SLIDE 11

RHD simulations

Do hydro step Convert density field to number density field Do Monte Carlo photoionization Update pressure based

  • n ionization

structure

π‘ˆ = π‘ˆ

π‘œπ‘¦πΌ + π‘ˆπ‘— 1 βˆ’ 𝑦𝐼

π‘ˆ

π‘œ = 500 K

π‘ˆπ‘— = 8,000 K

Vandenbroucke & Wood (2018)

Harder than it sounds!

Slide 11 of 16

slide-12
SLIDE 12

Example: photoionization in DIG

Ionising luminosity per source

Vandenbroucke & Wood (submitted to MNRAS)

Slide 12 of 16

slide-13
SLIDE 13

Example: ionized accretion torus

Sartorio et al. (2019)

Slide 13 of 16

slide-14
SLIDE 14

Example: turbulent clouds

Sartorio & Vandenbroucke, in prep.

Slide 14 of 16

slide-15
SLIDE 15

The future: PMacIonize

Vandenbroucke & Wood (2018)

CMacIonize 2.0

Slide 15 of 16

slide-16
SLIDE 16

Summary

  • Do you have sources of ionising radiation?
  • Do you want to know how they ionise their

surroundings?

  • Do you want to know how ionisation affects

the dynamics of their surroundings?

  • Consider using !

https://github.com/bwvdnbro/CMacIonize Vandenbroucke & Wood (2018)

Slide 16 of 16