Numerical Analysis of a Fast Finite Element Method for a Hidden-Memory Variable-Order Time-Fractional Diffusion Equation

Jinhong Jia,Hong Wang,Xiangcheng Zheng
DOI: https://doi.org/10.1007/s10915-022-01820-z
2022-04-04
Journal of Scientific Computing
Abstract:<p class="a-plus-plus">We investigate a fast finite element scheme to a hidden-memory variable-order time-fractional diffusion equation. Different from the traditional L1 methods, a fast approximation to the hidden-memory variable-order fractional derivative is derived to reduce the computational cost of generating coefficients from <span class="a-plus-plus inline-equation id-i-eq1"><span class="a-plus-plus equation-source format-t-e-x"><span class="mjpage"><svg xmlns:xlink="http://www.w3.org/1999/xlink" width="6.76ex" height="3.009ex" style="vertical-align: -0.838ex;" viewBox="0 -934.9 2910.5 1295.7" role="img" focusable="false" xmlns="http://www.w3.org/2000/svg"><g stroke="currentColor" fill="currentColor" stroke-width="0" transform="matrix(1 0 0 -1 0 0)"> <use xlink:href="#MJMATHI-4F" x="0" y="0"></use> <use xlink:href="#MJMAIN-28" x="763" y="0"></use><g transform="translate(1153,0)"> <use xlink:href="#MJMATHI-4E" x="0" y="0"></use> <use transform="scale(0.707)" xlink:href="#MJMAIN-32" x="1292" y="513"></use></g> <use xlink:href="#MJMAIN-29" x="2520" y="0"></use></g></svg></span></span></span> to <span class="a-plus-plus inline-equation id-i-eq2"><span class="a-plus-plus equation-source format-t-e-x"><span class="mjpage"><svg xmlns:xlink="http://www.w3.org/1999/xlink" width="11.456ex" height="2.843ex" style="vertical-align: -0.838ex;" viewBox="0 -863.1 4932.3 1223.9" role="img" focusable="false" xmlns="http://www.w3.org/2000/svg"><g stroke="currentColor" fill="currentColor" stroke-width="0" transform="matrix(1 0 0 -1 0 0)"> <use xlink:href="#MJMATHI-4F" x="0" y="0"></use> <use xlink:href="#MJMAIN-28" x="763" y="0"></use> <use xlink:href="#MJMATHI-4E" x="1153" y="0"></use><g transform="translate(2208,0)"> <use xlink:href="#MJMAIN-6C"></use> <use xlink:href="#MJMAIN-6F" x="278" y="0"></use> <use xlink:href="#MJMAIN-67" x="779" y="0"></use></g> <use xlink:href="#MJMATHI-4E" x="3654" y="0"></use> <use xlink:href="#MJMAIN-29" x="4542" y="0"></use></g></svg></span></span></span>, where <em class="a-plus-plus">N</em> refers to the number of time steps. We then develop different techniques from the analysis of L1 methods to prove error estimates for the corresponding fast fully-discrete finite element scheme. Furthermore, a fast divide and conquer algorithm is proposed to reduce the complexity of solving the linear systems from <span class="a-plus-plus inline-equation id-i-eq3"><span class="a-plus-plus equation-source format-t-e-x"><span class="mjpage"><svg xmlns:xlink="http://www.w3.org/1999/xlink" width="9.202ex" height="3.009ex" style="vertical-align: -0.838ex;" viewBox="0 -934.9 3962 1295.7" role="img" focusable="false" xmlns="http://www.w3.org/2000/svg"><g stroke="currentColor" fill="currentColor" stroke-width="0" transform="matrix(1 0 0 -1 0 0)"> <use xlink:href="#MJMATHI-4F" x="0" y="0"></use> <use xlink:href="#MJMAIN-28" x="763" y="0"></use> <use xlink:href="#MJMATHI-4D" x="1153" y="0"></use><g transform="translate(2204,0)"> <use xlink:href="#MJMATHI-4E" x="0" y="0"></use> <use transform="scale(0.707)" xlink:href="#MJMAIN-32" x="1292" y="513"></use></g> <use xlink:href="#MJMAIN-29" x="3572" y="0"></use></g></svg></span></span></span> to <span class="a-plus-plus inline-equation id-i-eq4"><span class="a-plus-plus equation-source format-t-e-x"><span class="mjpage"><svg xmlns:xlink="http://www.w3.org/1999/xlink" width="14.952ex" height="3.176ex" style="vertical-align: -0.838ex;" viewBox="0 -1006.6 6437.7 1367.4" role="img" focusable="false" xmlns="http://www.w3.org/2000/svg"><g stroke="currentColor" fill="currentColor" stroke-width="0" transform="matrix(1 0 0 -1 0 0)"> <use xlink:href="#MJMATHI-4F" x="0" y="0"></use> <use xlink:href="#MJMAIN-28" x="763" y="0"></use> <use xlink:href="#MJMATHI-4D" x="1153" y="0"></use> <use xlink:href="#MJMATHI-4E" x="2204" y="0"></use><g transform="translate(3259,0)"> <use xlink:href="#MJMAIN-6C"></use> <use xlink:href="#MJMAIN-6F" x="278" y="0"></use> <use xlink:href="#MJMAIN-67" x="779" y="0"></use> <use transform="scale(0.707)" xlink:href="#MJMAIN-32" x="1809" y="598"></use></g> <use xlink:href="#MJMATHI-4E" x="5159" y="0"></use> <use xlink:href="#MJMAIN-29" x="6048" y="0"></use></g></svg></span></span></span> where <em class="a-plus-plus">M</em> stands for the spatial degree of freedom. Numerical experiments are presented to substantiate the theoretical results.</p><svg xmlns="http://www.w3.org/2000/svg" style="display: none;"><defs id="MathJax_SVG_glyphs"><path stroke-width="1" id="MJMATHI-4F" d="M740 435Q740 320 676 213T511 42T304 -22Q207 -22 138 35T51 201Q50 209 50 244Q50 346 98 438T227 601Q351 704 476 704Q514 704 524 703Q621 689 680 617T740 435ZM637 476Q637 565 591 615T476 665Q396 665 322 605Q242 542 200 428T157 216Q157 126 200 73T314 19Q404 19 485 98T608 313Q637 408 637 476Z"></path><path stroke-width="1" id="MJMAIN-28" d="M94 250Q94 319 104 381T127 488T164 576T202 643T244 695T277 729T302 750H315H319Q333 750 333 741Q333 738 316 720T275 667T226 581T184 443T167 250T184 58T225 -81T274 -167T316 -220T333 -241Q333 -250 318 -250H315H302L274 -226Q180 -141 137 -14T94 250Z"></path><path stroke-width="1" id="MJMATHI-4E" d="M234 637Q231 637 226 637Q201 637 196 638T191 649Q191 676 202 682Q204 683 299 683Q376 683 387 683T401 677Q612 181 616 168L670 381Q723 592 723 606Q723 633 659 637Q635 637 635 648Q635 650 637 660Q641 676 643 679T653 683Q656 683 684 682T767 680Q817 680 843 681T873 682Q888 682 888 672Q888 650 880 642Q878 637 858 637Q787 633 769 597L620 7Q618 0 599 0Q585 0 582 2Q579 5 453 305L326 604L261 344Q196 88 196 79Q201 46 268 46H278Q284 41 284 38T282 19Q278 6 272 0H259Q228 2 151 2Q123 2 100 2T63 2T46 1Q31 1 31 10Q31 14 34 26T39 40Q41 46 62 46Q130 49 150 85Q154 91 221 362L289 634Q287 635 234 637Z"></path><path stroke-width="1" id="MJMAIN-32" d="M109 429Q82 429 66 447T50 491Q50 562 103 614T235 666Q326 666 387 610T449 465Q449 422 429 383T381 315T301 241Q265 210 201 149L142 93L218 92Q375 92 385 97Q392 99 409 186V189H449V186Q448 183 436 95T421 3V0H50V19V31Q50 38 56 46T86 81Q115 113 136 137Q145 147 170 174T204 211T233 244T261 278T284 308T305 340T320 369T333 401T340 431T343 464Q343 527 309 573T212 619Q179 619 154 602T119 569T109 550Q109 549 114 549Q132 549 151 535T170 489Q170 464 154 447T109 429Z"></path><path stroke-width="1" id="MJMAIN-29" d="M60 749L64 750Q69 750 74 750H86L114 726Q208 641 251 514T294 250Q294 182 284 119T261 12T224 -76T186 -143T145 -194T113 -227T90 -246Q87 -249 86 -250H74Q66 -250 63 -250T58 -247T55 -238Q56 -237 66 -225Q221 -64 221 250T66 725Q56 737 55 738Q55 746 60 749Z"></path><path stroke-width="1" id="MJMAIN-6C" d="M42 46H56Q95 46 103 60V68Q103 77 103 91T103 124T104 167T104 217T104 272T104 329Q104 366 104 407T104 482T104 542T103 586T103 603Q100 622 89 628T44 637H26V660Q26 683 28 683L38 684Q48 685 67 686T104 688Q121 689 141 690T171 693T182 694H185V379Q185 62 186 60Q190 52 198 49Q219 46 247 46H263V0H255L232 1Q209 2 183 2T145 3T107 3T57 1L34 0H26V46H42Z"></path><path stroke-width="1" id="MJMAIN-6F" d="M28 214Q28 309 93 378T250 448Q340 448 405 380T471 215Q471 120 407 55T250 -10Q153 -10 91 57T28 214ZM250 30Q372 30 372 193V225V250Q372 272 371 288T364 326T348 362T317 390T268 410Q263 411 252 411Q222 411 195 399Q152 377 139 338T126 246V226Q126 130 145 91Q177 30 250 30Z"></path><path stroke-width="1" id="MJMAIN-67" d="M329 409Q373 453 429 453Q459 453 472 434T485 396Q485 382 476 371T449 360Q416 360 412 390Q410 404 415 411Q415 412 416 414V415Q388 412 363 393Q355 388 355 386Q355 385 359 381T368 369T379 351T388 325T392 292Q392 230 343 187T222 143Q172 143 123 171Q112 153 112 133Q112 98 138 81Q147 75 155 75T227 73Q311 72 335 67Q396 58 431 26Q470 -13 470 -72Q470 -139 392 -175Q332 -206 250 -206Q167 -206 107 -175Q29 -140 29 -75Q29 -39 50 -15T92 18L103 24Q67 55 67 108Q67 155 96 193Q52 237 52 292Q52 355 102 398T223 442Q274 442 318 416L329 409ZM299 343Q294 371 273 387T221 404Q192 404 171 388T145 343Q142 326 142 292Q142 248 149 227T179 192Q196 182 222 182Q244 182 260 189T283 207T294 227T299 242Q302 258 302 292T299 343ZM403 -75Q403 -50 389 -34T348 -11T299 -2T245 0H218Q151 0 138 -6Q118 -15 107 -34T95 -74Q95 -84 101 -97T122 -127T170 -155T250 -167Q319 -167 361 -139T403 -75Z"></path><path stroke-width="1" id="MJMATHI-4D" d="M289 629Q289 635 232 637Q208 637 201 638T194 648Q194 649 196 659Q197 662 198 666T199 671T201 676T203 679T207 681T212 683T220 683T232 684Q238 684 262 684T307 683Q386 683 398 683T414 678Q415 674 451 396L487 117L510 154Q534 190 574 254T662 394Q837 673 839 675Q840 676 842 678T846 681L852 683H948Q965 683 988 683T1017 684Q1051 684 1051 673Q1051 668 1048 656T1045 643Q1041 637 1008 637Q968 636 957 634T939 623Q936 618 867 340T797 59Q797 55 798 54T805 50T822 48T855 46H886Q892 37 892 35Q892 19 885 5Q880 0 869 0Q864 0 828 1T736 2Q675 2 644 2T609 1Q592 1 592 11Q592 13 594 25Q598 41 602 43T625 46Q652 46 685 49Q699 52 704 61Q706 65 742 207T813 490T848 631L654 322Q458 10 453 5Q451 4 449 3Q444 0 433 0Q418 0 415 7Q413 11 374 317L335 624L267 354Q200 88 200 79Q206 46 272 46H282Q288 41 289 37T286 19Q282 3 278 1Q274 0 267 0Q265 0 255 0T221 1T157 2Q127 2 95 1T58 0Q43 0 39 2T35 11Q35 13 38 25T43 40Q45 46 65 46Q135 46 154 86Q158 92 223 354T289 629Z"></path></defs></svg>
mathematics, applied
What problem does this paper attempt to address?