![Page 1: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/1.jpg)
NUMERICAL GENERATION OF VECTOR POTENTIALS FOR GENERAL
RELATIVITY SIMULATIONSZACHARY SILBERMAN
MAGNETIC FIELDS IN THE UNIVERSE VI
OCTOBER 19, 2017Collaborators:Joshua Faber (advisor), RITZachariah Etienne, West Virginia UniversityIan Ruchlin, West Virginia UniversityThomas Adams, West Virginia University
![Page 2: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/2.jpg)
NUMERICAL RELATIVITY
HARM3D
• Initial inspiral
• Approximations
• Evolves the magnetic field B
IllinoisGRMHD
• Merger
• Full Relativity
• Evolves the magnetic vector potential A
2
Gammie, McKinney, & Tóth 2003,Etienne et al. 2012
![Page 3: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/3.jpg)
A STAGGERED CELL
3
• Different quantities are defined at different locations in the grid
• Cell centers, face centers, edge centers, vertices
![Page 4: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/4.jpg)
STAGGERED GRID
4
![Page 5: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/5.jpg)
CALCULATING THE CURL
• I use a centered finite-differencing stencil to calculate the curl:
• with equivalent expressions for the other components of the magnetic field
5
![Page 6: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/6.jpg)
CELL-BY-CELL GENERATION
• Calculate the 12 A-field values of a cell using
• the 6 B-field values in that cell
• any A-field values from previously determined cells
• Seek to maximize symmetry of solution
• In cases where multiple solutions exist, attempt to minimize cell-to-cell variation in A
6
![Page 7: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/7.jpg)
RIGHT NOW
• Interpolates the data into my staggered coordinates
• Removes any divergence in B, if any exists
• Builds the A field
• Transforms into Coulomb gauge:
• Performs smoothness tests
• Ensures that the curl of A is B
• It works!
7
![Page 8: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/8.jpg)
GLOBAL LINEAR ALGEBRA
• Treat the curl operator as a matrix
• Include Coulomb gauge conditions
• More symmetric than the cell-by-cell method
• Also works!
8
![Page 9: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/9.jpg)
ROTOR TEST
• Initial B-field along x
• Edge of disk moving at 0.995c
• t = 0.2
9
![Page 10: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/10.jpg)
ROTOR TEST
• Initial B-field along x
• Edge of disk moving at 0.995c
• t = 0.2
10
![Page 11: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/11.jpg)
ROTOR TEST
• Final time t=0.4
• Run 1:
• t=0 to t=0.4
• Run 2:
• t=0 to t=0.2
• B->A
• t=0.2 to t=0.4
11
![Page 12: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/12.jpg)
ROTOR TEST
• Final time t=0.4
• Run 1:
• t=0 to t=0.4
• Run 2:
• t=0 to t=0.2
• B->A
• t=0.2 to t=0.4
12
![Page 13: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/13.jpg)
TOLMAN-OPPENHEIMER-VOLKOFF EQUATION
• Structure of a body with:
• Spherical symmetry
• Isotropic material
• Static gravitational equilibrium
13
![Page 14: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/14.jpg)
TOV STARS:STABLE
14
• Final time: 4tdyn ≈ 11.2• Restart Time: tdyn ≈ 2.8
![Page 15: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/15.jpg)
TOV STARS:UNSTABLE
15
• Final time: 4tdyn ≈ 11.2• Restart Time: tdyn ≈ 2.8
![Page 16: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/16.jpg)
THE FUTURE
• More Evolution Tests
• Continue taking A-field data for other configurations and evolve in time using the numerical relativity codes
• Mesh Refinement
• Create multiple levels of data, where for example the grid spacing of one level is half that of another
• Use with HARM3D and IllinoisGRMHD
• Run full binary black hole simulations
16
![Page 17: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/17.jpg)
BONUS SLIDES
17
![Page 18: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/18.jpg)
ROTOR DENSITY
18
![Page 19: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/19.jpg)
TOV DENSITY:STABLE
19
![Page 20: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/20.jpg)
TOV DENSITY:UNSTABLE
20
![Page 21: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/21.jpg)
ROTOR VELOCITY
21
![Page 22: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/22.jpg)
TOV VELOCITY:STABLE
22
![Page 23: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/23.jpg)
TOV VELOCITY:UNSTABLE
23
![Page 24: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/24.jpg)
NO MONOPOLES
• The divergence-free nature of the magnetic field can be ensured by defining a vector potential. Maxwell’s equations imply:
• Getting from A to B is easy; getting from B to A is not!
24
![Page 25: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/25.jpg)
DIVERGENCE REMOVAL
• The divergence is calculated for a cell by finite differencing, using the centered stencil
• If this is non-zero, magnetic flux is symmetrically removed through the faces of the cell
• My A-field solver introduces no divergence by construction
25
Tóth 2000
![Page 26: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/26.jpg)
SINGLE CELL SYMMETRY CONSIDERATIONS
• Each edge takes values based on:
• The two neighboring faces
• The two opposite faces
• NOT the perpendicular faces!
26
![Page 27: NUMERICAL GENERATION OF VECTOR POTENTIALS FOR … · numerical generation of vector potentials for general relativity simulations zachary silberman magnetic fields in the universe](https://reader030.vdocument.in/reader030/viewer/2022041213/5e00d4d27e1f9d4e32686db9/html5/thumbnails/27.jpg)
GAUGE CONDITIONS
• My results are in an arbitrary gauge, dependent on the order I evaluate my cells
• Coulomb:
• There is a prescription to take an arbitrary gauge and produce a Coulomb gauge solution:
27