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.     

A LOCKING-FREE ANISOTROPIC NONCONFORMING FINITE ELEMENT FOR PLANAR LINEAR ELASTICITY PROBLEM

The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value prob-lem of planar linear elasticity. The optimal error estimates are obtained by using some novel approaches and techniques. The method proposed in this article is robust in the sense that the convergence estimates in the energy and L2-norms are independent of the Lame parameterλ.     

Boundary Identification for a Blast Furnace

In this paper,the authors discuss an inverse boundary problem for the axi- symmetric steady-state heat equation,which arises in monitoring the boundary corrosion for the blast-furnace.Measure temperature at some locations are used to identify the shape of the corrosion boundary. The numerical inversion is complicated and consuming since the wear-line varies during the process and the boundary in the heat problem is not fixed.The authors suggest a method that the unknown boundary can be represented by a given curve plus a small perturbation,then the equation can be solved with fixed boundary,and a lot of computing time will be saved. A method is given to solve the inverse problem by minimizing the sum of the squared residual at the measuring locations,in which the direct problems are solved by axi- symmetric fundamental solution method. The numerical results are in good agreement with test model data as well as industrial data,even in severe corrosion case.     

NATURAL BOUNDARY ELEMENT METHOD FOR THREE DIMENSIONAL EXTERIOR HARMONIC PROBLEM WITH AN INNER PROLATE SPHEROID BOUNDARY

In this paper, we study natural boundary reduction for Laplace equation with Dirichletor Neumann boundary condition in a three-dimensional unbounded domain, which is theoutside domain of a prolate spheroid. We express the Poisson integral formula and naturalintegral operator in a series form explicitly. Thus the original problem is reduced to aboundary integral equation on a prolate spheroid. The variational formula for the reducedproblem and its well-posedness are discussed. Boundary element approximation for thevariational problem and its error estimates, which have relation to the mesh size andthe terms after the series is truncated, are also presented. Two numerical examples arepresented to demonstrate the effectiveness and error estimates of this method.     

An Inverse Problem of Identifying the Radiative Coefficient in a Degenerate Parabolic Equation

The authors investigate an inverse problem of determining the radiative coefficient in a degenerate parabolic equation from the final overspecified data. Being different from other inverse coefficient problems in which the principle coefficients are assumed to be strictly positive definite, the mathematical model discussed in this paper belongs to the second order parabolic equations with non-negative characteristic form, namely, there exists a degeneracy on the lateral boundaries of the domain. Based on the optimal control framework, the problem is transformed into an optimization problem and the existence of the minimizer is established. After the necessary conditions which must be satisfied by the minimizer are deduced, the uniqueness and stability of the minimizer are proved. By minor modification of the cost functional and some a priori regularity conditions imposed on the forward operator, the convergence of the minimizer for the noisy input data is obtained in this paper. The results can be extended to more general degenerate parabolic equations.     

Efficient estimation for semiparametric varying- coefficient partially linear regression models with current status data

This article considers a semiparametric varying-coefficient partially linear regression model with current status data. The semiparametric varying-coefficient partially linear regression model which is a generalization of the partially linear regression model and varying-coefficient regression model that allows one to explore the possibly nonlinear effect of a certain covariate on the response variable. A Sieve maximum likelihood estimation method is proposed and the asymptotic properties of the proposed estimators are discussed. Under some mild conditions, the estimators are shown to be strongly consistent. The convergence rate of the estimator for the unknown smooth function is obtained and the estimator for the unknown parameter is shown to be asymptotically efficient and normally distributed. Simulation studies are conducted to examine the small-sample properties of the proposed estimates and a real dataset is used to illustrate our approach.     

GLOBAL EXISTENCE AND CONVERGENCE RATES OF SMOOTH SOLUTIONS FOR THE 3-D COMPRESSIBLE MAGNETOHYDRODYNAMIC EQUATIONS WITHOUT HEAT CONDUCTIVITY

In this paper, we are concerned with the global existence and convergence rates of the smooth solutions for the compressible magnetohydrodynamic equations without heat conductivity, which is a hyperbolic-parabolic system. The global solutions are obtained by combining the local existence and a priori estimates if H3-norm of the initial perturbation around a constant states is small enough and its L1-norm is bounded. A priori decay-in-time estimates on the pressure, velocity and magnetic field are used to get the uniform bound of entropy. Moreover, the optimal convergence rates are also obtained.     

A numerical algorithm for determination of a control parameter in two-dimensional parabolic inverse problems

Numerical solution of the parabolic partial differential equations with an unknown parameter play a very important role in engineering applications. In this study we present a high order scheme for determining unknown control parameter and unknown solution of two-dimensional parabolic inverse problem with overspecialization at a point in the spatial domain. In this approach, a compact fourth-order scheme is used to discretize spatial derivatives of equation and reduces the problem to a system of ordinary differential equations(ODEs).Then we apply a fourth order boundary value method to the solution of resulting system of ODEs. So the proposed method has fourth order of accuracy in both space and time components and is unconditionally stable due to the favorable stability property of boundary value methods. The results of numerical experiments are presented and some comparisons are made with several well-known finite difference schemes in the literature.Also we will investigate the effect of noise in data on the approximate solutions.     

非线性抛物型方程有限元法数值积分的有效性

Abstract. The effect of numerical integration in finite element methods applied to a class of nonlinear parabolic equations is considered and some sufficient conditions on the quadrature scheme to ensure that the order of convergence is unaltered in the presence of numerical integration are given. Optimal Lz and H1 estimates for the error and its time derivative are established.     

THE OPTIMAL CONVERGENCE ORDER OF THE DISCONTINUOUS FINITE ELEMENT METHODS FOR FIRST ORDER HYPERBOLIC SYSTEMS

In this paper, a discontinuous finite element method for the positive and symmetric, first-order hyperbolic systems （steady and nonsteady state） is constructed and analyzed by using linear triangle elements, and the O（h^2）-order optimal error estimates are derived under the assumption of strongly regular triangulation and the Ha-regularity for the exact solutions. The convergence analysis is based on some superclose estimates of the interpolation approximation. Finally, we discuss the Maxwell equations in a two-dimensional domain, and numerical experiments are given to validate the theoretical results.     

Determination of Unknown Boundary in the Composite Materials with Stefan-Boltzmann Conditions

Abstract The authors consider one specific kind of heat transfer problems in a threedimensional layered domain, with nonlinear Stefan-Boltzmann conditions on the boundaries as well as on the interfaces. To determine the unknown part of the boundary （or corrosion） by the Cauchy data on the reachable part is an important inverse problem in engineering. The mathematical model of this problem is introduced, the well-posedness of the forward problems and the uniqueness of the inverse problems are obtained.     

FINITE ELEMENT METHODS FOR SOBOLEV EQUATIONS

A new high-order time-stepping finite element method based upon the high-order numerical integration formula is formulated for Sobolev equations, whose computations consist of an iteration procedure coupled with a system of two elliptic equations. The optimal and superconvergence error estimates for this new method are derived both in space and in time. Also, a class of new error estimates of convergence and superconvergence for the time-continuous finite element method is demonstrated in which there are no time derivatives of the exact solution involved, such that these estimates can be bounded by the norms of the known data. Moreover, some useful a-posteriori error estimators are given on the basis of the superconvergence estimates.     

BOUNDARY ELEMENT ANALYSIS FOR A CLASS OF MULTI-DIMENSIONAL LINEAR PARABOLIC EQUATIONS

In this paper,by me as of beundary element method,we try to deal with the initial -boundary value problem for a class of linear parunolic equations,which is a linear heat conduction equation. We tresent a boundary integral equation for the solution to the problem and its variational formalation The well-posedness of the variational formulation is proved. And the error estimates for the approsutate solutions are provided. The results of this paper are more general than those of[1]     

SOURCE TERM IDENTIFICATION WITH DISCONTINUOUS DUAL RECIPROCITY APPROXIMATION AND QUASI-NEWTON METHOD FROM BOUNDARY OBSERVATIONS

This paper deals with discontinuous dual reciprocity boundary element method for solving an inverse source problem.The aim of this work is to determine the source term in elliptic equations for nonhomogenous anisotropic media,where some additional boundary measurements are required.An equivalent formulation to the primary inverse problem is established based on the minimization of a functional cost,where a regularization term is employed to eliminate the oscillations of the noisy data.Moreover,an efficient algorithm is presented and tested for some numerical examples.     

LOCAL GAUSSIAN-COLLOCATION SCHEME TO APPROXIMATE THE SOLUTION OF NONLINEAR FRACTIONAL DIFFERENTIAL EQUATIONS USING VOLTERRA INTEGRAL EQUATIONS

This work describes an accurate and effective method for numerically solving a class of nonlinear fractional differential equations.To start the method,we equivalently convert these types of differential equations to nonlinear fractional Volterra integral equations of the second kind by integrating from both sides of them.Afterward,the solution of the mentioned Volterra integral equations can be estimated using the collocation method based on locally supported Gaussian functions.The local Gaussian-collocation scheme estimates the unknown function utilizing a small set of data instead of all points in the solution domain,so the proposed method uses much less computer memory and volume computing in comparison with global cases.We apply the composite non-uniform Gauss-Legendre quadrature formula to estimate singular-fractional integrals in the method.Because of the fact that the proposed scheme requires no cell structures on the domain,it is a meshless method.Furthermore,we obtain the error analysis of the proposed method and demon-strate that the convergence rate of the approach is arbitrarily high.Illustrative examples clearly show the reliability and efficiency of the new technique and confirm the theoretical error estimates.     

ERROR ESTIMATES FOR THE RECURSIVE LINEARIZATION OF INVERSE MEDIUM PROBLEMS

This paper is devoted to the mathematical analysis of a general recursive linearization algorithm for solving inverse medium problems with multi-frequency measurements. Under some reasonable assumptions, it is shown that the algorithm is convergent with error estimates. The work is motivated by our effort to analyze recent significant numerical results for solving inverse medium problems. Based on the uncertainty principle, the recursive linearization allows the nonlinear inverse problems to be reduced to a set of linear problems and be solved recursively in a proper order according to the measurements. As an application, the convergence of the recursive linearization algorithm [Chen, Inverse Problems 13（1997）, pp.253-282] is established for solving the acoustic inverse scattering problem.     

A POSTERIORI ERROR ESTIMATES FOR LOCAL CO DISCONTINUOUS GALERKIN METHODS FOR KIRCHHOFF PLATE BENDING PROBLEMS

We derive some residual-type a posteriori error estimates for the local CO discontinuous Galerkin （LCDG） approximations （[31]） of the Kirchhoff bending plate clamped on the boundary. The estimator is both reliable and efficient with respect to the moment-field approximation error in an energy norm. Some numerical experiments are reported to demonstrate theoretical results.     

Iterative methods for a forward-backward heat equation in two-dimension

A finite difference method is introduced to solve the forward-backward heat equation in two space dimensions. In this procedure, the backward and forward difference scheme in two subdomains and a coarse-mesh second-order central difference scheme at the middle interface are used. Maximum norm error estimate for the procedure is derived. Then an iterative method based on domain decomposition is presented for the numerical scheme and the convergence of the given method is established. Then numerical experiments are presented to support the theoretical analysis.     

Discontinuous Galerkin finite element method for a nonlinear boundary value problem

In this paper, we investigate the a priori and a posteriori error estimates for the discontinuous Galerkin finite element approximation to a regularization version of the variational inequality of the second kind. We show the optimal error estimates in the DG-norm （stronger than the H1 norm） and the L2 norm, respectively. Furthermore, some residual-based a posteriori error estimators are established which provide global upper bounds and local lower bounds on the discretization error. These a posteriori analysis results can be applied to develop the adaptive DG methods.     

EXISTENCE OF GLOBAL SMOOTH SOLUTION TO INITIAL BOUNDARY VALUE PROBLEM FOR A VISCOELASTIC MODEL WITH RELAXATION

This paper is concerned with the initial boundary value problem for a vis-coelastic model with relaxation. Under the only assumption that the C^0-norm of theinitial data is small, without smallness hypothesis for the C^1-norm, the existence of theglobal smooth solution to the corresponding initial boundary value problem is proved.The analysis is based on some a priori estimates obtained by the “maximum principle” offirst-order quasilinear hyperbolic system.