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.     

Sharp estimates for approximations to a nonlinear backward heat equation

A nonlinear backward heat problem for an infinite strip is considered. The problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. In this paper, we use the Fourier regularization method to solve the problem. Some sharp estimates of the error between the exact solution and its regularization approximation are given.     

FOURIER REGULARIZATION FOR DETERMINING SURFACE HEAT FLUX FROM INTERIOR OBSERVATION BASED ON A SIDEWAYS PARABOLIC EQUATION

In this paper we consider a non-standard inverse heat conduction problem for determining surface heat flux from an interior observation which appears in some applied subjects. This problem is ill-posed in the sense that the solution (if it exists) does not depend continuously on the data. A Fourier method is applied to formulate a regularized approximation solution, and some sharp error estimates are also given.     

Mixed problem for a first-order partial differential equation with involution and periodic boundary conditions

The Fourier method is used to find a classical solution of the mixed problem for a first-order differential equation with involution and periodic boundary conditions. The application of the Fourier method is substantiated using refined asymptotic formulas obtained for the eigenvalues and eigenfunctions of the corresponding spectral problem. The Fourier series representing the formal solution is transformed using certain techniques, and the possibility of its term-by-term differentiation is proved. Minimal requirements are imposed on the initial data of the problem.     

Splitting method for an inverse source problem in parabolic differential equations: Error analysis and applications

In this work, we present a numerical method based on a splitting algorithm to find the solution of an inverse source problem with the integral condition. The source term is reconstructed by using the specified data and by employing the Lie splitting method, we decompose the equation into linear and nonlinear parts. Each subproblem is solved by the Fourier transform and then by combining the solutions of subproblems, the solution of the original problem is computed. Moreover, the framework of strongly continuous semigroup (or C₀-semigroup) is employed in error analysis of operator splitting method for the inverse problem. The convergence of the proposed method is also investigated and proved. Finally, some numerical examples in one, two, and three-dimensional spaces are provided to confirm the efficiency and capability of our work compared with some other well-known methods.     

Two regularization methods for a spherically symmetric inverse heat conduction problem

In this paper, we consider a spherically symmetric inverse heat conduction problem of determining the internal surface temperature of a hollow sphere from the measured data at a fixed location inside it. This is an ill-posed problem in the sense that the solution (if it exists) does not depend continuously on the data. A Tikhonov type’s regularization method and a Fourier regularization method are applied to formulate regularized solutions which are stably convergent to the exact ones with order optimal error estimates.     

Fourier Moment Method with Regularization for the Cauchy Problem of Helmholtz Equation

In this paper, we consider the reconstruction of the wave field in a bounded domain. By choosing a special family of functions, the Cauchy problem can be transformed into a Fourier moment problem. This problem is ill-posed. We propose a regularization method for obtaining an approximate solution to the wave field on the unspecified boundary. We also give the convergence analysis and error estimate of the numerical algorithm. Finally, we present some numerical examples to show the effectiveness of this method.     

Application of generalized separation of variables to solving mixed problems with irregular boundary conditions

One of the methods for solving mixed problems is the classical separation of variables (the Fourier method). If the boundary conditions of the mixed problem are irregular, this method, generally speaking, is not applicable. In the present paper, a generalized separation of variables and a way of application of this method to solving some mixed problems with irregular boundary conditions are proposed. Analytical representation of the solution to this irregular mixed problem is obtained.     

A characterisation of oscillations in the discrete two-dimensional convection-diffusion equation

It is well known that discrete solutions to the convection-diffusion equation contain nonphysical oscillations when boundary layers are present but not resolved by the discretisation. However, except for one-dimensional problems, there is little analysis of this phenomenon. In this paper, we present an analysis of the two-dimensional problem with constant flow aligned with the grid, based on a Fourier decomposition of the discrete solution. For Galerkin bilinear finite element discretisations, we derive closed form expressions for the Fourier coefficients, showing them to be weighted sums of certain functions which are oscillatory when the mesh Péclet number is large. The oscillatory functions are determined as solutions to a set of three-term recurrences, and the weights are determined by the boundary conditions. These expressions are then used to characterise the oscillations of the discrete solution in terms of the mesh Péclet number and boundary conditions of the problem.

     

A wavelets method for plane elasticity problem

A Neumann boundary value problem of plane elasticity problem in the exterior circular domain is reduced into an equivalent natural boundary integral equation and a Poisson integral formula with the DtN method. Using the trigonometric wavelets and Galerkin method, we obtain a fast numerical method for the natural boundary integral equation which has an unique solution in the quotient space. We decompose the stiffness matrix in our numerical method into four circulant and symmetrical or antisymmetrical submatrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform (FFT) and the inverse fast Fourier transform (IFFT) instead of the inverse matrix. Examples are given for demonstrating our method has good accuracy of our method even though the exact solution is almost singular.     

关于波动方程混合问题的特征线方法

传统的求解1维波动方程混合问题的方法是分离变量法，进而解出该问题的Fourier级数解，本文将用特征线方法给出该问题在求解区域内解的显式表达式。     

The fast Fourier transform and the numerical solution of one-dimensional boundary integral equations

Summary Here we present a fully discretized projection method with Fourier series which is based on a modification of the fast Fourier transform. The method is applied to systems of integro-differential equations with the Cauchy kernel, boundary integral equations from the boundary element method and, more generally, to certain elliptic pseudodifferential equations on closed smooth curves. We use Gaussian quadratures on families of equidistant partitions combined with the fast Fourier transform. This yields an extremely accurate and fast numerical scheme. We present complete asymptotic error estimates including the quadrature errors. These are quasioptimal and of exponential order for analytic data. Numerical experiments for a scattering problem, the clamped plate and plane estatostatics confirm the theoretical convergence rates and show high accuracy.     

Analysis of pseudo-periodic chronological series with irregularly time-spaced data in view to their prediction. II. Fourier analysis

The problem originates from the necessity to predict luminosities of large-amplitude variable stars that are to be observed by the astronomical satellite HIPPARCOS. The data have a specific character: they are unequally time-spaced and can be missing during a long time in comparison to the pseudo-period. So the classical method of time-series analysis must be adapted and new methods are to be searched. In this paper we present a numerical solution derived from a Fourier analysis.     

Spectral smoothed boundary methods: The role of external boundary conditions

The spectral smoothed boundary method (SSBM) is a recently proposed numerical method to approximate the solution of partial differential equations in irregular domains with no‐flux boundary conditions by means of Fourier spectral methods. In this article we explore the robustness and accuracy of the scheme under variations of the artificial boundary conditions that must be imposed on the boundary of the enlarged domain in which the problem is solved. As a test model, we present quantitative numerical results based on a problem of propagation of waves of electrical activity in cardiac tissue for which the method is relevant. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Dirichlet problem for a third-order equation of mixed type in a rectangular domain

For a third-order differential equation of parabolic-hyperbolic type, we suggest a method for studying the first boundary value problem by solving an inverse problem for a second-order equation of mixed type with unknown right-hand side. We obtain a uniqueness criterion for the solution of the inverse problem. The solution of the inverse problem and the Dirichlet problem for the original equation is constructed in the form of the sum of a Fourier series.     

关于一类具阻尼项的非线性波方程的Cauchy问题

The paper concerns with the existence, uniqueness and nonexistence of global solution to the Cauchy problem for a class of nonlinear wave equations with damping term. It proves that under suitable assumptions on nonlinear the function and initial data the above-mentioned problem admits a unique global solution by Fourier transform method. The sufficient conditions of nonexistence of the global solution to the above-mentioned problem are given by the concavity method.     

Approximation of an Inverse Initial Problem for a Biparabolic Equation

In this paper, we consider the problem of finding the initial distribution for the linear inhomogeneous and nonlinear biparabolic equation. The problem is severely ill-posed in the sense of Hadamard. First, we apply a general filter method to regularize the linear nonhomogeneous problem. Then, we also give a regularized solution and consider the convergence between the regularized solution and the sought solution. Under the a priori assumption on the exact solution belonging to a Gevrey space, we consider a generalized nonlinear problem by using the Fourier truncation method to obtain rigorous convergence estimates in the norms on Hilbert space and Hilbert scale space.     

A mollification method for ill-posed problems

Summary. A mollification method for a class of ill-posed problems is suggested. The idea of the method is very simple and natural: if the data are given inexactly then we try to find a sequence of ``mollification operators" which map the improper data into well-posedness classes of the problem (mollify the improper data). Within these mollified data our problem becomes well-posed. And when these facts are in hand we try to obtain error estimates and optimal or ``quasi-optimal" mollification parameters. The method is working not only for problems in Hilbert spaces, but also for problems in Banach spaces. Applications of the method to concrete problems, like numerical differentiation, parabolic equations backwards in time, the Cauchy problem for the Laplace equation, one- and multidimensional non-characteristic Cauchy problems for parabolic equations (in infinite or finite domains),... give us very sharp stability estimates of H\"older continuous type. In these cases the method is optimal in the sense that it gives the same order of H\"older continuous dependence on the data as for the regularized problems. Furthermore, the method may be implemented numerically using fast Fourier transforms. For the first time a uniform stability estimate of H\"older continuous type of the solution of the heat equation backwards in time in the space for all could be established by our mollification method. A new simple sharp pointwise estimate of H\"older type for the weak solution of a non-characteristic Cauchy problem for parabolic equations in a finite domain is established. Received June 25, 1993 / Revised version received February 18, 1994     

The Fourier-finite element method for the Poisson problem on a non-convex polyhedral cylinder

We study the Poisson problem with zero boundary datum in a (finite) polyhedral cylinder with a non-convex edge. Applying the Fourier sine series to the equation along the edge and by a corner singularity expansion for the Poisson problem with parameter, we define the edge flux coefficient and the regular part of the solution on the polyhedral cylinder. We present a numerical method for approximating the edge flux coefficient and the regular part and show the stability. We derive an error estimate and give some numerical experiments.     

A mixed problem for an inhomogeneous wave equation with a summable potential

A mixed problem for an inhomogeneous wave equation with fixed ends in the case of a summable potential is studied. Using the Krylov method for acceleration of the convergence of Fourier series, a classical solution under minimal conditions on the initial data and a generalized solution in the case of quadratic summable initial data and perturbing function are obtained.