Numerical method using cubic B-spline for the heat and wave equation

In the present paper, the cubic B-splines method is considered for solving one-dimensional heat and wave equations. A typical finite difference approach had been used to discretize the time derivative while the cubic B-spline is applied as an interpolation function in the space dimension. The accuracy of the method for both equations is discussed. The efficiency of the method is illustrated by some test problems. The numerical results are found to be in good agreement with the exact solution.     

A logarithmically improved regularity criterion for the surface quasi-geostrophic equation

In this paper, we consider the two-dimensional dissipative surface quasi-geostrophic equation, and establish a logarithmically improved regularity criterion. Consequently, our result extends the regularity criterion result of Dong and Pavlovié (2009). In addition, a logarithmically improved regularity criterion to the three-dimensional Navier–Stokes equations is also derived by the same arguments. Therefore, this result extends and improves many previous works.     

Multilevel difference schemes for the heat conduction equation and its application to the dirichlet problem in two and three dimensions

In this paper we have discussed a general method of deriving high order multilevel schemes for the heat conduction equation in higher dimensions. The application of the three level parabolic schemes for the solution of the steady state Dirichlet problem in two and three dimensions is discussed by introducing two parameters. The error bounds given by Hadjidimos are further refined which drastically reduces the theoretical estimate of the number of iteration cycles required for a given accuracy η. We have also defined a new set of iteration parameters. A numerical example is computed and it is seen that the present methods produce accurate results.     

The Bang-Bang principle of time optimal controls for the heat equation with internal controls

In this paper, we study a certain time optimal control problem governed by the internal controlled heat equation and with a closed ball centered at 0 as the target set. We derive, by making use of Pontryagin's maximum principle and a special kind of unique continuation property for solutions of the heat equation, the bang-bang property for time optimal controls of the problem.     

Continuous dependence of optimal control to controlled domain of actuator for heat equation

In this paper, we consider continuous dependence of the optimal control with respect to the actuator domain which is varying as open subset in the spatial domain for a multi-dimensional heat equation. Both time optimal control and norm optimal control problems are considered. The reason behind combining these two problems together is that these two problems are actually equivalent: The energy to be used to drive the system to target set in minimal time interval is actually the minimal energy of driving the system to target set in this minimal time interval, and visa versa. It is shown that both optimal control and optimal cost are continuous with respect to open controlled actuator domain under the Lebesgue measure.     

A new unconditionally stable ADI compact scheme for the two-space-dimensional linear hyperbolic equation

We formulate a new alternating direction implicit compact scheme of O(τ²+h ⁴) for the linear hyperbolic equation u _tt +2α u _t +β² u=u _xx +u _yy +f(x, y, t), 0<x, y<1, 0<t≤T, subject to appropriate initial and Dirichlet boundary conditions, where α>0 and β≥0 are real numbers. In this article, we show the method is unconditionally stable by the Von Neumann method. At last, numerical demonstrations are given to illustrate our result.     

A quasi-reversibility regularization method for the Cauchy problem of the Helmholtz equation

In this paper, a Cauchy problem for the Helmholtz equation is considered. It is known that such a problem is severely ill-posed, i.e. the solution does not depend continuously on the given Cauchy data. We propose a quasi-reversibility regularization method to solve it. Convergence estimates are established under two different a priori assumptions for an exact solution. Numerical results obtained by two different schemes show that our proposed methods work well.     

Domain decomposition algorithm based on the group explicit formula for the heat equation

A finite difference domain decomposition algorithm (DDA) for solving the heat equation in parallel is presented. In this procedure, interface values between subdomains are calculated by the group explicit formula, whereas interior values of subdomains are determined by the classical implicit scheme. The stability and convergence for this DDA are proved. The stability bound of the procedure is derived to be eight times that of the classical explicit scheme. Though the truncation error at the interface is O(τ?+?h), L ²-error is proved to be O(τ?+?h ²). Numerical examples confirm the second-order convergence and indicate that the stability condition is sharp. A comparison of the numerical errors of this procedure with other known methods is also included.     

A new filtering method for the Cauchy problem of the Laplace equation

In the present paper, we consider a Cauchy problem for the Laplace equation in a rectangle domain. A new filtering method is presented for approximating the solution of this problem, and the Hölder-type error estimates are obtained by the different parameter choice rules. Numerical illustration shows that the proposed method works effectively.     

Numerical study of a nonlinear heat equation for plasma physics

This paper is devoted to the numerical approximation of a nonlinear parabolic balance equation, which describes the heat evolution of a magnetically confined plasma in the edge region of a tokamak. The nonlinearity implies some numerical difficulties, in particular for the long-time behaviour approximation, when solved with standard methods. An efficient numerical scheme is presented in this paper, based on a combination of a directional splitting scheme and the implicit–explicit scheme introduced in Filbet and Jin [A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources, J. Comput. Phys. 229 (2010), pp. 7625–7648].     

Solvability of a non-linear Cauchy problem for an elliptic equation

ABSTRACT

We study a non-linear operator equation originating from a Cauchy problem for an elliptic equation. The problem appears in applications where surface measurements are used to calculate the temperature below the earth surface. The Cauchy problem is ill-posed and small perturbations to the used data can result in large changes in the solution. Since the problem is non-linear certain assumptions on the coefficients are needed. We reformulate the problem as an non-linear operator equation and show that under suitable assumptions the operator is well-defined. The proof is based on making a change of variables and removing the non-linearity from the leading term of the equation. As a part of the proof we obtain an iterative procedure that is convergent and can be used for evaluating the operator. Numerical results show that the proposed procedure converges faster than a simple fixed point iteration for the equation in the the original variables.     

Reconstruction of the corrosion boundary for the Laplace equation by using a boundary collocation method

In this paper, we consider the identification of a corrosion boundary for the two-dimensional Laplace equation. A boundary collocation method is proposed for determining the unknown portion of the boundary from the Cauchy data on a part of the boundary. Since the resulting matrix equation is badly ill-conditioned, a regularized solution is obtained by employing the Tikhonov regularization technique, while the regularization parameter is provided by the generalized cross-validation criterion. Numerical examples show that the proposed method is reasonable and feasible.     

Uniqueness and numerical scheme for the Robin coefficient identification of the time-fractional diffusion equation

We study an inverse problem of determining the Robin coefficient of fractional diffusion equation from a nonlocal boundary condition. Based on the property of Caputo fractional derivative, the uniqueness is proved. The numerical schemes for the direct problem and the inverse problem are developed. Three examples are given to show the effectiveness of the presented methods.     

Probability method for cryptanalysis of general multivariate modular linear equation

Finding the solution to a general multivariate modular linear equation plays an important role in cryptanalysis field. Earlier results show that obtaining a relatively short solution is possible in polynomial time. However, one problem arises here that if the equation has a short solution in given bounded range, the results outputted by earlier algorithms are often not the ones we are interested in. In this paper, we present a probability method based on lattice basis reduction to solve the problem. For a g...     

Numerical procedure for the determination of an unknown parameter in a parabolic equation

A numerical procedure for an inverse problem of determining unknown source parameter of one-dimensional parabolic equation subject to the specification of the solution at internal point along with the usual initial boundary conditions is considered. By using some transformation the problem is reformulated to a nonlocal parabolic problem. Some numerical examples using the proposed numerical procedure are presented.     

The truncation method for identifying the heat source dependent on a spatial variable

This paper deals with the ill-posed problem of determining an unknown source depending only on the spatial variable in the heat equation. In this problem, the truncation method is employed to get the regularization solution. Moreover, the Hölder type error estimate is also obtained between this regularization solution and the exact solution. Various typical examples are given to verify the efficiency and accuracy of the proposed method.     

Identification of variable spacial coefficients for a beam equation from boundary measurements

We consider a system described by the Euler-Bernoulli beam equation, with one end clamped and with torque input at the other end. The output function are the displacement and the angle velocity at the non-clamped end of the beam. We study the identification of the spatially variable coefficients in the beam equation, from input-output data. We show that both the density and the flexural rigidity of the beam (which are assumed to be of class C⁴) can be uniquely determined if the input and output functions are known for all positive times.