地形构造中地震波传播的非对称交错网格模拟

提出了一种新的三维空间对称交错网格差分方法,模拟地形构造中弹性波传播过程．通过具有二阶时间精度和四阶空间精度的不规则网格差分算子用来近似一阶弹性波动方程,引入附加差分公式解决非均匀交错网格的不对称问题．该方法无需在精细网格和粗糙网格间进行插值,所有网格点上的计算在同一次空间迭代中完成．使用精细不规则网格处理海底粗糙界面、 断层和空间界面等复杂几何构造, 理论分析和数值算例表明, 该方法不但节省了大量内存和计算时间, 而且具有令人满意的稳定性和精度．在模拟地形构造中地震波传播时,该方法比常规方法效率更高．     

求解流固耦合问题的一种四步分裂有限元算法

基于arbitrary Lagrangian Eulerian (ALE) 有限元方法,发展了一种求解流固耦合问题的弱耦合算法．将半隐式四步分裂有限元格式推广至求解ALE描述下的Navier-Stokes（N-S）方程,并在动量方程中引入迎风流线（streamline upwind/Petrov-Galerkin, SUPG）稳定项以消除对流引发的速度场数值振荡;采用Newmark-β法对结构方程进行时间离散;运用经典的Galerkin有限元法求解修正的Laplace方程以实现网格更新,每个计算步施加网格总变形量防止结构长时间、大位移运动时的网格质量恶化．运用上述算法对弹性支撑刚性圆柱体的流致振动问题进行了数值模拟,计算结果与已有结果相吻合,初步验证了该算法的正确性和有效性．     

炸药爆轰物质点法三维数值模拟

炸药爆轰以及多点起爆所产生的爆轰波汇聚问题很难应用有限元法进行模拟分析,尤其当网格发生畸变时,导致有限元法计算效率和数值精度严重下降,甚至无法得到正确结果.为此,该文应用显式积分算法的物质点法对炸药两点起爆和按时间序列的多点起爆的爆轰过程进行数值模拟,与炸药爆轰的理论计算结果相吻合.物质点法不但可以有效地避免网格畸变问题,而且为炸药爆轰的数值模拟提供了新的思路.     

Bakhvalov-Shishkin网格上求解边界层问题的差分进化算法

该文在Bakhvalov-Shishkin网格上求解具有左边界层或右边界层的对流扩散方程,并采用差分进化算法对Bakhvalov-Shishkin网格中的参数进行优化,获得了该网格上具有最优精度的数值解.对三个算例进行了数值模拟,数值结果表明:采用差分进化算法求解具有较高的计算精度和收敛性,特别是边界层的数值解精度明显...     

Lagrange中心型守恒格式

提出了Lagrange中心型守恒气体动力学格式．引入了当前时刻子网格密度与当前时刻网格声速产生的网格分片常数压力．初始网格密度乘以初始子网格体积得到子网格质量,这些子网格质量除以当前时刻子网格体积得到当前时刻子网格密度．应用网格分片常数压力,构造了满足动量守恒、总能量守恒的Lagrange中心型守恒气体动力学格式,格点速度以与网格面的数值通量相容的方式计算．对Saltzman活塞问题等进行了数值模拟,数值结果显示Lagrange中心型守恒气体动力学格式的有效性和精确性．     

多介质问题分析的前处理程序PreGenGrid

概述了多介质问题数值模拟前处理程序PreGenGrid的主要功能,着重介绍了该前处理程序网格自动生成功能和多介质问题多块网格生成后的无缝拼接技术,给出了多介质问题网格实例,验证了前处理程序PreGenGrid的特点．PreGenGrid是科学计算和复杂工程技术人员进行数值模拟的理想工具.     

抛物型积分微分方程的变网格有限元法

众所周知,对通常的抛物型边值问题的解法就是对空间区域采用有限元法,对时间区间采用差分方法,并且不同时刻采用相同的网格。1985年,梁国平在文[5]中提出了一种解抛物型初边值问题的新思想,即变网格有限元法。此后许多数值分析学者就相继采用这一思想来解决其它问题。本文我们采用变网格有限元方法来处理抛物型积分微分方程问题,给出     

Lagrange非结构网格高阶交错型守恒气体动力学格式

提出Lagrange(拉格朗日)非结构网格高阶交错型守恒气体动力学格式.用产生于当前时刻子网格密度和网格声速的子网格压力和MUSCL方法构造了高阶子网格力,利用高阶子网格力构造了高阶空间通量,借助时间中点通量的Taylor(泰勒)展开完成了高阶时间通量离散.研制了Lagrange非结构网格高阶交错型守恒气体动力学格式.对Saltzman活塞问题等进行了数值模拟,数值结果显示了Lagrange非结构网格高阶交错型守恒气体动力学格式的有效性和精确性.     

基于高阶空间精度格式和嵌套网格对旋翼尾迹捕捉的改进

发展了一种基于高阶迎风格式和嵌套网格捕捉直升机悬停旋翼涡尾迹的方法．无粘通量采用Roe Reimann求解器,使用改进的5阶加权基本无振荡(WENO)格式对交界面左右状态进行高阶插值,并与MUSCL插值进行比较．为便于捕捉尾迹和实施周期性边界条件,计算采用结构嵌套网格,其中高质量的旋翼网格完全嵌套于背景网格中．当解达到近似收敛后在桨尖涡分布区域对背景网格进行加密,如此经过3次得到优化的背景网格．考虑到WENO格式插值的特点,提出了搜索3层洞边界和人工外边界的方法以便插值的直接进行．用该方法对一跨音速和一亚音速悬停旋翼粘性流场进行了数值计算．数值结果表明:所发展方法对涡尾迹具有很高的捕捉能力;与MUSCL格式相比,WENO格式具有较低的数值耗散．     

用CT-VOF法在数值波试验槽中生成线性和非线性波

为渡水槽中波的模拟和传播提出了二维的数值模型．假设流动的流体为粘性、不可压缩的,并将Navier-Stokes方程和连续性方程作为控制方程．用标准的k-ε模型来模拟紊流流动;用交错网格的有限差分法,离散化Navier-Stokes方程;并用简化的标识和单元(SMAC)方法进行求解．使用活塞型波发生器生成并传播波;数值渡水槽的端部采用敞开式的边界条件．为了证明模型的有效性,进行了一些标准的试验,如顶盖驱动的方腔测试试验、单向的常速度场试验以及干燥河床上的溃坝试验．为了论证方法的性能及其精度,将所生成波的结果与已有波理论的结果进行比较．最后,采用群集技术(CT)生成网格,并提出最佳的网格生成条件．     

Pantograph-Catenary Dynamics: An Analysis of Models and Simulation Techniques

In this paper we analyze models and simulation techniques for the interaction of pantograph and catenary. Detailed models for catenary and pantograph and the propagation of waves are first investigated. Next, the semi discretization by the finite element method and the time integration are described. In this context numerical techniques like GGL-stabilization and superconvergent patch recovery are applied. The latter yields an error estimation for the finite element grid and shows the critical points of the system.     

不可压缩流动的并行数值方法

不可压缩流动的数值模拟是计算流体力学的重要组成部分. 基于有限元离散方法, 本文设计了不可压缩Navier-Stokes (N-S)方程支配流的若干并行数值算法. 这些并行算法可归为两大类: 一类是基于两重网格离散方法, 首先在粗网格上求解非线性的N-S方程, 然后在细网格的子区域上并行求解线性化的残差方程, 以校正粗网格的解; 另一类是基于新型完全重叠型区域分解技巧, 每台处理器用一局部加密的全局多尺度网格计算所负责子区域的局部有限元解. 这些并行算法实现简单, 通信需求少, 具有良好的并行性能, 能获得与标准有限元方法相同收敛阶的有限元解. 理论分析和数值试验验证了并行算法的高效性     

ERROR ESTIMATES FOR TWO-SCALE COMPOSITE FINITE ELEMENT APPROXIMATIONS OF NONLINEAR PARABOLIC EQUATIONS

We study spatially semidiscrete and fully discrete two-scale composite finite element method for approximations of the nonlinear parabolic equations with homogeneous Dirich-let boundary conditions in a convex polygonal domain in the plane.This new class of finite elements,which is called composite finite elements,was first introduced by Hackbusch and Sauter[Numer.Math.,75(1997),pp.447-472]for the approximation of partial differential equations on domains with complicated geometry.The aim of this paper is to introduce an efficient numerical method which gives a lower dimensional approach for solving par-tial differential equations by domain discretization method.The composite finite element method introduces two-scale grid for discretization of the domain,the coarse-scale and the fine-scale grid with the degrees of freedom lies on the coarse-scale grid only.While the fine-scale grid is used to resolve the Dirichlet boundary condition,the dimension of the finite element space depends only on the coarse-scale grid.As a consequence,the resulting linear system will have a fewer number of unknowns.A continuous,piecewise linear composite finite element space is employed for the space discretization whereas the time discretization is based on both the backward Euler and the Crank-Nicolson methods.We have derived the error estimates in the L∞(L2)-norm for both semidiscrete and fully discrete schemes.Moreover,numerical simulations show that the proposed method is an efficient method to provide a good approximate solution.     

Extended finite element modeling of deformable porous media with arbitrary interfaces

In this paper, an enriched finite element method is presented for numerical simulation of saturated porous media. The arbitrary discontinuities, such as material interfaces, are encountered via the extended finite element method (X-FEM) by enhancing the standard FEM displacements. The X-FEM technique is applied to the governing equations of porous media for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. In X-FEM, the material interfaces are represented independently of element boundaries and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of elements. Finally, several numerical examples are analyzed, including the dynamic analysis of the failure of lower San Fernando dam, to demonstrate the efficiency of the X-FEM technique in saturated porous soils.     

Discrete concepts for elliptic optimal control problems with constraints on the gradient of the state

We consider elliptic optimal control problems with constraints on the gradient of the state and propose two distinguish concepts for their discretization. The first concept uses piecewise linear, continuous finite element Ansatz functions for the state, while the second concept uses the lowest order Raviart–Thomas mixed finite element. In both cases variational discretization from [5] is used for the controls. We present optimal finite element error estimates for the numerical solutions and confirm our theoretical findings by a numerical experiment. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions

In this paper we present a finite element discretization of the Joule-heating problem. We prove existence of solution to the discrete formulation and strong convergence of the finite element solution to the weak solution, up to a sub-sequence. We also present numerical examples in three spatial dimensions. The first example demonstrates the convergence of the method in the second example we consider an engineering application.     

Two-level multiscale finite element methods for the steady navier-stokes problem

In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.     

A monotone scheme for Hamilton–Jacobi equations via the nonstandard finite difference method

A usual way of approximating Hamilton–Jacobi equations is to couple space finite element discretization with time finite difference discretization. This classical approach leads to a severe restriction on the time step size for the scheme to be monotone. In this paper, we couple the finite element method with the nonstandard finite difference method, which is based on Mickens' rule of nonlocal approximation. The scheme obtained in this way is unconditionally monotone. The convergence of the new method is discussed and numerical results that support the theory are provided. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical solution of the elastic body-plate problem by nonoverlapping domain decomposition type techniques

The purpose of this paper is to provide two numerical methods for solving the elastic body-plate problem by nonoverlapping domain decomposition type techniques, based on the discretization method by Wang. The first one is similar to an older method, but here the corresponding Schur complement matrix is preconditioned by a specific preconditioner associated with the plate problem. The second one is a ``displacement-force' type Schwarz alternating method. At each iteration step of the two methods, either a pure body or a pure plate problem needs to be solved. It is shown that both methods have a convergence rate independent of the size of the finite element mesh.

     

A Lagrange multiplier/fictitious domain method for the numerical simulation of incompressible viscous flow around moving rigid bodies: (I) case where the rigid body motions are known a priori

We describe in this Note a method for the numerical simulation of incompressible viscous flow around moving rigid bodies; we suppose the rigid body motions a priori known. The computational technique takes advantage of a time discretization by operator splitting à la Marchuk-Yanenko and of a finite element space discretization on a fixed mesh, to combine a Lagrange multiplier/fictitious domain treatment of the rigid body motions with an L²-projection technique, to force the incompressibility condition. The results of numerical experiments concerning flow around moving disks at Reynolds number of the order of 100 are presented.