Analytical Solutions of the Space-Time Fractional Derivative of Advection Dispersion Equation

Fractional advection-dispersion equations are used in groundwater hydrology to model the transport of passive tracers carried by fluid flow in porous medium. A space-time fractional advection-dispersion equation (FADE) is a generalization of the classical ADE in which the first-order space derivative is replaced with Caputo or Riemann-Liouville derivative of order <svg style="vertical-align:-3.56265pt;width:65.162498px;" id="M1" height="16.6625" version="1.1" viewBox="0 0 65.162498 16.6625" width="65.162498" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns="http://www.w3.org/2000/svg"> <g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path id="x30" d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105 q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" /></g><g transform="matrix(.017,-0,0,-.017,12.931,12.162)"><path id="x3C" d="M512 -3l-437 233v51l437 233v-58l-378 -200v-2l378 -199v-58z" /></g><g transform="matrix(.017,-0,0,-.017,27.634,12.162)"><path id="x1D6FD" d="M558 587q0 -32 -14 -61t-40 -53.5t-48.5 -41t-54.5 -36.5q144 -51 144 -174q0 -55 -43.5 -108t-104.5 -87q-77 -42 -131 -42q-31 0 -54 20t-31 47l11 18q48 -29 108 -29q79 0 119.5 43t40.5 109t-44.5 107.5t-119.5 50.5l22 47q34 1 65 21q96 61 96 157q0 42 -24 67.5 t-62 25.5q-24 0 -43.5 -9t-35 -29.5t-27 -44t-22.5 -63t-19.5 -75.5t-18.5 -91q-57 -294 -68 -380q-26 -190 -35 -200q-26 -31 -97 -37l-4 26q19 9 48 170l77 413q23 121 52.5 187.5t83.5 114.5q70 62 148 62q51 0 88.5 -34t37.5 -91z" /></g><g transform="matrix(.017,-0,0,-.017,42.236,12.162)"><path id="x2264" d="M531 71l-474 214v50l474 215v-56l-416 -184l416 -183v-56zM531 -40h-474v50h474v-50z" /></g><g transform="matrix(.017,-0,0,-.017,56.94,12.162)"><path id="x31" d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" /></g> </svg>, and the second-order space derivative is replaced with the Caputo or the Riemann-Liouville fractional derivative of order <svg style="vertical-align:-0.546pt;width:64.9375px;" id="M2" height="11.6" version="1.1" viewBox="0 0 64.9375 11.6" width="64.9375" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns="http://www.w3.org/2000/svg"> <g transform="matrix(.017,-0,0,-.017,.062,10.862)"><use xlink:href="#x31"/></g><g transform="matrix(.017,-0,0,-.017,12.931,10.862)"><use xlink:href="#x3C"/></g><g transform="matrix(.017,-0,0,-.017,27.634,10.862)"><path id="x1D6FC" d="M545 106q-67 -118 -134 -118q-24 0 -40 37.5t-30 129.5h-2q-47 -72 -103 -119.5t-108 -47.5q-47 0 -76 45.5t-29 119.5q0 113 85 204t174 91q47 0 70 -33.5t43 -119.5h3q32 47 80 140l55 13l10 -9q-47 -80 -138 -201q17 -99 27.5 -136t22.5 -37q23 0 69 61zM333 204 q-14 98 -31 149.5t-50 51.5q-49 0 -94 -70t-45 -164q0 -55 15.5 -86t40.5 -31q70 0 164 150z" /></g><g transform="matrix(.017,-0,0,-.017,42.015,10.862)"><use xlink:href="#x2264"/></g><g transform="matrix(.017,-0,0,-.017,56.719,10.862)"><path id="x32" d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" /></g> </svg>. We derive the solution of the new equation in terms of Mittag-Leffler functions using Laplace transfrom. Some examples are given. The results from comparison let no doubt that the FADE is better in prediction than ADE.