PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

Partition of Unity Finite Element Analysis of Nonlinear Transient Diffusion Problems Using p-Version Refinement

We propose a high-order enriched partition of unity finite element method for linear and nonlinear time-dependent diffusion problems.The solution of this class of problems often exhibits non-smooth features such as steep gradients and boundary layers which can be very challenging to recover using the conventional low-order finite element methods.A class of steady-state exponential functions has been widely used for enrichment and its performance to numerically solve these challenges has been demonstrated.However,these enrichment functions have been used only in context of the standard h-version refinement or the so-called q-version refinement.In this paper we demonstrate that the p-version refinement can also be a very attractive option in terms of the efficiency and the accuracy in the enriched partition of unity finite element method.First,the transient diffusion problem is integrated in time using a semi-implicit scheme and the semi-discrete problem is then integrated in space using the p-version enriched finite elements.Numerical results are presented for three test examples of timedependent diffusion problems in both homogeneous and heterogeneous media.The computed results show the significant improvement when using the p-version refined enriched approximations in the finite element analysis.In addition,these results support our expectations for a robust and high-order accurate enriched partition of unity finite element method.

GALERKIN MESHLESS METHODS BASED ON PARTITION OF UNITY QUADRATURE

Numerical quadrature is an important ingredient of Galerkin meshless methods. A new numerical quadrature technique, partition of unity quadrature （PUQ）,for Galerkin meshless methods was presented. The technique is based on finite covering and partition of unity. There is no need to decompose the physical domain into small cell. It possesses remarkable integration accuracy. Using Element-free Galerkin methods as example, Galerkin meshless methods based on PUQ were studied in detail. Meshing is always not required in the procedure of constitution of approximate function or numerical quadrature, so Galerkin meshless methods based on PUQ are “truly” meshless methods.

Partitioned Method of Insect Flapping Flight for Maneuvering Analysis

This study proposed a partitioned method to analyze maneuvering of insects during flapping flight.This method decomposed the insect flapping flight into wing and body subsystems and then coupled them via boundary conditions imposed on the wing’s base using one-way coupling.In the wing subsystem,the strong coupling of the flexible wings and surrounding fluid was accurately analyzed using the finite element method to obtain the thrust forces acting on the insect’s body.The resulting thrust forces were passed from the wing subsystem to the body subsystem,and then rigid body motion was analyzed in the body subsystem.The rolling,yawing,and pitching motions were simulated using the proposed method as follows:In the rolling simulation,the difference of the stroke angle between the right and left wings caused a roll torque.In the yawing simulation,the initial feathering angle in the right wing only caused a yaw torque.In the pitching simulation,the difference between the front-and back-stroke angles in both the right and left wings caused a pitch torque.All three torques generated maneuvering motion comparable with that obtained in actual observations of insect flight.These results demonstrate that the proposed method can adequately simulate the fundamental maneuvers of insect flapping flight.In the present simulations,the maneuvering mechanisms were investigated at the governing equation level,which might be difficult using other approaches.Therefore,the proposed method will contribute to revealing the underlying insect flight mechanisms.

Parallel hydrodynamic finite element model with an N-Best refining partition scheme

We enhance a robust parallel finite element model for coasts and estuaries cases with the use of N-Best refinement algorithms,in multilevel partitioning scheme.Graph partitioning is an important step to construct the parallel model,in which computation speed is a big concern.The partitioning strategy includes the division of the research domain into several semi-equal-sized sub-domains,minimizing the sum weight of edges between different sub-domains.Multilevel schemes for graph partitioning are divided into three phases:coarsening,partitioning,and uncoarsening.In the uncoarsening phase,many refinement algorithms have been proposed previously,such as KL,Greedy,and Boundary refinements.In this study,we propose an N-Best refinement algorithm and show its advantages in our case study of Xiamen Bay.Compared with original partitioning algorithm in previous models,the N-Best algorithm can speed up the computation by 1.9 times,and the simulation results are in a good match with the in-situ data.

Partition of Unity for a Class of Nonlinear Parabolic Equation on Overlapping Non-Matching Grids

A class of nonlinear parabolic equation on a polygonal domain Ω R2 is inves- tigated in this paper. We introduce a finite element method on overlapping non-matching grids for the nonlinear parabolic equation based on the partition of unity method. We give the construction and convergence analysis for the semi-discrete and the fully discrete finite element methods. Moreover, we prove that the error of the discrete variational problem has good approximation properties. Our results are valid for any spatial dimensions. A numerical example to illustrate the theoretical results is also given.

EXPANDABLE PARALLEL FINITE ELEMENT METHODS FOR LINEAR ELLIPTIC PROBLEMS

In this article,two kinds of expandable parallel finite element methods,based on two-grid discretizations,are given to solve the linear elliptic problems.Compared with the classical local and parallel finite element methods,there are two attractive features of the methods shown in this article:1)a partition of unity is used to generate a series of local and independent subproblems to guarantee the final approximation globally continuous;2)the computational domain of each local subproblem is contained in a ball with radius of O(H)(H is the coarse mesh parameter),which means methods in this article are more suitable for parallel computing in a large parallel computer system.Some a priori error estimation are obtained and optimal error bounds in both H^1-normal and L^2-normal are derived.Finally,numerical results are reported to test and verify the feasibility and validity of our methods.

CONVERGENCE ANALYSIS OF SOME FINITE ELEMENT PARALLEL ALGORITHMS FOR THE STATIONARY INCOMPRESSIBLE MHD EQUATIONS

By combination of iteration methods with the partition of unity method(PUM),some finite element parallel algorithms for the stationary incompressible magnetohydrodynamics(MHD)with different physical parameters are presented and analyzed.These algorithms are highly efficient.At first,a global solution is obtained on a coarse grid for all approaches by one of the iteration methods.By parallelized residual schemes,local corrected solutions are calculated on finer meshes with overlapping sub-domains.The subdomains can be achieved flexibly by a class of PUM.The proposed algorithm is proved to be uniformly stable and convergent.Finally,one numerical example is presented to confirm the theoretical findings.

Meshfree Method Numerical Simulation of Sheet Metal Forming Process

Metal forming plays an important role in manufacturing industry and is widely applied in industries.The tradi- tional finite element method(FEM)numerical simulation is commonly used to predict metal forming process.Conventional finite element analysis of metal forming processes often breaks down due to severe mesh distortion,therefore time-consuming remeshing is necessary.Meshfree methods have been developed since 1977 and can avoid this problem.This new generation of computational methods reduces time-consuming model generation and refinement effort,and its shape function has higher order connectivity than FEM’s.In this paper the velocity shape functions are developed from a reproducing kernel approximation that satisfies consistency conditions and is used to analyze metal tension rigid viscoplastic deforming and Magnesium Alloy(MB 15)sheet superplastic ten- sion forming.A meshfree method metal forming modeling program is set up,the partition of unity method is used to compute the integrations in weak form equations and penalty method is used to impose the essential boundary condition exactly.Metal forming examples,such as sheet metal superplastic tension forming and metal rigid viscoplastic tension forming,are analyzed to demon- strate the performance of mesh free method.

A novel virtual node method for polygonal elements

A novel polygonal finite element method （PFEM） based on partition of unity is proposed, termed the virtual node method （VNM）. To test the performance of the present method, numerical examples are given for solid mechanics problems. With a polynomial form, the VNM achieves better results than those of traditional PFEMs, including the Wachspress method and the mean value method in standard patch tests. Compared with the standard triangular FEM, the VNM can achieve better accuracy. With the ability to construct shape functions on polygonal elements, the VNM provides greater flexibility in mesh generation. Therefore, several fracture problems are studied to demonstrate the potential implementation. With the advantage of the VNM, the convenient refinement and remeshing strategy are applied.

OPTIMAL ERROR ESTIMATES OF THE PARTITION OF UNITY METHOD WITH LOCAL POLYNOMIAL APPROXIMATION SPACES

In this paper, we provide a theoretical method（PUFEM）, which belongs to the analysis of the partition of unity finite element family of meshfree methods. The usual error analysis only shows the order of error estimate to the same as the local approximations[12]. Using standard linear finite element base functions as partition of unity and polynomials as local approximation space, in l-d case, we derive optimal order error estimates for PUFEM interpolants. Our analysis show that the error estimate is of one order higher than the local approximations. The interpolation error estimates yield optimal error estimates for PUFEM solutions of elliptic boundary value problems.

A Quadratic Serendipity Finite Volume Element Method on Arbitrary Convex Polygonal Meshes

Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.

Superconvergence of Rectangular Mixed Finite Element Methods for Constrained Optimal Control Problem

We investigate the superconvergence properties of the constrained quadratic elliptic optimal control problem which is solved by using rectangular mixed finite element methods.We use the lowest order Raviart-Thomas mixed finite element spaces to approximate the state and co-state variables and use piecewise constant functions to approximate the control variable.We obtain the superconvergence of O(h^(1+s))(0<s≤1)for the control variable.Finally,we present two numerical examples to confirm our superconvergence results.

Methodology for Comparing Coupling Algorithms for Fluid-Structure Interaction Problems

The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.

On hp refinements of independent cover numerical manifold method——some strategies and observations

In this paper,strategies are provided for a powerful numerical manifold method(NMM)with h and p refinement in analyses of elasticity and elasto-plasticity.The new NMM is based on the concept of the independent cover,which gets rid of NMM's important defect of rank deficiency when using higher-order local approximation functions.Several techniques are presented.In terms of mesh generation,a relationship between the quadtree structure and the mathematical mesh is established to allow a robust h-refinement.As to the condition number,a scaling based on the physical patch is much better than the classical scaling based on the mathematical patch;an overlapping width of 1%–10%can ensure a good condition number for 2nd,3rd,and 4th order local approximation functions;the small element issue can be overcome after the local approximation on small patch is replaced by that on a regular patch.On numerical accuracy,local approximation using complete polynomials is necessary for the optimal convergence rate.Two issues that may damage the convergence rate should be prevented.The first is to approximate the curved boundary of a higher-order element by overly few straight lines,and the second is excessive overlapping width.Finally,several refinement strategies are verified by numerical examples.

基于无单元局部弱式法的2.5维直流电阻率正演

传统的数值模拟方法依赖于节点连接信息的单元,基于全局弱式的无单元Galerkin法虽然无需节点连接信息的单元剖分,但仍依赖于全局域内剖分的背景积分单元。文中基于无单元Galerkin法,利用单位分解积分法将全局域积分转化为节点局部域积分,从而不再需要全局域内剖分的背景积分单元,进一步减小了对单元的依赖,实现了更灵活和更高精度的2.5维直流电阻率无单元局部弱式法正演模拟。分别应用该算法、无单元Galerkin法及有限单元法对不同地电模型进行数值模拟及效果对比,证明了所提算法对直流电阻率正演模拟的正确性和有效性。相比于无单元Galerkin法和有限单元法,文中算法具有更强的灵活性和适应性、更高的精度。为了提高该算法对地形起伏模型的模拟精度,对地形起伏区域进行节点加密,并对地形以外的高斯点不做计算,可获得与无单元Galerkin法和有限单元法(FEM)基本一致的效果。

考虑系统参数和载荷不确定性的声振分析

为了准确描述随机性对声振系统响应的影响,本文基于单位分解有限元/有限元法和直接概率积分法实现了精确高效的随机声振分析。本文首先通过在声振耦合界面划分虚拟单元、实施坐标变换和引入精确积分方法,发展了新型单位分解有限元/有限元方法。然后将新方法的单元内捕捉波形信息的能力强、系统矩阵规模小、易于处理复数矩阵等优势和直接概率积分法需要重复计算的特点相结合,构建了精确高效的随机声振分析方案。数值算例分析表明,所提方案得到的统计特性均和蒙特卡洛模拟吻合很好,其能够进行可靠的随机响应分析。

FINITE ELEMENT APPROXIMATION OF A NONLINEAR STEADY-STATE HEAT CONDUCTION PROBLEM

Examines a nonlinear partial differential equation of elliptic type with the homogeneous Dirichlet boundary conditions; Proof of the comparison and maximum principles; Approximation of the finite element; Introduction of a discrete analogue of the maximum principle for linear elements.

FINITE VOLUME ELEMENT METHOD WITH LAGRANGIAN CUBIC FUNCTIONS

这份报纸在三角形的网孔上为泊松方程建立一个新有限体积元素计划。试用函数空间在三角形的分区上作为 Lagrangian 立方的有限元素空格被花，并且测试函数空格在双分区上被定义为 piecewise 常数空格。在关于三角形的网孔的一些弱状况下面，作者证明僵硬矩阵是一致地积极的明确并且是 O 的集中率(h 3 ) 在 H 1 标准。一些数字实验证实理论考虑。