Partition of Unity Methods for Heterogeneous Domain Decomposition
|
|
- Anna Phelps
- 5 years ago
- Views:
Transcription
1 Partition of Unity Methods for Heterogeneous Domain Decomposition Gabriele Ciaramella and Martin J. Gander 1 Heterogeneous problems and partition of unity decomposition We are interested in solving linear PDEs of the form L (u) = f in, u = g on, (1) where is a bounded domain in R d with d = 1,2, L is a linear (elliptic) differential operator, f and g are the data, and u is the solution to (1). The weak form of (1) with a Hilbert space (V,, ) of functions v : R is a(u,v) = l(v) v V 0, with u = g on, (2) where V 0 := {v V : v = 0 on }, a : V V R is the bilinear form corresponding to the operator L, and l : V R is the linear functional induced by f. We assume that (2) has a unique solution u {v V : v = g on }, and that u is heterogeneous, behaving very differently in different parts of. Typical examples are advection-diffusion problems, where there are advection dominated and diffusion dominated regions (subdomains), and the boundaries in between are not clearly defined, see [8, 10] and references therein. Apart from the χ-method [6, 1], there are no methods to determine such subdomain decompositions, and our goal is to present and study a new such method. We thus introduce (see [18, 11]) Definition 1 (Membership function). Let R d be a set. A membership function ϕ is a map ϕ : [0,1], and its support S is S := {x : ϕ(x) 0}. Given two membership functions ϕ 1,ϕ 2 : [0,1] that form a partition of unity on, ϕ 1 (x) + ϕ 2 (x) = 1 for all x, their supports provide then a domain decompo- G. Ciaramella Universität Konstanz, Germany, gabriele.ciaramella@uni-konstanz.de M. J. Gander Université de Genève, Genève, Switzerland martin.gander@unige.ch 1
2 2 Gabriele Ciaramella and Martin J. Gander sition = suppϕ 1 suppϕ 2. We introduce the approximation u dd := ϕ 1 u 1 +ϕ 1 u 2 u, where u 1 and u 2 represent two different possible behaviors of u, and we assume that u dd = g on. We proceed as follows to define the spaces that u dd, u 1 and u 2 are to be sought in: first, we introduce two approximate problems, L 1 (u 1 ) = f 1 in, B 1 (u 1, g) = g, and L 2 (u 2 ) = f 2 in, B 2 (u 2, g) = 0. (3) Here L j are approximation operators of L, f j are approximations of f, and B j are operators to define the boundary conditions of (3), see Section 3 for concrete examples. The function g represents a control and belongs to an appropriate Hilbert space W. Notice that g is different from the actual boundary data g: the latter is defined on, while we will define the former only on a subset of. We assume 1 that (3) (left) is uniquely solvable in V for any g W and (3) (right) has a unique solution u 2 V. To reformulate (3) (left), we introduce two operators A : V V and B g : W V, such that (3) (left) becomes Au 1 = B g g+ f 1. Notice that B g represents the boundary conditions of (3) (left) and takes into account also g. This problem is formally solved by u 1 = A 1 B g g + A 1 f 1, where A 1 is well defined if (3) (left) is well posed. Now, we define the spaces V 1 := {v V : v = A 1 (B g q + f 1 ), q W}, and V 2 := {u 2 }. Here V 1 represents the space of all possible solutions to the first problem in (3) generated by all the possible (control) functions in W, while V 2 is a singleton containing only the unique solution u 2 to (3) (right). Finally, we use the definition of a partition of unity method space (PUM-space [2, 13]) V PUM := ϕ 1 V 1 +ϕ 2 V 2 V, where ϕ 1,ϕ 2 are membership functions. V PUM, V 1 and V 2 are the spaces that the approximations u dd, u 1 and u 2 have to be sought in. In particular, for the approximation u dd the functions ϕ 1, ϕ 2 and g have to be computed. These are defined as solutions to optimal control problems, as described in Section 2. Here we need to remark that our approach could be computationally expensive. However, it is motivated by applications in astrophysics governed by hyperbolic equations like the Boltzmann equation. In many cases, like for supernova explosion, physical phenomena are modeled using two different (limiting) regimes. However, this would require an a-priori knowledge of the transition regime; see, e.g. [8, 3, 11] and references therein. This is exactly the role of the partition of unity functions obtained by our computational framework. In practice, one could use our computationally expensive approach to obtain the partition of unity functions for one representative case and then reuse them (as approximations) in a domain decomposition fashion to compute approximate solutions of other cases of interest. 1 This specific approximation is motivated by asymptotic expansion techniques providing in general two problems, one that is uniquely determined and a second one that is determined up to some constants for asymptotic matching [15].
3 Partition of Unity Methods for Heterogeneous Domain Decomposition 3 2 Optimal control approaches To compute ϕ 1, ϕ 2 and g, we embed the PUM formulation into an optimal control framework. We begin by inserting u dd into (2) and obtain the bounded linear functional r : V R defined by r(v) := a(ϕ 1 u 1 + ϕ 2 u 2,v) l(v), where v V. In the case that v = ϕ 1 w and v = ϕ 2 w with w V, we get the functionals r j (w) := r(ϕ j w) = a(ϕ 1 u 1 + ϕ 2 u 2,ϕ j w) l(ϕ j w) w V, for j = 1,2. Since w V r j (w), j = 1,2, are bounded linear functionals, they are elements in V, and by the Riesz representation theorem [7], there exist R 1 and R 2 in V such that R j,v = a(ϕ 1 u 1 + ϕ 2 u 2,ϕ j v) l(ϕ j v) v V 0, j = 1,2, (4) where we used V 0, since u dd is exact on and thus R 1 and R 2 must vanish there. Now, we define ϕ := ϕ 1 with ϕ 2 = 1 ϕ, and recall that r j V = R j V. Minimizing the norms of the residuals R j V leads to the optimal control problem min R 1,R 2,u 1,g,ϕ J(R 1,R 2,g,ϕ) := 1 2 R 1 2 V R 2 2 V + α 2 ϕ 2 V + β 2 g 2 W s.t. R 1,v = a(ϕu 1 + (1 ϕ)u 2,ϕv) l(ϕv) v V 0, R 2,v = a(ϕu 1 + (1 ϕ)u 2,(1 ϕ)v) l((1 ϕ)v) v V 0, Au 1 = B g g + f 1, g W, u 2 V 2, ϕ V, 0 ϕ 1 a.e. in, (5) where α,β > 0 are two regularization parameters used to tune the cost of ϕ and g, and f 1 is the same approximation to f introduced in (3). Solving (5) by an iterative procedure [5, 17] requires at each iteration to solve the two equations (4) for R 1 and R 2, and (3) for u 1. A less expensive optimal control problem is obtained by summing (4) for j = 1,2, and we obtain with R := R 1 + R 2 R,v = a(ϕu 1 + (1 ϕ)u 2,v) l(v) v V 0, (6) which is a Petrov-Galerkin type equation that we could have obtained directly applying a Petrov-Galerkin method to (2) using V PUM and V as trial and test spaces. Using (6), we get the less expensive optimal control problem min R,u 1,g,ϕ J(R,g,ϕ) := 1 2 R 2 V + α 2 ϕ 2 V + β 2 g 2 W s.t. R,v = a(ϕu 1 + (1 ϕ)u 2,v) l(v) v V 0, Au 1 = B g g + f 1, g W, u 2 V 2, ϕ V, 0 ϕ 1 a.e. in. (7)
4 4 Gabriele Ciaramella and Martin J. Gander A D M 2 1 B E C Fig. 1 Example of a boundary decomposition = Optimal control for elliptic boundary-layer problems As main test cases we consider elliptic problems of the form L (u) := µ u + a u + cu = f in, u = g on, (8) where is a bounded domain in R d, for d = 1,2, g C( ), f is sufficiently smooth, and the components of a are assumed to be strictly-positive. The assumption on a is restrictive, but it simplifies the presentation below and can be relaxed. The corresponding weak problem is to find a u { v H 1 () v = g on } such that a(u,v) := µ u v + a uv + cuvdx = f vdx =: l(v) v H0 1 (). We also assume that is such that the boundary can be decomposed into = 1 2, where the intersection 1 2 has a non-zero measure, as illustrated in Figure 1. To obtain u dd = ϕu 1 + (1 ϕ)u 2 u, we define Γ := \ 1 and introduce the operator L 1 := µ +c. Then, as in (3), for any choice of the control g H 1 0 (Γ ) the corresponding approximate problem for u 1 is µ u 1 v+cu 1 vdx = 0 v H0 1 (), { } u 1 w H 1 () w = g on 1, w = g + g on Γ, τ(w) C( ), where τ is the trace operator on. Notice that we have chosen f 1 = 0. As before, we introduce the operator A : H 1 () H 1 () defined as Au,v H 1,H 1 := µ u v + cuvdx for all v H1 0 (), and the operator B g : H 1 0 (Γ ) H 1 () such that v (B g g)(v) is a bounded linear functional in H 1 (). The operator B g represents the Dirichlet boundary conditions of (9). Au 1 = B g g is then equivalent to (9). The corresponding set V 1 is given by V 1 = {v H 1 () : Av = B g q for any q H 1 0 (Γ )}. Now, consider the operator L 2 := a + c and f 2 = f. The problem for u 2 is then L 2 (u 2 ) = a u 2 + cu 2 = f in, u 2 = g on 2, (10) (9)
5 Partition of Unity Methods for Heterogeneous Domain Decomposition 5 which we assume uniquely solvable in H 1 () C(). Notice that (10) is a pure advection problem and the subset 2 is given as the set of points where the characteristic curves enter the domain. This is the main assumption we make on 2 for the problem (10) to be well posed. The set V 2 contains only the solution to (10), i.e. V 2 = {u 2 }. The approximation u dd u is then obtained as u dd = ϕ 1 u 1 + ϕ 2 u 2, where the membership functions ϕ 1 = ϕ,ϕ 2 = 1 ϕ H 1 () form a partition of unity, and ϕ is such that {1} if x \ 2, ϕ(x) [0,1] if x 1 2, (11) {0} if x Γ, with τ(ϕ) C( ). Notice that this definition of ϕ makes u dd exact on the boundary, τ(u dd ) = τ(ϕ 1 u 1 + ϕ 2 u 2 ) = g. In what follows, we study the control problem (7) ((5) would have a similar structure) to optimize ϕ and g for computing the approximation u dd to the solution to (8). In particular, we first show well-posedness, and then we derive the firstorder optimality system. We consider directly a 2-dimensional problem (d = 2), since the analysis of the 1-dimensional version is simpler and relies on the same arguments. To define our optimal control problem, as in (7), we consider the cost functional J(R,g,ϕ) := 2 1 R 2 H 1 () + α 2 ϕ 2 H 1 () + β 2 g 2 H 1. Now, we introduce (Γ ) the control-to-state maps g u 1 (g) and (g,ϕ) R(u 1 (g),ϕ), where u 1 (g) and R(u 1 (g),ϕ) solve (9) and R,v H 1 () = µ u dd v + a u dd v + cu dd v f vdx v H 1 0 (). (12) Notice that the left-hand side of (12), that is R,v H 1 () = R v + Rvdx, is of a similar form to the left-hand side in (9). These maps are well defined according to the lemmas below and allow us to define the reduced cost functional J(g,ϕ) := J(R(u 1 (g),ϕ),g,ϕ) and the optimal control problem min g,ϕ J(g,ϕ) s.t. 0 ϕ(x) 1 in and (11) holds. (13) For well-posedness of this optimization problem, we need four Lemmas: Lemma 1. Let z H 1 ( ) with R 2 convex and Lipschitz. Then the problem µ u 1 v + cu 1 vdx = 0 v H0 1 () (14) with u 1 = z on is uniquely solvable by u 1 H 1 () C(), and there exists a positive constant c such that u 1 H 1 () c z H 1 ( ). Proof. To show that there exists a unique u 1 C(), we define w as the harmonic extension of z in. Recalling the embedding H 1 C for one-dimensional domains, we have that z C( ). Therefore, since is a Lipschitz domain,
6 6 Gabriele Ciaramella and Martin J. Gander w C 2 () C(); see, e.g., [12]. Now, consider the problem µ v + cv = cw in with v = 0 on. Since is convex, Theorems in [14] ensure that this problem is uniquely solved by v H 2 () H0 1(). Since R2, the Sobolev embedding H 2 () C() [7] ensures that v C(). Noticing that the function w + v solves (14), u 1 C() and is unique by the linearity of (14). Next, we show that u 1 H 1 () with u 1 H 1 () c z H 1 ( ). Consider the trace operator τ : H 1 () H 1/2 ( ). Since is a Lipschitz domain, by [16, Theorem 3.37, page 102] this operator has a bounded right-inverse τ 1 : H 1/2 ( ) H 1 (). Now, we define w := τ 1 z and note that w H 1 (). So, if we decompose u 1 as u 1 = w + ṽ, then ṽ must solve in a weak sense the problem µ ṽ + cṽ = ( µ w + cw) in with ṽ = 0 on. By the Lax-Milgram theorem we have that the unique solution is ṽ H0 1() and there exists a constant C such that ṽ H 1 () C w H 1 (). Therefore, u 1 H 1 () and using the decomposition u 1 = w+ṽ we get u 1 H 1 () (1 +C) w H 1 () = (1 +C) τ 1 z H 1 () K z H 1 ( ), for some positive constant K, where we used the boundedness of τ 1 [16]. Lemma 2. Let ϕ H 1 () such that 0 ϕ(x) 1 a.e. in. Then for any function v H 1 () C() it holds that ϕv H 1 (). Proof. An application of Theorem 1 in [9, page 247] shows that (vϕ) = v ϕ + ϕ v. Then a simple estimate of the norm (vϕ) L 2 () allows us to obtain the result. Lemma 3. Let {z n } n be a sequence that converges weakly in H 1 ( ) to a weak limit ẑ H 1 ( ), i.e. z n ẑ in H 1 ( ). Define the sequence {u 1,n } n by u 1,n := u 1 (z n ), where u 1 (z n ) solves (14) with u 1 = z n on. Then there exists a subsequence that converges weakly in H 1 () and strongly in L 2 () to the limit û 1 = u 1 (ẑ) H 1 (), i.e., û 1 in H 1 () and û 1 in L 2 (). Proof. Since the sequence {z n } n converges weakly in H 1 ( ), it is bounded in the norm H 1 ( ). By Lemma 1, we have that u 1,n H 1 () c z n H 1 ( ) K, for some positive constant K, and the sequence u 1,n is bounded in H 1 (). Since H 1 () is reflexive, there exists a weakly convergent subsequence û 1 in H 1 (). Now, from (14), we have that for any v H0 1() µ v + c vdx µ û 1 v + cû 1 vdx. Moreover, the weak convergence z n j ẑ and the continuity of the trace operator τ : H 1 () H 1/2 ( ) [16, Theorem 3.37] implies that z n j = τ( ) τ(û 1 ) = ẑ, weakly in H 1/2 ( ). Therefore, û 1 = u 1 (ẑ). We conclude by recalling the Sobolev compact embedding H 1 () L 2 (); see, e.g., [7]. Lemma 4. Let {u 1,n } n be the sequence defined in Lemma 3 such that û 1 (weakly) in H 1 (). Consider a sequence {ϕ n } n in H 1 () such that 0 ϕ n (x) 1
7 Partition of Unity Methods for Heterogeneous Domain Decomposition 7 and ϕ n ϕ (weakly) in H 1 () with 0 ϕ(x) 1. Then there exist two subsequences {ϕ n j } j and { } j such that ϕ n j ϕ and û 1 (strongly) in L 2 (), and for any v H 1 0 () (ϕ n j ) v + ϕ n j vdx ( ϕû 1 ) v + ϕû 1 vdx. Proof. The existence of the subsequences {ϕ n j } j and { } j such that ϕ n j ϕ and û 1 (strongly) in L 2 () follows from the fact that ϕ n ϕ (weakly in H 1 ()), Lemma 3, and the Sobolev (compact) embedding H 1 () L 2 () [7]. Now, recalling Lemma 1 and according to the proof of Lemma 2 it holds that ( ϕ n j ) = ϕ n j + ϕ n j. Therefore, to treat the products of sequences û 1.n j ϕ n j, ϕ n j, and ϕ n j, we use [7, Theorem ] to obtain for any v H0 1 () that (ϕ n j ) v+ϕ n j vdx= ϕ n j v+ϕ n j v+ϕ n j vdx û 1 ϕ v+ ϕ û 1 v+û 1 ϕvdx= ( ϕû 1 ) v+ ϕû 1 v+û 1 ϕvdx. We are now ready to prove that (13) is well posed. Theorem 1. Let α,β > 0, then there exists a solution to problem (13). Proof. Consider a minimizing sequence {(R n,ϕ n,u 1,n,g n )} n, where g n is extended by zero on. Since J is coercive in ϕ and g we have the bounds ϕ n H 1 () c and g n H 1 ( ) c, for two positive constants c,c ; see, e.g., [17]. The reflexivity of H 1 () and H 1 ( ) ensures the existence of weakly convergent subsequences: ϕ n j ϕ in H 1 () and g n j ĝ in H 1 ( ). By the Sobolev (compact) embedding H 1 () L 2 () [7], the sequence {ϕ n j } j converges strongly in L 2 () to ϕ. Since the set {v L 2 () : 0 v(x) 1 a.e. in } is (weakly) closed in L 2 () [17], we have 0 ϕ(x) 1. Consider now the sequence {u 1,n } n and the corresponding subsequence = u 1 (g n j ). By Lemma 3, we have that û 1 = u 1 (ĝ) weakly in H 1 () and û 1 = u 1 (ĝ) strongly in L 2 (). Consider the sequence {R n } n. Since R n satisfies R n,v H 1 () = µ u dd,n v + a u dd,n v + cu dd,n v f vdx v H0 1 (), where u dd,n = ϕ n u 1,n + (1 ϕ n )u 2, from the Lax-Milgram theorem we have that R n H 1 () K( u 1,n H 1 (), ϕ n H 1 () ), where the constant K depends on u 1,n H 1 () and ϕ n H 1 (), which are bounded. Therefore, R n is bounded as well, and by Lemma 4, one can show that R n j R = R(û 1, ϕ) weakly in H 1 (). Now, the weaklower semi-continuity of J implies the claim [17, 4]. To obtain the first-order optimality system, we rely on the Lagrange multiplier approach and work in the reduced space of solutions of constraint and adjoint equa-
8 8 Gabriele Ciaramella and Martin J. Gander tions; see, e.g., [5, 17]. We first recall the control-to-state maps g u 1 (g) and (g,ϕ) R(u 1 (g),ϕ) and the reduced cost functional J(g,ϕ). Then we notice that its derivatives, for δg H0 1(Γ ) and δϕ H1 0 (), are D g J(g,ϕ)(δg)= βg+r g,δg H 1 (Γ ), D ϕ J(w,ϕ)(δϕ)= αϕ +R ϕ,δϕ H 1 (). (15) Here R g is the solution of the problem R g,δg H 1 0 (Γ ) = B gδg,λ H 1,H1, (16) where, H 1,H 1 : H 1 () H0 1() R denotes the duality pairing, and R ϕ is the Riesz representative of the linear functional δϕ µ [ (u 1 u 2 )δϕ ] Rdx + a [ (u 1 u 2 )δϕ ] R + c(u 1 u 2 )δϕ Rdx. In (16), λ H0 1 () is a Lagrange multiplier that solves the adjoint equation λ v + cλ vdx = µ (vϕ) R + a (vϕ)r + cvϕ Rdx, (17) for all v H0 1 (). Therefore, the first-order optimality system is given by (9), (12), (17) and (16) together with the conditions [4, 17] D g J(g,ϕ)(δg) = 0, for all δg H0 1 (Γ ), and for any arbitrary θ > 0 ϕ = P Vad (ϕ θ ( αϕ + R ϕ ) ), where P Vad is the projection onto V ad := {v H 1 () : 0 v(x) 1 a.e. in }. 4 Numerical experiments We present now numerical experiments for the one-dimensional elliptic problem µ xx u x u = 1 in (0,1), with u(0) = 0, u(1) = 0, (18) for given µ = 0.01, computing u dd = ϕ 1 u 1 + ϕ 2 u 2, with µ xx u 1 = 0 in (0,1), u 1 (0) = 0, u 1 (1) = g, and x u 2 = 1 in [0,1), u 2 (1) = 0. We solve both the PUM and Petrov-Galerkin optimality systems discretized by linear finite-elements with a projected-lbfgs method with stopping tolerance
9 Partition of Unity Methods for Heterogeneous Domain Decomposition Petrov PUM Petrov PUM Exact Petrov PUM ϕ and 1 ϕ u and udd J(φ n,g n ) x x n (iteration) Fig. 2 Comparison of the Petrov and PUM approaches: Left: partition of unity functions ϕ and 1 ϕ. Middle: exact solution and approximations. Right: Decay of the cost functional. on the (relative) residual norm. The regularization parameters are α = β = In Figure 2 (left) we see that the ϕ and 1 ϕ obtained by the two approaches are very similar, and catch well the boundary layer on the left. The small bumps in the right part (close to x = 1) are due numerical effects and we checked that they disappear for smaller tolerances. In Figure 2 (middle) the exact solution is compared with the two approximations u dd, and we see good agreement. In Figure 2 (right), we show the decay of the cost functional with respect to the number of iterations, and we see that the Petrov-Galerkin approach converges a bit faster. References 1. Y. Achdou and O. Pironneau. The χ-method for the Navier-Stokes equations. IMA journal of numerical analysis, 13(4): , I. Babuska and J. M. Melenk. The partition of unity method. International Journal of Numerical Methods in Engineering, 40: , H. Berninger, E. Frnod, M. Gander, M. Liebendrfer, and J. Michaud. Derivation of the isotropic diffusion source approximation (idsa) for supernova neutrino transport by asymptotic expansions. SIAM Journal on Mathematical Analysis, 45(6): , A. Borzì, G. Ciaramella, and M. Sprengel. Formulation and Numerical Solution of Quantum Control Problems. Society for Industrial and Applied Mathematics, Philadelphia, PA, A. Borzì and V. Schulz. Computational Optimization of Systems Governed by Partial Differential Equations. SIAM, Philadelphia, F. Brezzi, C. Canuto, and A. Russo. A self-adaptive formulation for the Euler-Navier Stokes coupling. Comput. Methods Appl. Mech. Eng., 73: , P. G. Ciarlet. Linear and Nonlinear Functional Analysis with Applications. SIAM, Philadelphia, P. Degond and S. Jin. A smooth transition model between kinetic and diffusion equations. SIAM J. Numerical Analysis, 42(6): , L. C. Evans. Partial differential equations. Graduate studies in mathematics. American Mathematical Society, Providence (R.I.), M. J. Gander, L. Halpern, and V. Martin. A new algorithm based on factorization for heterogeneous domain decomposition. Numerical Algorithms, 73(1): , M.J. Gander and J. Michaud. Fuzzy domain decomposition: a new perspective on heterogeneous DD methods. In Domain Decomposition Methods in Science and Engineering XXI, pages Springer, D. Gilbarg and N. S. Trudinger. Elliptic partial differential equations of second order. Grundlehren der mathematischen Wissenschaften. Springer-Verlag, Berlin, New York, 1983.
10 10 Gabriele Ciaramella and Martin J. Gander 13. M. Griebel and M. A. Schweitzer. A particle-partition of unity method for the solution of elliptic, parabolic, and hyperbolic pdes. SIAM Journal on Scientific Computing, 22(3): , P. Grisvard. Elliptic Problems in Nonsmooth Domains. Monographs and studies in mathematics 24. Pitman Advanced Publishing Program, Boston, London, Melbourne, M.H. Holmes. Introduction to Perturbation Methods. Texts in Applied Mathematics. Springer New York, W.C.H. McLean. Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press, F. Tröltzsch. Optimal Control of Partial Differential Equations: Theory, Methods and Applications. Grad. Stud. Math American Mathematical Society, Providence, RI, L. A. Zadeh. Fuzzy sets. Information and Control, 8: , 1965.
Postprint.
http://www.diva-portal.org Postprint This is the accepted version of a chapter published in Domain Decomposition Methods in Science and Engineering XXI. Citation for the original published chapter: Gander,
More informationChapter 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 informationNumerical 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 informationThe Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control
The Navier-Stokes problem in velocity-pressure formulation :convergence and Optimal Control A.Younes 1 A. Jarray 2 1 Faculté des Sciences de Tunis, Tunisie. e-mail :younesanis@yahoo.fr 2 Faculté des Sciences
More 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 informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 informationChapter 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 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 informationFINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET OPTIMAL CONTROL PROBLEMS
Numerical Functional Analysis and Optimization, 28(7 8):957 973, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0163-0563 print/1532-2467 online DOI: 10.1080/01630560701493305 FINITE ELEMENT APPROXIMATION
More informationGlowinski Pironneau method for the 3D ω-ψ equations
280 GUERMOND AND QUARTAPELLE Glowinski Pironneau method for the 3D ω-ψ equations Jean-Luc Guermond and Luigi Quartapelle 1 LIMSI CNRS, Orsay, France, and Dipartimento di Fisica, Politecnico di Milano,
More information1 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 informationAdvection Diffusion Problems with Pure Advection Approximation in Subregions
Advection Diffusion Problems with Pure Advection Approximation in Subregions M. J. Gander, L. Halpern 2, C. Japhet 2, and V. Martin 3 Université de Genève, 2-4 rue du Lièvre, CP 64, CH-2 Genève. Switzerland.
More informationA Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains
A Remark on the Regularity of Solutions of Maxwell s Equations on Lipschitz Domains Martin Costabel Abstract Let u be a vector field on a bounded Lipschitz domain in R 3, and let u together with its divergence
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationWeak 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 informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationTheory 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 informationFact 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 informationSTOKES 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 informationand finally, any second order divergence form elliptic operator
Supporting Information: Mathematical proofs Preliminaries Let be an arbitrary bounded open set in R n and let L be any elliptic differential operator associated to a symmetric positive bilinear form B
More informationContents: 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 informationTechnische Universität Graz
Technische Universität Graz Stability of the Laplace single layer boundary integral operator in Sobolev spaces O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2016/2 Technische
More informationLeast-Squares Finite Element Methods
Pavel В. Bochev Max D. Gunzburger Least-Squares Finite Element Methods Spri ringer Contents Part I Survey of Variational Principles and Associated Finite Element Methods 1 Classical Variational Methods
More informationON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS
Chin. Ann. Math.??B(?), 200?, 1 20 DOI: 10.1007/s11401-007-0001-x ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS Michel CHIPOT Abstract We present a method allowing to obtain existence of a solution for
More informationLecture 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 informationA NOTE ON THE INITIAL-BOUNDARY VALUE PROBLEM FOR CONTINUITY EQUATIONS WITH ROUGH COEFFICIENTS. Gianluca Crippa. Carlotta Donadello. Laura V.
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume, Number, Xxxx XXXX pp. A NOTE ON THE INITIAL-BOUNDARY VALUE PROBLEM FOR CONTINUITY EQUATIONS WITH ROUGH COEFFICIENTS Gianluca
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationLectures 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 informationThe Dirichlet-to-Neumann operator
Lecture 8 The Dirichlet-to-Neumann operator The Dirichlet-to-Neumann operator plays an important role in the theory of inverse problems. In fact, from measurements of electrical currents at the surface
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 informationTraces, 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 information5. The Indirect Approach To Domain Decomposition
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c2003 DDM.org 5. The Indirect Approach To Domain Decomposition
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.
ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS Emerson A. M. de Abreu Alexandre N. Carvalho Abstract Under fairly general conditions one can prove 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 informationOptimal Control of Partial Differential Equations I+II
About the lecture: Optimal Control of Partial Differential Equations I+II Prof. Dr. H. J. Pesch Winter semester 2011/12, summer semester 2012 (seminar), winter semester 2012/13 University of Bayreuth Preliminary
More informationRegularity 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 informationON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume, Number, Pages S -9939(XX- ON THE DYNAMICAL SYSTEMS METHOD FOR SOLVING NONLINEAR EQUATIONS WITH MONOTONE OPERATORS N. S. HOANG AND A. G. RAMM (Communicated
More informationChapter 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 informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More 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 informationUniversité de Metz. Master 2 Recherche de Mathématiques 2ème semestre. par Ralph Chill Laboratoire de Mathématiques et Applications de Metz
Université de Metz Master 2 Recherche de Mathématiques 2ème semestre Systèmes gradients par Ralph Chill Laboratoire de Mathématiques et Applications de Metz Année 26/7 1 Contents Chapter 1. Introduction
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
More informationSecond 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 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 informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationWeak 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 informationREGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS
REGULAR LAGRANGE MULTIPLIERS FOR CONTROL PROBLEMS WITH MIXED POINTWISE CONTROL-STATE CONSTRAINTS fredi tröltzsch 1 Abstract. A class of quadratic optimization problems in Hilbert spaces is considered,
More informationCONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM
CONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM IAN KNOWLES AND MARY A. LARUSSA Abstract. For the elliptic equation (p(x) v) + λq(x)v = f, x R n, the problem of determining when one or more
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationNONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian UNDER NONHOMOGENEOUS NEUMANN BOUNDARY CONDITION
Electronic Journal of Differential Equations, Vol. 2016 (2016), No. 210, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu NONLINEAR FREDHOLM ALTERNATIVE FOR THE p-laplacian
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationu = f in Ω, u = q on Γ. (1.2)
ERROR ANALYSIS FOR A FINITE ELEMENT APPROXIMATION OF ELLIPTIC DIRICHLET BOUNDARY CONTROL PROBLEMS S. MAY, R. RANNACHER, AND B. VEXLER Abstract. We consider the Galerkin finite element approximation of
More informationSPECTRAL REPRESENTATIONS, AND APPROXIMATIONS, OF DIVERGENCE-FREE VECTOR FIELDS. 1. Introduction
SPECTRAL REPRESENTATIONS, AND APPROXIMATIONS, OF DIVERGENCE-FREE VECTOR FIELDS GILES AUCHMUTY AND DOUGLAS R. SIMPKINS Abstract. Special solutions of the equation for a solenoidal vector field subject to
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationRegularization on Discrete Spaces
Regularization on Discrete Spaces Dengyong Zhou and Bernhard Schölkopf Max Planck Institute for Biological Cybernetics Spemannstr. 38, 72076 Tuebingen, Germany {dengyong.zhou, bernhard.schoelkopf}@tuebingen.mpg.de
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationWEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS
WEAK LOWER SEMI-CONTINUITY OF THE OPTIMAL VALUE FUNCTION AND APPLICATIONS TO WORST-CASE ROBUST OPTIMAL CONTROL PROBLEMS ROLAND HERZOG AND FRANK SCHMIDT Abstract. Sufficient conditions ensuring weak lower
More informationDivergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints
Divergence Boundary Conditions for Vector Helmholtz Equations with Divergence Constraints Urve Kangro Department of Mathematical Sciences Carnegie Mellon University Pittsburgh, PA 15213 Roy Nicolaides
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationDomain 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 informationDiscreteness of Transmission Eigenvalues via Upper Triangular Compact Operators
Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research
More informationA SHORT PROOF OF INCREASED PARABOLIC REGULARITY
Electronic Journal of Differential Equations, Vol. 15 15, No. 5, pp. 1 9. ISSN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu A SHORT PROOF OF INCREASED
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationBIHARMONIC WAVE MAPS INTO SPHERES
BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.
More informationGeneralized Korn s inequality and conformal Killing vectors
Generalized Korn s inequality and conformal Killing vectors arxiv:gr-qc/0505022v1 4 May 2005 Sergio Dain Max-Planck-Institut für Gravitationsphysik Am Mühlenberg 1 14476 Golm Germany February 4, 2008 Abstract
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 informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationDISCONTINUOUS GALERKIN FINITE ELEMENT DISCRETIZATION FOR STEADY STOKES FLOWS WITH THRESHOLD SLIP BOUNDARY CONDITION
DISCONTINUOUS GALRKIN FINIT LMNT DISCRTIZATION FOR STADY STOKS FLOWS WITH THRSHOLD SLIP BOUNDARY CONDITION J.K. Djoko Department of Mathematics and Applied Mathematics, University of Pretoria, Private
More informationInfinite-Dimensional Dynamical Systems in Mechanics and Physics
Roger Temam Infinite-Dimensional Dynamical Systems in Mechanics and Physics Second Edition With 13 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vii ix GENERAL
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationCONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM
CONDITIONAL WELL-POSEDNESS FOR AN ELLIPTIC INVERSE PROBLEM IAN KNOWLES AND MARY A. LARUSSA Abstract. For the elliptic equation (p(x) v) + λq(x)v = f, x R n, the problem of determining when one or more
More informationA CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS 1. INTRODUCTION
Trends in Mathematics Information Center for Mathematical Sciences Volume 4, Number 2, December 21, Pages 39 44 A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS
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 informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationPDE Constrained Optimization selected Proofs
PDE Constrained Optimization selected Proofs Jeff Snider jeff@snider.com George Mason University Department of Mathematical Sciences April, 214 Outline 1 Prelim 2 Thms 3.9 3.11 3 Thm 3.12 4 Thm 3.13 5
More informationA MIXED VARIATIONAL FRAMEWORK FOR THE RADIATIVE TRANSFER EQUATION
SpezialForschungsBereich F 32 Karl Franzens Universita t Graz Technische Universita t Graz Medizinische Universita t Graz A MIXED VARIATIONAL FRAMEWORK FOR THE RADIATIVE TRANSFER EQUATION H. Egger M. Schlottbom
More informationON THE EXISTENCE OF TRANSMISSION EIGENVALUES. Andreas Kirsch1
Manuscript submitted to AIMS Journals Volume 3, Number 2, May 29 Website: http://aimsciences.org pp. 1 XX ON THE EXISTENCE OF TRANSMISSION EIGENVALUES Andreas Kirsch1 University of Karlsruhe epartment
More informationTHE STOKES SYSTEM R.E. SHOWALTER
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
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES
Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 29, 2004, 167 176 ON BOUNDEDNESS OF MAXIMAL FUNCTIONS IN SOBOLEV SPACES Piotr Haj lasz and Jani Onninen Warsaw University, Institute of Mathematics
More informationTechnische Universität Graz
Technische Universität Graz A note on the stable coupling of finite and boundary elements O. Steinbach Berichte aus dem Institut für Numerische Mathematik Bericht 2009/4 Technische Universität Graz A
More informationAffine covariant Semi-smooth Newton in function space
Affine covariant Semi-smooth Newton in function space Anton Schiela March 14, 2018 These are lecture notes of my talks given for the Winter School Modern Methods in Nonsmooth Optimization that was held
More informationCompressible Euler equations with second sound: asymptotics of discontinuous solutions. Beixiang Fang Reinhard Racke
Universität Konstanz Compressible Euler equations with second sound: asymptotics of discontinuous solutions Beixiang Fang Reinhard Racke Konstanzer Schriften in Mathematik Nr. 306, August 0 ISSN 430-3558
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationOn Two Nonlinear Biharmonic Evolution Equations: Existence, Uniqueness and Stability
On Two Nonlinear Biharmonic Evolution Equations: Existence, Uniqueness and Stability Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We study the following two nonlinear evolution equations with a fourth
More informationThe 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 informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 155 A posteriori error estimates for stationary slow flows of power-law fluids Michael Bildhauer,
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationIntroduction. Christophe Prange. February 9, This set of lectures is motivated by the following kind of phenomena:
Christophe Prange February 9, 206 This set of lectures is motivated by the following kind of phenomena: sin(x/ε) 0, while sin 2 (x/ε) /2. Therefore the weak limit of the product is in general different
More informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More 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 informationNumerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods
Numerical Analysis of Higher Order Discontinuous Galerkin Finite Element Methods Contents Ralf Hartmann Institute of Aerodynamics and Flow Technology DLR (German Aerospace Center) Lilienthalplatz 7, 3808
More informationt y n (s) ds. t y(s) ds, x(t) = x(0) +
1 Appendix Definition (Closed Linear Operator) (1) The graph G(T ) of a linear operator T on the domain D(T ) X into Y is the set (x, T x) : x D(T )} in the product space X Y. Then T is closed if its graph
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More information