1. Introduction. Consider nonstationary Navier Stokes equations with a nonslip boundary condition
|
|
- Eugenia Phelps
- 5 years ago
- Views:
Transcription
1 SIAM J. SCI. COMPUT. Vol. 3, No., pp c 8 Society for Industrial and Applied Mathematics LONG TIME NUMERICAL SOLUTION OF THE NAVIER STOKES EQUATIONS BASED ON A SEQUENTIAL REGULARIZATION FORMULATION PING LIN, JIAN-GUO LIU, AND XILIANG LU Abstract. The sequential regularization method is a reformulation of the unsteady Navier Stokes equations from the viewpoint of constrained dynamical systems or the approximate Helmholtz Hodge projection. In this paper we study the long time behavior of the sequential regularization formulation. We give a uniform-in-time estimate between the solution of the reformulated system and that of the Navier Stokes equations. We also conduct an error analysis for the temporal discrete system and show that the error bound is independent of time. A couple of long time flow examples are computed to demonstrate this method. Key words. Navier Stokes equations, sequential regularization, iterative penalty method, long time solution, constrained dynamical system, approximate projection AMS subject classifications. 65M, 76D5 DOI..37/ Introduction. Consider nonstationary Navier Stokes equations with a nonslip boundary condition.).).3) u t +u )u = ν u p + f, divu =, u Ω =, u t= = u, in a smooth bounded domain Ω R over a long time. In estimating the error of numerical methods the Gronwall inequality is used in dealing with the nonlinear convection term and the pressure p term. Then a factor expct) C > ) usually appears in the error bound. This is a difficulty in that usual analysis techniques do not work in studying long time approximations of the Navier Stokes equations. Another computational difficulty of the Navier Stokes equations is to maintain the incompressibility constraint divu = over a long time. Various techniques have been developed to overcome some of the above difficulties. The dissipative term ν u is usually used to control the nonlinear convection term and the pressure term. Then the Reynolds number has to be restricted to be quite low in order to have enough dissipation. The projection method may be a widely used method in dealing with the incompressibility constraint. However, in the case of a no-slip boundary only normal velocity is enforced in the projection formulation. Received by the editors October 3, 6; accepted for publication in revised form) May 8, 8; published electronically October 5, 8. Department of Mathematics, National University of Singapore, Science Drive, Singapore 7543, Singapore. Current address: Division of Mathematics, University of Dundee, 3 Perth Road, Dundee DD 4HN, Scotland UK plin@maths.dundee.ac.uk). Department of Mathematics and Institute for Physical Sciences and Technology, University of Maryland, College Park, MD jliu@math.umd.edu). Department of Mathematics, National University of Singapore, Science Drive, Singapore 7543, Singapore. Current address: Johann Radon Institute for Computational and Applied Mathematics RICAM), Austrian Academy of Sciences, Altenbergerstrasse 69, A-44 Linz, Austria lu xiliang@hotmail.com). 398
2 LONG TIME NUMERICAL SOLUTION 399 The tangential velocity is not enforced in the formulation and may cause a slip error in a long time simulation. Another method is to write the Navier Stokes equations to a first-order system and then derive a new formulation based on the least squares minimization of its residual. However, quite a few new variables are introduced in the formulation which complicate the problem and increase the computational cost. Other formulations are available in dealing with the above difficulties, for example, a class of pseudocompressibility methods, which include the penalty method, the artificial compressibility method, etc. cf. [3,, 7,, 8]). But the small parameter associated with these methods introduces extra stiffness into the system, and the incompressibility constraint is not guaranteed over a long time. In [, 7, ], a sequential regularization method SRM) is introduced and analyzed to maintain the incompressibility constraint and to solve the unsteady Navier Stokes equations. In this paper we would like to study the long time behavior of the method because, in this method, many of the above issues are resolved or improved. The incompressibility constraint is uniformly under control in time, while there is no extra stiffness introduced due to the approximate Hodge projection to be explained in next section) or the built-in Baumgarte stabilization [, 7] for the incompressibility constraint. There is no slip error since the no-slip boundary condition is fully enforced. The pressure term is eliminated/replaced by a damping term, and then there is no need to control it in the analysis. The method or formulation reads: given an initial guess p and a positive constant α, for s =,, 3,..., find the solution u s,p s ) of the following system:.4).5).6) u s ) t +u s )u s = ν u s p s + f, divu s ) t + divαu s = ɛp s p s ), u s Ω =, u s t= = u. The method is developed from the viewpoint of constrained dynamical systems cf. [5, 7]). It is a combination of the penalty method cf. [4, 5]) and the Baumgarte time stabilization cf. [, 7]) for the constraint. It may be seen as an extension of the augmented Lagrangian method or the iterative penalty method for unsteady Navier Stokes equations cf. [3]). It may also be seen as an approximate Helmholtz Hodge projection, which will be explained later. Unlike the penalty method the parameter ɛ is not necessarily very small, and thus the reformulated system is more stable or less stiff see [7, ] or the convergence estimate 3.) later). It approximates the incompressibility constraint better than the penalty method. Roughly speaking we expect from.5) and 3.) that the divergence of the velocity u s is of Oɛ s ) and the bound may be independent of the time. From the regularity point of view, the method is a more natural formulation than the penalty formulation for the unsteady Navier Stokes equations see [] for more details). This method also decouples the velocity and the pressure. We can eliminate p s from the system, solve an equation only with unknown u s, and then recover p s from.5). When we eliminate p s from system.4).6), we obtain an equation only with the unknown u s. Let us omit the iteration index s and rewrite f p s as f; the equation can be written in the following form:.7).8) u t ɛ divu t + αu) ν u +u )u = f, u Ω =, u t= = a. In [], we consider.4).6) in a finite time interval [,T] and obtain the existence and regularity of the solution as well as the error estimates of the SRM.
3 4 PING LIN, JIAN-GUO LIU, AND XILIANG LU For some flow problems, we need to compute the solution over a very long time, for instance, the steady state solution of driven cavity flow. It is necessary to study the long time behavior of this formulation and extend the results in [] over a long time. This paper will be organized as follows. We will relate the SRM to an approximate Helmholtz Hodge projection and give another explanation/understanding of the method in section. Some preliminaries such as notations and assumptions will be given there as well. In sections 3 and 4, we establish energy inequalities and error estimates for the continuous equations and the temporal discrete equations, respectively. A couple of long time flow examples will be computed in section 5.. Approximate Helmholtz Hodge projection and the SRM. The SRM can be explained as an approximate Helmholtz Hodge projection. Let P : L Ω, R N ) H Ω)) denote the Helmholtz Hodge projection onto vector fields that are divergence-free and have zero normal component on the boundary. In prin- we may define Pu = u φ, where φ satisfies φ = divu, u n ciple = u n, and φdx =. One may apply P to both sides of.) and obtain Ω.) u t + Pu u f) =νp u. The dissipation in.) appears degenerate due to the fact that P annihilates gradients. So the analysis of.) is usually restricted to a divergence-free space. However, for the stability of numerical computation, some additional damping effects to the divergence need to be added [9]. In this formulation, solutions formally satisfy t divu) =. Consequently the divergence-free condition.) needs to be imposed only on the initial data. However, when a numerical perturbation applies, weak instability to the divergence-free condition may occur. In addition, the projection P may not be easily implemented numerically. To overcome these difficulties with the projection, we can do the following. We introduce an approximate/desingularized projection operator P ɛ by P ɛ u := u φ ɛ, where φ ɛ satisfies φ ɛ ɛφ ɛ = divu and φɛ n = u n. We can then verify.) I div ) P ɛ u = u, P ɛ u n Ω =. ɛ Thus we can denote Pɛ =I div ɛ ), which is a local operator asking no additional boundary conditions. Now, replacing P by P ɛ and taking the inverse operator yields an approximate Navier Stokes equation: or u t + P ɛ u u f) =νp ɛ u I div ) u t + u u = ν u + f. ɛ It is then the SRM.7) if adding a damping term α ɛ divu to the lefthand side of above equation accordingly to adding P α ɛ divu = to.)). It is.7), which will be discussed in this paper. From.5) the divergence of u will not drift too much away from zero as the time goes. Therefore, this formulation has better stability. The divergence of the projected solution is uniformly controlled due to the following identity obtained from.)): P ɛ u + ɛ divp ɛu + ɛ divp ɛu = u.
4 LONG TIME NUMERICAL SOLUTION 4 Here and in what follows we used notations and m as the norms for L and H m, respectively. Unlike the penalty method, there is no additional time step restriction. From the definition of P and P ɛ earlier, we can have P P ɛ )u Cɛ u, which justifies the approximation to the Navier Stokes equation as well. The approximate Helmholtz projection can also be viewed as a regularization to the Helmholtz projection. The Fourier symbol for P is represented by ˆP ξ = I ξ ξ, while the symbol for P ɛ is a desingularized one: ˆP ɛ = I ξ ξ ξ + ɛ. Simple computation shows that the symbol of its inverse is given by ˆP ɛ = I + ξ ξ ɛ, which again agrees with Pɛ = I div ɛ. It is known that the solution of the Navier Stokes equation is regular up to an arbitrary time interval in a two-dimensional domain cf. [6, 7]). More precisely, consider two conditions:.3).4) u N, u N, f,t) dt N, divu =, f,t) N, <t< f t dt N. If condition.3) holds, then there exists a constant M which depends on N such that.5) u + <t< u + p + u t ) dt M. If conditions.3) and.4) hold, then there exists a constant M which depends on N and N such that.6) u + u t + p ) M. <t< For the sequential regularized reformulation, we have similar results for.7).8); i.e., if we assume that.3) holds, then there exists a constant M such that.7) u + u t + <t< ɛ divu t + ) ɛ divu + u ds M. If.3) and.4) hold, then there exists a constant M such that.8) <t< ɛ divu t + ) ɛ divu + u + u t M. Since the pressure can be recovered by divu t +αdivu, the same bound for pressure follows. The key part in the proof is Lemma 3., which is similar to Lemma 4. in []
5 4 PING LIN, JIAN-GUO LIU, AND XILIANG LU but holds uniformly in time. Based on these results, we can extend the error estimation in [] to any time interval. Therefore, this sequential regularization formulation is suitable for long time computations. The numerical results in the last section will also demonstrate this claim. Next we introduce a few operators for the treatment of the nonlinear convection term. Let One can easily check Bu, v) =u )v, Bu, v) =Bu, v)+ divu)v, bu, v, w) =Bu, v), w), bu, v, w) = Bu, v), w). Therefore, bu, v, w) = {bu, v, w) bu, w, v)} u, v, w H Ω). bu, v, v) = u, v H Ω). For the trilinear form b or b), we can prove the following inequality by combination of integration of parts, the Hölder inequality, and the Sobolev inequality see, for instance, [7]) u v w.9) bu, v, w) C u v w u v w. Moreover, if Ω R, we have the following:.) bu, v, w) C u u v v w v H H, u H, w L, where the trilinear form b can be replaced by b. In what follows we will replace B by B. Hence.7).8) become.).) u t ɛ divu t + αu) ν u + Bu, u) =f, u Ω =, u t= = u. We remark here that we can also replace B by ω u based on the identity u )u = ω u + u, where ω is the vorticity u and the second term on the right-hand side may be absorbed into the pressure see [, 9]). Our analysis may be carried out for this form as well. 3. A priori estimations. We first assume that the solution of.).) exists and then establish an a priori estimation. We will give a basic energy estimation <t< u + u dt and a strong solution estimation <t< u + u dt in Theorem 3.3 and then obtain a higher-order regularity estimation <t< u in Theorem 3.4. Since we consider an arbitrary long time interval [,T], we need to improve the solution estimates in [] such that the choice of ɛ and the bound of solution do not depend on the time length T. The following two Lemmas 3. and 3. play an essential role in obtaining a time-independent estimation.
6 LONG TIME NUMERICAL SOLUTION 43 Lemma 3.. Let X be a Banach space, let f :, ) X be a vector value function which satisfies ft) dt < +, and let α be a positive constant; then we have [ t ] fs) e αs t) ds dt 3.) fs) ds, 3.) <t<t [ t α ] fs) e αs t) ds α ft). <t<t Proof. We apply the Cauchy Schwarz inequality and change the integration order to obtain [ t [ t fs) e ds] αs t) dt = e αt fs) e αs/ e ds] αs/ dt and α = α [ t fs) e ds] αs t) <t<t <t<t e αt [ t ][ t ] fs) e αs ds e αs ds dt t e αt fs) e αs dsdt fs) e αs e αt dtds α fs) ds <t<t <t<t <t<t e αt [ t α e αt α e αt α t t s ][ t ] fs) e αs ds e αs ds fs) e αs ds fs) <s<t = α ft). <t<t ) fs) e αs ds <s<t Lemma 3.. Define an operator Au = ɛ divu t +αu) ν u, where the initial data u is divergence-free. Then there exist a ɛ and a constant C, only depending on the domain Ω, α, and ν, such that ɛ ɛ, we have t ɛ divu t + ) t 3.3) ɛ divu + u dt C Au dt, 3.4) <t<t ɛ divu t + ɛ divu + u ) C <t<t Au. Proof. Define w = Au, p = ɛ divu t + αu), and g = divu. First we solve g from the ODE g t + αg = ɛp
7 44 PING LIN, JIAN-GUO LIU, AND XILIANG LU with initial condition g) = divu =. The solution is gt) = ɛ t ps)e αs t) ds. Let u,p) satisfy the nonhomogeneous Stokes equation ν u + p = w, divu = g. Then using the estimation for the nonhomogeneous Stokes equation cf. [6]), we have 3.5) u + p C w + g ), where the constant C only depends on the domain Ω and ν. We now work out a few inequalities between g and p by using Lemma 3.: [ t ] 3.6) g ɛ ps)e αs t) ds dt ɛ p dt, 3.7) <t<t gt) ɛ <t<t [ t α ] p e αs t) ds ɛ α pt). <t<t Now we are ready to prove the results of this lemma. Take the integration from to T for the inequality 3.5), and we have u + p ) dt C w + g ) dt C w ɛ dt + C α Take the for both sides of the inequality 3.5), and we have <t<t u + p ) C <t<t w + g ) dt Define ɛ = Then we notice that <t<t <t<t C <t<t α C, and then for all ɛ ɛ, we have ɛ divu ds = ɛ divu t ds = ɛ divu = ɛ divu t = <t<t w + C ɛ α u + p ) dt C w dt, u + p ) C w. <t<t <t<t <t<t ɛ g ds α ɛ g t ds = ɛ g α ɛ g t p ds, p dt. p. <t<t ɛ ɛ p + α g ds 4 p, <t<t <t<t p ds, ɛ ɛ p + α g 4 p. <t<t Combining the above inequalities together, we can easily conclude the lemma.
8 LONG TIME NUMERICAL SOLUTION 45 Remark 3.. Comparing with a similar lemma in [], the improvement here is that the choice of ɛ and the constant of the bound do not depend on the time interval T, which implies that we may fix ɛ in advance and then compute the solution to any arbitrarily large time. Next we will consider the modified SRMs.).). Theorem 3.3. For.).), we assume.3). Then there exist two constants M and M, which only depend on the domain and N, such that <t< u + ɛ divu ) + u + <t< ɛ divu + u ) ds M, u t + ɛ divu t + ɛ divu + u Proof. Multiplying u for both sides of.), we have d dt u + d ɛ dt divu + α ɛ divu + ν u =f, u). Using the Young inequality and integrating for both sides, we can obtain u + ɛ divu + t ) ɛ divu + u ds C u + t ) ds M. ) f ds, where the constant C only depends on the domain Ω, α, and ν. Since the above inequality holds true for any time t and C does not depend on t, it is indeed a uniform-in-time estimation. We denote M = C u + f ds) and simplify the inequality 3.8) u + ɛ ) divu + <t< ) ɛ divu + u ds M. Then let us fix a time interval,t), where T is arbitrarily chosen. Recall the definition of operator A in Lemma 3.. We multiply Au on both sides of.), ɛ divu t + α d ɛ dt divu + ν d dt u + Au + bu, u,au) =f,au). Since we consider the D case, we have bu, u,au) C u / u u / Au C δ u u Au + δ u ). Define yt) = ɛ divut) + ut) and kt) =C δ u u, and obtain ɛ divu t + dy dt + Au 4 f + δ u + kt)yt). Then we can obtain by the technique of the Gronwall lemma letting Kt) = t ks)ds) e Kt) ɛ divu t + de Kt) y + e Kt) Au δ u dt ) 4e Kt) f, yt )+ e KT ) Kt) ) Au δ u ) T dt e KT ) 4 f dt + y).
9 46 PING LIN, JIAN-GUO LIU, AND XILIANG LU Then noticing that e KT ) Kt) is monotonically decreasing, choosing δ to be a constant such that t Au ds t δ u ds is true for any t we can choose δ = C, where C is the constant in Lemma 3.), defining Gt) = t Au δ u )ds, and doing integration by parts, we have e KT ) Kt) Au δ u ) dt = e KT ) Kt) Gt) T +.5 Au δ u ) dt Au dt. Gt)e KT ) Kt) K t)dt Let M = e CM ft) dt + u) ), where C only depends on the choice of δ. Hence T 3.9) ɛ divu + u + Au dt M. Since the choice of δ is uniform for t, ) and due to 3.8), the above inequality is uniform-in-time as well. Then we multiply u t on both sides of.) and obtain u t +Au, u t )+ bu, u, u t )=f, u t ). By applying the Young inequality, Lemma 3., and inequality 3.9), we can prove this theorem. The error between the solutions of the sequential regularization formulation and the Navier Stokes equations can be estimated uniformly in time by using Theorem 3.3 and following the technique in [7, ]. We simply state the result here: 3.) u u s + <t< u t u s ) t + u u s + p p s ) )dt Mɛ) s p p ) dt, where M is a constant which depends on Ω, ν, α, u, f, and p. The proof is exactly the same as Theorem.5 in [], by noticing that all of the constants in the proof are uniform in time. Theorem 3.4. For.).), we assume.3).4). Then there exists a constant M 3, which depends on N and N, such that <t< ɛ divu t + ɛ divu + u + u t ) M 3. Proof. Multiplying u t on both sides of.), we have 3.) u t + α ) ɛ divu t ɛ divu, u t ν u, u t )+ bu, u, u t )=f, u t ). Choosing the time t =, applying the Young inequality, and noting that divu) =, u) H, and f L L ), we have 3.) u t ) + ɛ divu t) C.
10 LONG TIME NUMERICAL SOLUTION 47 Then differentiating.) with respect to time t, we have 3.3) u tt ɛ divu tt + αu t ) ν u t + Bu, u t )+ Bu t, u) =f t. Multiplying u t on both sides of 3.3) yields d dt u t + d ɛ dt divu t + α ɛ divu t +ν u t + bu t, u, u t )+ bu, u t, u t )=f t, u t ). Using but, u, u t ) C u u t + ν 4 u t, bu, ut, u t )=, f t, u t ) ν f t + ν 4 u t, and defining yt) = u t + ɛ divu t, we have dy dt + ɛ divu t + u t C 3 u yt)+ f t ). Applying the Gronwall inequality, we obtain t ) t yt)+ ɛ divu t + u t dt C 4 ) f t dt + y). Using inequality 3.) to control the right-hand side of the above inequality, we have the estimation for u t, i.e., 3.4) u t + t ) ɛ divu t + ɛ divu t + u t dt C 4 t, ). Then multiplying Au on both sides of.), we have 3.5) u t,au)+ Au + bu, u,au) =f,au). Since we have u t,au) u t + 4 Au, f,au) f + 4 Au, bu, u,au) C δ u 6 + δ u + 4 Au, 4 Au u t + C δ u 6 + δ u + f. We can choose δ = C, where C is the constant in Lemma 3., and then take the for both sides of the above inequality. Using Lemma 3. and noticing that u, f, and u t are bounded, we have 3.6) Au M 3, <t< where the constant M 3 is uniform-in-time. Finally applying Lemma 3., we obtain the estimation of each term in the lemma.
11 48 PING LIN, JIAN-GUO LIU, AND XILIANG LU 4. Time discretization. Let us consider the time discretization of.) by the semi-implicit scheme: u n+ u n u n+ ɛ div u n ) + αu n+ 4.) 4.) ν u n+ + B u n, u n+) = f n+, u Ω =, u = u, tn+ where t n = n and f n = t n fτ)dτ. Lemma 4.. Let a n and b n be two sequences in a Banach space which satisfy + αh)a n a n = hb n, a =, where α and h are two positive constants. We have the inequalities a n α n= N n= b n, a n n N α b n n N for all integers N>. Proof. First we express a n in terms of b i by induction: a n = By the Holder inequality Hence a n = h +αh b n + +αh a n = h + αh) n+ h + αh) n+ h α + αh) n+ a n n= = n= i= n i= n=i i= n i= n b i + αh) i i= n b i + αh) i i= n b i + αh) i. i= hb i + αh) n+ i. n i= + αh) i h α + αh) n+ b i + αh) i h α + αh) n+ b i + αh) i α b i n= = N α b i, i= ) n +αh
12 LONG TIME NUMERICAL SOLUTION 49 and a n h α + αh) n+ α b i. i N n i= + αh) i b i i n Then this lemma follows by taking the for both sides of the inequality. Lemma 4.. Assume that the time step is. Let N be any positive integer. Define Au n+ = ɛ div un+ u n + αu n+ ) ν u n+, n =,, 3,..., and u = u). Then there exist a constant ɛ, which depends on Ω, α and ν, such that, when ɛ ɛ, we have u n + ɛ divun + ɛ div un u n ) C Au n, n= n= u n + n N ɛ divun + ɛ div un u n ) C Au n, n N where the constant C does not depend on the choice of ɛ and N. Proof. Define w n = Au n, p n = ɛ div un u n + αu n ), and g n = divu n ; then and ν u n + p n = w n, divu n = g n, + α)g n g n = ɛp n. We have the inequality for the Stokes equations see 3.5)) 4.3) u n + p n C w n + g n ). Since u n satisfies the homogeneous Dirichlet boundary condition, we can choose the Banach space as H L with the norm f, and notice that u) is divergence-free; then we apply Lemma 3., Sum up inequality 4.3), n= n= g n ɛ α p n, n= g n ɛ n N α b n. n N u n + p n ) C w n + g n ) n= n= C w n + C ɛ ) α pn,
13 4 PING LIN, JIAN-GUO LIU, AND XILIANG LU and take the of 4.3), u n + p n ) C Let ɛ = n= n N C α C. Then ɛ ɛ, we have n= n N n N w n + g n ) w n + C ɛ α u n + p n ) C w n, n N n= u n + p n ) C w n. n N ) p n. n N Notice that ɛ divun = ɛ gn α p n, n= n= n= ) u n ɛ div u n = α ɛ divun + p n N 4 p n, n N ɛ n N ɛ divun = n= n N ) u n div u n = n N ɛ gn n N α ɛ divun + p n 4 n= α pn, n N p n. Combining all of the above inequalities together, we have the lemma. Lemma 4.3. For 4.), if we assume.3).4), then there exist constants M, which depends on N, and M, which depends on N and N, such that u n + u i + divu i + div ui u i n + u i u i ) M, u n M. n Proof. We will follow the proof of Lemmas 3.3 and 4.5 of []. First we notice that conditions.3).4) imply f n N, f n N, f n+ f n n< N. n= Multiplying u n+ on both sides of 4.) and using identity b a, b) = b a + b a ), we obtain 4.4) u n + ) ɛ divun + ɛ divun + u n C. Then we multiply Au n+ on both sides of 4.), define y n = ɛ divun + u n, k n = u n u n+, and we have y n+ y n + Au n+ δ u n+ + C δ n= k n y n + f n+ ).
14 LONG TIME NUMERICAL SOLUTION 4 We let δ be suitably small such that δ N u.5 N Aun holds for all N one choice is δ = C, where C is constant in Lemma 4.) and denote Kn =Π n +C δ k i ), a n =.5 Au n, b n =.5 Au n δ u ), and d n = C δ f n. Then we can rewrite the above inequality as y n+ + a n+ + b n+ Kn K n yn + d n+. By eliminating y, y,...,y n from above the inequalities, we obtain Using identity n+ n+ y n+ + b i K i = n+ K n n+ a i + K n i= b i + b i n+ K i Kn y + K n n K i ) i K i i= j= b j d i. + b K and noticing that K n is monotone increasing, i j= bj, we conclude that 4.5) ɛ divun + u + Au n M. Multiplying un+ u n on both sides of 4.), applying Lemma 4. and the above inequality, and noticing that the choice of N does not affect the constant on the righthand side, we have 4.6) u n + n u n + divu n + div un u n + u n u n ) M. Finally, the proof of the second inequality u M can be similarly done cf. Lemma 4.5 of []). In addition the existence and uniqueness of the solution to.).) can be shown after we obtain a priori estimates for the time discrete solution i.e., Lemma 4.3). We ignore the proof here since the steps are similar to those in [] for the finite time case. Define the error function e n = u n ut n ); then e n satisfies e n+ e n e n+ 4.7) ɛ grad div e n ) + αe n+ ν e n+ + Bu n, u n+ ) But n+ ), ut n+ )) = r n+, where initial data e = and the remainder term r n+ = tn+ [ η t n ) u tt η) ] t n ɛ grad divu ttη) dη. To obtain the full first-order error estimation, we need the remainder to be bounded; i.e., 4.8) r n C.
15 4 PING LIN, JIAN-GUO LIU, AND XILIANG LU Theorem 4.4. estimates of e n, Assume conditions.3).4) and 4.8); we can obtain error e n + n< ) ɛ diven + e n C. Proof. Since the proof is almost same as Theorem 4.6 of [], we give a sketch here. Due to the bilinear property of B, we have B u n, u n+) But n+ ), ut n+ )) = B u n u n+, u n+) + B u n+, e n+) + B e n+, ut n+ ) ). Then multiplying e n+ on both sides of 4.7), we obtain ) e n + n< ɛ diven + e n C. Then we multiply Ae n+ on both sides of 4.7) and notice ut) and u n to be uniformly bounded; we can conclude the result. Remark 4.. By the same proof of Lemma 4.3 in [], we can extend the result to a large time interval; i.e., with the same condition we will simplify the condition ) g = f) + ν u Bu, u ) H later) of Lemma 4.3 in [], we have 4.9) u t + <t< ɛ divu tt + u t + u tt dt C. From the above estimation and Remark 4. in [], we know that the condition ) and inequality 4.8) are true when the global compatibility condition holds cf. [6]). For simplicity, we just assume this global compatibility holds. Now we can compare the time discrete solution for the sequential regularization formulation with the solution for the Navier Stokes equations. We denote s as the SRM index and n as the time index see the detail of this implementation in section 5.). We combine inequality 3.) with Theorem 4.4 and obtain the following estimation. Theorem 4.5. Assume conditions.3).4) and the global compatibility condition at t =. Let u i s and u be the solutions to the time discrete sequential regularization reformulation and the original Navier Stokes equation, respectively. Then we have n u n s ut n ) + ɛ div u i s ut i ) ) + u i s ut i ) ) C + Mɛ ) ) s. Proof. We only need to check if all of the conditions for inequality 3.) and Theorem 4.4 are satisfied. Since the sequential regularization formulation replaces the external force f by f p s at the s + )st step of the SRM iteration, we need to check the following inequalities: 4.) 4.) p s + <t< p s dt C s, p s ) t dt C s,
16 LONG TIME NUMERICAL SOLUTION 43 where C s is a constant which does not depend on ɛ. The case of s = is not a problem since we can choose p freely. When s =,, 3,..., we have p s = p s ɛ divu s) t + αdivu s ) from.5). Hence we just need to show that 4.) 4.3) ɛ divu s) t + ɛ divu s dt C s, ) <t< ɛ divu s) t + ɛ divu s + ɛ divu s) tt + ɛ divu s) t dt C s. By mathematical induction, inequality 4.) is simply a corollary of Theorem 3.3, and inequality 4.3) is a corollary of Theorem 3.4 and estimation 4.9). Remark 4.. The analysis in this chapter is based on a semi-implicit scheme. In practice we are also interested in the fully explicit treatment of the nonlinear convection term. If u n is uniformly bounded for n, the theorem in this chapter is also true for the fully explicit treatment of the nonlinear convection term. Without any specific assumption, u n only remains bounded in a short time interval see Lemma 3.5 in []). But in practice, the velocity of many flow problems is stably away from blow up. And since the fully explicit treatment is easier to implement, we will use it in the next chapter. For other schemes, for example, the Crank Nicolson scheme, there is no essential difference for analysis, and the accepted order can be achieved by assuming enough regularity of the solution see the discussion in [, ]). 5. Numerical examples. In this section, we will present results of numerical experiments of long time simulation of two-dimensional viscous flows. 5.. Implementation. Let V h be the finite element space; the full discrete scheme of SRM can be represented as follows: given an initial guess p n, where n =,,,..., for s =,,..., we solve u n+ s,p n+ s ) from the system u n+ s u n ) s, v + div un+ s u n ) s, divv + α divu n+ s, divv ) ɛ ɛ +ν u n+ s, v ) + b u n s, u n s, v) = fs n+, v ) p n+ s, divv) v V h, p n+ s = div un+ s u n ) s + αdivu n+ s. ɛ The above scheme is as follows: starting from s = for each n, find u n and p n at time t n, and then repeat this procedure to find u n, p n, and so on. Since we need to store the values of u n and p n for all n time steps before the next SRM iteration, the storage may be very expensive since we compute the solution to a very long time. It is not what we do in practice. This difficulty can be overcome by an equivalent implementation. We make an observation: the solution u n+ s,p n+ s ) only depends on u n s and p n+ s. Now fix the number of iterations in advance according to the accuracy requirement numerical experiments indicate that s = 3 or 4 may be enough in most cases). Since the initial guess p n and the initial data u s are known, we can rearrange the order of the computation as follows. Starting from n =, find u,p ) by using information p and u, and find u,p ) by using information p and u,..., up to u s,p s). Save all the data u s,p s), and move to n =. Compute u,p ) by using p and u up
17 44 PING LIN, JIAN-GUO LIU, AND XILIANG LU Fig. 5.. Implementation instruction. to u s,p s). We can repeat this procedure at all time steps until the terminal time is reached cf. [8]). The procedure is illustrated in Figure 5.. The original SRM is to do the computation row by row, and the practical one is to do the computation column by column. It solves the storage problem since we only need to store u n,p n ),..., u n s,p n s ) with s = 3 or 4 for computations at the next time step n +. At every time step and SRM iteration, we need to solve a PDE of the form u +α divu ν u = g. ɛ Due to the noncommutativity of div and, the fast Poisson solver does not apply here. We will solve it by a direct solver. Since the equation does not depend on the time step n and SRM index s, we can discrete this PDE into matrix form and do the LU or the Choleski factorization at the initial time. This factorization at the later time would be the same as that of the initial time, and thus we simply do backward and forward substitutions in computations of all other time steps except the initial time. This implementation could save a lot of computational time. In our numerical examples, we will use the P finite element for the spatial discretization. Let Ω be a convex polygon and T h be a quasi-regular triangulation of Ω, where h represents the maximal diameters of the triangles in the triangulation. Define V h = { v C, v = on Ω, v Ti is a polynomial of degree }. 5.. Cavity flow. The first example is lid-driven cavity flow in a unit square. To find out the steady state solution, we use time-dependent Navier Stokes equations and compute them to a long time. The external force is f =, and the boundary condition is u =, ) at the top side of the square u = at the other three sides. We choose α =, =.5, ɛ =., s = 3, and the P element on a 3748 triangular mesh. We do the computation for Reynolds numbers Re = and Re = 5,
18 LONG TIME NUMERICAL SOLUTION e Fig. 5.. Streamline, Re = e 6 e 5 5e 5.. e e 5 e 5 e e 5..5 e 6 Fig Left eddy, Re =. Fig Right eddy, Re =. respectively. We will compare the computational results with the benchmark solution in [4] as well. Re =. The solution approaches the steady state around T =. Figure 5. is the streamline profile. There are three eddies in the streamline figure. We zoom in on the left and the right eddies in Figures 5.3 and 5.4 and depict the value of the velocity u at vertical line x =.5 and v at horizontal line y =.5 in Figures 5.5 and 5.6, respectively. Re = 5. When Reynolds number is larger, we must compute the solution to a longer time T to ensure that the flow approaches steady state. Meanwhile we need to choose a smaller time step to maintain the stability. It increases the computational cost. We can use the steady state solution at a smaller Reynolds number as our initial condition to reduce the computational time. Figure 5.7 is the streamline profile. There are four eddies in the streamline figure. Figures 5.8, 5.9, and 5. are zoomed figures of the left bottom eddy, right bottom eddy, and left top eddy, respectively. Figures 5. and 5. depict the value of u at vertical line x =.5 and v at horizontal line y =.5, respectively Flow past a cylinder. The second example is flow past a cylinder. The domain is unbounded. It is not possible to do numerical computation in an infinite domain, so we must do a truncation. The computational domain is a rectangle excluding a disk the cross section of the cylinder) inside. The x coordinate of the rectangle is from to 5, the y coordinate is from 4 to 4, and the disk is centered at 4, ) with diameter. The boundary condition can be represented as follows: the inflow at the left side has a constant velocity which is parallel to the x-axis u =, )); the flow at the upper and lower boundaries is parallel to the x-axis, and its velocity
19 46 PING LIN, JIAN-GUO LIU, AND XILIANG LU Re= Present work) Re= Ghia et al., 98).. Re= Present work) Re= Ghia et al., 98) Fig u.5,y), Re =. Fig vx,.5), Re =..9 Eddy TL Eddy BL Eddy BR e 5 e 6 e Fig Streamline, Re = 5. Fig Right eddy, Re = e 5 e e 6 5e 5 e Fig Top eddy, Re = 5. Fig. 5.. Left eddy, Re = Re=5 Present work) Re=5 Ghia al et., 98) Re=5 Present work) Re=5 Ghia al et., 98) Fig. 5.. u.5,y), Re = 5. Fig. 5.. vx,.5), Re = 5.
20 LONG TIME NUMERICAL SOLUTION 47 4 solid wall, u,v)=,) 4 4 r=.5 y u/ n=, v= outflow, u/ n=, v/ n= inflow, u,v)=,) x u/ n=, v= 4 5 Fig Flow past disk: geometry and boundary conditions Fig Re = 4, T = Fig Re =, T = 35. does not change along the y-direction u n =,v = ); the boundary of the disk is a solid wall with a no-slip boundary condition u = ); the outflow at the right side is imposed by a natural boundary condition homogeneous Neumann boundary condition, u n = v n = ). Figure 5.3 is a diagram for the computational domain and the boundary conditions. It is well known that the Reynolds number plays a crucial role in this problem. We will compute the solution with two different Reynolds numbers Re = 4 and Re = ). Re = 4. In this case, the flow approaches steady state after a certain time. Figure 5.4 depicts the streamline after the flow reaches its steady state. There are two symmetric eddies behind the disk. Re =. When Re =, the flow does not have a steady state. Figures depict the streamlines at different times. From these figures, we can see how the flow separates. All of our computational results are pretty close to existing results.
21 48 PING LIN, JIAN-GUO LIU, AND XILIANG LU Fig Re =, T = Fig Re =, T = 6. REFERENCES [] U. Ascher and P. Lin, Sequential regularization methods for higher index DAEs with constraint singularities: The linear index- case, SIAM J. Numer. Anal., ), pp [] A. J. Chorin, Numerical solution of the Navier Stokes equations, Math. Comp., 968), pp [3] M. Fortin and R. Glowinski, Augmented Lagrangian Methods: Applications to the Numerical Solution of Boundary-Value Problems, Stud. Math. Appl. 5, North-Holland, Amsterdam, 983. [4] U. Ghia, K. N. Ghia, and C. T. Shin, High-Re solutions for incompressible flow using the Navier Stokes equations and a mutigrid method, J. Comput. Phys., 48 98), pp [5] O. Gonzalez, J. H. Maddocks, and R. L. Pego, Multi-multiplier ambient-space formulations of constrained dynamical systems, with an application to elastodynamics, Arch. Ration. Mech. Anal., 57 ), pp [6] 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 98), pp [7] P. Lin, A sequential regularization method for time-dependent incompressible Navier Stokes equations, SIAM J. Numer. Anal., ), pp [8] P. Lin, X. Chen, and M. Ong, Finite element methods based on a new formulation for the non-stationary incompressible Navier Stokes equations, Int. J. Numer. Methods Fluids, 46 4), pp [9] J.-G. Liu, J. Liu, and R. Pego, Stability and convergence of efficient Navier Stokes solvers via a commutator estimate, Comm. Pure Appl. Math., 6 7), pp [] J.-G. Liu and W. C. Wang, Energy and helicity preserving schemes for hydro- and magnetohydro-dynamics flows with symmetry, J. Comput. Phys., 4), pp [] X. Lu, Error Analysis for Incompressible Viscous Flow Based on a Sequential Regularization Formulation, Ph.D. thesis, Department of Mathematics, National University of Singapore, Singapore, 6. [] X. Lu, P. Lin, and J.-G. Liu, Analysis of a sequential regularization method for the unsteady Navier Stokes equations, Math. Comp., 77 8), pp [3] A. Prohl, Projection and Quasi-Compressibility Methods for Solving the Incompressible Navier Stokes Equations, Teubner, Stuttgart, 997.
22 LONG TIME NUMERICAL SOLUTION 49 [4] R. Rannacher, On the numerical solution of the incompressible Navier Stokes equations, ZAMM Z. Angew. Math. Mech., ), pp [5] J. Shen, On error estimates of the penalty method for unsteady Navier Stokes equations, SIAM J. Numer. Anal., 3 995), pp [6] R. Temam, Navier Stokes Equations, North Holland, Amsterdam, 977. [7] R. Temam, Navier Stokes Equations and Nonlinear Functional Analysis, nd ed., CBMS-NSF Regional Conf. Ser. Appl. Math., Philadelphia, 995. [8] R. Temam, Une méthode d approximation de la solution des équations de Navier Stokes, Bull. Soc. Math. France, ), pp [9] T. E. Tezduyar, J. Liou, and D. K. Ganjoo, Incompressible flow computations based on the vorticity-stream function and velocity-pressure formulations, Comput. & Structures, 35 99), pp
ANALYSIS OF A SEQUENTIAL REGULARIZATION METHOD FOR THE UNSTEADY NAVIER-STOKES EQUATIONS
MATHEMATICS OF COMPUTATION Volume 77, Number 263, July 28, Pages 467 494 S 25-578(8)287-5 Article electronically published on February, 28 ANALYSIS OF A SEQUENTIAL REGULARIZATION METHOD FOR THE UNSTEADY
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 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 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 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 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 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 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 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 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 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 informationFast Solvers for Unsteady Incompressible Flow
ICFD25 Workshop p. 1/39 Fast Solvers for Unsteady Incompressible Flow David Silvester University of Manchester http://www.maths.manchester.ac.uk/~djs/ ICFD25 Workshop p. 2/39 Outline Semi-Discrete Navier-Stokes
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 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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
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 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 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 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 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 informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationLocal null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain
Local null controllability of the N-dimensional Navier-Stokes system with N-1 scalar controls in an arbitrary control domain Nicolás Carreño Université Pierre et Marie Curie-Paris 6 UMR 7598 Laboratoire
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 informationAnalysis and Computation of Navier-Stokes equation in bounded domains
1 Analysis and Computation of Navier-Stokes equation in bounded domains Jian-Guo Liu Department of Physics and Department of Mathematics Duke University with Jie Liu (National University of Singapore)
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 information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
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 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 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 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 informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
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 informationDIRECTION OF VORTICITY AND A REFINED BLOW-UP CRITERION FOR THE NAVIER-STOKES EQUATIONS WITH FRACTIONAL LAPLACIAN
DIRECTION OF VORTICITY AND A REFINED BLOW-UP CRITERION FOR THE NAVIER-STOKES EQUATIONS WITH FRACTIONAL LAPLACIAN KENGO NAKAI Abstract. We give a refined blow-up criterion for solutions of the D Navier-
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 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 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 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 Review of Solution Techniques for Unsteady Incompressible Flow
Zeist 2009 p. 1/57 A Review of Solution Techniques for Unsteady Incompressible Flow David Silvester School of Mathematics University of Manchester Zeist 2009 p. 2/57 PDEs Review : 1966 1999 Update : 2000
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
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 informationON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS
Submitted to Math. Comp. ON THE ERROR ESTIMATES FOR THE ROTATIONAL PRESSURE-CORRECTION PROJECTION METHODS J.L. GUERMOND 1 AND JIE SHEN 2 Abstract. In this paper we study the rotational form of the pressure-correction
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 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 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 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 informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationPost-processing of solutions of incompressible Navier Stokes equations on rotating spheres
ANZIAM J. 50 (CTAC2008) pp.c90 C106, 2008 C90 Post-processing of solutions of incompressible Navier Stokes equations on rotating spheres M. Ganesh 1 Q. T. Le Gia 2 (Received 14 August 2008; revised 03
More informationVector Penalty-Projection Methods for the Solution of Unsteady Incompressible Flows
Author manuscript, published in "5th International Symposium on Finite Volumes for Complex Applications, Aussois : France 2008)" Vector Penalty-Projection Methods for the Solution of Unsteady Incompressible
More informationOn a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions
Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability
More informationFinite Element Decompositions for Stable Time Integration of Flow Equations
MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationImplementation of 3D Incompressible N-S Equations. Mikhail Sekachev
Implementation of 3D Incompressible N-S Equations Mikhail Sekachev Navier-Stokes Equations The motion of a viscous incompressible fluid is governed by the Navier-Stokes equations u + t u = ( u ) 0 Quick
More informationOn the local existence for an active scalar equation in critical regularity setting
On the local existence for an active scalar equation in critical regularity setting Walter Rusin Department of Mathematics, Oklahoma State University, Stillwater, OK 7478 Fei Wang Department of Mathematics,
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 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 informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationPartially Penalized Immersed Finite Element Methods for Parabolic Interface Problems
Partially Penalized Immersed Finite Element Methods for Parabolic Interface Problems Tao Lin, Qing Yang and Xu Zhang Abstract We present partially penalized immersed finite element methods for solving
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationOptimization and control of a separated boundary-layer flow
Optimization and control of a separated boundary-layer flow Journal: 2011 Hawaii Summer Conferences Manuscript ID: Draft lumeetingid: 2225 Date Submitted by the Author: n/a Contact Author: PASSAGGIA, Pierre-Yves
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
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 informationA High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations
A High-Order Discontinuous Galerkin Method for the Unsteady Incompressible Navier-Stokes Equations Khosro Shahbazi 1, Paul F. Fischer 2 and C. Ross Ethier 1 1 University of Toronto and 2 Argonne National
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
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 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 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 informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationThe 2D Magnetohydrodynamic Equations with Partial Dissipation. Oklahoma State University
The 2D Magnetohydrodynamic Equations with Partial Dissipation Jiahong Wu Oklahoma State University IPAM Workshop Mathematical Analysis of Turbulence IPAM, UCLA, September 29-October 3, 2014 1 / 112 Outline
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 informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
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 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 informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
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 informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
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 informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationBlow-up of solutions for the sixth-order thin film equation with positive initial energy
PRAMANA c Indian Academy of Sciences Vol. 85, No. 4 journal of October 05 physics pp. 577 58 Blow-up of solutions for the sixth-order thin film equation with positive initial energy WENJUN LIU and KEWANG
More informationFEniCS Course. Lecture 6: Incompressible Navier Stokes. Contributors Anders Logg André Massing
FEniCS Course Lecture 6: Incompressible Navier Stokes Contributors Anders Logg André Massing 1 / 11 The incompressible Navier Stokes equations u + u u ν u + p = f in Ω (0, T ] u = 0 in Ω (0, T ] u = g
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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
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 informationTHE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS
ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 007 014, March 2009 002 THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS Y. CHARLES LI Abstract. Nadirashvili presented a
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationFine Grid Numerical Solutions of Triangular Cavity Flow
Published in : The European Physical Journal - Applied Physics (2007) Eur. Phys. J. Appl. Phys. 2007; Vol 38: pp 97-105 Fine Grid Numerical Solutions of Triangular Cavity Flow Ercan Erturk 1 and Orhan
More informationJMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman
JMBC Computational Fluid Dynamics I Exercises by A.E.P. Veldman The exercises will be carried out on PC s in the practicum rooms. Several (Matlab and Fortran) files are required. How these can be obtained
More informationWeak Solutions to Nonlinear Parabolic Problems with Variable Exponent
International Journal of Mathematical Analysis Vol. 1, 216, no. 12, 553-564 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.216.6223 Weak Solutions to Nonlinear Parabolic Problems with Variable
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationVANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE
VANISHING VISCOSITY IN THE PLANE FOR VORTICITY IN BORDERLINE SPACES OF BESOV TYPE ELAINE COZZI AND JAMES P. KELLIHER Abstract. The existence and uniqueness of solutions to the Euler equations for initial
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationVector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy
Vector penalty-projection methods for outflow boundary conditions with optimal second-order accuracy Philippe Angot, Rima Cheaytou To cite this version: Philippe Angot, Rima Cheaytou. Vector penalty-projection
More informationSolving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Solving Distributed Optimal Control Problems for the Unsteady Burgers Equation in COMSOL Multiphysics Fikriye Yılmaz 1, Bülent Karasözen
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 information