scholarly journals Applied Mathematics of Space-time & Space+time: Problems in General Relativity and Cosmology

2021 ◽  
Author(s):  
◽  
Celine Cattoen

<p>Cosmography is the part of cosmology that proceeds by making minimal dynamic assumptions. That is, one does not assume the Friedmann equations (Einstein equations) unless and until absolutely necessary. On the other hand, cosmodynamics is the part of cosmology that relates the geometry to the density and pressure using the Friedmann equations. In both frameworks, we consider the amount of information and the nature of the constraints we can obtain from the Hubble flow in a FLRW universe. Indeed, the cosmological parameters contained in the Hubble relation between distance and redshift provide information on the behaviour of the universe (expansion, acceleration etc...). In the first framework, it is possible to concentrate more directly on the observational situation in a model-independent manner. We perform a number of inter-related cosmographic fits to supernova datasets, and pay particular attention to the extent to which the choice of distance scale and manner of representing the redshift scale affect the cosmological parameters. In the second framework, we use the class of w-parameter models which has become increasingly popular in the last decade. We explore the extent to which a constraint on the w-parameter leads to useful and non-trivial constraints on the Hubble flow in terms of cosmological parameters H(z), density p(z), density parameter O(z), distance scales d(z), and lookback time T(z). On another front, Numerical Relativity has experienced many breakthroughs since 2005, with full inspiral-merger-ringdown simulations now possible. One of the main goals is to provide very accurate templates of gravitational waves for ground-based and space-based interferometers. We explore the potential of a very recent and accurate numerical method, the Spectral Element Method (SEM), for Numerical Relativity, by treating a singular Schwarszchild black hole evolution as a test case. Spectral elements combine the theory of spectral and pseudo-spectral methods for high order polynomials and the variational formulation of finite elements and the associated geometric flexibility. We use the BSSN formulation of the Einstein equations with the method of the moving punctures. After applying the variational formulation to the BSSN system, we present several possible weak forms of this system and its spectral element discretization in space. We use a Runge-Kutta fourth order time discretization. The accuracy of high order methods can deteriorate in the presence of discontinuities or sharp gradients. We show that we can treat the element that contains the puncture with a filtering method to avoid artificial and spurious oscillations. These might form and propagate into the domain coming from discontinuous initial data from the BSSN system.</p>

2021 ◽  
Author(s):  
◽  
Celine Cattoen

<p>Cosmography is the part of cosmology that proceeds by making minimal dynamic assumptions. That is, one does not assume the Friedmann equations (Einstein equations) unless and until absolutely necessary. On the other hand, cosmodynamics is the part of cosmology that relates the geometry to the density and pressure using the Friedmann equations. In both frameworks, we consider the amount of information and the nature of the constraints we can obtain from the Hubble flow in a FLRW universe. Indeed, the cosmological parameters contained in the Hubble relation between distance and redshift provide information on the behaviour of the universe (expansion, acceleration etc...). In the first framework, it is possible to concentrate more directly on the observational situation in a model-independent manner. We perform a number of inter-related cosmographic fits to supernova datasets, and pay particular attention to the extent to which the choice of distance scale and manner of representing the redshift scale affect the cosmological parameters. In the second framework, we use the class of w-parameter models which has become increasingly popular in the last decade. We explore the extent to which a constraint on the w-parameter leads to useful and non-trivial constraints on the Hubble flow in terms of cosmological parameters H(z), density p(z), density parameter O(z), distance scales d(z), and lookback time T(z). On another front, Numerical Relativity has experienced many breakthroughs since 2005, with full inspiral-merger-ringdown simulations now possible. One of the main goals is to provide very accurate templates of gravitational waves for ground-based and space-based interferometers. We explore the potential of a very recent and accurate numerical method, the Spectral Element Method (SEM), for Numerical Relativity, by treating a singular Schwarszchild black hole evolution as a test case. Spectral elements combine the theory of spectral and pseudo-spectral methods for high order polynomials and the variational formulation of finite elements and the associated geometric flexibility. We use the BSSN formulation of the Einstein equations with the method of the moving punctures. After applying the variational formulation to the BSSN system, we present several possible weak forms of this system and its spectral element discretization in space. We use a Runge-Kutta fourth order time discretization. The accuracy of high order methods can deteriorate in the presence of discontinuities or sharp gradients. We show that we can treat the element that contains the puncture with a filtering method to avoid artificial and spurious oscillations. These might form and propagate into the domain coming from discontinuous initial data from the BSSN system.</p>


Author(s):  
Lonny L. Thompson

Abstract A Computational Structural Acoustics (CSA) capability for solving scattering, radiation, and other problems related to the acoustics of submerged structures has been developed by employing some of the recent algorithmic trends in Computational Fluid Dynamics (CFD), namely time-discontinuous Galerkin Least-Squares finite element methods. Traditional computational methods toward simulation of acoustic radiation and scattering from submerged elastic bodies have been primarily based on frequency domain formulations. These classical time-harmonic approaches (including boundary element, finite element, and finite difference methods) have been successful for problems involving a limited range of frequencies (narrow band response) and scales (wavelengths) that are large compared to the characteristic dimensions of the elastic structure. Attempts at solving large-scale structural acoustic systems with dimensions that are much larger than the operating wavelengths and which are complex, consisting of many different components with different scales and broadband frequencies, has revealed limitations of many of the classical methods. As a result, there has been renewed interest in new innovative approaches, including time-domain approaches. This paper describes recent advances in the development of a new class of high-order accurate and unconditionally stable space-time methods for structural acoustics which employ finite element discretization of the time domain as well as the usual discretization of the spatial domain. The formulation is based on a space-time variational equation for both the acoustic fluid and elastic structure together with their interaction. Topics to be discussed include the development and implementation of higher-order accurate non-reflecting boundary conditions based on the exact impedance relation through the. Dirichlet-to-Neumann (DtN) map, and a multi-field representation for the acoustic fluid based on independent pressure and velocity potential variables. Numerical examples involving radiation and scattering of acoustic waves are presented to illustrate the high-order accuracy achieved by the new methodology for CSA.


2018 ◽  
Vol 16 (01) ◽  
pp. 1850093 ◽  
Author(s):  
Chaoxu Pei ◽  
Mark Sussman ◽  
M. Yousuff Hussaini

A space-time discontinuous Galerkin spectral element method is combined with two different approaches for treating problems with discontinuous solutions: (i) adding a space-time dependent artificial viscosity, and (ii) tracking the discontinuity with space-time spectral accuracy. A Picard iteration method is employed to solve nonlinear system of equations derived from the space-time DG spectral element discretization. Spectral accuracy in both space and time is demonstrated for the Burgers’ equation with a smooth solution. For tests with discontinuities, the present space-time method enables better accuracy at capturing the shock strength in the element containing shock when higher order polynomials in both space and time are used. The spectral accuracy of the shock speed and location is demonstrated for the solution of the inviscid Burgers’ equation obtained by the tracking method.


2011 ◽  
Vol 3 (4) ◽  
pp. 448-469 ◽  
Author(s):  
Zakaria Belhachmi ◽  
Andreas Karageorghis

AbstractIn this paper, we study the numerical solution of the Stokes system in deformed axisymmetric geometries. In the azimuthal direction the discretization is carried out by using truncated Fourier series, thus reducing the dimension of the problem. The resulting two-dimensional problems are discretized using the spectral element method which is based on the variational formulation in primitive variables. The meridian domain is subdivided into elements, in each of which the solution is approximated by truncated polynomial series. The results of numerical experiments for several geometries are presented.


2010 ◽  
Vol 19 (03) ◽  
pp. 233-244
Author(s):  
RIZWAN UL HAQ ANSARI ◽  
P. K. SURESH

We considered an anisotropic brane embedded in a five-dimensional bulk space–time. The Einstein equations and Friedmann equations are obtained. The corresponding Friedmann equations are modified with quadratic corrections to energy density. The anisotropic DGP model also admits a self-accelerating phase in the late universe. The acceleration conditions are obtained and they are consistent with the braneworld scenario. The acceleration condition for DGP model, in the late universe, is dependent on the crossover scale rc.


Author(s):  
Diogo B. Pitz ◽  
John W. Chew ◽  
Olaf Marxen ◽  
Nicholas J. Hills

A high-order numerical method is employed to investigate flow in a rotor/stator cavity without heat transfer and buoyant flow in a rotor/rotor cavity. The numerical tool used employs a spectral element discretization in two dimensions and a Fourier expansion in the remaining direction, which is periodic and corresponds to the azimuthal coordinate in cylindrical coordinates. The spectral element approximation uses a Galerkin method to discretize the governing equations, but employs high-order polynomials within each element to obtain spectral accuracy. A second-order, semi-implicit, stiffly stable algorithm is used for the time discretization. Numerical results obtained for the rotor/stator cavity compare favorably with experimental results for Reynolds numbers up to Re1 = 106 in terms of velocities and Reynolds stresses. The buoyancy-driven flow is simulated using the Boussinesq approximation. Predictions are compared with previous computational and experimental results. Analysis of the present results shows close correspondence to natural convection in a gravitational field and consistency with experimentally observed flow structures in a water-filled rotating annulus. Predicted mean heat transfer levels are higher than the available measurements for an air-filled rotating annulus, but in agreement with correlations for natural convection under gravity.


2006 ◽  
Vol 03 (01) ◽  
pp. 81-141 ◽  
Author(s):  
PIOTR T. CHRUŚCIEL ◽  
SZYMON ŁȨSKI

The study of Einstein equations leads naturally to Cauchy problems with initial data on hypersurfaces which closely resemble hyperboloids in Minkowski space-time, and with initial data with polyhomogeneous asymptotics, that is, with asymptotic expansions in terms of powers of ln r and inverse powers of r. Such expansions also arise in the conformal method for analysing wave equations in odd space-time dimension. In recent work it has been shown that for non-linear wave equations, or for wave maps, polyhomogeneous initial data lead to solutions which are also polyhomogeneous provided that an infinite hierarchy of corner conditions holds. In this paper we show that the result is true regardless of corner conditions.


Sign in / Sign up

Export Citation Format

Share Document