Generalized Newtonian Fluid Flow through a Porous Medium
|
|
- Antony Barrett
- 5 years ago
- Views:
Transcription
1 Generalized Newtonian Fluid Flow through a Porous Medium V.J. Ervin Hyesuk Lee A.J. Salgado April 24, 204 Abstract We present a model for generalized Newtonian fluid flow through a porous medium. In the model the dependence of the fluid viscosity on the velocity is replaced by a dependence on a smoothed (locally averaged) velocity. With appropriate assumptions on the smoothed velocity, existence of a solution to the model is shown. Two examples of smoothing operators are presented in the appendix. A numerical approximation scheme is presented and an a priori error estimate derived. A numerical example is given illustrating the approximation scheme and the a priori error estimate. Key words. Darcy equation, Generalized Newtonian fluid AMS Mathematics subject classifications. 65N30, 75D03, 76A05, 76M0 Introduction Of interest in this article is the modeling and approximation of generalized Newtonian fluid flow through a porous medium. Darcy s modeling equations for a steady-state fluid flow through a porous medium,, are ν eff K u + p = 0, in, (.) u = 0, in. (.2) where u and p denote the velocity and pressure of the fluid, respectively. K(x) in (.) represents the permeability of the medium at x, which is assumed to be a symmetric, positive definite tensor. As our investigations are not concerned with K, we assume that K is of the form k(x)i where k(x) is a Lipschitz continuous, positive, bounded and bounded away from zero, scalar function. ν eff in (.) represents the effective viscosity of the fluid. Department of Mathematical Sciences, Clemson University, Clemson, SC 29634, USA. (vjervin@clemson.edu) Department of Mathematical Sciences, Clemson University, Clemson, SC , U.S.A. (hklee@mail.clemson.edu) Partially supported by the NSF under grant no. DMS Department of Mathematics, University of Tennessee, Knoxville, TN 37996, U.S.A. (asalgad@utk.edu)
2 In the case of a Newtonian fluid we have that ν eff is a positive constant. For a generalized Newtonian fluid ν eff is a function of u. Two such examples are Power Law Model: ν eff ( u ) = c ν u r 2, Cross Model: ν eff ( u ) = ν + ν 0 ν, (.3) + c ν u 2 r where c ν, ν 0, ν and r are fluid dependent constants. For shear thinning fluids < r < 2. (In modeling the viscosity of shear thinning fluids the Power Law model suffers the criticism that as u 0 ν eff.) For the case of a Newtonian fluid (.), (.2) are well studied. The two standard approaches in analyzing (.), (.2) are: (i) study (.), (.2) as a mixed formulation problem for u and p (either (u, p) H div () L 2 (), or (u, p) L 2 () H ()), or (ii) use (.2) to eliminate u in (.) to obtain a generalized Laplace s equation for p. For generalized Newtonian fluids, with ν eff = ν eff ( u ), assumptions are required on ν eff in order to establish existence and uniqueness of solutions. Typical assumptions are uniform continuity of ν eff ( u )u and strong monotonicity of ν eff ( u ) [7, 8, 0], i.e., there exists C > 0 such that ν eff ( u )u ν eff ( v )v C u v, u, v IR d, (.4) (ν eff ( u )u ν eff ( v )v) (u v) C (u v) (u v), u, v IR d. (.5) A more general setting where the fluid rheology is defined implicitly has been analyzed in [5, 6]. The case where the fluid viscosity depends on the shear rate and pressure has been studied in [3, 2]. For both of these cases additional structure beyond (.4) and (.5) is required in order to establish existence and uniqueness of a solution. A nonlinear Darcy fluid flow problem, with a permeability dependent upon the pressure was investigated by Azaïez, Ben Belgacem, Bernardi, and Chorfi [2], and Girault, Murat, and Salgado []. For a Lipschitz continuous permeability function, bounded above and bounded away from zero, existence of a solution (u, p) L 2 () H () was established. Important in handling the nonlinear permeability function, in establishing existence of a solution, was the property that p H (). In [2] the authors also investigated a spectral numerical approximation scheme for the nonlinear Darcy problem, assuming an axisymmetric domain. A convergence analysis for the finite element discretization of that problem was given in []. Our interest in this paper is in relaxing the assumptions (.4) and (.5). Specifically, our interest is assuming that ν eff ( ) is only Lipschitz continuous and both bounded above and bounded away from zero. However, relaxing the conditions (.4) and (.5) requires us to make an additional assumption regarding the argument of ν eff ( ). In order to obtain a modeling system of equations for which a solution can be shown to exist, we replace u in ν eff ( u ) by a smoothed velocity, u s. The approach of regularizing the model with the introduction of u s is, in part, motivated by the fact that the Darcy fluid flow equations can be derived by averaging, e.g. volume averaging [6], homogenization [], or mixture theory [4]. Presented in the Appendix are two smoothing operators for u. One is a local averaging operator, whereby u s (x) is obtained by averaging u in a neighborhood of x. The second smoothing operator, which is nonlocal, computes u s (x) using a differential filter applied to u. That is, u s is given by the solution to an elliptic differential equation whose right hand side is u. For establishing the existence 2
3 of a solution to (.)-(.2), the key property of the smoothing operators is that they transform a weakly convergent sequence in L 2 () into a sequence which converges strongly in L (). For the mathematical analysis of this problem it is convenient to have homogeneous boundary conditions. This is achieved by introducing a suitable change of variables. For example, assuming = Γ in Γ Γ out, in the case the specified boundary conditions are u ( n) = g in on Γ in, u n = 0 on Γ, p = p out on Γ out, we introduce functions b(x) and p b (x) defined on satisfying b = 0, in, b n = g in, on Γ in, b t i = 0, on Γ in, b = 0, on \Γ in, where t i, i =,..., (d ) denotes an orthogonal set of tangent vectors on Γ in. p b = 0, in, p b = p out, on Γ out, p b n = 0, on \Γ out. (In case the pressure is specified on the inflow boundary Γ in, then b = 0, and the definition of p b is appropriately modified.) With the change of variables: u = u 0 + b and p = p 0 + p b, and subsequent relabeling u 0 = u, p 0 = p and f = p b we obtain the following system of modeling equations: where β( u s + b ) = ν eff ( u s + b ) k. β( u s + b )u + β( u s + b )b + p = f, in, (.6) u = 0, in, (.7) u n = 0, on Γ in Γ, (.8) p = 0, on Γ out, (.9) In the next section we show that, under suitable assumptions on β( ) and u s, there exists a unique solution to (.6)-(.9). An approximation scheme is presented in Section 3, and an a priori error estimate derived. A numerical example illustrating the approximation scheme and the a priori error estimate is presented in Section 4. 2 Existence and Uniqueness In this section we investigate the existence and uniqueness of solutions to the nonlinear system equations (.6)-(.9). We assume that IR d, d = 2 or 3, is a convex polyhedral domain and for vectors in IR d denotes the Euclidean norm. Throughout, we use C to denote a generic nonnegative constant, independent of the mesh parameter h, whose actual value may change from line to line in the analysis. 3
4 We make the following assumptions on β( ) and u s. Assumptions on β( ) Aβ : β( ) : IR + IR +, Aβ2 : 0 < β min β(s) β max, s IR +, Aβ3: β is Lipschitz continuous, β(s ) β(s 2 ) C β s s 2. Assumptions on u s Au s : For u L 2 (), u s L () C s u L 2 (), Au s 2: For {u n } n= L2 (), with u n converging weakly to u L 2 (), then {u s n} n= u s in L (), Au s 3: The mapping u u s is linear. Weak formulation of (.6)-(.9) Let X = {v H div () : v n = 0, on Γ in Γ}. We use (f, g) := f g d, and f := (f, f) /2 converges to to denote the L 2 inner product and the L 2 norm over, respectively, for both scalar and vector valued functions. Additionally, we introduce the norm v X = ( ) /2 ( v v + v v) d. Remark: For v H div () it follows that v n H /2 ( ). For the interpretation of the condition v n = 0 on Γ in Γ see [9, 5]. We restate (.6)-(.9) as: Given b, f L 2 (), find (u, p) X L 2 (), such that for all v X and q L 2 () (β( u s + b )u, v) (p, v) + (β( u s + b )b, v) = (f, v), (2.) (q, u) = 0. (2.2) For the spaces X and L 2 () we have the following inf-sup condition inf sup (q, v) c 0 > 0. (2.3) q L 2 () v X q v X We begin by establishing boundedness of any solution to (2.)-(2.2). Lemma 2. Any solution (u, p) X L 2 () to (2.)-(2.2) satisfies u X + p C ( b + f ). (2.4) Proof: From (2.2) and that X L 2 () we have that any solution u to (2.)-(2.2) satisfies u = 0. (2.5) 4
5 With the choice v = u, q = p, subtracting (2.2) from (2.), and using assumption Aβ2 yields (β( u s + b )u, u) = (β( u s + b )b, u) + (f, u), β min u 2 β max b u + f u. (2.6) Combining (2.5) and (2.6) we obtain the stated bound for u. The estimate for p is obtained using the inf-sup condition (2.3). p (p, v) sup c 0 v X = sup v X c 0 v X ( β( u s + b )u + β( u s + b )b + f ) c 0 c 0 (β max ( u + b ) + f ), from which the stated bound follows. (β( u s + b )u, v) + (β( u s + b )b, v) (f, v) v X Define Z = { v X : (q, v) = 0, q L 2 () }. Because of the inf-sup condition (2.3), the weak formulation (2.)-(2.2) can be equivalently stated as: Given b, f L 2 ()), find u Z, such that for all v Z (β( u s + b )u, v) + (β( u s + b )b, v) = (f, v). (2.7) Remark: For v Z, v X = v, as v = 0. To establish the existence of a solution to (2.7) we use the Leray-Schauder fixed point theorem. To do this we show that a solution to (2.7) is a fixed point of a compact mapping Φ. Theorem 2. For β( ) and u s satisfying assumptions Aβ Aβ3 and Au s Au s 2, respectively, there exists a solution u to (2.7). Proof: Let Φ : Z Z be defined by Φ(u) = w, where w satisfies (β( u s + b )w, v) + (β( u s + b )b, v) = (f, v). (2.8) That Φ is well defined follows from Aβ2 and the Lax-Milgram theorem. To show that Φ is a compact operator, let {u n } n= denote a bounded sequence in Z. From {u n} n= we can extract a subsequence, which we again denote as {u n } n=, such that {u n} n= converges weakly to u Z. For w n = Φ(u n ), using (2.8) (β( u s + b )w, v) (β( u s n + b )w n, v) = (β( u s + b )b, v) + (β( u s n + b )b, v) (β( u s n + b ) (w w n ), v) = ((β( u s + b ) β( u s n + b )) w, v) ((β( u s + b ) β( u s n + b )) b, v). 5
6 With v = w w n, and using Aβ2 and Aβ3 β min w w n 2 C β ( u s + b u s n + b ) w w w n + C β ( u s + b u s n + b ) b w w n C β u s u s n w w w n + C β u s u s n b w w n C β d u s u s n L () w w w n + C β d u s u s n L () b w w n w w n X = w w n C β d u s u s β n L () ( w + b ), min from which, with Au s 2, we can conclude that Φ is a compact operator. For r = βmax β min ( b + f ), from Lemma 2. we have that Φ(u) r, u Z. Then, applying the Leray-Schauder fixed point theorem [7] we obtain that there exists a u Z such that u = Φ(u). Under small data conditions we have the following theorem guaranteeing uniqueness of solutions to (2.7). Theorem 2.2 { With the stated assumptions Aβ Aβ3 and Au s Au s 2, and the condition that } b max β min /β max, β min /(C β d Cs ), if a solution u to (2.7) exists satisfying { } β min β u < max, min b, (2.9) β max C β d Cs then there is no other solution to (2.7). Proof: Suppose that both u and w Z satisfy (2.7), i.e., together with (2.7) we have that (β( w s + b )w, v) + (β( w s + b )b, v) (f, v) = 0, v Z. (2.0) With v = u w, subtracting (2.0) from (2.7) and using the bounds for β( ) we obtain (β( w s + b )(u w), (u w)) + ((β( u s + b ) β( w s + b ))u, (u w)) + ((β( u s + b ) β( w s + b ))b, (u w)) = 0 (2.) β min u w 2 β max u u w + β max b u w (β min β max ( u + b )) u w 0 (2.2) w = u, provided u < β min β max b. (2.3) Alternatively, from (2.), using Aβ3 and Au s, β min u w 2 C β d u s w s L () u u w + C β d u s w s L () b u w C β d Cs ( u + b ) u w 2 (β min C β d Cs ( u + b )) u w 2 0 β min w = u, provided u < b. C β d Cs 6
7 3 Finite Element Approximation In this section we investigate the finite element approximation to (u, p) satisfying (.6)-(.9). Let T h be a triangulation of made of triangles (in IR 2 ) or tetrahedrons (in IR 3 ). computational domain is defined by = K Th K. We assume that there exist constants c, c 2 such that Thus, the c h h K c 2 ρ K, where h K is the diameter of triangle (tetrahedron) K, ρ K is the diameter of the greatest ball (sphere) included in K, and h = max K Th h K. For k IN, let P k (A) denote the space of polynomials on A of degree no greater than k, and RT k (T h ) the (Piola) affine transformation of the Raviart-Thomas elements of order k on the unit triangle. We define the finite element spaces X h, X s h and Q h as follows. Note that as X h X h := {RT k (T h ) X}, (3.) X s h := { v X C 0 () : v K P l (K), K T h }, (3.2) Q h := { q L 2 () : q K P k (K), K T h }. (3.3) Additionally, let Z h := {v X h : (q, v) = 0, q Q h }. (3.4) Q h, for v Z h we have that v = 0, thus v X = v. For X h and Q h defined in (3.) and (3.3), the following discrete inf-sup condition is satisfied inf q Q h (q, v) sup c 0 > 0. (3.5) v X h q Q v X With X h, Z h, Q h defined above, we have the following approximation properties [4, 3]. For u Z H k+ () and p H k+ () inf v Z h u v X = inf v Z h u v C inf v X h u v = C h k+ u H k+ (), (3.6) inf q Q h p q C h k+ p H k+ (). (3.7) The approximation scheme we investigate is: Given b, f L 2 (), determine (u h, p h ) X h Q h, satisfying (β( u s h + b )u h, v) (p h, v) + (β( u s h + b )b, v) = (f, v), v X h (3.8) (q, u h ) = 0, q Q h. (3.9) Regarding u s h, note that applying a smoother to a function v X h (typically) does not result in v s Xh s. Therefore, we let ũs h Hl+ () C 0 () denote the result of the smoother applied to u h, and define u s h (x) = I hũ s h (x), (3.0) 7
8 where I h : C 0 () Xh s denotes an interpolation operator. We assume that the smoothed velocity ũ s h dependent on ũ s h, C ũ s such that h is sufficiently regular such that there exists a constant The precise dependence of Cũs h ũ s h I hũ s h L () Cũs h h l+. (3.) on ũ s h will depend on the particular smoother used. The existence, uniqueness, and boundedness of the solutions (u n h, pn h ) to (3.8)-(3.9) are established in a completely analogous manner as for the continuous problem. Corollary 3. (See Lemma 2..) Any solution (u, p) X h Q h to (3.8)-(3.9) satisfies u h X + p h C ( b + f ). (3.2) Corollary 3.2 (See Theorem 2..) For β( ) and u s h satisfying assumptions Aβ Aβ3 and Aus Au s 2, respectively, there exists a solution (u h, p h ) to (3.8)-(3.9). Proof: The existence of u h is established as that for u in Theorem 2.. The existence of p h then follows from the discrete inf-sup condition (3.5). In the next lemma we present the a priori error estimate for the approximation given by (3.8)-(3.9). Lemma 3. For (u, p) H k+ () X H k+ () satisfying (2.)-(2.2), (u h, p h ) satisfying (3.8)- (3.9), and u satisfying the small data condition C β d Cs ( u + b ) < β min, (3.3) and assuming that Cũs h given in (3.) is bounded by a constant C u, we have that there exists C > 0 such that ( u u h X + p p h C h k+ u H k+ () + h k+ p H k+ () + C u h l+). (3.4) Remark: The condition (3.3) guarantees uniqueness of the solution to (3.8)-(3.9), see Theorem 2.2. Proof: We have that the solutions u h and u to (3.8)-(3.9) and (2.)-(2.2), respectively, satisfy the following equations for all v Z h : and (β( u s h + b )u h, v) + (β( u s h + b )b, v) = (f, v), (3.5) (β( u s h + b )u, v) + (β( us h + b )b, v) = (f, v) ((β( us + b ) β( u s h + b )) u, v) ((β( u s + b ) β( u s h + b )) b, v). (3.6) 8
9 With e = u u h, subtracting equations (3.5) and (3.6) we obtain (β( u s h + b )e, v) = ((β( us + b ) β( u s h + b )) u, v) ((β( u s + b ) β( u s h + b )) b, v), v Z h. (3.7) For U Z h, let e = (u U) + (U u h ) := Λ + E. Then, for v = E, (3.7) becomes (β( u s h + b )E, E) = (β( us h + b )Λ, E) ((β( us + b ) β( u s h + b )) u, E) ((β( u s + b ) β( u s h + b )) b, E). (3.8) Next we bound each of the terms in (3.8). (β( u s h + b )E, E) β min E 2. (3.9) (β( u s h + b )Λ, E) β max Λ E ɛ E 2 + 4ɛ β2 max Λ 2. (3.20) ((β( u s + b ) β( u s h + b )) u, E) (β( us + b ) β( u s h + b )) u E C β d u s u s h L () u E C β d ( u s ũ s h L () + ũ s h us h L ()) u E C β d ( Cs u u h + ũ s h I hũ s h L ()) u E C β d ( Cs ( Λ + E ) + ũ s h I hũ s h L ()) u E C β d Cs u E 2 + ɛ 2 E 2 + 2ɛ 2 C 2 β d u 2 ( C 2 s Λ 2 ) + ũ s h I hũ s h 2 L () (3.2). A similar bound to that given in (3.2) holds for the third term on the right hand side of (3.8). Combining the estimates (3.9)-(3.2) with (3.8) we have ) (β min ɛ C β d Cs ( u + b ) 2ɛ 2 ( 4 ɛ β 2 max + 2 ɛ 2 C 2 β d C2 s ( u 2 + b 2 ) E 2 ) Λ ɛ 2 C 2 β d ( u 2 + b 2 ) ) ũ s h I hũ s h 2 L (). (3.22) Hence, in view of the stated hypothesis (3.3), there exists C > 0 such that E C ( Λ + ũ s h I hũ s h L ()). Finally, from the triangle inequality and (3.6) we have that ( u u h X = u u h Λ + E C h k+ u H k+ () + C u h l+). (3.23) To obtain the error estimate for the pressure, let P Q h. Then, from (3.5) we have that there exists v X h such that c 0 P p h (P p h, v) v X = (P, v) (p h, v) v X. 9
10 Using (3.8) and (2.) we obtain c 0 v X P p h (P, v) + (f, v) (β( u s h + b )u h, v) (β( u s h + b )b, v) = (P, v) (p, v) + (β( u s + b )u, v) + (β( u s + b )b, v) (β( u s h + b )u h, v) (β( u s h + b )b, v) = (P p, v) + (β( u s h + b ) (u u h), v) + ((β( u s + b ) β( u s h + b )) u, v) + ((β( u s + b ) β( u s h + b )) b, v) c 0 P p h P p + β max u u h ( + C β d Cs u u h + ũ s h I hũ h L ()) ( u + b ). Using the triangle inequality, (2.4), (3.7), (3.) and (3.23) we obtain the stated estimate for p p h. Remark: The L () norm used for the term (ũ s h I hũ s h ) and the L2 () norm used for u and b in (3.22) may be interchanged, assuming that the functions u and b are sufficiently regular. 4 Numerical Computations In this section we present a numerical example to demonstrate the numerical approximation scheme (3.8)-(3.9), and investigate the a priori error estimate (3.4). Let = (, ) (0, ), β(s) = v + (v 0 v )/( + ks 2 r ), with parameters v =, v 0 = 5, k =, and r = /2. (β( ) represents the Cross model for the effective viscosity for a generalized Newtonian fluid.) The true solution u and p are taken to be [ ] sin(πx) cos(πy) u(x, y) =, p(x, y) = xy. (4.) cos(πx) sin(πy) For this choice of u, u 0, hence a right hand side function is added to (3.9). The boundary conditions used are u n along {} (0, ), (, ) {}, { } (0, ), with p = 0 weakly imposed along (, ) {0}. A computation mesh corresponding to mesh parameter h = /4 is presented in Figure 4.. Plots of β( u ), u and p are given in Figures 4.2, 4.3 and 4.4, respectively. Example. For u s h, the interpolate of ũs h (the smoothed function of u h), we compute a continuous, piecewise quadratic, velocity by taking a simple average of u h at the nodal points of u s h. Computations were performed using RT 0 discp 0, RT discp, and RT 2 discp 2 elements for the velocity and pressure. (By RT k we are referring to Raviart-Thomas elements of degree k, and discp k refers to the space of discontinuous scalar functions which are polynomials of degree less that or equal to k on each triangle in the triangulation.) The results, together with the experimental convergence rates are presented in Table 4.. The experimental convergence rates are consistent with those predicted by (3.4) for l = 2. (Regarding the O(h 4 ) experimental convergence rate for the pressure using RT 2 discp 2 elements, note that the true solution for the pressure lies in the discp 2 approximation space.) Example 2. In order to investigate the dependence of the approximation on the interpolant of the smoother, 0
11 Figure 4.: Computational mesh for h = /4. Figure 4.2: Plot of β( u ). Velocity Figure 4.3: Plot of the velocity flow field u. Figure 4.4: Plot of the pressure function p. in this case we take u s h to be a continuous, piecewise linear function, obtained by taking a simple average of ũ s h at the vertices of the triangles in the triangulations. The results obtained using RT discp, and RT 2 discp 2 approximating elements are presented in Table 4.2. In this case (l = ) we observe optimal convergence for RT discp (and RT 0 discp 0, results not included). However, the experimental convergence rates for the RT 2 discp 2 approximation is limited to 2 for the velocity and pressure, consistent with (3.4). A Example of a local smoothing function In this section we give an example of a local smoothing function which satisfies properties Au s and Au s 2 presented in Section 2. The smoothing function is a simple averaging operator. We use the term domain to refer to an open connected set in IR n. For simplicity we present the case for a scalar function u(x). smoother is simply applied to each of the coordinate functions. For a vector valued function the
12 h u u h L 2 () Cvg. rate (u u h ) L 2 () Cvg. rate p p h L 2 () Cvg. rate X h = RT 0 Q h = discp 0 / E E E-2.29 / E E E-2.0 /8.790E E E-2.08 /0.433E E E-2.05 /2.95E E E-2 Pred X h = RT Q h = discp / E E E / E E E /8.456E E E / E E E / E E E-04 Pred X h = RT 2 Q h = discp 2 /4 6.66E E E /6.905E E E / E E E / E E E /2 2.32E E E-05 Pred Table 4.: Example, u s h a quadratic interpolant of ũs h. h u u h L 2 () Cvg. rate (u u h ) L 2 () Cvg. rate p p h L 2 () Cvg. rate X h = RT Q h = discp / E E E-2.99 /6 3.50E E E /8.808E E E /0.70E E E /2 8.69E E E-3 Pred X h = RT 2 Q h = discp 2 / E E E-2.7 /6.727E E E / E E E /0 6.30E E E / E E E-3 Pred Table 4.2: Example 2, u s h a linear interpolant of ũs h. Let denote a bounded domain in IR n and L() the Lebesgue measurable sets in. Let δ > 0 denote the (fixed) volume measure over which we average a function to obtain its smoothed value. 2
13 For x the typical averaging volume which comes to mind is B(x, r δ ), where B(x, r δ ) denotes the ball centered at x of radius r δ having volume δ. As δ is fixed the difficulty in using B(x, r δ ) arises for points whose distance from is less that r δ. This requires us to consider averaging volumes other than balls. Namely, for each point x we associate a domain V (x) having a volume of δ. We require that the association of x with V (x) be continuous. This continuity is formally described in the next paragraph. Let ν denote the Lebesgue measure in IR n. For S, S 2 L(), introduce the metric d(s, S 2 ) defined by d(s, S 2 ) := ν(s S 2 ), where S S 2 := (S \S 2 ) (S 2 \S ). (A.) Now, let V : L() satisfy: (i) V (x) is a domain with ν(v (x)) = δ for all x, and (ii) d(v (x), V (y)) = ν(v (x) V (y)) C V x y for all x, y, where C V a fixed constant. For convenience we denote the domain V (x) as V x. Definition: Local Smoothing Operator For u L 2 (), define u s as u s (x) = δ We have the following properties for u s (x). V x u(z) d. (A.2) Lemma A. For u L 2 (), u s defined by (A.2) satisfies the following properties. (i) u s L () δ /2 u L 2 (). (ii) u s : IR is uniformly continuous. (iii) Suppose that {u n } n= L2 () and that u n converges weakly to u L 2 (). Then {u s n} n= converges to u s in L (). Proof: Let S L 2 () denote the characteristic function of the domain S. From (A.2), for x u s (x) = u(z) d = Vx u(z) d δ V x δ ( ) /2 ( ) /2 ( Vx ) 2 d u(z) 2 d δ which establishes (i). = δ /2 u L 2 (), For x, y, u s (x) u s (y) δ = δ ( Vx Vy u(z) d ( Vx Vy ) 2 d ) /2 ( = δ u L 2 () d(v (x), V (y)) /2 = C/2 V δ u L 2 () x y /2, ) /2 u(z) 2 d 3
14 which establishes the uniform continuity of u s. As u s is bounded on then u s can be continuously extended to. To establish (iii), as {u n } converges weakly, let sup n u n = M <. In addition, for ɛ > 0, σ = ( ) 2, ɛ/(6 M C /2 V ) let {zi } N i= denote a σ-net of, i.e., for all x there exists an i x {, 2,..., N} such that x z ix < σ. Now, u s n(x) u s (x) = = (u n (y) u(y)) d V x (u n (y) u(y)) d + (u n (y) u(y)) d V zix V x\v zix (u n (y) u(y)) d V zix + u n (y) u(y) d. V x V zix (A.3) Since {u n } converges weakly to u in L 2 (), for all w L 2 () there exists N w such that for n > N w (u n u) w d < ɛ 3. (A.4) Let N = max i=,2,...,n { N Vzi }. Then, for n > N (u n (y) u(y)) d V zix = For the second term on the right hand side of (A.3) we have (u n (y) u(y)) Vzi d < ɛ 3. V x V zix u n (y) u(y) d ( V x V zix u n (y) u(y) 2 d 2M ν(v x V zix ) /2 ) /2 ( 2M C /2 V x z ix /2 2M C /2 V σ /2 V x V zix d ) /2 = ɛ 3. (A.5) Thus, from (A.3)-(A.5) it follows that for all x, for n > N u s n(x) u s (x) < 2 3 ɛ, i.e., us n u s L () < 2 3 ɛ < ɛ. A. Regularity of u s (for u L ()) If, in place of u L 2 (), we have u L () then u s defined by (A.2) is a H () function. To establish this regularity result we begin by citing a characterization of the W,p (IR n ) function space. 4
15 Theorem A. ( [8], Theorem 2..6) Let < p <. Then u W,p (IR n ) if and only if u L p (IR n ) and ( u(x + h) u(x) p /p IR n h dx) = h u(x + h) u(x) L p (IR n ) (A.6) remains bounded for all h IR n. Theorem A.2 If u L () then, for u s defined by (A.2), u s H (). Proof: In order to apply Theorem A. we need to define an extension of u to IR n. Let { u(x), x ũ(x) = 0, x, and Ṽ : IRn L(IR n ) denote an extension of V satisfying properties (i) and (ii) (with replaced by IR n ), and additionally that there exists constants C > 0 and C 2 0 such that (iii) diameter(ṽ (z)) C for all z IR n, and (iv) sup z IR n inf y Ṽ (z) z y C 2. Let B denote the bounded set, B := {x IR n : inf y x y < + C + C 2 } support(ũ s ). Note that for x IR n \ B and h <, ũ s (x + h) = 0. Now, for h, ũ s (x + h) ũ s (x) IR n h For h <, ũ s (x + h) ũ s (x) IR n h 2 2 dx dx = = = ( ) 2 h 2 IR (ũs (x + h)) 2 dx + n IR (ũs (x)) 2 dx n 4 h 2 IR (ũs (x)) 2 dx 4 n h 2 ũs 2 L ( B ) ν( B) 4 ν( B ) ũ 2 L (IR n = 4 ν( ) B ) u 2 L (). (A.7) h 2 ũ s (x + h) ũ s (x) 2 dx B ( h 2 Ṽx+h (z) Ṽx (z) δ B h 2 h 2 ũ(z) B δ 2 ũ 2 L ( B ) B ( B ) dz Ṽx+h (z) Ṽx (z) d(ṽx+h, Ṽx) 2 dx δ 2 ũ 2 L ( B ) B h 2 δ 2 ũ 2 L ( B ) C2 V h 2 ν( B ) 2 dx 2 dz) dx = δ 2 C2 V ν( B ) ũ 2 L ( B ) = δ 2 C2 V ν( B ) u 2 L (). (A.8) From (A.7) and (A.8), together with Theorem A., we obtain that ũ s H (IR n ). As u s = ũ s, it then follows that u s H (). 5
16 B Example of a differential smoothing function As an alternative to the local averaging filter discussed in Section A, in this section we present a differential smoothing filter. Let X s = H 0 () = {v H () : v = 0 on } X. (B.9) Definition: Differential Smoothing Operator For u L 2 (), define u s X s as ( u s, v) = (u s, v), v X s. (B.0) The well posedness of u s follows from an application of the Lax-Milgram theorem. Next we show that this smoothing operation satisfies properties Au s and Au s 2 presented in Section 2. Lemma B.2 For u L 2 (), u s defined by (B.0) satisfies the following properties. (i) u s L () C u L 2 (). (ii) Suppose that {u n } n= L2 (), and that u n converges weakly to u L 2 (). The {u s n} converges to u s in L (). Proof: From (B.0) we have that u s X s, and as u L 2 (), from the shift theorem (together with a sufficiently smooth ), it follows that u s H 2 () X s, with u s H 2 () C u. (B.) Using the embedding of H 2 () in L () we establish (i). Let W : L 2 () H 2 () X s, W(u) := u s, denote the filter mapping. Then from (B.) W is a bounded (linear) transformation from L 2 () H 2 () X s. Let W : ( H 2 () X s) L 2 () denote the adjoint operator of W. (The existence of W follows immediately from the Riesz Representation Theorem.) Now, for η ( H 2 () X s) u s n u s, η H 2,(H 2 ) = W(u n) W(u), η H 2,(H 2 ) = W(u n u), η H 2,(H 2 ) = (u n u, W (η)) 0 as n, as u n converges weakly in L 2 () to u. Hence as H 2 () X s is compactly embedded in L () X s, then u s n converges to u s strongly in L (). References [] G. Allaire. Homogenization of the Navier-Stokes equations with a slip boundary condition. Comm. Pure Appl. Math., 44(6):605 64, 99. 6
17 [2] M. Azaïez, F. Ben Belgacem, C. Bernardi, and N. Chorfi. Spectral discretization of Darcy s equations with pressure dependent porosity. Appl. Math. Comput., 27(5): , 200. [3] S.C. Brenner and L.R. Scott. The mathematical theory of finite element methods, volume 5 of Texts in Applied Mathematics. Springer-Verlag, New York, second edition, [4] F. Brezzi and M. Fortin. Mixed and hybrid finite element methods. Springer-Verlag, New York, 99. [5] M. Bulíček, P. Gwiazda, J. Málek, and A. Świerczewska-Gwiazda. On steady flows of incompressible fluids with implicit power-law-like rheology. Adv. Calc. Var., 2(2):09 36, [6] L. Diening, C. Kreuzer, and E. Süli. Finite element approximation of steady flows of incompressible fluids with implicit power-law-like rheology. SIAM J. Numer. Anal., 5(2):984 05, 203. [7] V.J. Ervin, E.W. Jenkins, and S. Sun. Coupled generalized nonlinear Stokes flow with flow through a porous medium. SIAM J. Numer. Anal., 47(2): , [8] V.J. Ervin, E.W. Jenkins, and S. Sun. Coupling nonlinear Stokes and Darcy flow using mortar finite elements. Appl. Numer. Math., 6():98 222, 20. [9] J. Galvis and M. Sarkis. Non matching mortar discretization analysis for the coupling Stokes Darcy equations. Electron. Trans. Numer. Anal., 26: , [0] G.N. Gatica, R. Oyarzúa, and F.-J. Sayas. A twofold saddle point approach for the coupling of fluid flow with nonlinear porous media flow. IMA J. Numer. Anal., 32(3): , 202. [] V. Girault, F. Murat, and A. Salgado. Finite element discretization of Darcy s equations with pressure dependent porosity. M2AN Math. Model. Numer. Anal., 44(6):55 9, 200. [2] A. Hirn, M. Lanzendörfer, and J. Stebel. Finite element approximation of flow of fluids with shear-rate- and pressure-dependent viscosity. IMA J. Numer. Anal., 32(4): , 202. [3] J. Málek, J. Nečas, and K.R. Rajagopal. Global analysis of the flows of fluids with pressuredependent viscosities. Arch. Ration. Mech. Anal., 65(3): , [4] K. R. Rajagopal. On a hierarchy of approximate models for flows of incompressible fluids through porous solids. Math. Models Methods Appl. Sci., 7(2):25 252, [5] L. Traverso, T.N. Phillips, and Y. Yang. Mixed finite element methods for groundwater flow in heterogeneous aquifers. Computers & Fluids, 88:60 80, 203. [6] S. Whitaker. Flow in porous media I. A theoretical derivation of Darcy s law. Transp. Porous Media, :3 25, 986. [7] E. Zeidler. Applied functional analysis, volume 08 of Applied Mathematical Sciences. Springer- Verlag, New York, 995. Applications to mathematical physics. [8] W.P. Ziemer. Weakly differentiable functions, volume 20 of Graduate Texts in Mathematics. Springer-Verlag, New York, 989. Sobolev spaces and functions of bounded variation. 7
FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationCoupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements
Coupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements V.J. Ervin E.W. Jenkins S. Sun Department of Mathematical Sciences Clemson University, Clemson, SC, 29634-0975, USA. Abstract We study
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 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 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 informationApproximation of Axisymmetric Darcy Flow
Approximation of Axisymmetric Darcy Flow V.J. Ervin Department of Mathematical Sciences Clemson University Clemson, SC, 29634-0975 USA. January 6, 2012 Abstract In this article we investigate the numerical
More informationMathematical analysis of the stationary Navier-Stokes equations
Mathematical analysis of the Department of Mathematics, Sogang University, Republic of Korea The 3rd GCOE International Symposium Weaving Science Web beyond Particle Matter Hierarchy February 17-19, 2011,
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 informationA Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements
W I S S E N T E C H N I K L E I D E N S C H A F T A Domain Decomposition Method for Quasilinear Elliptic PDEs Using Mortar Finite Elements Matthias Gsell and Olaf Steinbach Institute of Computational Mathematics
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 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 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 informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
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 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 informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
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 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 informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
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 informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 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 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 informationProblem of Second grade fluids in convex polyhedrons
Problem of Second grade fluids in convex polyhedrons J. M. Bernard* Abstract This article studies the solutions of a three-dimensional grade-two fluid model with a tangential boundary condition, in a polyhedron.
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
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 informationNumerical Simulations of Fluid Pressure in the Human Eye
Numerical Simulations of Fluid Pressure in the Human Eye T.R. Crowder V.J. Ervin Department of Mathematical Sciences Clemson University Clemson, SC, 2963-0975 USA. January 22, 2013 Abstract In this article
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 informationStenberg s sufficiency criteria for the LBB condition for Axisymmetric Stokes Flow
Stenberg s sufficiency criteria for the LBB condition for Axisymmetric Stokes Flow V.J. Ervin E.W. Jenkins Department of Mathematical Sciences Clemson University Clemson, SC, 29634-0975 USA. Abstract In
More informationNonlinear Darcy fluid flow with deposition
Nonlinear Darcy fluid flow wit deposition V.J. Ervin Hyesuk Lee J. Ruiz-Ramírez Marc 23, 2016 Abstract In tis article we consider a model of a filtration process. Te process is modeled using te nonlinear
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
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 informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationVariable Exponents Spaces and Their Applications to Fluid Dynamics
Variable Exponents Spaces and Their Applications to Fluid Dynamics Martin Rapp TU Darmstadt November 7, 213 Martin Rapp (TU Darmstadt) Variable Exponent Spaces November 7, 213 1 / 14 Overview 1 Variable
More informationNumerical Methods for the Navier-Stokes equations
Arnold Reusken Numerical Methods for the Navier-Stokes equations January 6, 212 Chair for Numerical Mathematics RWTH Aachen Contents 1 The Navier-Stokes equations.............................................
More informationThe Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control
The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control A.Younes 1 A. Jarray 2 1 Faculté des Sciences de Tunis, Tunisie. e-mail :younesanis@yahoo.fr 2 Faculté des Sciences
More 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 informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
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 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 informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More informationMEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EQUATIONS AND THEIR APPLICATIONS. Yalchin Efendiev.
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX MEYERS TYPE ESTIMATES FOR APPROXIMATE SOLUTIONS OF NONLINEAR ELLIPTIC EUATIONS AND THEIR APPLICATIONS
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More informationNumerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell
J.S.Howell Research Statement 1 Numerical Analysis and Computation of Fluid/Fluid and Fluid/Structure Interaction Problems Jason S. Howell Fig. 1. Plot of the magnitude of the velocity (upper left) and
More informationAnalysis of the Oseen-Viscoelastic Fluid Flow Problem
Analysis of the Oseen-Viscoelastic Fluid Flow Problem Vincent J. Ervin Hyesuk K. Lee Louis N. Ntasin Department of Mathematical Sciences Clemson University Clemson, South Carolina 9634-0975 December 30,
More informationOn Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations
On Smoothness of Suitable Weak Solutions to the Navier-Stokes Equations G. Seregin, V. Šverák Dedicated to Vsevolod Alexeevich Solonnikov Abstract We prove two sufficient conditions for local regularity
More informationUniform inf-sup condition for the Brinkman problem in highly heterogeneous media
Uniform inf-sup condition for the Brinkman problem in highly heterogeneous media Raytcho Lazarov & Aziz Takhirov Texas A&M May 3-4, 2016 R. Lazarov & A.T. (Texas A&M) Brinkman May 3-4, 2016 1 / 30 Outline
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationSELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY
Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS
More informationICES REPORT A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams
ICS RPORT 15-17 July 2015 A Generalized Mimetic Finite Difference Method and Two-Point Flux Schemes over Voronoi Diagrams by Omar Al-Hinai, Mary F. Wheeler, Ivan Yotov The Institute for Computational ngineering
More informationExistence of global weak solutions to implicitly constituted kinetic models of incompressible homogeneous dilute polymers
1 / 31 Existence of global weak solutions to implicitly constituted kinetic models of incompressible homogeneous dilute polymers Endre Süli Mathematical Institute, University of Oxford joint work with
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationStrong Stabilization of the System of Linear Elasticity by a Dirichlet Boundary Feedback
To appear in IMA J. Appl. Math. Strong Stabilization of the System of Linear Elasticity by a Dirichlet Boundary Feedback Wei-Jiu Liu and Miroslav Krstić Department of AMES University of California at San
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 informationUniform estimates for Stokes equations in domains with small holes and applications in homogenization problems
Uniform estimates for Stokes equations in domains with small holes and applications in homogenization problems Yong Lu Abstract We consider the Dirichlet problem for the Stokes equations in a domain with
More informationHARNACK S INEQUALITY FOR GENERAL SOLUTIONS WITH NONSTANDARD GROWTH
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 37, 2012, 571 577 HARNACK S INEQUALITY FOR GENERAL SOLUTIONS WITH NONSTANDARD GROWTH Olli Toivanen University of Eastern Finland, Department of
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
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 informationA note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations
A note on discontinuous Galerkin divergence-free solutions of the Navier-Stokes equations Bernardo Cockburn Guido anschat Dominik Schötzau June 1, 2007 Journal of Scientific Computing, Vol. 31, 2007, pp.
More informationHOMEOMORPHISMS OF BOUNDED VARIATION
HOMEOMORPHISMS OF BOUNDED VARIATION STANISLAV HENCL, PEKKA KOSKELA AND JANI ONNINEN Abstract. We show that the inverse of a planar homeomorphism of bounded variation is also of bounded variation. In higher
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 On the divergence constraint in mixed finite element methods for incompressible
More informationINSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES. Note on the fast decay property of steady Navier-Stokes flows in the whole space
INSTITUTE OF MATHEMATICS THE CZECH ACADEMY OF SCIENCES Note on the fast decay property of stea Navier-Stokes flows in the whole space Tomoyuki Nakatsuka Preprint No. 15-017 PRAHA 017 Note on the fast
More informationA Product Property of Sobolev Spaces with Application to Elliptic Estimates
Rend. Sem. Mat. Univ. Padova Manoscritto in corso di stampa pervenuto il 23 luglio 2012 accettato l 1 ottobre 2012 A Product Property of Sobolev Spaces with Application to Elliptic Estimates by Henry C.
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 informationON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationMORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS
MORTAR MULTISCALE FINITE ELEMENT METHODS FOR STOKES-DARCY FLOWS VIVETTE GIRAULT, DANAIL VASSILEV, AND IVAN YOTOV Abstract. We investigate mortar multiscale numerical methods for coupled Stokes and Darcy
More informationSUBELLIPTIC CORDES ESTIMATES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 000-9939XX0000-0 SUBELLIPTIC CORDES ESTIMATES Abstract. We prove Cordes type estimates for subelliptic linear partial
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 informationarxiv: v1 [math.na] 19 Dec 2017
Arnold-Winther Mixed Finite Elements for Stokes Eigenvalue Problems Joscha Gedicke Arbaz Khan arxiv:72.0686v [math.na] 9 Dec 207 Abstract This paper is devoted to study the Arnold-Winther mixed finite
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 informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationConvergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions
BIT Numer Math (2010) 50: 781 795 DOI 10.1007/s10543-010-0287-z Convergence analysis of finite element approximations of the Joule heating problem in three spatial dimensions Michael J. Holst Mats G. Larson
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationFIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION PROBLEM
SIAM J. NUMER. ANAL. c 1998 Society for Industrial Applied Mathematics Vol. 35, No. 1, pp. 320 335, February 1998 016 FIRST-ORDER SYSTEM LEAST SQUARES (FOSLS) FOR PLANAR LINEAR ELASTICITY: PURE TRACTION
More informationEverywhere differentiability of infinity harmonic functions
Everywhere differentiability of infinity harmonic functions Lawrence C. Evans and Charles K. Smart Department of Mathematics University of California, Berkeley Abstract We show that an infinity harmonic
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More informationA posteriori error estimates for non conforming approximation of eigenvalue problems
A posteriori error estimates for non conforming approximation of eigenvalue problems E. Dari a, R. G. Durán b and C. Padra c, a Centro Atómico Bariloche, Comisión Nacional de Energía Atómica and CONICE,
More informationA connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes
A connection between grad-div stabilized FE solutions and pointwise divergence-free FE solutions on general meshes Sarah E. Malick Sponsor: Leo G. Rebholz Abstract We prove, for Stokes, Oseen, and Boussinesq
More informationPhilippe. Functional Analysis with Applications. Linear and Nonlinear. G. Ciarlet. City University of Hong Kong. siajtl
Philippe G Ciarlet City University of Hong Kong Linear and Nonlinear Functional Analysis with Applications with 401 Problems and 52 Figures siajtl Society for Industrial and Applied Mathematics Philadelphia
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 informationACM/CMS 107 Linear Analysis & Applications Fall 2017 Assignment 2: PDEs and Finite Element Methods Due: 7th November 2017
ACM/CMS 17 Linear Analysis & Applications Fall 217 Assignment 2: PDEs and Finite Element Methods Due: 7th November 217 For this assignment the following MATLAB code will be required: Introduction http://wwwmdunloporg/cms17/assignment2zip
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationA Stabilized Finite Element Method for the Darcy Problem on Surfaces
arxiv:1511.03747v1 [math.na] 12 Nov 2015 A Stabilized Finite Element Method for the Darcy Problem on Surfaces Peter Hansbo Mats G. Larson October 8, 2018 Abstract We consider a stabilized finite element
More informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationNumerical Analysis of a Higher Order Time Relaxation Model of Fluids
Numerical Analysis of a Higher Order Time Relaxation Model of Fluids Vincent J. Ervin William J. Layton Monika Neda Abstract. We study the numerical errors in finite element discretizations of a time relaxation
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
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 informationNonlinear elliptic systems with exponential nonlinearities
22-Fez conference on Partial Differential Equations, Electronic Journal of Differential Equations, Conference 9, 22, pp 139 147. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationEnergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations of the Navier-Stokes equations
INTRNATIONAL JOURNAL FOR NUMRICAL MTHODS IN FLUIDS Int. J. Numer. Meth. Fluids 19007; 1:1 [Version: 00/09/18 v1.01] nergy norm a-posteriori error estimation for divergence-free discontinuous Galerkin approximations
More informationSobolev Spaces. Chapter Hölder spaces
Chapter 2 Sobolev Spaces Sobolev spaces turn out often to be the proper setting in which to apply ideas of functional analysis to get information concerning partial differential equations. Here, we collect
More information