A note on a Cauchy problem for the Laplace equation: Regularization and error estimates

In this paper, a Cauchy problem for the Laplace equation is investigated. Based on the fundamental solution to the elliptic equation, we propose to solve this problem by the truncation method, which generates well-posed problem. Then the well posedness of the proposed regularizing problem and convergence property of the regularizing solution to the exact one are proved. Error estimates for this method are provided together with a selection rule for the regularization parameter. The numerical results show that our proposed numerical methods work effectively. This work extends to earlier results in Qian et al. (2008) [14] and Hao et al. (2009) [5].     

缺陷对结构临界载荷的影响

利用普适开折理论研究了多重缺陷对工程结构临界载荷的影响.对音叉问题给出了由缺陷引起的结构临界载荷改变量的下界估计公式及相应的简单、有效的数值计算方法.     

An iterative method for the reconstruction of a stationary flow

In this article, an iterative algorithm based on the Landweber‐Fridman method in combination with the boundary element method is developed for solving a Cauchy problem in linear hydrostatics Stokes flow of a slow viscous fluid. This is an iteration scheme where mixed well‐posed problems for the stationary generalized Stokes system and its adjoint are solved in an alternating way. A convergence proof of this procedure is included and an efficient stopping criterion is employed. The numerical results confirm that the iterative method produces a convergent and stable numerical solution. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Modified regularization method for the Cauchy problem of the Helmholtz equation

In this paper, the Cauchy problem for the Helmholtz equation is investigated. It is known that such problem is severely ill-posed. We propose a modified regularization method to solve it based on the solution given by the method of separation of variables. Convergence estimates are presented under two different a-priori bounded assumptions for the exact solution. Finally, numerical examples are given to show the effectiveness of the proposed numerical method.     

The numerical treatment of perturbed bifurcation problems in ordinary differential equations

A new numerical method is presented to analyze perturbations of bifurcations of the solutions of nonlinear boundary value problems. The perturbations may result from imperfections or other inhomogeneities in the corresponding scientific problem. The nonisolated solutions are calculated in dependence of the perturbation parameter. Therefore, it is possible to determine the singular solution as well as a solution branch through this nonisolated solution simultaneously. Standard procedures of numerical analysis are applicable.     

Three-dimensional numerical simulation of the inverse problem of thermal convection using the quasi-reversibility method

Inverse (time-reverse) simulation of three-dimensional thermoconvective flows is considered for a highly viscous incompressible fluid with temperature-dependent density and viscosity. The model of the fluid dynamics is described by the Stokes equations, the incompressibility and heat balance equations subject to the appropriate initial and boundary conditions. To solve the problem backward in time, the quasi-reversibility method is applied to the heat balance equation. The numerical solution is based on the introduction of a two-component vector potential for the velocity of the medium, on the application of the finite element method with a special tricubic spline basis for computing this potential, and on the application of the splitting method and the method of characteristics for computing the temperature. The numerical algorithm is designed to be executed on parallel computers. The proposed numerical algorithm is used to reconstruct the evolution of diapiric structures in the Earth’s upper mantle. The computational efficiency of the algorithm is analyzed on the basis of the appropriate functionals of residuals.     

A Fourth-Order Modified Method for the Cauchy Problem of the Modified Helmholtz Equation

This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0 < x≤ 1, y ∈ R. The Cauchy data at x = 0 is given and the solution is then sought for the interval 0 < x ≤1. This problem is highly ill-posed and the solution (if it exists) does not depend continuously on the given data. In this paper, we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution. Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.     

一类各向异性外问题的非重叠型区域分解算法

In this paper, based on the natural integral operator on elliptic boundary, a nonoverlapping domain decomposition method is presented for a kind of anisotropic elliptic problem with constant coefficients in an exterior domain, and the convergence of the method is analyzed. The choice of the relaxition factor is discussed.Some numerical examples are given. Theoretical analysis as well as numerical examples show that our method is performance.     

双调和方程柯西问题的修正的Tikhonov正则化方法

本文研究了双调和方程柯西问题,这类是不适定的,即问题的解(如果存在)不连续依赖于测量数据.首先在精确解的先验假设下给出问题的条件稳定性结果.接着利用修正的Tikhonov正则化方法求解此不适定问题.在先验和后验正则化参数选取规则下,给出正则解和精确解之间的误差估计式.最后给出几个数值例子验证此正则化方法求解此类反问题的有效性.     

Extremal Eigenvalues of the Sturm-Liouville Problems with Discontinuous Coefficients

In this paper, an extremal eigenvalue problem to the Sturm-Liouville equations with discontinuous coefficients and volume constraint is investigated. Liouville transformation is applied to change the problem into an equivalent minimization problem. Finite element method is proposed and the convergence for the finite element solution is established. A monotonic decreasing algorithm is presented to solve the extremal eigenvalue problem. A global convergence for the algorithm in the continuous case is proved. A few numerical results are given to depict the efficiency of the method.     

Analytical solution for the evolution of exit point along the sloping seepage face

This paper developed an analytical solution for the problem of exit point evolution on the seepage face in the unconfined aquifer with sloping interface. A theoretical model for the groundwater drawdown problem in a half‐infinite aquifer with a sloping boundary is built in accordance with the linearized one‐dimensional Boussinesq equation and the Neumann boundary condition at the seepage point. The homotopy analysis method is then adopted for solving this dynamic boundary problem. By constructing two continuous deformations, the original problem could be converted into a group of subproblems with the same physical essence and similar mathematical solutions. To compare this analytical solution, a numerical model based on the finite volume method is developed, which employs adaptive grids to settle the dynamic boundary condition. The comparisons show that the analytical solution agrees with the numerical model well. The results are useful for the quantification of various hydrological problems. The methodology applied in this study is referential for other dynamic boundary problems as well.     

Identification of a time-dependent source term for a time fractional diffusion problem

In this paper, we study an inverse problem of identifying a time-dependent term of an unknown source for a time fractional diffusion equation using nonlocal measurement data. Firstly, we establish the conditional stability for this inverse problem. Then two regularization methods are proposed to for reconstructing the time-dependent source term from noisy measurements. The first method is an integral equation method which formulates the inverse source problem into an integral equation of the second kind; and a prior convergence rate of regularized solutions is derived with a suitable choice strategy of regularization parameters. The second method is a standard Tikhonov regularization method and formulates the inverse source problem as a minimizing problem of the Tikhonov functional. Based on the superposition principle and the technique of finite-element interpolation, a numerical scheme is proposed to implement the second regularization method. One- and two-dimensional examples are carried out to verify efficiency and stability of the second regularization method.     

Eigenvalues and eigenfunctions of the Laplace operator on an equilateral triangle

A boundary value problem for the Laplace equation with Dirichlet and Neumann boundary conditions on an equilateral triangle is transformed to a problem of the same type on a rectangle. This enables us to use, e.g., the cyclic reduction method for computing the numerical solution of the problem. By the same transformation, explicit formulae for all eigenvalues and all eigenfunctions of the corresponding operator are obtained.     

A Herglotz wavefunction method for solving the inverse Cauchy problem connected with the Helmholtz equation

This paper is concerned with the Cauchy problem connected with the Helmholtz equation. On the basis of the denseness of Herglotz wavefunctions, we propose a numerical method for obtaining an approximate solution to the problem. We analyze the convergence and stability with a suitable choice of regularization method. Numerical experiments are also presented to show the effectiveness of our method.     

The a posteriori Fourier method for solving the Cauchy problem for the Laplace equation with nonhomogeneous Neumann data

In the present paper, the Cauchy problem for the Laplace equation with nonhomogeneous Neumann data in an infinite “strip” domain is considered. This problem is severely ill-posed, i.e., the solution does not depend continuously on the data. A conditional stability result is given. A new a posteriori Fourier method for solving this problem is proposed. The corresponding error estimate between the exact solution and its regularization approximate solution is also proved. Numerical examples show the effectiveness of the method and the comparison of numerical effect between the a posteriori and the a priori Fourier method are also taken into account.     

Recovery of the unknown coefficients in a two-dimensional hyperbolic equation

In this paper, an inverse boundary value problem for a two-dimensional hyperbolic equation with overdetermination conditions is studied. To investigate the solvability of the original problem, we first consider an auxiliary inverse boundary value problem and prove its equivalence to the original problem in a certain sense. We then use the Fourier method to reduce such an equivalent problem to a system of integral equations. Furthermore, we prove the existence and uniqueness theorem for the auxiliary problem by the contraction mappings principle. Based on the equivalency of these problems, the existence and uniqueness theorem for the classical solution of the original inverse problem is proved. Some discussions on the numerical solutions for this inverse problem are presented including some numerical examples.     

A Spectral Element/Laguerre Coupled Method to the Elliptic Helmholtz Problem on the Half Line

A Legendre spectral element/Laguerre coupled method is proposed to numerically solve the elliptic Helmholtz problem on the half line. Rigorous analysis is carried out to establish the convergence of the method. Several numerical examples are provided to confirm the theoretical results. The advantage of this method is demonstrated by a numerical comparison with the pure Laguerre method.     

Finite element approximation with numerical integration for differential eigenvalue problems

Error estimates of the finite element method with numerical integration for differential eigenvalue problems are presented. More specifically, refined results on the eigenvalue dependence for the eigenvalue and eigenfunction errors are proved. The theoretical results are illustrated by numerical experiments for a model problem.     

A general filter regularization method to solve the three dimensional Cauchy problem for inhomogeneous Helmholtz-type equations: Theory and numerical simulation

In this paper, we solve the Cauchy problem for an inhomogeneous Helmholtz-type equation with homogeneous Dirichlet and Neumann boundary condition. The proposed problem is ill-posed. Up to now, most investigations on this topic focus on very specific cases, and with Dirichlet boundary condition. Recently, we solve this problem in 2D for an inhomogeneous modified Helmholtz equation (2012). This work is a continuous expansion of our previous results. Herein we introduce a general filter regularization (GFR) method, and then from the GFR we deduce two concrete filters, which are a foundation to implement a numerical procedure. In addition, we develop a numerical model for solving this problem in three dimensional region. The proposed filter method has been verified by numerical experiments.