A SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS

Size: px
Start display at page:

Download "A SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS"

Transcription

1 INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS HONGFEI FU, HONGING RUI, HUI GUO, JIANSONG ZHANG, AND JIAN HOU Abstract. In tis paper, we propose a splitting least-squares mixed finite element metod for te approximation of elliptic optimal control problem wit te control constrained by pointwise inequality. By selecting a properly least-squares minimization functional, we derive equivalent two independent, symmetric and positive definite weak formulation for te primal state variable and its flux. Ten, using te first order necessary and also sufficient optimality condition, we deduce anoter two corresponding adjoint state equations, wic are bot independent, symmetric and positive definite. Also, a variational inequality for te control variable is involved. For te discretization of te state and adjoint state equations, eiter RT mixed finite element or standard C finite element can be used, wic is not necessary subject to te Ladyzenkaya-Babuska-Brezzi condition. Optimal a priori error estimates in corresponding norms are derived for te control, te states and adjoint states, respectively. Finally, we use some numerical examples to validate te teoretical analysis. Key words. Optimal control, splitting least-squares, mixed finite element metod, positive definite, a priori error estimates.. Introduction Optimal control problems are playing an increasingly important role in modern scientific and engineering numerical simulations. Nowadays, finite element metod seems to be te most widely used numerical metod in practical computation. Te readers are referred to, for example, Refs. [, 2, 3, 4, 5] for systematic introductions of finite element metods and optimal control problems. In tis paper, we are interested in te following convex quadratic optimal control problem wit te control constrained by pointwise inequality: min u U ad J y, σ, u = 2 subject to 2 and 3 y y d 2 + Ω Ω σ σ d 2 + γ divσ + cy = f + Bu, in Ω, σ + A y =, in Ω, y =, on Ω, ξ ux ξ 2, a.e. in Ω U. Ω U u 2 Here γ > is a constant, Ω and Ω U Ω are two bounded domain in R 2, wit Lipscitz boundaries Ω and Ω U. A precise formulation of tis problem including a functional analytic setting is given in te next section. Received by te editors June 22, 25 and, in revised form, December 3, Matematics Subject Classification. 49K2, 49M25, 65N5, 65N3. Corresponding autor: Hongxing Rui. Tis researc was supported by te NSF of Cina 2485, 9336, 57367, 4588, te NSF of Sandong Province BS23NJ, and te Fundamental Researc Funds for te Central Universities 4C227, 5C84A. 6

2 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 6 Classical mixed finite element metods see Refs. [6, 7, 8, 9, ] ave been proved effectively for solving elliptic equations and fluid problems. Tey ave an advantage of approximating te unknown scalar variable and its diffusive flux simultaneously. Besides, tese metods can approximate te unknown variable and its flux to a same order of accuracy. Recently, tere are some researc articles on tese metods for solving optimal control problems, see Refs. [, 2, 3], for example. However, it is well known tat tese metods usually produce a symmetric but indefinite system for elliptic equations. Tus te popular conjugate gradient CG or algebraic multi-grid AMG solvers can not be used for te solution of linear algebraic equation systems. To conquer tese difficulties appeared in using classical mixed finite element metods, least-squares mixed finite element metod, for first-order elliptic mixed system in unknown variable y and unknown velocity flux σ, was introduced by Pelivanov et al. [4]. It is well known tat te least-squares mixed finite element metod as two typical advantages: First, it is not subjected to te Ladyzenkaya- Babuska-Brezzi consistency condition, so te coice of finite element spaces becomes flexible; Second, it results in a symmetric and positive definite system, wic can be solved using tose solvers suc as CG and AMG quickly. Te idea of splitting least-squares was first proposed by Rui et al. in [5] for a reaction-diffusion equation, were by selecting a properly least-squares functional, te autors derived two independent, symmetric and positive definite equations, respectively, for te unknown state variable y and its flux σ. Ten it is applied to solve linear and nonlinear parabolic equations [6], sobolev equations [7], pseudo-parabolic equations [8], and nonlinear convection-diffusion equations [9] and so on. In tis paper, we apply te splitting least-squares mixed finite element metod for te discretization of elliptic optimal control problem. Pointwise inequality constraints on te control variable are considered. We derive optimal a priori error estimates, respectively, for te optimal control u in L 2 Ω U -norm, wic is approximated by piecewise constant or piecewise linear discontinuous elements; for te primal state y and adjoint state z bot in L 2 Ω-norm and H Ω-norm, wic are approximated by standard piecewise linear C finite elements; for te flux state σ and adjoint state ω in Hdiv; Ω-norm, wic are approximated by te lowewt-order RT mixed finite elements or standard piecewise linear C finite elements. Here, te Ladyzenkaya-Babuska-Brezzi consistency condition for te discretization spaces of y and σ is not needed. Tis paper is organized as follows. In Sect. 2, we introduce te optimal control problem and derive te continuous optimality conditions based on te idea of leastsquares. In Sect. 3, a splitting least-squares mixed finite element approximation to te continuous optimal control problem is proposed, and ten we derive te corresponding discrete optimality conditions. In Sect. 4, some a priori error estimates for te states, adjoint states and control are derived under control constrained by pointwise inequality. In Sect. 5, we conduct some numerical experiments to observe te convergence beavior of te numerical sceme. In te last section, some concluding remarks are given. In te following, we employ te standard notations W m,p Ω for Sobolev spaces on Ω wit norm W m,p Ω and seminorm W m,p Ω. We set W m,p Ω = v W m,p Ω : v = on Ω}. For p = 2, we denote H m Ω = W m,2 Ω and H m Ω = W m,2 Ω. In addition, C denotes a general positive constant wic is independent of te spatial mes parameters.

3 62 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU 2. Optimality system based on splitting least-squares In tis section, we first describe te optimal control problem under consideration, and ten derive te continuous optimality system based on te idea of splitting least-squares. Let Hdiv; Ω = τ L 2 Ω 2 : divτ L 2 Ω}, endowed wit te norm given by τ Hdiv;Ω = τ 2 L 2 Ω 2 + divτ 2 L 2 Ω /2. To formulate te optimal control problem, in te rest we sall take te state spaces V = H Ω, W = Hdiv; Ω, te control space and te admissible control set 4 U = L 2 Ω U, U ad = u U : ξ u ξ 2, a.e. in Ω U }, were te bounds ξ, ξ 2 R fulfill ξ < ξ 2. For te optimal control problem -2, we need te following assumptions. H- Te desired states and source function satisfy te following regularity: y d H Ω, σ d H Ω 2, f L 2 Ω. H-2 A = Ax = a i,j x 2 2 is a bounded symmetric and positive definite matrix, i.e., tere are positive constants α and β suc tat α 2 T A β 2, R 2. H-3 c = cx is positive definite and bounded, tat is, tere exist positive constants c and c 2 suc tat < c c c 2. H-4 B is a linear continuous operator from U = L 2 Ω U to L 2 Ω, suc tat Bu, v = B v, u U C u L2 Ω U v L2 Ω, u U, v V, were B is te adjoint operator of B. Te classical mixed weak formula for problems -3 is to find y, σ, u L 2 Ω Hdiv; Ω L 2 Ω U suc tat 5 J y, σ, u = min u U ad J y, σ, u and 6 A σ, τ y, divτ =, τ Hdiv; Ω, divσ, v + cy, v = f + Bu, v, v L 2 Ω. It is clear tat te two state variables σ and y in problem 6 are coupled eac oter, and te discretization of 6 will lead to a saddle-point problem, wic is symmetric but indefinite. To derive a symmetric and positive definite weak formulation for problem 2, we define for a given control u te least-squares functional Fv, τ as follows: 7 Fv, τ = 2 c /2 divτ + cv f Bu 2 L 2 Ω + 2 A /2 τ + A v 2 L 2 Ω 2.

4 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 63 Ten te solution to 2 is corresponding to te least-squares minimization problem: 8 Fy, σ = min Fv, τ. v,τ V W Furtermore, it follows from te first order necessary and sufficient minimization condition, problem 8 becomes: find y, σ V W suc tat A y, v + cy, v = f + Bu, v, v V, 9 c divσ, divτ + A σ, τ = c f + Bu, divτ, τ W. Remark 2.. It is clear tat te two equations in problem 9 are split, and tus can be solved independently. Besides, it is easy to ceck tat tey are also symmetric, coercive and continuous. For te given control set U ad, we now restate te optimal control problem -3 as follows: OCP J y, σ, u = min u U ad J y, σ, u suc tat y, σ, u V W U and A y, v + cy, v = f + Bu, v, v V, c divσ, divτ + A σ, τ = c f + Bu, divτ, τ W. Here te inner product in L 2 Ω or L 2 Ω 2 is indicated by,. Since te objective functional J y, σ, u is convex, it ten follows from Ref. [2] tat te optimal control problem OCP as a unique solution y, σ, u V W U ad. Furtermore, y, σ, u is te solution of OCP if and only if tere is a pair of adjoint state z, ω V W, suc tat y, σ, z, ω, u V W 2 U ad satisfies te following optimality conditions: OCP-OPT A y, v + cy, v = f + Bu, v, v V, c divσ, divτ + A σ, τ = c f + Bu, divτ, τ W. 2 A z, v + cz, v = y y d, v, v V, c divω, divτ + A ω, τ = σ σ d, τ, τ W. 3 γu B z B c divω, u u U, u U ad U, were, U denotes te inner product in U = L 2 Ω U. Remark 2.2. It is clear tat te two variables z and ω in problem 2 are also independent, so tey can be solved separately, too. Besides, te two equations in problem 2 are also symmetric, coercive and continuous. Remark 2.3. [2]. Inequality 3 is equivalent to te following: γu B z B c divω, u = ξ ; 4 γu B z B c divω, u = ξ 2 ; γu B z B c divω =, ξ < u < ξ 2. Define a pointwise projection operator P Uad from U to U ad tat P Uad fx, t = max ξ, minξ 2, fx, t. Ten, te optimality condition 3 can be expressed as 5 u = P Uad γ B z + c divω.

5 64 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU 3. Splitting least-squares mixed finite element approximation In tis section, we consider te approximation of optimal control problem OCP based on te splitting least-squares mixed finite element metod. Let T y and T σ be two families of regular triangulations of Ω, wic could be identical or not. Here y and σ denote te largest diameters of elements in T y and T σ, respectively. Construct two finite element spaces V W V W wit te following approximate properties see Refs. [, 8, 9, ]: Tere exist two integers k k and l suc tat 6 inf v V v v L2 Ω + y v v L2 Ω 2} Cl+ inf τ W τ τ L2 Ω 2 Ck+ σ τ H k+ Ω 2, inf τ W divτ τ L2 Ω C k σ τ H k + Ω 2, y v H l+ Ω, for any v V H l+ Ω and τ W H k+ Ω 2. Here k = k + wen W is cosen as RT mixed finite elements in Ref. [8], and k = k wen W is any one of te oter mixed elements in Refs. [9, ] and standard C elements in Ref. []. Let T u be a family of regular triangulations of Ω U, were u denotes te largest diameters of elements in T u. Associate wit T u is anoter finite element space U U = L 2 Ω U consist of piecewise r-order polynomials on eac element of T u. Here tere is no requirement for te continuity, due to te limited regularity of te optimal control. Let Uad be a closed convex set in U, tat is, 7 U ad = u U : ξ u ξ 2, a.e. in Ω U }. Ten te corresponding splitting least-squares mixed finite element approximation of OCP, wic will be labeled as OCP, is defined as follows: 8 J y, σ, u = min J y, σ, u u Uad were y, σ, u V W U satisfies A y, v + cy, v = f + Bu, v, v V, c divσ, divτ + A σ, τ = c f + Bu, divτ, τ W. It is again well known see, e.g., Ref. [2] tat te optimal control problem OCP as a unique solution y, σ, u V W Uad, and tat a triplet y, σ, u is te solution of OCP if and only if tere is a pair of adjoint state z, ω V W, suc tat y, σ, z, ω, u V W 2 Uad satisfies te following discrete optimality conditions: OCP-OPT 9 A y, v + cy, v = f + Bu, v, v V, c divσ, divτ + A σ, τ = c f + Bu, divτ, τ W. 2 A z, v + cz, v = y y d, v, v V, c divω, divτ + A ω, τ = σ σ d, τ, τ W. 2 γu B z B c divω, u u U, u U ad U. To obtain a priori error estimates for te above numerical scemes 9-2, in te end of tis section, we introduce auxiliary variables y u, σ u V W

6 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 65 and z u, ω u V W, associated wit te optimal control u as follows: A y u, v + cy u, v = f + Bu, v, 22 and 23 c divσ u, divτ + A σ u, τ = c f + Bu, divτ, A z u, v + cz u, v = y u y d, v, c divω u, divτ + A ω u, τ = σ u σ d, τ, for v V and τ W. For simplicity, in te following we denote ỹ := y u, z := z u, σ := σ u, ω := ω u. 4. A priori error estimates In tis section, we are able to derive some a priori error estimates for te optimal control problem OCP-OPT and its splitting least-squares mixed finite element approximation OCP-OPT. Due to te weak regularity of te optimal control problem, we consider piecewise constant elements r = or piecewise linear discontinuous elements r = for te approximation of te optimal control u, and piecewise linear C elements l = are adopted for bot te approximations of te states y and z. Wile for te discretizations of te flux states σ and ω, we use eiter te lowest-order RT mixed finite elements k =, k = or standard piecewise linear C elements k = k =, were in bot cases te Ladyzenkaya- Babuska-Brezzi consistency condition is not satisfied. Let Ω U = τ U : τ U Ω U, ξ < u τu < ξ 2 }, 24 Ω a U = τ U : τ U Ω U, u τu ξ, or u τu ξ 2 }, Ω b U = Ω U \ Ω U Ω a U. It is easy to ceck tat te tree parts do not intersect on eac oter, and Ω U = Ω U Ωa U Ωb U. In tis paper, we assume tat u and T U are regular suc tat measω b U C u see Refs. [2, 22]. Moreover, set Ω U = x Ω U, ξ < ux < ξ 2 }. Ten it is easy to see tat Ω U Ω U. Before deriving te main error estimates for te splitting least-squares mixed finite element approximation of optimal control problem governed by elliptic equations, some interpolation and projection estimates are prepared witout proof. Lemma 4.. []. Let P be te L 2 -projection from U = L 2 Ω U to U suc tat for any u U 25 u P u, ϕ U =, ϕ U. Ten we ave for any u H Ω U. u P u L2 Ω U C u u H Ω U Lemma 4.2. Let I be te standard Lagrange interpolation operator defined in Ref. []. Ten tere is a constant C > suc tat u I u W m,p Ω C 2 m u u W 2,p Ω for u W 2,p Ω, p >, and m = or.

7 66 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU Lemma 4.3. Let R y V be defined as te elliptic projection of y V suc tat 26 A y R y, v + cy R y, v =, w V. It ten follows from Ref. [2] tat y R y L 2 Ω + y y R y H Ω C 2 y y H 2 Ω. Lemma 4.4. [5]. From te approximate property 6 of finite element spaces we know tat for any σ W H 2 Ω 2, tere exists a vector-valued function Q σ W suc tat σ Q σ L2 Ω 2 Ck+ σ σ H k+ Ω 2, divσ Q σ L2 Ω C σ σ H2 Ω 2, were k = for te lowest-order RT mixed finite elements, and k = for piecewise linear C elements. In te following, we wis to demonstrate tree main lemmas wic can be collected to obtain te main conclusion of tis paper. Lemma 4.5. Let y, σ, ω, z and ỹ, σ, ω, z be te solutions of 9-2 and 22-23, respectively. Ten te following estimates old 27 y ỹ H Ω + z z H Ω C u u L 2 Ω U, 28 σ σ Hdiv;Ω + ω ω Hdiv;Ω C u u L 2 Ω U. Proof. It follows from 9 and 22 tat 29 A y ỹ, v + cy ỹ, v = Bu u, v, c divσ σ, divτ + A σ σ, τ = c Bu u, divτ, for v V and τ W. Note tat te two equations in problem 29 are split. By selecting v = y ỹ, τ = σ σ, respectively, and using te assumptions H-2-H-4, we can easily prove tat 3 y ỹ H Ω C u u L 2 Ω U, 3 σ σ Hdiv;Ω C u u L 2 Ω U. Similarly, one can obtain from 2 and 23 tat 32 A z z, v + cz z, v = ỹ y, v, c divω ω, divτ + A ω ω, τ = σ σ, τ, for v V and τ W. Ten wit v = z z and τ = ω ω in 32, and again using te assumptions H-2-H-4, we similarly derive 33 z z H Ω C y ỹ L2 Ω C u u L 2 Ω U, 34 ω ω Hdiv;Ω C σ σ L2 Ω 2 C u u L 2 Ω U, wic ends te proof of Lemma 4.5.

8 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 67 Lemma 4.6. Let y, σ, z, ω, u and y, σ, z, ω, u be te solutions of OCP- OPT and OCP-OPT, respectively. Let U be te piecewise constant element space r =. Assume tat u H Ω U, z H Ω, and divω H Ω. Ten we ave 35 u u L2 Ω U C u + z z L2 Ω + divω ω L2 Ω, were z, ω is defined in 23, and te constant C depends on some norms of te solution z, ω, u, but is independent of te mes parameters u. Let U be te piecewise linear discontinuous element space r =. Furtermore, assume tat u W, Ω U H 2 Ω U and divω W, Ω. Ten we ave 36 u u L 2 Ω U C, were Ω U is defined above, z W, Ω, 3 2 u + z z L 2 Ω + divω ω L 2 Ω were te constant C depends on some norms of te solution z, ω, u, but is independent of te mes parameter u. Proof. Let u I U ad be an approximation of u. Note tat U U and U ad U ad. Recalling te variational inequalities 3 and 2 imply tat 37 Hence, we can deduce tat γu B z B c divω, u u U, γu B z B c divω, u u I U. γ u u 2 L 2 Ω U = γu, u u U γu, u u U = γu B z B c divω, u u U +γu B z B c divω, u u I U, 38 +γu B z B c divω, u I u U +B z z, u u U + B c divω ω, u u γu u, u I u U + γu B z B c divω, u I u U U +B z z, u I u U + B c divω ω, u I u U +B z z, u I u U + B c div ω ω, u I u U +B z z, u u U + B c divω ω, u u U +B z z, u u U + B c div ω ω, u u Te rigt-and side of 38 can be bounded in a standard way. We ere just pay special attentions on te last two terms of te rigt-and side of 38. By selecting v = z z, τ = ω ω in 29 and v = y ỹ, τ = σ σ in 32 we ave 39 and B z z, u u U = z z, Bu u = A y ỹ, z z + cy ỹ, z z = y ỹ, y ỹ = y ỹ 2 L 2 Ω, B c div ω ω, u u U = c div ω ω, Bu u U. 4 = c divσ σ, divω ω + A σ σ, ω ω

9 68 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU = σ σ, σ σ = σ σ 2 L 2 Ω 2. Tus we obtain from 39-4, Lemma 4.5 and Caucy-Scwarz inequality tat γ u u 2 L 2 Ω U γu B z B c divω, u I u U + Cδ u u I 2 L 2 Ω U +Cδ z z 2 L 2 Ω + Cδ divω ω 2 L 2 Ω 2 4 +Cδ z z 2 L 2 Ω + Cδ div ω ω 2 L 2 Ω + Cδ u u 2 L 2 Ω U γu B z B c divω, u I u U + C u u I 2 L 2 Ω U +C z z 2 L 2 Ω + C divω ω 2 L 2 Ω + Cδ u u 2 L 2 Ω U, were δ is an arbitrary small positive number. Tanks to te definition of te control finite element space U, in te following we sall discuss two different cases for estimates of te first two terms on te rigtand side of 4. First let us consider te case tat U is te piecewise constant element space. Let u I be defined as te L2 -projection P u see, Lemma 4.. It is clear tat P u U ad, and 42 u P u L 2 Ω U C u u H Ω U. Moreover, if u H Ω U, z H Ω, and divω H Ω, we conclude from te definition of P in 25 tat γu B z B c divω, P u u U 43 = Θ P Θ P u u Ω U divω 2 H Ω + z 2 H Ω + u 2 H Ω U, C 2 u were Θ := γu B z B c divω. Ten 35 follows from 4-43 wit Cδ = γ/2. Next we consider te case tat U is te piecewise linear element space. Let u I = I u U be defined as te standard Lagrange linear interpolation of u see, Lemma 4.2. It is clear tat I u Uad, and for u W, Ω U H 2 Ω U we ave u I u L 2 Ω U C 2 u u H 2 Ω U, u I u L Ω b U C u u W, Ω b U. Note tat I u = u on Ω a U, ten it follows tat u I u 2 L 2 Ω U = + + u I u 2 Ω U Ω a U Ω b U 44 C 4 u u 2 H 2 Ω U + + C2 u u 2 W, Ω b U measωb U C 3 u u 2 H 2 Ω U + u 2 W, Ω U. Moreover, it follows from 4 tat γu B z B c divω = on Ω U, and we conclude from te definition of Ω b U tat for any element τ U Ω b U, tere is a x suc

10 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 69 tat ξ < ux < ξ 2, and ence γu B z B c divω x =. Terefore, for any τ U Ω b U we ave Tus γu B z B c divω L τ U = γu B z B c divω γu B z B c divω x L τ U C u γu B z B c divω W, τ U. 45 γu B z B c divω, I u u U = + + γu B z B c divω I u u Ω U Ω a U Ω b U = + + γu B z B c divω I u u Ω b U γu B z B c divω L Ω b U I u u L Ω b U measω b U C 3 u u 2 W, Ω U + z 2 W, Ω + divω 2 W, Ω. Ten 36 is proved by inserting into 4 wit Cδ = γ/2. Lemma 4.7. Let y, σ, z, ω and ỹ, σ, z, ω be te solutions of -2 and 22-23, respectively. Assume tat problem 2 is H 2 -regular, and te solutions y, z H Ω H 2 Ω and σ, ω H 2 Ω 2. Ten te following estimates old 46 y ỹ Hs Ω + z z Hs Ω C 2 s y, s =,, 47 σ σ Hdiv;Ω + ω ω Hdiv;Ω C σ, were te constant C depends on some norms of y, σ, z, ω, but is independent of te mes parameters y and σ. Proof. In te following, we divide te proof into two parts. Part I. Note tat σ, ỹ is te least-squares mixed finite element solution of σ, y. Terefore, it follows immediately from Ref. [5] combined wit Lemmas 4.3 and 4.4 tat 48 y ỹ Hs Ω C 2 s y y H2 Ω, s =,, and 49 σ σ Hdiv;Ω C σ σ H2 Ω 2. Part II. To derive corresponding error estimates for z z Hs Ω and ω ω Hdiv;Ω in 46-47, we write z z = z R z + R z z := ρ + π, ω ω = ω Q ω + Q ω ω := ζ + η. Here R z V and Q ω W are defined as in Lemmas , and of course ρ and ζ are bounded in te desired way. In te following, we first estimate π = R z z. By subtracting te first equation in 23 wit tat of 2, we obtain for v V, 5 A z z, v + cz z, v = ỹ y, v,

11 62 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU wic implies 5 A π, v + cπ, v = ỹ y, v, v V. Let v = π = R z z. It ten follows from te standard finite element error estimates and 48 tat 52 R z z H Ω C y ỹ L2 Ω C 2 y y H2 Ω. 53 Similarly, we can obtain anoter error equation c divω ω, divτ + A ω ω, τ = σ σ, τ, τ W. Using te interpolation operator Q, equation 53 can also be written as 54 c divη, divτ + A η, τ = σ σ, τ c divζ, divτ A ζ, τ, τ W. Let τ = η = Q ω ω in 54. It is easy to see tat 55 c 2 divq ω ω 2 L 2 Ω + β Q ω ω 2 L 2 Ω 2 σ σ L 2 Ω 2 Q ω ω L 2 Ω 2 +c divω Q ω L2 Ω divq ω ω L2 Ω +α ω Q ω L 2 Ω 2 Q ω ω L 2 Ω 2 C σ σ L 2 Ω 2 + divω Q ω L 2 Ω + ω Q ω L 2 Ω 2 Q ω ω Hdiv;Ω. Following from Lemma 4.4 and te proved result 49, we deduce 56 Q ω ω Hdiv;Ω C σ σ L2 Ω 2 + divω Q ω L2 Ω + ω Q ω L2 Ω 2 C σ τ H2 Ω 2. τ=σ,ω Incorporating te above results 52, 56 wit te well-known estimates for ρ in Lemma 4.3 and ζ in Lemma 4.4, we derive 57 z z Hs Ω C 2 s y v H2 Ω, s =,, and 58 v=y,z ω ω Hdiv;Ω C σ τ=σ,ω τ H2 Ω 2. Tus Lemma 4.7 is proved. Collecting te bounds given by Lemmas , we ave te following main result. Teorem 4.8. Suppose tat y, σ, z, ω, u and y, σ, z, ω, u are te solutions of OCP-OPT and OCP-OPT, respectively. Assume tat all conditions of Lemmas are valid. Ten for r, s =,, we ave 59 u u L 2 Ω U C + r 2 u + 2 y + σ,

12 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 62 6 y y H s Ω + z z H s Ω C + r 2 u + 2 s y + σ, 6 σ σ Hdiv;Ω + ω ω Hdiv;Ω C + r 2 u + 2 y + σ, were te constant C depends on some norms of y, σ, z, ω, u, but is independent of te mes parameters y, σ and u. Proof. It follows from Lemmas tat 62 u u L 2 Ω U C + r 2 u + z z L2 Ω + divω ω L2 Ω C + r 2 u + 2 y + σ. Moreover, it follows from Lemma 4.5, Lemma 4.7 and 62 tat 63 y y H s Ω + z z H s Ω y ỹ Hs Ω + z z Hs Ω + y ỹ Hs Ω + z z Hs Ω C + r 2 u + 2 s y + σ, 64 σ σ Hdiv;Ω + ω ω Hdiv;Ω σ σ Hdiv;Ω + ω ω Hdiv;Ω + σ σ Hdiv;Ω + ω ω Hdiv;Ω C + r 2 u + 2 y + σ. Tus, Teorem 4.8 follows immediately from Numerical experiments In tis section, we conduct some numerical examples to ceck te efficiency of te spitting least-squares mixed finite element sceme 9-2. Let te spatial domain Ω =, 2 and consider te problem: 65 subject to 66 2 y y d 2 + Ω Ω σ σ d 2 + diva y + cy = f + u, in Ω, y =, Ω u 2 on Ω, wit σ = A y. In fact, te adjoint state z satisfies te following equation: diva z + cz = y yd, in Ω, 67 z =, on Ω. For simplicity we sall adopt te same mes partitions for te states and control in te numerical tests, i.e., Ty = Tσ = Tu. To solve te optimal control problem numerically, we use te C++ software package: AFEpack, it is available at ttp://dsec.pku.edu.cn/ rli/. As for constrained optimal control problems, people pay more attention on te states and te control. Terefore, in te following numerical examples, we mostly center on te primal state variable y, wic is approximated by piecewise linear C elements; te flux state variable σ, wic is

13 622 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU approximated by piecewise linear C elements or te lowest-order RT elements; and te control variable u, wic is discretized using piecewise constant elements or piecewise linear discontinuous elements. Furtermore, in te following tests, we use E u, E y, E y, and E σ to represent te error norms u u L 2 Ω, y y L 2 Ω, y y H Ω, and σ σ Hdiv;Ω, respectively. Example. For te first example, te corresponding analytical solutions for problems wit a constant diffusion coefficient matrix A = I and a constant reaction coefficient c = are given by yx = sinπx sinπx 2, σx = A y, zx = 2 sinπx sinπx 2, ωx = sinπx 2, sinπx, ux = max, min., z + c divω. Te source function f, te desired states y d and σ d can be determined using te above functions. In tis test, we consider te lowest-order RT elements for te approximation of te flux state σ, and piecewise constant elements for tat of te control u. Table displays te errors and convergence orders for te control u, te primal state y, and te flux state σ in corresponding norms. We can see tat te convergence orders for te control u in L 2 Ω-norm, te state y in H Ω-norm, and σ in Hdiv; Ω-norm are almost first-order, wile te convergence order for te state y in L 2 Ω-norm is second-order, wic are consistent wit Teorem 4.8 very well. Figure sows te approximate profiles of te control u and te primal state y. Te two components of te numerical flux state σ are presented in Figure 2. It can be observed tat te numerical figures are agreement wit te analytical solutions very well. Table. σ: te lowest-order RT elements, u: piecewise constant elements. E u Order E y Order E y Order E σ Order E E E E E E E E E E E E E-3..92E E E-. Example 2. For te second example, te corresponding analytical solutions for problems are yx = x x x 2 x 2, σx = A y, zx = 6x x x 2 x 2, ωx = x 2 x 2, x x, ux = max, min.5, z + c divω, were A = I and c =. Te source function f and te desired states y d and σ d can also be determined using te above functions.

14 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP y_ u_ Figure. Te approximate control u left and te approximate state y rigt for Example sigma_2, sigma_, Figure 2. Te approximate flux state σ for Example. In tis example, we consider piecewise linear C elements for te approximation of te flux state σ, and piecewise constant elements for tat of te control u. We list te corresponding numerical errors and convergence orders for te control u, te primal state y and te flux state σ in Table 2. Te convergence orders for te control, te primal state, and te flux state are also well matced wit te teoretical analysis in Teorem 4.8. We also sow te approximate profiles of te control u, te primal state y and te flux state σ in Figures 3-4, respectively. It can be observed tat te numerical results are also agreement wit te analytical solutions very well. Table 2. σ: piecewise linear C elements, u: piecewise constant elements. Eu Order Ey Order Ey Order Eσ E E E-2 6.9E E E E E E E E E E-3..35E E E-3.8 Order Example 3. For te tird example, we consider a reaction-dominated diffusion problem. In tis case, te solution typically as exponential boundary layers on all sides of Ω or internal layers in Ω. We take te analytical solutions for problems

15 624 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU y_ u_ Figure 3. Te approximate control u left and te approximate state y rigt for Example sigma_2, sigma_, Figure 4. Te approximate flux state σ for Example as follows: yx = exp x.52 + x2.52 /ε, σx = ε y, zx = 2 exp x.52 + x2.52 /ε, ωx =,, ux = max, min.5, z + c divω, were ε = 2 and c =. Te source function f and te desired states yd and σd can also be determined using te above functions. In tis example, we consider te lowest-order RT elements for te approximation of te flux state σ, and piecewise linear discontinuous elements for tat of te control u. Some corresponding numerical errors and convergence orders for te control u, te primal state y and te flux state σ are presented in Table 3. It 3 can be observed tat an nearly O 2 -order convergence for te control is derived, wic is consistent wit Teorem 4.8 for te case r =. Te approximate profiles of te control u, te primal state y and te flux state σ are displayed in Figures 5-6, respectively. All tese results also confirms te teoretical analysis. 6. Concluding remarks In tis paper, we derive a priori error estimates for te splitting least-squares mixed finite element discretization of optimal control problem governed by elliptic

16 SPLITTING LEAST-SQUARES MFE FOR ELLIPTIC OCP 625 Table 3. σ: te lowest-order RT elements; u: piecewise linear discontinuous elements. Eu Order Ey Order Ey Order Eσ E E E-.9923E E E E E E E E E E E E E-2. Order y_ u_ Figure 5. Te approximate control u left and te approximate state y rigt for Example sigma_2 sigma_ Figure 6. Te approximate flux state σ for Example 3. equations, were pointwise inequality constraints on te control variable are considered. Tis splitting mixed metod leads to uncoupled positive definite systems for te primal state variable and its flux, also for te adjoint state and flux state variables. Finally, numerical experiments are addressed to verify te teoretical analysis. Acknowledgments Te autors would like to tank te editor and te anonymous referee for teir valuable comments and suggestions on an earlier version of tis paper. Te first autor was partially supported by te Cina Scolarsip Council. References [] P. G. Ciarlet, Te Finite Element Metod for Elliptic Problems, SIAM, Piladelpia, 22. [2] J. L. Lions, Optimal Control of Systems Governed by Partial Differential Equations, SpringerVerlag, Berlin, 97.

17 626 H. FU, H. RUI, H. GUO, J. ZHANG, AND J. HOU [3] W. Liu and N. an, Adaptive Finite Element Metods for Optimal Control Governed by PDEs, Series in Information and Computational Science 4, Science Press, Beijing, 28. [4] M. Hinze, R. Pinnau, M. Ulbric, and S. Ulbric, Optimization wit PDE Constraints, Matematical Modelling: Teory and Applications 23, Springer, Neterlands, 29. [5] D. Tiba, Lectures on te Optimal Control of Elliptic Equations, University of Jyvaskyla Press, Finland, 995. [6] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Metods, Springer-Verlag, 99. [7] J. Douglas Jr. and J. E. Roberts, Global estimates for mixed metods for second order elliptic equations, Mat. Comput., [8] P. A. Raviart and J. M. Tomas, A mixed finite element metod for 2nd order elliptic problems, Lecture Notes in Matematics, vol. 66, Springer, Berlin, 977, [9] F. Brezzi, J. Douglas Jr., M. Fortin and L. D. Marini, Efficient rectangular mixed finite elements in two and tree space variables, RAIRO Model. Mat. Numer. Anal., [] F. Brezzi, J. Douglas Jr. and L. D. Marini, Two family of mixed finite elements for second order elliptic problems, Numer. Mat., []. Cen and W. Liu, Error estimates and superconvergence of mixed finite elements for quadratic optimal control, Int. J. Numer. Anal. Model., [2]. Cen and Z. Lu, Error estimates of fully discrete mixed finite element metods for semilinear quadratic parabolic optimal control problem, Comput. Metods Appl. Mec. Engrg., [3] J. Zou,. Cen and. Dai, Superconvergence of triangular mixed finite elements for optimal control problems wit an integral constraint, Appl. Mat. Comput., [4] A. Pelivanov, G. F. Carey and D. Lazarov, Least-squares mixed finite elements for secondorder elliptic problems, SIAM J. Numer. Anal., [5] H. Rui, S. Kim and S. D. Kim, A remark on least-squares mixed element metods for reactiondiffusion problems, J. Comput. Appl. Mat., [6] H. Rui, S. D. Kim and S. Kim, Split least-squares finite element metods for linear and nonlinear parabolic problems, J. Comput. Appl. Mat., [7] F. Gao and H. Rui, A split least-squares caracteristic mixed finite element metod for Sobolev equations wit convection term, Mat. Comput. Simulation, [8] H. Guo, A remark on split least-squares mixed element procedures for pseudo-parabolic equations, Appl. Mat. Comput., [9] J. Zang and H. Guo, A split least-squares caracteristic mixed element metod for nonlinear nonstationary convection-diffusion problem, Int. J. Comput. Mat., [2] M. F. Weeler, A priori L 2 error estimates for Galerkin approximations to parabolic partial differential equations, SIAM J. Numer. Anal., [2] R. Becker and B. Vexler, Optimal control of te convection-diffusion equation using stabilized finite element metods, Numer. Mat., [22] C. Meyer and A. Rösc, Superconvergence properties of optimal control problems, SIAM J. Control Optim., College of Science, Cina University of Petroleum, Qingdao 26658, Cina ongfeifu@upc.edu.cn Scool of Matematics, Sandong University, Jinan 25, Cina xrui@sdu.edu.cn College of Science, Cina University of Petroleum, Qingdao 26658, Cina sdug@63.com College of Science, Cina University of Petroleum, Qingdao 26658, Cina jszang@upc.edu.cn College of Petroleum Engineering, Cina University of Petroleum, Qingdao 26658, Cina oujian@upc.edu.cn

1. Introduction. We consider the model problem: seeking an unknown function u satisfying

1. Introduction. We consider the model problem: seeking an unknown function u satisfying A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced

More information

Preconditioning in H(div) and Applications

Preconditioning in H(div) and Applications 1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition

More information

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for

More information

AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS

AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous

More information

Parallel algorithm for convection diffusion system based on least squares procedure

Parallel algorithm for convection diffusion system based on least squares procedure Zang et al. SpringerPlus (26 5:69 DOI.86/s464-6-3333-8 RESEARCH Parallel algoritm for convection diffusion system based on least squares procedure Open Access Jiansong Zang *, Hui Guo 2, Hongfei Fu 3 and

More information

Overlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions

Overlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related

More information

arxiv: v1 [math.na] 17 Jul 2014

arxiv: v1 [math.na] 17 Jul 2014 Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order

More information

LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS

LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO

More information

arxiv: v1 [math.na] 20 Jul 2009

arxiv: v1 [math.na] 20 Jul 2009 STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed

More information

A Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems

A Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

A Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems

A Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:

More information

Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods

Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business

More information

A Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems

A Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems

More information

Different Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method

Different Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles

More information

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of

More information

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

Priority Program 1253

Priority Program 1253 Deutsce Forscungsgemeinscaft Priority Program 53 Optimization wit Partial Differential Equations K. Deckelnick, A. Günter and M. Hinze Finite element approximation of elliptic control problems wit constraints

More information

Computers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations

Computers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares

More information

Dedicated to the 70th birthday of Professor Lin Qun

Dedicated to the 70th birthday of Professor Lin Qun Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,

More information

ch (for some fixed positive number c) reaching c

ch (for some fixed positive number c) reaching c GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan

More information

Analysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem

Analysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for

More information

Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS

Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating

More information

arxiv: v1 [math.na] 27 Jan 2014

arxiv: v1 [math.na] 27 Jan 2014 L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates

More information

AN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS

AN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES

More information

b i (x) u + c(x)u = f in Ω,

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

A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION

A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable

More information

A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS

A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a

More information

1. Introduction. Consider a semilinear parabolic equation in the form

1. Introduction. Consider a semilinear parabolic equation in the form A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic

More information

Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations

Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY

More information

arxiv: v1 [math.na] 9 Sep 2015

arxiv: v1 [math.na] 9 Sep 2015 arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on

More information

The Laplace equation, cylindrically or spherically symmetric case

The Laplace equation, cylindrically or spherically symmetric case Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

Adaptive methods for control problems with finite-dimensional control space

Adaptive methods for control problems with finite-dimensional control space Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy

More information

CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES

CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL SINGULARLY PERTURBED DIFFERENTIAL EQUATIONS ON EQUIDISTRIBUTING MESHES INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume, Number, Pages 57 74 c 5 Institute for Scientific Computing and Information CONVERGENCE ANALYSIS OF FINITE ELEMENT SOLUTION OF ONE-DIMENSIONAL

More information

MANY scientific and engineering problems can be

MANY scientific and engineering problems can be A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial

More information

Subdifferentials of convex functions

Subdifferentials of convex functions Subdifferentials of convex functions Jordan Bell jordan.bell@gmail.com Department of Matematics, University of Toronto April 21, 2014 Wenever we speak about a vector space in tis note we mean a vector

More information

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0

Math 212-Lecture 9. For a single-variable function z = f(x), the derivative is f (x) = lim h 0 3.4: Partial Derivatives Definition Mat 22-Lecture 9 For a single-variable function z = f(x), te derivative is f (x) = lim 0 f(x+) f(x). For a function z = f(x, y) of two variables, to define te derivatives,

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

arxiv: v1 [math.na] 12 Mar 2018

arxiv: v1 [math.na] 12 Mar 2018 ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure

More information

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms

More on generalized inverses of partitioned matrices with Banachiewicz-Schur forms More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,

More information

Department of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801

Department of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801 RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused

More information

Effect of Numerical Integration on Meshless Methods

Effect of Numerical Integration on Meshless Methods Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Key words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness

Key words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent

More information

MAT Calculus for Engineers I EXAM #1

MAT Calculus for Engineers I EXAM #1 MAT 65 - Calculus for Engineers I EXAM # Instructor: Liu, Hao Honor Statement By signing below you conrm tat you ave neiter given nor received any unautorized assistance on tis eam. Tis includes any use

More information

Research Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation

Research Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation Advances in Numerical Analysis Volume 204, Article ID 35394, 8 pages ttp://dx.doi.org/0.55/204/35394 Researc Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic

More information

Finite Difference Method

Finite Difference Method Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According

More information

Linearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization

Linearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems

More information

Downloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see

Downloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT

More information

c 2004 Society for Industrial and Applied Mathematics

c 2004 Society for Industrial and Applied Mathematics SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares

More information

A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions

A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical

More information

MIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp

MIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a

More information

A SHORT INTRODUCTION TO BANACH LATTICES AND

A SHORT INTRODUCTION TO BANACH LATTICES AND CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,

More information

On convergence of the immersed boundary method for elliptic interface problems

On convergence of the immersed boundary method for elliptic interface problems On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many

More information

INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS WITH IMPLICIT TIME-INTEGRATION TECHNIQUES FOR NONLINEAR PARABOLIC EQUATIONS

INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS WITH IMPLICIT TIME-INTEGRATION TECHNIQUES FOR NONLINEAR PARABOLIC EQUATIONS INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS WITH IMPLICIT TIME-INTEGRATION TECHNIQUES FOR NONLINEAR PARABOLIC EQUATIONS LUNJI SONG,, GUNG-MIN GIE 3, AND MING-CHENG SHIUE 4 Scool of Matematics and Statistics,

More information

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS

ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods

More information

SOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL

SOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL Electronic Journal of Differential Equations, Vol. 7(7, No., pp.. ISSN: 7-669. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu (login: ftp SOLUTIONS OF FOURTH-ORDER PARTIAL

More information

A Finite Element Primer

A Finite Element Primer A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.

More information

Finite Element Methods for Linear Elasticity

Finite Element Methods for Linear Elasticity Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu

More information

ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS

ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis

More information

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014 Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.

More information

L -Estimates of Rectangular Mixed Methods for Nonlinear Constrained Optimal Control Problem

L -Estimates of Rectangular Mixed Methods for Nonlinear Constrained Optimal Control Problem BULLETIN of the MALAYSIAN MATHEMATICAL SCIENCES SOCIETY http:/math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 37(1) (2014), 271 284 L -Estimates of Rectangular Mixed Methods for Nonlinear Constrained

More information

Continuity. Example 1

Continuity. Example 1 Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)

More information

arxiv: v1 [math.na] 20 Nov 2018

arxiv: v1 [math.na] 20 Nov 2018 An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018

More information

ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS. 1. Introduction. We consider the parabolic PDE of the form,

ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS. 1. Introduction. We consider the parabolic PDE of the form, ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS K. CHRYSAFINOS AND NOEL. J. WALKINGTON Abstract. We analyze te classical discontinuous Galerkin metod for a general parabolic

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

Recent Advances in Time-Domain Maxwell s Equations in Metamaterials

Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Recent Advances in Time-Domain Maxwell s Equations in Metamaterials Yunqing Huang 1, and Jicun Li, 1 Hunan Key Laboratory for Computation and Simulation in Science and Engineering, Xiangtan University,

More information

Weak Solutions for Nonlinear Parabolic Equations with Variable Exponents

Weak Solutions for Nonlinear Parabolic Equations with Variable Exponents Communications in Matematics 25 217) 55 7 Copyrigt c 217 Te University of Ostrava 55 Weak Solutions for Nonlinear Parabolic Equations wit Variable Exponents Lingeswaran Sangerganes, Arumugam Gurusamy,

More information

Variational Localizations of the Dual Weighted Residual Estimator

Variational Localizations of the Dual Weighted Residual Estimator Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)

More information

arxiv: v1 [math.na] 28 Apr 2017

arxiv: v1 [math.na] 28 Apr 2017 THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular

More information

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER* EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

arxiv: v1 [math.na] 19 Mar 2018

arxiv: v1 [math.na] 19 Mar 2018 A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne

More information

A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION

A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous

More information

Combining functions: algebraic methods

Combining functions: algebraic methods Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)

More information

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006

Math 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2

More information

arxiv: v1 [math.na] 18 Jul 2015

arxiv: v1 [math.na] 18 Jul 2015 ENERGY NORM ERROR ESTIMATES FOR FINITE ELEMENT DISCRETIZATION OF PARABOLIC PROBLEMS HERBERT EGGER arxiv:150705183v1 [matna] 18 Jul 2015 Department of Matematics, TU Darmstadt, Germany Abstract We consider

More information

Technische Universität Graz

Technische Universität Graz Tecnisce Universität Graz Boundary element metods for Diriclet boundary control problems Günter Of, Tan Pan Xuan, Olaf Steinbac Bericte aus dem Institut für Numerisce Matematik Berict 2009/2 Tecnisce

More information

Priority Program 1253

Priority Program 1253 Deutsce Forscungsgemeinscaft Priority Program 1253 Optimization wit Partial Differential Equations K. Deckelnick and M. Hinze Convergence and error analysis of a numerical metod for te identification of

More information

DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS

DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS Journal o Applied Analysis Vol. 14, No. 2 2008, pp. 259 271 DIFFERENTIAL POLYNOMIALS GENERATED BY SECOND ORDER LINEAR DIFFERENTIAL EQUATIONS B. BELAÏDI and A. EL FARISSI Received December 5, 2007 and,

More information

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations

Parameter Fitted Scheme for Singularly Perturbed Delay Differential Equations International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department

More information

University Mathematics 2

University Mathematics 2 University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at

More information

CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS

CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7, Number 1, Pages 1 9 c 010 Institute for Scientific Computing and Information CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION

More information

On convexity of polynomial paths and generalized majorizations

On convexity of polynomial paths and generalized majorizations On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal

More information

arxiv: v1 [math.na] 29 Jan 2019

arxiv: v1 [math.na] 29 Jan 2019 Iterated total variation regularization wit finite element metods for reconstruction te source term in elliptic systems Micael Hinze a Tran Nan Tam Quyen b a Department of Matematics, University of Hamburg,

More information

An approximation method using approximate approximations

An approximation method using approximate approximations Applicable Analysis: An International Journal Vol. 00, No. 00, September 2005, 1 13 An approximation metod using approximate approximations FRANK MÜLLER and WERNER VARNHORN, University of Kassel, Germany,

More information

Smoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations

Smoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations 396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt

More information

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means

More information

CONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA

CONVERGENCE ANALYSIS OF YEE SCHEMES FOR MAXWELL S EQUATIONS IN DEBYE AND LORENTZ DISPERSIVE MEDIA INTRNATIONAL JOURNAL OF NUMRICAL ANALYSIS AND MODLING Volume XX Number 0 ages 45 c 03 Institute for Scientific Computing and Information CONVRGNC ANALYSIS OF Y SCHMS FOR MAXWLL S QUATIONS IN DBY AND LORNTZ

More information

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,

NUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example, NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing

More information

EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD

EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1379 1396 S 0025-5718(99)01093-5 Article electronically publised on February 24, 1999 EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD

More information

H(div) conforming and DG methods for incompressible Euler s equations

H(div) conforming and DG methods for incompressible Euler s equations H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible

More information

FINITE ELEMENT ANALYSIS OF OPTIMAL CONTROL PROBLEM GOVERNED BY STOKES EQUATIONS WITH L 2 -NORM STATE-CONSTRAINTS * 1. Introduction

FINITE ELEMENT ANALYSIS OF OPTIMAL CONTROL PROBLEM GOVERNED BY STOKES EQUATIONS WITH L 2 -NORM STATE-CONSTRAINTS * 1. Introduction Journal of Computational Matematics Vol.29, No.5, 2011, 589 604. ttp://www.global-sci.org/jcm doi:10.4208/jcm.1103-m3514 FINITE ELEMENT ANALYSIS OF OPTIMAL CONTROL PROBLEM GOVERNED BY STOKES EQUATIONS

More information

Order of Accuracy. ũ h u Ch p, (1)

Order of Accuracy. ũ h u Ch p, (1) Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical

More information

Pre-Calculus Review Preemptive Strike

Pre-Calculus Review Preemptive Strike Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

arxiv: v1 [math.dg] 4 Feb 2015

arxiv: v1 [math.dg] 4 Feb 2015 CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola

More information