Novel Second-Order Accurate Implicit Numerical Methods for the Riesz Space Distributed-Order Advection-Dispersion Equations

<div>Exact solutions and numerical solutions with <svg height="12.6917pt" id="M308" style="vertical-align:-3.40337pt" version="1.1" viewbox="-0.0498162 -9.28833 82.7981 12.6917" width="82.7981pt" 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="M558 454L532 470C512 451 497 442 465 442C421 442 376 450 318 450C146 450 24 313 24 164C24 46 91 -12 183 -12C317 -12 439 105 439 247C439 306 428 346 382 385L383 389C414 384 450 379 485 379C517 379 539 412 558 454ZM357 273C357 191 313 28 205 28C149 28 111 79 111 160C111 266 175 396 280 396C298 396 328 386 341 357C353 331 357 315 357 273Z" id="g50-240"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="576" vert-adv-y="576"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.166,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,22.784,0)"><path d="M525 275C525 381 464 451 380 451C338 451 289 439 241 405C158 347 127 277 106 152C51 -168 37 -216 24 -242L28 -257C42 -256 92 -247 117 -222C126 -213 134 -193 150 -69L166 32C184 10 211 -12 245 -12C276 -12 332 8 384 46C458 100 525 176 525 275ZM436 260C436 131 366 41 267 41C246 41 208 51 175 82C181 125 187 162 195 200C220 322 257 400 328 400C398 400 436 338 436 260Z" id="g50-239"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="549" vert-adv-y="549"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,33.599,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,45.218,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><g transform="matrix(.013,0,0,-0.013,51.549,0)"><path d="M381 704H317L48 -163H112L381 704Z" id="g50-48"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="429" vert-adv-y="429"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,57.131,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><g transform="matrix(.013,0,0,-0.013,63.462,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,69.793,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,76.124,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> and <svg height="11.4781pt" id="M309" style="vertical-align:-2.18977pt" version="1.1" viewbox="-0.0498162 -9.28833 75.1515 11.4781" width="75.1515pt" 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="M487 458L455 460C436 434 425 433 397 433C331 433 279 437 225 437C103 437 52 374 24 338L49 312C100 363 152 379 229 374C210 308 184 213 163 117C158 94 156 74 156 58C156 6 180 -12 211 -12C261 -12 313 27 348 68L332 99C307 76 285 64 267 64C252 64 236 86 250 148L301 373C320 373 344 369 364 369C399 369 435 370 450 372C458 386 470 411 487 458Z" id="g50-241"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="499" vert-adv-y="499"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,10.166,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,21.784,0)"><path d="M509 90L493 118C465 91 427 70 418 70S409 77 415 103L469 327C494 429 466 451 436 451C412 451 387 444 356 427C307 400 222 342 160 264H158L246 667C253 701 252 710 243 710C228 710 169 680 87 673L83 644H117C155 644 160 640 149 588L24 -4L32 -12C57 -4 86 4 108 9L143 187C198 275 324 386 375 386C394 386 399 367 385 311L331 87C315 22 322 -12 351 -12C387 -12 452 24 509 90Z" id="g50-105"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="527" vert-adv-y="527"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,32.31,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,43.929,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><g transform="matrix(.013,0,0,-0.013,50.26,0)"><path d="M381 704H317L48 -163H112L381 704Z" id="g50-48"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="429" vert-adv-y="429"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,55.842,0)"><path d="M290 377C321 398 342 415 358 430C378 450 389 473 389 502C389 578 329 635 238 635H237C184 635 137 610 109 578L64 515L88 493C112 529 154 573 208 573S303 542 303 482C303 409 233 370 141 341L149 308C165 313 190 319 215 319C272 319 341 283 341 193C342 98 292 43 222 43C163 43 122 72 96 94C88 101 79 100 70 94C61 87 47 73 46 60C44 47 48 37 62 23C76 10 118 -12 165 -12C238 -12 430 62 430 223C430 297 379 359 290 375V377Z" id="g50-52"></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,62.173,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,68.504,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> at <svg height="9.01192pt" id="M310" style="vertical-align:-0.2063904pt" version="1.1" viewbox="-0.0498162 -8.80553 39.6416 9.01192" width="39.6416pt" 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="M655 676H628C608 656 600 650 571 650H156C118 650 114 651 102 676H77C63 615 45 548 24 485L60 486C81 531 97 562 114 580C134 603 148 611 252 611H308L217 127C201 41 196 38 105 32L98 0H384L393 32C300 38 288 41 304 127L393 611H466C553 611 571 604 582 580C592 559 597 530 595 485L630 489C636 548 645 624 655 676Z" id="g50-85"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="673" vert-adv-y="673"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,12.064,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,23.683,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><g transform="matrix(.013,0,0,-0.013,30.014,0)"><path d="M117 -12C153 -12 178 12 178 50C178 84 153 110 118 110C84 110 59 84 59 50C59 12 84 -12 117 -12Z" id="g50-47"></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,33.095,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>

Advances in Mathematical Physics

fig2

Figure 2

Figure 2: Novel Second-Order Accurate Implicit Numerical Methods for the Riesz Space Distributed-Order Advection-Dispersion Equations 