Blo lowin wing u up a a p polytr lytropic sta ic star Justin - - PowerPoint PPT Presentation

▶
blo lowin wing u up a a p polytr lytropic sta ic star
SMART_READER_LITE
LIVE PREVIEW

Blo lowin wing u up a a p polytr lytropic sta ic star Justin - - PowerPoint PPT Presentation

Nov 22, 2023 36 likes •250 views

Blo lowin wing u up a a p polytr lytropic sta ic star Justin Lietz Code Justin Lietz Code Kaitlin Cook Kaitlin Cook Analysis Analysis Sherwood Richers Inputs Sherwood Richers Inputs Motivations Toy model of

slide-1
SLIDE 1

Blo lowin wing u up a a p polytr lytropic sta ic star

Justin Lietz – Code Justin Lietz – Code Kaitlin Cook – Kaitlin Cook – Analysis Analysis Sherwood Richers – Inputs Sherwood Richers – Inputs

slide-2
SLIDE 2

Motivations

  • Toy model of a star + toy model of a shock

= toy supernova.

  • We can then follow the nucleosynthesis.
  • Acts as a good way to combine the tools from

last week.

slide-3
SLIDE 3

Aims

Lane Emden star Lagrangian VH1 Sedov blast Dimensionalise XNet

slide-4
SLIDE 4

Key Challenges

  • Working out units
  • Recasting VH1 in Lagrangian coordinates
  • Writing a wrapper script to take output from

VH1 to XNet

  • Displaying the results in a meaningful way
slide-5
SLIDE 5

Polytropic Stars: Lane-Emden Solutions

Poisson Equation: Polytropic EOS: Lane-Emden: Plot of a solution to Lane-Emden n = 3/2 Hydrostatic Equilibrium: Energy of the Blast = 1049 Ergs

slide-6
SLIDE 6

Sedov Blasts

  • The Sedov blast - a

large amount of energy in a small region.

  • Very high pressure in

the first zone, smaller pressures in all of the

  • ther zones (from star

model)

slide-7
SLIDE 7

Lagrangian Coordinates

  • Eulerian coordinates – zones transfer matter
  • Lagrangian coordinates – zones follow matter

→ Follow nucleosynthesis without advecting species

slide-8
SLIDE 8

VH1 Lagrangian

  • ppmlr.f90

Setup Eulerian Grid Map to Lagrangian Evolve Remap to Eulerian Map to Lagrangian Evolve Remap to Eulerian Setup Lagrangian Grid Evolve

slide-9
SLIDE 9

Lagrangian vs Eulerian

slide-10
SLIDE 10

Connecting VH1 to Xnet

  • VH1 output looks familiar to a lot of you by now:
slide-11
SLIDE 11

Hydro Results – 1 Msol + 1049 ergs

slide-12
SLIDE 12

10 Msol + 1049 ergs

slide-13
SLIDE 13

1 Msol Blast at different radii

slide-14
SLIDE 14

Nucleosynthesis 1 Msol + 1049 ergs 0.03 Rsol – Peak T = 3.3 GK

slide-15
SLIDE 15

Nucleosynthesis 1 Msol + 1049 ergs 0.10 Rsol– Peak T = 0.04 GK

slide-16
SLIDE 16

Nucleosynthesis 1 Msol + 1049 ergs 0.17 Rsol – Peak T = 0.001 GK

slide-17
SLIDE 17

50/50 n/p progenitor

slide-18
SLIDE 18

4He progenitor

slide-19
SLIDE 19

56Fe progenitor

slide-20
SLIDE 20

Abundances at t = 105

56Fe initial 4He initial

n+p initial

slide-21
SLIDE 21

With a few more days?

  • Add in gravity
  • Realistic EOS
  • Realistic initial conditions for hydro
  • Further variation of initial abundances
  • Larger Network
  • Vary blast energy
  • Nucleosynthesis on all zones
  • Feed the Nucleosynthesis energy into the hydro