scholarly journals A MODEL PROBLEM FOR THE INITIAL-BOUNDARY VALUE FORMULATION OF EINSTEIN'S FIELD EQUATIONS

2005 ◽  
Vol 02 (02) ◽  
pp. 397-435 ◽  
Author(s):  
OSCAR REULA ◽  
OLIVIER SARBACH

In many numerical implementations of the Cauchy formulation of Einstein's field equations one encounters artificial boundaries which raises the issue of specifying boundary conditions. Such conditions have to be chosen carefully. In particular, they should be compatible with the constraints, yield a well posed initial-boundary value formulation and incorporate some physically desirable properties like, for instance, minimizing reflections of gravitational radiation. Motivated by the problem in General Relativity, we analyze a model problem, consisting of a formulation of Maxwell's equations on a spatially compact region of space–time with timelike boundaries. The form in which the equations are written is such that their structure is very similar to the Einstein–Christoffel symmetric hyperbolic formulations of Einstein's field equations. For this model problem, we specify a family of Sommerfeld-type constraint-preserving boundary conditions and show that the resulting initial-boundary value formulations are well posed. We expect that these results can be generalized to the Einstein–Christoffel formulations of General Relativity, at least in the case of linearizations about a stationary background.

2011 ◽  
Vol 152 (3) ◽  
pp. 473-496 ◽  
Author(s):  
DAVID A. SMITH

AbstractWe study initial-boundary value problems for linear evolution equations of arbitrary spatial order, subject to arbitrary linear boundary conditions and posed on a rectangular 1-space, 1-time domain. We give a new characterisation of the boundary conditions that specify well-posed problems using Fokas' transform method. We also give a sufficient condition guaranteeing that the solution can be represented using a series.The relevant condition, the analyticity at infinity of certain meromorphic functions within particular sectors, is significantly more concrete and easier to test than the previous criterion, based on the existence of admissible functions.


2007 ◽  
Vol 04 (04) ◽  
pp. 587-612 ◽  
Author(s):  
ALEXANDER M. ALEKSEENKO

A well-posed initial-boundary value problem is formulated for the model problem of the vector wave equation subject to the divergence-free constraint. Existence, uniqueness and stability of the solution is proved by reduction to a system evolving the constraint quantity statically, i.e. the second time derivative of the constraint quantity is zero. A new set of radiation-controlling constraint-preserving boundary conditions is constructed for the free evolution problem. Comparison between the new conditions and the standard constraint-preserving boundary conditions is made using the Fourier–Laplace analysis and the power series decomposition in time. The new boundary conditions satisfy the Kreiss condition and are free from the ill-posed modes growing polynomially in time.


2005 ◽  
Vol 02 (04) ◽  
pp. 839-883 ◽  
Author(s):  
OLIVIER SARBACH ◽  
MANUEL TIGLIO

Outer boundary conditions for strongly and symmetric hyperbolic formulations of 3D Einstein's field equations with a live gauge condition are discussed. The boundary conditions have the property that they ensure constraint propagation and control in a sense made precise in this article the physical degrees of freedom at the boundary. We use Fourier–Laplace transformation techniques to find necessary conditions for the well posedness of the resulting initial-boundary value problem and integrate the resulting three-dimensional nonlinear equations using a finite-differencing code. We obtain a set of constraint-preserving boundary conditions which pass a robust numerical stability test. We explicitly compare these new boundary conditions to standard, maximally dissipative ones through Brill wave evolutions. Our numerical results explicitly show that in the latter case the constraint variables, describing the violation of the constraints, do not converge to zero when resolution is increased while for the new boundary conditions, the constraint variables do decrease as resolution is increased. As an application, we inject pulses of "gravitational radiation" through the boundaries of an initially flat spacetime domain, with enough amplitude to generate strong fields and induce large curvature scalars, showing that our boundary conditions are robust enough to handle nonlinear dynamics. We expect our boundary conditions to be useful for improving the accuracy and stability of current binary black hole and binary neutron star simulations, for a successful implementation of characteristic or perturbative matching techniques, and other applications. We also discuss limitations of our approach and possible future directions.


Author(s):  
А.М. Слиденко ◽  
В.М. Слиденко

Приводится анализ механических колебаний элементов ударного устройства с помощью модели стержневого типа. Ударник и инструмент связаны упругими и диссипативными элементами, которые имитируют их взаимодействие. Аналогично моделируется взаимодействие инструмента с рабочей средой. Сформулирована начально-краевая задача для системы двух волновых уравнений с учетом переменных поперечных сечений стержней. Площади поперечных сечений определяются параметрическими формулами при сохранении объемов стержней. Параметрические формулы позволяют получать различного вида зависимости площади поперечного сечения стержня от его длины. Начальные условия отражают физическую картину взаимодействия инструмента с ударником и рабочей средой. Краевые условия описывают контактные взаимодействия ударника с инструментом и последнего с рабочей средой. В качестве модельной задачи рассматривается соударение ударника и инструмента через элемент большой жесткости. Начально-краевая задача исследуется разностным методом. Проводится сравнение решений задачи, полученных с помощью двухслойной и трехслойной разностных схем. Такие схемы реализованы в общей компьютерной программе в системе Mathcad. Показано, что при вычислениях распределения нормальных напряжений по длине стержня лучшими свойствами относительно устойчивости обладает двухслойная схема The article gives the analysis of mechanical vibrations of the impact device elements using the model of the rod type. The hammer and the tool are connected by elastic and dissipative elements that simulate their interaction. The interaction of the tool with the processing medium is simulated in a similar way. An initial boundary-value problem is formulated for a system of two wave equations taking into account the variable cross sections of the rods. Cross-sectional areas are determined by parametric formulas maintaining the volume of the rods. Parametric formulas allow one to obtain various dependence types of the cross-sectional area of the rod on its length. The initial and boundary conditions reflect the physical phenomenon of the tool interaction with the processing medium, and also describe the contact interactions of the hammer with the tool. The impacting of the hammer and the tool through an element of high rigidity is considered as a model problem. To control the limiting values, the solution of the model problem by the Fourier method is used. The initial-boundary-value problem is investigated by the difference method. A comparison of solutions obtained for the two-layer and three-layer difference schemes is given. Such schemes are realized in a common computer program in the Mathcad. It is shown that the two-layer scheme has the best properties in relation to stability while calculating the distribution of normal voltage along the length of the rod


2020 ◽  
Vol 75 (8) ◽  
pp. 713-725 ◽  
Author(s):  
Guenbo Hwang

AbstractInitial-boundary value problems for the one-dimensional linear advection–dispersion equation with decay (LAD) are studied by utilizing a unified method, known as the Fokas method. The method takes advantage of the spectral analysis of both parts of Lax pair and the global algebraic relation coupling all initial and boundary values. We present the explicit analytical solution of the LAD equation posed on the half line and a finite interval with general initial and boundary conditions. In addition, for the case of periodic boundary conditions, we show that the solution of the LAD equation is asymptotically t-periodic for large t if the Dirichlet boundary datum is periodic in t. Furthermore, it can be shown that if the Dirichlet boundary value is asymptotically periodic for large t, then so is the unknown Neumann boundary value, which is uniquely characterized in terms of the given asymptotically periodic Dirichlet boundary datum. The analytical predictions for large t are compared with numerical results showing the excellent agreement.


Sign in / Sign up

Export Citation Format

Share Document