Global smooth solutions to relativistic Euler-Poisson equations with repulsive force

In this paper, we are concerned with the global existence of smooth solutions for the one dimen- sional relativistic Euler-Poisson equations： Combining certain physical background, the relativistic Euler-Poisson model is derived mathematically. By using an invariant of Lax＇s method, we will give a sufficient condition for the existence of a global smooth solution to the one-dimensional Euler-Poisson equations with repulsive force.     

一维相对论振子运动的研究

该文运用椭圆积分的理论，给出了一维相对论振子运动的解析解以及振子的振动周期. 指出相对论振子的振动周期不但与振子的固有性质有关，而且还与振子的振幅有关.     

A relaxation method of an alternating iterative MFS algorithm for the Cauchy problem associated with the two‐dimensional modified Helmholtz equation

We propose two algorithms involving the relaxation of either the given Dirichlet data or the prescribed Neumann data on the over‐specified boundary in the case of the alternating iterative algorithm of Kozlov et al. (USSR Comput Math Math Phys 31 (1991), 45–52) applied to the Cauchy problem for the two‐dimensional modified Helmholtz equation. The two mixed, well‐posed and direct problems corresponding to every iteration of the numerical procedure are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is selected according to the generalized cross‐validation criterion. The iterative MFS algorithms with relaxation are tested for Cauchy problems associated with the modified Helmholtz equation in two‐dimensional geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the method. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Simplified immersed interface methods for elliptic interface problems with straight interfaces

In this article, we propose simplified immersed interface methods for elliptic partial/ordinary differential equations with discontinuous coefficients across interfaces that are few isolated points in 1D, and straight lines in 2D. For one‐dimensional problems or two‐dimensional problems with circular interfaces, we propose a conservative second‐order finite difference scheme whose coefficient matrix is symmetric and definite. For two‐dimensional problems with straight interfaces, we first propose a conservative first‐order finite difference scheme, then use the Richardson extrapolation technique to get a second‐order method. In both cases, the finite difference coefficients are almost the same as those for regular problems. Error analysis is given along with numerical example. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 188–203, 2012     

On well‐posedness of nonclassical problems for elliptic equations

In the present paper, we consider nonclassical problems for multidimensional elliptic equations. A finite difference method for solving these nonlocal boundary value problems is presented. Stability, almost coercive stability and coercive stability for the solutions of first and second orders of approximation are obtained. The theoretical statements for the solutions of these difference schemes are supported by numerical examples for the two‐dimensional elliptic equations. Copyright © 2014 John Wiley & Sons, Ltd.     

Asymptotic stability of the relativistic Maxwellian

Global classical solutions near the relativistic Maxwellian are constructed for the relativistic Boltzmann equation in both a periodic box and the whole space. For both cases, we are able to get the non‐negativity of the global solution under less restriction than in (Publ. Res. Inst. Math. Sci. 1993; 29 :301–347) on the scattering kernel. Copyright © 2006 John Wiley & Sons, Ltd.     

Steger–Warming flux vector splitting method for special relativistic hydrodynamics

This paper discusses the properties of the rotational invariance and hyperbolicity in time of the governing equations of the ideal special relativistic hydrodynamics and proves for the first time that the ideal relativistic hydrodynamical equations satisfy the homogeneity property, which is the footstone of the Steger–Warming flux vector splitting method [J. L. Steger and R. F. Warming, J. Comput. Phys., 40(1981), 263–293]. On the basis of this remarkable property, the Steger–Warming flux vector splitting (SW‐FVS) is given. Two high‐resolution SW‐FVS schemes are also given on the basis of the initial reconstructions of the solutions and the fluxes, respectively. Several numerical experiments are conducted to validate the performance of the SW‐FVS method. Copyright © 2013 John Wiley & Sons, Ltd.     

Preconditioning trace coupled 3d‐1d systems using fractional Laplacian

Multiscale or multiphysics problems often involve coupling of partial differential equations posed on domains of different dimensionality. In this work, we consider a simplified model problem of a 3d‐1d coupling and the main objective is to construct algorithms that may utilize standard multilevel algorithms for the 3d domain, which has the dominating computational complexity. Preconditioning for a system of two elliptic problems posed, respectively, in a three‐dimensional domain and an embedded one dimensional curve and coupled by the trace constraint is discussed. Investigating numerically the properties of the well‐defined discrete trace operator, it is found that negative fractional Sobolev norms are suitable preconditioners for the Schur complement of the system. The norms are employed to construct a robust block diagonal preconditioner for the coupled problem.     

Quasi-exact solution of the problem of relativistic bound states in the (1+1)-dimensional case

We investigate the problem of bound states for bosons and fermions in the framework of the relativistic configurational representation with the kinetic part of the Hamiltonian containing purely imaginary finite shift operators e^±ihd/dx instead of differential operators. For local (quasi)potentials of the type of a rectangular potential well in the (1+1)-dimensional case, we elaborate effective methods for solving the problem analytically that allow finding the spectrum and investigating the properties of wave functions in a wide parameter range. We show that the properties of these relativistic bound states differ essentially from those of the corresponding solutions of the Schrödinger and Dirac equations in a static external potential of the same form in a number of fundamental aspects both at the level of wave functions and of the energy spectrum structure. In particular, competition between ? and the potential parameters arises, as a result of which these distinctions are retained at low-lying levels in a sufficiently deep potential well for ? ? 1 and the boson and fermion energy spectra become identical.     

Gravity and Lorentz Force

The relativistic covariant equations are obtained for the relativistic Newton gravity law proposed by Poincaré. These equations are solved for the case where the mass of one gravitating body is equal to zero.     

Delta Shocks and Vacuum States in Vanishing Pressure Limits of Solutions to the Relativistic Euler Equations

The Riemann problems for the Euler system of conservation laws of energy and momentum in special relativity as pressure vanishes are considered. The Riemann solutions for the pressureless relativistic Euler equations are obtained constructively. There are two kinds of solutions, the one involves delta shock wave and the other involves vacuum. The authors prove that these two kinds of solutions are the limits of the solutions as pressure vanishes in the Euler system of conservation laws of energy and momentum in special relativity.     

The Robin problem for the scalar Oseen equation

We study the Robin problem for the scalar Oseen equation in an open n‐dimensional set with compact Ljapunov boundary. We prescribe two types of Robin boundary conditions, and prove the unique solvability of these problems as well as a representation formula for the solution in form of a scalar Oseen single layer potential. Moreover, we prove the maximum principle for the solution to the Robin problem of the scalar Oseen equation. Copyright © 2013 John Wiley & Sons, Ltd.     

Uniform approximation of singularly perturbed reaction‐diffusion problems by the finite element method on a Shishkin mesh

We consider the numerical approximation of singularly perturbed reaction‐diffusion problems over two‐dimensional domains with smooth boundary. Using the h version of the finite element method over appropriately designed piecewise uniform (Shishkin) meshes, we are able to uniformly approximate the solution at a quasi‐optimal rate. The results of numerical computations showing agreement with the analysis are also presented. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 89–111, 2003     

Inf‐sup stability of Petrov‐Galerkin immersed finite element methods for one‐dimensional elliptic interface problems

In this article, we analyze the Petrov‐Galerkin immersed finite element method (PG‐IFEM) when applied to one‐dimensional elliptic interface problems. In the PG‐IFEM (T. Hou, X. Wu and Y. Zhang, Commun. Math. Sci., 2 (2004), 185‐205, and S. Hou and X. Liu, J. Comput. Phys., 202 (2005), 411‐445), the classic immersed finite element (IFE) space was taken as the trial space while the conforming linear finite element space was taken as the test space. We first prove the inf‐sup condition of the PG‐IFEM and then show the optimal error estimate in the energy norm. We also show the optimal estimate of the condition number of the stiffness matrix. The results are extended to two dimensional problems in a special case.     

Global existence of smooth solution to 2D relativistic membrane equation in asymptotically flat FLRW space‐time

This paper is concerned with the 2D relativistic membrane equation evolving in curved space‐time, whose metric is prescribed as a small perturbation to the flat Friedmann‐Lemaître‐Robertson‐Walker (FLRW) metric with time‐dependent coefficients. Thanks to the partial “null structure” of the membrane equation and the properties of the background metric, we could prove the global stability of a class of time‐dependent solutions by weighted energy estimate.     

On global symmetric solutions to the relativistic Vlasov–Poisson equation in three space dimensions

A collisionless plasma is modelled by the Vlasov–Maxwell system. In the presence of very large velocities, relativistic corrections are meaningful. When magnetic effects are ignored this formally becomes the relativistic Vlasov–Poisson equation. The initial datum for the phase space density ƒ₀(x, v) is assumed to be sufficiently smooth, non‐negative and cylindrically symmetric. If the (two‐dimensional) angular momentum is bounded away from zero on the support of ƒ₀(x, v), it is shown that a smooth solution to the Cauchy problem exists for all times. Copyright © 2001 John Wiley & Sons, Ltd.     

Vanishing pressure limits of Riemann solutions to the isentropic relativistic Euler system for Chaplygin gas

The Riemann solutions to the isentropic relativistic Euler system for Chaplygin gas with a small parameter are considered. Unlike the polytropic or barotropic gas cases, we find that firstly, as the parameter decreases to a certain critical number, the two-shock solution converges to a delta shock wave solution of the same system. Moreover, as the parameter goes to zero, that is, the pressure vanishes, the solution is nothing but the delta shock wave solution to the zero-pressure relativistic Euler system. Meanwhile, the two-rarefaction wave solution tends to the vacuum solution to the zero-pressure relativistic system, and the solution containing one rarefaction wave and one shock wave tends to the contact discontinuity solution to the zero-pressure relativistic system as pressure vanishes.     

Derivation of the Hydrodynamic Equations in the Framework of the Bogoliubov Functional Hypothesis

We generalize the Bogoliubov functional hypothesis to the case of multiparticle interaction depending on both the coordinates and momenta of particles. We illustrate this with the examples of two weakly relativistic models: the Darwin model in the theory of charged particles and the Fock model in the general theory of relativity. For these models, based on the chain of the BBGKY equations, we calculate weakly relativistic corrections to the classical transport coefficients and find the conditions under which there is no bijective relation between the parameters of the local equilibrium distribution and the hydrodynamic variables.     

Discontinuous Galerkin method for solving the black‐oil problem in porous media

We introduce a new algorithm for solving the three‐component three‐phase flow problem in two‐dimensional and three‐dimensional heterogeneous media. The oil and gas components can be found in the liquid and vapor phases, whereas the aqueous phase is only composed of water component. The numerical scheme employs a sequential implicit formulation discretized with discontinuous finite elements. Capillarity and gravity effects are included. The method is shown to be accurate and robust for several test problems. It has been carefully designed so that calculation of appearance and disappearance of phases does not require additional steps.     

Determination of the Robin coefficient in a nonlinear boundary condition for a steady‐state problem

In this paper, a constant heat transfer coefficient present in a nonlinear Robin‐type boundary condition associated with an elliptic equation is reconstructed uniquely from a single boundary energy measurement. Two types of such boundary energy measurement are considered, and solvability theorems for the solution of the resulting nonlinear inverse problems are provided. Further, one‐dimensional numerical results are presented and discussed. Copyright © 2008 John Wiley & Sons, Ltd.