Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Operator-Based Local Discontinuous Galerkin Method Compatible With the BSSN Formulation of the Einstein Equations
58
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Discontinuous Galerkin Finite Element (DGFE) methods offer a mathematically beautiful, computationally efficient, and efficiently parallelizable way to solve hyperbolic partial differential equations. These properties make them highly desirable for numerical calculations in relativistic astrophysics and many other fields. The BSSN formulation of the Einstein equations has repeatedly demonstrated its robustness. The formulation is not only stable but allows for puncture-type evolutions of black hole systems. To-date no one has been able to solve the full (3+1)-dimensional BSSN equations using DGFE methods. This is partly because DGFE discretization often occurs at the level of the equations, not the derivative operator, and partly because DGFE methods are traditionally formulated for manifestly flux-conservative systems. By discretizing the derivative operator, we generalize a particular flavor of DGFE methods, Local DG methods, to solve arbitrary second-order hyperbolic equations. Because we discretize at the level of the derivative operator, our method can be interpreted as either a DGFE method or as a finite differences stencil with non-constant coefficients.
rate research
Read More
We present a strongly hyperbolic first-order formulation of the Einstein equations based on the conformal and covariant Z4 system (CCZ4) with constraint-violation damping, which we refer to as FO-CCZ4. As CCZ4, this formulation combines the advantages of a conformal and traceless formulation, with the suppression of constraint violations given by the damping terms, but being first order in time and space, it is particularly suited for a discontinuous Galerkin (DG) implementation. The strongly hyperbolic first-order formulation has been obtained by making careful use of first and second-order ordering constraints. A proof of strong hyperbolicity is given for a selected choice of standard gauges via an analytical computation of the entire eigenstructure of the FO-CCZ4 system. The resulting governing partial differential equations system is written in non-conservative form and requires the evolution of 58 unknowns. A key feature of our formulation is that the first-order CCZ4 system decouples into a set of pure ordinary differential equations and a reduced hyperbolic system of partial differential equations that contains only linearly degenerate fields. We implement FO-CCZ4 in a high-order path-conservative arbitrary-high-order-method-using-derivatives (ADER)-DG scheme with adaptive mesh refinement and local time-stepping, supplemented with a third-order ADER-WENO subcell finite-volume limiter in order to deal with singularities arising with black holes. We validate the correctness of the formulation through a series of standard tests in vacuum, performed in one, two and three spatial dimensions, and also present preliminary results on the evolution of binary black-hole systems. To the best of our knowledge, these are the first successful three-dimensional simulations of moving punctures carried out with high-order DG schemes using a first-order formulation of the Einstein equations.
In this paper, we develop an oscillation free local discontinuous Galerkin (OFLDG) method for solving nonlinear degenerate parabolic equations. Following the idea of our recent work [J. Lu, Y. Liu, and C.-W. Shu, SIAM J. Numer. Anal. 59(2021), pp. 1299-1324.], we add the damping terms to the LDG scheme to control the spurious oscillations when solutions have a large gradient. The $L^2$-stability and optimal priori error estimates for the semi-discrete scheme are established. The numerical experiments demonstrate that the proposed method maintains the high-order accuracy and controls the spurious oscillations well.
In this paper, an energy-based discontinuous Galerkin method for dynamic Euler-Bernoulli beam equations is developed. The resulting method is energy-dissipating or energy-conserving depending on the simple, mesh-independent choice of numerical fluxes. By introducing a velocity field, the original problem is transformed into a first-order in time system. In our formulation, the discontinuous Galerkin approximations for the original displacement field and the auxiliary velocity field are not restricted to be in the same space. In particular, a given accuracy can be achieved with the fewest degrees of freedom when the degree for the approximation space of the velocity field is two orders lower than the degree of approximation space for the displacement field. In addition, we establish the error estimates in an energy norm and demonstrate the corresponding optimal convergence in numerical experiments.
We adopt a reference-metric approach to generalize a covariant and conformal version of the Z4 system of the Einstein equations. We refer to the resulting system as ``fully covariant and conformal, or fCCZ4 for short, since it is well suited for curvilinear as well as Cartesian coordinates. We implement this fCCZ4 formalism in spherical polar coordinates under the assumption of spherical symmetry using a partially-implicit Runge-Kutta (PIRK) method and show that our code can evolve both vacuum and non-vacuum spacetimes without encountering instabilities. Our method does not require regularization of the equations to handle coordinate singularities, nor does it depend on constraint-preserving outer boundary conditions, nor does it need any modifications of the equations for evolutions of black holes. We perform several tests and compare the performance of the fCCZ4 system, for different choices of certain free parameters, with that of BSSN. Confirming earlier results we find that, for an optimal choice of these parameters, and for neutron-star spacetimes, the violations of the Hamiltonian constraint can be between 1 and 3 orders of magnitude smaller in the fCCZ4 system than in the BSSN formulation. For black-hole spacetimes, on the other hand, any advantages of fCCZ4 over BSSN are less evident.
Gravitational wave emission from extreme mass ratio binaries (EMRBs) should be detectable by the joint NASA-ESA LISA project, spurring interest in analytical and numerical methods for investigating EMRBs. We describe a discontinuous Galerkin (dG) method for solving the distributionally forced 1+1 wave equations which arise when modeling EMRBs via the perturbation theory of Schwarzschild blackholes. Despite the presence of jump discontinuities in the relevant polar and axial gravitational master functions, our dG method achieves global spectral accuracy, provided that we know the instantaneous position, velocity, and acceleration of the small particle. Here these variables are known, since we assume that the particle follows a timelike geodesic of the Schwarzschild geometry. We document the results of several numerical experiments testing our method, and in our concluding section discuss the possible inclusion of gravitational self-force effects.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
