Partial Differential Equations - Modelling and ... - ResearchGate

More documents

Recommendations

Info

228 R. Glowinski and D.C. Sorensen If these discrete boundary conditions are used to eliminate the unknowns u I+1j ,u 0j ,u iI+1 and u i0 , we obtain the following discrete eigenproblem (in R N , N = I 2 ): If 2 ≤ i, j ≤ I − 1, 2[(Rρ −1 + cos θ j ) −1 + Rρ −1 + cos θ j cos(h/2)]u ij − (Rρ −1 + cos θ j ) −1 (u i+1j + u i−1j ) − (Rρ −1 + cos(θ j + h/2))u ij+1 − (Rρ −1 + cos(θ j − h/2))u ij−1 = λρ 2 (Rρ −1 + cos θ j )h 2 u ij . (5) If i =1and2≤ j ≤ I − 1, 2[(Rρ −1 + cos θ j ) −1 + Rρ −1 + cos θ j cos(h/2)]u 1j − (Rρ −1 + cos θ j ) −1 (u 2j + u Ij ) − (Rρ −1 + cos(θ j + h/2))u 1j+1 − (Rρ −1 + cos(θ j − h/2))u 1j−1 = λρ 2 (Rρ −1 + cos θ j )h 2 u 1j . (6) If i = j =1, 2[(Rρ −1 + cos h) −1 + Rρ −1 + cos h cos(h/2)]u 11 − (Rρ −1 + cos h) −1 (u 21 + u I1 ) − (Rρ −1 +cos(3h/2))u 12 − (Rρ −1 + cos(h/2))u 1I = λρ 2 (Rρ −1 + cos h)h 2 u 11 . (7) If i =1andj = I, 2[(Rρ −1 +1) −1 + Rρ −1 + cos(h/2)]u 1I − (Rρ −1 +1) −1 (u 2I + u II ) − (Rρ −1 + cos(h/2))u 11 − (Rρ −1 + cos(h/2))u 1I−1 = λρ 2 (Rρ −1 +1)h 2 u 1I . (8) If i = I and 2 ≤ j ≤ I − 1, 2[(Rρ −1 + cos θ j ) −1 + Rρ −1 + cos θ j cos(h/2)]u Ij − (Rρ −1 + cos θ j ) −1 (u 1j + u I−1j ) − (Rρ −1 + cos(θ j + h/2))u Ij+1 − (Rρ −1 + cos(θ j − h/2))u Ij−1 = λρ 2 (Rρ −1 + cos θ j )h 2 u Ij . (9) If i = I and j =1, 2[(Rρ −1 + cos h) −1 + Rρ −1 + cos h cos(h/2)]u I1 − (Rρ −1 + cos h) −1 (u 11 + u I−11 ) − (Rρ −1 +cos(3h/2))u I2 − (Rρ −1 + cos(h/2))u II = λρ 2 (Rρ −1 + cos h)h 2 u I1 . (10)
Eigenvalues of the Laplace–Beltrami Operator on the Surface of a Torus 229 If i = I and j = I, 2[(Rρ −1 +1) −1 + Rρ −1 + cos(h/2)]u II − (Rρ −1 +1) −1 (u 1I + u I−1I ) − (Rρ −1 + cos(h/2))u I1 − (Rρ −1 + cos(h/2))u II−1 = λρ 2 (Rρ −1 +1)h 2 u II . (11) If 2 ≤ i ≤ I − 1andj =1, 2[(Rρ −1 + cos(h)) −1 + Rρ −1 + cos(h) cos(h/2)]u i1 − (Rρ −1 + cos(h)) −1 (u i+11 + u i−11 ) − (Rρ −1 +cos(3h/2))u i2 − (Rρ −1 + cos(h/2))u iI = λρ 2 (Rρ −1 + cos(h))h 2 u i1 . (12) If 2 ≤ i ≤ I − 1andj = I, 2[(Rρ −1 +1) −1 + Rρ −1 + cos(h/2)]u iI − (Rρ −1 +1) −1 (u i+1I + u i−1I ) − (Rρ −1 + cos(h/2))u i1 − (Rρ −1 + cos(h/2))u iI−1 = λρ 2 (Rρ −1 +1)h 2 u iI . (13) If 2 ≤ i, j ≤ I − 1, 2[(Rρ −1 + cos θ j ) −1 + Rρ −1 + cos θ j cos(h/2)]u ij − (Rρ −1 + cos θ j ) −1 (u i+1j + u i−1j ) − (Rρ −1 + cos(θ j + h/2))u ij+1 − (Rρ −1 + cos(θ j − h/2))u ij−1 = λρ 2 (Rρ −1 + cos θ j )h 2 u ij . (14) These finite difference formulas generate an approximation to the problem (3) in the form of a symmetric generalized eigenvalue problem Ax = λDx, (15) with A sparse and symmetric positive semi-definite and with D positive definite and diagonal (independent of the ordering of the variables). We used Matlab to solve the problem (15) to obtain approximations to eigenvalues and corresponding eigenfunctions of (2). 5 Numerical Experiments The Matlab function eigs which is based upon ARPACK [Sor92, LSY98] was used to perform the numerical calculation of eigenvalues and corresponding eigenvectors. In all cases, we computed the 15 lowest (algebraically smallest) eigenvalues of the generalized eigenvalue problem (15) using the shift-invert option with shift σ = −.0001. Since the eigenvalues are real and non-negative, the eigenvalues closest to the origin are enhanced with this transformation and thus easily computed with a Krylov method.
Page 1 and 2:
Partial Differential Equations
Page 3 and 4:
Partial Differential Equations Mode
Page 5 and 6:
Dedicated to Olivier Pironneau
Page 7 and 8:
VIII Preface computers has been at
Page 9 and 10:
Contents List of Contributors .....
Page 11 and 12:
List of Contributors Yves Achdou UF
Page 13 and 14:
List of Contributors XV Claude Le B
Page 15 and 16:
Discontinuous Galerkin Methods Vive
Page 17 and 18:
Discontinuous Galerkin Methods 5 2
Page 19 and 20:
Discontinuous Galerkin Methods 7 g
Page 21 and 22:
Discontinuous Galerkin Methods 9 (
Page 23 and 24:
Page 25 and 26:
Page 27 and 28:
Discontinuous Galerkin Methods 15 W
Page 29 and 30:
Let a h and b h denote the bilinear
Page 31 and 32:
Discontinuous Galerkin Methods 19 t
Page 33 and 34:
Discontinuous Galerkin Methods 21 l
Page 35 and 36:
Table 1. Primal DG for transport Di
Page 37 and 38:
Discontinuous Galerkin Methods 25 [
Page 39 and 40:
Mixed Finite Element Methods on Pol
Page 41 and 42:
Mixed FE Methods on Polyhedral Mesh
Page 43 and 44:
Page 45 and 46:
Page 47 and 48:
4 Hybridization and Condensation Mi
Page 49 and 50:
Page 51 and 52:
is symmetric and positive definite,
Page 53 and 54:
with some coefficient α ∈ R wher
Page 55 and 56:
44 E.J. Dean and R. Glowinski so fa
Page 57 and 58:
46 E.J. Dean and R. Glowinski 2 A L
Page 59 and 60:
48 E.J. Dean and R. Glowinski S:T=
Page 61 and 62:
50 E.J. Dean and R. Glowinski minim
Page 63 and 64:
52 E.J. Dean and R. Glowinski 6 On
Page 65 and 66:
54 E.J. Dean and R. Glowinski Fig.
Page 67 and 68:
56 E.J. Dean and R. Glowinski and
Page 69 and 70:
58 E.J. Dean and R. Glowinski 7 Num
Page 71 and 72:
60 E.J. Dean and R. Glowinski Fig.
Page 73 and 74:
62 E.J. Dean and R. Glowinski Assum
Page 75 and 76:
Higher Order Time Stepping for Seco
Page 77 and 78:
u n+1 h Optimal Higher Order Time D
Page 79 and 80:
Optimal Higher Order Time Discretiz
Page 81 and 82:
Page 83 and 84:
Page 85 and 86:
Page 87 and 88:
Page 89 and 90:
Page 91 and 92:
Page 93 and 94:
Page 95 and 96:
Page 97 and 98:
Page 99 and 100:
Page 101 and 102:
Page 103 and 104:
96 I. Sazonov et al. To provide a p
Page 105 and 106:
98 I. Sazonov et al. In the first s
Page 107 and 108:
100 I. Sazonov et al. Fig. 1. An ex
Page 109 and 110:
102 I. Sazonov et al. H z 1 exact F
Page 111 and 112:
104 I. Sazonov et al. (a) (b) Fig.
Page 113 and 114:
106 I. Sazonov et al. Scattering Wi
Page 115 and 116:
108 I. Sazonov et al. 6.4 Scatterin
Page 117 and 118:
110 I. Sazonov et al. (a) (b) Fig.
Page 119 and 120:
112 I. Sazonov et al. [MHP96] K. Mo
Page 121 and 122:
114 R. Sanders and A.M. Tesdall I R
Page 123 and 124:
116 R. Sanders and A.M. Tesdall imp
Page 125 and 126:
118 R. Sanders and A.M. Tesdall alo
Page 127 and 128:
120 R. Sanders and A.M. Tesdall (a)
Page 129 and 130:
122 R. Sanders and A.M. Tesdall D C
Page 131 and 132:
124 R. Sanders and A.M. Tesdall 8.6
Page 133 and 134:
126 R. Sanders and A.M. Tesdall 0.3
Page 135 and 136:
128 R. Sanders and A.M. Tesdall [TR
Page 137 and 138:
132 S. Lapin et al. Ω R γ Ω 2 Γ
Page 139 and 140:
134 S. Lapin et al. ∫ ∂ 2 ∫
Page 141 and 142:
136 S. Lapin et al. Ω R γ Fig. 3.
Page 143 and 144:
138 S. Lapin et al. 4 Energy Inequa
Page 145 and 146:
140 S. Lapin et al. 5 Numerical Exp
Page 147 and 148:
142 S. Lapin et al. Fig. 6. Contour
Page 149 and 150:
144 S. Lapin et al. Fig. 9. Obstacl
Page 151 and 152:
Domain Decomposition and Electronic
Page 153 and 154:
Domain Decomposition Approach for C
Page 155 and 156:
Page 157 and 158:
Page 159 and 160:
Page 161 and 162:
Page 163 and 164:
Page 165 and 166:
Page 167 and 168:
Page 169 and 170:
Numerical Analysis of a Finite Elem
Page 171 and 172:
Page 173 and 174:
Page 175 and 176:
Page 177 and 178: Numerical Analysis of a Finite Elem
Page 181 and 182: so that |u| 2 1,ω h ≤ 2 Numerica
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: 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 243 and 244: Reduced-Order Modelling of Dispersi
Page 255 and 256: Calibration of Lévy Processes with
Page 263 and 264: We have proved Calibration of Lévy
Page 271 and 272: Note that p ∗ satisfies Calibrati
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
show all

Partial Differential Equations - Modelling and ... - ResearchGate

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?