An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations

The present paper deals with the numerical solution of the coupled Schrdinger-KdV equations using the elementfree Galerkin(EFG) method which is based on the moving least-square approximation.Instead of traditional mesh oriented methods such as the finite difference method(FDM) and the finite element method(FEM),this method needs only scattered nodes in the domain.For this scheme,a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method.In numerical experiments,the results are presented and compared with the findings of the finite element method,the radial basis functions method,and an analytical solution to confirm the good accuracy of the presented scheme.     

方位电阻率成像仪（ARI）的三维有限元模拟

为了研究方位电阻率成像仪的定量解释方法，开发了一套能模拟方位电阻率成像仪在三维非均匀地层模拟中响应的三维有限元程序，在简化的地层模型中，将该程序与解析方法和二维有限差分方法进行了对比，结果表明此程序的误差不超过3％，在此程序得到验证后，用该程序考查了方位电阻率成在一些典型三维非均匀地层模型中的响应，结果表明此程序能为方位电阻率成像仪的定量解释提供支持。     

微波元件与加速结构的有限元分析

研究了有限元方法在二维均匀结构和轴对称加速结构中的应用,采用了带宽优化技术和子空间求解特征值方法,并给出了部分例子,结果表明该方法精度高、速度快.这些方法可直接用于三维程序、并行计算.     

平面应变问题广义有限元法

结合广义有限元和理性有限元,针对平面应变问题提出一种广义四边形单元.该单元考虑泊松效应,以节点位移自由度约束弹性力学平面应变方程的半解析解,构造单元位移模式的附加项,较准确地反映真实位移场,提高了单元的计算精度.首先推导广义单元及其等参单元的形函数,之后设计分片试验和数值算例来验证单元的精度,数值算例表明,在规则和非规则网格下新单元计算精度均优于传统有限元和广义有限元.新单元精度高且易于程序实现,可应用于实际工程的结构分析.     

分数阶对流扩散方程的特征有限元方法

讨论非线性分数阶对流扩散方程的特征有限元方法.利用特征线法和分数阶有限元框架,构建一种基于特征方向的全离散有限元格式.模拟物理问题,并在数值上与常规有限元格式进行比较,计算结果表明：该方法能准确地捕捉到控制方程的精确解,即使是在对流效应占优时,也具有稳定性好和逼近精度高等特征.     

超声波检测计算机模拟和仿真的研究及应用现状

随着高性能计算机的广泛应用，超声检测模拟和仿真的研究方兴未艾。本文对这些研究所选择的数学模型和方法进行了综合评述，探讨了Gaussiam模型、声线模型、衍射理论、有限积分技术等的应用，涉及的方法包括各种解析方法及有限差分、有限元、边界元等数值方法。     

二维多流体域耦合声场的光滑有限元解法

针对声学有限元分析中四节点等参单元计算精度低,对网格质量敏感的问题,将光滑有限元法引入到多流体域耦合声场的数值分析中,提出了二维多流体域耦合声场的光滑有限元解法。该方法在Helmholtz控制方程与多流体域耦合界面的声压/质点法向速度连续条件的基础上,得到二维多流体耦合声场的离散控制方程,并采用光滑有限元的分区光滑技术将声学梯度矩阵形函数导数的域内积分转换形函数的域边界积分,避免了雅克比矩阵的计算。以管道二维多流体域耦合内声场为数值分析算例,研究结果表明,与标准有限元相比,对单元尺寸较大或扭曲严重的四边形网格模型,光滑有限元的计算精度更高。因此光滑有限元能很好地应用于大尺寸单元或扭曲严重的网格模型下二维多流体域耦合声场的预测,具有良好的工程应用前景。     

响应面法的杆式超声电机有限元模型修正

现有杆式超声电机的有限元模型大都采用连续复合材料的简化结构形式,忽略了实际超声电机各部件间的接触特性,导致计算结果与实验结果具有较大偏差。针对该问题,提出了一种采用响应面法的杆式超声电机有限元模型修正方法,来获得高精度的有限元模型。该方法考虑了超声电机螺栓预紧力及各部件接触界面的法向接触刚度和摩擦系数,筛选出显著影响电机工作模态的参数,建立响应面模型替代超声电机的有限元模型实现快速计算结构响应的目的,并以实验模态分析结果为目标对模型进行修正。修正结果表明,修正后的模型模态频率的平均误差由修正前的1.20%降到0.21%,模型精度得到明显改善,表明以响应面的有限元模型修正方法对杆式超声电机以及类似夹心式压电振子的设计具有应用价值。     

光滑粒子流体动力学有限元法接触算法研究

 耦合光滑粒子流体动力学（SPH）方法和有限元法（FEM）,是研究冲击动力学问题的一种有效途径。为解决SPH粒子和有限单元间的接触问题,提出了SPH-FEM接触算法。该算法是在有限元节点处设置背景粒子,采用SPH接触算法的思想,计算施加在SPH粒子和有限元节点上的接触力,并且以外力的形式分别加入到SPH动量方程和有限元动力学方程中。利用SPH-FEM接触算法,对两杆撞击以及圆柱形钢弹正冲击钢板发生的冲塞破坏过程进行了三维数值模拟,靶板采用含损伤的Johnson-Cook模型和Grüneisen状态方程,模拟结果与实验结果吻合较好。     

A coupled edge-/face-based smoothed finite element method for structural-acoustic problems

The edge-based smoothed finite element method (ES-FEM) and the face-based smoothed finite element method (FS-FEM) developed recently have shown great efficiency in solving solid mechanics problems with triangular and tetrahedral meshes. In this paper, a coupled ES-/FS-FEM model is extended to solve the structural-acoustic problems consisting of a plate structure interacting with the fluid medium. Three-node triangular elements and four-node tetrahedral elements are used to discretize the two-dimensional (2D) plate and three-dimensional (3D) fluid, respectively, as they can be generated easily and even automatically for complicated geometries. The field variable in each element is approximated using the linear shape functions, which is exactly the same as that in the standard FEM. The gradient field of the problem is obtained particularly using the gradient smoothing operation over the edge-based and face-based smoothing domains in 2D and 3D, respectively. The gradient smoothing technique can provide a proper softening effect to the model, effectively solve the problems caused by the well-known “overly-stiff” phenomenon existing in the standard FEM, and hence significantly improve the accuracy of the solution for the coupled systems. Intensive numerical studies have been conducted to verify the effectiveness of the coupled ES-/FS-FEM for structural-acoustic problems.     

Finite element model updating of a rod-type ultrasonicmotor based on response surface method

The conventional finite element model(FEM) of a rod-type ultrasonic motor is usually simplified by means of continuous composite structure. Because the actual contact characteristics between the parts of the ultrasonic motor is ignored, there is bigger error between the calculated values and experimental results. Aiming at solving problem, a new modeling method of a rod-type ultrasonic motor is presented to obtain a high-accuracy FEM. The bolt pretension and the normal contact stiffness and friction coefficient of the contact surface of ultrasonic motor are all considered in this method, and the significant parameters of working mode of the motor are selected by the response surface method, and the goal of calculating the structural response rapidly is realized by building the response surface model to replace the FEM. The result of finite element model updating shows that the average error of modal frequencies of updated model drops to 0.21% from 1.20%. The accuracy of FEM is obviously improved, which indicates that the FEM updating based on response surface method is of great application value on the design for a rod-type ultrasonic motor.     

Free Mesh Method: fundamental conception, algorithms and accuracy study

The finite element method (FEM) has been commonly employed in a variety of fields as a computer simulation method to solve such problems as solid, fluid, electro-magnetic phenomena and so on. However, creation of a quality mesh for the problem domain is a prerequisite when using FEM, which becomes a major part of the cost of a simulation. It is natural that the concept of meshless method has evolved. The free mesh method (FMM) is among the typical meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, especially on parallel processors. FMM is an efficient node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm for the finite element calculations. In this paper, FMM and its variation are reviewed focusing on their fundamental conception, algorithms and accuracy.     

A finite-element method using dispersion reduced spline elements for room acoustics simulation

This paper presents a finite element method (FEM) using hexahedral 27-node spline acoustic elements (Spl27) with low numerical dispersion for room acoustics simulation in both the frequency and time domains, especially at higher frequencies. Dispersion error analysis in one dimension is performed to increase the accuracy of FEM using Spl27 by modifying the numerical integration points of element stiffness and mass matrices. The basic accuracy and efficiency of the FEM using the improved Spl27, which uses modified integration points, are presented through numerical experiments using benchmark problems in both the frequency and time domains, revealing that FEM using the improved Spl27 in both domains provides more accurate results than the conventional method does, and with fewer degrees of freedom. Moreover, the effectiveness of FEM using the improved Spl27 over that using hexahedral 27-node Lagrange elements is shown for time domain analysis of the sound field in a practical sized room.     

标准有限元方法模拟大气压流注等离子体中数值振荡与伪扩散问题研究与改进

基于商用软件的等离子体数值模拟方案大多是基于标准Galerkin格式的有限元方法,在计算大气压流注放电时,存在数值振荡和伪扩散的问题,计算结果可靠性不高.针对针-板体流注放电结构,构建了修正的COMSOL Multiphysics等离子体模型.修正模型与COMSOL Multiphysics内置等离子体模块的区别主要在于增加了基于弱形式的人工稳定项.利用修正模型复现了经典的雪崩-流注转捩算例,并与经典纳秒秒冲过电压针板放电实验结果进行了对比验证.结果表明,引入修正后的商用软件模型能够克服数值振荡和伪扩散缺陷,确保等离子体数值计算的准确性和可复现性.      

Finite element approximation in quantum field theory

Lattices are conventionally used to regulate quantum field theory. However, lattices have a number of drawbacks, one of the most noteworthy being the awkwardness which arises in treating the first difference operator associated with the fermion kinetic energy. We propose here an alternative to the conventional lattice formulation of quantum field theories. This alternative is based on the method of finite elements, which has been a powerful tool in continuum mechanics and fluid dynamics for many years. We believe that this method may be advantageous in quantum field theory because it treats bosons and fermions on a nearly equal footing. This is possible because in the finite element method one does not introduce a finite difference approximation for derivatives. In this paper, we adapt the finite element technique to quantum field theory, and apply it to several elementary field theory models.     

基于大规模并行计算的三维多群中子扩散方程有限差分方法

三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点;然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法:采用中心有限差分格式离散中子扩散方程;基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算;采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。     

Morphological evolution and migration of void in bi-piezoelectric interface based on nonlocal phase field method

This paper reports the result of investigation into the morphological evolution and migration of void in bi-piezoelectric material interface by utilizing nonlocal phase field model and finite element method (FEM), where the small scale effect containing the long-range forces among atoms is considered. The nonlocal elastic strain energy and the nonlocal electric energy around the void are firstly calculated by the finite element method. Then based on the finite difference method (FDM), the thermodynamic equilibrium equation containing the surface energy and anisotropic diffusivity is solved to simulate the morphological evolution and migration of elliptical void in bi-piezoelectric films interface. Results show that the way of load condition plays a significant role in the evolution process, and the boundary of void's long axis gradually collapses toward the center of ellipse. In addition, the evolutionary speed of left boundary gradually decreases with scale effect coefficient growth. This work can provide references for the safety evaluation of piezoelectric materials in micro electro mechanical system.     

基于大规模并行计算的三维多群中子扩散方程有限差分方法

三维多群中子扩散方程的精确、高效求解是核动力堆芯设计及燃料管理的基础。应用有限差分方法求解该方程具有简便、精确、成熟的优点;然而,该方法的计算量和存储量均较大,极大地限制了它的计算规模和应用范围。本文基于大规模并行计算,研究三维多群中子扩散方程有限差分方法：采用中心有限差分格式离散中子扩散方程;基于MPI并行编程模型,采用空间区域分解的方式实现大规模并行计算;采用多群多区域耦合PGMRES算法进行并行加速。在集群服务器上开发了ParaFiDi程序,并采用IAEA3D,PHWR等多个基准题对该程序进行验证。数值结果表明,ParaFiDi程序具有较高的计算精度和计算效率。     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

STOCHASTIC PERTURBATION APPROACH TO ENGINEERING STRUCTURE VIBRATIONS BY THE FINITE DIFFERENCE METHOD

The main idea of the paper is to introduce the second order perturbation second probabilistic moment analysis in the context of the finite difference method (FDM) modelling of vibrations. The approach can be successfully applied in all those engineering analyses where FDM modelling of engineering structures vibrations is still useful and, at the same time, some structural parameters are random variables or fields. The general advantage of the stochastic finite difference method (SFDM) proposed is the relatively easy extension of the existing deterministic results of the classical elastodynamics on the random or stochastic case. However, similarly to stochastic boundary or finite element methods, the approach proposed has its limitations on the second order random uncertainties measures of input random variables.