ترغب بنشر مسار تعليمي؟ اضغط هنا

Constructing of constraint preserving scheme for Einstein equations

46   0   0.0 ( 0 )
 نشر من قبل Takuya Tsuchiya
 تاريخ النشر 2016
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We propose a new numerical scheme of evolution for the Einstein equations using the discrete variational derivative method (DVDM). We derive the discrete evolution equation of the constraint using this scheme and show the constraint preserves in the discrete level. In addition, to confirm the numerical stability using this scheme, we perform some numerical simulations by discretized equations with the Crank-Nicolson scheme and with the new scheme, and we find that the new discretized equations have better stability than that of the Crank-Nicolson scheme.



قيم البحث

اقرأ أيضاً

We derive the discretized Maxwells equations using the discrete variational derivative method (DVDM), calculate the evolution equation of the constraint, and confirm that the equation is satisfied at the discrete level. Numerical simulations showed t hat the results obtained by the DVDM are superior to those obtained by the Crank-Nicolson scheme. In addition, we study the two types of the discretized Maxwells equations by the DVDM and conclude that if the evolution equation of the constraint is not conserved at the discrete level, then the numerical results are also unstable.
Exact black hole solutions in the Einstein-Maxwell-scalar theory are constructed. They are the extensions of dilaton black holes in de Sitter or anti de Sitter universe. As a result, except for a scalar potential, a coupling function between the scal ar field and the Maxwell invariant is present. Then the corresponding Smarr formula and the first law of thermodynamics are investigated.
The linear Einstein-Boltzmann equations describe the evolution of perturbations in the universe and its numerical solutions play a central role in cosmology. We revisit this system of differential equations and present a detailed investigation of its mathematical properties. For this purpose, we focus on a simplified set of equations aimed at describing the broad features of the matter power spectrum. We first perform an eigenvalue analysis and study the onset of oscillations in the system signaled by the transition from real to complex eigenvalues. We then provide a stability criterion of different numerical schemes for this linear system and estimate the associated step-size. We elucidate the stiffness property of the Einstein-Boltzmann system and show how it can be characterized in terms of the eigenvalues. While the parameters of the system are time dependent making it non-autonomous, we define an adiabatic regime where the parameters vary slowly enough for the system to be quasi-autonomous. We summarize the different regimes of the system for these different criteria as function of wave number $k$ and scale factor $a$. We also provide a compendium of analytic solutions for all perturbation variables in 6 limits on the $k$-$a$ plane and express them explicitly in terms of initial conditions. These results are aimed to help the further development and testing of numerical cosmological Boltzmann solvers.
In this work a static solution of Einstein-Cartan (EC) equations in 2+1 dimensional space-time is given by considering classical spin-1/2 field as external source for torsion of the space-time. Here, the torsion tensor is obtained from metricity cond ition for the connection and the static spinor field is determined as the solution of Dirac equation in 2+1 spacetime with non-zero cosmological constant and torsion. The torsion itself is considered as a non-dynamical field.
Theories of gravity other than general relativity (GR) can explain the observed cosmic acceleration without a cosmological constant. One such class of theories of gravity is f(R). Metric f(R) theories have been proven to be equivalent to Brans-Dicke (BD) scalar-tensor gravity without a kinetic term. Using this equivalence and a 3+1 decomposition of the theory it has been shown that metric f(R) gravity admits a well-posed initial value problem. However, it has not been proven that the 3+1 evolution equations of metric f(R) gravity preserve the (hamiltonian and momentum) constraints. In this paper we show that this is indeed the case. In addition, we show that the mathematical form of the constraint propagation equations in BD-equilavent f(R) gravity and in f(R) gravity in both the Jordan and Einstein frames, is exactly the same as in the standard ADM 3+1 decomposition of GR. Finally, we point out that current numerical relativity codes can incorporate the 3+1 evolution equations of metric f(R) gravity by modifying the stress-energy tensor and adding an additional scalar field evolution equation. We hope that this work will serve as a starting point for relativists to develop fully dynamical codes for valid f(R) models.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا