Numerical Evolution of Strong Brill Waves
Dae-Il (Dale) Choi
USRA, NASA/GSFC
Collaborators:
Dr. D. Brown and Ms. L. Lowe (NCSU)
and
Members of Goddard Numerical Relativity Group
Fri Apr 4 17:24:38 EST 2003
Initial Data
- Generated using elliptic solver AMRMG developed by Brown and Lowe.
Feb. 21, 2003 Presentation [PDF] by David Brown
Refer to a following webpage for various initial value problems in GR to study using AMRMG code.
- ADM Mass for various amplitude M_ADM
- \Psi(x,y,z=0) for A=5,h=0.125
- Convergence factor(x,y,z=0) for Octant symm. runs with A=1, h=0.25,0.125,0.0625
- Convergence factor(x,y,z=0) for full grid runs with A=1, h=0.5,0.25,0.125
- Convergence factor(x,y,z=0) for full grid runs with A=1, h=0.25,0.125,0.0625
- Conformal factor(x,y,z=0) for full grid run with A=1, h=0.5 (Low resol)
- Conformal factor(x,y,z=0) for full grid run with A=1, h=0.125 (High resol)
-
- Evolution of Strong Teukolsky wave (A=0.001) as a evolution code test: metric functions are 2nd-order convergent (e.g. gzz [PS])
Evolution
- Unigrid---- Brill parameter [A=3, M_ADM=0.26]
- Code is 2nd-order convergent! (Hamiltonian Constraint Violation is 2nd-order convergent (e.g. gzz [PDF])
- Lapse function (ln(alpha) at x=y=z=0) vs time [PDF] -- disperses away
- Unigrid---- Brill parameter [A=4.7, M_ADM=]
- Unigrid---- Brill parameter [A=4.75, M_ADM=]
- Lapse function vs x (9 time shot) [PDF]
- Unigrid---- Brill parameter [A=4.8, M_ADM=0.643]
- Lapse function vs x (in 9 time shot) [PDF]
- Lapse function(0,0,0) vs Time [PDF]
- Fixed Mesh Refinement---- Brill parameter [A=3, M_ADM=0.26]
- Adaptive Mesh Refinement---- First tests
[A=4MPEG]
[A=4.7MPEG]
[A=4.75MPEG]
[A=4.8MPEG]
[A=5(solid)MPEG]
[A=5(wire)MPEG]
[A=4,4.7,4.8 (lapse at the origin)JPEG]
Talks
- Presentation at APS April Meeting (Philadelphia, PA) APR 8, 2003 [PS]
- FIGURES for the post presentation at AGWS (College Park, MD) APR 24-26, 2003
References
- Matt's IMA NR workshop Minneapolis, MN June 26, 2002 talk