scholarly journals Constraint damping in first-order evolution systems for numerical relativity

2007 ◽  
Vol 76 (4) ◽  
Author(s):  
Robert Owen
1993 ◽  
Vol 04 (04) ◽  
pp. 883-907 ◽  
Author(s):  
C. BONA ◽  
J. MASSÓ

The construction of numerical solutions of Einstein's General Relativity equations is formulated as an initial-value problem. The space-plus-time (3 + 1) decomposition of the spacetime metric tensor is used to discuss the structure of the field equations. The resulting evolution system is shown to depend in a crucial way on the coordinate gauge. The mandatory use of singularity avoiding coordinate conditions (like maximal slicing or similar gauges) is explained. A brief historical review of Numerical Relativity is included, showing the enormous effort in constructing codes based in these gauges, which lead to non-hyperbolic evolution systems, using "ad hoc" numerical techniques. A new family of first order hyperbolic evolution systems for the vacuum Einstein field equations in the harmonic slicing gauge is presented. This family depends on a symmetric 3 × 3 array of parameters which can be used to scale the dynamical variables in future numerical applications.


Author(s):  
Anatoly Kochubei

AbstractWe describe a class of evolution systems of linear partial differential equations with the Caputo-Dzhrbashyan fractional derivative of order α ∈ (0, 1) in the time variable t and the first order derivatives in spatial variables x = (x 1, …, x n), which can be considered as a fractional analogue of the class of hyperbolic systems. For such systems, we construct a fundamental solution of the Cauchy problem having exponential decay outside the fractional light cone {(t,x) : |t -α| ≤ 1}.


2002 ◽  
Vol 66 (8) ◽  
Author(s):  
C. Bona ◽  
T. Ledvinka ◽  
C. Palenzuela

2016 ◽  
Vol 16 (3) ◽  
pp. 409-428 ◽  
Author(s):  
Willy Dörfler ◽  
Stefan Findeisen ◽  
Christian Wieners

AbstractWe introduce a space-time discretization for linear first-order hyperbolic evolution systems using a discontinuous Galerkin approximation in space and a Petrov–Galerkin scheme in time. We show well-posedness and convergence of the discrete system. Then we introduce an adaptive strategy based on goal-oriented dual-weighted error estimation. The full space-time linear system is solved with a parallel multilevel preconditioner. Numerical experiments for the linear transport equation and the Maxwell equation in 2D underline the efficiency of the overall adaptive solution process.


2010 ◽  
Vol 82 (12) ◽  
Author(s):  
C. Bona ◽  
C. Bona-Casas ◽  
C. Palenzuela

1997 ◽  
Vol 56 (6) ◽  
pp. 3405-3415 ◽  
Author(s):  
C. Bona ◽  
J. Massó ◽  
E. Seidel ◽  
J. Stela

Sign in / Sign up

Export Citation Format

Share Document