泊松方程及平面弹性问题有限元方法中求高阶导数的提取法

在许多有限元计算中经常在求得近似解后还要求得到近似的解的导数.如在弹性计算中,如何从计算得到的位移近似解较好地计算应力早已被研究多年.如果计算中包含直接对近似解求导数,必然会丧失部分精度,得不到满意的结果.特别,若近似解为分片常数函数,则根本无法从直接求导数得到应力的近似值.Babuska和 Miller提出了所谓“提取法”,即利用推导出来的提取公式来求解的导数的近似值,以得到与近似解本身同     

泊松方程及平面弹性问题有限元方法中求高阶导数的提取法

在许多有限元计算中经常在求得近似解后还要求得到近似的解的导数.如在弹性计算中,如何从计算得到的位移近似解较好地计算应力早已被研究多年.如果计算中包含直接对近似解求导数,必然会丧失部分精度,得不到满意的结果.特别,若近似解为分片常数函数,则根本无法从直接求导数得到应力的近似值.Babuska和 Miller提出了所谓“提取法”,即利用推导出来的提取公式来求解的导数的近似值,以得到与近似解本身同     

四阶抛物偏微分方程的H1-Galerkin混合元方法及数值模拟

到目前为止, H¹-Galerkin 混合有限元方法研究的问题仅局限于二阶发展方程. 然而对于高阶发展方程, 特别是重要的四阶发展方程问题的研究却没有出现. 本文首次提出四阶发展方程的H¹-Galerkin 混合有限元方法, 为了给出理论分析的需要, 我们考虑四阶抛物型发展方程. 通过引进三个适当的中间辅助变量, 形成四个一阶方程组成的方程组系统, 提出四阶抛物型方程的H¹-Galerkin 混合有限元方法. 得到了一维情形下的半离散和全离散格式的最优收敛阶误差估计和多维情形的半离散格式误差估计, 并采用迭代方法证明了全离散格式的稳定性. 最后, 通过数值例子验证了提出算法的可行性. 在一维情况下我们能够同时得到未知纯量函数、一阶导数、负二阶导数和负三阶导数的最优逼近解, 这一点是以往混合元方法所不能得到的.     

二阶时滞Volterra微积分方程数值研究（英文）

本文研究二阶时滞Volterra微积分方程收敛问题.利用勒让德谱方法,获得方程的精确解与近似解及精确导数与近似导数误差在指定范数空间呈指数收敛结果,推广了二阶Volterra方程的结果.     

无界区域上 Stokes 方程组的混合有限元方法

本文讨论无界区域上 Stokes 方程组边值问题的有限元近似解.为了克服区域的无界性所造成的困难,本文采用“局部化”技巧,首先将问题化为一个等价的有界区域上的边值问题,然后求解这个等价问题的混合有限元近似解,最后给出了有限元近似解的误差分析.     

定常Navier-Stokes方程的一种高效稳定有限元方法

提出了二维定常Navier-Stokes(N-S)方程的一种两层稳定有限元方法.该方法基于局部高斯积分技术,通过不满足inf-sup条件的低次等阶有限元对N-S方程进行有限元求解.该方法在粗网格上解定常N-S方程,在细网格上只需解一个Stokes方程.误差分析和数值试验都表明:两层稳定有限元方法与直接在细网格上采用的传统有限元方法得到的解具有同阶的收敛性,但两层稳定有限元方法节省了大量的工作时间.     

关于轴对称情形下粘性流体Stokes方程的有限元分析

§1 引言 本文研究轴对称Stokes方程组的有限元解法,由于问题的轴对称性,相应的有限元可以看为立体环状单元,这种环状单元虽早已被工程师们用于力学计算,然而收敛性尚未经过仔细分析。本文从柱坐标形式的Stokes方程组出发,根据问题的轴对称性,得到一个退化了的方程组,然后直接对速度和压力的轴向、径向分量进行有限元近似,讨论近似解的收敛性与误差估计。     

拟线性伪双曲型方程的变网格有限元解法

其中Ω为R~m中有界凸域。这类方程有其实际的物理背景,例如,动物神经系统中生物电讯息传播过程,就可以用一维伪双曲型方程来描述。对于问题(1.1)解的正则性,已有不少作者作了研究,且提出了求解的数值方法。本文提出了变网格有限元格式来求其近似解。这类方法已被不少作者采用,这样我们可以在不同的时间层采用不同的网格,使计算结果更好。采用本文提出的方法,不但可以求得精确解x的近似值,而且可以求得u_t的近似值。在网格变动次数M=0(1/h)时,有限元解对精确解的     

软物质准晶广义流体动力学的一个应用——圆柱绕流的近似解

文中报道了笔者建议的软物质准晶广义流体动力学的一个应用——软物质准晶圆柱绕流的近似解.人们熟知,在普通流体动力学中, 二维圆柱绕流问题遇到很大的困难,Stokes求解它,未能成功,这就是著名Stokes佯谬.为了克服这一困难, Oseen分析了原因不在边界条件的提法,也不在Stokes的求解,问题出在Navier-Stokes方程, 他对方程进行了修改, 得到了二维绕流问题的有物理意义的近似分析解.本文借助于Oseen的方法讨论12次对称软物质准晶广义流体动力学二维绕流问题.由于问题比普通流体动力学复杂得多,严格的求解,在目前的条件下是根本不可能的.笔者提出一种近似方法——交替程序去构造其零级近似解,并且把该结果用于软物质准晶的位错问题.     

带有阻尼项的Stokes方程的有限元分析

本文研究了带有阻尼项的定常Stokes方程,证明了弱解的存在唯一性得到了有限元逼近问题适定性,给出了有限元逼近误差,并提出了求解逼近解的迭代算法.数值算例表明算法是正确的和有效的.     

Semi-discrete stabilized finite element methods for Navier–Stokes equations with nonlinear slip boundary conditions based on regularization procedure

Based on the pressure projection stabilized methods, the semi-discrete finite element approximation to the time-dependent Navier–Stokes equations with nonlinear slip boundary conditions is considered in this paper. Because this class of boundary condition includes the subdifferential property, then the variational formulation is the Navier–Stokes type variational inequality problem. Using the regularization procedure, we obtain a regularized problem and give the error estimate between the solutions of the variational inequality problem and the regularized problem with respect to the regularized parameter \({\varepsilon}\), which means that the solution of the regularized problem converges to the solution of the Navier–Stokes type variational inequality problem as the parameter \({\varepsilon\longrightarrow 0}\). Moreover, some regularized estimates about the solution of the regularized problem are also derived under some assumptions about the physical data. The pressure projection stabilized finite element methods are used to the regularized problem and some optimal error estimates of the finite element approximation solutions are derived.     

On the convergence of finite element solutions to the interface problem for the Stokes system

The Stokes system with a discontinuous coefficient (Stokes interface problem) and its finite element approximations are considered. We firstly show a general error estimate. To derive explicit convergence rates, we introduce some appropriate assumptions on the regularity of exact solutions and on a geometric condition for the triangulation. We mainly deal with the MINI element approximation and then consider P1-iso-P2/P1 element approximation. Results are expected to give an instructive remark in numerical analysis for two-phase flow problems.     

A multilevel correction method for Stokes eigenvalue problems and its applications

In this paper, a new multilevel correction scheme is proposed to solve Stokes eigenvalue problems by the finite element method. This new scheme contains a series of correction steps, and the accuracy of eigenpair approximation can be improved after each step. In each correction step, we only need to solve a Stokes problem on the corresponding fine finite element space and a Stokes eigenvalue problem on the coarsest finite element space. This correction scheme can improve the efficiency of solving Stokes eigenvalue problems by the finite element method. As applications of this multilevel correction method, a multigrid method and an adaptive finite element technique are introduced for Stokes eigenvalue problems. Some numerical results are given to validate our schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

Infinite Element Approximation to Axial Symmetric Stokes Flow

We considered in [1] the finite element approximation to axial symmetric Stokes flow in a bounded domain. The problem for the flow passing an obstacle in an unbounded domain is also frequently encountered. In this paper, we are going to give approximate solutions for this problem by an approach stated in [2]. An iterative method is used to calculate the combined stiffness matrix.     

Superconvergence of a stabilized approximation for the Stokes eigenvalue problem by projection method

This paper presents a superconvergence result based on projection method for stabilized finite element approximation of the Stokes eigenvalue problem. The projection method is a postprocessing procedure that constructs a new approximation by using the least squares method. The paper complements the work of Li et al. (2012), which establishes the superconvergence result of the Stokes equations by the stabilized finite element method. Moreover, numerical tests confirm the theoretical analysis.     

An optimal order convergence for a weak formulation of the compressible Stokes system with inflow boundary condition

Summary. We formulate the compressible Stokes system given in (1.1) into a (new) weak formulation (2.1). A finite element method for this is presented. Existence and uniqueness of the finite element method is shown. An optimal error estimate for the numerical approximation is obtained. Numerical examples are given, showing its efficiency and rates of convergence of the approximate solutions that results from the discrete problem (3.1). Received October 20, 1996 / Revised version received January 21, 1999 / Published online: April 20, 2000     

非定常Stokes方程的全离散稳定化有限元格式

本文研究二维非定常Stokes方程全离散稳定化有限元方法.首先给出关于时间向后一步Euler半离散格式,然后直接从该时间半离散格式出发,构造基于两局部高斯积分的稳定化全离散有限元格式,其中空间用P_1—P_1元逼近,证明有限元解的误差估计.本文的研究方法使得理论证明变得更加简便,也是处理非定常Stokes方程的一种新的途径.     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS I： SPATIAL DISCRETIZATION

In this article we consider a two-level finite element Galerkin method using mixed finite elements for the two-dimensional nonstationary incompressible Navier-Stokes equations. The method yields a $H^1$-optimal velocity approximation and a $L_2$-optimal pressure approximation. The two-level finite element Galerkin method involves solving one small, nonlinear Navier-Stokes problem on the coarse mesh with mesh size $H$, one linear Stokes problem on the fine mesh with mesh size $h << H$. The algorithm we study produces an approximate solution with the optimal, asymptotic in $h$, accuracy.     

Convergence of a finite element/ALE method for the Stokes equations in a domain depending on time

We consider the approximation of the unsteady Stokes equations in a time dependent domain when the motion of the domain is given. More precisely, we apply the finite element method to an Arbitrary Lagrangian Eulerian (ALE) formulation of the system. Our main results state the convergence of the solutions of the semi-discretized (with respect to the space variable) and of the fully-discrete problems towards the solutions of the Stokes system.     

A finite element,filtered eddy-viscosity method for the Navier–Stokes equations with large Reynolds number

The direct numerical simulation of the Navier–Stokes system in turbulent regimes is a formidable task due to the disparate scales that have to be resolved. Turbulence modeling attempts to mitigate this situation by somehow accounting for the effects of small-scale behavior on that at large-scales, without explicitly resolving the small scales. One such approach is to add viscosity to the problem; the Smagorinsky and Ladyzhenskaya models and other eddy-viscosity models are examples of this approach. Unfortunately, this approach usually results in over-dampening at the large scales, i.e., large-scale structures are unphysically smeared out. To overcome this fault of simple eddy-viscosity modeling, filtered eddy-viscosity methods that add artificial viscosity only to the high-frequency modes were developed in the context of spectral methods. We apply the filtered eddy-viscosity idea to finite element methods based on hierarchical basis functions. We prove the existence and uniqueness of the finite element approximation and its convergence to solutions of the Navier–Stokes system; we also derive error estimates for finite element approximations.