A regularization method for solving the Cauchy problem for the Helmholtz equation

In this paper, we investigate a Cauchy problem associated with Helmholtz-type equation in an infinite “strip”. This problem is well known to be severely ill-posed. The optimal error bound for the problem with only nonhomogeneous Neumann data is deduced, which is independent of the selected regularization methods. A framework of a modified Tikhonov regularization in conjunction with the Morozov’s discrepancy principle is proposed, it may be useful to the other linear ill-posed problems and helpful for the other regularization methods. Some sharp error estimates between the exact solutions and their regularization approximation are given. Numerical tests are also provided to show that the modified Tikhonov method works well.     

A regularization method for the cauchy problem of the modified Helmholtz equation

In the present paper, an iteration regularization method for solving the Cauchy problem of the modified Helmholtz equation is proposed. The a priori and a posteriori rule for choosing regularization parameters with corresponding error estimates between the exact solution and its approximation are also given. The numerical example shows the effectiveness of this method. Copyright © 2014 John Wiley & Sons, Ltd.     

Backward heat equation with time dependent variable coefficient

Backward heat equation with time dependent variable coefficient is severely ill‐posed in the sense of Hadamard, so we need regularization. In this paper, we consider Backward heat equation with time dependent variable coefficient, and by small perturbing, we obtain an approximation problem. We show this approximation problem is well‐posed with small parameter. Also, we show this approximation system converges to the original problem when parameter goes to zero. Here, we use modified‐quasi boundary value method to regularize this problem. Copyright © 2016 John Wiley & Sons, Ltd.     

Stable sequential Lagrange principles in the inverse final observation problem for the system of Maxwell equations in the quasistationary magnetic approximation

We justify the possibility of using stable, with respect to errors in the input data, algorithms of dual regularization and iterative dual regularization for solving the inverse final observation problem for the system of Maxwell equations in the quasistationary magnetic approximation under general conditions on the coefficients, which is treated as an optimal control problem for the differential equation describing the magnetic field intensity with an operator equality constraint. We state a classical parametric Lagrange principle and stable Lagrange principles in sequential form for the posed problem. We present a stopping rule for the iterative process for the stable sequential Lagrange principle in iterative form in the case of finite fixed error in the input data.     

The numerical inversion of the Laplace transform in reproducing kernel Hilbert spaces for ill-posed problems

In this paper we have converted the Laplace transform into an integral equation of the first kind of convolution type, which is an ill-posed problem, and used a statistical regularization method to solve it. The method is applied to three examples. It gives a good approximation to the true solution and compares well with the method given by Rodriguez.     

Lavrentiev regularization + Ritz approximation = uniform finite element error estimates for differential equations with rough coefficients

We consider a parametric family of boundary value problems for a diffusion equation with a diffusion coefficient equal to a small constant in a subdomain. Such problems are not uniformly well-posed when the constant gets small. However, in a series of papers, Bakhvalov and Knyazev have suggested a natural splitting of the problem into two well-posed problems. Using this idea, we prove a uniform finite element error estimate for our model problem in the standard parameter-independent Sobolev norm. We also study uniform regularity of the transmission problem, needed for approximation. A traditional finite element method with only one additional assumption, namely, that the boundary of the subdomain with the small coefficient does not cut any finite element, is considered.

One interpretation of our main theorem is in terms of regularization. Our FEM problem can be viewed as resulting from a Lavrentiev regularization and a Ritz-Galerkin approximation of a symmetric ill-posed problem. Our error estimate can then be used to find an optimal regularization parameter together with the optimal dimension of the approximation subspace.

     

An FEM approximation for a fourth-order variational inequality of second kind

A fourth-order variational inequality of the second kind arising in a plate frictional bending problem is considered. By using regularization method, the original problem can be formulated as a differentiable variational equation, and the corresponding discrete FEM variational equation is presented afterwards. Abstract error estimates and error estimates of the approximation are derived in terms of energy norm and L^2-norm.     

On the Stochastic Wave Equation with Nonlinear Damping

We discuss an initial boundary value problem for the stochastic wave equation with nonlinear damping. We establish the existence and uniqueness of a solution. Our method for the existence of pathwise solutions consists of regularization of the equation and data, the Galerkin approximation and an elementary measure-theoretic argument. We also prove the existence of an invariant measure when the equation has pure nonlinear damping.     

On a numerical solution of the Cauchy problem for the Laplace equation by the fundamental solutions method

We discuss a numerical method to solve a Cauchy problem for the Laplace equation in the two-dimensional annular domain. We consider the case that the Cauchy data is given on an arc. We develop an approximation method based of the fundamental solutions method using the least squares method with Tikhonov regularization. The effectiveness of our method is examined by a numerical experiment. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Small mass asymptotic for the motion with vanishing friction

We consider the small mass asymptotic (Smoluchowski–Kramers approximation) for the Langevin equation with a variable friction coefficient. The friction coefficient is assumed to be vanishing within certain region. We introduce a regularization for this problem and study the limiting motion for the 1-dimensional case and a multidimensional model problem. The limiting motion is a Markov process on a projected space. We specify the generator and the boundary condition of this limiting Markov process and prove the convergence.     

The Coupling of Boundary Integral and Finite Element Methods for the Navier-Stokes Equations in an Exterior Domain

In this paper, a technique of coupling variational formulation of FEM and BIE (boundary integral equation) is used to deal with stationary Navier-Stokes equations in an unbounded domain. We discuss well-posedness for the coupling variational problem, the regularization method and FEM-BEM approximation. Finally, operator splitting and optimal control techniques are used to treat the difficulty of nonlinearity and constraints in computer implementation.     

Convergence of a numerical scheme of the type of the vortex loop method on a closed surface with an approximation to the surface shape

We consider a linear integral equation with a hypersingular integral treated in the sense of the Hadamard finite value. This equation arises when solving the Neumann boundary value problem for the Laplace equation with the use of the representation of the solution in the form of a double layer potential. We study the case in which an exterior or interior boundary value problem is solved in a domain whose boundary is a smooth closed surface and the integral equation is written out on that surface. For the numerical solution of the integral equation, the surface is approximated by spatial polygons whose vertices lie on the surface. We construct a numerical scheme for solving the integral equation on the basis of such an approximation to the surface with the use of quadrature formulas of the type of the method of discrete singularities with regularization. We prove that the numerical solutions converge to the exact solution of the hypersingular integral equation uniformly on the grid.     

带损伤弹性反问题的数值分析

考虑一类由椭圆性方程和热传导方程共同来刻画的准静态弹性模型,通过给定观测值来反演边界的牵引力.首先构造一个凸目标泛函,并引入Tikhonov正则化方法,使之极小化得到一个稳定的近似解.再用有限元离散求解,导出误差估计.最后,用数值例子说明算法的可行性和有效性.     

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.     

On the convergence of the vortex loop method with regularization for the Neumann boundary value problem on a plane screen

We consider a linear integral equation with a supersingular integral treated in the sense of the Hadamard finite value, which arises in the solution of the Neumann boundary value problem for the Laplace equation with the representation of the solution in the form of a doublelayer potential. We consider the case in which the exterior boundary value problem is solved outside a plane surface (a screen). For the integral operator in the above-mentioned equation, we suggest quadrature formulas of the vortex loop method with regularization, which provide its approximation on the entire surface when using an unstructured partition. In the problem in question, the derivative of the unknown density of the double-layer potential, as well as the errors of quadrature formulas, has singularities in a neighborhood of the screen edge. We construct a numerical scheme for the integral equation on the basis of the suggested quadrature formulas and prove an estimate for the norm of the inverse matrix of the resulting system of linear equations and the uniform convergence of the numerical solutions to the exact solution of the supersingular integral equation on the grid.     

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.     

Optimal control for systems described by hyperbolic equations with strong nonlinearity

An optimization control problem for a hyperbolic equation is considered. The system is nonlinear with respect to the state derivative. The regularization technique for the state equation is applied. The necessary conditions of optimality for the regularized control problem are proved. It uses the extended differentiability of the control-state mapping for the regularized equation. The convergence of the regularization method is proved. Thus the optimal control for the regularized problem with a small enough regularization parameter can be chosen as an approximate solution of the initial optimization problem.     

Tikhonov regularization method for a backward problem for the time-fractional diffusion equation

This paper is devoted to solve a backward problem for a time-fractional diffusion equation with variable coefficients in a general bounded domain by the Tikhonov regularization method. Based on the eigenfunction expansion of the solution, the backward problem for searching the initial data is changed to solve a Fredholm integral equation of the first kind. The conditional stability for the backward problem is obtained. We use the Tikhonov regularization method to deal with the integral equation and obtain the series expression of solution. Furthermore, the convergence rates for the Tikhonov regularized solution can be proved by using an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule. Two numerical examples in one-dimensional and two-dimensional cases respectively are investigated. Numerical results show that the proposed method is effective and stable.     

Adaptivity with relaxation for ill-posed problems and global convergence for a coefficient inverse problem

A new framework of the functional analysis is developed for the finite element adaptive method (adaptivity) for the Tikhonov regularization functional for some ill-posed problems. As a result, the relaxation property for adaptive mesh refinements is established. An application to a multidimensional coefficient inverse problem for a hyperbolic equation is discussed. This problem arises in the inverse scattering of acoustic and electromagnetic waves. First, a globally convergent numerical method provides a good approximation for the correct solution of this problem. Next, this approximation is enhanced via the subsequent application of the adaptivity. Analytical results are verified computationally. Bibliography: 30 titles. Illustration: 2 figures.     

A fractional Tikhonov method for solving a Cauchy problem of Helmholtz equation

In this paper, we study a fractional Tikhonov regularization method (FTRM) for solving a Cauchy problem of Helmholtz equation in the frequency domain. On the one hand, the FTRM retains the advantage of classical Tikhonov method. On the other hand, our method can prevent the effect of oversmoothing of classical Tikhonov method and conveniently control the amount of damping. The convergence error estimates between the exact solution and its regularization approximation are constructed. Several interesting numerical examples are provided, which validate the effectiveness of the proposed method.