Finite element approximation of viscoelastic fluid flow: Existence of approximate solutions and error bounds

Summary We study a finite element approximation of viscoelastic fluid flow obeying an Oldroyd B type constitutive law. The approximate stress, velocity and pressure are respectivelyP ₁ discontinuous,P ₂ continuous,P ₁ continuous. We use the method of Lesaint-Raviart for the convection of the extra stress tensor. We suppose that the continuous problem admits a sufficiently smooth and sufficiently small solution. We show by a fixed point method that the approximate problem has a solution and we give an error bound.This work has been supported in part by the GDR CNRS 901 Rhéologie der polymères fondus.     

A type of new posteriori error estimators for stokes problems

1. IntroductionIn the numerical approximation of PDE, it is often very importals to detect regionswhere the accuracy of the numerical solution is degraded by local singularities of the solutionof the continuous problem such as the singularity near the re-entrant corller. An obviousremedy is to refine the discretization in the critical regions, i.e., to place more gridpointswhere the solution is less regular. The question is how to identify these regions automdticallyand how to determine a goo…     

Numerical analysis to discontinuous Galerkin methods for the age structured population model of marine invertebrates

In this article we consider the age structured population growth model of marine invertebrates. The problem is a nonlinear coupled system of the age‐density distribution of sessile adults and the abundance of larvae. We propose the semidiscrete and fully‐discrete discontinuous Galerkin schemes to the nonlinear problem. The DG method is well suited to approximate the local behavior of the problem and to easily take the locally refined meshes with hanging nodes adaptively. The simple communication pattern between elements makes the DG method ideal for parallel computation. The global existence of the approximation solution is proved for the nonlinear approximation system by using the broken Sobolev spaces and the Schauder's fixed point theorem, and error estimates are obtained for both the semidiscrete scheme and the fully‐discrete scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Variational and numerical analysis of a quasistatic viscoelastic contact problem with normal compliance and friction

In this paper, a class of generalized evolution variational inequalities arising in quasistatic friction contact problem for viscoelastic materials is introduced and studied. Under some suitable assumptions, we obtain an existence and uniqueness theorem of the solution for the generalized evolution variational inequalities by using Banach’s fixed point theorem. Moreover, we study two numerical approximation schemes of the problem: semidiscrete scheme and fully discrete scheme. For both schemes, we prove the existence of the solution and derive the error estimations.     

Stokes问题的杂交-混合有限元分析

本文发展Stokes问题的一个四变量杂交-混合变分方程：应力-速度-压力-拉格朗日乘子．然后发展其有限元方法：对应四变量分别用间断型Raviart—Thomas最低阶元，分片常数元，连续线性元和连续线性元的迹空间．我们获得了稳定性和最优误差界．通过后处理办法，我们得到一个适合于计算的速度-压力格式，该格式可视为“Mini”元方法的一个变形(本文格式中引入了局部投影算子)．然而，本文格式关于压力具有“超收敛”结果：得到了压力关于H^1-范的误差界O(h)．     

A Dynamic Frictional Contact Problem with Normal Damped Response

We consider a mathematical model which describes the frictional contact between a viscoelastic body and a reactive foundation. The process is assumed to be dynamic and the contact is modeled with a general normal damped response condition and a local friction law. We present a variational formulation of the problem and prove the existence and uniqueness of the weak solution, using results on evolution equations with monotone operators and a fixed point argument. We then introduce and study a fully discrete numerical approximation scheme of the variational problem, in terms of the velocity variable. The numerical scheme has a unique solution. We derive error estimates under additional regularity assumptions on the data and the solution.     

A Proximal Bundle Method Based on Approximate Subgradients

In this paper a proximal bundle method is introduced that is capable to deal with approximate subgradients. No further knowledge of the approximation quality (like explicit knowledge or controllability of error bounds) is required for proving convergence. It is shown that every accumulation point of the sequence of iterates generated by the proposed algorithm is a well-defined approximate solution of the exact minimization problem. In the case of exact subgradients the algorithm behaves like well-established proximal bundle methods. Numerical tests emphasize the theoretical findings.     

Stabilization method for continuous approximation of transient convection problem

The aim of this work is to study a new finite element (FE) formulation for the approximation of nonsteady convection equation. Our approximation scheme is based on the Streamline Upwind Petrov Galerkin (SUPG) method for space variable, x, and a modified of the Euler implicit method for time variable, t. The most interest for this scheme lies in its application to resolve by continuous (FE) method the complex of viscoelastic fluid flow obeying an Oldroyd‐B differential model; this constituted our aim motivation and allows us to treat the constitutive law equation, which expresses the relation between the stress tensor and the velocity gradient and includes tensorial transport term. To make the analysis of the method more clear, we first study, in this article this modified method for the advection equation. We point out the stability of this new method and the error estimate of the approximation solution is discussed. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A priori and a posteriori error analyses in the study of viscoelastic problems

In this work, the numerical approximation of a viscoelastic problem is studied. A fully discrete scheme is introduced by using the finite element method to approximate the spatial variable and an Euler scheme to discretize time derivatives. Then, two numerical analyses are presented. First, a priori estimates are proved from which the linear convergence of the algorithm is derived under suitable regularity conditions. Secondly, an a posteriori error analysis is provided extending some preliminary results obtained in the study of the heat equation. Upper and lower error bounds are obtained.     

Taylor-Galerkin algorithms for viscoelastic flow: Application to a model problem

Coupled and decoupled Taylor-Galerkin algorithms are considered for viscoelastic flow and a model problem—transient startup Poiseuille flow in a channel under a fixed pressure gradient. All algorithms reproduce the steady-state solutions and are stable at high elasticity numbers (E). For a fixed mesh, the coupled and decoupled versions (TGC and TGD) give exceptional time-accuracy at low elasticity numbers [to within O(1%) at E = 1] and reasonable accuracy at high elasticity numbers [to within O(10%) at E = 10, 100]. By definition, the decoupled false-transient scheme (TGF), which uses different time scales for velocity and stress time stepping, provides a poor transient history. Where the main requirement is to compute a steady-state algorithm efficiency is crucial. The TGF scheme attains a steady state between six to eight times faster than does the TGC scheme, and the latter is over twice as fast as the TGD form. © 1994 John Wiley & Sons, Inc.     

On the convergence of inexact two-step Newton-like algorithms using recurrent functions

Using our new concept of recurrent functions, we approximate a locally unique solution of a nonlinear equation by an inexact two-step Newton-like algorithm in a Banach space setting. Our semilocal analysis provides tighter error bounds than before, and in many interesting cases, weaker sufficient convergence conditions. Applications including the solution of a nonlinear Chandrasekhar-type integral equation appearing in radiative transfer, and a two point boundary value problem with a Green kernel are also provided in this study.     

A V-cycle multigrid method for a viscoelastic fluid flow satisfying an Oldroyd-B-type constitutive equation

A V-cycle multigrid method is developed for a time-dependent viscoelastic fluid flow satisfying an Oldroyd-B-type constitutive equation in two-dimensional domains. Also existence, uniqueness, and error estimates of an approximate solution are discussed. The approximate stress, velocity, and pressure are, respectively, σ _k -discontinuous, u _k -continuous, and p _k -continuous.     

A normal compliance contact problem in viscoelasticity: An a posteriori error analysis and computational experiments

In this work, the numerical approximation of a viscoelastic contact problem is studied. The classical Kelvin-Voigt constitutive law is employed, and contact is assumed with a deformable obstacle and modelled using the normal compliance condition. The variational formulation leads to a nonlinear parabolic variational equation. An existence and uniqueness result is recalled. Then, a fully discrete scheme is introduced, by using the finite element method to approximate the spatial variable and the implicit Euler scheme to discretize time derivatives. A priori error estimates recently proved for this problem are recalled. Then, an a posteriori error analysis is provided, extending some preliminary results obtained in the study of the heat equation and other parabolic equations. Upper and lower error bounds are proved. Finally, some numerical experiments are presented to demonstrate the accuracy and the numerical behaviour of the error estimates.     

Discretization of the dynamic problem of elasticity theory with a discontinuous solution

A discretization algorithm for initial boundary-value problems is developed for systems of two linear equations of hyperbolic type with discontinuous solutions. A Crank—Nicholson scheme is constructed for discretization of the Cauchy problem and error bounds are obtained for the approximate solution. A model example is solved.Institute of Cybernetics of the Ukrainian Academy of Sciences. Kiev University. Translated from Vychislitel'naya i Prikladnaya Matematika, No. 75, pp. 75–83, 1991.     

Analysis and finite element simulations of a second‐order fluid model in a bounded domain

This article is concerned with the equations governing the steady motion of a viscoelastic incompressible second‐order fluid in a bounded domain. A new proof of existence and uniqueness of strong solutions is given. In addition, using appropriate finite element methods to approximate a coupled equivalent problem, sharp error estimates are obtained using a fixed point argument. The method is applied to the two‐dimensional lid‐driven cavity problem, at low Reynolds number and in a certain range of values of the viscoelastic parameters, to analyze the combined effects of inertia and viscoelasticity on the flow. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Nonlinear eigenvalue Neumann problems with discontinuities

In this paper we study nonlinear eigenvalue problems with Neumann boundary conditions and discontinuous terms. First we consider a nonlinear problem involving the p-Laplacian and we prove the existence of a solution for the multivalued approximation of it, then we pass to semilinear problems and we prove the existence of multiple solutions. The approach is based on the critical point theory for nonsmooth locally Lipschitz functionals.     

A unified framework for a posteriori error estimation for the Stokes problem

In this paper, a unified framework for a posteriori error estimation for the Stokes problem is developed. It is based on $[H^1_0(\Omega )]^d$ -conforming velocity reconstruction and $\underline{\varvec{H}}(\mathrm{div},\Omega )$ -conforming, locally conservative flux (stress) reconstruction. It?gives guaranteed, fully computable global upper bounds as well as local lower bounds on the energy error. In order to apply this framework to a given numerical method, two simple conditions need to be checked. We show how to do this for various conforming and conforming stabilized finite element methods, the discontinuous Galerkin method, the Crouzeix–Raviart nonconforming finite element method, the mixed finite element method, and a general class of finite volume methods. The tools developed and used include a new simple equilibration on dual meshes and the solution of local Poisson-type Neumann problems by the mixed finite element method. Numerical experiments illustrate the theoretical developments.     

ITERATIVE METHODS FOR THE BOUNDARY VALUE PROBLEM OF A THIRD ORDER DIFFERENCE EQUATION

ＩＴＥＲＡＴＩＶＥＭＥＴＨＯＤＳＦＯＲＴＨＥＢＯＵＮＤＡＲＹＶＡＬＵＥＰＲＯＢＬＥＭＯＦＡＴＨＩＲＤＯＲＤＥＲＤＩＦＦＥＲＥＮＣＥＥＱＵＡＴＩＯＮＷａｎｇＰｅｉｇｕａｎｇ（王培光）（ＨｅｂｅｉＵｎｉｖｅｒｓｉｔｙ，河北大学，邮编：０７１００２）＆Ｌｕ...     

A combined hybrid mixed element method for incompressible miscible displacement problem with local discontinuous Galerkin procedure

In this article, we propose a combined hybrid discontinuous mixed finite element method for miscible displacement problem with local discontinuous Galerkin method. Here, to obtain more accurate approximation and deal with the discontinuous case, we use the hybrid mixed element method to approximate the pressure and velocity, and use the local discontinuous Galerkin finite element method for the concentration. Compared with other combined methods, this method can improve the efficiency of computation, deal with the discontinuous problem well and keep local mass balance. We study the convergence of this method and give the corresponding optimal error estimates in L^∞(L²) for velocity and concentration and the super convergence in L^∞(H¹) for pressure. Finally, we also present some numerical examples to confirm our theoretical analysis.     

One Approach to the Comparison of Error Bounds at a Point and on a Set in the Solution of Ill-Posed Problems

The approximate solution of ill-posed problems by the regularization method always involves the issue of estimating the error. It is a common practice to use uniform bounds on the whole class of well-posedness in terms of the modulus of continuity of the inverse operator on this class. Local error bounds, which are also called error bounds at a point, have been studied much less. Since the solution of a real-life ill-posed problem is unique, an error bound obtained on the whole class of well-posedness roughens to a great extent the true error bound. In the present paper, we study the difference between error bounds on the class of well-posedness and error bounds at a point for a special class of ill-posed problems. Assuming that the exact solution is a piecewise smooth function, we prove that an error bound at a point is infinitely smaller than the exact bound on the class of well-posedness.