Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathe
|
|
- Alexandrina Cummings
- 5 years ago
- Views:
Transcription
1 Convergence of A Galerkin Method for 2-D Discontinuous Euler Flows Jian-Guo Liu 1 Institute for Physical Science andtechnology and Department of Mathematics University of Maryland College Park, MD 2742 houping Xin 2 Courant Institute, New York University and IMS and Dept. of Math., The Chinese University of Hong Kong Shatin, N.T., Hong Kong Key Words: discontinuous Galerkin Method, discontinuous ows, incompressible Euler equations, strong convergence. AMS(MOS) subject classication: 65M6, 76M2. Abstract: We prove the convergence of a discontinuous Galerkin Method approximating the 2-D incompressible Euler equations with discontinuous initial vorticity:! 2 L 2 (). Furthermore, when! 2 L 1 (), the whole sequence is shown to be strongly convergent. This is the rst convergence result in numerical approximations of this general class of discontinuous ows. Some important ows such asvortex patches belong to this class. x1. Introduction. Numerical simulation of 2-D discontinuous incompressible ows is of considerable interests in both theoretical analysis and applications. It is believed that the Lagrangian methods such as vortex methods [5,9], or the ones based on contour 1 jliu@math.umd.edu. The research was supported in part by NSF grant DMS xinz@cims.nyu.edu. The research was supported in part by the heng Ge Ru Foundation, NSF Grant DMS , and DOE Grant De-FG2-88ER
2 dynamics [1,15] give preferable treatments for such ows especially for inviscid inter-facial ows. However, the convergence of such methods poses great diculties. Past eorts concentrate on either special ows (see [4,12,13]), or require heavy machinery (such as large derivation [14]) and yield much weak convergence results [2,3,14]. However, for more complicated ows (such as a ow mixing), such front-tracking methods are impossible to implement. Thus, grid-based methods such as nite dierence and nite elements are called for. Yet, the convergence of such methods is unknown as we know of [1, 12]. Recently, a discontinuous Galerkin method was proposed in [12] which has the main advantages that the energy is conserved even for upwind type numerical uxes, and amusingly, thenumerical enstropy is non-increasing in time. The main observation of this paper is to point out that the boundness of energy and enstropy are sucient condition for strong convergence for a class of discontinuous initial data! 2 L 2 including vortex patches. In particularly, our results imply that the discontinuous Galerkin methods in [12] do converge for such ows. x2. A Discontinuous Galerkin Method. The 2-D incompressible Euler equation in vorticity stream-function formulation reads: t! +(r? r)! = r? =(;@ x ) =! with no-ow boundary condition (2:1b) = and initial condition (2:1c)!j t= =! (x) 2 L 2 () where R 2 is a simply-connected domain with a C 2 boundary, or piecewise C 2 boundary with convex corners. Assume that is equipped with a quasi-uniform triangulation T h = fkg consisting of polygons K of maximum size (diameter) h. Denote h = [K. The vorticity! is approximated by! h in a discontinuous nite element space Vh k o = nv : v j K 2 P k (K) 8K 2T h, while the stream function is approximated by h in a 2
3 continuous one W h k = V h k \ C ( h ). Here P k (K) denotes the set of all polynomials of degree at most k on the cell K. In the following, we will also use the notations that hi stands for the standard integration over the whole domain h, while an integral over a sub-domain K is denoted by hi K. The semi-discrete discontinuous Galerkin method in [12] can be described by looking for! h 2 Vh k and h 2 W h k such that X (2:2a) h@ t! h v h i K ;h! h uh rv h i K + huh n c! h v ; h i e = e2@k 8v h 2 V k h (2:2b) ;hr h r' h i h = h! h ' h i h 8' h 2 W k h : where e is a cell boundary and n is its unit out-normal. We now explain the notations used in (2.2). First the velocity eld is given by uh = r? h. Note that even through both! h and test function v h may be discontinuous across the cell boundaries, yet the velocity eld possesses continuous normal component across each cell boundary due to the denition of the nite element space W h k. Thus the numerical ux in (2.2a) can be dened as follows: Denote by v ; h (v+ h ) the value of v h from the inside (outside) of the element K, then the upwind ux is set to be (2:3) c! h = 8 >< >:! ; h if uh n! + h if uh n < : It should be remarked that for smooth ows, we could use a central ux dened by (2:3 ) c! h = 1 2! h + +!; h : However, the up-wind uxes (2.3) are preferred since the main concerns here are discontinuous ows. The rst important property of this scheme is the conservation of (no numerical dissipation in) energy (2:4) kr? h( t)k L 2 (h ) = kr? h( )k L 2 (h ) for the upwind ux (2.3), which canbe veried directly by taking v h = h in (2.2a), summing up the resulting equations over all K in the triangulation, and using (2.2b) and the 3
4 continuity of the normal velocity across the cell boundaries. Next, taking v h =! h,integrating by parts for the second term in (2.2a), and summing up for all K and estimating the terms involving cell boundaries by using (2.3) and the continuity of the normal component of the velocity eld across the cell boundaries, one can show that the enstropy decays in the sense that (2:5) k! h ( t)k L 2 (h ) k! h ( )k L 2 (h ) k! k L 2 (h ) where the initial data! h ( ) is taken as the L 2 projection of! and hence is uniformly bounded in L 2. Furthermore, taking ' h = h in (2.2b), one derives the fact that (2:6) kr h k 2 L 2 ( h ) = ;h! h h i h : Our main observation in this paper is the fact that these three simple properties, (2.4)- (2.6), yield a strong convergence. To prove and state such a result, one needs some timeregularity estimate rst. x3. Time Regularity Estimate. In this section, we will prove the following lemma about the time regularity for the approximate solutions constructed by the discontinuous Galerkin method. Lemma 1. (time regularity) It holds that (3:1) k@ t! h k L 1 ([ T ) W ;2 q ()) + k@ t h k L 1 ([ T ) Lq ()) C for any 1 q<2 where! h and h denote respectively their natural extension or restriction from h to. Proof: We rst show that (3:2) k@ t! h k L 1 ([ T ) L 2 (h )) C h 2 For any smooth function v 2 C 1 () with zero extension to the outside, let v h be the L 2 projection of v in the space of V k h. Now taking v h as a test function in (2.2a), one gets (3:3) h@ t! h vi h = h@ t! h v h i h = X K h! h uh rv h i K ; X K 4 X e2@k huh n c! h v ; h i e I 1 + I 2
5 I 1 can be estimated directly as ji 1 jk! h k L 1 X k kuhk L 2 (K) krv h k L 2 (K) Ck! h k L 1 (h )kuhk L 2 (h )krv h k L 2 (h ) Using the inverse inequality, the L 2 estimate for uh and! h in (2.4) and (2.5), and the fact that v h is the L 2 projection of v, one has (3:4) ji 1 j C h 2 k! hk L 2 (h )kuhk L 2 (h )kv h k L 2 (h ) C h 2 kvk L 2 ( h ) Next, we proceed to estimate I 2. Since uh n is continuous and c! h take same value from both sides of a cell boundary, the contribution of the two sides will gives the jump of v h. Therefore, one can obtain (3:5) ji 2 j X K X e2@k hjuh n c! h (v + h ; v; h )ji e Thanks to the quasi-uniformly regularity in the triangulation, one can show (3:6) Hence, X X K e2@k k! h k L 2 (e)kv h k L 2 (e) C h X K k! h k L 2 (K)kv h k L 2 (K) (3:7) ji 2 j C h 2 kv hk L 2 (h ) C h 2 kvk L 2 ( h ) Combining (3.4) with (3.7) shows (3:8) jh@ t! h vi h j C h 2 kvk L 2 ( h ) Which yields the desired estimate (3.2). Since the natural extension t! h to from h gives equivalent norms. Hence we have also shown (3:2 ) k@ t! h k L 1 ([ T ) L 2 ()) C h 2 Next, let I h v be the piecewise linear interpolation of v in Vh k. Then one can decompose I 1 in (3.3) as (3:9) I 1 = X K h! h uh ri h vi K + X K h! h uh r(v h ;I h v)i K I 11 + I 12 5
6 I 11 is bounded by (3:1) ji 11 j X K k! h k L 2 (K) kuhk L 2 (K) kri h vk L 1 k! h k L 2 (h )kuhk L 2 (h )kri h vk L 1 (h ) Due to the inequality one obtains from (2.4), (2.5), and (3.1) that kri h vk L 1 (h ) Ckvk W 1 1 () (3:11) ji 11 jckvk W 1 1 () One can estimate I 12 similarly. Indeed, ji 12 j X K k! h k L 2 (K) kuhk L 2 (K) kr(v h ;I h v)k L 1 (h ) Ckr(v h ;I h v)k L 1 (h ) Using the inverse inequality kr(v h ;I h v)k L 1 (h ) C h 2 kv h ;I h vk L 2 (h ) and (3:12) kv h ;I h vk L 2 (h ) kv ;I h vk L 2 (h ) which holds true since v h is the L 2 projection of v, one gets (3:13) ji 12 j C h 2 kv ;I hvk L 2 (h ) This, together with the standard estimate for the interpolation, leads to (3:14) ji 12 jckvk H 2 () It follows from (3.11) and (3.14) that we have obtained an estimate on I 1 in (3.3) independent ofh. Now we can also derive anh-independent estimate on I 2 in (3.3) as follows. Noting that I h v is continuous at the cell boundary, one can insert it into the right hand side of (3.5) to obtain (3:15) ji 2 j X K X hjuh n c! h (v h ;I h v)ji e Ckuhk L 1 X X e2@k K e2@k 6 k! h k L 2 (e) kv h ;I h vk L 2 (e)
7 As in (3.6), one has (3:16) ji 2 j C h 2 kv h ;I h vk L 2 ( h ) C h 2 kv ;I hvk L 2 ( h ) Ckvk H 2 () Collecting all the estimates (3.11), (3.14), and (3.16), we arrive at (3:17) jh@ t! h vi h jc(kvk W 1 1 () + kvk H 2()) Ckvk W 2 p () for any p>2 : To complete the proof of the rst part of (3.1), we need to estimate the dierence for the left hand term between and h. First, jh@ t! h vi nh jk@ t! h k L 2 (nh )kvk L 1 (nh ) q jn h j Note the simple estimate jn h jch 2,and kvk L 1 (nh ) Ch 2 kvk W 1 1 () since v vanishes on the we can obtain from these and (3.2') that (3:18) jh@ t! h vi nh jchkvk W 1 1 () Chkvk W 2 q () for any 1 q<2. Consequently, (3:19) k@ t! h k L 1 ([ T ) W ;2 q ()) C for any 1 q<2 This proves the rst part of (3.1). Next, we prove the second part of inequality in (3.1). For f 2 L p () with zero extension to the outside, we let solve the following problem (3:2) ; = f in = and let h 2 W k h be the nite element solution: (3:21) hr h r' h i h = hf ' h i h for any ' h 2 W k h t h 2 W k h,wehave (3:22) h@ t h fi h = hr@ t h r h i h = ;h@ t! h h i h 7
8 where we have used (3.21) and (2.2b) after taking a time derivative. Rewrite (3.22) as (3:23) h@ t h fi h = ;h@ t! h i h ;h@ t! h ( h ; )i h Let p>2 be the dual number of q: 1=p +1=q = 1, one gets (3:24) jh@ t h fi h j k@ t! h k L 1 ([ T ) W ;2 q ())kk W 2 p () + k@ t! h k L 1 ([ T ) L 2 (h ))k h ; k L 2 (h ) Ckk W 2 p () + C h 2 k h ; k L 2 (h ) Since the domain is either C 2 or piecewise C 2 with convex corners, the following the elliptic regularity and L 2 estimate is true: (3:25) kk W 2 p () Ckfk L p () and (3:26) k h ; k L 2 (h ) Ch 2 kk H 2 () Thus, ( ) imply that Or jh@ t h fi h jckfk Lp () (3:27) k@ t h k L 1 ([ T ) Lq ( h )) C for any 1 q<2 Since natural extension of h to from h gives equivalent norms, we therefore have shown (3:27) k@ t h k L 1 ([ T ) Lq () C for any 1 q<2 This gives the second part of (3.1). The proof of the lemma is completed. x4. A Uniqueness Theorem. In this section, we generalize Yudovich's uniqueness theorem to show that weak solutions with corresponding vorticity inl 1 are unique in the wider class where the vorticities are in L 2. This generalization will be used in next section to obtain a stronger convergence theorem. The precise statement is given by: 8
9 Theorem 1: Assume that! 2 L 1. Then the weak solution (4:1)! 2 L 1 ([ T) L 1 ()) \ Lip ([ T) W ;2 r ()) to (2.1) is unique in the space L 1 ([ T) L q ()) \ Lip ([ T) W ;2 r ()) where q>4=3 and 1 r<2. Proof: Suppose that the initial boundary value problem (2.1) for the 2D Euler equations has two weak solutions with same initial vorticity! 2 L 1 and the following regularities: (4:2)! 1 2 L 1 ([ T) L 1 ()) \ Lip ([ T) W ;2 r ()) and (4:3)! 2 2 L 1 ([ T) L q ()) \ Lip ([ T) W ;2 r ()) q>4=3 and1 r<2. It suces to show that! 1 =! 2. Denote by 1 and 2 the stream functions in (2.1) corresponding to! 1 and! 2 respectively. Then by the elliptic regularity, we have (4:4) 1 2 L 1 ([ T) W 2 p ()) \ Lip ([ T) L r ()) and (4:5) 2 2 L 1 ([ T) W 2 q ()) \ Lip ([ T) L r ()) We denote the correspond velocities by u 1 = r? 1 and u 2 = r? 2. Then one can rewrite the Euler equations (2.1a) in a distribution sense: (4:6) r? (@ t u 1 + u 1 ru 1 )= 2D and (4:7) r? (@ t u 2 + u 2 ru 2 )= 2D Set (4:8) u = u 1 ; u 2 = 1 ; 2 9
10 and subtract (4.7) from (4.6) to get (4:9) r? t u + u 2 ru + uru 1 )= 2D It follows from (4.4) and (4.5) that (4:1) 2 L 1 ([ T) W 2 q ()) \ Lip ([ T) L r ()) Therefore, one can take as a test function in (4.9) to obtain integration by parts, and use fact that u = r? and vanishes on the boundary, wehave (4:11) u(@ t u + u 2 ru + uru 1 ) dx = where one has used the fact that u = r? and vanishes on the boundray. Dene (4:12) E(t) = juj 2 dx Due to the regularity assumption ( ), one can show that (4:13) Using the fact (4:14) d dt E(t) =2 u ut dx uu 2 ru dx = and the equation (4.9), we have (4:15) d dt E(t) 2 juj 2 jru 1 j dx Using the classical potential estimate (4:16) kru 1 ( t)k L p Cpk! 1 ( t)k L 1 for all 1 <p<1, where C is an constant independent onp, u 1 and! 1, together with the fact that (4:18) k! 1 ( t)k L 1 k! k L 1 1
11 one shows by the Holder inequality that d (p;1)=p dt E(t) juj dx 2p=(p;1) kru 1 k L p Cp juj 2p=(p;1) dx (p;1)=p where C is independent onp. The right hand side above can be further estimated as juj 2p=(p;1) dx = On the other hand, (juj 2 ) (p;2)=(p;1) (juj 4 ) 1=(p;1) = kuk 2(p;2)=(p;1) kuk L 4 Ckuk q=(4(q;1)) W kuk (3q;4)=(4(q;1)) 1 q L 2 L 2 kuk 4=(p;1) L 4 Since kuk W 1 q (4:19) Therefor (4:2) is bounded, we thus have shown that d E(t) CpE(t)1;(5q;4)=(4p(q;1)) dt d E(t) (5q;4)=(4p(q;1)) C dt Now one can conclude that E(t). Indeed, taking an interval [ T ] with the property that CT ( 1 ), one obtains from (4.2) and E() = that 2 (4:21) E(t) ( 1 2 )4p(q;1))=(5q;4)! as p tends to innity So E(t) for t 2 [ T ]. Repeating these arguments we conclude that E(t) =for all t<t. This completes the proof. x5. Main Convergence Theorem. Finally, we are able to state and prove our main convergence theorem. Theorem 2: Let R 2 be a simply connected domain with C 2 boundary (or piecewise smooth C 2 boundary with convex corners), and equipped with a quasi-uniform triangulation. Suppose that the initial vorticity! belongs to L 2 (). Let (! h h ) 2 Vh k W h k be the approximate solutions generated by the discontinuous Galerkin method (2.2). Then there exists a convergent subsequence of (! h h ) (for which we will still use the same notation for simplicity) such that (5:1)! h *! (star weakly) in L 1 ([ T) L 2 ()) \ Lip ([ T) W ;2 q ()) 11
12 for any 1 q<2and (5:2) h! (strongly) in L 2 ([ T) H 1 ()) and the limiting functions! have the properties that (5:3)! 2 L 1 ([ T) L 2 ()) \ Lip ([ T) W ;2 q ()) and (5:4) 2 L 2 ([ T) H 1 ( )) \ Lip ([ T) L q ()) for any 1 q<2, and for any 2 C 1 (( t) ) (5:5) T (! t +!r? r) dxdt = (5:6) =! in D In other words, (! ) isaweak solution to the Euler equation (2.1) with initial data!. Furthermore, if the initial data! 2 L 1 (), then the whole sequence of (! h h ) will converge to the unique solution of the Euler equations and the limiting vorticity! is bounded in L 1 ([ T) L 1 ()). Proof: As in the Lemma 1, we extend! h and h to form h naturally. First, (2.5) and (3.1) show that there is a subsequence of! h (for which we still use the same notation) such that (5:7)! h *! (star weakly) in L 1 ([ T) L 2 ()) \ Lip ([ T) W ;2 q ()) and the limiting function! satises (5.3). Next it follows from (2.5), (2.6), and the Poincare inequality that (5:8) k h k L 1 ([ T ) H 1 ) C: This, together with (3.1), shows that there there is a subsequence of h such that (5:9) h * (star weakly) in L 1 ([ T) H 1 ()) \ Lip ([ T) L q ()) 12
13 and 2 L 1 ([ T) H 1 ()) \ Lip ([ T) L q ()) Since the distance h is of order O(h 2 )and h vanishes h,wehave (5:1) k h k L 1 (@) Ch 2 kr h k L 1 (h ) Chkr h k L 2 (h ) Ch! Hence satises (5.4). We now show that h converges strongly. First, It follows from (5.8), (3.1) and the Lions-Aubin lemma that (5:11) h! strongly in L 2 ([ T) ) which, together with (2.6) and (5.7), yields (5:12) T T kr h k 2 L 2 dt = ; h! h h i dt!; T h! i dt: To obtain the strong convergence that (5:13) r h!r strongly in L 2 ([ T) ) it suces to show that T T kr h k 2 L 2 dt! kr k 2 L 2 dt which is a direct consequence of (5:14) ; T h! i dt = T kr k 2 L 2 dt: that will be veried below. Indeed, for any 2 C 1 (( T) ), taking ' h = I h in (2.2b) yields (5:15) ; T hr h ri h i dt = T h! h I h i dt where I h is the interpolation operator in W h k. Using the strong convergence [6], (5:16) I h! strongly in L 1 ([ T) L 2 ()) and L 1 ([ T) H 1 ()) 13
14 and weak convergences (5.7) and (5.9), one shows from (5.15) that (5:17) ; T hr ri dt = T h!i dt This immediately gives (5.14) since C 1 (( T) ) is dense in L 2 ([ T) H 1 ()). As a sequence of (5.13), one concludes that (5:18) uh! u strongly in L 2 ([ T) ) where u = r? Now we show that the limit functions (! ) are indeed a weak solution to the Euler equations. To this end, one can take v h to be I h for any 2 C 1 (( t) ) in (2.2a), sum over all the cells, and integrate in time to obtain (5:19) T (h@ t! h I h i;h! h uh ri h i) dt = where we have used the facts that the upwind uxes c! h are the same for the adjacent elements, both I h and the normal component of the velocity eld are continuous across the interior cell boundaries, and u n = on the exterior cell boundaries. Since is compactly supported, t I h = I t,wecanintegrate by parttoget (5:2) T (h! h I t i + h! h uh ri h i) dt = Using the weak convergence of! h in (5.7) and strong convergence of I h, one shows that T h! h I t i dt! T h!@ t i dt Similarly, it follows from the weak convergence of! h in (5.7), strong convergence of uh in (5.18) and strong convergence of r h that T h! h uh ri h i dt! T h! u ri dt: Hence (5:21) T h!@ t i dt + T h! u ri dt =: 14
15 This gives (5.5). Finally, (5.6) follows from (2.2b) by taking ' h = I h. Thus we have proved that (!) isaweak solution to the Euler equations (2.1). In the case that the initial data! 2 L 1 (), then the Cauchy problem for the Euler equations has a solutions! 2 L 1 ([ T) ), and from Theorem 1 we know that this solution is unique in the class of (5.3) and (5.4). Therefor every convergent subsequence has the same limit. As a consequence, the whole sequence of (! h h ) converges to the unique solution. This completes the proof of the theorem. Acknowledgment The authors would like to thank Professors Qiang Du and Jun ou for their careful reading and many useful comments on the rst draft of this paper, and thank Professor John Osborn for some helpful discussions. References [1] G.R. Baker and M.J. Shelley, On the connection between thin vortex layers and vortex sheets, J. FluidMech., 215 (199), 161{194. [2] J.T. Beale, The approximation of weak solutions to the 2-D Euler equations by vortex elements, Multidimensional Hyperbolic Problems and Computations, Edited by J. Glimm and A. Majda, IMA 29 (1991), 23{37. [3] Y. Brenier and G.H. Cottet, Convergence of particle methods with random rezoning for the 2-D Euler and Navier-Stokes equations, SIAM J. Numer. Anal., 32 (1995) 18{197. [4] R. Caisch and J. Lowengrub, Convergence of the vortex method for vortex sheets, SIAM J. Numer. Anal., 26 (1989), 16{18. 15
16 [5] A.J. Chorin, Numerical study of slightly viscous ow, J. Fluid Mech., 57 (1973), 785{ 796. [6] P. Ciarlet, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, [7] B. Cockburn and C.-W. Shu, TVB Runge-Kutta local projection discontinuous Galerkin nite element method for scalar conservation laws II: general framework, Math. Comp., 52 (1989) 411{435. [8] R. DiPerna and A. Majda, Concentrations in regularizations for 2-D incompressible ow, Comm. Pure Appl. Math., 4 (1987), 31{345. [9] R. Krasny, Computation of vortex sheet roll-up in the Tretz plane, J. Fluid Mech., 184 (1987), 123{155. [1] D. Levy and E. Tadmor, Non-oscillatory central schemes for the incompressible 2-D Euler equations, Math Research Letters, 4 (1997) 1{2. [11] J.-G. Liu and C.-W. Shu, A high order discontinuous Galerkin method for incompressible ows, preprint. (1998) [12] J.-G. Liu and. Xin, Convergence of vortex methods for weak solutions to the 2-D Euler equations with vortex sheets data, Comm. Pure Appl. Math., 48 (1995) 611{628. [13] S. Schochet, The weak vorticity formulation of the 2-D Euler equations and concentration-cancellation, Comm. Partial Dierential Equations 2 (1995) 177{114. [14] S. Schochet, The point-vortex method for periodic weak solutions of the 2-D Euler equations, Comm. Pure Appl. Math., 49 (1996) 911{965. [15] R. Temam, Navier-Stokes Equations, Theory and Numerical Analysis, North-Holland, New Yrok, (1979) 16
17 [16] N. abusky, M.H. Hughes, K.V. Roberts, Contour dynamics for Euler equations in two dimensions, J. Comput. Phys., 3 (1979) pp
Pointwise 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 information2 Ben Schweizer As Hocking describes in [5], a discussion focused on the following question: Can one assume that the dynamic contact angle is a consta
A well-posed model for dynamic contact angles Ben Schweizer February 27, 1998 Abstract We consider uid systems with a free boundary and with a point of contact of the free boundary with a solid wall. We
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 informationAn explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
Journal of Computational and Applied Mathematics 156 (2003) 319 343 www.elsevier.com/locate/cam An explicit nite element method for convection-dominated compressible viscous Stokes system with inow boundary
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 informationElement diameter free stability parameters. for stabilized methods applied to uids
Element diameter free stability parameters for stabilized methods applied to uids by Leopoldo P. Franca Laboratorio Nacional de Computac~ao Cientica (LNCC/CNPq) Rua Lauro Muller 455 22290 Rio de Janeiro,
More informationMODIFIED STREAMLINE DIFFUSION SCHEMES FOR CONVECTION-DIFFUSION PROBLEMS. YIN-TZER SHIH AND HOWARD C. ELMAN y
MODIFIED STREAMLINE DIFFUSION SCHEMES FOR CONVECTION-DIFFUSION PROBLEMS YIN-TZER SHIH AND HOWARD C. ELMAN y Abstract. We consider the design of robust and accurate nite element approximation methods for
More informationIt is known that Morley element is not C 0 element and it is divergent for Poisson equation (see [6]). When Morley element is applied to solve problem
Modied Morley Element Method for a ourth Order Elliptic Singular Perturbation Problem Λ Wang Ming LMAM, School of Mathematical Science, Peking University Jinchao u School of Mathematical Science, Peking
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 informationMath 497 R1 Winter 2018 Navier-Stokes Regularity
Math 497 R Winter 208 Navier-Stokes Regularity Lecture : Sobolev Spaces and Newtonian Potentials Xinwei Yu Jan. 0, 208 Based on..2 of []. Some ne properties of Sobolev spaces, and basics of Newtonian potentials.
More informationEntropic Schemes for Conservation Laws
CONSTRUCTVE FUNCTON THEORY, Varna 2002 (B. Bojanov, Ed.), DARBA, Sofia, 2002, pp. 1-6. Entropic Schemes for Conservation Laws Bojan Popov A new class of Godunov-type numerical methods (called here entropic)
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
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 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 informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationA mixed nite volume element method based on rectangular mesh for biharmonic equations
Journal of Computational and Applied Mathematics 7 () 7 3 www.elsevier.com/locate/cam A mixed nite volume element method based on rectangular mesh for biharmonic equations Tongke Wang College of Mathematical
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 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 informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationon! 0, 1 and 2 In the Zienkiewicz-Zhu SPR p 1 and p 2 are obtained by solving the locally discrete least-squares p
Analysis of a Class of Superconvergence Patch Recovery Techniques for Linear and Bilinear Finite Elements Bo Li Zhimin Zhang y Abstract Mathematical proofs are presented for the derivative superconvergence
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationOn the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data
On the limiting behaviour of regularizations of the Euler Equations with vortex sheet initial data Monika Nitsche Department of Mathematics and Statistics University of New Mexico Collaborators: Darryl
More informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationarxiv: v1 [math.ap] 5 Nov 2018
STRONG CONTINUITY FOR THE 2D EULER EQUATIONS GIANLUCA CRIPPA, ELIZAVETA SEMENOVA, AND STEFANO SPIRITO arxiv:1811.01553v1 [math.ap] 5 Nov 2018 Abstract. We prove two results of strong continuity with respect
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationRearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D
Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, U.K. R.J. Douglas, Isaac Newton
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationAbstract. A front tracking method is used to construct weak solutions to
A Front Tracking Method for Conservation Laws with Boundary Conditions K. Hvistendahl Karlsen, K.{A. Lie, and N. H. Risebro Abstract. A front tracking method is used to construct weak solutions to scalar
More informationWALTER LITTMAN BO LIU kind of problem is related to the deformation of a membrane with strings attached to several sides of the polygon. In other word
THE REGULARITY AND SINGULARITY OF SOLUTIONS OF CERTAIN ELLIPTIC PROBLEMS ON POLYGONAL DOMAINS Walter Littman Bo Liu School of Math., University of Minnesota, Minneapolis, MN55455 Abstract. The regularity
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More informationFormulation of the problem
TOPICAL PROBLEMS OF FLUID MECHANICS DOI: https://doi.org/.43/tpfm.27. NOTE ON THE PROBLEM OF DISSIPATIVE MEASURE-VALUED SOLUTIONS TO THE COMPRESSIBLE NON-NEWTONIAN SYSTEM H. Al Baba, 2, M. Caggio, B. Ducomet
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationThe Incompressible Navier Stokes Flow in Two Dimensions with Prescribed Vorticity
The Incompressible Navier Stokes Flow in Two Dimensions with Prescribed Vorticity Sagun Chanillo, Jean Van Schaftingen, and Po-Lam Yung To Dick Wheeden in friendship and appreciation Abstract We study
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 informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationConvergence of Vortex Methods for Weak Solutions to the 2-D Euler Equations with Vortex Sheet Data
Convergence of Vortex Methods for Weak Solutions to the 2-D Euler Equations with Vortex Sheet Data JIAN-GUO LIU Temple University AND ZHOUPING XIN Courant Institute Abstract We prove the convergence of
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 informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationMIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS
MIXED FINITE ELEMENT METHODS FOR PROBLEMS WITH ROBIN BOUNDARY CONDITIONS JUHO KÖNNÖ, DOMINIK SCHÖTZAU, AND ROLF STENBERG Abstract. We derive new a-priori and a-posteriori error estimates for mixed nite
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
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 informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
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 informationTHE LUBRICATION APPROXIMATION FOR THIN VISCOUS FILMS: REGULARITY AND LONG TIME BEHAVIOR OF WEAK SOLUTIONS A.L. BERTOZZI AND M. PUGH.
THE LUBRICATION APPROXIMATION FOR THIN VISCOUS FILMS: REGULARITY AND LONG TIME BEHAVIOR OF WEAK SOLUTIONS A.L. BERTOI AND M. PUGH April 1994 Abstract. We consider the fourth order degenerate diusion equation
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 informationGlobal and Localised A Posteriori Error Analysis in the maximum norm for finite element approximations of a convection-diffusion problem
FINITE ELEMENT CENTER PREPRINT 000 11 Global and Localised A Posteriori Error Analysis in the maximum norm for finite element approximations of a convection-diffusion problem Mats Boman Chalmers Finite
More informationMARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann
More informationStrong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators
Journal of Scientific Computing, Vol. 8, No., February 3 ( 3) Strong Stability Preserving Properties of Runge Kutta Time Discretization Methods for Linear Constant Coefficient Operators Sigal Gottlieb
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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More informationConvergence of The Multigrid Method With A Wavelet. Abstract. This new coarse grid operator is constructed using the wavelet
Convergence of The Multigrid Method With A Wavelet Coarse Grid Operator Bjorn Engquist Erding Luo y Abstract The convergence of the two-level multigrid method with a new coarse grid operator is studied.
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 informationContents Introduction Energy nalysis of Convection Equation 4 3 Semi-discrete nalysis 6 4 Fully Discrete nalysis 7 4. Two-stage time discretization...
Energy Stability nalysis of Multi-Step Methods on Unstructured Meshes by Michael Giles CFDL-TR-87- March 987 This research was supported by ir Force Oce of Scientic Research contract F4960-78-C-0084, supervised
More informationBOUNDARY CONTROLLABILITY IN PROBLEMS JOHN E. LAGNESE. linear hyperbolic systems having piecewise constant coecients in
ESAIM: Control, Optimisation and Calculus of Variations URL: http://www.emath.fr/cocv/,, BOUNDARY CONTROLLABILITY IN PROBLEMS OF TRANSMISSION FOR A CLASS OF SECOND ORDER HYPERBOLIC SYSTEMS JOHN E. LAGNESE
More informationk=6, t=100π, solid line: exact solution; dashed line / squares: numerical solution
DIFFERENT FORMULATIONS OF THE DISCONTINUOUS GALERKIN METHOD FOR THE VISCOUS TERMS CHI-WANG SHU y Abstract. Discontinuous Galerkin method is a nite element method using completely discontinuous piecewise
More informationLow Froude Number Limit of the Rotating Shallow Water and Euler Equations
Low Froude Number Limit of the Rotating Shallow Water and Euler Equations Kung-Chien Wu Department of Pure Mathematics and Mathematical Statistics University of Cambridge, Wilberforce Road Cambridge, CB3
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 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 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 informationEntropy stable schemes for degenerate convection-diffusion equations
Entropy stable schemes for degenerate convection-diffusion equations Silvia Jerez 1 Carlos Parés 2 ModCompShock, Paris 6-8 Decmber 216 1 CIMAT, Guanajuato, Mexico. 2 University of Malaga, Málaga, Spain.
More informationThe Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-19 Wien, Austria The Negative Discrete Spectrum of a Class of Two{Dimentional Schrodinger Operators with Magnetic
More informationAn Optimal Finite Element Mesh for Elastostatic Structural. Analysis Problems. P.K. Jimack. University of Leeds. Leeds LS2 9JT, UK.
An Optimal Finite Element Mesh for Elastostatic Structural Analysis Problems P.K. Jimack School of Computer Studies University of Leeds Leeds LS2 9JT, UK Abstract This paper investigates the adaptive solution
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 informationDecay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients
South Asian Journal of Mathematics 2012, Vol. 2 2): 148 153 www.sajm-online.com ISSN 2251-1512 RESEARCH ARTICLE Decay rate of the compressible Navier-Stokes equations with density-dependent viscosity coefficients
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 informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
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 informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationc1992 Society for Industrial and Applied Mathematics THE CONVERGENCE RATE OF APPROXIMATE SOLUTIONS FOR NONLINEAR SCALAR CONSERVATION LAWS
SIAM J. NUMER. ANAL. Vol. 9, No.6, pp. 505-59, December 99 c99 Society for Industrial and Applied Mathematics 00 THE CONVERGENCE RATE OF APPROXIMATE SOLUTIONS FOR NONLINEAR SCALAR CONSERVATION LAWS HAIM
More informationFINITE TIME BLOW-UP FOR A DYADIC MODEL OF THE EULER EQUATIONS
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 357, Number 2, Pages 695 708 S 0002-9947(04)03532-9 Article electronically published on March 12, 2004 FINITE TIME BLOW-UP FOR A DYADIC MODEL OF
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More information2 Jian Shen fundamental law in porous media ow. f is the source term. In general, the conductivity term is discontinuous for heterogeneous reservoirs
MIXED FINITE ELEMENT METHODS ON DISTORTED RECTANGULAR GRIDS 1 Jian Shen Communicated by R. E. Ewing, R. D. Lazarov and J. M. Thomas. Institute for Scientic Computation Texas A & M University College Station,
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 informationNumerical computation of three-dimensional incompressible Navier Stokes equations in primitive variable form by DQ method
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2003; 43:345 368 (DOI: 10.1002/d.566) Numerical computation of three-dimensional incompressible Navier Stokes equations
More informationANNALES DE L I. H. P., SECTION C
ANNALES DE L I. H. P., SECTION C KAISING TSO Remarks on critical exponents for hessian operators Annales de l I. H. P., section C, tome 7, n o 2 (1990), p. 113-122
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 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 informationOn the Regularity of Weak Solutions to the Magnetohydrodynamic Equations
On the Regularity of Weak Solutions to the Magnetohydrodynamic Equations Cheng HE (Institute of Applied Mathematics, Academy of Mathematics and System Science, Chinese Academy of Sciences, Beijing, 18,
More informationquantum semiconductor devices. The model is valid to all orders of h and to rst order in the classical potential energy. The smooth QHD equations have
Numerical Simulation of the Smooth Quantum Hydrodynamic Model for Semiconductor Devices Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-84 Abstract
More informationFinite Difference Schemes for Incompressible Flows in the Velocity Impulse Density Formulation
JOURNAL OF COMPUTATIONAL PHYSICS 130, 67 76 (1997) ARTICLE NO. CP965537 Finite Difference Schemes for Incompressible Flows in the Velocity Impulse Density Formulation Weinan E 1 Courant Institute of Mathematical
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 informationRening the submesh strategy in the two-level nite element method: application to the advection diusion equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2002; 39:161 187 (DOI: 10.1002/d.219) Rening the submesh strategy in the two-level nite element method: application to
More informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
More informationThe Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations
The Role of Convection and Nearly Singular Behavior of the 3D Navier-Stokes Equations Thomas Y. Hou Applied and Comput. Mathematics, Caltech PDE Conference in honor of Blake Temple, University of Michigan
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 informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationNonlinear Analysis. A regularity criterion for the 3D magneto-micropolar fluid equations in Triebel Lizorkin spaces
Nonlinear Analysis 74 (11) 5 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A regularity criterion for the 3D magneto-micropolar fluid equations
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More information