ANALYSIS OF A SEQUENTIAL REGULARIZATION METHOD FOR THE UNSTEADY NAVIER-STOKES EQUATIONS
|
|
- Nathaniel Chase
- 5 years ago
- Views:
Transcription
1 MATHEMATICS OF COMPUTATION Volume 77, Number 263, July 28, Pages S (8)287-5 Article electronically published on February, 28 ANALYSIS OF A SEQUENTIAL REGULARIZATION METHOD FOR THE UNSTEADY NAVIER-STOKES EQUATIONS XILIANG LU, PING LIN, AND JIAN-GUO LIU Abstract. The incompressibility constraint makes Navier-Stokes equations difficult. A reformulation to a better posed problem is needed before solving it numerically. The sequential regularization method (SRM) is a reformulation which combines the penalty method with a stabilization method in the context of constrained dynamical systems and has the benefit of both methods. In the paper, we study the existence and uniqueness for the solution of the SRM and provide a simple proof of the convergence of the solution of the SRM to the solution of the Navier-Stokes equations. We also give error estimates for the time discretized SRM formulation.. Introduction We consider the unsteady Navier-Stokes equations with homogeneous Dirichlet boundary condition: (.) u t +(u grad)u = ν u gradp + f, (.2) divu =, (.3) u Ω =, u t= = u, in a bounded domain Ω and a time interval t T. Here u(x,t)represents the velocity of a viscous incompressible fluid, p(x,t) represents the pressure, f the prescribed external force, u (x) the prescribed initial velocity, and with no-slip velocity boundary condition. A lot of methods have been proposed for the numerical solution of the Navier- Stokes equation. The equation is an infinite-dimensional constrained dynamical system (see [9,, 7]), where the incompressibility condition (.2) can be seen as a constraint and the pressure p as the Lagrange multiplier. The special role of the pressure p leads to computational difficulties, and a direct discretization to the system is not recommended. A reformulation or regularization is thus needed to obtain a better behaved problem, where a direct discretization may then be applied. An important structural property of a reformulation is the stability of the associated constraint or invariant set in the dynamical system context. A simple stabilization formulation is the Baumgarte stabilization (see [4]), which is also called an ambientspace (impetus-striction) formulation in [7]. The method basically replaces (.2) Received by the editor July 25, 26 and, in revised form, June 4, Mathematics Subject Classification. Primary 65M2; Secondary 76D5. Key words and phrases. Navier-Stokes equations, iterative penalty method, implicit parabolic PDE, error estimates, constrained dynamical system, stabilization method. The research was supported by several grants at the Department of Mathematics, National University of Singapore. c 28 American Mathematical Society Reverts to public domain 28 years from publication 467
2 468 XILIANG LU, PING LIN, AND JIAN-GUO LIU by α (divu) t + α 2 divu =. Under this formulation the invariant set divu = is exponentially attractive if α,α 2 >. However, when coupled with the momentum equation (.) a Poisson equation for the pressure has to be solved and it is well known that an artificial boundary condition for the pressure has to be imposed. Also, variables u and p are coupled. On the other hand, another reformulation, the penalty method (cf. [5, 7, ]), is simple and, in particular, it does not require the boundary condition for the pressure, and the variables u and p are uncoupled. However, it introduces a penalty parameter ɛ, which has to be very small to maintain the accuracy of the penalized system. This small parameter ɛ gives a sharp initial layer in the pressure and causes some inaccuracy in the pressure near the beginning of time if the initial pressure is not zero. Also, the error bound of a spatial discretization will be inversely proportional to the small parameter ɛ (cf. [4]) in general, which could be useless as ɛ is required to be very small in the formulation. In [], a sequential regularization method (SRM) is proposed and analyzed, based on methods of dealing with differential algebraic equations (see, e.g. [2, 3]). The method is a combination of the penalty method and Baumgarte stabilization and can be seen as an iterative penalty method or a modified augmented Lagrangian method for nonstationary problems (cf. [6]). The method is defined as follows: with p (x,t) an initial guess, given a small parameter ɛ and two nonnegative constants α, α 2,fors =, 2,...,solvetheproblem (.4) (.5) (.6) (u s ) t +(u s grad)u s = ν u s gradp s + f, div(α (u s ) t + α 2 u s )=ɛ(p s p s ), u s Ω =, u s t= = u. If we choose α =,α 2 =,s =andp =, the SRM is exactly the penalty method. We are interested in the SRM because it keeps the benefits of the penalty method, but unlike the penalty method, the parameter ɛ is not necessarily very small, and thus the reformulated system is more stable or less stiff (see [] and the convergence estimate later in Section 2). Furthermore, it approximates the divergence condition better than the penalty method. We can simply see from (.5) that the divergence of the velocity is of O(ɛ s ) and this bound is independent of the time t if p s p s = O(ɛ s ), which we can show later. The SRM may also avoid the initial layer in the pressure variable (see some explanation in [2, 3]). Since it does not require ɛ to be very small, it would be fine even if the error bound of the finite element method is inversely proportional to ɛ. So any standard finite elements can possibly be used in numerical computation. In addition, under this formulation, a fully explicit method can possibly be designed (see [2]). In this paper we will take α =andα 2 = α>. The domain Ω will be assumed to be smooth and bounded in R 2 (a remark will be given for the 3D case). At each SRM iteration, after eliminating p s from (.4), (.5), omitting the index s, and redefining the right-hand side as f, weneedtosolveapdeoftheform (.7) (.8) u t ɛ grad div(u t + αu) ν u +(u grad)u = f, u Ω =, u t= = u.
3 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 469 This PDE is an implicit parabolic equation, the operator I ɛ grad div is a degenerate elliptic operator, and we will show the existence of the solution to this PDE in this paper. From the formulation point of view we may find that the SRM is a more natural formulation than the penalty formulation for the unsteady Navier-Stokes equations. To see this let us consider the time-independent Stokes equations with a perturbation in the incompressibility condition (i.e., a steady-state inhomogeneous Stokes equation) ν u + gradp = f, divu = g, u Ω =. If the basic compatibility holds ( Ω g =),wehave (.9) u + gradp C ( f + gradg ), where C only depends on the domain Ω and ν (see [8]). When we try to estimate the error for the penalty method (or iterative penalty method) for the timeindependent problem, this estimation plays a crucial role. When the time is taken into account, i.e. considering u t ν u + gradp = f, divu = g, u Ω =, u t= = u, where the basic compatibility holds ( Ω g =,divu = g()), we can obtain (see [3]) (.) sup u 2 + t T ( u 2 + gradp 2 )dt C( u 2 + f 2 + gradg 2 + g t 2 (H ) dt), where (H ) is the dual space of H. If we apply the penalty method directly to the time-dependent Navier-Stokes equation, then the last term on the right-hand side of the above inequality will depend on the time derivative of the pressure (the extra requirement for the pressure reads p H (,T,(H ) ). It is not trivial to obtain the estimation of p t for the Navier-Stokes equations. In general the global compatibility condition is required (see [9]). After iterations, things will be even worse, where the regularity for higher order derivatives of p with respect to t may be needed. From this point of view, we may not expect a simple iterative penalty scheme to significantly improve its accuracy. However, if we modify the inhomogeneous time-dependent Stokes equation to what the SRM is based on, then u t ν u + gradp = f, divu t + αdivu = g, u Ω =, u t= = u. If Ω g =anddivu =,wehave(see[]) (.) sup u 2 + t T C( u 2 + ( u 2 + gradp 2 )dt f 2 + gradg 2 dt),
4 47 XILIANG LU, PING LIN, AND JIAN-GUO LIU where the constant C only depends on ν, α, the domain Ω and the time length T. The Navier-Stokes system is an index-2 partial differential algebraic system (see []) and ill-posed in time in the sense that the solution depends on g t if the equation (.2) is perturbed by a function g. From estimation (.), we can see that the solution of the sequential regularization reformulation is well posed (or more stable in time) since it does not depend on g t. Using this estimation, we can obtain an error bound for the SRM similar to that for the time-independent Navier-Stokes equations. The paper is organized as follows. In 2, we describe the sequential regularization formulation. A modified convergence theorem is proved using a simpler method than in []. In 3, we prove the existence of the solution of equation (.7), (.8). In 4, we discretize the regularized problems in the temporal direction and analyze the convergence of a number of semi-discrete schemes. Numerical experiments which show the SRM works well may be found in [,, 2, 4]. 2. Preliminaries, the sequential regularization method and its convergence We first describe some notation and assumptions. As usual, we use L p (Ω), or simply L p, to denote the space of pth-power integrable functions in Ω, and u L p =( u i p dx) p Ω i= its norm, where u =(u,,u n ). We denote the inner product in L 2 by (, ) and let L 2. H m includes functions which are square integrable up to mth derivatives with norm u H m =( D α u 2 ) 2. α m For simplicity, let m H m. For a normed linear space B with norm B and a sufficiently regular function g :[α, β] B, we define β g L 2 ([α,β];b) =( g(,t) 2 Bdt) 2 and g L ([α,β];b) = sup g(,t) B. α α t β If [α, β] =[,T], we simplify the notation as g L 2 (B) and g L (B), respectively. We list some widely used inequalities here, which may be found in many analysis textbooks (see [5]). We define C as a generic constant which does not depend on the choice of functions. Poincaré s inequality: (2.) u C gradu, u H or, more generally, for u H, ( ) (2.2) u C gradu + udx. Ω Hölder s inequality: (2.3) fgh C f L p g L q h L r where p, q, r>and p + q + r =. Ω
5 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 47 Sobolev s inequality: (2.4) u L q C u m, if 2 m n > and q = 2 m n,wheren is the dimension of the domain Ω; (2.5) u L C u m, if 2 m n <, where n is the dimension of the domain Ω; (2.6) u L 4 C u 2 u 2, where the domain is in R 2. Young s inequality: (2.7) ab ɛa p + c ɛ b q, where <p<, ɛ>, p + q =,c ɛ is a constant which only depends on ɛ. Gronwall s inequality in differential form. Let y(t) be a nonnegative, absolutely continuous function in [,t] thatsat- isfies, for almost every t, the differential inequality: (2.8) y (t) a(t)y(t)+b(t), where a(t) andb(t) are nonnegative, summable functions in [,t]. Then we have: (2.9) y(t) e t a(s)ds [y() + t ] b(s)ds. Discrete Gronwall s inequality. Let y n, a n, b n and c n be nonnegative sequences, satisfying (2.) y n + a i y + (b i y i + c i ), i= with b i <. Then we have (2.) y n + a i e C n i= bi ( c i + y ), i= where C =max i n ( b i ). We are interested in the case that the problem (.), (.2), (.3) has a unique solution and the solution satisfies (2.2) u 2 + i= i= ( u t 2 + u p 2 ) dt M, where the pressure p is up to an arbitrary constant which is determined by (2.3) p(x, ) dx =. Ω Hence, some basic compatibility condition has to be assumed (cf. [9]): (2.4) u Ω =, divu.
6 472 XILIANG LU, PING LIN, AND JIAN-GUO LIU Furthermore, we assume (2.5) sup t [,T ] f M, u 2 M, where M is a positive constant. Eliminating p s from equations (.4), (.5) (with α =,α 2 = α), we can reformulate the SRM equation (cf. []) as follows: with p (x,t) an initial guess, for s =, 2,, solve the problem (u s ) t ɛ grad((divu s) t + αdivu s )+(u s grad)u s (2.6) (2.7) (2.8) = ν u s gradp s + f, u s Ω =, u s t= = u, p s = p s ɛ ((divu s) t + αdivu s ). We take p satisfying (2.3). Then it is easy to see that p s satisfies (2.3) for all s. To simplify the nonlinear term, we introduce a few operators and inequalities (see [8, 9]): B(u, v) =(u grad)v, B(u, v) =B(u, v)+ 2 (divu)v, b(u, v, w) =(B(u, v), w), b(u, v, w) =( B(u, v), w). We can easily check that b(u, v, w) = 2 {b(u, v, w) b(u, w, v)}, u, v, w H (Ω). Therefore, b(u, v, v) =, u, v H (Ω). For trilinear forms b (or b), we can prove the following inequality by a combination of integration of parts, Hölder s and Sobolev s inequalities (see for instance [9]). Lemma 2.. Assume Ω R n. Then the trilinear form b (or b) isdefinedasa bounded linear functional on H m H m2+ H m 3,wherem i and (2.9) m + m 2 + m 3 n 2 if m i n 2, (2.2) m + m 2 + m 3 > n 2 if m i = n 2 for some i. Lemma 2.2. As a corollary of Lemma 2., we choose proper m i to get: u v w (2.2) b(u, v, w) C u 2 v w u v 2 w. Moreover, if Ω R 2, we have: (2.22) b(u, v, w) C u 2 u 2 v 2 v 2 2 w, v H 2 H, u, w H, (2.23) b(u, v, w) C u 2 u 2 v w 2 w 2, v, u, w H ; if Ω R 3, we have: b(u, v, w) C u v 2 v 2 (2.24) 2 w, v H 2 H, u, w H, where the trilinear form b can be replaced with b.
7 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 473 We will use the modified nonlinear convection term B(u, u) to replace the original one B(u, u) in the remaining parts. The Navier-Stokes equations and SRM equations are modified respectively as (2.25) (2.26) (2.27) and (2.28) (2.29) (2.3) u t + B(u, u) =ν u gradp + f, divu =, u Ω =, u t= = u, (u s ) t + B(u s, u s )=ν u s gradp s + f, div((u s ) t + αu s )=ɛ(p s p s ), u s Ω =, u s t= = u. We should mention that the systems (.)-(.3) and (2.25)-(2.27) have the same solution. Now we start to analyze the SRM (2.28)-(2.3). We have the following two lemmas (see []), and these lemmas are also true when we replace B by B. Lemma 2.3. For the solution {u s,p s } of (2.28), (2.29), (2.3), we have the following estimates: u s 2 + C M,s e C M 2,s, ( ɛ 2 (divu s) t 2 + ɛ 2 divu s 2 + (u s ) t 2 + u s p s 2 ) dt where M,s = u 2 + ( f 2 + p s 2 ) dt, andc is a generic constant which does not depend on f, u and p s. The second lemma is for a linear auxiliary problem: (2.3) (2.32) (2.33) w t + B(U, w)+ B(w, V) =ν w gradp + f, (divw) t + αdivw = g, w Ω =, w t= = u, where U, V, g and u are given functions, u satisfies the compatibility conditions (2.4) and g satisfies (2.3). Lemma 2.4. For the solution of (2.3)-(2.33), ifu and V satisfy (2.34) 2 + then we have the following estimate: (2.35) 2 2 dt M 2, w 2 + ( w w t 2 + p 2 ) dt ) T C 2 ( + M 2 )e C 2M 2 ( u 2 + ( f 2 + g 2 ) dt where C 2 does not depend on f, g and u.
8 474 XILIANG LU, PING LIN, AND JIAN-GUO LIU For Navier-Stokes equations, we have a classical energy bound (see [8]): (2.36) u 2 + ( u p 2 )dt C 3 M 3 e C 3M 2 3, where M 3 = u 2 + f 2 dt and C 3 does not depend on f, u. We will prove the convergence results by combining the above lemmas. Let C =max{c,c 2,C 3 },M 4 =2+ u 2 + M 5 = CM 4 e CM2 4,M= C(M5 +)e CM2 5, and ɛ satisfy (2.37) 4Mɛ 2 ( + 2 ( p 2 + p 2 )dt). ( f 2 +2 p 2 )dt, Theorem 2.5 (Convergence). Let u and p be the solution of problem (2.25), (2.26), (2.27) and u s and p s the solution of problem (2.28), (2.29), (2.3). If we choose our initial guess p to satisfy p 2 dt 2+2 p 2 dt (e.g. p =), and ɛ to satisfy (2.37), then for s =, 2,, we have the following error estimate: (2.38) u u s 2 + ( u u s 2 2+ (u u s ) t 2 + p p s 2 ) dt (4Mɛ 2 ) s p p 2 dt. We have seen most parts of this estimate in [], where the technique of asymptotic expansion is used in the proof. But here we are going to provide a much simpler proof following an idea in []. Here we also give clear criteria that p and ɛ should satisfy for convergence. Proof. Subtracting (2.25)-(2.27) (replacing (2.26) by divu t + αdivu = ) from (2.28)-(2.3) and denoting by e s = u s u, h s = p s p for s =, 2,, we have (2.39) (e s ) t + B(u s, e s )+ B(e s, u) =ν e s grad(h s ), (2.4) (dive s ) t + αdive s = ɛ(h s h s ), (2.4) e s Ω =, e s t= =. Clearly, from estimation (2.36), we have (2.42) u 2 + u 2 2dt M 5. We will prove the following claims by induction; the convergence is the corollary of claim (2.45). Claim: For s =, 2, 3,..., wehave (2.43) (2.44) (2.45) u s 2 + u s dt M 5, 4Mɛ 2 h s 2 dt, e s 2 + ( e s (e s ) t 2 + h s 2 )dt 4Mɛ 2 h s 2 dt.
9 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 475 Step : Case s =. Since our initial guess p satisfies p 2 dt 2+2 p 2 dt, hence M, M 4 (see Lemma 2.3 for the definition of M,s ). Claim (2.43) follows from Lemma 2.3 and the definition of ɛ implies claim (2.44). Estimates (2.42) and (2.43) together with Lemma 2.4 imply (2.46) e 2 + ( e (e ) t 2 + h 2 )dt Mɛ 2 h h 2 dt. Applying the Cauchy inequality, and noticing the choice of ɛ, wehave Mɛ 2 h h 2 dt 2Mɛ2 ( h 2 + h 2 )dt (2.47) h 2 2 dt +2Mɛ 2 h 2 dt. Inequalities (2.46) and (2.47) yield claim (2.45). Step 2: Assume that claims (2.43) and (2.45) are true for s = k and let s = k+. Claims (2.45) and (2.44) for s = k imply h k 2 dt 4Mɛ 2 h k 2 dt. This and give M,k+ = u 2 + p k 2 dt = h k + p 2 dt 2 ( f 2 + p k 2 )dt 2+ u 2 + ( p 2 + h k 2 )dt ( f 2 +2 p 2 )dt = M 4. Claim (2.43) is a corollary of Lemma 2.3 by taking s = k +. The proof of claim (2.45) is similar to the case of s =. (2.42) and (2.43) together with Lemma 2.4 imply e k+ 2 + ( e k (e k+ ) t 2 + h k+ 2 )dt Mɛ 2 h k h k+ 2 dt. Applying the Cauchy inequality, and noticing the choice of ɛ, wehave Mɛ 2 h k h k+ 2 dt 2Mɛ 2 ( h k 2 + h k+ 2 )dt h k+ 2 2 dt +2Mɛ 2 h k 2 dt. The above two inequalities yield claim (2.45). Claim (2.45) immediately implies 4Mɛ 2 h k+ 2 dt 4Mɛ 2 (4Mɛ 2 h k 2 dt) 4Mɛ 2. So the claim (2.44) is also true for s = k +. Remark 2.. As a corollary, we can obtain p s p s L2 (H ) = O(ɛ s )from inequality (2.38); hence the divergence of the velocity is of O(ɛ s ).
10 476 XILIANG LU, PING LIN, AND JIAN-GUO LIU Remark 2.2. From Theorem 2.5, the choice of ɛ should satisfy (2.37). In practice, it is not necessarily very small. A good approximation can be ensured by iterations as long as ɛ = o() (see numerical results in [, 2]). 3. Existence of the strong solution When we apply SRM to Navier-Stokes equations, we need to solve equations (2.28)-(2.3) at each iteration step. By eliminating p s from equations (2.28)-(2.3), omitting the index s and denoting the right-hand side as a new function f, we obtain equations which are similar to (.7), (.8): (3.) u t ɛ grad div(u t + αu) ν u + B(u, u) =f, (3.2) u Ω =, u t= = u, where divu =. We will establish the existence and uniqueness of the solution of equation (3.), (3.2) in this section. Following standard ideas we need to obtain energy estimates of a semi-discretized problem and then pass to the limit of the discrete solution. There are two ways to construct the discrete solution, i.e. temporally discrete or spatially discrete solution. If we use the spatially discrete solution, it is difficult to obtain enough energy or regularity estimates (since we do not have enough boundary conditions to ensure the integration by parts). So we turn to time discretization (see also [3] for using the temporally discrete solution). We first establish the energy estimations for the time discrete scheme, then let the time step approach. The limit function is the strong solution of problem (3.), (3.2). To take a semi-discrete scheme in time, we first use a semi-implicit scheme, i.e. an explicit-implicit scheme for the nonlinear convention term and an implicit scheme for the remaining. The system reads: (3.3) u n+ u n ɛ grad u n div(un+ + αu n+ ) ν u n+ + B(u n, u n+ )=f n+ with the homogeneous Dirichlet boundary condition, where u = u, f n+ = f(x,t n+ )andt n = n. Most of our results are obtained for this semi-implicit scheme. But in practice, we are also interested in a fully explicit treatment of the nonlinear convection term (thus the stiffness matrix will not change in time). The corresponding time-discrete scheme reads: (3.4) un+ u n ɛ grad u n div(un+ +αu n+ ) ν u n+ + B(u n, u n )=f n+ satisfying the same boundary and initial conditions. If we only consider a short time solution (Ω can be in either 2D or 3D), we have a similar estimation for this scheme. We will also consider a finite time solution in a remark. Clearly, equations (3.3) and (3.4) are linear second order elliptic equations for any fixed n. The existence and uniqueness of the solutions are standard. Moreover, we have stronger energy estimations for both equations. Lemma 3.. Assumethatthetimeintervalis[,T]. Let N = T, and define Au n+ = ɛ grad div( un+ u n +αu n+ ) ν u n+, n =,, 2, 3... and divu =.
11 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 477 Then there exists a constant ɛ which does not depend on u n, such that when ɛ ɛ, we have N ( u n ɛ 2 grad divun 2 + ɛ 2 grad div un u n N 2 ) C Au n 2, sup ( u n 2 + n N ɛ grad divun + ɛ grad div un u n ) C sup Au n, n N where the constant C does not depend on the choice of ɛ and u n. Proof. Define w n = Au n, p n = ɛ div( un u n + αu n ). Then (3.5) ν u n + gradp n = w n, (3.6) div( un u n + αu n )= ɛp n. Solving equation (3.6), and defining g n =divu n,wehave (3.7) ( + α)g n g n = ɛp n (3.8) gradg n gradg + ɛ gradp i. Since u is divergence free, the first term on the right-hand side of above equation is. Then we apply the Cauchy-Schwarz inequality, gradg n 2 ɛ 2 2 gradp i 2 (3.9) ɛ 2 2 n ɛ 2 T i= i= gradp i 2 gradp i 2. Using the inequality (.9) for the inhomogeneous Stokes equations (see [8]), we have (3.) u n gradp n 2 C ( w n 2 + gradg n 2 ). Defining the constant ɛ =min( 2C, T 2 2C T ), taking the summation from to n in inequality (3.) and using inequality (3.9), we have (3.) ( u i gradp i 2 ) C 2 w i 2 for any ɛ ɛ. From the definition of p i, we obtain the first conclusion in the lemma. To obtain the uniform estimation, we take sup in inequality (3.8), that is, sup gradg n sup ɛ gradp i n N n N i= ɛn gradp n. sup n N Then taking sup at inequality (3.) and choosing ɛ ɛ, we complete the second part of this lemma.
12 478 XILIANG LU, PING LIN, AND JIAN-GUO LIU Lemma 3.2. For equation (3.3), if we assume that N f n 2 is bounded, u L 2 and divu =, we have the following estimation: u n 2 + ɛ divun 2 + u i 2 C. Proof. Multiplying by u n+ on both sides of equations (3.3), we have: 2 ( un+ 2 u n 2 + u n+ u n 2 ) + 2ɛ ( divun+ 2 divu n 2 + div(u n+ u n ) 2 ) + α ɛ divun+ 2 ν gradu n+ 2 + b(u n, u n+, u n+ )=(f n+, u n+ ). Since b(u n, u n+, u n+ )=and(f n+, u n+ ) ν 2 gradun+ 2 + C f n+ 2,we can then define y n = u n 2 + +α ɛ divu n 2 to obtain the inequality y n+ y n + gradu n+ 2 C 2 f n+ 2. Taking the summation from to n in the above inequality yields (3.2) y n y + gradu i 2 C 2 f i 2. Since y = u 2 is bounded, we conclude the lemma. With more regularity assumptions for f and u, we can have a higher regularity estimation. Lemma 3.3. For equation (3.3), if we assume that N f n 2 is bounded, u H H2, divu =and ɛ ɛ, we have the following estimation: u n+ 2 + { ɛ (3.3) 2 grad div ui+ u i 2 + ɛ 2 grad divui+ 2 + u i+ 2 + ui+ u i 2 } C. Proof. In Lemma 3.2, we have proved that sup u n and n ui 2 are bounded, which will be used in the proof without so mentioning. Multiplying by un+ u n on both sides of equation (3.3), we have un+ u n 2 + ν 2 ( gradun+ 2 gradu n 2 + grad(u n+ u n ) 2 ) + ɛ div( un+ u n ) 2 + α 2ɛ ( divun+ 2 divu n 2 + div(u n+ u n ) 2 ) + b(u n, u n+, un+ u n )=(f n+, un+ u n ). Since b(u n, u n+, un+ u n ) C u n 2 u n 2 u n+ 2 u n+ 2 2 un+ u n u n 4 un+ 2 + C 2 u n u n u n+ u n+ 2 u n 4 un+ 2 + δ u n C 2 δ un 2 u n+ 2
13 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 479 and we have (3.4) (f n+, un+ u n ) u n 4 un+ 2 + f n+ 2, un+ u n 2 + ɛ ( divun+ 2 divu n 2 ) + ( gradun+ 2 gradu n 2 ) If we rewrite the equation (3.3) as C 3 f n+ 2 + δ u n C δ u n 2 u n+ 2. ν u n+ + gradp n+ = f n+ un+ u n B(u n, u n+ ), divu n+ = g n+, ( + α)g n+ g n = ɛp n+, by applying estimation (.9) for the nonhomogeneous Stokes equation, and using Young s inequality and Sobolev s inequality, we have (3.5) u n gradpn+ C ( f n+ 2 + un+ u n 2 + gradg n+ 2 + δ u n C δ u n 2 u n+ 2 ). Then (3.4) + 2C (3.5) implies 2 un+ u n 2 + ( u n gradp n+ 2 )+ 2C ɛ ( divun+ 2 divu n 2 ) + ( gradun+ 2 gradu n 2 ) (C 3 +) f n+ 2 + gradg n+ 2 +δ u n C δ u n 2 u n+ 2. Multiplying by gradg n+ on both sides of the recursive equation (+α)gradg n+ gradg n = ɛgradp n+,wehave 2 ( gradgn+ 2 gradg n 2 + grad(g n+ g n ) 2 )+α gradg n+ 2 = ɛ(gradp n+, gradg n+ ) ɛ2 2α gradpn+ 2 + α 2 gradgn+ 2. Therefore, ( gradgn+ 2 gradg n 2 )+α gradg n+ 2 ɛ2 α gradpn+ 2. Denoting y n = ɛ gradg n 2, z n = 2 ɛ divun 2 + gradu n 2,then 2 un+ u n 2 + u n α2 y n+ + α (y n+ y n )+ 2C 2C 2C (zn+ z n ) (C 3 +) f n+ 2 + ɛ 2 y n+ + δ u n C δ u n 2 zn+. If ɛ is less than some critical value ɛ (e.g. δ such that α 2C ), we can choose a properly small un+ u n 2 + u n (yn+ + z n+ y n z n ) C 4 ( f n+ 2 + u n 2 z n+ ).
14 48 XILIANG LU, PING LIN, AND JIAN-GUO LIU Applying the discrete Gronwall s inequality yields u n+ 2 + { ɛ 2 grad divui+ 2 + u i+ 2 + ui+ u i 2 } C. The last term n ɛ grad div ui+ u i 2 2 can be obtained from the identity g n+ g n = p n+ α ɛ ɛ gn+, where n pi 2 is bounded by (3.5). Now with the estimation (3.3), we can pass to the limit and then prove the existence of the strong solution. Define the strong solution to equation (3.)-(3.2) as follows. Problem 3.. For f L 2 (L 2 )andu H H 2 with divergence free, find u to satisfy u L 2 (H H 2 ), u t L 2 (L 2 ), divu t L 2 (H ), and (3.6) u t ɛ grad div(u t + αu) ν u + B(u, u) =f, (3.7) u t= = u. Theorem 3.4. There exists a unique solution u to problem 3.. Proof. We introduce the time semi-discrete solution which depends on the time step, and define u n to be a function of spatial variables satisfying the time semi-discrete equations (3.3). Then define u and U to be a piecewise constant or linear interpolation of u n, respectively, u (,t)=u n+, t [n, (n +)), (3.8) U (,t)=u n + t n (u n+ u n ), t [n, (n +)). Then we can rewrite equation (3.3) as (3.9) t U ɛ grad div( t U +u ) ν u + B(u (,t ), u )=f, where f is a constant interpolation of f n+ in the time step. Clearly, u is bounded in L 2 (H H 2 )andu is bounded in L 2 (H H 2 ) H (L 2 ). Passing to the limit of subsequence j,wehave w u j u L 2 (H H 2 ), w U j u L 2 (H H2 ) H (L 2 ). Define Q =[,T] Ω, and L 2 (Q) =L 2 (L 2 (Ω)), H (Q) ={u : u L 2 (Q), gradu L 2 (Q),u t L 2 (Q)}. Then we have L 2 (H H 2 ) H (L 2 ) H (Q) L 2 (Q), where the second embedding is compact (see [5]); hence U j u L 2 (L 2 ). Then we estimate the difference between u and U. By the definition of u and U, (3.2) u (,t) U (,t) u n+ u n, t [n, (n +)).
15 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 48 Hence, (3.2) N u U 2 dt u n+ u n 2. Then by Lemma 3.3, we get u U 2 L 2 (Q) C2. Hence, u and u are the same element in L 2 (Q). Thus we also obtain that u j u L 2 (L 2 ). The next step is to check that u is indeed a strong solution. Choosing a test function v to be smooth enough in Q and taking the inner product with v for equation (3.9), we can pass to the limit weakly in L 2 (Q) for all terms in (3.9) except the nonlinear term. Consider the nonlinear term B(u (,t ), u ). Since grad u converges to grad u weakly and u converges to u strongly in L 2 (Q), we conclude that B(u (,t ), u )convergesto B(u, u) in the sense of distribution on Q (see [8]). After passing to the limit, we obtain u satisfying equation (3.6). Then we need to verify that u satisfies the initial condition u() = u. There is a standard embedding result, L 2 (H H 2 ) H (L 2 ) C(H ) (see [5]), and then the map u u() is continuous from C(H )toh. From above, the initial condition is satisfied automatically. The last step is to check the uniqueness of the solution. Assume that we have two solutions u and v. Define w = u v. Then we have the equation of w: (3.22) w t ɛ graddiv(w t + αw) ν w + B(u, w)+ B(w, v) = with homogeneous Dirichlet boundary condition and zero initial condition. Multiplying the above equation by w t and applying Sobolev s and Young s inequalities, we have d dt ( ɛ divw 2 + gradw 2 ) C( u v 2 2) gradw 2. Since u and v are in L 2 (H 2 )andfromw() = and Gronwall s inequality we have gradw =, hence the uniqueness of solution follows. As we mentioned before, if we only look for a short time solution, the fully explicit scheme for the nonlinear convection term may be analyzed accordingly as well. Moreover, the argument works for the 3D case. The following lemma is similar to Lemma 3.3. Lemma 3.5. Assume that N f n 2 is bounded, u H H 2, divu =, and ɛ is less than a positive constant ɛ. Then there exists a positive constant T, whenever <n T =min(t,t ), and the solution to equation (3.4) satisfies (3.23) u n+ 2 + ( ɛ 2 grad div ui+ u i 2 + ɛ 2 grad divui+ 2 + u i+ 2 + ui+ u i 2 ) C.
16 482 XILIANG LU, PING LIN, AND JIAN-GUO LIU Proof. The proof will be slightly different from Lemma 3.3 since we do not have an a priori estimation for the solution of scheme (3.4). Multiplying by un+ u n on both sides of equation (3.4), similarly to Lemma 3.3, but noticing b(u n, u n, un+ u n ) C u n 3 2 u n 2 2 un+ u n u n 4 un+ 2 + δ u n C δ u n 6, we have (3.24) un+ u n 2 + ɛ ( divun+ 2 divu n 2 ) + ( gradun+ 2 gradu n 2 ) C 2 f n+ 2 + δ u n C δ u n 6. Then we rewrite the equation (3.4) as ν u n+ + gradp n+ = f n+ un+ u n B(u n, u n ), divu n+ = g n+, ( + α)g n+ g n = ɛp n+. By applying estimation (.9) for the nonhomogeneous Stokes equation, Young s and Sobolev s inequalities, we obtain (3.25) u n gradp n+ C ( f n+ 2 + un+ u n + gradg n+ 2 + δ u n C δ u n 6 ). 2 Then (3.24) + 2C (3.25) implies un+ u n 2 + ( u n gradp n+ 2 ) 2 2C + ɛ ( divun+ 2 divu n 2 )+ ( gradun+ 2 gradu n 2 ) (C 3 +) f n+ 2 + gradg n+ 2 + δ u n C δ u n 6. Multiplying by gradg n+ on both sides of the recursive equation for gradg n, ( + α)gradg n+ gradg n = ɛgradp n+, we have 2 ( gradgn+ 2 gradg n 2 + grad(g n+ g n ) 2 )+α gradg n+ 2 = ɛ(gradp n+, gradg n+ ) ɛ2 2α gradpn+ 2 + α 2 gradgn+ 2 ; hence ( gradgn+ 2 gradg n 2 )+α gradg n+ 2 ɛ2 α gradpn+ 2.
17 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 483 Define y n = ɛ 2 gradg n 2, z n = ɛ divun 2 + gradu n 2.Then 2 un+ u n 2 + u n α2 y n+ + α (y n+ y n )+ 2C 2C 2C (zn+ z n ) (C 3 +) f n+ 2 + ɛ 2 y n+ + δ u n C δ (z n ) 3. α If ɛ is less than some critical value ɛ (e.g. 2C ), we can choose a properly small δ (e.g. δ = 4C ), and denote d n = y n + z n + u n 2 2 ; hence un+ u n 2 + u n (dn+ d n ) C 4 ( f n+ 2 +(z n ) 3 ). Summing for n yields ( ui+ u i 2 + u i+ 2 2)+d n+ d + C 4 ( f i+ 2 +(d n ) 3 ). Denoting by w n = d + C 4 N f i 2 + C 4 n (dn ) 3, and noticing d n+ w n from the above inequality, we then have w n+ w n = C 4 (d n+ ) 3 C 4 (w n ) 3. Now consider a concave function F (x) = 2x. Clearly F (x) is negative and 2 monotonically increasing in (, ), and F (x) = x. By the convexity of F, 3 F (w n+ ) F (w n ) F (w n )(w n+ w n )= wn+ w n (w n ) 3 C 4. Hence F (w n ) F (w )+nc 4. SinceF (w ) is negative, we let T be the critical time that satisfies C = F (w )+C 4 T < andt =min(t,t ). As long as n T,wehavew n F (C ) and hence the estimation u n+ 2 + ( ɛ 2 grad divui+ 2 + u i+ 2 + ui+ u i 2 ) C. The last term n ɛ grad div ui+ u i 2 2 C can be obtained by the same argument as in Lemma 3.3. Remark 3.. Lemma 3.5 states the stability of the fully explicit treatment of the nonlinear convection term. With this energy estimation, we can prove the existence of a strong solution in a short time for equation (3.), (3.2) when Ω R 3. The proof will be exactly the same as that of Lemma 3.5, just replacing T by T. 4. Error estimation of the time discretization In this section, we will estimate the error between the SRM solution and the semi-implicit time discrete solution. Hence, combining it with Theorem 2.5 we can obtain an error estimate between the solution of the semi-implicit scheme and the solution of the Navier-Stokes equations (.)-(.2). Subtracting equation (3.) from equation (3.3), denoting the error function by e n = u n u(t n ), where t n = n, we obtain the equation that the error e n satisfies: (4.) e n+ e n ɛ grad e n div(en+ + αe n+ ) ν e n+ + B(u n, u n+ ) B(u(t n+ ), u(t n+ )) = r n+,
18 484 XILIANG LU, PING LIN, AND JIAN-GUO LIU tn+ where the remainder term r n+ = t n (η t n )[u tt (η) ɛ grad divu tt(η)]dη. Since = t n+ t n,wehave r n+ 2 tn+ tn+ 2 η t n 2 dη u tt (η) t n t n ɛ graddivu tt(η) 2 dη = tn+ (4.2) u tt (η) 3 t n ɛ grad divu tt(η) 2 dη. Hence we need an energy estimation for u tt(η) ɛ grad divu tt(η) 2 dη to control the truncation error. 4.. A priori estimation. Lemma 4.. Define an operator Au = ɛ grad div(u t + αu) ν u, wherethe time interval is [,T]. Then there exists a constant ɛ independent of u and when ɛ ɛ, we have t ɛ 2 grad divu t 2 + t ɛ 2 graddivu 2 + u 2 dt C( graddivu() 2 + Au 2 dt), sup ( t [,T ] ɛ grad divu t + grad divu + u ) C( graddivu() + sup Au ), ɛ t [,T ] where the constant C does not depend on ɛ and u. Proof. The proof is similar to the discrete case, Lemma 3.. Let w = Au, p = ɛ div(u t + αu) andg =divu. Firstly we solve for g from the ODE (4.3) g t + αg = ɛp with initial condition g() = divu(). The solution is (4.4) g(t) =e αt g() ɛ t p(s)e α(s t) ds. Now (u,p) satisfies the inhomogeneous Stokes equations u + gradp = w, divu = g. Then using the estimation (.9), we have u 2 + gradp C ( w + gradg ). Hence u 2 + gradp C ( graddivu() + w + ɛ t gradp(s) ds). Choosing ɛ to be small enough (ɛ ɛ = 2C T ), we can obtain (4.5) sup ( u 2 + gradp ) C( graddivu() + sup w ) t [,T ] t [,T ] and (4.6) t u gradp 2 ds C( graddivu() 2 + t w 2 ds). From the above inequalities and the definition of p, we can conclude the lemma. The following lemmas give a higher regularity estimate for the solution of equation (3.)-(3.2).
19 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 485 Lemma 4.2. Let f L (L 2 ), f t L 2 (H ), u H 2 H, divu =and ɛ ɛ. Then we have: sup ( u t + u 2 + t [,T ] ɛ grad divu t + grad divu ) C. ɛ Proof. Multiplying by u t on both sides of equation (3.), we have (4.7) u t 2 + ɛ divu t 2 ( α ɛ grad divu, u t) (ν u, u t )+ b(u, u, u t )=(f, u t ). Notice that and ( α ɛ grad divu, u t) α ɛ grad divu u t 2, (ν u, u t ) ν u u t 2, b(u, u, u t ) C u u 2 u t C 2 u u t 2, (f, u t ) f u t 2. Hence u t 2 + ɛ divu t 2 C 3 ( ɛ grad divu 2 + u f 2 ). Choosing the time t = in the above inequality, and noting that divu =, u H 2 and f L (L 2 ), we have (4.8) u t () 2 + ɛ divu t() 2 C 4. Then taking the derivative of equation (3.) with respect to t, wehave (4.9) u tt ɛ grad div(u tt + αu t ) ν u t + B(u, u t )+ B(u t, u) =f t. Multiplying by u t on both sides of the above equation yields d 2 dt u t 2 + d 2ɛ dt divu t 2 + α ɛ divu t 2 + ν gradu t 2 + b(u t, u, u t )+ b(u, u t, u t )= f t, u t. Using b(u t, u, u t ) C u 2 u t u t b(u, u t, u t )=, C 5 u 2 2 u t 2 + ν 4 gradu t 2, f t, u t ν f t 2 + ν 4 gradu t 2, we have d dt u t 2 + d ɛ dt divu t 2 + ɛ divu t 2 + gradu t 2 C 6 ( u 2 2 u t 2 + f t 2 ). Since we have already had u L 2 (H 2 ) (Lemma 2.3), applying Gronwall s inequality, we obtain u t 2 + t ɛ divu t 2 + ( α t ɛ divu t 2 + u t 2 )dt C 7 ( f t 2 dt + u t () 2 + ɛ divu t() 2 ).
20 486 XILIANG LU, PING LIN, AND JIAN-GUO LIU Applying the inequality (4.8) to control the right-hand side of the above inequality, we have the estimation for u t 2, i.e. (4.) u t 2 + ɛ divu t 2 + t ( ɛ divu t 2 + u t 2 )dt C 8. Then multiplying by Au on both sides of equation (3.), we have (4.) (u t,au)+ Au 2 + b(u, u,au) =(f,au). Since (u t,au) u t Au 2, (f,au) f Au 2, b(u, u,au) C u 3 2 u 2 2 Au C 9 4δ u 6 + δ u Au 2, we have 4 Au 2 u t 2 + C 9 4δ u 6 + δ u f 2. Taking sup for both sides of the above inequality, using Lemma 4., choosing a small δ ( 8C ) and noticing that u, f and u t are bounded, we have (4.2) sup t [,T ] Au C. Finally applying Lemma 4. and noting divu =, we obtain the estimation for all terms in the lemma. Lemma 4.3. Assume all conditions of Lemma 4.2. Moreover, f t L 2 (L 2 ), g = f() + ν u B(u, u ) H and divg =. Then we have u t 2 + ɛ 2 grad divu tt 2 + ɛ 2 grad divu t 2 + u t u tt 2 C. Proof. Defining h = u t () and letting t = in equation (3.), we have (4.3) h graddivh = g. ɛ From the regularity result for the second order elliptic equation, we have v C (Ω), w C (Ω), s.t. w w = v. ɛ Since h H(div) {u : u L 2, divu L 2 } from (4.8), we have (divh)vdx = (divh)(w ɛ w)dx = (divh)wdx (divh) ɛ divgradwdx Ω Ω Ω Ω = h gradwdx + graddivh gradwdx = (h graddivh) gradwdx Ω Ω ɛ Ω ɛ = g gradwdx = (divg)wdx =. Ω Ω Hence divh = in the distribution sense, and h H(div) implies divh =. Then h = g H ; hence u t () is bounded. Multiplying by Au t on both sides of equation (4.9), we have (u tt,au t )+ Au t 2 + b(u t, u,au t )+ b(u, u t,au t )=(f t,au t ).
21 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 487 Using (u tt,au t )= ɛ divu tt 2 + α d 2ɛ dt divu t 2 + ν d 2 dt gradu t 2, b(u t, u,au t ) C u t u 2 Au t we have 4 Au t 2 + C 2 u t 2, b(u, u t,au t ) 4 Au t 2 + C 2 u t 2, (f t,au t ) 4 Au t 2 + f t 2, ɛ divu tt 2 + ɛ Applying Gronwall s inequality, we obtain d dt divu t 2 + d dt u t 2 + Au t 2 C 3 ( u t 2 + f t 2 ). u t 2 + Au t 2 C 4. Multiplying by u tt on both sides of equation (4.9), we have u tt 2 +(Au t, u tt )+ b(u t, u, u tt )+ b(u, u t, u tt )=(f t, u tt ). Using b(u t, u, u tt ) C u t u 2 u tt 8 u tt 2 + C 5 u t 2, b(u, u t, u tt ) 8 u tt 2 + C 5 u t 2, (f t, u tt ) 4 u tt 2 + f t 2, (u tt,au t ) 4 u tt 2 + Au t 2, and combining all the above inequalities together with Lemma 4., we can obtain the results of this lemma. Remark 4... In [9], we know that for Navier-Stokes equations, if ( u tt 2 + u t 2 2)dt is bounded, then the following over-determined Neumann problem (4.4) q =div(f() B(u, u )) in Ω, (4.5) q Ω = u + f() B(u, u ) Ω has solutions. For the sequential regularization formulation (2.28)-(2.3), we have a similar property. If equations (4.4)-(4.5) are solvable, we can choose the initial guess of pressure p such that p () is the solution. Let s =, and rewrite equations (2.28)-(2.3) by eliminating p from the systems: (u ) t ɛ div((u ) t + αu ) ν u + B(u, u )=f p, u Ω =, u t= = u. Let g = f() + ν u B(u, u ) p (). Since p () satisfies equations (4.4)- (4.5), we have g H and divg =. Then following the proof of Lemma 4.3, we
22 488 XILIANG LU, PING LIN, AND JIAN-GUO LIU have div(u ) t () = ; hence p () = p (). We can then repeat this procedure to have estimations of u s for all s. 2. Since the equations (4.4)-(4.5) are over-determined, we need the compatibility condition of u and f() to fulfil the existence of a solution. This is nonlocal and virtually uncheckable for given data. In the absence of such a compatibility condition, we will discuss another type of estimation as follows. Lemma 4.4. Assume all the conditions as in Lemma 4.2. Moreover, f t L 2 (L 2 ), and we have (4.6) t u t 2 + t( Au t 2 + u tt 2 )dt C. As a consequence of inequality (4.6), ifwefixt (,T), t t, (4.7) t u t 2 + ( t ɛ 2 grad divu tt 2 + ɛ 2 grad divu t 2 + u t u tt 2 )dt C(t ). Proof. Multiplying by Au t on both sides of equation (4.9) and estimating the trilinear form by using Young s inequality, we have ɛ divu tt 2 + d ɛ dt divu t 2 + d dt u t 2 + Au t 2 C ( u t 2 + f t 2 ). Multiplying by t on both sides of the inequality, we obtain d dt t( ɛ divu t 2 + u t 2 )+t Au t 2 (C +)( u t 2 + ɛ divu t 2 + t f t 2 ). Integrating both sides and recalling the inequality (4.), we have ( ɛ divu t 2 + u t 2 )dt C, t( ɛ divu t 2 + u t 2 )+ t τ Au t 2 dτ C 3. Then multiplying by u tt on both sides of equation (4.9), we have (4.8) u tt 2 +(Au t, u tt )+b(u t, u, u tt )+b(u, u t, u tt )=(f t, u tt ). Multiplying by t on both sides of the equation, integrating from to t, and noting tb(u t, u, u tt ) C 4 u t u 2 u tt 8 t u tt 2 + C 5 t u t 2, tb(u, u t, u tt ) 8 t u tt 2 + C 5 t u t 2, t(f t, u tt ) 4 t u tt 2 + t f t 2, t(u tt,au t ) 4 t u tt 2 + t Au t 2, we obtain t τ u tt 2 dτ C.
23 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 489 Combining all the above inequalities we obtain the inequality (4.6). To obtain the inequality (4.7), we notice that graddivu t (t ) Cɛ (Lemma 4.2). Then applying Lemma 4. and replacing the initial time by t,weobtain t >t, t t ( ɛ 2 grad divu tt 2 + ɛ 2 grad divu t 2 + u t 2 2)dt C(ɛ 2 + t t Au t 2 dt). Then the conclusion is straightforward from inequality (4.6). Remark 4.2. Lemma 4.3 and Lemma 4.4 are two upper bounds of the truncation error u tt ɛ grad divu tt 2 dt. With the assumption of global compatibility, we have the energy estimation up to the initial time. Without such an assumption, the regularity remains true from a fixed small time t (i.e. t is a little bit away from ). We will mainly use Lemma 4.3 in the next section and give a remark about the result related to Lemma Error estimation. Lemma 4.5. Let f L (L 2 ), f t L 2 (L 2 ), u H 2 H, divu =and ɛ ɛ (as in Lemma 3.). Then for the semi-implicit scheme (3.3), we have the estimation (4.9) u n 2 C. Proof. Recall that we have already had a few inequalities in the previous section (see Lemma 3.3). Taking the backward difference quotient for equation (3.3), we have ( div)(un+ ɛ grad u n ν un+ u n un u n ) α ɛ grad div un+ u n + B(u n, u n+ ) B(u n, u n ) = f n+ f n. Then multiplying by un+ u n on both sides, we have ( un+ u n α ɛ div un+ u n 2 un u n 2 + un+ u n ɛ ( div un+ u n 2 div un u n un u n 2 ) 2 + div( un+ u n un u n ) 2 ) 2 + ν grad un+ u n 2 + ( b(u n, u n+, un+ u n ) b(u n, u n, un+ u n )) = ( f n+ f n, un+ u n ). Since b is a trilinear form, we have ( b(u n, u n+, un+ u n ) b(u n, u n, un+ u n )) = b(u n, un+ u n, un+ u n ) b( un u n, u n, un+ u n ).
24 49 XILIANG LU, PING LIN, AND JIAN-GUO LIU Since and b( un u n ( f n+ f n b(u n, un+ u n, u n, un+ u n, un+ u n, un+ u n )=, ) C u n 2 un+ u n un u n δ un+ u n 2 + C δ un 2 2 un u n 2, ) δ un+ u n we can choose δ to be properly small, i.e. Defining y n+ = un+ u n (4.2) 2 + ɛ div un+ u n 2 + n+ δ f f n 2, δ un+ u n 2 ν 4 grad un+ u n 2. 2,wehave (yn+ y n )+ ɛ div un+ u n 2 + un+ u n 2 C 2 ( u n 2 2(y n+ + y n )+ f n+ f n 2 ). Hence, f n+ f n 2 tn+ t n for n, weget N n= f t 2 dt due to f n+ f n = f n+ f n 2 f t 2 dt. tn+ t n f t dt. Summing Since we already had that N n= un 2 2 is bounded, summing (4.2) for n and applying the discrete Gronwall s inequality, we obtain y n C 3 ( f t 2 dt + y ). We still need to estimate y = u u 2 + ɛ div u u 2. Multiplying by un+ u n on both sides of equation (3.3) and letting n =,wehave y + α ɛ (divu, div u u )+ν(gradu, grad u u )+ b(u, u, u u )=(f, u u ). Noting divu =(sinceu = u )and(gradu, grad u u )= grad(u u ) 2 ( u, u u ), ( u, u u ) 4 y + C 4 u 2 2, (f, u u ) 4 y + f 2, b(u, u, u u ) C u 2 u u u 4 y + C u 2 2 u 2, and that u 2 and u are bounded, we obtain y C 5. Hence, y n C 6. Multiplying by Au n+ on both sides of equation (3.3), where the operator A is defined in Lemma 3., we have ( un+ u n,au n+ )+ Au n+ 2 + b(u n, u n+,au n+ )=(f n+,au n+ ).
25 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 49 Since and ( un+ u n,au n+ ) 4 Aun+ 2 + y n, (f n+,au n+ ) 4 Aun+ 2 + f n+ 2 b(u n, u n+,au n+ ) C 7 u n u n+ 2 u n+ 2 2 Au n+ 4 Aun+ 2 + δ u n C 8 δ, choosing a properly small δ and using Lemma 3. and the boundedness of y n,we can conclude that Au n is uniformly bounded. Therefore, u n 2 C from Lemma 3.. Theorem 4.6. For the semi-implicit scheme, defining the error function e n = u n u(t n ), then with the same assumption as in Lemma 4.3, we have N (4.2) e n 2 + ( ɛ 2 grad diven + e n 2 2) C 2. Proof. The error function e n satisfies the equation: e n+ e n (4.22) ɛ grad e n div(en+ + αe n+ ) ν e n+ + B(u n, u n+ ) B(u(t n+ ), u(t n+ )) = r n+. The difference of the nonlinear terms can be rewritten as: B(u n, u n+ ) B(u(t n+ ), u(t n+ )) = B(u n u n+, u n+ )+ B(u n+, e n+ )+ B(e n+, u(t n+ )). Multiplying by e n+ on both sides of equation (4.22), we now estimate the nonlinear part term by term: b(u n u n+, u n+, e n+ ) C u n u n+ u n+ 2 e n+ δ e n+ 2 + C 2 δ u n u n+ 2, b(u n+, e n+, e n+ )=, b(e n+, u(t n+ ), e n+ ) δ e n+ 2 + C 2 δ en+ 2. Choosing a properly small δ, i.e. δ e n+ 2 4 graden+ 2, and defining y n = e n 2 + ɛ diven 2,wehave (yn+ y n )+ ɛ diven+ 2 + e n+ 2 C 3 ( r n+ 2 + u n+ u n 2 + y n+ ). Noting that r n 2 C 2 (see inequality (4.2) and Lemma 4.4) and u n+ u n 2 C 2 (see Lemma 3.3) and applying the discrete Gronwall s inequality we obtain y n + ( ɛ divei 2 + e i 2 ) C 4 (y + 2 ). Since u = u(), we then have y = ; hence N (4.23) e n 2 + e n 2 C 2.
26 492 XILIANG LU, PING LIN, AND JIAN-GUO LIU To obtain a stronger error estimation, we need to multiply by Ae n+ on both sides of equation (4.22) and have ( en+ e n,ae n+ )+ Ae n+ 2 + (4.24) b(u n u n+, u n+,ae n+ ) + b(u n+, e n+,ae n+ )+ b(e n+, u(t n+ ),Ae n+ )=(r n+,ae n+ ). We now deal with the above equation term by term: ( en+ e n,ae n+ )= ɛ div en+ e n 2 + α 2ɛ ( diven+ 2 dive n 2 + div(e n+ e n ) 2 )+ ν 2 ( graden+ 2 grade n 2 + grad(e n+ e n ) 2 ), b(u n u n+, u n+,ae n+ ) 8 Aen+ 2 + C 5 u n u n+ 2 u n+ 2, b(u n+, e n+,ae n+ ) 8 Aen+ 2 + C 5 e n+ 2, b(e n+, u(t n+ ),Ae n+ ) 8 Aen+ 2 + C 5 e n+ 2, (r n+,ae n+ ) 8 Aen r n+ 2. Defining y n = ɛ diven + e n 2 and using the above inequalities, we have (y n+ y n ) (4.25) + Ae n+ 2 C 6 ( u n u n+ 2 + e n+ 2 + r n+ 2 ). From the inequality (4.2), we obtain u n u n+ 2 C 2. We have just derived e n 2 C 2 in (4.23). From condition r n 2 C 2, applying the discrete Gronwall s inequality, using e = and Lemma 3., we have N (4.26) e n 2 + ( ɛ 2 grad diven 2 + e n 2 2) C 2. We apply the semi-implicit scheme to the equations (2.28)-(2.3). Let the solution of the time discrete scheme be u n s and the solution of the Navier-Stokes equations be (u,p). If the initial value u and f() satisfy the global compatibility condition such that the overdetermined Neumann problem (4.4)-(4.5) has a solution, then we can choose our initial guess p () to solve this equation with Ω p ()dx =. Combining Theorem 2.5 and Theorem 4.6, we have our main theorem. Theorem 4.7. Let f L (L 2 ), f t L 2 (L 2 ), u H 2 H, divu =, ɛ be less than a constant ɛ, the initial guess satisfy p 2 dt 2+2 p 2 dt and p () solve equation (4.4)-4.5. Then we have u n s u(t n ) 2 + ( ɛ 2 graddiv(ui s u(t i )) 2 + u i s u(t i ) 2 2) C( 2 +(Mɛ 2 ) s ). Remark 4.3. We use r n 2 C 2 in the proof. Without a global compatibility assumption, this may not hold in general. But we should have N n r n 2 C(n ) 2, where the constant C may depend on n. In this case, after a small fixed time, the error estimation remains true.
27 ANALYSIS OF SEQUENTIAL REGULARIZATION METHOD 493 Remark 4.4. In practice we are also interested in the fully explicit treatment of the nonlinear convection term, i.e. the scheme (3.4). If we assume a priori estimation of u n, i.e. u n is uniformly bounded, we can have a similar result as in Theorem 4.7. This assumption of uniform boundedness of u n is technically difficult to be generally verified for the scheme (3.4) up to a given time (see Lemma 3.5 for the boundedness of the short time solution u n ). But it is often true in physical or practical examples. We can also consider other schemes in time, for example, the fully implicit scheme or the Crank-Nicolson scheme as in [4]. The results will be similar; i.e., with some global compatibility assumption, we have error estimations of Theorem 4.7 up to the initial time; without it the results remain true after a small fixed time. References [] C. Amrouche and V. Girault, Decomposition of vector spaces and application to the Stokes problem in arbitrary dimension, Czechoslovak Mathematical J. 44 (994), No. 9, pp MR25794 (95c:359) [2] U. Ascher and P. Lin, Sequential regularization methods for higher index DAEs with constraint singularities: Linear index-2 case, SIAM J. Numer. Anal. 33 (996) MR4855 (98c:6533) [3] U. M. Ascher and P. Lin, Sequential regularization methods for nonlinear higher index DAEs, SIAM J. Sci. Comput. 8 (997), no., 6-8. MR (97k:658) [4] J. Baumgarte, Stabilization of constraints and integrals of motion in dynamical systems, Comp. Meth. in Appl. Mech. and Eng. (972), -6. MR39628 (52:2448) [5] L. C. Evans, Partial Differential Equations. American Mathematical Society, 998. MR (99e:35) [6] M. Fortin and R. Glowinski, Augmented Lagrangian Methods: Applications to the Numerical Solution of Boundary-Value Problems, Studies in Mathematics and Its Applications 5, North-Holland, 983. MR72472 (85a:494) [7] O. Gonzalez, J.H. Maddocks and R.L. Pego, Multi-multiplier ambient-space formulations of constrained dynamical systems, with an application to elastodynamics, Arch. Rational Mech. Anal. 57 (2) MR8374 (22k:7427) [8] K. Gustafson, Theory and computation of periodic solutions of autonomous partial differential equation boundary value problems, with application to the driven cavity problem, Math. Comput. Modelling 22 (995), No. 9, MR365 (96h:7623) [9] J.G. Heywood and R. Rannacher, Finite element approximation of the nonstationary Navier- Stokes problem, I. Regularity of solutions and second-order error estimates for spatial discretization, SIAM J. Numer. Anal. 9 (982), MR6552 (83d:6526) [] P. Lin, A sequential regularization method for time-dependent incompressible Navier-Stokes equations, SIAM J. Numer. Anal. 34 (997), no.3, 5-7. MR453 (98b:65) [] P. Lin, Regularization Methods for Differential Equations and Their Numerical Solution, Ph.D. thesis, Department of Mathematics, University of British Columbia, 995. [2] P. Lin, Q. Guo and X. Chen, A fully explicit method and its parallelization for incompressible flow computation, Comput. Methods Appl. Mech. Engrg. 92/22-24 (23), [3] J. G. Liu, J. Liu, R. L. Pego, Stability and convergence of efficient Navier-Stokes solvers via a commutator estimate, Comm. Pure Appl. Math. 6 (27), pp MR [4] Lu Xiliang, Error Analysis for Incompressible Viscous Flow Based on a Sequential Regularization Formulation, Ph.D. thesis, Department of Mathematics, National University of Singapore, 26. [5] R. Rannacher, On the numerical solution of the incompressible Navier-Stokes equations, Z. Angew. Math. Mech. 73 (993), pp MR24222 (94h:7673) [6] J. Shen, Hopf bifurcation of the unsteady regularized driven cavity flows, J. Comp. Phy. 95 (99), [7] J. Shen, On error estimates of the penalty method for the unsteady Navier-Stokes equations, SIAM J. Numer. Anal. 32 (995), pp MR (96d:6553)
1. Introduction. Consider nonstationary Navier Stokes equations with a nonslip boundary condition
SIAM J. SCI. COMPUT. Vol. 3, No., pp. 398 49 c 8 Society for Industrial and Applied Mathematics LONG TIME NUMERICAL SOLUTION OF THE NAVIER STOKES EQUATIONS BASED ON A SEQUENTIAL REGULARIZATION FORMULATION
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationSECOND-ORDER FULLY DISCRETIZED PROJECTION METHOD FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 3 (06), pp. 9 0. ISSN: 07-669. URL: http://ejde.math.txstate.edu or
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES
MATHEMATICS OF COMPUTATION Volume, Number, Pages S 5-578XX- ERROR ESTIMATES FOR SEMI-DISCRETE GAUGE METHODS FOR THE NAVIER-STOKES EQUATIONS : FIRST-ORDER SCHEMES RICARDO H. NOCHETTO AND JAE-HONG PYO Abstract.
More informationGlobal regularity of a modified Navier-Stokes equation
Global regularity of a modified Navier-Stokes equation Tobias Grafke, Rainer Grauer and Thomas C. Sideris Institut für Theoretische Physik I, Ruhr-Universität Bochum, Germany Department of Mathematics,
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
More informationPROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
MATHEMATICS OF COMPUTATION Volume 71, Number 237, Pages 27 47 S 0025-5718(01)01313-8 Article electronically published on May 14, 2001 PROJECTION METHOD III: SPATIAL DISCRETIZATION ON THE STAGGERED GRID
More informationA SEQUENTIAL REGULARIZATION METHOD FOR TIME-DEPENDENT INCOMPRESSIBLE NAVIER STOKES EQUATIONS
SIAM J. NUMER. ANAL. c 1997 Society for Industrial and Applied Mathematics Vol. 34, No. 3, pp. 151 171, June 1997 1 A SEQUENTIAL REGULARIZATION METHOD FOR TIME-DEPENDENT INCOMPRESSIBLE NAVIER STOKES EQUATIONS
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More informationKonstantinos Chrysafinos 1 and L. Steven Hou Introduction
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS
More informationNew Discretizations of Turbulent Flow Problems
New Discretizations of Turbulent Flow Problems Carolina Cardoso Manica and Songul Kaya Merdan Abstract A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is
More informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationhal , version 6-26 Dec 2012
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS ABDEHAFID YOUNSI Abstract. In this paper, we give a new regularity criterion on the uniqueness results of weak solutions for the 3D Navier-Stokes equations
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More information1 Introduction. J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods
J.-L. GUERMOND and L. QUARTAPELLE 1 On incremental projection methods 1 Introduction Achieving high order time-accuracy in the approximation of the incompressible Navier Stokes equations by means of fractional-step
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationParameter Dependent Quasi-Linear Parabolic Equations
CADERNOS DE MATEMÁTICA 4, 39 33 October (23) ARTIGO NÚMERO SMA#79 Parameter Dependent Quasi-Linear Parabolic Equations Cláudia Buttarello Gentile Departamento de Matemática, Universidade Federal de São
More informationBLOW-UP AND EXTINCTION OF SOLUTIONS TO A FAST DIFFUSION EQUATION WITH HOMOGENEOUS NEUMANN BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 016 (016), No. 36, pp. 1 10. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu BLOW-UP AND EXTINCTION OF SOLUTIONS TO A FAST
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN
ON THE STRONG SOLUTIONS OF THE INHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATIONS IN A THIN DOMAIN XIAN LIAO Abstract. In this work we will show the global existence of the strong solutions of the inhomogeneous
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationA Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations
A Two-Grid Stabilization Method for Solving the Steady-State Navier-Stokes Equations Songul Kaya and Béatrice Rivière Abstract We formulate a subgrid eddy viscosity method for solving the steady-state
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
More informationGlobal well-posedness of the primitive equations of oceanic and atmospheric dynamics
Global well-posedness of the primitive equations of oceanic and atmospheric dynamics Jinkai Li Department of Mathematics The Chinese University of Hong Kong Dynamics of Small Scales in Fluids ICERM, Feb
More informationThe Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control
The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control A.Younes 1 A. Jarray 2 1 Faculté des Sciences de Tunis, Tunisie. e-mail :younesanis@yahoo.fr 2 Faculté des Sciences
More informationThe Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:
Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply
More informationOn Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations
On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations G. Seregin, V. Šverák Dedicated to Vsevolod Alexeevich Solonnikov Abstract We prove two sufficient conditions for local regularity
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationDECAY AND GROWTH FOR A NONLINEAR PARABOLIC DIFFERENCE EQUATION
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 133, Number 9, Pages 2613 2620 S 0002-9939(05)08052-4 Article electronically published on April 19, 2005 DECAY AND GROWTH FOR A NONLINEAR PARABOLIC
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationOn the uniqueness of heat flow of harmonic maps and hydrodynamic flow of nematic liquid crystals
On the uniqueness of heat flow of harmonic maps and hydrodynamic flow of nematic liquid crystals Fanghua Lin Changyou Wang Dedicated to Professor Roger Temam on the occasion of his 7th birthday Abstract
More informationCONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND JIAN-GUO LIU
MATHEMATICS OF COMPUTATION Volume 69, Number 232, Pages 135{1407 S 0025-571(00)0124-5 Article electronically published on March 24, 2000 CONVERGENCE OF GAUGE METHOD FOR INCOMPRESSIBLE FLOW CHENG WANG AND
More informationA posteriori regularity of the three-dimensional Navier-Stokes equations from numerical computations
A posteriori regularity of the three-dimensional Navier-Stokes equations from numerical computations Sergei I. Chernyshenko, Aeronautics and Astronautics, School of Engineering Sciences, University of
More informationCONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS
J. Korean Math. Soc. 34 (1997), No. 3, pp. 515 531 CONVERGENCE OF FINITE DIFFERENCE METHOD FOR THE GENERALIZED SOLUTIONS OF SOBOLEV EQUATIONS S. K. CHUNG, A.K.PANI AND M. G. PARK ABSTRACT. In this paper,
More informationGENERIC SOLVABILITY FOR THE 3-D NAVIER-STOKES EQUATIONS WITH NONREGULAR FORCE
Electronic Journal of Differential Equations, Vol. 2(2), No. 78, pp. 1 8. ISSN: 172-6691. UR: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) GENERIC SOVABIITY
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationSome recent results for two extended Navier-Stokes systems
Some recent results for two extended Navier-Stokes systems Jim Kelliher UC Riverside Joint work with Mihaela Ignatova (UCR), Gautam Iyer (Carnegie Mellon), Bob Pego (Carnegie Mellon), and Arghir Dani Zarnescu
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationPDE Constrained Optimization selected Proofs
PDE Constrained Optimization selected Proofs Jeff Snider jeff@snider.com George Mason University Department of Mathematical Sciences April, 214 Outline 1 Prelim 2 Thms 3.9 3.11 3 Thm 3.12 4 Thm 3.13 5
More informationProjection Method with Spatial Discretizations
Projection Method with Spatial Discretizations Weinan E School of Mathematics Institute for Advanced Study Princeton, NJ 08540 and Jian-Guo Liu Courant Institute of Mathematical Sciences New York, NY 10012
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.
ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS Emerson A. M. de Abreu Alexandre N. Carvalho Abstract Under fairly general conditions one can prove that
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationSTOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN
Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationPartial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces
Partial Differential Equations, nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Shih-Hsin Chen, Yung-Hsiang Huang 7.8.3 Abstract In these exercises always denote an open set of with smooth boundary. As
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationA Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains
A Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains Martin Costabel Abstract Let u be a vector field on a bounded Lipschitz domain in R 3, and let u together with its divergence
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationSPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM. Mirjana Stojanović University of Novi Sad, Yugoslavia
GLASNIK MATEMATIČKI Vol. 37(57(2002, 393 403 SPLINE COLLOCATION METHOD FOR SINGULAR PERTURBATION PROBLEM Mirjana Stojanović University of Novi Sad, Yugoslavia Abstract. We introduce piecewise interpolating
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationREMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID
REMARKS ON THE VANISHING OBSTACLE LIMIT FOR A 3D VISCOUS INCOMPRESSIBLE FLUID DRAGOŞ IFTIMIE AND JAMES P. KELLIHER Abstract. In [Math. Ann. 336 (2006), 449-489] the authors consider the two dimensional
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationOn Two Nonlinear Biharmonic Evolution Equations: Existence, Uniqueness and Stability
On Two Nonlinear Biharmonic Evolution Equations: Existence, Uniqueness and Stability Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We study the following two nonlinear evolution equations with a fourth
More informationThe Navier-Stokes Equations with Time Delay. Werner Varnhorn. Faculty of Mathematics University of Kassel, Germany
The Navier-Stokes Equations with Time Delay Werner Varnhorn Faculty of Mathematics University of Kassel, Germany AMS: 35 (A 35, D 5, K 55, Q 1), 65 M 1, 76 D 5 Abstract In the present paper we use a time
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationON SINGULAR PERTURBATION OF THE STOKES PROBLEM
NUMERICAL ANALYSIS AND MATHEMATICAL MODELLING BANACH CENTER PUBLICATIONS, VOLUME 9 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 994 ON SINGULAR PERTURBATION OF THE STOKES PROBLEM G. M.
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationPartial regularity for suitable weak solutions to Navier-Stokes equations
Partial regularity for suitable weak solutions to Navier-Stokes equations Yanqing Wang Capital Normal University Joint work with: Quansen Jiu, Gang Wu Contents 1 What is the partial regularity? 2 Review
More informationModels. A.Dunca and Y.Epshteyn
On the Adams-Stolz Deconvolution LES Models A.Dunca and Y.Epshteyn Abstract We consider a family of LES models with arbitrarily high consistency error O(δ 2N+2 ) for N = 1, 2, 3.. that are based on the
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationIntegro-differential equations: Regularity theory and Pohozaev identities
Integro-differential equations: Regularity theory and Pohozaev identities Xavier Ros Oton Departament Matemàtica Aplicada I, Universitat Politècnica de Catalunya PhD Thesis Advisor: Xavier Cabré Xavier
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationFourth Order Convergence of Compact Finite Difference Solver for 2D Incompressible Flow
Fourth Order Convergence of Compact Finite Difference Solver for D Incompressible Flow Cheng Wang 1 Institute for Scientific Computing and Applied Mathematics and Department of Mathematics, Indiana University,
More informationDecay in Time of Incompressible Flows
J. math. fluid mech. 5 (23) 231 244 1422-6928/3/3231-14 c 23 Birkhäuser Verlag, Basel DOI 1.17/s21-3-79-1 Journal of Mathematical Fluid Mechanics Decay in Time of Incompressible Flows Heinz-Otto Kreiss,
More informationUNIFORM DECAY OF SOLUTIONS FOR COUPLED VISCOELASTIC WAVE EQUATIONS
Electronic Journal of Differential Equations, Vol. 16 16, No. 7, pp. 1 11. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu UNIFORM DECAY OF SOLUTIONS
More informationu( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)
M ath 5 2 7 Fall 2 0 0 9 L ecture 4 ( S ep. 6, 2 0 0 9 ) Properties and Estimates of Laplace s and Poisson s Equations In our last lecture we derived the formulas for the solutions of Poisson s equation
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) Key Concepts Current Semester 1 / 25 Introduction The purpose of this section is to define
More informationON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM
Indian J. Pure Appl. Math., 44(3): 343-354, June 2013 c Indian National Science Academy ON THE SOLVABILITY OF A NONLINEAR PSEUDOPARABOLIC PROBLEM S. Mesloub and T. Hayat Mathematics Department, College
More informationA RECURRENCE THEOREM ON THE SOLUTIONS TO THE 2D EULER EQUATION
ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 001 006, March 2009 001 A RECURRENCE THEOREM ON THE SOLUTIONS TO THE 2D EULER EQUATION Y. CHARLES LI Abstract. In this article, I will prove
More informationVANISHING VISCOSITY LIMIT FOR THE 3D NONHOMOGENEOUS INCOMPRESSIBLE NAVIER-STOKES EQUATION WITH SPECIAL SLIP BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 169, pp. 1 13. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu VANISHING VISCOSITY LIMIT FOR THE 3D NONHOMOGENEOUS
More information1. Introduction Since the pioneering work by Leray [3] in 1934, there have been several studies on solutions of Navier-Stokes equations
Math. Res. Lett. 13 (6, no. 3, 455 461 c International Press 6 NAVIER-STOKES EQUATIONS IN ARBITRARY DOMAINS : THE FUJITA-KATO SCHEME Sylvie Monniaux Abstract. Navier-Stokes equations are investigated in
More information