双曲型积分微分方程一个新H1-Galerkin混合元格式

在半离散格式下,本文针对一类双曲型积分微分方程,研究了一个新的H1-Galerkin混合有限元方法.该方法不需要满足离散的LBB条件,而且网格剖分不需要满足正则性条件.利用单元的特殊性质,在不需要使用Rita-Vrolterra投影,而是直接使用插值的情况下,得到了与传统混合有限元方法相同的误差估计,并且得到了超逼近性质.最后,通过使用插值后处理技巧,还得到了相应的超收敛结果.     

奇型非线性抛物问题的有限元误差估计

引入带权的Sobolev空间,讨论了奇型非线性抛物问题的有限元方法,在a(u),f(u)均满足Lipschitz条件下,证明了相应椭圆投影算子Rhu与u之间误差估计式,并在一定的假设下,给出了奇型非线性抛物问题的广义解u及半离散问题的有限元解uh误差估计式。     

双曲型积分微分方程混合元法的误差估计

基于Raviart-Thomas空间Vh×Wh,本文研究了双曲型积分微分方程初边值问题混合元方法的L2,L∞误差估计。给出了未知函数u,ut和乱utt伴随速度P,散度divP逼近解的最优阶L2误差估计。还得到了逼近u及P的拟最优阶L∞误差估计。     

Navier-Stokes方程的集中质量非协调有限元法

本文通过所谓的速度-压力型公式讨论了Navier-Stokes方程的集中质量非协调有限元法(半离散情形)。首先给出了所讨论方程的集中质量非协调有限元逼近格式,其次对所讨论方程的真解与逼近格式的解之间的误差进行了分析,最后利用Navier-Stokes投影算子及其性质,得到了在确定模意义下的速度、压力误差估计,且某些误差估计能达到最优。     

非线性双曲型积分微分方程的各向异性非协调有限元逼近

在各向异性网格下,讨论了一类非线性双曲型积分微分方程的一个矩形非协调有限元方法逼近,给出了半离散格式下的有限元解的收敛性分析和误差估计。在精确解适当光滑的前提下,利用新的技巧和精细估计得到了其超逼近性质。同时利用插值后处理技术导出了整体超收敛结果。本文的结论表明传统有限元分析中对网格的正则性要求和对Ritz-Volterra投影的依赖不是必要的,从而进一步扩展了非协调有限元方法的应用范围。     

四阶拟线性椭圆方程的有限元误差估计

针对一类四阶拟线性椭圆方程,本文给出了它的协调有限元逼近。当网格参数h足够小时,得到了有限元逼近解与真解之间的误差估计,并且这些误差估计是最优的。最后,通过数值实验验证了理论分析的准确性。证明方法可以类似地应用到某些二阶拟线性椭圆方程的有限元逼近。     

抛物问题的质量集中非协调有限元法

主要讨论了一类抛物问题的质量集中非协调有限元方法。首先,我们给出了所讨论问题的质量集中非协调有限元Crank-Nicolson全离散逼近格式。其次,对讨论问题的解与所给出逼近格式的解之间的误差估计进行了分析研究。最后利用椭圆投影算子,我们得到了关于L2模和能量模方面的一些误差估计式。     

非线性双曲型方程动边界问题的全离散有限元格式及数值分析

研究具有变动边界的三维区域上的非线性双曲型方程的初边值问题．提出一类垒离散有限元逼近格式．并表明了其稳定性．通过进行空问变量代换、引入椭圆投影，以及采用其它非线性做分方程先验误差估计技巧,得到了最优阶的L^1模和日^1模收敛结果。     

关于Stokes问题的一个各向异性非协调有限元的逼近方法

本文构造了一个可用于Stokes问题的各向异性非协调混合有限元,并通过引入新的证明技巧,在各向异性网格下得到了该混合元方法的最优误差估计。     

抛物型积分微分方程各向异性非协调有限元分析

各向异性有限元方法的显著的优点之一就是可以用较少的自由度得到与传统有限元正则剖分时同样的估计结果.然而,在这种情况下,Sobolev空间上的:Bramble-Hilbert引理在插值误差分析中不能直接应用,而且对于非协调元来说其传统边界估计技巧也不再适用.本文证明了一个非协调单元具有各向异性特征,并将它应用到研究抛物积分微分方程半离散格式下的Galerkin逼近.利用单元的特殊性,验证了Ritz-Volterra投影与有限元插值是相同的.在解适当光滑时,通过引入一些新的技巧,得到了与传统方法相同的收敛误差估计和超逼近性质.最后,通过构造适当的插值后处理算子,得到了各向异性网格下的整体超收敛结果.该文的结果对进一步探索和设计数值的自适应算法是有帮助的.     

Coupled and decoupled stabilized mixed finite element methods for nonstationary dual-porosity-Stokes fluid flow model

In this paper, we propose and analyze two stabilized mixed finite element methods for the dual-porosity-Stokes model, which couples the free flow region and microfracture-matrix system through four interface conditions on an interface. The first stabilized mixed finite element method is a coupled method in the traditional format. Based on the idea of partitioned time stepping, the four interface conditions, and the mass exchange terms in the dual-porosity model, the second stabilized mixed finite element method is decoupled in two levels and allows a noniterative splitting of the coupled problem into three subproblems. Due to their superior conservation properties and convenience of the computation of flux, mixed finite element methods have been widely developed for different types of subsurface flow problems in porous media. For the mixed finite element methods developed in this article, no Lagrange multiplier is used, but an interface stabilization term with a penalty parameter is added in the temporal discretization. This stabilization term ensures the numerical stability of both the coupled and decoupled schemes. The stability and the convergence analysis are carried out for both the coupled and decoupled schemes. Three numerical experiments are provided to demonstrate the accuracy, efficiency, and applicability of the proposed methods.     

A mixed finite element for the nonlinear analysis of in-plane loaded masonry walls

A mixed membrane eight-node quadrilateral finite element for the analysis of masonry walls is presented. Assuming that a nonlinear and history-dependent 2D stress-strain constitutive law is used to model masonry material, the element derivation is based on a Hu-Washizu variational statement, involving displacement, strain, and stress fields as primary variables. As the behavior of masonry structures is often characterized by strain localization phenomena, due to strain softening at material level, a discontinuous, piecewise constant interpolation of the strain field is considered at element level, to capture highly nonlinear strain spatial distributions also within finite elements. Newton's method of solution is adopted for the element state determination problem. For avoiding pathological sensitivity to the finite element mesh, a novel algorithm is proposed to perform an integral-type nonlocal regularization of the constitutive equations in the present mixed formulation. By the comparison with competing serendipity displacement-based formulation, numerical simulations prove high performances of the proposed finite element, especially when coarse meshes are adopted.     

On the analysis of a mixed mode bending sandwich specimen for debond fracture characterization

The mixed mode bending specimen originally developed for mixed mode delamination fracture characterization of unidirectional composites has been extended to the study of debond propagation in foam cored sandwich specimens. The compliance and strain energy release rate expressions for the mixed mode bending sandwich specimen are derived based on a superposition analysis of solutions for the double cantilever beam and cracked sandwich beam specimens by applying a proper kinematic relationship for the specimen deformation combined with the loading provided by the test rig. This analysis provides also expressions for the global mode mixities. An extensive parametric analysis to improve the understanding of the influence of loading conditions, specimen geometry and mechanical properties of the face and core materials has been performed using the derived expressions and finite element analysis. The mixed mode bending compliance and energy release rate predictions were in good agreement with finite element results. Furthermore, the numerical crack surface displacement extrapolation method implemented in finite element analysis was applied to determine the local mode mixity at the tip of the debond.     

特征值问题迭代伽略金法与Rayleigh商加速

该文讨论特征值问题非协调有限元和混合有限元的加速计算方法。基于迭代伽略金法和Rayleigh商加速技巧,我们建立了特征值问题Wilson非协调有限元和Ciarlet-Raviart混合有限元的加速计算方案。这些新方案把在细网格上解一个特征值问题简化为在粗网格上解一个特征值问题和在细网格上解一个线性方程。文中证明了新方案的计算结果仍然保持了渐近最优精度阶,并用数值实验验证了理论结果。     

Volumetric coupling approaches for multiphysics simulations on non‐matching meshes

In finite element analysis of volume coupled multiphysics, different meshes for the involved physical fields are often highly desirable in terms of solution accuracy and computational costs. We present a general methodology for volumetric coupling of different meshes within a monolithic solution scheme. A straightforward collocation approach is compared to a mortar‐based method for nodal information transfer. For the latter, dual shape functions based on the biorthogonality concept are used to build the projection matrices, thus further reducing the evaluation costs. We give a detailed explanation of the integration scheme and the construction of dual shape functions for general first‐order and second‐order Langrangian finite elements within the mortar method, as well as an analysis of the conservation properties of the projection operators. Moreover, possible incompatibilities due to different geometric approximations of curved boundaries are discussed. Numerical examples demonstrate the flexibility of the presented mortar approach for arbitrary finite element combinations in two and three dimensions and its applicability to different multiphysics coupling scenarios. Copyright © 2016 John Wiley & Sons, Ltd.     

多孔介质中两相可混溶驱动问题的特征-混合动态有限元方法(英文)

数学上,多孔介质中一种不可压流体对另一不可压流体的相溶驱动由两个耦合的非线性偏微分方程组成,其中一个是关于压力的椭圆方程,另一个是关于浓度的抛物方程。本文用特征有限元方法结合动态有限元空间来逼近浓度,而压力和达西速度则由混合元方法来同时逼近。通过采用负模估计,我们给出了收敛性分析与误差估计。     

A Computational Study of Stabilized, Low-order C 0 Finite Element Approximations of Darcy Equations

We consider finite element methods for the Darcy equations that are designed to work with standard, low order C ⁰ finite element spaces. Such spaces remain a popular choice in the engineering practice because they offer the convenience of simple and uniform data structures and reasonable accuracy. A consistently stabilized method [20] and a least-squares formulation [18] are compared with two new stabilized methods. The first one is an extension of a recently proposed polynomial pressure projection stabilization of the Stokes equations [5,13]. The second one is a weighted average of a mixed and a Galerkin principles for the Darcy problem, and can be viewed as a consistent version of the classical penalty stabilization for the Stokes equations [8]. Our main conclusion is that polynomial pressure projection stabilization is a viable stabilization choice for low order C ⁰ approximations of the Darcy problem.     

Stabilized low‐order finite elements for strongly coupled poromechanical problems

Finite element solutions of poromechanical problems often exhibit oscillating pore pressures in the limits of low permeability, fast loading rates, coarse meshes, and/or small time step sizes. To suppress completely the pore pressure oscillations, a stabilized finite element scheme with a better performance on monotonicity is proposed for modeling compressible fluid‐saturated porous media. This method, based on the polynomial pressure projection technique, allows the use of linear equal‐order interpolation for both displacement and pore pressure fields, which is more straightforward for both code development and maintenance compared to others. By employing the discrete maximum principle, a proper stabilization parameter is deduced, which is efficient to guarantee the monotonicity and optimal in theory in the 1‐dimensional case. An appealing feature of the method is that the stabilization parameter is evaluated in terms of the properties of porous material only, while no mesh or time step size is involved. Through comparing the numerical simulations with the analytical benchmarks, the efficiency of the proposed stabilization scheme is confirmed.     

A dynamic variational multiscale method for viscoelasticity using linear tetrahedral elements

In this article, we develop a dynamic version of the variational multiscale (D‐VMS) stabilization for nearly/fully incompressible solid dynamics simulations of viscoelastic materials. The constitutive models considered here are based on Prony series expansions, which are rather common in the practice of finite element simulations, especially in industrial/commercial applications. Our method is based on a mixed formulation, in which the momentum equation is complemented by a pressure equation in rate form. The unknown pressure, displacement, and velocity are approximated with piecewise linear, continuous finite element functions. To prevent spurious oscillations, the pressure equation is augmented with a stabilization operator specifically designed for viscoelastic problems, in that it depends on the viscoelastic dissipation. We demonstrate the robustness, stability, and accuracy properties of the proposed method with extensive numerical tests in the case of linear and finite deformations.     

Finite element Lagrange multiplier formulation for size-dependent skew-symmetric couple-stress planar elasticity

We develop a variational principle based on recent advances in couple-stress theory and the introduction of an engineering mean curvature vector as energy conjugate to the couple stresses. This new variational formulation provides a base for developing a couple-stress finite element approach. By considering the total potential energy functional to be not only a function of displacement, but of an independent rotation as well, we avoid the necessity to maintain C¹ continuity in the finite element method that we develop here. The result is a mixed formulation, which uses Lagrange multipliers to constrain the rotation field to be compatible with the displacement field. Interestingly, this formulation has the noteworthy advantage that the Lagrange multipliers can be shown to be equal to the skew-symmetric part of the force-stress, which otherwise would be cumbersome to calculate. Creating a new consistent couple-stress finite element formulation from this variational principle is then a matter of discretizing the variational statement and using appropriate mixed isoparametric elements to represent the domain of interest. Finally, problems of a hole in a plate with finite dimensions, the planar deformation of a ring, and the transverse deflection of a cantilever are explored using this finite element formulation to show some of the interesting effects of couple stress. Where possible, results are compared to existing solutions to validate the formulation developed here.