A Fourier Continuation Method for the Solution of Elliptic Eigenvalue Problems in General Domains

<div>Convergence in the first 10 eigenvalues for the unit disk with the variable coefficient <svg height="15.1252pt" id="M267" style="vertical-align:-3.4033pt" version="1.1" viewbox="-0.0498162 -11.7219 111.846 15.1252" width="111.846pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M563 583C563 664 501 710 433 710C381 710 334 692 285 648C212 582 178 506 149 350L74 -56C49 -190 36 -222 24 -226L30 -257C52 -256 95 -247 121 -229C131 -222 134 -206 165 -21L219 303C254 510 296 665 397 665C447 665 481 626 481 572C481 512 446 457 390 420C367 405 351 400 323 397L302 349C398 338 466 285 466 197C466 117 419 45 306 45C268 45 224 59 198 74L186 50C197 20 228 -14 274 -14C318 -14 359 3 401 26C478 68 554 144 554 225C554 322 493 367 412 396C490 446 563 500 563 583Z" id="g50-224"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,7.632,0)"><path d="M309 -142C211 -61 151 99 151 271S210 601 309 683L288 710C141 611 75 451 75 272V271C75 90 141 -72 288 -169L309 -142Z" id="g50-41"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,12.286,0)"><path d="M545 403C545 425 527 451 498 451C448 451 400 404 312 286L290 341C262 412 246 451 219 451C182 451 138 404 92 345L112 323C152 368 169 378 181 378C192 378 201 359 216 321L257 214C184 115 142 69 113 69C102 69 93 73 88 79S78 89 68 89C47 89 24 61 24 40C24 8 49 -12 75 -12C125 -12 175 31 271 175L312 64C330 15 357 -12 382 -12C421 -12 475 35 515 98L496 121C466 88 439 62 420 62C402 62 384 92 363 147L325 247C349 280 374 309 397 333C422 361 446 381 462 381C471 381 481 374 488 366C493 360 499 357 505 357C521 357 545 380 545 403Z" id="g50-121"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="569" vert-adv-y="569"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,19.684,0)"><path d="M98 134C72 134 46 117 46 90C46 73 55 65 60 64C95 55 124 32 124 -4C124 -42 95 -68 44 -89L57 -123C122 -104 194 -60 194 22C194 94 136 134 98 134Z" id="g50-45"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,24.97,0)"><path d="M563 395C563 428 543 451 519 451C500 451 480 438 472 423C468 415 465 403 470 394C477 381 482 368 482 347C482 272 386 136 342 74H340C334 163 323 257 312 337C302 410 290 451 257 451C214 451 160 385 127 338L145 310C169 341 203 373 209 373C219 373 225 369 232 325C249 215 269 62 271 -15C220 -80 135 -165 18 -231L26 -257L142 -230C254 -105 279 -74 342 13C394 85 484 211 526 286C553 334 563 365 563 395Z" id="g50-122"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,32.603,0)"><path d="M283 271C283 451 220 610 70 710L48 683C146 603 207 443 207 271S149 -59 48 -142L70 -169C220 -67 283 90 283 271Z" id="g50-42"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="358" vert-adv-y="358"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,40.932,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,52.551,0)"><path d="M545 403C545 425 527 451 498 451C448 451 400 404 312 286L290 341C262 412 246 451 219 451C182 451 138 404 92 345L112 323C152 368 169 378 181 378C192 378 201 359 216 321L257 214C184 115 142 69 113 69C102 69 93 73 88 79S78 89 68 89C47 89 24 61 24 40C24 8 49 -12 75 -12C125 -12 175 31 271 175L312 64C330 15 357 -12 382 -12C421 -12 475 35 515 98L496 121C466 88 439 62 420 62C402 62 384 92 363 147L325 247C349 280 374 309 397 333C422 361 446 381 462 381C471 381 481 374 488 366C493 360 499 357 505 357C521 357 545 380 545 403Z" id="g50-121"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="569" vert-adv-y="569"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,59.948,-5.872)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,67.834,0)"><path d="M556 236V289H337V504H275V289H56V236H275V-4H337V236H556Z" id="g54-36"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,78.718,0)"><path d="M563 395C563 428 543 451 519 451C500 451 480 438 472 423C468 415 465 403 470 394C477 381 482 368 482 347C482 272 386 136 342 74H340C334 163 323 257 312 337C302 410 290 451 257 451C214 451 160 385 127 338L145 310C169 341 203 373 209 373C219 373 225 369 232 325C249 215 269 62 271 -15C220 -80 135 -165 18 -231L26 -257L142 -230C254 -105 279 -74 342 13C394 85 484 211 526 286C553 334 563 365 563 395Z" id="g50-122"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,86.351,-5.872)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,94.238,0)"><path d="M556 236V289H337V504H275V289H56V236H275V-4H337V236H556Z" id="g54-36"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,105.121,0)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z" id="g50-50"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg> is shown with respect to the finest discretization; <svg height="8.80314pt" id="M268" style="vertical-align:-0.2063999pt" version="1.1" viewbox="-0.0498162 -8.59674 64.6691 8.80314" width="64.6691pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M870 650H622L613 619C694 615 708 606 708 559C708 531 702 477 686 396L636 136H633L281 650H127L120 619C180 616 198 611 212 592C227 572 232 558 220 503L166 254C148 169 139 131 126 89C112 43 85 35 33 32L24 0H280L286 32C208 37 189 44 189 93C189 124 198 179 212 250L267 526H271L622 -8H655L732 396C747 477 762 532 770 556C785 600 798 616 862 619L870 650Z" id="g50-79"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="888" vert-adv-y="888"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,14.754,0)"><path d="M556 331V384H55V331H556ZM556 140V193H55V140H556Z" id="g54-34"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,26.373,0)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,32.704,0)"><path d="M158 548H390L417 615L410 623H122L83 318C105 326 143 337 185 337C296 337 350 275 350 188C350 116 308 42 225 42C164 42 122 74 100 93C90 101 82 99 72 92C60 82 51 68 50 59C48 46 52 38 66 24C82 9 125 -12 172 -12C225 -11 292 15 346 59C408 108 437 166 437 226C437 309 371 397 242 397C214 397 170 382 133 369L158 548Z" id="g50-54"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,39.035,0)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z" id="g50-51"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,45.366,0)"><path d="M253 635C142 635 71 555 71 471C71 401 113 353 178 316C144 294 106 267 91 252C70 231 46 202 46 157C46 51 132 -12 241 -12C328 -12 442 52 442 168C442 257 379 305 311 343C354 373 380 398 388 407C406 428 417 458 417 487C417 568 348 635 253 635ZM240 600C288 600 340 563 340 481C340 423 313 387 279 360C207 394 148 427 148 499C148 549 182 600 240 600ZM252 23C188 23 129 73 129 163C129 217 162 266 211 298C289 260 360 217 360 144C360 69 312 23 252 23Z" id="g50-57"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,51.697,0)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,58.028,0)"><path d="M245 635C92 635 37 457 37 312C37 149 91 -12 244 -12C395 -12 449 166 449 312C449 469 395 635 245 635ZM243 598C332 598 358 454 358 312C358 173 334 26 245 26C158 26 128 174 128 313S152 598 243 598Z" id="g50-49"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg>.</div>

Mathematical Problems in Engineering

fig11

Figure 11

Figure 11: A Fourier Continuation Method for the Solution of Elliptic Eigenvalue Problems in General Domains 