Existence and uniqueness for an equation of mixed type in a rectangle

A boundary value problem is considered for an equation of mixed type in a rectangle. An energy method is used to prove uniqueness for appropriate boundary conditions. Similar considerations for the adjoint problem allow one to prove the existence of a generalized solution to the problem. It is also shown that the energy method leads naturally to a finite element method. The method is described for a special case and convergence is shown for the resulting numerical procedure.     

Local existence and non-relativistic limits of shock solutions to a multidimensional piston problem for the relativistic Euler equations

The multidimensional piston problem is a special initial-boundary value problem. The boundary conditions are given in two conical surfaces: one is the boundary of the piston, and the other is the shock whose location is to be determined later. In this paper, we are concerned with spherically symmetric piston problem for the relativistic Euler equations. A local shock front solution with the state equation p = a ² ρ, a is a constant and has been established by the Newton iteration. To overcome the difficulty caused by the free boundary, we introduce a coordinate transformation to fix it and employ the linear iteration scheme to establish a sequence of approximate solutions to the auxiliary problems by iteration. In each step, the value of the solution of the previous problem is taken as the data to determine the solution of the next problem. We obtain the existence of the original problem by establishing the convergence of these sequences. Meanwhile, we establish the convergence of the local solution as c → ∞ to the corresponding solution of the classical non-relativistic Euler equations.     

Numerical Computation of Connecting Orbits in Delay Differential Equations

We discuss the numerical computation of homoclinic and heteroclinic orbits in delay differential equations. Such connecting orbits are approximated using projection boundary conditions, which involve the stable and unstable manifolds of a steady state solution. The stable manifold of a steady state solution of a delay differential equation (DDE) is infinite-dimensional, a problem which we circumvent by reformulating the end conditions using a special bilinear form. The resulting boundary value problem is solved using a collocation method. We demonstrate results, showing homoclinic orbits in a model for neural activity and travelling wave solutions to the delayed Hodgkin–Huxley equation. Our numerical tests indicate convergence behaviour that corresponds to known theoretical results for ODEs and periodic boundary value problems for DDEs.     

Radial basis functions method for solving of a non-local boundary value problem with Neumann’s boundary conditions

In this paper, the problem of solving the two-dimensional diffusion equation subject to a non-local condition involving a double integral in a rectangular region is considered. The solution of this type of problems are complicated. Therefore, a simple meshless method using the radial basis functions is constructed for the non-local boundary value problem with Neumann’s boundary conditions. Numerical examples are included to demonstrate the reliability and efficiency of this method. Also N_e and root mean square errors are obtained to show the convergence of the method.     

Convergence estimates of a projection-difference method for an operator-differential equation

This article investigates the projection-difference method for a Cauchy problem for a linear operator-differential equation with a leading self-adjoint operator A(t) and a subordinate linear operator K(t) in Hilbert space. This method leads to the solution of a system of linear algebraic equations on each time level; moreover, the projection subspaces are linear spans of eigenvectors of an operator similar to A(t). The convergence estimates are obtained. The application of the developed method for solving the initial boundary value problem is given.     

On an upwind difference scheme for strongly degenerate parabolic equations modelling the settling of suspensions in centrifuges and non-cylindrical vessels

We prove the convergence of an explicit monotone finite difference scheme approximating an initial-boundary value problem for a spatially one-dimensional quasilinear strongly degenerate parabolic equation, which is supplied with two zero-flux boundary conditions. This problem arises in a model of sedimentation–consolidation processes in centrifuges and vessels with varying cross-sectional area. We formulate the definition of entropy solution of the model in the sense of Kružkov and prove the convergence of the scheme to the unique BV entropy solution of the problem. The scheme and the model are illustrated by numerical examples.     

Convergence of an operator splitting method on a bounded domain for a convection-diffusion-reaction system

We solve a convection-diffusion-sorption (reaction) system on a bounded domain with dominant convection using an operator splitting method. The model arises in contaminant transport in groundwater induced by a dual-well, or in controlled laboratory experiments. The operator splitting transforms the original problem to three subproblems: nonlinear convection, nonlinear diffusion, and a reaction problem, each with its own boundary conditions. The transport equation is solved by a Riemann solver, the diffusion one by a finite volume method, and the reaction equation by an approximation of an integral equation. This approach has proved to be very successful in solving the problem, but the convergence properties where not fully known. We show how the boundary conditions must be taken into account, and prove convergence in L1,loc of the fully discrete splitting procedure to the very weak solution of the original system based on compactness arguments via total variation estimates. Generally, this is the best convergence obtained for this type of approximation. The derivation indicates limitations of the approach, being able to consider only some types of boundary conditions. A sample numerical experiment of a problem with an analytical solution is given, showing the stated efficiency of the method.     

A discrete collocation method for Fredholm integro-differential equations with weakly singular kernels

A fully discrete version of a piecewise polynomial collocation method is constructed to solve initial or boundary value problems of linear Fredholm integro-differential equations with weakly singular kernels. Using an integral equation reformulation and special graded grids, optimal global convergence estimates are derived. For special values of parameters an improvement of the convergence rate of elaborated numerical schemes is established. Some of our theoretical results are illustrated by numerical experiments.     

L'équation de Laplace avec conditions aux limites discontinues : convergence d'une discrétisation par éléments spectraux

We prove an existence result for the Laplace equation in a disk with discontinuous Dirichlet boundary conditions: 0 on part of the boundary and 1 on its complement. The problem is discretized by the mortar spectral element method and a convergence estimate is derived which is confirmed numerically.     

Existence and regularity of the solution of a mixed boundary value problem for the Keldysh equation with a nonlinear absorption term

The Keldysh equation is a more general form of the classic Tricomi equation from fluid dynamics. Its well-posedness and the regularity of its solution are interesting and important. The Keldysh equation is elliptic in y>0 and is degenerate at the line y=0 in R². Adding a special nonlinear absorption term, we study a nonlinear degenerate elliptic equation with mixed boundary conditions in a piecewise smooth domain—similar to the potential fluid shock reflection problem. By means of an elliptic regularization technique, a delicate a priori estimate and compact argument, we show that the solution of a mixed boundary value problem of the Keldysh equation is smooth in the interior and Lipschitz continuous up to the degenerate boundary under some conditions. We believe that this kind of regularity result for the solution will be rather useful.     

Order reduction in the integral equations of the null-field method

We reduce the system of equations of the null-field method in the problem of diffraction on a permeable scatterer to a single equation. For a special representation of the field inside the scatterer, we prove the equivalence of that equation to the original boundary value problem of diffraction. We state conditions of the convergence of the approximate solution to the exact one. We prove the completeness and closedness of a system of plane waves.     

Boundary Value Problems with a Turning Point

A boundary value problem involving a second order differential equation on an interval containing a single turning point is considered. Asymptotic approximations to solutions of the differential equation are obtained by the comparison equation method. An essential feature of the work is that domains of validity are restricted so that all approximations are “complete” in the sense of Olver. Asymptotic solutions of the boundary value problem are studied in the special cases where resonance is possible.     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．     

Diffraction of water waves due to a presence of two headlands

This paper is concerned with a boundary value problem for the Helmholtz equation on a horizontal infinite strip with obstacles. The derivation of Helmholtz equation from shallow water equations is given and the boundary value problem with an arbitrary shap of headland is stated. The boundary conditions are of the general Neumann type, and thus we use the finite difference method in numerical solution. Helholtz equation is replaced by the five-points formula and for the points close to the boundary, Taylors expansions are made useful with non-uniform spacing. For solving the resulting system of linear equations, the “Mathematica” package is used. The graphs show the velocity potential contours in the cases, of semielliptic, semicircular and narrow headland. Also, we discuss the problem in the presence of two headlands.     

Convergence of the finite element method for an elliptic equation with strong degeneration

Under consideration are the questions of the numerical solution by the finite element method (FEM) of the first boundary value problem for an elliptic equation with degeneration on a part of the boundary. The weak and strong variational statements are posed in the function spaces with the coordinated weights that correspond to the problem. Using the method of the multiplicative extraction of singularities for the finite element method that utilizes piecewise linear elements, we prove that the convergence of the approximate solutions to the exact solution in the weighted norm is not worse than in the case of an elliptic equation without degeneration.     

Computational methods of solving the boundary value problems for the loaded differential and Fredholm integro‐differential equations

The article presents a new general solution to a loaded differential equation and describes its properties. Solving a linear boundary value problem for loaded differential equation is reduced to the solving a system of linear algebraic equations with respect to the arbitrary vectors of general solution introduced. The system's coefficients and right sides are computed by solving the Cauchy problems for ordinary differential equations. Algorithms of constructing a new general solution and solving a linear boundary value problem for loaded differential equation are offered. Linear boundary value problem for the Fredholm integro‐differential equation is approximated by the linear boundary value problem for loaded differential equation. A mutual relationship between the qualitative properties of original and approximate problems is obtained, and the estimates for differences between their solutions are given. The paper proposes numerical and approximate methods of solving a linear boundary value problem for the Fredholm integro‐differential equation and examines their convergence, stability, and accuracy.     

On an iterative method for the solution of discretized elliptic problems with imperfect contact condition

Iterative method for the solution of discretized elliptic boundary value problem with imperfect contact condition on an internal boundary is constructed. It is proved that the method is robust with respect to the coefficient K appearing in the imperfect contact condition, i.e., its convergence rate is independent of the coefficient K. Numerical tests confirming the theoretical result are presented.     

无界区域Stokes 问题非重叠型区域分解算法及其收敛性

本文研究无界区域Stokes方程外问题的利用有限元法和自然边界归化的非蕈叠型区域分解算法,此方法对无界区域Stokes问题非常有效.给出连续和离散情形的D-N算法及其收敛性分析,得到算法收敛的充要条件及充分条件,并得到最优的松弛因子和压缩因子,最后给出数值算例予以验证.     

Existence and regularity of solution of mixed boundary value problem of a Busemann-equation with nonlinear absorb term

The Busemann-equation is a classical equation coming from fluid dynamics. The well-posed problem and regularity of solution of Busemann-equation with nonlinear term are interesting and important. The Busemann-equation is elliptic in y>0 and is degenerate at the line y=0 in R². With a special nonlinear absorb term, we study a nonlinear degenerate elliptic equation with mixed boundary conditions in a piecewise smooth domain. By means of elliptic regularization technique, a delicate prior estimate and compact argument, we show that the solution of mixed boundary value problem of Busemann-equation is smooth in the interior and Lipschitz continuous up to the degenerate boundary on some conditions. The result is better than the result of classical boundary degenerate elliptic equation.     

On the convergence of a numerical method for a hypersingular integral equation on a closed surface

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises in the solution of the Neumann boundary value problem for the Laplace equation with a representation of a solution in the form of a double-layer potential. We consider the case in which the interior or exterior boundary value problem is solved in a domain; whose boundary is a smooth closed surface, and an integral equation is written out on that surface. For the integral operator in that equation, we suggest quadrature formulas like the method of vortical frames with a regularization, which provides its approximation on the entire surface for the use of a nonstructured partition. We construct a numerical scheme for the integral equation on the basis of suggested quadrature formulas, prove an estimate for the norm of the inverse matrix of the related system of linear equations and the uniform convergence of numerical solutions to the exact solution of the hypersingular integral equation on the grid.