Partial Differential Equations - Modelling and ... - ResearchGate
Partial Differential Equations - Modelling and ... - ResearchGate Partial Differential Equations - Modelling and ... - ResearchGate
240 P. Neittaanmäki and D. Tiba The discussion in Section 2 shows the existence of infinitely many admissible pairs [u, y] for the constrained control problem (15)–(18). (Here g is fixed satisfying the necessary smoothness properties.) In case g and Ω g ⊂ D are variable and unknown, we say that (15)–(18) is the variational fixed domain (in D!) formulation of the Neumann boundary value problem. One can write the optimality conditions that give a system of equations equivalent with (15)–(18) and extend the Neumann problem from Ω g to D. We introduce the penalized control problem, for ε>0, as follows (here [g ≡ 0] denotes ∂Ω g ): subject to Above, − d∑ i,j=1 Min u∈L 2 (D) { ∫ 1 u 2 dx + 1 ∫ } F (y ε ) 2 dσ 2 D 2ε [g≡0] (19) ( ) ∂ ∂y ε a ij + a 0 y ε = f +(1− H(g))u in D, (20) ∂x j ∂x i F (y) = ∂y ε ∂n A =0 on∂D. (21) d∑ i,j=1 a ij ∂y ∂x j ∇g · e i and the problem (19)–(21), which is unconstrained, remains a coercive and strictly convex control problem. That is, we have the existence and the uniqueness of the approximating optimal pair [u ε ,y ε ] ∈ L 2 (D) × H 2 (D) (if∂D is smooth enough). Proposition 1. We have |F (y ε )| L2 (∂Ω g) ≤ Cε 1 2 , (22) u ε → û strongly in L 2 (D), (23) y ε → ŷ strongly in H 2 (D), (24) where C is a constant independent of ε>0 and [û, ŷ] ∈ L 2 (D) × H 2 (D) is the unique optimal pair of (15)–(18). Proof. As in Section 2, by the trace theorem, we may choose ỹ ∈ H 2 (D \ Ω g ) ∂ỹ with the property that ∂n A =0in∂(D \ Ω g )andỹ may be extended to the solution of (2) inside Ω g . We can compute ũ ∈ L 2 (D \ Ω g ) by (20) and extend it by 0 inside Ω g .Then[ũ, ỹ] is an admissible pair for the control problem (19)–(21) and, by the optimality of [u ε ,y ε ], we get
Shape Optimization Problems with Neumann Boundary Condition 241 ∫ 1 u 2 ε dx + 1 ∫ F (y ε ) 2 dσ ≤ 1 ∫ ũ 2 dx (25) 2 D 2ε [g≡0] 2 D since F (ỹ) =0in∂Ω g . The inequality (25) gives (22) and {u ε } bounded in L 2 (D). By (20), (21), {y ε } is bounded in H 2 (D) and, on a subsequence, we have y ε → ŷ, u ε → û weakly in H 2 (D), respectively in L 2 (D), where [û, ŷ] again satisfy (20), (21). Moreover, one can pass to the limit in (22) with ε → 0, to see that F (ŷ) =0in∂Ω g . This shows that [û, ŷ] is an admissible pair for the original state constrained control problem (15)–(18). For any admissible pair [µ, z] ∈ L 2 (D) × H 2 (D) of (15)–(18), we have F (z) =0on∂Ω g and the inequality (25) is valid with ũ replaced by µ and we infer ∫ 1 u 2 ε dx ≤ 1 ∫ µ 2 dx. 2 D 2 D The weak lower semicontinuity of the norm gives ∫ 1 2 D(û) 2 dx ≤ 1 ∫ µ 2 dx, 2 D that is, the pair [û, ŷ] is, in fact, the unique optimal pair of (15)–(18) and we also have ∫ ∫ lim u 2 ε dx = (û) 2 dx. ε→0 D D Then u ε → û strongly in L 2 (D) andy ε → ŷ strongly in H 2 (D) by the strong convergence criterion in uniformly convex spaces. The convergence is valid without taking subsequences due to the uniqueness of [û, ŷ]. Remark 2. One can further regularize H in (20), by replacing it with a mollification H ε of the Yosida approximation H ε of the maximal monotone extension of H. Remark 3. One may take in D even null Dirichlet boundary conditions instead of (16). Similar distributed controllability properties (approximate or exact) maybeestablishedinverymuchthesameway. To write shortly, we consider the case of the Laplace operator. The penalized and regularized problem is the following: Min u∈L 2 (D) { ∫ 1 u 2 dx + 1 ∫ [∇y ·∇g] 2 dσ 2 D 2ε [g≡0] −∆y + y = f +(1− H ε (g))u in D, y =0 on∂D. Here, the control u ensures the “transfer” from Dirichlet to Neumann (null) conditions on ∂Ω g and all the results are similar as for the Neumann–Neumann case. } ,
- Page 187 and 188: Numerical Analysis of a Finite Elem
- Page 189 and 190: 188 A. Bonito et al. of the model a
- Page 191 and 192: 190 A. Bonito et al. Fig. 1. The sp
- Page 193 and 194: 192 A. Bonito et al. 3 1 16 4 1 1 1
- Page 195 and 196: 194 A. Bonito et al. ∫ v n+1 h
- Page 197 and 198: 196 A. Bonito et al. −pn +2µD(v)
- Page 199 and 200: 198 A. Bonito et al. each of its pa
- Page 201 and 202: 200 A. Bonito et al. The normal vec
- Page 203 and 204: 202 A. Bonito et al. with initial c
- Page 205 and 206: 204 A. Bonito et al. Fig. 8. Jet bu
- Page 207 and 208: 206 A. Bonito et al. References [AM
- Page 209 and 210: 208 A. Bonito et al. [Set96] J. A.
- Page 211 and 212: 210 J. Hao et al. due to shear flow
- Page 213 and 214: 212 J. Hao et al. The backward reac
- Page 215 and 216: 214 J. Hao et al. where u and p den
- Page 217 and 218: 216 J. Hao et al. * * * * * * * * *
- Page 219 and 220: 218 J. Hao et al. and solve for V n
- Page 221 and 222: 220 J. Hao et al. Table 2. The calc
- Page 223 and 224: 222 J. Hao et al. References [ASS80
- Page 225 and 226: Computing the Eigenvalues of the La
- Page 227 and 228: Eigenvalues of the Laplace-Beltrami
- Page 229 and 230: Eigenvalues of the Laplace-Beltrami
- Page 231 and 232: Eigenvalues of the Laplace-Beltrami
- Page 233 and 234: A Fixed Domain Approach in Shape Op
- Page 235 and 236: Shape Optimization Problems with Ne
- Page 237: Shape Optimization Problems with Ne
- Page 241 and 242: Shape Optimization Problems with Ne
- Page 243 and 244: Reduced-Order Modelling of Dispersi
- Page 245 and 246: Reduced-Order Modelling of Dispersi
- Page 247 and 248: Reduced-Order Modelling of Dispersi
- Page 249 and 250: Reduced-Order Modelling of Dispersi
- Page 251 and 252: Reduced-Order Modelling of Dispersi
- Page 253 and 254: Reduced-Order Modelling of Dispersi
- Page 255 and 256: Calibration of Lévy Processes with
- Page 257 and 258: Calibration of Lévy Processes with
- Page 259 and 260: Calibration of Lévy Processes with
- Page 261 and 262: Calibration of Lévy Processes with
- Page 263 and 264: We have proved Calibration of Lévy
- Page 265 and 266: Calibration of Lévy Processes with
- Page 267 and 268: Calibration of Lévy Processes with
- Page 269 and 270: Calibration of Lévy Processes with
- Page 271 and 272: Note that p ∗ satisfies Calibrati
- Page 273 and 274: Calibration of Lévy Processes with
- Page 275 and 276: 280 S. Ikonen and J. Toivanen the p
- Page 277 and 278: 282 S. Ikonen and J. Toivanen Merto
- Page 279 and 280: 284 S. Ikonen and J. Toivanen For H
- Page 281 and 282: 286 S. Ikonen and J. Toivanen and {
- Page 283 and 284: 288 S. Ikonen and J. Toivanen 1.6 1
- Page 285 and 286: 290 S. Ikonen and J. Toivanen 8 Con
- Page 287: 292 S. Ikonen and J. Toivanen [Mer7
Shape Optimization Problems with Neumann Boundary Condition 241<br />
∫<br />
1<br />
u 2 ε dx + 1 ∫<br />
F (y ε ) 2 dσ ≤ 1 ∫<br />
ũ 2 dx (25)<br />
2 D 2ε [g≡0] 2 D<br />
since F (ỹ) =0in∂Ω g .<br />
The inequality (25) gives (22) <strong>and</strong> {u ε } bounded in L 2 (D). By (20), (21),<br />
{y ε } is bounded in H 2 (D) <strong>and</strong>, on a subsequence, we have y ε → ŷ, u ε →<br />
û weakly in H 2 (D), respectively in L 2 (D), where [û, ŷ] again satisfy (20),<br />
(21). Moreover, one can pass to the limit in (22) with ε → 0, to see that<br />
F (ŷ) =0in∂Ω g . This shows that [û, ŷ] is an admissible pair for the original<br />
state constrained control problem (15)–(18). For any admissible pair [µ, z] ∈<br />
L 2 (D) × H 2 (D) of (15)–(18), we have F (z) =0on∂Ω g <strong>and</strong> the inequality<br />
(25) is valid with ũ replaced by µ <strong>and</strong> we infer<br />
∫<br />
1<br />
u 2 ε dx ≤ 1 ∫<br />
µ 2 dx.<br />
2 D 2 D<br />
The weak lower semicontinuity of the norm gives<br />
∫<br />
1<br />
2 D(û) 2 dx ≤ 1 ∫<br />
µ 2 dx,<br />
2 D<br />
that is, the pair [û, ŷ] is, in fact, the unique optimal pair of (15)–(18) <strong>and</strong> we<br />
also have<br />
∫ ∫<br />
lim u 2 ε dx = (û) 2 dx.<br />
ε→0<br />
D<br />
D<br />
Then u ε → û strongly in L 2 (D) <strong>and</strong>y ε → ŷ strongly in H 2 (D) by the strong<br />
convergence criterion in uniformly convex spaces. The convergence is valid<br />
without taking subsequences due to the uniqueness of [û, ŷ].<br />
Remark 2. One can further regularize H in (20), by replacing it with a mollification<br />
H ε of the Yosida approximation H ε of the maximal monotone extension<br />
of H.<br />
Remark 3. One may take in D even null Dirichlet boundary conditions instead<br />
of (16). Similar distributed controllability properties (approximate or exact)<br />
maybeestablishedinverymuchthesameway.<br />
To write shortly, we consider the case of the Laplace operator. The penalized<br />
<strong>and</strong> regularized problem is the following:<br />
Min<br />
u∈L 2 (D)<br />
{ ∫<br />
1<br />
u 2 dx + 1 ∫<br />
[∇y ·∇g] 2 dσ<br />
2 D 2ε [g≡0]<br />
−∆y + y = f +(1− H ε (g))u in D,<br />
y =0<br />
on∂D.<br />
Here, the control u ensures the “transfer” from Dirichlet to Neumann (null)<br />
conditions on ∂Ω g <strong>and</strong> all the results are similar as for the Neumann–Neumann<br />
case.<br />
}<br />
,