有限元计算中疏密网格过渡方法研究

工程计算中出于节省计算量的目的,往往需要在一个有限元模型中布置粗细不同的网格。为保证计算结果的准确性,必须保证网格突变情况下的位移协调问题。本文工作之一是在强天驰界面过渡单元的基础上,引入虚拟节点和子单元,在子单元中应用节理元思想,提出了基于最小势能原理的弹簧节理单元法。简化了积分运算,避免了精度要求极高的坐标转换,从而提高了方法的精度和实用性;二是提出了基于位移约束的主从自由度法,简便实用,只需简单的矩阵运算即可实现。两种方法均实现了不同尺寸网格间位移的协调性和刚度的匹配,从而使之满足有限元收敛准则,且生成的刚度阵具有对称性及带状性。算例证明两种方法精度良好,并可方便地应用于求解大规模工程问题。     

CGNS API和FVM在非结构混合网格计算中的应用

用CGNS API（CFD General Notation System Application Programming Interface）作为非结构混合网格求解器的前处理和后处理,用FVM（Finite Volume Method）作为偏微分方程求解方法．在前处理过程中,用hash表法对内部网格面和边界网格面进行编号,并计算出相应的几何信息,以满足FVM求解器的需要．从FVM求解器计算出来的各种场信息可以写入原来的CGNS文件,该文件可以被许多专业商业后处理软件（如Tecplot,Fluent,CFX等）读取和进行可视化;对于求解器,用基于网格中心的FVM及SIMPLEC（Semi Implicit Methodfor Pressure Linked Equation Consistent）方法求解压力速度耦合．最后给出两个说明算例．     

热传导方程的一类无网格方法

构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.     

A stabilized extremum-preserving scheme for nonlinear parabolic equation on polygonal meshes

In this paper, a stabilized extremum-preserving scheme is introduced for the nonlinear parabolic equation on polygonal meshes. The so-called harmonic averaging points located at the interface of heterogeneity are employed to define the auxiliary unknowns and can be interpolated by the cell-centered unknowns. This scheme has only cell-centered unknowns and possesses a small stencil. A stabilized term is constructed to improve the stability of this scheme. The stability analysis of this scheme is obtained under standard assumptions. Numerical results illustrate that the scheme satisfies the extremum principle with anisotropic full tensor coefficient problems and has optimal convergence rate in space on distorted meshes.     

Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition

Fully discrete discontinuous Galerkin methods with variable mesh- es in time are developed for the fourth order Cahn-Hilliard equation arising from phase transition in materials science. The methods are formulated and analyzed in both two and three dimensions, and are proved to give optimal order error bounds. This coupled with the flexibility of the methods demonstrates that the proposed discontinuous Galerkin methods indeed provide an efficient and viable alternative to the mixed finite element methods and nonconforming (plate) finite element methods for solving fourth order partial differential equations.

     

CONVERGENCE OF THE NONCONFORMING WILSON''''S BRICK FOR ARBITRARY HEXAHEDRAL MESHES

The nonconforming Wilson's brick classically is restricted to regular hexahedral meshes. Lesaint and Zlamal[6] relaxed this constraint for the two-dimensional analonue of this element In this paper we extend their results to three dimensions and prove that and where u is the exact solution, u_h is the approximate solution and is the usual norm for the Sobolev space H~1(?).     

一个Crouzeix-Raviart型有限元方法关于二阶椭圆问题超逼近的一个新证明

In this paper,a new proof of superclose of a Crouzeix-Raviart type finite element is given for second order elliptic boundary value problem by Bramble-Hilbert lemma on anisotropic meshes.     

THE APPROXIMATE COMPUTATION OF HYPERSINGULAR INTEGRALS ON INTERVAL

In this paper some numerical methods for computing hypersingular integrals on interval are given. Using geometric meshes, these methods lead 10 an exponential convergence in the range of engineering compulation. A numerical example shows their effeclivity and accuracy.     

Large‐eddy simulation with complex 2‐D geometries using a parallel finite‐element/spectral algorithm

A parallel stabilized finite‐element/spectral formulation is presented for incompressible large‐eddy simulation with complex 2‐D geometries. A unique discretization scheme is developed consisting of a streamline‐upwind Petrov–Galerkin/Pressure‐Stabilized Petrov–Galerkin (SUPG/PSPG) finite‐element discretization in the 2‐D plane with a collocated spectral/pseudospectral discretization in the out‐of‐plane direction. This formulation provides an efficient approach for solving 3‐D flows over arbitrary 2‐D geometries. Utilizing this discretization and through explicit temporal treatment of the non‐linear terms, the system of equations for each Fourier mode is decoupled within each time step. A novel parallelization approach is then taken, where the computational work is partitioned in Fourier space. A validation of the algorithm is presented via comparison of results for flow past a circular cylinder with published values for Re=195, 300, and 3900. Copyright © 2003 John Wiley & Sons, Ltd.     

How Accurate is the Streamline Diffusion Finite Element Method?

We investigate the optimal accuracy of the streamline diffusion finite element method applied to convection-dominated problems. For linear/bilinear elements the theoretical order of convergence given in the literature is either for quasi-uniform meshes or for some uniform meshes. The determination of the optimal order in general was an open problem. By studying a special type of meshes, it is shown that the streamline diffusion method may actually converge with any order within this range depending on the characterization of the meshes.