Least-Squares Finite Element Method for the Steady Upper-Convected Maxwell Fluid

Shaoling Zhou^1,2*, Lei Hou¹

¹Department of Mathematics, Shanghai University, Shanghai, China.
²School of Science, Hebei University of Engineering, Handan, China.

DOI: 10.4236/apm.2015.55024                

Abstract

In this paper, a least-squares finite element method for the upper-convected Maxell (UCM) fluid is proposed. We first linearize the constitutive and momentum equations and then apply a least-squares method to the linearized version of the viscoelastic UCM model. The L² least-squares functional involves the residuals of each equation multiplied by proper weights. The corresponding homogeneous functional is equivalent to a natural norm. The error estimates of the finite element solution are analyzed when the conforming piecewise polynomial elements are used for the unknowns.

Keywords

Upper-Convected Maxwell Fluid, Least-Squares Finite Element Method, Viscoelastic Fluid Model

1. Introduction

In recent years, there has been an increased interest in the least-squares finite element method for the approximation of partial differential equations, see e.g. [1] -[6] . This technique is attractive because the linear systems generated by the discretization are symmetric and positive definite, thus the algebraic system can be solved by fast direct or iterative algorithms. Moreover, in contrast to the mixed finite element method, the inf-sup or LBB type of conditions is naturally satisfied. However, without the weights in the least-squares functional, this method results in poor numerical solutions even for simple problems. In [7] , Bochev and Gunzburger pointed that the weighted least-squares method was optimal for the velocity-pressure-stress formulation of the Stokes equations. The weighted least-squares method has also been used to solve other viscoelastic problems, such as the Oldroyd-B, Carreau, and Phan-Thien-Tanner models [1] [3] [8] .

In the viscoelastic fluids of the differential type, the constitutive equations consist of an algebraic tensorial relationship between the stress tensor and the rate of deformation tensor. The upper-convected Maxwell fluid [9] is the simplest, if not the easiest, representative of that class and has served as a model fluid for developing numerical techniques. The purpose of this paper is to present a finite element method for the upper-convected Maxwell fluid which is one of the most used viscoelastic models. The nonlinear model is first approached by linearizing the equations and a weighted least-squares finite element method is applied to solve the linear equations. Error estimates of the finite element solutions to the linear system are derived.

2. Governing Equations

Assume that is a bounded and connected domain in with Lipschitz boundary. We consider the steady incompressible flows governed by the conservation equations for mass and momentum

(2.1)

(2.2)

where denotes the velocity vector, the constant density, the pressure, the extra-stress tensor and the body force.

For the upper-convected Maxwell model, the extra-stress tensor satisfies the following constitutive equation

(2.3)

where is the constant viscosity, the relaxation time and the standard strain rate tensor. The subscript (1) denotes the upper-convected material derivative

and

To simplify our analysis, homogeneous boundary conditions are assumed on. The results in this article can be extended to nonhomogeneous boundary conditions easily. Collecting (2.1)-(2.3), we obtain the steady UCM model

(2.4)

3. Mathematical Notation and Preliminaries

Throughout the paper, we use the standard notation and definition for the Sobolev spaces, , with inner products and norms denoted by and. For, we write, and. As usual, denotes the closure of with respect to the norm and denotes the space of squ- are integrable functions with zero mean

.

The spaces with positive values of is defined as the dual space of with the following norm

(3.1)

where stands for the duality pairing between and.

We use the following approximation

to linearize the equations in (2.4). Moreover we assume that and the approximation satisfies

. (3.2)

We introduce the replacement rules

which result in the linearized system

(3.3)

where

The velocity, the pressure and the extra-stress tensor belong to their respective spaces

and let.

Based on [7] , we define the weighted least-squares functional for the linearized system (3.3)

(3.4)

Now we show that the homogeneous least-squares functional of (3.4) is equivalent to the norm

.

Theorem 1. There exist positive constants and, which depend on, and, such that

(3.5)

hold for any.

Proof. The upper bound in (3.5) follows easily from the triangle inequality and (3.2). For the lower bound, we will show that

(3.6)

where

.

Using the Green’s formula and Cauchy-Schwarz inequality, we obtain

for any. By using Lemma 2.1 in [10] as

we obtain

which implies that

. (3.7)

Similarly, we have

(3.8)

By the arguments similar to Theorem 4.1 in [11] , we obtain

(3.9)

Combining (3.7)-(3.9) yields (3.6).

From the inequality, we establish

(3.10)

Note that

and, using (3.10),

where for and chosen sufficiently small. This completes the proof.

However the least squares functional is not practical. The negative order Sobolev norm

leads to difficulties in the assembly of the linear algebraic equations. In [7] , Bochev and Gunzburger used the weighted norm where denotes some parameter of the finite element space instead of the norm. Hence we will consider the mesh dependent functional in which the residuals of each equation in -norm are multiplied by proper mesh dependent weights.

4. Finite Element Approximations

We assume that the domain is a polygon, and is a triangulation of made of triangular elements with. Thus, the computational domain is given by

.

We assume that the partition is regular and satisfies the inverse assumption. Let denote the space of polynomials on of degree less or equal to. We define the finite element spaces for the approximation of as follows

.

Let be a finite dimensional subspace of X with the following approximation properties:

(3.11)

where and. The space admits the property

(3.12)

with and. The properties hold for finite element spaces consisting of continuous piecewise polynomials based on quasi-uniform triangulations [11] [12] .

The mesh dependent least squares functional is defined by the weighted sum in -norms of the residuals of the equations in (3.3)

(3.13)

where is a positive constant. The least squares finite element problem is to minimize this functional over X^h: seek such that

.

The minimizer of (3.13) necessarily satisfies the Euler-Lagrange equation given by

(3.14)

where

and the double-dot product is defined as

Based on Theorem 1, we establish the ellipticity of the functional in Theorem 2.

Theorem 2. For any, there exist positive constants and, which depend on, and, such that

, (3.15)

for.

Proof. The first inequality in (3.15) is straightforward from Theorem 1. To prove the upper bound, we assume that the spaces and satisfy the following inverse inequalities

and

.

From the triangle inequality, we obtain

This completes the proof of the theorem.

By virtue of Theorem 2 and the Lax-Milgram theorem, we establish the following theorem.

Theorem 3. For any, the functional (3.13) has the unique minimizer out of, i.e., there exists a unique solution satisfies the Euler-Lagrange equation (3.14).

Now we derive error estimates for the least-squares finite element solution which satisfies (3.14).

Theorem 4. Assume that is the solution to (3.3), then the least-squares finite element solution satisfies

, (3.16)

for.

Proof. From Theorem 2, we obtain the following bound

.

Combining the properties (3.11) and (3.12), we have

This completes the proof of the theorem.

5. Conclusion

In this paper, we have proposed and analyzed a weighted least-squares method for the approximate solution of the upper-convected Maxwell fluid. The weights in our least-squares functional involve mesh dependent weight and mass conservation constant. The homogenous functional is shown to be equivalent to a natural norm. A prior error estimate is given for the finite element solutions. An adaptive least-squares finite element method for this viscoelastic fluid model will be discussed in the future.

Acknowledgements

The authors’ work is supported by the National Science Foundation of China (No. 11271247) and the Natural Science Foundation of Hebei Province (No. G2013402063).

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Chen, T.F., Lee, H. and Liu, C.C. (2013) Numerical Approximation of the Oldroyd-B Model by the Weighted Least- Squares/Discontinuous Galerkin Method. Numerical Methods for Partial Differential Equations, 29, 531-548. http://dx.doi.org/10.1002/num.21719
[2]	Cai, Z. and Ku, J. (2006) The L2 Norm Error Estimates for the Div Least-Squares method. SIAM Journal on Numerical Analysis, 44, 1721-1734. http://dx.doi.org/10.1137/050636504
[3]	Zhou, S.L. and Hou, L. (2015) Decoupled Algorithm for Solving Phan-Thien-Tanner Viscoelastic Fluid by Finite Element Method. Computer & Mathematics with Applications, 69, 423-437. http://dx.doi.org/10.1016/j.camwa.2015.01.006
[4]	Cai, Z., Lazarov, R. and Manteuffel, T.A. and McCormick, S.F. (1994) First-Order System Least Squares for Second- Order Partial Differential Equations: Part I. SIAM Journal on Numerical Analysis, 31, 1785-1799. http://dx.doi.org/10.1137/0731091
[5]	Cai, Z., Lee, B. and Wang, P. (2004) Least-Squares Methods for Incompressible Newtonian Fluid Flow: Linear Stationary Problems. SIAM Journal on Numerical Analysis, 42, 843-859. http://dx.doi.org/10.1137/S0036142903422673
[6]	Lee, H.C. and Chen, T.F. (2015) Adaptive Least-Squares Finite Element Approximations to Stokes Equations. Journal of Computational and Applied Mathematics, 280, 396-412. http://dx.doi.org/10.1016/j.cam.2014.11.041
[7]	Bochev, P.B. and Gunzburger, M.D. (1995) Least-Squares Methods for the Velocity-Pressure-Stress Formulation of the Stokes Equations. Computer Methods in Applied Mechanics and Engineering, 126, 267-287. http://dx.doi.org/10.1016/0045-7825(95)00826-M
[8]	Lee, H.C. (2014) An Adaptively Refined Least-Squares Finite Element Method for Generalized Newtonian Fluid Flows Using the Carreau Model. SIAM Journal on Scientific Computing, 36, A193-A218. http://dx.doi.org/10.1137/130912682
[9]	Fan, Y., Tanner, R.I. and Phan-Thien, N. (1999) Galerkin/Least-Square Finite-Element Methods for Steady Viscoelastic Flows. Journal of Non-Newtonian Fluid Mechanics, 84, 233-256. http://dx.doi.org/10.1016/S0377-0257(98)00154-2
[10]	Cai, Z., Manteuffel, T.A. and McCormich, S.F. (1995) First-Order System Least Squares for Velocity-Vorticity-Pres- sure from of the Stokes Equations, with Application to Linear Elasticity. Electronic Transactions on Numerical Analysis, 3, 150-159.
[11]	Cai, Z. and Westphal, C.R. (2009) An Adaptive Mixed Least-Squares Finite Element Method for Viscoelastic Fluids of Oldroyd Type. Journal of Non-Newtonian Fluid Mechanics, 159, 72-80. http://dx.doi.org/10.1016/j.jnnfm.2009.02.004
[12]	Braess, D. (2007) Finite Elements: Theory, Fast Solvers, and Applications in Solid Mechanics. Cambridge University Press, Cambridge. http://dx.doi.org/10.1017/CBO9780511618635