High order numerical approximation of minimal surfaces

We present an algorithm for finding high order numerical approximations of minimal surfaces with a fixed boundary. The algorithm employs parametrization by high order polynomials and a linearization of the weak formulation of the Laplace–Beltrami operator to arrive at an iterative procedure to evolve from a given initial surface to the final minimal surface. For the steady state solution we measure the approximation error in a few cases where the exact solution is known. In the framework of parametric interpolation, the choice of interpolation points (mesh nodes) is directly affecting the approximation error, and we discuss how to best update the mesh on the evolutionary surface such that the parametrization remains smooth. In our test cases we may achieve exponential convergence in the approximation of the minimal surface as the polynomial degree increases, but the rate of convergence greatly differs with different choices of mesh update algorithms. The present work is also of relevance to high order numerical approximation of fluid flow problems involving free surfaces.     

A meshless Galerkin method with moving least square approximations for infinite elastic solids

Combining moving least square approximations and boundary integral equations, a meshless Galerkin method, which is the Galerkin boundary node method （GBNM）, for twoand three-dimensional infinite elastic solid mechanics problems with traction boundary conditions is discussed. In this numerical method, the resulting formulation inherits the symmetry and positive definiteness of variational problems, and boundary conditions can be applied directly and easily. A rigorous error analysis and convergence study for both displacement and stress is presented in Sobolev spaces. The capability of this method is illustrated and assessed by some numerical examples.     

Coupling of finite element and boundary integral methods for electromagnetic scattering in a two-layered medium

Consider a time-harmonic electromagnetic plane wave incident on an inhomogeneity embedded in a two-layered medium. In this paper, a method of coupling of finite element and boundary integral equation methods is presented for the solutions of electromagnetic scattering in both transverse electric and magnetic polarization cases. The well-posedness of the continuous and discrete problems, as well as optimal error estimates for the coupled variational approximations, are obtained. Numerical results are included to illustrate the accuracy with the optimal convergence property of the proposed method and to show the wave features in a two-layered medium.     

一类各项异性半线性椭圆方程自然边界元与有限元耦合法

将冯康和余德浩提出的自然边界归化方法用于研究一类半线性椭圆方程外区域问题，提出一种自然边界元与有限元的耦合算法、针对某一类半线性椭圆方程，应用变分原理，研究其弱解性及Galerkin逼近，得到有限元解的误差估计及收敛阶O(h^n)，最后给出相应数值例子。     

Solving unsteady Schr?dinger equation using the improved element-free Galerkin method

By employing the improved moving least-square (IMLS) approximation, the improved element-free Galerkin (IEFG) method is presented for the unsteady Schrödinger equation. In the IEFG method, the two-dimensional (2D) trial function is approximated by the IMLS approximation, the variation method is used to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. Because the number of coefficients in the IMLS approximation is less than in the moving least-square (MLS) approximation, fewer nodes are needed in the entire domain when the IMLS approximation is used than when the MLS approximation is adopted. Then the IEFG method has high computational efficiency and accuracy. Several numerical examples are given to verify the accuracy and efficiency of the IEFG method in this paper.     

解弹性力学第二类边界积分方程的求积法与分裂外推

利用Sidi奇异求积公式,提出了解曲边多角形域上线性弹性力学第二类边界积分方程的求积法,即离散矩阵的每个元素的生成只需赋值不需计算任何奇异积分.通过估计离散矩阵的特征值和利用Anselone聚紧收敛理论,证明了近似解的收敛性;同时得到了误差的多参数渐近展开式;通过并行地解粗网格上的离散方程,利用分裂外推获得了高精度近似解和后验误差.     

Theoretical and Mathematical Physics

A new effective method for solving three-dimensional problems of electromagnetic wave diffraction by impedance bodies with irregular boundaries is proposed. The method offers a high rate of convergence. Examples of solving the problems of wave scattering by bodies of revolution are given, and results illustrating the rate of convergence of a computational algorithm for bodies of various shapes are presented. The impedance approximation is shown to be valid for simulation of scattering characteristics of bodies with an insulating coating even when the boundary has irregularities and the refractive index of the coating is not too high. Various ways of characterizing “black bodies” and the results of studying their scattering characteristics are discussed.     

基于对偶混合变分原理的Signorini问题的数值模拟

基于Signorini问题的对偶混合变分形式,提出了一种非协调有限元逼近格式,证明了离散的B-B条件,获得了Raviart-Thomas(k=0)有限元逼近的误差界O(h^3/4),并且Uzawa型算法对协调与非协调有限元逼近格式进行了数值求解.根据数值结果的分析和比较,表明应用非协调有限元逼近格式求解更有效.     

Covolume-upwind finite volume approximations for linear elliptic partial differential equations

In this paper, we study covolume-upwind finite volume methods on rectangular meshes for solving linear elliptic partial differential equations with mixed boundary conditions. To avoid non-physical numerical oscillations for convection-dominated problems, nonstandard control volumes (covolumes) are generated based on local Peclet’s numbers and the upwind principle for finite volume approximations. Two types of discretization schemes with mass lumping are developed with use of bilinear or biquadratic basis functions as the trial space respectively. Some stability analyses of the schemes are presented for the model problem with constant coefficients. Various examples are also carried out to numerically demonstrate stability and optimal convergence of the proposed methods.     

An accurate, stable and efficient domain-type meshless method for the solution of MHD flow problems

The aim of the present paper is the development of an efficient numerical algorithm for the solution of magnetohydrodynamics flow problems for regular and irregular geometries subject to Dirichlet, Neumann and Robin boundary conditions. Toward this, the meshless point collocation method (MPCM) is used for MHD flow problems in channels with fully insulating or partially insulating and partially conducting walls, having rectangular, circular, elliptical or even arbitrary cross sections. MPC is a truly meshless and computationally efficient method. The maximum principle for the discrete harmonic operator in the meshfree point collocation method has been proven very recently, and the convergence proof for the numerical solution of the Poisson problem with Dirichlet boundary conditions have been attained also. Additionally, in the present work convergence is attained for Neumann and Robin boundary conditions, accordingly. The shape functions are constructed using the Moving Least Squares (MLS) approximation. The refinement procedure with meshless methods is obtained with an easily handled and fully automated manner. We present results for Hartmann number up to 10⁵${10}^5$. The numerical evidences of the proposed meshless method demonstrate the accuracy of the solutions after comparing with the exact solution and the conventional FEM and BEM, for the Dirichlet, Neumann and Robin boundary conditions of interior problems with simple or complex boundaries.     

铁热稠密物质状态方程的相对论Hartree-Fock-Slater计算

 本文在球元胞、中心力场、相对论Fermi统计近似下，用平均原子模型计算铁热稠密物质在任意温度和密度下的单电子能级和状态方程。这种算法以Zink的解析拟合势作为初始势，求解满足Wigner-Seitz边界条件的Dirac-Slater径向波函数方程。通过简并度是密度的连续函数，考虑了大密度的压致电离效应。这样从另一个角度考虑了电子的能带结构。     

An improved boundary element-free method （IBEFM） for two-dimensional potential problems

The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker δ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker δ function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.     

Elastic-wave identification of penetrable obstacles using shape-material sensitivity framework

This study deals with elastic-wave identification of discrete heterogeneities (inclusions) in an otherwise homogeneous “reference” solid from limited-aperture waveform measurements taken on its surface. On adopting the boundary integral equation (BIE) framework for elastodynamic scattering, the inverse query is cast as a minimization problem involving experimental observations and their simulations for a trial inclusion that is defined through its boundary, elastic moduli, and mass density. For an optimal performance of the gradient-based search methods suited to solve the problem, explicit expressions for the shape (i.e. boundary) and material sensitivities of the misfit functional are obtained via the adjoint field approach and direct differentiation of the governing BIEs. Making use of the message-passing interface, the proposed sensitivity formulas are implemented in a data-parallel code and integrated into a nonlinear optimization framework based on the direct BIE method and an augmented Lagrangian whose inequality constraints are employed to avoid solving forward scattering problems for physically inadmissible (or overly distorted) trial inclusion configurations. Numerical results for the reconstruction of an ellipsoidal defect in a semi-infinite solid show the effectiveness of the proposed shape-material sensitivity formulation, which constitutes an essential computational component of the defect identification algorithm.     

改进的多重网格法重建含遮拦物的干涉波前

阐述了由单张干涉图求取相位主值图像及最小二乘相位恢复。推导出基于梯度拟合的泊松方程及其离散化形成后,详细讨论了多重网格法对此方程求解的原理和过程。结合路径无关与路径相关算法的优点,设计了自适应最优路径法作为多重网格法的预处理,改进迭代的初始条件,提高了重建精度,同时大大加快了收敛速度。实验结果证明改进的多重网格法对含遮挡物的干涉图能获得很好的波前重建效果。     

Development of arbitrary-order multioperators-based schemes for parallel calculations.: Part 2: Families of compact approximations with two-diagonal inversions and related multioperators

One-parameter families of compact approximations to grid functionals with inverses of two-point operators and their properties are described. As particular examples, interpolation/extrapolations operators, quadratures formulas and approximations to derivatives are presented. Using operators from the families with fixed parameters values as basis operators, their linear combinations providing formally arbitrary-order approximations (multioperators) are constructed. Numerical illustrations are presented. Special emphasis is placed on first derivatives discretizations in the context of conservation laws. As an example, a highly accurate tenth-order scheme is outlined and tested against the Burgers’ equation. It is shown how extrapolation multioperators can be used to create boundary closures.     

Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method

The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503–1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710–1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.     

基于高精度自适应hp-FEM的随钻电阻率测井电场数值模拟

采用高精度自适应hp有限元(FEM)算法对随钻电阻率测井电场分布情况进行模拟,能够根据地层模型的实际情况和误差指示自动选择合适的细化方式和计算策略,不需要对整个求解域进行计算就能使有用信号呈指数速率准确地收敛到其真值附近.计算结果显示,获得的接收线圈感应电动势的相位差和幅度比曲线和地层分层情况吻合较好,有助于实际随钻电阻率测井资料的解释.     

A new complex variable meshless method for transient heat conduction problems

<正>In this paper,based on the improved complex variable moving least-square(ICVMLS) approximation,a new complex variable meshless method(CVMM) for two-dimensional(2D) transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations,and the essential boundary conditions are imposed by the penalty method.As the transient heat conduction problems are related to time,the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization.Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained.In order to demonstrate the applicability of the proposed method,numerical examples are given to show the high convergence rate,good accuracy,and high efficiency of the CVMM presented in this paper.     

The complex variable meshless local Petrov—Galerkin method of solving two-dimensional potential problems

Based on the complex variable moving least-square(CVMLS) approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin(CVMLPG) method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square(MLS) approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin(MLPG) method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.     

Some approximations in the linear dynamic equations of thin cylinders

Theoretical analysis is performed on the linear dynamic equations of thin cylindrical shells to find the error committed by making the Donnell assumption and the neglect of in-plane inertia. At first, the effect of these approximations is studied on a shell with classical simply supported boundary condition. The same approximations are then investigated for other boundary conditions from a consistent approximate solution of the eigenvalue problem. The Donnell assumption is valed at frequencies high compared with the ring frequencies, for finite length thin shells. The error in the eigenfrequencies from omitting tangential inertia is appreciable for modes with large circumferential and axial wave lengths, independent of shell thickness and boundary conditions.