570 likes | 687 Views
Numerical comparison between maximal and geodesic slice in Schwarzschild black hole. 2006. 04. 20 ~21 Hyung Won Lee , Kyung Yee Kim Inje Institute for Mathematical Sciences School of Computer Aided Science, Inje University. Contents. Introduction
E N D
Numerical comparison between maximal and geodesic slice in Schwarzschild black hole 2006. 04. 20 ~21 Hyung Won Lee , Kyung Yee Kim Inje Institute for Mathematical Sciences School of Computer Aided Science, Inje University
Contents • Introduction • Gauge Conditions* Geodesic slice * Maximal slice • The Initial Value Problem • The Boundary condition • The Finite difference Scheme • Results
Solving Einstein’s Equations Numerically To learn radiated gravitational energy at infinity To understand dynamical processes near singularity What is Numerical Relativity? Introduction
Introduction Component to need the numerical construction
Introduction • Schwarzschild Black hole: • Vacuum, Spherically symmetric • Stationary and Static solution • Time reversal symmetry • Invariant on time translation • Asymptotically flat solution • Geometric mass:
Introduction Schwarzschild black hole line element: Transformation to isotropic coordinate
Introduction Isotropic coordinate • Invariant isometry transformation: • Allowing us to apply the isometry boundary condition. • Identical under mapping from the outside of throat to the inside of throat w.r.t. . • The form of map J which identifies two asymptotically flat sheet through the throat (isometry boundary condition):
: shift vector Introduction The 3+1 ADM (Arnowitt, Deser and Misner) Formalism : Lapse function To determine the proper time measured by an observer falling normal to the slice:
C B A Introduction The 3+1 ADM (Arnowitt, Deser and Misner) Formalism
Introduction The 3+1 ADM (Arnowitt, Deser and Misner) Formalism The Extrinsic Curvature • To describe the curvature of 3D slice in the 4D space-time in which it is embedded. • To be defined via the Lie derivative of the 3-metric with respect to the future-pointing unit normal vector to the 3-surface.
Introduction The 3+1 ADM (Arnowitt, Deser and Misner) Formalism The Einstein Equation: • Vacuum Einstein Equation • Hamilton Constraint Equation: • Momentum Constraint Equation: • Evolution Equation :
Introduction The 3+1 ADM (Arnowitt, Deser and Misner) Formalism Constraint equation Evolution equation Definition of
Contents • Introduction • Gauge Conditions* Geodesic slice * Maximal slice • The Initial Value Problem • The Boundary condition • The Finite difference Scheme • Results
C B A The Gauge Conditions
The Gauge Condition :Geodesic slice geodesic slice : radially infalling free particle t( coordinate time) r =0 r =2m r =r0
t=0 r =0 r =2m r =2m at crash t>0 r =2m t>0 t=0 0<r <2m r =2m at crash r =0 The Gauge Condition : Geodesic Slice geodesic slice:
t= The Gauge Condition : Maximal Slice • Maximal slice: t=0 r =2m r =0 r =2m t= t>0 r =2m t>0 t=0 0<r <2m r =2m r =3/2m
The Gauge Condition : Maximal Slice • Condition:
Contents • Introduction • Gauge Conditions* Geodesic slice * Maximal slice • The Initial Value Problem • The Boundary condition • The Finite difference Scheme • Results
To be determined 3-metric and to satisfying the constraint equation. • Using York’s conformal decomposition method • Insert 3-metric and extrinsic curvature to determine in Hamilton constraint equation : The Initial Value problem
The Initial Value problem Metric: Spherical symmetry line of element 3- metric: York’s conformal decomposition method
The Initial Value problem four dynamic variable two gauge variable
Contents • Introduction • Gauge Conditions* Geodesic slice * Maximal slice • The Initial Value Problem • The Boundary condition • The Finite difference Scheme • Results
The Boundary Condition The Inner boundary condition: • Black hole ‘s throat • Inner- most surface is locally an areal minimum. • For Schwarzschild black hole, the throat is the event horizon. • Isometry boundary condition. • Invariant :
The Boundary Condition The outer boundary condition: • far away from the black hole • Asymptotically flat at infinite region. The Inner boundary condition The Outer boundary condition
Contents • Introduction • Gauge Conditions* Geodesic slice * Maximal slice • The Initial Value Problem • The Boundary condition • The Finite difference Scheme • Results
The Numerical Algorithm Finite Difference method:
The Numerical Algorithm: Evolution Equation: • ½ time slice estimation of :
The Numerical Algorithm Grid:
The Numerical Algorithm: Evolution Equation for geodesic slice:
The Numerical Algorithm: To apply the finite difference method
The Numerical Algorithm: To apply the finite difference method
The Numerical Algorithm: Runge - Kutta method, shooting method
Contents • Introduction • Gauge Conditions* Geodesic slice * Maximal slice • The Initial Value Problem • The Boundary condition • The Finite difference Scheme • Results
Result Isotropic gauge: Static solution Code test result: static solution
Result T=0 T= 400(4M)
Result: geodesic slice The result Crashed around time:
Result Horizon [figure] metric a T = 100, 150, 200, 250
Result [Figure] Solution of Schwarzschild black hole
Result Schwarzchild Black hole Simulation Result Fig. a-variation with time
Result Schwarzchild Black hole Simulation Result Fig. b-variation with time
Result Schwarzchild Black hole Simulation Result Fig. -variation with time
Result Schwarzchild Black hole Simulation Result Fig. - variation with time
Result Schwarzchild Black hole Simulation Result
Result: Maximal slice [0M,5M,10M,20M,30M,40M,50M,60M,80M,90M] [gr-qc0307036 Fig9.Result]
Result: Maximal slice [0M,5M,10M,20M,30M,40M,50M,60M,80M,90M]