Coupling Non-Linear Stokes and Darcy Flow using Mortar Finite Elements
|
|
- Paul Page
- 5 years ago
- Views:
Transcription
1 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, , USA. Abstract We study a system composed of a non-linear Stokes flow in one subdomain coupled with a nonlinear porous medium flow in another subdomain. Special attention is paid to the mathematical consequence of the shear-dependent fluid viscosity for the Stokes flow and the velocity-dependent effective viscosity for the Darcy flow. Motivated by the physical setting, we consider the case where only flow rates are specified on the low and outflow boundaries in both subdomains. We recast the coupled Stokes-Darcy system as a reduced matching problem on the interface using a mortar space approach. We prove a number of properties of the nonlinear interface operator associated with the reduced problem, which directly yield the existence, uniqueness and regularity of a variational solution to the system. We further propose and analyze a numerical algorithm based on mortar finite elements for the interface problem and conforming finite elements for the subdomain problems. Optimal a priori error estimates are established for the interface and subdomain problems, and a number of compatibility conditions for the finite element spaces used are discussed. Numerical simulations are presented to illustrate the algorithm and to compare two treatments of the defective boundary conditions. Key words. Generalized non-linear Stokes flow; coupled Stokes and Darcy flow; defective boundary condition AMS Mathematics subject classifications. 65N30 Introduction The problem of approximating coupled Stokes Darcy flow has received considerable attention in the mathematics community over the past ten years. Many of the applications considered use a Newtonian fluid in both the Stokes and Darcy regions, where the motivating problem is often a coupled surface water / groundwater model see, for instance, [6, 4, 8, 22, 24, 27, 28]. Our interest in coupled flows arises from filtration applications, which continues to be an active area of research see, for instance, [3, 2, 20, 26, 23]. The purpose of the filtration mechanism can be the removal of particulates [23, 3, 26, 29] or impurities [20], but in all cases the availability of an accurate and efficient simulation tool would aid in the design and assessment of effective filtration devices. vjervin@clemson.edu, lea@clemson.edu, shuyu@clemson.edu.
2 Our particular focus is on the effective removal of debris particles from a molten polymer [29, 0, 5], which is a non-newtonian fluid. Earlier work on this filtration problem focused only on the flow in the porous, or Darcy, region [, 2, 25, 29, 30], but analysis of the fully coupled problem is essential to accurately simulate the transport of the suspended particles into the filter. The fully coupled problem for non-newtonian Stokes and Darcy flows was initially analyzed in []. Motivated by numerical implementation considerations, we have recast this coupled problem as a matching problem on the interface between the domains. This approach naturally gives rise to a parallel algorithm for the subproblems in the Stokes domain Ω f and the Darcy domain Ω p, and it also allows one to combine existing codes for Stokes and Darcy simulations to solve the coupled problem. Key ideas from [7, 22, 4, 7] are used in the new formulation. In the literature, two basic approaches have been used to analyze the coupled problem considered herein: Stokes Darcy coupling and Stokes Laplace coupling. In Stokes Darcy coupling, the velocity and pressure are resolved in both the fluid flow and the porous media domains. For the Stokes Laplace formulation the velocity and pressure are the unknowns in the fluid flow domain and the pressure is the only unknown in the porous media domain. In [7] Discacciati and Quateroni studied the coupled Stokes Laplace formulation. After showing existence and uniqueness for the coupled problem they reformulated the problem, using the Steklov-Poincaré operator, as an interface problem for the interfacial pressure. Parallel and serial implementations for the numerical approximation of the interface problem were reported in [8, 9]. There has been considerable work done on efficient numerical solution algorithms for the Stokes Laplace formulation. We refer the reader to [4, 9, 24] and the references therein. Using a Lagrange multiplier technique Layton, Schieweck, and Yotov [22] see [28] also introduced and studied the Stokes Darcy coupling approach. Galvis and Sarkis in [4] extended this approach and showed how the Lagrange multiplier space could be more appropriately defined. In [3] they investigated efficient preconditioning strategies for the Stokes Darcy formulation. Our interest herein is the Stokes Darcy coupling approach. Denoting by u f, p f, u p, p p the velocity and pressure in the fluid flow domain Ω f and the porous media domain Ω p, respectively, we have the following boundary conditions holding along the interface Γ: u f n f + u p n p = 0, Conservation of mass. p f σ f n f n f = p p, Balance of interfacial pressure.2 and a boundary condition for the tangential stress in Ω f on Γ. The boundary condition for the tangential stress becomes a natural boundary condition for the variational formulation on Ω f. In.2 σ f denotes the extra stress tensor in Ω f. For the coupled problem under consideration, fully described in Section 2 below, given an interfacial pressure p I, the Stokes and Darcy problems are independently solvable in Ω f and Ω p, respectively, to yield u f p I, p f p I, u pp I, p pp I. Thus, we investigate the problem of determining λ such that Λλ := u f λ n f + u pλ n p = 0..3 The matching problem lies on the interface Γ. For the discrete approximation we introduce a mortar space for the representation of the interfacial pressure of Γ [7, 6]. The mortar space discretization is independent of the partitions of Ω f and Ω p. However a compatibility condition for the mortar space discretization and the partitions of Ω f and Ω p is required see Assumptions 4.4,
3 For the modeling equations describing the fluid flow we consider the case of a shear thinning fluid, where the fluid s viscosity is a nonlinear function of the magnitude of the deformation tensor. Additionally we assume that only flow rates are specified for the low and outflow boundaries i.e. defective boundary conditions. The paper is organized as follows. In Section 2 we formally describe the modeling equations and assumptions on the non-linear functions modeling the fluid s viscosity. Variational formulations for the Stokes problem in Ω f, and Darcy problem in Ω p, are given in Section 3 where the interfacial pressure, p I, is treated as a known quantity for both problems. Existence and uniqueness of the problems is verified and solutions are shown to depend continuously on p I. The matching problem of determining p I such that. is satisfied is then formulated and shown to have a unique solution. In Section 4 we present the finite element approximation scheme and analysis. The analysis for the discrete problem follows the same pattern as the continuous problem. First the discrete Stokes and Darcy problems are separately considered. Existence and uniqueness of the discrete approximations is verified and error estimates between the continuous solutions and discrete approximations are derived. The discrete matching problem is then formulated and shown to have a unique solution. A combined error estimate is then given for the solution between the coupled Stokes Darcy problem and its discrete approximation. In Section 5 a numerical example is given for two different treatments of the defective boundary conditions. 2 Modeling Equations For convenience, we work with Hilbert spaces instead of general Sobolev spaces for our weak formulations. This requires us to assume that the viscosity is bounded from above. A general formulation involving the spaces W,r Ω f and W /r,r Γ can be found in []. We would like to mention that there exists a certain physical limit for the viscosity, so our assumption is actually physically reasonable. The power law, for example, is an approximation of the reality, and the implication from the power law that the viscosity goes to inity when the deformation goes to zero is not physically correct as the power law does not apply to the case of a zero strain tensor. For convenience of implementation, we use Einstein s notation in addition to the vector notation. All repeated subscript indices imply summation over all dimensions except the subscript m for the tangential vector t mi, where we list the summation of m explicitly. Let Ω R n, n = 2 or 3, denote the flow domain of interest. Let Ω f Ω and Ω p Ω be bounded Lipschitz domains for the non-linear Stokes flow and non-linear Darcy flow, respectively. In this work, we consider only two subdomains, and assume Ω = Ω f Ω p Γ I, where Γ I := Ω f Ω p, denoting a smooth interface between the subdomains. 3
4 We first consider the non-linear Stokes flow in Ω f. p σ ij = f i, x i x j in Ω f, 2.4 u i = 0, x i in Ω f, 2.5 σ ij = 2ν f Du Du ij, in Ω f, 2.6 u i = 0, on Γ f,d, 2.7 σ ij n j pn i = τ B,i, on Γ f,n, 2.8 α S u i + σ ij n j t mi = 0, on Γ I, m =,, d, 2.9 u i n i ds = q f,k ds = q f,k measγ f,k, k =, 2,, K f. 2.0 Γ f,k Γ f,k Note that u B and τ B denote vector functions on the domain boundary, where Ω f = Γ f,d Γ f,n Γ f,f Γ I this is understood as Ω f = Γ f,d Γ f,n Γ f,f Γ I and similar notational simplification applied to the rest of this paper, Γ f,f := K f k= Γ f,k and K f is the number of defective boundaries in the Stokes flow domain. The boundaries Γ f,d, Γ f,n, Γ f,f, and Γ I are pairwise disjoint, and their definitions are clear from the following modeling equations for Stokes flow. We use D to denote the symmetric gradient operator. In other words, Du is the deformation tensor Du ij = ui 2 x j + u j x i. The viscosity ν f depends nonlinearly on Du. We assume that ν f is bounded from above and from below; that is, there exist ν min > 0 and ν max such that ν min ν f Du ν max for all u. We assume that measγ f,n > 0 and measγ I > 0. On Γ I, we impose the Beavers-Joseph-Saffman BJS condition 2.9. On each Γ f,k, we impose a defective boundary condition, where the averaged flow rate per area is imposed; or equivalently, the integral of the flow rate q f,k is specified on Γ f,k. Note that we assume that α S is a scalar for simplicity of presentation, but all results in this paper can be easily extended to treat a full tensor α S,ij. The Darcy flow in Ω p is modeled using the equations: ν p u R ij u j + p = f i, x i in Ω p, 2. u i = 0, x i in Ω p, 2.2 u i n i = 0, on Γ p,n, 2.3 p = p B, on Γ p,d, 2.4 u i n i ds = q p,k ds = q p,k measγ p,k, k =, 2,, K p, 2.5 Γ p,k Γ p,k where Ω p = Γ p,n Γ p,d Γ p,f Γ I. We assume that the permeability tensor K associated with the porous medium is symmetric positive definite, and we denote its inverse by the flow resistivity tensor R; that is, R = K. We also assume that ν p u is bounded from above by ν max and from below by ν min > 0. Similar to the Stokes region, we have imposed a defective boundary condition on the Darcy region boundary Γ p,k, with q p,k R being the imposed averaged flow rate; q p R Kp may be also viewed as a piecewise constant function defined on Γ p,f = K p k= Γ p,k. 4
5 To couple the two flow systems, we impose conservation of mass and balance of the normal forces across the interface Γ I : [u i n i ] ΓI = 0, [p σ ij n i n j ] ΓI = 0, where we extend σ ij by zero to Ω p for notational ease; that is σ ij Ωp = 0. We will use p f := p Ωf, u f := u Ωf, f f := f Ωf, p p := p Ωp, u p := u Ωp, f p := f Ωp to emphasize the corresponding variables applied to a specific region, although we will simply write p, u and f when the region we refer to is clear from the context. On the interface Γ I, n f = n Ωf = n p = n Ωp denotes the unit normal vector pointing from Ω f toward Ω p. When we restrict our attention to the Stokes or Darcy region, we simply write n f or n p as n. Throughout the paper, we use C to denote a generic positive constant, and ɛ to denote a fixed positive constant that may be chosen arbitrarily small. We assume that the nonlinear functions ν f Du Du and ν p u u are uniformly continuous with regard to Du and u, respectively; that is, there exists a constant C > 0 such that ν f T T ν f S S C T S, T R d d, S R d d, ν p u u ν p v v C u v, u R d, v R d. In addition, we assume that ν f and ν p are strictly monotone; namely, there exists a constant C > 0 such that: ν f T T ν f S S : T S CT S : T S, T R d d, S R d d, ν p u u ν p v v u v Cu v u v, u R d, v R d. 3 Variational Formulations with Mortar Spaces We first consider the Stokes flow and the Darcy flow problems separately as if they were two independent processes. This corresponds to the modularity of their code implementation, where it is desirable to implement and test individual codes for separate Stokes flow and Darcy flow before numerically coupling them. The variational formulations for the fully coupled problem is presented and analyzed in []. 3. Stokes flow We now restrict our attention to the Stokes flow problem, by assuming the interface pressure p I is given. The equations to be solved include together with the balance of the normal forces across the interface Γ I : σ ij n i n j p = p I, on Γ I. { Let X f = v : v H Ω f } d, v Γf,D = 0. 5
6 To derive the weak formulation, we multiply 2.4 by a vector function v X f, integrate over Ω f and apply the divergence theorem Green s theorem to obtain p f i, v i Ωf = σ ij v i dx x i x j Noting that = = Ω f Ω f Ω f σ ij Dv ij p v i x i dx 2ν f Du Du ij Dv ij p v i x i σ ij n j pn i v i ds Ω f dx τ B,i v i ds Γ f,n σ ij n j pn i v i ds σ ij n j pn i v i ds, Γ I Γ f,f d m= t mi t mj + n i n j = δ ij, we can treat the BJS slip boundary condition as σ ij n j pn i v i ds = Γ I σ ij n j pn i δ ik v k ds Γ I = = = σ ij n j pn i t mi t mk v k ds σ ij n j pn i n i n k v k ds Γ I Γ I σ ij n j t mi v k t mk ds σ ij n i n j p n k v k ds Γ I Γ I α S u i t mi v k t mk ds + p I n k v k ds. Γ I Γ I d m= d d m= m= As for the defective boundary conditions on Γ f,f, a total flow rate is specified on each of its pieces Γ f,k, but the traction vector on the boundary is not imposed. We follow Ervin et al. s approach [] and use a Lagrange multiplier method for treating the defective boundary conditions on Γ f,f which for a sufficiently smooth function is equivalent to the condition that the traction vector is a constant vector normal to the boundary surface. Thus the traction integral over Γ f,f is replaced by K f σ ij n j pn i v i ds Γ f,f k= β f,k v i n i ds = β f v i n i ds, 3.6 Γ f,k Γ f,f where in the last equality the Lagrange multiplier β f R K f is to be viewed as a piecewise constant function on Γ f,f, with β f n i representing the traction vector on the boundary. 6
7 Combining the above results, we have f i, v i Ωf = 2ν f Du Du ij Dv ij dx p v i dx Ω f Ω f x i τ B,i v i ds β f v i n i ds Γ f,n Γ f,f + d m= The averaged flow rates are imposed by α S u i t mi v k t mk ds + p I v i n i ds. Γ I Γ I K f K f γ f u i n i ds = γ f,k q f,k ds = γ f,k q f,k measγ f,k, 3.7 Γ f,f Γ f,k k= where γ f is a piecewise constant function on Γ f,f. The divergence-free equation 2.5 is weakly imposed using u i q = 0, q L 2 Ω f. x i We now introduce where τ B, v Γf,N Ω f Ω f a f u, v := 2ν f Du Du ij Dv ij dx + Ω f v i b f v, p := pdx, x i k= d m= K f b f,b v, β f := β f v i n i ds = β f,k v i n i ds, Γ f,f Γ f,k k= l f v := f i, v i Ωf + τ B, v Γf,N p I, v n ΓI, denotes a duality pairing between p I, v n ΓI a duality pairing between H 2 Γ I and H 2 Γ I. The weak formulation for the Stokes flow is: Given f X f, τ B p I H 2 Γ I, determine u, p, β f X f L 2 Ω f R K f Γ I α S u i t mi v k t mk ds, H d d, 2 Γ f,n and H 2 Γ f,n and such that H 2 Γ f,n d, qf R K f, and a f u, v b f v, p b f,b v, β f = l f v, v X f, 3.8 b f u, q + b f,b u, γ f = b f,b q f n f, γ f, q, γ f L 2 Ω f R K f. 3.9 Again, as remarked above, the term q f n f is to be interpreted for q f a piecewise constant function on Γ f,n. 7
8 We first establish an -sup condition before presenting our theorem on the existence, uniqueness and regularity of a solution. Let Z f,b := { v X f : b f,b v, β f = 0, β f R K } f, { } Z f,0 := v X f : v Γf,F = 0, Z f := { v X f : b f v, q + b f,b v, β f = 0, β f R K f, q L 2 Ω f }. Lemma. There exists a constant C > 0 such that β f R Kf, q L 2 Ω f b f v, q + b f,b v, β f sup C v X f v Xf q L 2 Ω f + β f Proof. Taking advantage of the fact that R K f is finite-dimensional, we easily conclude b f,b v, β f sup C > β f R K f v X f v Xf β f That is, for each j =,, K f, we find a v j X f such as v j k j Γ = 0 and f,k Γ f,j v j i n i ds 0. These v j s must exist because Γ k,j s are pair-wise disjointed. It is then a trivial exercise to construct a function v as a linear combination of these v j s for a given β f R K f to achieve b f,b v, β f C v Xf β f. It is well known that an -sup condition holds for b f v, q on Z f,0 L 2 Ω f. Since Z f,0 Z f,b, we know there exists an -sup condition for b f v, q on Z f,b L 2 Ω f. We now recall the fact that the continuity of the bilinear forms and individual -sup conditions imply the combined sup condition [] see also Theorem 8 and conclude the existence of an -sup condition for b f v, q + b f,b v, β f. Theorem. There exists a unique solution u, p, β f X f L 2 Ω f R K f satisfying In addition, there exists a constant C > 0 such that u Xf + d α S u i t mi 2 L 2 Γ I + p L 2 Ω f + β f m= C f X + τ B f H d + p I 2 Γ f,n H + q f 2 Γ I Proof. The existence and uniqueness of a solution follows from the continuity and strict monotonicity of a f, on Z f Z f, together with the -sup condition From the weak form 3.8, the -sup condition 3.20, and the assumption that ν f Du is 8
9 bounded from above, we have p L 2 Ω f + β f C sup v X f b f v, p + b f,b v, β f v Xf a f u, v l f v = C sup v X f v Xf C u Xf + f X f + τ B H 2 Γ f,n d + p I H 2 Γ I We choose v = u, q = p, and γ f = β f, and add together to obtain a f u, u = l f u + b f,b q f n f, β f = f i, u i Ωf + τ B, u Γf,N p I, u n ΓI + β f q f ds Γ f,f C u Xf f X f + τ B H 2 Γ f,n d + p I H 2 Γ I + C β f q f. Due to the assumption that measγ f,d > 0 and u B = 0, we know from Korn s inequality: which implies the strict positivity of a f u, u: We then have u 2 X f + u Xf = u H Ω f d C Du L 2 Ω f d d, a f u, u = 2ν Du Du ij Du ij dx + Ω f d m= d 2ν min Du 2 L 2 Ω f d d + C u 2 X f + α S u i t mi 2 L 2 Γ I C d m= The theorem follows from 3.23 and m= α S u i t mi 2 L 2 Γ I. d m= α S u i t mi 2 L 2 Γ I Γ S α S u i t mi u k t mk ds f 2 X + τ B 2 f H d + p I 2 2 Γ f,n H 2 + q f 2 + ɛ β f Γ I Remark: See [] For sufficiently smooth data, f, τ B, p I and solution u, p, with σn f = s n n f + s T, where s n = σn f n f and s T = σn f s n n f, 9
10 and τ B, v Γf,N := τ B,i v i ds, Γ f,n p I, v n f ΓI := p I v i n i ds, Γ I the unique solution of satisfies p + s n = β f,k and s T = 0 on Γ f,k, k =, 2,..., K f Thus the variational form corresponds to the boundary value problem with the additional constraint The reduced Stokes problem on the interface We let f X f, τ B H 2 Γ f,n d, and qf R K F be fixed given data as before, but we consider p I H 2 Γ I L 2 Γ I as a variable input. From Theorem, we know that there exists a unique solution u, p, β f X f L 2 Ω f R K F of that is a function of p I. We denote the solution as u p I, p p I, βf p I. Before studying the relationship between u p I and p I, we need a lemma, which we use below in the proof of Theorem 2: Lemma 2. If measγ f,n > 0, then there exists a constant C > 0 such that for all λ H 2 Γ I : sup v X f λ, v n ΓI v Xf C sup v Z f λ, v n ΓI v Xf. As a result, we have C λ H 2 Γ I sup v Z f λ, v n ΓI v Xf λ H 2 Γ I. Proof. We first note two -sup conditions: λ H 2 Γ I sup v X f λ, v n ΓI v Xf λ H 2 Γ I C > 0, 3.26 and where β f R K f b f,b v, β f sup C > 0, 3.27 v Z Γ I v Xf β f f { } Z Γ I f := v X f : λ, v n ΓI = 0, λ H 2 ΓI. Inequality 3.26 is well known, and 3.27 can be established by taking advantage of the finite dimensionality of R K f as we did for Lemma. Theorem 8 and imply λ H 2 Γ I, β f R K f λ, v n ΓI + b f,b v, β f sup C > v X f v Xf λ H + β f 2 Γ I 0
11 As measγ f,n > 0, we can show the following -sup condition by applying an argument similar to the one used to prove Lemma : q L 2 Ω f sup v Z b f,b,γ I f b f v, q v Xf q L 2 Ω f C > 0, 3.29 where { } Z b f,b,γ I f := v X f : λ, v n ΓI + b f,b v, β f = 0, λ H 2 ΓI, β f R K f. Now Theorem 8 and imply λ H 2 Γ I, q L 2 Ω f, β f R K f Applying Theorem 8 again, we see that 3.30 leads to λ, v n ΓI + b f v, q + b f,b v, β f sup C > v X f v Xf λ H + q 2 Γ I L 2 Ω f + β f λ H 2 Γ I sup v Z f λ, v n ΓI v Xf λ H 2 Γ I C > 0, or equivalently, λ H 2 Γ I C sup v Z f λ, v n ΓI v Xf, λ H 2 ΓI. Theorem 2. There exists a constant C > 0 such that for all λ H 2 Γ I and µ H 2 Γ I, If measγ f,n > 0, we further have u λ u µ Xf C λ µ H 2 Γ I. 3.3 C λ µ H 2 Γ I u λ u µ Xf C λ µ H 2 Γ I Proof. To show 3.3, we pick u = u λ, v = u µ, p = p λ, and β f = βf λ in 3.8, and pick u = u µ, q = p λ, and γ f = βf λ in 3.9, and then add them together to obtain The strict monotonicity of a f, then implies a f u λ, u µ = l f u µ + b f,b qf n f, β f λ. C u λ u µ 2 X f a f u λ, u λ u µ a f u µ, u λ u µ = λ µ, u λ u µ n f ΓI which yields the desired bound in 3.3. C u λ u µ Xf λ µ H 2 Γ I,
12 To show the first inequality in 3.32, we consider a f u λ, v + a f u µ, v + b f v, p λ p µ + b f,b v, β f λ βf µ = λ µ, v n f ΓI. For v Z f X f, we have λ µ, v n f ΓI = a f u λ, v + a f u µ, v C u λ u µ Xf v Xf Consequently, using Lemma 2 and 3.33 λ µ H 2 Γ I = sup q H 2 Γ I λ µ, q ΓI q H 2 ΓI = sup v X f λ µ, v n f ΓI v n H 2 ΓI C sup v Z f λ µ, v n f ΓI v n H 2 ΓI C u λ u µ Xf. If the normal velocity component q I H 2 Γ I is given on the interface Γ I, rather than the pressure, we can solve for the pressure by seeking λ H 2 Γ I such that with the form A f, defined by A f λ, µ = µ, q I ΓI µ H 2 ΓI, 3.34 A f λ, µ := µ, u f λ n f ΓI. We now prove a few properties of A f,. Theorem 3. A f λ, µ is a nonlinear functional of λ and a linear functional of µ. A f, is continuous on H 2 Γ I H 2 Γ I. In addition, A f, is monotone on H 2 Γ I H 2 Γ I. If measγ f,n > 0, A f, is strictly monotone on H 2 Γ I H 2 Γ I. Proof. It is clear that A f λ, µ is linear in µ. Its continuity follows from Theorem 2. We now show the strict monotonicity assuming measγ f,n > 0. From the weak formulation 3.8, we know A f λ, µ = µ, u f λ n f ΓI = a f u µ, u λ b f u λ, p µ b f,b u λ, βf µ f, u λ Ωf τ B, u f λ Γ f,n = a f u µ, u λ b f,b qf n f, βf µ f, u λ Ωf τ B, u f λ Γ f,n. We then obtain A f λ, λ µ A f µ, λ µ = A f λ, λ A f λ, µ A f µ, λ + A f µ, µ = a f u λ, u λ u µ a f u µ, u λ u µ C u λ u µ 2 H Ω f d C λ µ 2 H 2, Γ I 2
13 where we used Theorem 2 to obtain the last inequality above. If measγ f,n = 0, the last inequality above may fail to hold, and we then have only monotonicity not strict monotonicity. Theorem 4. If measγ f,n > 0, then there is a unique solution λ H 2 Γ I to the reduced problem 3.34 for any given q I H 2 Γ I. Proof. This theorem follows directly from the continuity and the strict monotonicity of A f,. 3.3 Darcy flow We next restrict our attention to the Darcy flow problem, assuming that the pressure is specified on the interface. Specifically, we consider together with: p = p I, on Γ I. Multiplying 2. by a smooth vector function v, with v n p Γp,N = 0, integrate over Ω p and apply the divergence theorem we obtain f i, v i Ωp = ν p u R ij u j + p v i dx 3.35 Ω p x i = ν p u R ij u j v i dx p v i dx + pv i n i dx Ω p Ω p x i Ω p = ν p u R ij u j v i dx p v i dx + p B v i n i ds Ω p Ω p x i Γ p,d + pv i n i ds + p I v i n i ds. Γ p,f Γ I In part, to incorporate the specified flow rate constraints into the formulation, for each Γ p,k we replace Γ p,k pv i n i ds by Γ p,k pv i n i ds equivalently β p,k v i n i ds, β p,k R, k =,..., K p, Γ p,k pv i n i ds Γ p,f β p v i n i ds, Γ p,f 3.36 where β p R K p in 3.36 is interpreted as a piecewise constant function on Γ p,f. The averaged flow rates are imposed by K p K p γ p u i n i ds = γ p,k q p,k ds = γ p,k q p,k measγ p,k. Γ p,f Γ p,k k= Similar to the previous treatment, the divergence-free equation 2.2 is weakly imposed using u i q = 0, q L 2 Ω p. Ω p x i 3 k=
14 For u Hdiv, Ω p, u n p H /2 Ω p. For p H /2 Γ I duality pairing does not define u n p, p I ΓI, as u n p acts on functions in H /2 Ω p. Following the work of Galvis and Sarkis [4], Lemma 2., given Γ s Ω p, r H /2 Γ s, let E /2 Γ s r H /2 Ω p denote the extension of r to Ω p. Then, for f H /2 Ω p, we denote f, r Γs := f, E /2 Γ s r Ωp. Also, as given in [4], for f H /2 Ω p, f Γs = 0 is defined as f, E /2 00,Γ s w Ωp = 0, for all w H /2 00 Γ s, where E /2 00,Γ s w denotes the extension by 0 of w to Ω p \Γ s. We now introduce X p = { } v : v Hdiv, Ω p, v n Γp,N = 0, a p u, v := ν u R ij u j v i dx, Ω p v i b p v, p := pdx, x i Ω p b p,b v, β p := v n p, β p Γp,F, l p v := f i, v i Ωp v n p, p B Γp,D v n p, p I ΓI. The weak formulation for the Darcy flow is: Given f X p, p B H 2 Γ p,d, q p R K p, and p I H 2 Γ I, determine u, p, β p X p L 2 Ω p R K p such that a p u, v b p v, p b p,b v, β p = l p v, v X p, 3.37 b p u, q + b p,b u, γ p = b p,b q p n p, γ p, q, γ p L 2 Ω p R K p Before presenting our results, we introduce Z p,b := { v X p : b f,b v, β p = 0, β p R K } p, { } Z p,0 := v X p : v n Γp,F = 0, Z p := { v X p : b p v, q + b p,b v, β p = 0, β p R Kp, q L 2 Ω f }. Lemma 3. There exists a constant C > 0 such that β p R Kp, q L 2 Ω f b p v, q + b p,b v, β p sup C v X p v Xp q L 2 Ω f + β p Proof. Similar to our previous arguments, we can utilize the finite dimensionality of R Kp the following -sup condition: to obtain β p R K p b p,b v, β p sup v X p v Xp β p C > 0. 4
15 It is well known that an -sup condition holds for b p v, q on Z p,0 L 2 Ω p. Since Z p,0 Z p,b, we know there exists an -sup condition for b p v, q on Z p,b L 2 Ω p. Theorem 8 then implies the existence of an -sup condition for b p v, q + b p,b v, β p. Theorem 5. There exists a unique solution u, p, β p X p L 2 Ω p R Kp satisfying In addition, there exists a constant C > 0 such that u Xp + p L 2 Ω p + β p C f X p + p B + p H 2 Γp,D I + q H 2 ΓI p Proof. The existence and uniqueness of the solution follows from the continuity and strict monotonicity of a p, on Z p Z p and the -sup condition From the weak formulation 3.37, the -sup condition 3.39 and the assumption that ν p u and R are both bounded from above, we have p L 2 Ω p + β b p v, p + b p,b v, β p p C sup v X p v Xp a p u, v l p v = C sup v X p v Xp C u Xp + f X p + p B + p H 2 Γp,D I H ΓI We pick q = u and γ p = 0 in 3.38 to obtain u 2 L 2 Ω p = b p u, u = 0, which implies that u Xp = u Hdiv,Ωp = u L 2 Ω p d. To see the regularity inequality, we pick v = u, q = p, and γ p = β p, and add together to obtain a p u, u = l p u + b p,b q p n p, β p = f i, u i Ωp u n p, p B Γp,D u n p, p I ΓI q p n p, β p Γp,F C f 2Xp + p B 2H 2 + p I 2H Γp,D 2 + q p 2 + ɛ u 2 ΓI X p + β p From the assumptions that the viscosity ν p u is bounded from below and that R is strictly positive definite, we have a p u, u = ν p u R ij u j v i dx C u L 2 Ω p d Ω p Now the theorem follows from Remark: See [] For sufficiently smooth data, f, p B, p I and solution u, p the unique solution of satisfies p = β p,k on Γ p,k, k =, 2,..., K p Thus the variational form corresponds to the boundary value problem with the additional constraint
16 3.4 The reduced Darcy problem on the interface Let f X p, p B H 2 Γ p,d, and q p R K p be fixed given data as before. For each p I H 2 Γ I H 2 Γ I, we have a unique solution u p I, p p I, βpp I X p L 2 Ω p R Kp that satisfies Theorem 6. There exists a constant C > 0such that If measγ p,d > 0, we further have u λ u µ Xp C λ µ H 2 ΓI C λ µ H 2 Γ I u λ u µ Xp C λ µ H 2 ΓI Proof. To show the first inequality in 3.46, we consider a p u λ, v + a p u µ, v + b p v, p λ p µ + b p,b v, β p λ β pµ For v Z p X p, C u λ u µ Xp v Xp = C u λ u µ L 2 Ω p d v L 2 Ω p d a p u λ, v + a p u µ, v = v n p, λ µ ΓI. = v n p, λ µ ΓI. Thus, with the assumption that measγ p,d > 0 employing a result analogous to Lemma 2 for X p and Z p, we have λ µ H 2 ΓI = sup q H 2 Γ I q, λ µ ΓI q H 2 Γ I C sup v X p v n p, λ µ ΓI v n H 2 Ω p C sup v Z p v n p, λ µ ΓI v n H 2 Γ I C u λ u µ Xp. Similar to the proof for Theorem 2, 3.45 and the second inequality in 3.46 can be shown by using the weak formulation and the strict monotonicity of a p,. One difference here is that we need to use the fact that the Darcy velocity is divergence free. If p I is an unknown but the normal component of Darcy velocity q I H 2 Γ I is given instead, we can solve for the pressure on Γ I by seeking λ H 2 Γ I such that A p λ, µ = q I µds, µ H 2 ΓI, 3.47 Γ I 6
17 with the form A p, defined by We now prove a few properties of A p,. A p λ, µ := u pλ n p, µ ΓI. Theorem 7. A p λ, µ is a nonlinear functional of λ and a linear functional of µ. A p, is continuous on H 2 Γ I H 2 Γ I. Moreover, A p, is monotone on H 2 Γ I H 2 Γ I. If measγ p,d > 0, A p, is strictly monotone on H 2 Γ I H 2 Γ I. Remark. Unlike the space on Γ I for Stokes flow, λ and µ here need to be defined in H 2 Γ I, not H 2 Γ I. Proof. The continuity follows from Theorem 6. To see the strict monotonicity under the assumption of measγ p,d > 0, we note A p λ, λ µ A p µ, λ µ = a p u λ, u λ u µ a p u µ, u λ u µ C u λ u µ 2 L 2 Ω p d = C u λ u µ 2 Hdiv,Ω p C λ µ 2 H 2 Γ I, where we have used Theorem 6 for the last inequality above. similarly as that in Theorem 3. The rest of this theorem follows Theorem 8. If measγ p,d > 0, then there is a unique solution λ H 2 Γ I to the reduced problem 3.47 for any given q I H 2 Γ I. Proof. This theorem follows directly from the continuity and the strict monotonicity of the form A p,. 3.5 Coupled system For convenience, in the following analysis we assume that measγ p,d > 0, which as we will see, gives the uniqueness of the pressure solution. We will remark later on the modification of the analysis for the case of measγ f,n = measγ p,d = 0. We couple the fluid flow through the two regions using continuity of the flux on Γ I. We define From the continuity of the normal velocity A λ, µ := A f λ, µ + A p λ, µ. [u n] ΓI = u f n f + u p n p ΓI = 0, we know that the coupled system can be formulated as: Determine λ H 2 Γ I such that We now list a few properties of A,. A λ, µ = 0, µ H 2 ΓI
18 Theorem 9. A λ, µ is a nonlinear functional of λ and a linear functional of µ. Under the assumption that measγ p,d > 0, A, is continuous and strictly monotone on H 2 Γ I H 2 Γ I. Proof. It is clear that A λ, µ is linear in µ. Continuity follows from the continuity of its parts: Af Ap A λ, µ A λ, µ λ, µ A f λ, µ + λ, µ A p λ, µ C λ λ λ µ H 2 Γ I H + C 2 λ µ Γ I H 2 ΓI H 2 Γ I C λ λ µ H 2 ΓI H 2, Γ I and A λ, µ A λ, µ A f λ, µ A f λ, µ + A p λ, µ A p λ, µ C λ H 2 ΓI µ µ H 2 Γ I. The strict monotonicity also follows from the individual monotonicities. That is, A f λ, λ µ A f µ, λ µ 0 and A p λ, λ µ A p µ, λ µ C λ µ 2 H together imply A λ, λ µ 2 Γ I A µ, λ µ C λ µ 2 H 2 Γ I. Theorem 0. We assume that measγ p,d > 0. For any given q I H 2 Γ I, there is a unique solution λ H 2 Γ I to the reduced problem Proof. This theorem follows directly from the continuity and the strict monotonicity of the form A,. Theorem. We assume that measγ p,d > 0. Let p I = λ be the solution to the reduced interface problem 3.48, and u f p I, p f p I, β f p I, u pp I, p pp I, and β pp I be the subdomain solutions to , and We have the following regularity result: p I + H 2 u ΓI f p I Xf + u p p I Xp + p p I L 2 Ω + β f p I + β p p I C f f X + f p f X + τ B p H 2 d + p B + q Γ f,n H 2 Γp,D f + q p. Proof. From the weak formulation 3.48, we have A p I, p I = 0. Using this together with the continuity and strict monotonicity of the form A,, we conclude p I 2 H 2 Γ I = p I 0 2 H 2 Γ I C A p I, p I 0 A 0, p I 0 = CA 0, p I = C u f 0 n f, p I ΓI C u p0 n p, p I ΓI ɛ p I 2 H 2 Γ I + C u f 0 n f 2 H 2 Γ I + C u p 0 n p 2 H 2 Γ I ɛ p I 2 H 2 Γ I + C u f 0 2 X f + C u p 0 2 X p. 8
19 Consequently, we have p I H 2 ΓI C u f 0 Xf + C u p0 Xp. Theorems and 5 imply u f 0 Xf C f f X + τ B f H d + q f 2 Γ f,n and u p 0 Xp C f p X + p B + q p H 2 Γp,D p. This Theorem follows by another application of Theorems and 5. Remark. If measγ f,n = measγ p,d = 0, the pressure solutions in the two stand-alone subdomain problems are still unique. But the pressure solution to the interface problem 3.48 is no longer unique; instead it is unique up to an additive constant. As a result, the two subdomain pressure solutions to the coupled system are unique up to an additive constant. In particular, we see that A f, and A p,, and consequently A,, lose their strict monotonicity, which leads to the nonuniqueness of solutions. To see that the pressure solution actually exists uniquely up to an additive constant under measγ f,n = measγ p,d = 0, we replace the original spaces for p f, p p and p I by L 2 Ω f /R, L 2 Ω p /R, and H 2 Γ I /R, respectively. With the new spaces, one can show that A f, and A p,, and consequently A,, recover their strict monotonicity under measγ f,n = measγ p,d = 0. All other theorems above then follow under the modified pressure spaces. 4 Finite Element Approximations We discretize the coupled nonlinear Stokes-Darcy system 3.48, , and with finite element approximations. We use conforming approximating spaces: X f,h X f, M f,h L 2 Ω f, X p,h X p, M p,h L 2 Ω p, L h H 2 ΓI. Here X f,h and M f,h denote velocity and pressure spaces typically used for Stokes fluid flow approximations, for example, the Taylor-Hood spaces. X p,h and M p,h denote typical velocity and pressure spaces in mixed finite element methods for Darcy flow, such as the Raviart-Thomas spaces. In this section, we again assume that measγ p,d > 0 for the uniqueness of the pressure solution. We first list a few assumptions on our finite element spaces. Assumption. We assume that there exists an -sup condition for the Stokes flow approximation: β f R Kf, q h M f,h b f v h, q h + b f,b v h, β f sup C > v h X f,h v h Xf q h L 2 Ω f + β f Remark. The above condition holds for the Taylor-Hood spaces on a quasi-uniform mesh of triangles or tetrahedra see [2, ]. 9
20 Assumption 2. We assume that there exists an -sup condition for the Darcy flow approximation: β p R Kp, q h M p,h b p v h, q h + b p,b v h, β p sup C > v h X p,h v h Xp q h L 2 Ω + β p p Remark. The above assumption holds for the Raviart-Thomas spaces on a quasi-uniform mesh of triangles or tetrahedra see [2, ]. Assumption 3. Compatibility condition of the mixed spaces for Darcy flow We assume that there exists a swapped -sup condition for the Darcy flow approximation: v h X p,h b p v h, q h sup q h M p,h v h L 2 Ω p q h L 2 Ω p C > Remark. The above assumption holds whenever X p,h M p,h, which is satisfied by commonly used mixed finite element spaces for example, the Raviart-Thomas spaces. Analogous to the analysis for the continuous problems, we define two null spaces: Z f,h := { v h X f,h : b f v h, q h + b f,b v h, γ f = 0, γ f R K f, q h M f,h }, Z p,h := { v h X p,h : b p v h, q h + b p,b v h, γ p = 0, γ p R K p, q h M p,h }, and introduce two affine sets: Z A f,h := { v h X f,h : b f v h, q h + b f,b v h, γ f = b f,b q f n f, γ f, γ f R K f, q h M f,h }, Z A p,h := { v h X p,h : b p v h, q h + b p,b v h, γ p = b f,b q p n p, γ p, γ p R K p, q h M p,h }. Assumption 4. Mortar compatibility condition for Stokes flow We assume that there exists an -sup condition between the mortar space and the Stokes flow approximation spaces: v h n f, λ h ΓI sup λ h L h v h Z f,h v h Xf λ h H 2 Γ I C Assumption 5. Mortar compatibility condition for Darcy flow We assume that there exists an -sup condition between the mortar space and the Darcy flow approximation spaces: v h n p, λ h ΓI sup λ h L h v h Z p,h v h Xp λ h H 2 ΓI C Remark. Roughly speaking, 4.52 specifies that the normal-velocity component restriction of Z f,h on Γ I should be at least as dense as the mortar space L h. Similarly, 4.53 specifies that the normalvelocity component restriction of Z p,h on Γ I should be at least as dense as the mortar space L h. We remark that it is sufficient but not necessary to require both mortar compatibility conditions 4.52 and 4.53 for the existence and uniqueness of a solution to the coupled system; either one of the two conditions is sufficient. We assume both compatibility conditions here for convenience. Following from the above -sup assumptions and Theorem 8, we have the combined compatibility conditions: 20
21 Lemma 4. There exist a constant C > 0 such that λ h L h, q h M f,h, β f R K f λ h L h, q h M p,h, β p R Kp λ h, v h n f ΓI + b f v h, q h + b f,b v h, β f sup C > 0, 4.54 v h X f,h v h Xf λ h H + q h 2 Γ I L 2 Ω f + β f v h n p, λ h ΓI + b p v h, q h + b p,b v h, β p sup C > v h X p,h v h Xp λ h + q H 2 ΓI h L 2 Ω p + β p 4. Stokes flow approximation We first restrict out attention to the finite element discretization of Stokes flow, assuming that p I,h L h H 2 Γ I H 2 Γ I is given. Define l f,h v := f i, v i Ωf + τ B,i v i ds Γ f,n p I,h v i n i ds. Γ I The weak formulation is: Given f X f, τ N u h, p h, β f,h X f,h M f,h R K f such that H 2 Γ f,n d, qf R K F, and p I,h L h, determine a f u h, v h b f v h, p h b f,b v h, β f,h = l f,h v h, v h X f,h, 4.56 b f u h, q h + b f,b u h, γ f = b f,b q f n f, γ f,, q h, γ f M f,h R K f Theorem 2. There exists a unique solution u h, p h, β f,h X f,h M f,h R K f satisfying In addition, there exists a constant C > 0 such that u u h Xf + p p h L 2 Ω f + β f β f,h C u û h Xf + p p h L û h X f,h p h M 2 Ω f f,h + C p I,h p I H 2 Γ I Proof. The existence and uniqueness of the solution follows from the continuity and strict monotonicity of a f, on Z f,h Z f,h, together with the discrete -sup condition Subtracting from , respectively, we obtain the following error equations: a f u h, v h = a f u, v h + b f v h, p h p + b f,b v h, β f,h β f + l f,h v h l f v h, v h X f,h, 4.59 b f u h u, q h + b f,b u h u, γ f = 0, q h, γ f M f,h R K f We note that the right-hand side of 4.59 contains the error propagated from the interface approximation: l f,h v h l f v h = p I,h p I, v h n ΓI. 4.6 We pick û h Z A f,h X f,h, p h M f,h, β f,h R K f, and note that b f u h û h, p h p h + b f,b u h û h, β f,h β f,h =
22 With v h = u h û h Z f,h X f,h, using 4.59, 4.6 and 4.62, we obtain a f u h, u h û h a f û h, u h û h = a f u, u h û h a f û h, u h û h + b f u h û h, p h p + b f,b u h û h, β f,h β f p I,h p I, u h û h n ΓI Because of the strict monotonicity, we have u h û h Xf ɛ u h û h 2 X f + C u û h 2 X f + C p p h 2 L 2 Ω f + C β f β 2 f,h + C pi,h p I 2 H 2. Γ I a f u h, u h û h a f û h, u h û h C u h û h 2 X f. Since û h, p h, and β f,h can be arbitrarily chosen from their corresponding spaces, we have C β f β f,h = C u û h û h Z A Xf + f,h p p h L p h M 2 Ω f + f,h + C p I,h p I H 2 Γ I u û h û h Z A Xf + p p h L p f,h h M 2 Ω f f,h β f,h R K f + C p I,h p I H 2 Γ I We proceed to lift û h from Z A f,h to X f,h in the above estimate. From the -sup condition 4.49, we know that there exists an operator Π h : X f X f,h such that b f u Π h u, q h + b f,b u Π h u, γ f = 0, q h, γ f M f,h R K f and Π h u Xf C u Xf. For any given ũ h X f,h, we now define û h := ũ h Π h ũ h u. Using 4.60 and that fact that u satisfies 3.9, we easily verify that û h Z A f,h. We then have u û h Xf u ũ h Xf + ũ h û h Xf = u ũ h Xf + Π h ũ h u Xf + C u ũ h Xf, which implies u û h û h Z A Xf + C u ũ h Xf ũ f,h h X f,h The triangle inequality together with 4.63 and 4.64 yield u u h Xf C u û h Xf + p p h L û h X f,h p h M 2 Ω f f,h + C p I,h p I H 2 Γ I
23 To get the error estimate for p h and β f,h, we use the -sup condition 4.49 and the error equation Picking p h M f,h and β f,h R K f, we have v h, β f,h β f,h p h p h L 2 Ω f + βf,h β b f v h, p h p h + b f,b f,h C sup v h X f,h v h Xf implying b f v h, p h p + b f,b v h, β f,h β f C sup v h X f,h v h Xf v h, β f β f,h b f v h, p p h + b f,b + C sup v h X f,h v h Xf a f u h, v h a f u, v h + p I,h p I, v h n ΓI = C sup v h X f,h v h Xf v h, β f β f,h b f v h, p p h + b f,b + C sup v h X f,h v h Xf p p h L 2 Ω p + β f β f,h C u u h Xf + C C u u h Xf + C p I,h p I H 2 Γ I + C p p h L 2 Ω f + C βf β f,h, + C β f,h R K f The estimate 4.58 follows from 4.65 and p h M f,h p p h L 2 Ω f β f β f,h + C pi,h p I H 2 Γ I = C u u h Xf + C p h M f,h p p h L 2 Ω f + C p I,h p I H 2 Γ I Darcy flow approximation We now consider the finite element discretization for the Darcy flow, assuming that p I,h L h H 2 Γ I is given. We define l p,h v := f i, v i Ωp v h n p, p B Γp,D v h n p, p I,h ΓI. The weak formulation is: Given f X p, p B H 2 Γ p,d, q p R Kp, and p I,h L h, determine u h, p h, β p,h X p,h M p,h R K p such that a p u h, v h b p v h, p h b p,b v h, β p,h = l p,h v h, v h X p,h, 4.67 b p u h, q h + b p,b u h, γ p = b p,b q p n p, γ p, q h, γ p M p,h R Kp
24 Theorem 3. There exists a unique solution u h, p h, β p,h X p,h M p,h R Kp satisfying In addition, there exists a constant C > 0 such that u u h Xp + p p h L 2 Ω p + β p β p,h C u û h Xp + û h X p,h p h M p,h p p h L 2 Ω p + C p I,h p I H 2 ΓI Proof. The existence and uniqueness of the solution follows from the continuity and strict monotonicity of a p, on Z p,h Z p,h, together with the discrete -sup condition Subtracting from , respectively, we obtain the error equations: a p u h, v h a p u, v h b p v h, p h p b p,b v h, β p,h β p = v h n p, p I,h p I ΓI, v h X p,h, 4.70 b p u h u, q h + b p,b u h u, γ p = 0, q h, γ p M p,h R K p. 4.7 We now take û h Z A p,h X p,h, p h M p,h, β p,h R K p, and set v h = u h û h Z p,h X p,h in 4.70 to obtain a p u h, u h û h a p û h, u h û h = a p u, u h û h a p û h, u h û h + b p u h û h, p h p + b p,b u h û h, β p,h β p u h û h n p, p I,h p I ΓI Due to the strict monotonicity of a p,, we have ɛ u h û h 2 X p + C u û h 2 L 2 Ω p d + C p p h 2 L 2 Ω p + C β p β 2 p,h + C pi,h p I 2 H 2. Γ I a p u h, u h û h a p û h, u h û h C u h û h 2 L 2 Ω p d. Since the û h, p h, and β p,h can be arbitrarily chosen from their corresponding spaces, we have u u h L 2 Ω p u d û h L 2Ωp d + û h u h L 2Ωp d C C û h Z A p,h u û h û h Z A L 2 Ω p d + p,h + C p I,h p I H 2 ΓI + ɛ û h Z A p,h p p h L p h M 2 Ω + p p,h u h û h Xp u û h û h Z A Xp + p p h L p p,h h M 2 Ω p p,h β p,h R K p β p β p,h + C p I,h p I H 2 ΓI + ɛ u u h Xp Now we use the swapped -sup condition 4.5 and 4.68 to obtain u h L 2 Ω p C sup b p u h, q h q h M p,h q h L 2 Ω p 24 = 0,
25 which then implies As a consequence of 4.72 and 4.73, we have u u h Xp C u h u L 2 Ω p = u û h û h Z A Xp + p p h L p p,h h M 2 Ω p p,h + C p I,h p I H 2 ΓI Using the -sup condition 4.50 and the fact that u satisfies 3.38, we now lift û h from Z A p,h to X p,h by applying an argument similar to the one we have used before for lifting the velocity in the Stokes flow region to obtain: Now 4.74 and 4.75 yield u u h Xp C u û h û h Z A Xp C u ũ h Xp ũ p,h h X p,h u û h Xp + û h X p,h p p h L p h M 2 Ω p + C p I,h p I p,h H 2 ΓI To get the error estimate for p h and β p,h, we use the -sup condition 4.50 and the error equation 4.70 by picking p h M p,h and β p,h R K p similar to what we have done for Stokes flow: v h, β p,h β p,h p h p h L 2 Ω p + βp,h β b p v h, p h p h + b p,b p,h C sup v h X p,h v h Xp implying b p v h, p h p + b p,b v h, β p,h β p C sup v h X p,h v h Xp v h, β p β p,h b p v h, p p h + b p,b + C sup v h X p,h v h Xp a p u h, v h a p u, v h + v h n p, p I,h p I ΓI = C sup v h X p,h v h Xp v h, β p β p,h b p v h, p p h + b p,b + C sup v h X p,h v h Xp p p h L 2 Ω p + β p β p,h C u u h Xp + C C u u h Xp + C p I,h p I H 2 ΓI + C p p h L 2 Ω f + C βp β p,h, + C β p,h R K p The estimate 4.69 follows from 4.76 and p h M p,h p p h L 2 Ω p β p β p,h + C pi,h p I H 2 ΓI = C u u h Xp + C p h M p,h p p h L 2 Ω p + C p I,h p I H 2 ΓI
26 4.3 Discretization on the interface For each p I,h L h H 2 Γ I H 2 Γ I, we have a unique finite element solution u f,h p I,h, p f,h p I,h, βf,h p I,h X f,h M f,h R K f satisfying , and a unique finite element solution u p,h p I,h, p p,h p I,h, βp,h p I,h X p,h M p,h R K p satisfying We now define A f,h λ h, µ h := µ h, u f,h λ h n f ΓI, A p,h λ h, µ h := u p,h λ h n p, µ h ΓI, A h λ h, µ h := A f,h λ h, µ h + A p,h λ h, µ h. The discretized interface problem for the coupled system can be formulated as: Determine λ h L h H 2 Γ I such that A h λ h, µ h = 0, µ L h Theorem 4. There exists a constant C such that C λ h µ h H 2 Γ I u f,h λ h u f,h µ h Xf C λ h µ h H 2 Γ I, 4.79 C λ h µ h H 2 ΓI u p,h λ h u p,h µ h Xp C λ h µ h H 2 ΓI Proof and 4.80 can be shown similarly as we did for Theorems 2 and 6, except that we now need to replace the continuous -sup conditions by their discrete counterparts, i.e. the mortar compatibility conditions 4.52 and 4.53 for the lower bounds. We point out that the proof for the upper bound part of 4.80 uses the property that u h is divergence free in the Darcy region. We now prove a few properties of A h,. Theorem 5. A h λ h, µ h is a nonlinear functional of λ h and a linear functional of µ h. A h, is continuous and strictly monotone on L h L h. Proof. It is clear that A h λ, µ is linear in µ. To see the continuity, we apply Theorem 4: Af,h Ap,h A h λ h, µ h A h λh, µ h λ h, µ h A f,h λh, µ h + λ h, µ h A p,h λh, µ h C u f,h λ h u f,h λ h µ h Xf H 2 Γ I + C u p,h λ h u p,h λ h µ h Xp H 2 ΓI C λ h λ H h µ h 2 Γ I H + C λh λ H 2 h µ Γ I 2 h ΓI H 2 ΓI C λ h λ H h µ 2 h, ΓI H 2 ΓI 26
27 and A h λ h, µ h A h λ h, µ h A f,h λ h, µ h A f,h λ h, µ h + A p,h λ h, µ h A p,h λ h, µ h C λ h H 2 ΓI µ h µ h H 2 ΓI. Like its continuous counterpart in Theorem 3, we have the following strict monotonicity A f,h λ h, λ h µ h A f,h µ h, λ h µ h = a f u f,h λ h, u f,h λ h u f,h µ h a f u f,h µ h, u f,h λ h u f,h µ h C u f,h λ h u f,h µ h 2 X f C λ h µ h 2 H 2 Γ I. Similarly, we have A p,h λ h, λ h µ h A p,h µ h, λ h µ h C λ h µ h 2 H 2 Γ I. The strict monotonicity of A h, then follows directly from these two inequalities. Theorem 6. For any given q I H 2 Γ I, there is a unique solution λ h L h to the reduced problem Proof. This theorem directly follow from the continuity and the strict monotonicity of the form A h,. Theorem 7. Let p I,h = λ h be the solution to the discretized interface problem Let u f,h, p f,h, β f,h and u p,h, p p,h, β p,h be the subdomain solutions to and supplied with p I,h. Then, there exists a constant C such that p I p I,h H 2 u + f u f,h ΓI Xf + p f p f,h L 2 Ω f + β f β f,h + u p u p,h Xp + p p p p,h L 2 Ω p + β p β p,h C p I p I,h + C p I,h L h H 2 u f û h ΓI û h X Xf + f,h + C u p û h Xp + p p p h û h X p,h p h M L 2 Ω p. p,h Proof. Let p I,h L h. 4.78, we have p f p h p h M L 2 Ω f f,h From the weak formulation 3.48 and its finite dimensional counterpart A p I, p I,h p I,h = 0, A h p I,h, p I,h p I,h = 0. These two orthogonality conditions together with the strict monotonicity of A, yield C p I,h p I,h 2 H 2 Γ I A p I,h, p I,h p I,h A p I,h, p I,h p I,h = A p I,h, p I,h p I,h A h p I,h, p I,h p I,h + A p I, p I,h p I,h A p I,h, p I,h p I,h
28 The first two terms on the right-hand side of 4.8 contain the error propagated from the two subdomains and can be bounded using Theorems 2 and 3: A p I,h, p I,h p I,h A h p I,h, p I,h p I,h = p I,h p I,h, u f p I,h u f,h p I,h n f ΓI u pp I,h u p,h p I,h n p, p I,h p I,h ΓI C u f p I,h u f,h p I,h Xf p I,h p I,h H 2 + C u Γ I p p I,h u p,h p I,h Xp p I,h p I,h H 2 ΓI C u f û h û h X Xf + p f p h f,h p h M L 2 Ω f p I,h p I,h f,h H 2 Γ I + C u p û h Xp + p p p h û h X p,h p h M L 2 Ω p p I,h p I,h p,h H 2 ΓI ɛ p I,h p I,h 2 H 2 + C u f û h 2 Γ I û h X X f + p f p h 2 f,h p h M L 2 Ω f f,h + C u p û h 2 X û h X p + p p p h 2 p,h p h M L 2 Ω p. p,h The last two terms in 4.8 contain the error from the interface discretization: A p I, p I,h p I,h A p I,h, p I,h p I,h = p I,h p I,h, u f p I u f p I,h n f ΓI u pp I u pp I,h n p, p I,h p I,h ΓI C u f p I u f p I,h Xf p I,h p I,h H 2 Γ I + C u p p I u p p I,h Xp p I,h p I,h H 2 ΓI C p I p I,h H 2 Γ I p I,h p I,h H 2 Γ I + C p I p I,h H 2 ΓI p I,h p I,h H 2 ΓI C p I p I,h 2 H 2 Γ I + ɛ p I,h p I,h 2 H 2 Γ I. This theorem follows by first noting the fact that p I,h can be chosen arbitrarily from L h and then applying a triangle inequality and Theorems 2 and 3. 5 Numerical Examples We consider a coupled Stokes-Darcy system on the domain [0, 2] [0, ], representing a twodimensional version of an industrial filtering application where a non-newtonian fluid passes through a filter to remove unwanted particulates. Here the flow of the fluid through the channel [0, ] [0, ] is coupled with its flow in the porous medium [, 2] [0, ]. We impose a defective boundary condition on the left boundary by specifying an low flux of one, and we impose a defective boundary condition on the right boundary of an outflow flux of one. Along the top and bottom boundaries, we impose a no-slip boundary condition for the Stokes flow and a no-flow boundary condition for the Darcy flow. We assume a Cross model for the fluid viscosity ν f in the Stokes region and another Cross model for the effective viscosity ν p in the Darcy region: ν f Du = ν f, + ν f,0 ν f, + K f Du 2 r, f ν p u = ν p, + ν p,0 ν p, + K p u 2 r p, 28
Stochastic 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 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 informationGeneralized Newtonian Fluid Flow through a Porous Medium
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
More informationApproximation of the Stokes-Darcy system by optimization
Approximation of the Stokes-Darcy system by optimization V.J. Ervin E.W. Jenkins Hyesuk Lee June 2, 202 Abstract A solution algorithm for the linear/nonlinear Stokes-Darcy coupled problem is proposed and
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 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 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 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 informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
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 informationMixed Hybrid Finite Element Method: an introduction
Mixed Hybrid Finite Element Method: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola
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 informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More informationA Two-grid Method for Coupled Free Flow with Porous Media Flow
A Two-grid Method for Coupled Free Flow with Porous Media Flow Prince Chidyagwai a and Béatrice Rivière a, a Department of Computational and Applied Mathematics, Rice University, 600 Main Street, Houston,
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 informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
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 informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More 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 informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationNotes for Expansions/Series and Differential Equations
Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated
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 informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationDivergence-conforming multigrid methods for incompressible flow problems
Divergence-conforming multigrid methods for incompressible flow problems Guido Kanschat IWR, Universität Heidelberg Prague-Heidelberg-Workshop April 28th, 2015 G. Kanschat (IWR, Uni HD) Hdiv-DG Práha,
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 informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
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 informationAppendix A Functional Analysis
Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)
More informationHybridized DG methods
Hybridized DG methods University of Florida (Banff International Research Station, November 2007.) Collaborators: Bernardo Cockburn University of Minnesota Raytcho Lazarov Texas A&M University Thanks:
More informationFINITE 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 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 informationCoupling Stokes and Darcy equations: modeling and numerical methods
Coupling Stokes and Darcy equations: modeling and numerical methods Marco Discacciati NM2PorousMedia 24 Dubrovnik, September 29, 24 Acknowledgment: European Union Seventh Framework Programme (FP7/27-23),
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 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 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 informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
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 informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationNAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS
NAVIER-STOKES EQUATIONS IN THIN 3D DOMAINS WITH NAVIER BOUNDARY CONDITIONS DRAGOŞ IFTIMIE, GENEVIÈVE RAUGEL, AND GEORGE R. SELL Abstract. We consider the Navier-Stokes equations on a thin domain of the
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
More informationOn the Ladyzhenskaya Smagorinsky turbulence model of the Navier Stokes equations in smooth domains. The regularity problem
J. Eur. Math. Soc. 11, 127 167 c European Mathematical Society 2009 H. Beirão da Veiga On the Ladyzhenskaya Smagorinsky turbulence model of the Navier Stokes equations in smooth domains. The regularity
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 informationOutline: 1 Motivation: Domain Decomposition Method 2 3 4
Multiscale Basis Functions for Iterative Domain Decomposition Procedures A. Francisco 1, V. Ginting 2, F. Pereira 3 and J. Rigelo 2 1 Department Mechanical Engineering Federal Fluminense University, Volta
More informationCartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components
Cartesian Tensors Reference: Jeffreys Cartesian Tensors 1 Coordinates and Vectors z x 3 e 3 y x 2 e 2 e 1 x x 1 Coordinates x i, i 123,, Unit vectors: e i, i 123,, General vector (formal definition to
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 informationIMA Preprint Series # 2143
A BASIC INEQUALITY FOR THE STOKES OPERATOR RELATED TO THE NAVIER BOUNDARY CONDITION By Luan Thach Hoang IMA Preprint Series # 2143 ( November 2006 ) INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY
More informationc 2003 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 40, No. 6, pp. 195 18 c 003 Society for Industrial and Applied Mathematics COUPLING FLUID FLOW WITH POROUS MEDIA FLOW WILLIAM J. LAYTON, FRIEDHELM SCHIEWECK, AND IVAN YOTOV Abstract.
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. # 02 Conservation of Mass and Momentum: Continuity and
More informationIntroduction to finite element exterior calculus
Introduction to finite element exterior calculus Ragnar Winther CMA, University of Oslo Norway Why finite element exterior calculus? Recall the de Rham complex on the form: R H 1 (Ω) grad H(curl, Ω) curl
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationSolving the curl-div system using divergence-free or curl-free finite elements
Solving the curl-div system using divergence-free or curl-free finite elements Alberto Valli Dipartimento di Matematica, Università di Trento, Italy or: Why I say to my students that divergence-free finite
More informationDivergence-free or curl-free finite elements for solving the curl-div system
Divergence-free or curl-free finite elements for solving the curl-div system Alberto Valli Dipartimento di Matematica, Università di Trento, Italy Joint papers with: Ana Alonso Rodríguez Dipartimento di
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
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 informationA Multiscale Mortar Multipoint Flux Mixed Finite Element Method
A Multiscale Mortar Multipoint Flux Mixed Finite Element Method Mary F. Wheeler Guangri Xue Ivan Yotov Abstract In this paper, we develop a multiscale mortar multipoint flux mixed finite element method
More informationRobust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems
Robust Monolithic - Multigrid FEM Solver for Three Fields Formulation Rising from non-newtonian Flow Problems M. Aaqib Afaq Institute for Applied Mathematics and Numerics (LSIII) TU Dortmund 13 July 2017
More informationA Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 UNCORRECTED PROOF
1 A Hybrid Discontinuous Galerkin Method 2 for Darcy-Stokes Problems 3 Herbert Egger 1 and Christian Waluga 2 4 1 Center of Mathematics, Technische Universität München, Boltzmannstraße 3, 85748 5 Garching
More informationAllen Cahn Equation in Two Spatial Dimension
Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic
More informationAnalysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods
Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp. 830-844 DOI: 10.408/aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded
More informationON THE QUASISTATIC APPROXIMATION IN THE STOKES-DARCY MODEL OF GROUNDWATER-SURFACE WATER FLOWS
ON THE QUASSTATC APPROXMATON N THE STOKES-DARCY MODEL OF GROUNDWATER-SURFACE WATER FLOWS Marina Moraiti Department of Mathematics University of Pittsburgh Pittsburgh, PA 156, U.S.A. Tel: +1 (41) 889-9945
More informationNONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS
NONSTANDARD NONCONFORMING APPROXIMATION OF THE STOKES PROBLEM, I: PERIODIC BOUNDARY CONDITIONS J.-L. GUERMOND 1, Abstract. This paper analyzes a nonstandard form of the Stokes problem where the mass conservation
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
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 informationMULTIPLE SOLUTIONS FOR THE p-laplace EQUATION WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 2006(2006), No. 37, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) MULTIPLE
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationThe oblique derivative problem for general elliptic systems in Lipschitz domains
M. MITREA The oblique derivative problem for general elliptic systems in Lipschitz domains Let M be a smooth, oriented, connected, compact, boundaryless manifold of real dimension m, and let T M and T
More informationR. M. Brown. 29 March 2008 / Regional AMS meeting in Baton Rouge. Department of Mathematics University of Kentucky. The mixed problem.
mixed R. M. Department of Mathematics University of Kentucky 29 March 2008 / Regional AMS meeting in Baton Rouge Outline mixed 1 mixed 2 3 4 mixed We consider the mixed boundary value Lu = 0 u = f D u
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 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 informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part IV Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More information10 Transfer Matrix Models
MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More 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 informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationProject Description: MODELING FLOW IN VUGGY MEDIA. Todd Arbogast, Mathematics
Project Description: MODELING FLOW IN VUGGY MEDIA Todd Arbogast, Mathematics Steve Bryant, Petroleum & Geosystems Engineering Jim Jennings, Bureau of Economic Geology Charlie Kerans, Bureau of Economic
More informationh(x) lim H(x) = lim Since h is nondecreasing then h(x) 0 for all x, and if h is discontinuous at a point x then H(x) > 0. Denote
Real Variables, Fall 4 Problem set 4 Solution suggestions Exercise. Let f be of bounded variation on [a, b]. Show that for each c (a, b), lim x c f(x) and lim x c f(x) exist. Prove that a monotone function
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 informationSolvability of the Incompressible Navier-Stokes Equations on a Moving Domain with Surface Tension. 1. Introduction
Solvability of the Incompressible Navier-Stokes Equations on a Moving Domain with Surface Tension Hantaek Bae Courant Institute of Mathematical Sciences, New York University 251 Mercer Street, New York,
More informationCartesian Tensors. e 2. e 1. General vector (formal definition to follow) denoted by components
Cartesian Tensors Reference: Jeffreys Cartesian Tensors 1 Coordinates and Vectors z x 3 e 3 y x 2 e 2 e 1 x x 1 Coordinates x i, i 123,, Unit vectors: e i, i 123,, General vector (formal definition to
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 informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationA FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS. 1. Introduction
A FAMILY OF MULTISCALE HYBRID-MIXED FINITE ELEMENT METHODS FOR THE DARCY EQUATION WITH ROUGH COEFFICIENTS CHRISTOPHER HARDER, DIEGO PAREDES 2, AND FRÉDÉRIC VALENTIN Abstract. We aim at proposing novel
More informationOn Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma. Ben Schweizer 1
On Friedrichs inequality, Helmholtz decomposition, vector potentials, and the div-curl lemma Ben Schweizer 1 January 16, 2017 Abstract: We study connections between four different types of results that
More informationA u + b u + cu = f in Ω, (1.1)
A WEIGHTED H(div) LEAST-SQUARES METHOD FOR SECOND-ORDER ELLIPTIC PROBLEMS Z. CAI AND C. R. WESTPHAL Abstract. This paper presents analysis of a weighted-norm least squares finite element method for elliptic
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationOn a variational inequality of Bingham and Navier-Stokes type in three dimension
PDEs for multiphase ADvanced MATerials Palazzone, Cortona (Arezzo), Italy, September 17-21, 2012 On a variational inequality of Bingham and Navier-Stokes type in three dimension Takeshi FUKAO Kyoto University
More informationA NUMERICAL APPROXIMATION OF NONFICKIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXX, 1(21, pp. 75 84 Proceedings of Algoritmy 2 75 A NUMERICAL APPROXIMATION OF NONFICIAN FLOWS WITH MIXING LENGTH GROWTH IN POROUS MEDIA R. E. EWING, Y. LIN and J. WANG
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 information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More information