A Minimal Distance Constrained Adaptive Mesh Algorithm with Application to the Dual Reciprocity Method

We present a new technique for generating error equi-distributing meshes that satisfy both local quasi-uniformity and a preset minimal mesh spacing. This is first done in the one-dimensional case by extending the Kautsky and Nichols method and then in the two-dimensional case by generalizing the tensor product methods to alternating curved line equi-distributions. With the new meshing approach, we have achieved better accuracy in approximation using interpolatory radial basis functions. Furthermore, improved accuracy in numerical results has been obtained when the interpolatory strategy is applied to the dual reciprocity boundary element method for solving a class of linear and nonhomogeneous partial differential equations.     

Fourier analysis on trapezoids with curved sides

It is well known that smooth periodic functions can be expanded into Fourier series and can be approximated by trigonometric polynomials. The purpose of this paper is to do Fourier analysis for smooth functions on planar domains. A planar domain can often be divided into some trapezoids with curved sides, so first we do the Fourier analysis for smooth functions on trapezoids with curved sides. We will show that any smooth function on a trapezoid with curved sides can be expanded into Fourier sine series with simple polynomial factors, and so it can be well approximated by a combination of sine polynomials and simple polynomials. Then we consider the Fourier analysis on the global domain. Finally, we extend these results to the three-dimensional case.     

Approximation by Piecewise Constant Functions on a Square

In this paper we consider several algorithms for approximating functions defined on the unit square I = [0,1]² and ranging in $\mathbb{R}^2 $ . We use functions of zeroth-order Lagrange spline type as the approximation apparatus. They differ from the standard Lagrange splines on the plane by the rule for choosing grid lines according to which the spline is constructed; namely, a set of one-dimensional splines is used instead of a family of parallel lines determining the interpolation nodes.     

A film coating on a rough surface of an elastic body

A solution of the plane problem of the theory of elasticity for a film–substrate composite is solved by a perturbation method for a substrate with a rough surface. An algorithm for calculating any approximation, which ultimately leads to the solution of the same Fredholm equation of the second kind, is given. Formulae for calculating the right-hand side of this equation, which depends on all the preceding approximations, are derived. An exact solution of the integral equation in the form of Fourier series, whose coefficients are expressed in quadratures, is given in the case of a substrate with a periodically curved surface. The stresses on the flat surface of the film and on the interfacial surface are found in a first approximation as functions of the form of bending of the surface, the mean thickness of the film and the ratio of Young's moduli of the film and the substrate. It is shown, in particular, that the greatest stress concentration on the film surface occures on a protrusion of the softer substrate. ©2013     

The Use of Parabolic Arcs in Matching Curved Boundaries in the Finite Element Method

If isoparametric coordinates are used to deal with curved boundariesin the finite element method, the original boundary is implicitlyreplaced by a series of parabolic or cubic arcs. The equationsof these arcs involve parameters which are the coordinates ofpoints on the curved side, and a simple procedure is outlinedfor choosing these parameters in such a way that each arc isa parabola which passes through four points of the originalcurve thus ensuring a good approximation to it.     

Overcoming the boundary effects in surface spline interpolation

Surface spline interpolation when the domain is all of R^d isknown to converge much faster to the data function f than inthe case when the domain is the unit ball. This difference isunderstood to be due to boundary effects which, as will be shown,also affect the size of the surface spline's coefficients. Wepropose a modified form of surface spline interpolation which,to a great extent, overcomes these boundary effects. This modifiedsurface spline interpolant uses only the values of f at thegiven interpolation points.     

Fast deconvolution with approximated PSF by RSTLS with antireflective boundary conditions

The problem of reconstructing signals and images from degraded ones is considered in this paper. The latter problem is formulated as a linear system whose coefficient matrix models the unknown point spread function and the right hand side represents the observed image. Moreover, the coefficient matrix is very ill-conditioned, requiring an additional regularization term. Different boundary conditions can be proposed. In this paper antireflective boundary conditions are considered. Since both sides of the linear system have uncertainties and the coefficient matrix is highly structured, the Regularized Structured Total Least Squares approach seems to be the more appropriate one to compute an approximation of the true signal/image. With the latter approach the original problem is formulated as an highly nonconvex one, and seldom can the global minimum be computed. It is shown that Regularized Structured Total Least Squares problems for antireflective boundary conditions can be decomposed into single variable subproblems by a discrete sine transform. Such subproblems are then transformed into one-dimensional unimodal real-valued minimization problems which can be solved globally. Some numerical examples show the effectiveness of the proposed approach.     

A strip cut in a transversely isotropic elastic solid

The three-dimensional problems of a strip cut in a transversely isotropic elastic space, when the isotropy planes are perpendicular to the plane of the cut, are investigated using the asymptotic methods developed by Aleksandrov and his coauthors. Two cases of the location of the strip cut are considered: along the first axis of a Cartesian system of coordinates (Problem A) or along the second axis (Problem B). Assuming that the normal load, applied to the sides of the cut (normal separation friction) can be represented by a Fourier series, one-dimensional integral equations of problems A and B are obtained, the symbols of the kernels of which are independent of the number of the term of the Fourier series. A closed solution of the problem is derived for a special approximation of the kernel symbol. Regular and singular asymptotic methods are also used to solve the integral equations by introducing a dimensionless geometrical parameter, representing the ratio of the period of the applied wavy normal load to the thickness of the cut strip. The normal stress intensity factor on the strip boundary is calculated using the three methods of solving the integral equations indicated.     

Cubic spline for a class of singular two-point boundary value problems

In this paper we have presented a method based on cubic splines for solving a class of singular two-point boundary value problems. The original differential equation is modified at the singular point then the boundary value problem is treated by using cubic spline approximation. The tridiagonal system resulting from the spline approximation is efficiently solved by Thomas algorithm. Some model problems are solved, and the numerical results are compared with exact solution.     

Numerical calculation of fourier integrals with cubic splines

The most used formula for calculation of Fourier integrals is Filon's formula which is based on approximation of the function by a quadratic in each double interval. In order to obtain a better approximation we use the cubic spline fit. The method is not restricted to equidistant points, but the final formulas are only derived in this case. Test computations show that the spline formula may be superior to Filon's formula.     

Interpolation and approximation by monotone cubic splines

We study the reconstruction of a function defined on the real line from given, possibly noisy, data values and given shape constraints. Based on two abstract minimization problems characterization results are given for interpolation and approximation (in the euclidean norm) under monotonicity constraints. We derive from these results Newton-type algorithms for the computation of the monotone spline approximant.     

OnG 2 continuous cubic spline interpolation

In this paper the interpolation byG ² continuous planar cubic Bézier spline curves is studied. The interpolation is based upon the underlying curve points and the end tangent directions only, and could be viewed as an extension of the cubic spline interpolation to the curve case. Two boundary, and two interior points are interpolated per each spline section. It is shown that under certain conditions the interpolation problem is asymptotically solvable, and for a smooth curvef the optimal approximation order is achieved. The practical experiments demonstrate the interpolation to be very satisfactory. Supported in prat by the Ministry of Science and Technology of Slovenjia, and in part by the NSF and SF of National Educational Committee of China.     

THE GLOBAL ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE 3D FLOW AROUND A SUBMERGED BODY

We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the up-stream artificial boundary and two vertical sides are introduced as the downstream arti-ficial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.     

W_2~m空间中样条插值算子与线性泛函的最佳逼近

In this paper,the convergency of spline interpolation operators is obtained,these spline operators are determined by linear differential operators and constraint functionals.The errors of the interpolating spline with EHB fanctionals are estimated.The best approximation of linear functionals on W2^n spaces are investigated,which let to a useful computational method for the approximation solution of higher order linear differential equations with multipoint boundary value conditions.     

A uniformly convergent hybrid scheme for singularly perturbed system of reaction-diffusion Robin type boundary-value problems

This paper deals with the study on system of reaction diffusion differential equations for Robin or mixed type boundary value problems (MBVPs). A cubic spline approximation has been used to obtain the difference scheme for the system of MBVPs, on a piecewise uniform Shishkin mesh defined in the whole domain. It has been shown that our proposed scheme, i.e., central difference approximation for outer region with cubic spline approximation for inner region of boundary layers, leads to almost second order parameter uniform convergence whereas the standard method i.e., the forward-backward approximation for mixed boundary conditions with central difference approximation inside the domain leads to almost first order convergence on Shishkin mesh. Numerical results are provided to show the efficiency and accuracy of these methods.     

Approximate grid solution of a nonlocal boundary value problem for Laplace’s equation on a rectangle

A nonlocal boundary value problem for Laplace’s equation on a rectangle is considered. Dirichlet boundary conditions are set on three sides of the rectangle, while the boundary values on the fourth side are sought using the condition that they are equal to the trace of the solution on the parallel midline of the rectangle. A simple proof of the existence and uniqueness of a solution to this problem is given. Assuming that the boundary values given on three sides have a second derivative satisfying a Hölder condition, a finite difference method is proposed that produces a uniform approximation (on a square mesh) of the solution to the problem with second order accuracy in space. The method can be used to find an approximate solution of a similar nonlocal boundary value problem for Poisson’s equation.     

长波在弯曲管道中的传播——Ⅰ变截面弯管中长波传播的基本分析

本文研究长波在三维变截面弯管中的传播问题.通过建立正交曲线坐标系,以波数k和管道横截面的特征半径a的乘积ka作为小参数,对波动方程进行无量纲处理,用正则摄动法,把三维的Helmholtz方程化为二维的Laplace(或Poisson)方程和一维的Webster方程.并分析了管道的几何参数(横截面面积、管道中心线的曲率和挠度)对复速度势渐近展开的各阶项的影响.文中指出,横截面面积的变化首先影响浙近解的零阶项.在横截面的形状具有某种对称性时,管道中心线的曲率首先影响渐近解的二阶项,而挠度首先影响渐近解的三阶项.最后,给出了长波在弯曲圆管中传播的实例.     

Integral representation of slowly growing equidistant splines

In this paper we consider polynomial splines S(x) with equidistant nodes which may grow as O (|x|^s). We present an integral representation of such splines with a distribution kernel. This representation is related to the Fourier integral of slowly growing functions. The part of the Fourier exponentials herewith play the so called exponential splines by Schoenberg. The integral representation provides a flexible tool for dealing with the growing equidistant splines. First, it allows us to construct a rich library of splines possessing the property that translations of any such spline form a basis of corresponding spline space. It is shown that any such spline is associated with a dual spline whose translations form a biorthogonal basis. As examples we present solutions of the problems of projection of a growing function onto spline spaces and of spline interpolation of a growing function. We derive formulas for approximate evaluation of splines projecting a function onto the spline space and establish therewith exact estimations of the approximation errors.     

Time-dependent laser ignition of a combustible stagnant boundary layer

A spatially one-dimensional model for the ignition of a combustiblegas layer adjacent to a plane solid surface is considered. Theeffect of an incident laser beam from the gas side is to raisethe surface temperature; the other boundary is taken to be apoor thermal conductor in the form of an inert gaseous medium.It is assumed that the exothermic chemical reaction within thelayer has a negligible reactant consumption. Three examplesfor the effect of the laser on the solid surface are considered:(a) a large instantaneous temperature rise, (b) a temperaturejump at t=0 which is a linear funtion of time t>0, (c) asurface temperature variation which is a linear function oftime. For (a) and (b) conditions for criticality are obtained.For appropriate states of the system, times to ignition for(a), (b), and (c) have been determined. It is shown that thetheoretical results are in reasonably good agreement with experiment.The kind of physical situation envisaged might occur in a mineshaftwhen a stagnant combustible layer is subject to intense transientlight. Although an idealization, our analysis determines conditionsunder which thermal runaway and subsequent explosion could takeplace.     

Quadratic spline interpolation with coinciding interpolation and spline grids

In this paper the quadratic spline interpolation with coinciding interpolation and spline grids for continuous functions is considered. The theorems mainly concern error estimations which allow to formulate a convergence statement. To get such results it is assumed that the function to be interpolated is suitably smooth or possesses a special behavior. A best approximation property and a statement about the solution of boundary value problems using quadratic spline functions are added.