THE STOKES SYSTEM R.E. SHOWALTER

Size: px
Start display at page:

Download "THE STOKES SYSTEM R.E. SHOWALTER"

Transcription

1 THE STOKES SYSTEM R.E. SHOWALTER Contents 1. Stokes System 1 Stokes System 2 2. The Weak Solution of the Stokes System 3 3. The Strong Solution 4 4. The Normal Trace 6 5. The Mixed Problem 7 6. The Navier-Stokes System 8 1. Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the density and pressure are related by a state equation (1.1a) ρ = s(p) in which the constitutive function s( ) is non-decreasing and characterizes the type of fluid. The conservation of mass of fluid is expressed by (1.1b) t ρ + (ρv) = ρ g(x), and the conservation of momentum has the form (1.1c) t (ρv) + (v )(ρv) (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in. Here f(x) is the mass-distributed force density over and g(x) is the mass-distributed fluid source. These three nonlinear equations comprise the system for a general compressible fluid. For the calculation of the acceleration of the fluid element, the displacement u(x, t) of that element along with the points must be considered. The momentum of the small subdomain B travelling with the fluid is ρv(x + u(x, t), t) dx, and its derivative B is given by the chain rule with v(x, t) = t u(x, t) as ( t ρv(x + u(x, t), t) + j ρv(x + u(x, t), t) v j (x + u(x, t), t)) dx. B Date: May 26, Key words and phrases. basic notions of mechanics, diffusion, fluid flow. 1

2 2 R.E. SHOWALTER Thus, the momentum equation (1.1c) contains the additional term (v )ρv = v j j ρv for convective acceleration. If the fluid is incompressible, the density is constant and the corresponding system is the Navier-Stokes system ρ t v µ 1 v + ρ(v )v + p = ρf + (λ 1 + µ 1 ) g, v = g in, v = 0 on 0, p n + 2µ 1 ε(v, n) = g on 1, for a viscous incompressible fluid. Note that the quadratic nonlinearity arises from the geometry of the motion, and it is not based on any independent assumptions. Stokes System. In order to obtain a linear model to approximate the solutions of the system (1.1), we consider small oscillations about a rest state at which v = 0 and, hence, the functions ρ 0 (x), p 0 (x), f(x), g(x) satisfy ρ 0 = s(p 0 ), g(x) = 0, p 0 = ρ 0 f(x). ρ The quantity c 0 (x) 1 (p ρ 0 (x) p 0(x)) denotes the compressibility of the fluid at the rest state. Using the chain rule with the state equation (1.1a) yields ρ 0 (x) = ρ p p 0(x) = ρ 0 (x)c 0 (x)ρ 0 (x)f(x). Introduce a small parameter ε > 0 to characterize the size of the oscillations and the deviations from the rest state, and consider the corresponding asymptotic expansions Again from the chain rule we obtain ρ = ρ 0 (x) + ερ 1 (x, t) + O(ε 2 ), p = p 0 (x) + εp 1 (x, t) + O(ε 2 ), v = εv 1 (x, t) + O(ε 2 ), g = εg 1 (x) + O(ε 2 ). t ρ 1 (x, t) = ρ 0 (x)c 0 (x) t p 1 (x, t), ρ 1 (x, t) = ρ 0 (x)c 0 (x)p 1 (x, t). Conservation of mass (1.1b) implies to first order in ε that t ρ 1 + (ρ 0 v 1 ) = ρ 0 (x) g 1 (x). From conservation of momentum (1.1c) we get ε t ((ρ 0 +ερ 1 )v 1 )+(εv 1 )ερv 1 ε(λ 1 +µ 1 ) ( v 1 ) εµ 1 v 1 + (p 0 +εp 1 ) = (ρ 0 +ερ 1 )f(x), and then the definitions of rest state and compressibility give the linear system (1.2a) (1.2b) (1.2c) ρ 0 (x) t v 1 (λ 1 + µ 1 ) ( v 1 ) µ 1 v 1 + p 1 = c 0 (x)ρ 0 (x)f(x) p 1, c 0 (x) t p 1 + v 1 + c 0 (x)ρ 0 (x)f(x) v 1 = g 1 (x), v 1 = 0 on 0, λ 1 ( v 1 )n + 2µ 1 ε(v 1, n) p n = 0 on 1 for small variations of a compressible fluid. Here the two sets 0, 1 comprise a partition of the boundary =. In the incompressible case, c 0 (x) = 0, we obtain the Stokes

3 STOKES 3 system, ρ 0 (x) t v 1 µ 1 v 1 + p 1 = (λ 1 + µ 1 ) g 1 (x), v 1 = g 1 (x), v 1 = 0 on 0, 2µ 1 ε(v 1, n) p n = 0 on The Weak Solution of the Stokes System The (slow) flow of an incompressible homogeneous fluid is described by its pressure p(x, t) and velocity v(x, t). The (evolutionary) Stokes system is to find such a pair of functions on the smoothly bounded region in IR n for t > 0 which satisfy the initialboundary-value problem ρ 0 (x) t v µ v + p = f, v = 0 in IR +, v = 0 on IR +, v(0) = v 0 in, where = is the boundary of. If we separate variables, i.e., look for a solution in the form v(x, t) = e λt u(x) for some number λ, we are led to the stationary Stokes system (2.3a) (2.3b) ρ 0 (x)λu µ u + p = f, u = 0 in, u = 0 on, for the pair u(x), p(x). We focus on this system, but the results apply as well to the evolutionary system. Remark 2.1. For a pair of functions u H 1 (), q H 1 (), we have u q dx = u q dx + u n q ds Formally, this shows that Ker( ) = Rg( ) so we expect Rg( ) = Ker( ) up to closure in appropriate spaces. Eventually, we will need to construct carefully realizations of the gradient and divergence operators. Now, if the pair u(x), p(x) is a solution to the stationary Stokes system (2.3), then for every w H 1 () we have λ ρ 0 (x)u w dx + µ u i w i dx ( u i n)w i ds ( w)p dx + (w n)p ds = f w dx. We define the space V 0 = {w H 1 () : w = 0 in, w = 0 on }. The weak form of the problem is now to find (2.4) u V 0 : λ ρ 0 (x)u w dx + µ u i w i dx = f w dx for all w V 0.

4 4 R.E. SHOWALTER Define a continuous bilinear form a(, ) on V 0 by a(u, w) = (λρ 0 (x)u w + µ u i w i ) dx, u, w V 0, and the linear functional f( ) on V 0 by f(w) = f w dx, w V 0, where f L 2 () is given. Note that the principle part of a(, ) is a double sum n u i w i u i w i = x j x j which contains all first-order derivatives. Thus, the space V 0 is appropriate, and we have the estimate n ( ) 2 wi a(w, w) dx. x j This gives the ellipticity condition i,j=1 i,j=1 (2.5) a(w, w) c 0 w 2, for all w V, for some constant c 0 > 0. We noted already that a(, ) is continuous in both variables with respect to the V 0 norm, and likewise the functional f( ) is continuous. We have the following result immediately. Theorem 2.1. There is exactly one weak solution of the stationary Stokes system (2.6) u V 0 : a(u, w) = f(w) for all w V 0. There remains the issue of the sense in which the weak solution of (2.4) satisfies the strong formulation (2.3). 3. The Strong Solution With the space V 0 = {w H 1 () : w = 0 in, w = 0 on }, the weak form of the stationary Stokes problem is given by u V 0 : λ ρ 0 (x)u w dx + µ u i w i dx = f w dx for all w V 0. We need to characterize the annihilator of the space V 0, the kernel of the divergence operator in H 1 0(). This will be used to show the equivalence of the weak form of the Stokes equation with the strong formulation, and it is related to the characterization of the range of the gradient operator in H 1 (). For background, we begin with the following profound result on the annihilator of V {v C 0 () = D() : v = 0}. Theorem 3.1 (DeRham). Let be a domain in IR n. Then f = p in D () if and only if f(v) = 0 for all v V. A related result is the following.

5 STOKES 5 Theorem 3.2. Let be a bounded domain in IR n with Lipschitz boundary. If p D () satisfies p H 1 (), then p L 2 (), and we have the estimate (3.7) p L 2 ()/IR C p H 1 (). Theorem 3.2 was proved by Magenes & Stampacchia if has a C 1 boundary and by Neĉas for the case of a Lipschitz boundary. Related results were obtained by Deny & Lions. Recall that the space L 2 ()/IR is just the quotient of L 2 () with the constant functions. An estimate equivalent to (3.7) is ) (3.8) p L 2 () C ( p(x) dx + p H 1 (). Furthermore, these estimates are equivalent to the statement that the gradient operator : L 2 () H 1 () has closed range. Now we define the space V 1 to be the closure in H 1 0() of the space V. Theorem 3.3. Let be a bounded domain in IR n with Lipschitz boundary. Then V 1 = V 0. Proof. First we check that V 1 V 0 : if u = lim u m in H 1 0(), with u m V, then we have u m u = 0, so u V 0. Next we show that V 1 = V 0. Let T V 0 and assume further that the restriction to V 1 vanishes. Since V 0 is a closed subspace of H 1 0() = H0() 1 n, T has a continuous linear extension to all of H 1 0(), namely, t H 1 (). But t(v) = 0 for all v V, so from Theorem 3.1 and Theorem 3.2 we obtain t = p for some p L 2 (). Therefore, T (v) = t(v) = (p, v) = 0 for all v V 0, and this shows that T = 0. Since this is true for every T V 0 for which T V1 = 0, we have V 1 = V 0 as desired. Corollary 3.4. The annihilator V a 0 is { p : p L 2 ()} in H 1 (), and V is dense in V 0. Corollary 3.5. If u is a solution of the Stokes problem, then u H 1 () and there is a p L 2 () for which ρ 0 (x)λu µ u + p = f, u = 0 in, u = 0 on. Finally, we show that from Theorem 3.2 we can obtain a restricted form of Theorem 3.1 that is sufficient for our purposes above. Thus, our Theorem 3.3 is independent of Theorem 3.1. Theorem 3.6 (Luc Tartar). Let be a bounded domain in IR n with Lipschitz boundary. Then f H 1 () and f(v) = 0 for all v V if and only if f = p, p L 2 (). Proof. Let there be an increasing sequence of open sets for which m m+1 and m =. For each m 1 the gradient operator m L(L 2 (), H 1 ()) has closed range. Thus, its adjoint m L(H 1 0(), L 2 ()) is divergence, and its annihilator in H 1 () is Ker( m) a = Rg( m ).

6 6 R.E. SHOWALTER Let f H 1 () and f(v) = 0 for all v V. Then for each u Ker( m), the zeroextension to is denoted by ũ, and it has compact support in. By a regularization argument, ũ is the limit in H 1 0() of functions from V, so we have ũ V 0 and f(ũ) = 0. Hence, the restriction of f to m is in Ker( m) a = Rg( m ), so f = p m on m for a p m L 2 ( m ). Since m m+1 we have p m+1 p m = c on m, and we can modify p m+1 by a constant so that p m+1 = p m on m. Thus, f = p with p L 2 loc (), and then by Theorem 3.2 we obtain p L 2 (). 4. The Normal Trace Recall that the trace map γ : H 1 () L 2 () is the restriction to the boundary of each function u in H 1 (). The kernel of this continuous and linear map is the subspace H 1 0() of those functions which vanish on the boundary, and it is characterized as the closure of C 0 () in H 1 (). We denote the range of γ by B, and it follows that this is a complete space with the quotient norm ψ B inf{ w H 1 : w H 1 () with γ(w) = ψ}. Note that B is the space of all boundary values of functions in H 1 (), and we can identify L 2 () with its dual and obtain B L 2 () B. The corresponding quotient map of H 1 ()/H0() 1 onto B is an isomorphism, and so also is the dual map γ of B onto H0() 1 a, the indicated annihilator of H0() 1 in the dual H 1 (). Now we define the space W() = {u L 2 () : u L 2 ()} with the graph norm u 2 W() = u 2 L 2 + u 2 L 2 It follows easily that this space is complete, and it can be shown that the smooth functions are dense in this space. Furthermore W() is much larger than the space H 1 () = H 1 () n of vector-valued functions, and the injection H 1 () W() is continuous. The trace map has a natural meaning on H 1 (), namely, componentwise, but we shall need to extend it to all of W(). This is the point of the following result. Theorem 4.1. There exists a unique map γ n continuous and linear from W() onto B such that, for all functions u H 1 (), γ n (u) = γ(u) n is the normal trace, i.e., the generalized restriction to the boundary of u n, and this map satisfies the extended Stokes formula (4.9) (u, w) L 2 + ( u, w) L 2 = γ n (u)(γ(w)), u W(), w H 1 (). Proof. Let u W() and define f H 1 () by f(w) = ( uw + u w) dx, w H 1 (). For each w H0() 1 we have f(w) = 0. (This follows first for w C0 and then by continuity for all w H0().) 1 That is, f H0() 1 a, so from above we see there is a g B for which f = γ (g) = g γ. It follows that for each ψ B and for each w H 1 () with γ(w) = ψ, the value of (4.10) g(ψ) = ( uw + u w) dx

7 STOKES 7 is independent of the choice of w, and we have g(ψ) u W() w H 1 (). By taking the infimum over all such w we obtain g(ψ) u W() ψ B. It follows that this continuous linear functional on B satisfies g B u W(). We denote the dependence of this g on u by setting g = γ n (u), and this defines γ n L(W(), B ) as desired. Finally, we note that for u H 1 () and smooth w on we have γ n (u)(γ(w)) = (uw) dx = γ(u) nγ(w) ds, and the smooth functions are dense in H 1 (), so we have γ n (u) = u n in L 2 (). It remains to show that γ n maps onto B. Let h B, and then solve the Neumann problem for a p H 1 () with p p = 0 in and p = h on. That is, n p H 1 () : (p q + p q) dx = h(γ(q)) for all q H 1 (). This has exactly one solution, and we then set u = p L 2 (). It follows that u = p L 2 (), and so u W() with h = γ n (u) on. Remark 4.1. By solving an appropriate Dirichlet problem, we can show that the functional (4.10) satisfies g B = u W(). Remark 4.2. Note that it is only the normal component of the boundary trace that retains any meaning for functions from W(). 5. The Mixed Problem We shall consider the stationary Stokes system λu µ u + p = f, u = 0 in, u = 0 on 0, µ u n p n = g on 1, for the pair u(x), p(x). Here the two sets 0, 1 comprise a partition of the boundary =, and we are given the functions f L 2 () and g L 2 ( 1 ). Now, if the pair u(x), p(x) is a solution to the stationary Stokes system, then for every w H 1 () we have λu w dx + µ p ( w) dx + u i w i dx µ ( u i n)w i ds (w n)p ds = f w dx. We define the space V = {w H 1 () : w = 0 in, w = 0 on 0 }. Then our solution is chosen from this space, and if we also choose our test function w V above, it follows that λu w dx + µ u i w i dx µ ( u i n)w i ds 1 + (w n)p ds = f w dx. 1

8 8 R.E. SHOWALTER Thus, the weak form of the problem is now to find u V : (λu w + µ u i w i ) dx = f w dx + g w ds for all w V. 1 As before, we obtain the following result. Theorem 5.1. There is exactly one weak solution u to the stationary Stokes system. Note also that the pressure p is determined up to a constant. Conversely, let s assume that u is a solution of the weak form of the stationary Stokes system. Then from the inclusion u V we obtain u H 1 (), u = 0 in, and u = 0 on 0 in the sense of boundary trace. Furthermore, by taking test functions w V 0 = V H 1 0(), we find that λu µ u f V a 0, the indicated annihilator of V 0 in H 1 (), and so it is a gradient, and we have λu µ u f = p, p L 2 (). It follows that for each i = 1, 2,..., n we have µ u i p e i L 2 () and (µ u i p e i ) is the i-th component of µ u p = λu f L 2 (). This shows that we have each µ u i p e i W() and so there is a well-defined normal trace on the boundary. We use this equation to substitute for f in the weak form to obtain (µ u i w i + (µ u p) w) dx = g w ds for all w V. 1 Since w = 0, the i-th component of the integrand on the left side is given by [(µ u i pe i )w i ]. The generalized Stokes theorem implies that its integral over is given by γ n (µ u i pe i ) w i ds, 1 and so from above we obtain 1 ( µ u n p n It is in this sense that we have ) w dx = g w ds for all w V. 1 µ u n p n = g on 1. Remark 5.1. Because of the special properties of the space B, the condition γ n (w) = 0 on 1 is somewhat different from the condition γ n (w) = 0 for all w V. 6. The Navier-Stokes System We shall consider the stationary Navier-Stokes system (6.11) λu µ u + (u )u + p = f, u = 0 in, u = 0 on, for the pair u H 1 (), p L 2 (). The boundary condition is meant in the sense of trace on =, and we are given the function f L 2 (). This is distinguished from the Stokes system by the nonlinear momentum term (u )u = u j j u, and the corresponding problem is the Navier-Stokes system.

9 STOKES 9 Now, if the pair u(x), p(x) is a solution to the Navier-Stokes system, then for every w H 1 () we have λu w dx + µ u i w i dx µ ( u i n)w i ds + (u )u w dx p ( w) dx + (w n)p ds = f w dx. We define the space V 0 = {w H 1 () : w = 0 in, w = 0 on }. Then our solution belongs to this space, and if we also choose our test function w V 0 above, it follows that λu w dx + µ u i w i dx + (u )u w dx = f w dx. Thus, the weak form of the problem is now to find (6.12) u V 0 : (λu w + µ u i w i ) dx + (u )u w dx = f w dx for all w V 0. Next we confirm that the additional momentum terms above are meaningful, so we define the function (6.13) b(u, v, w) = (u )v w dx for appropriate u, v, w. To this end, we obtain the following consequences of the Sobolev imbedding Theorem 4.3 from our textbook. For each u H 1 0, we have These lead to the following. u L q C() u H 1, 1 q <, if n = 2, u L 6 C() u H 1, if n = 3, u L 4 C() u H 1, if n = 4. Lemma 6.1. If 2 n 4, then the trilinear form b(u, v, w) is bounded on H 1 0() 3. Proof. We use Hölder s inequality repeatedly. For n = 2: u j j v i w i dx C u j L 4 j v i L 2 w i L 4. For n 3: u j j v i w i dx C u j L 2n/(n 2) j v i L 2 w i L n, since n 2 2n + 1 n = 1 2.

10 10 R.E. SHOWALTER Conversely, let s assume that u is a solution of the weak Navier-Stokes system (6.12). Then from the inclusion u V 0 we obtain u H 1 (), u = 0 in, and u = 0 on in the sense of boundary trace. Furthermore, by taking test functions w V 0, we find that λu µ u + (u )u f V a 0, the indicated annihilator of V 0 in H 1 (), and so it is a gradient, and we have λu µ u + (u )u f = p, p L 2 (). Thus, u is a solution of the original system (6.11), and so these two formulations are equivalent. We develop some additional properties of the trilinear form b(u, v, w). Lemma 6.2. For u V 0 and v H 1 0() we have b(u, v, v) = 0. Proof. u j ( j v i ) v i dx = 1 u j j (v i ) 2 dx = 1 ( j u j )(v i ) 2 dx 2 2 so it follows by summing on all these that b(u, v, v) = 1 ( u) v 2 IR n dx = 0. Corollary 6.3. For u V 0 and v, w H 1 0() we have b(u, v, w) = b(u, w, v). Proof. 0 = b(u, v + w, v + w) = b(u, v, w) + b(u, w, v) We shall prove there is a weak solution of the Navier-Stokes system (6.12) as a direct application of our basic existence Theorem 2.1 of our textbook. To do so, we begin by defining an operator A : V 0 V 0 by A(u)(w) = λu w dx + µ u i w i dx + b(u, u, w), u, w V 0. We have shown in Lemma 6.1 that A( ) is bounded, and we check from above that it is V 0 -coercive: by Corollary 6.3 we have A(u)(u) µ u i u i dx + b(u, u, u) c 0 u 2, w V 0, for some constant c 0 > 0. It remains to establish some form of continuity of A( ). Lemma 6.4. The map u b(u, u, ) is weakly continuous from V 0 to V 0, that is, if u m w u in V 0 then b(u m, u m, w) b(u, u, w) for each w V 0. Proof. By the Rellich-Kandorochov Theorem 4.2 of the text, we have strong convergence u m u in L 2 (), and so for each w V we obtain b(u m, u m, w) = b(u m, w, u m ) = u m,j ( j w i ) u m,i dx b(u, w, u) = b(u, u, w). Since the sequence {b(u m, u m, )} is bounded in V 0, the same holds for each w V 0.

11 STOKES 11 Note that the preceding result depended on the observation that for each w V the bilinear function b(, w, ) is continuous on L 2 (). Corollary 6.5. The operator A : V 0 V 0 is weakly continuous. Theorem 6.6. There exists a weak solution u to the stationary Navier Stokes system (6.12). Proof. We have shown that A : V 0 V 0 is weakly continuous, hence, of type M, bounded, and coercive, so it is onto. Finally we note that the operator A is not monotone. We actually used the more general form of our Theorem 2.1. Furthermore, we have not said anything about uniqueness of the solution. Department of Mathematics, Oregon State University, Corvallis, OR address: show@math.oregonstate.edu

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.

1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω. 1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the

More information

Weak Formulation of Elliptic BVP s

Weak Formulation of Elliptic BVP s Weak Formulation of Elliptic BVP s There are a large number of problems of physical interest that can be formulated in the abstract setting in which the Lax-Milgram lemma is applied to an equation expressed

More information

LECTURE # 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) 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 information

Mixed exterior Laplace s problem

Mixed exterior Laplace s problem Mixed exterior Laplace s problem Chérif Amrouche, Florian Bonzom Laboratoire de mathématiques appliquées, CNRS UMR 5142, Université de Pau et des Pays de l Adour, IPRA, Avenue de l Université, 64000 Pau

More information

SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS

SPECTRAL 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 information

1 Constrained Optimization

1 Constrained Optimization 1 Constrained Optimization Let B be an M N matrix, a linear operator from the space IR N to IR M with adjoint B T. For (column vectors) x IR N and y IR M we have x B T y = Bx y. This vanishes for all y

More information

SHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction

SHARP 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 information

Numerical Methods for the Navier-Stokes equations

Numerical 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 information

Chapter 2 Finite Element Spaces for Linear Saddle Point Problems

Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark

More information

Second Order Elliptic PDE

Second Order Elliptic PDE Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic

More information

Solution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0

Solution Sheet 3. Solution Consider. with the metric. We also define a subset. and thus for any x, y X 0 Solution Sheet Throughout this sheet denotes a domain of R n with sufficiently smooth boundary. 1. Let 1 p

More information

1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e.,

1. Subspaces A subset M of Hilbert space H is a subspace of it is closed under the operation of forming linear combinations;i.e., Abstract Hilbert Space Results We have learned a little about the Hilbert spaces L U and and we have at least defined H 1 U and the scale of Hilbert spaces H p U. Now we are going to develop additional

More information

LECTURE 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, 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 information

The 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 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 information

Variational Formulations

Variational 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 information

Appendix A Functional Analysis

Appendix 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 information

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis

More information

Existence of Weak Solutions to a Class of Non-Newtonian Flows

Existence of Weak Solutions to a Class of Non-Newtonian Flows Existence of Weak Solutions to a Class of Non-Newtonian Flows 1. Introduction and statement of the result. Ladyzhenskaya [8]. Newtonian : Air, other gases, water, motor oil, alcohols, simple hydrocarbon

More information

Theory of PDE Homework 2

Theory of PDE Homework 2 Theory of PDE Homework 2 Adrienne Sands April 18, 2017 In the following exercises we assume the coefficients of the various PDE are smooth and satisfy the uniform ellipticity condition. R n is always an

More information

CONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence

CONVERGENCE 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 information

On a Nonlocal Elliptic System of p-kirchhoff-type Under Neumann Boundary Condition

On a Nonlocal Elliptic System of p-kirchhoff-type Under Neumann Boundary Condition On a Nonlocal Elliptic System of p-kirchhoff-type Under Neumann Boundary Condition Francisco Julio S.A Corrêa,, UFCG - Unidade Acadêmica de Matemática e Estatística, 58.9-97 - Campina Grande - PB - Brazil

More information

Lectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai.

Lectures on. Sobolev Spaces. S. Kesavan The Institute of Mathematical Sciences, Chennai. Lectures on Sobolev Spaces S. Kesavan The Institute of Mathematical Sciences, Chennai. e-mail: kesh@imsc.res.in 2 1 Distributions In this section we will, very briefly, recall concepts from the theory

More information

TRANSPORT IN POROUS MEDIA

TRANSPORT IN POROUS MEDIA 1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem

More information

Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.

Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological

More information

Lecture Notes on PDEs

Lecture Notes on PDEs Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential

More information

MATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY

MATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY MATH 263: PROBLEM SET 1: BUNDLES, SHEAVES AND HODGE THEORY 0.1. Vector Bundles and Connection 1-forms. Let E X be a complex vector bundle of rank r over a smooth manifold. Recall the following abstract

More information

The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control

The 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 information

Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems

Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter

More information

EXAMINATION SOLUTIONS

EXAMINATION SOLUTIONS 1 Marks & (a) If f is continuous on (, ), then xf is continuously differentiable 5 on (,) and F = (xf)/x as well (as the quotient of two continuously differentiable functions with the denominator never

More information

Fact Sheet Functional Analysis

Fact Sheet Functional Analysis Fact Sheet Functional Analysis Literature: Hackbusch, W.: Theorie und Numerik elliptischer Differentialgleichungen. Teubner, 986. Knabner, P., Angermann, L.: Numerik partieller Differentialgleichungen.

More information

Weak Convergence Methods for Energy Minimization

Weak Convergence Methods for Energy Minimization Weak Convergence Methods for Energy Minimization Bo Li Department of Mathematics University of California, San Diego E-mail: bli@math.ucsd.edu June 3, 2007 Introduction This compact set of notes present

More information

MAXIMIZATION AND MINIMIZATION PROBLEMS RELATED TO A p-laplacian EQUATION ON A MULTIPLY CONNECTED DOMAIN. N. Amiri and M.

MAXIMIZATION AND MINIMIZATION PROBLEMS RELATED TO A p-laplacian EQUATION ON A MULTIPLY CONNECTED DOMAIN. N. Amiri and M. TAIWANESE JOURNAL OF MATHEMATICS Vol. 19, No. 1, pp. 243-252, February 2015 DOI: 10.11650/tjm.19.2015.3873 This paper is available online at http://journal.taiwanmathsoc.org.tw MAXIMIZATION AND MINIMIZATION

More information

MAT 578 FUNCTIONAL ANALYSIS EXERCISES

MAT 578 FUNCTIONAL ANALYSIS EXERCISES MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.

More information

Homogenization of the compressible Navier Stokes equations in domains with very tiny holes

Homogenization of the compressible Navier Stokes equations in domains with very tiny holes Homogenization of the compressible Navier Stokes equations in domains with very tiny holes Yong Lu Sebastian Schwarzacher Abstract We consider the homogenization problem of the compressible Navier Stokes

More information

Simple Examples on Rectangular Domains

Simple Examples on Rectangular Domains 84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation

More information

Chapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation

Chapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f

More information

Equations paraboliques: comportement qualitatif

Equations paraboliques: comportement qualitatif Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Equations paraboliques: comportement qualitatif par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 25/6 1 Contents

More information

Analysis in weighted spaces : preliminary version

Analysis in weighted spaces : preliminary version Analysis in weighted spaces : preliminary version Frank Pacard To cite this version: Frank Pacard. Analysis in weighted spaces : preliminary version. 3rd cycle. Téhéran (Iran, 2006, pp.75.

More information

THE L 2 -HODGE THEORY AND REPRESENTATION ON R n

THE L 2 -HODGE THEORY AND REPRESENTATION ON R n THE L 2 -HODGE THEORY AND REPRESENTATION ON R n BAISHENG YAN Abstract. We present an elementary L 2 -Hodge theory on whole R n based on the minimization principle of the calculus of variations and some

More information

ASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS

ASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS ASYMPTOTIC BEHAVIOUR OF NONLINEAR ELLIPTIC SYSTEMS ON VARYING DOMAINS Juan CASADO DIAZ ( 1 ) Adriana GARRONI ( 2 ) Abstract We consider a monotone operator of the form Au = div(a(x, Du)), with R N and

More information

Domain Perturbation for Linear and Semi-Linear Boundary Value Problems

Domain Perturbation for Linear and Semi-Linear Boundary Value Problems CHAPTER 1 Domain Perturbation for Linear and Semi-Linear Boundary Value Problems Daniel Daners School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia E-mail: D.Daners@maths.usyd.edu.au

More information

Euler Equations: local existence

Euler Equations: local existence Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u

More information

ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen

ON 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 information

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang

ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of

More information

Some Notes on Elliptic Regularity

Some Notes on Elliptic Regularity Some Notes on Elliptic Regularity Here we consider first the existence of weak solutions to elliptic problems of the form: { Lu = f, u = 0, (1) and then we consider the regularity of such solutions. The

More information

Regularity Theory a Fourth Order PDE with Delta Right Hand Side

Regularity Theory a Fourth Order PDE with Delta Right Hand Side Regularity Theory a Fourth Order PDE with Delta Right Hand Side Graham Hobbs Applied PDEs Seminar, 29th October 2013 Contents Problem and Weak Formulation Example - The Biharmonic Problem Regularity Theory

More information

Uniform 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 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 information

A metric space X is a non-empty set endowed with a metric ρ (x, y):

A metric space X is a non-empty set endowed with a metric ρ (x, y): Chapter 1 Preliminaries References: Troianiello, G.M., 1987, Elliptic differential equations and obstacle problems, Plenum Press, New York. Friedman, A., 1982, Variational principles and free-boundary

More information

Variable Exponents Spaces and Their Applications to Fluid Dynamics

Variable 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 information

The Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation:

The Dirichlet s P rinciple. In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: Oct. 1 The Dirichlet s P rinciple In this lecture we discuss an alternative formulation of the Dirichlet problem for the Laplace equation: 1. Dirichlet s Principle. u = in, u = g on. ( 1 ) If we multiply

More information

Uniform 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 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 information

Sobolevology. 1. Definitions and Notation. When α = 1 this seminorm is the same as the Lipschitz constant of the function f. 2.

Sobolevology. 1. Definitions and Notation. When α = 1 this seminorm is the same as the Lipschitz constant of the function f. 2. Sobolevology 1. Definitions and Notation 1.1. The domain. is an open subset of R n. 1.2. Hölder seminorm. For α (, 1] the Hölder seminorm of exponent α of a function is given by f(x) f(y) [f] α = sup x

More information

Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...

Functional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability... Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................

More information

Sobolev Spaces. Chapter 10

Sobolev Spaces. Chapter 10 Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p

More information

NECESSARY OPTIMALITY CONDITIONS FOR AN OPTIMAL CONTROL PROBLEM OF 2D

NECESSARY OPTIMALITY CONDITIONS FOR AN OPTIMAL CONTROL PROBLEM OF 2D NECESSARY OPTIMALITY CONDITIONS FOR AN OPTIMAL CONTROL PROBLEM OF 2D g-navier-stokes EQUATIONS Nguyen Duc Loc 1 Abstract Considered here is the optimal control problem of 2D g-navier-stokes equations.

More information

ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM

ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,

More information

Math The Laplacian. 1 Green s Identities, Fundamental Solution

Math The Laplacian. 1 Green s Identities, Fundamental Solution Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external

More information

INCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS (II) Luan Thach Hoang. IMA Preprint Series #2406.

INCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS (II) Luan Thach Hoang. IMA Preprint Series #2406. INCOMPRESSIBLE FLUIDS IN THIN DOMAINS WITH NAVIER FRICTION BOUNDARY CONDITIONS II By Luan Thach Hoang IMA Preprint Series #2406 August 2012 INSTITUTE FOR MATHEMATICS AND ITS APPLICATIONS UNIVERSITY OF

More information

Variational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation

Variational Theory of Solitons for a Higher Order Generalized Camassa-Holm Equation International Journal of Mathematical Analysis Vol. 11, 2017, no. 21, 1007-1018 HIKAI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.710141 Variational Theory of Solitons for a Higher Order Generalized

More information

Solutions of Selected Problems

Solutions of Selected Problems 1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative

More information

i=1 α i. Given an m-times continuously

i=1 α i. Given an m-times continuously 1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable

More information

From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray

From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park

More information

INF-SUP CONDITION FOR OPERATOR EQUATIONS

INF-SUP CONDITION FOR OPERATOR EQUATIONS INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent

More information

On 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 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 information

Mathematical analysis of the stationary Navier-Stokes equations

Mathematical 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 information

Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms

Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz

More information

Finite Elements. Colin Cotter. February 22, Colin Cotter FEM

Finite 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 information

PROPERTIES OF CAPACITIES IN VARIABLE EXPONENT SOBOLEV SPACES

PROPERTIES OF CAPACITIES IN VARIABLE EXPONENT SOBOLEV SPACES PROPERTIES OF CAPACITIES IN VARIABLE EXPONENT SOBOLEV SPACES PETTERI HARJULEHTO, PETER HÄSTÖ, AND MIKA KOSKENOJA Abstract. In this paper we introduce two new capacities in the variable exponent setting:

More information

EXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON

EXISTENCE 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 information

EXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS

EXISTENCE OF SOLUTIONS FOR A RESONANT PROBLEM UNDER LANDESMAN-LAZER CONDITIONS Electronic Journal of Differential Equations, Vol. 2008(2008), No. 98, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) EXISTENCE

More information

Laplace s Equation. Chapter Mean Value Formulas

Laplace 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 information

HI CAMBRIDGE n S P UNIVERSITY PRESS

HI CAMBRIDGE n S P UNIVERSITY PRESS Infinite-Dimensional Dynamical Systems An Introduction to Dissipative Parabolic PDEs and the Theory of Global Attractors JAMES C. ROBINSON University of Warwick HI CAMBRIDGE n S P UNIVERSITY PRESS Preface

More information

A PROPERTY OF SOBOLEV SPACES ON COMPLETE RIEMANNIAN MANIFOLDS

A PROPERTY OF SOBOLEV SPACES ON COMPLETE RIEMANNIAN MANIFOLDS Electronic Journal of Differential Equations, Vol. 2005(2005), No.??, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A PROPERTY

More information

The oblique derivative problem for general elliptic systems in Lipschitz domains

The 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 information

STOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN

STOKES PROBLEM WITH SEVERAL TYPES OF BOUNDARY CONDITIONS IN AN EXTERIOR DOMAIN Electronic Journal of Differential Equations, Vol. 2013 2013, No. 196, pp. 1 28. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STOKES PROBLEM

More information

COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS

COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS Dynamic Systems and Applications 22 (203) 37-384 COMBINED EFFECTS FOR A STATIONARY PROBLEM WITH INDEFINITE NONLINEARITIES AND LACK OF COMPACTNESS VICENŢIU D. RĂDULESCU Simion Stoilow Mathematics Institute

More information

The Helmholtz Equation

The Helmholtz Equation The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness

More information

IMA Preprint Series # 2143

IMA 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 information

THE HODGE DECOMPOSITION

THE HODGE DECOMPOSITION THE HODGE DECOMPOSITION KELLER VANDEBOGERT 1. The Musical Isomorphisms and induced metrics Given a smooth Riemannian manifold (X, g), T X will denote the tangent bundle; T X the cotangent bundle. The additional

More information

Notes on Sobolev Spaces A. Visintin a.a

Notes on Sobolev Spaces A. Visintin a.a Notes on Sobolev Spaces A. Visintin a.a. 2017-18 Contents: 1. Hölder spaces. 2. Regularity of Euclidean domains. 3. Sobolev spaces of positive integer order. 4. Sobolev spaces of real integer order. 5.

More information

TD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle

TD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class

More information

Non-stationary Friedrichs systems

Non-stationary Friedrichs systems Department of Mathematics, University of Osijek BCAM, Bilbao, November 2013 Joint work with Marko Erceg 1 Stationary Friedrichs systems Classical theory Abstract theory 2 3 Motivation Stationary Friedrichs

More information

RNDr. Petr Tomiczek CSc.

RNDr. Petr Tomiczek CSc. HABILITAČNÍ PRÁCE RNDr. Petr Tomiczek CSc. Plzeň 26 Nonlinear differential equation of second order RNDr. Petr Tomiczek CSc. Department of Mathematics, University of West Bohemia 5 Contents 1 Introduction

More information

On the p-laplacian and p-fluids

On the p-laplacian and p-fluids LMU Munich, Germany Lars Diening On the p-laplacian and p-fluids Lars Diening On the p-laplacian and p-fluids 1/50 p-laplacian Part I p-laplace and basic properties Lars Diening On the p-laplacian and

More information

Lectures on. Weak Solutions of Elliptic Boundary Value Problems

Lectures on. Weak Solutions of Elliptic Boundary Value Problems Lectures on Weak Solutions of Elliptic Boundary Value Problems S. Kesavan The Institute of Mathematical Sciences, CIT Campus, Taramani, Chennai - 600 113. e mail: kesh@imsc.res.in 2 1 Some abstract variational

More information

Navier-Stokes equations in thin domains with Navier friction boundary conditions

Navier-Stokes equations in thin domains with Navier friction boundary conditions Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu

More information

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT

SEMILINEAR ELLIPTIC EQUATIONS WITH DEPENDENCE ON THE GRADIENT Electronic Journal of Differential Equations, Vol. 2012 (2012), No. 139, pp. 1 9. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu SEMILINEAR ELLIPTIC

More information

9. Boundary value problems in a constant-coefficient case

9. Boundary value problems in a constant-coefficient case 9.1 9. Boundary value problems in a constant-coefficient case 9.1. Boundary maps for the half-space. We have considered various realizations in Section 4.4 of the Laplacian and similar operators, defined

More information

Chapter 2. General concepts. 2.1 The Navier-Stokes equations

Chapter 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 information

Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains

Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains Equilibria with a nontrivial nodal set and the dynamics of parabolic equations on symmetric domains J. Földes Department of Mathematics, Univerité Libre de Bruxelles 1050 Brussels, Belgium P. Poláčik School

More information

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

More information

Getting started: CFD notation

Getting started: CFD notation PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =

More information

New estimates for the div-curl-grad operators and elliptic problems with L1-data in the whole space and in the half-space

New estimates for the div-curl-grad operators and elliptic problems with L1-data in the whole space and in the half-space New estimates for the div-curl-grad operators and elliptic problems with L1-data in the whole space and in the half-space Chérif Amrouche, Huy Hoang Nguyen To cite this version: Chérif Amrouche, Huy Hoang

More information

A Computational Approach to Study a Logistic Equation

A Computational Approach to Study a Logistic Equation Communications in MathematicalAnalysis Volume 1, Number 2, pp. 75 84, 2006 ISSN 0973-3841 2006 Research India Publications A Computational Approach to Study a Logistic Equation G. A. Afrouzi and S. Khademloo

More information

Konstantinos Chrysafinos 1 and L. Steven Hou Introduction

Konstantinos Chrysafinos 1 and L. Steven Hou Introduction Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ANALYSIS AND APPROXIMATIONS OF THE EVOLUTIONARY STOKES EQUATIONS WITH INHOMOGENEOUS

More information

On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions

On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability Boundary Conditions Proceedings of the 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August -, 5 pp36-41 On a Suitable Weak Solution of the Navier Stokes Equation with the Generalized Impermeability

More information

Partial Differential Equations

Partial Differential Equations Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,

More information

Basic Principles of Weak Galerkin Finite Element Methods for PDEs

Basic 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 information

GRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS

GRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS LE MATEMATICHE Vol. LI (1996) Fasc. II, pp. 335347 GRAND SOBOLEV SPACES AND THEIR APPLICATIONS TO VARIATIONAL PROBLEMS CARLO SBORDONE Dedicated to Professor Francesco Guglielmino on his 7th birthday W

More information