Numerical Solution of the 1D Advection-Diffusion Equation Using Standard and Nonstandard Finite Difference Schemes

<table>Plot of <svg height="11.3125" id="M273" style="vertical-align:-0.10919pt" version="1.1" viewbox="0 0 29.825001 11.3125" width="29.825001" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.113)"><path d="M631 18l-4 -26q-24 0 -35 1q-57 3 -92.5 23t-69.5 68q-21 30 -101 160q-13 22 -30.5 31t-54.5 9h-32v-159q0 -63 14 -77.5t74 -19.5v-28h-262v28q62 5 76 19.5t14 77.5v402q0 62 -13 76t-73 19v28h248q109 0 159 -34q67 -42 67 -131q0 -115 -127 -172q27 -52 89 -144
q49 -71 81 -104q31 -36 72 -47zM212 316h48q73 0 108 29q55 42 55 124q0 76 -42 112t-104 36q-41 0 -54 -10q-11 -7 -11 -44v-247z" id="x52"></path></g><g transform="matrix(.017,-0,0,-.017,10.619,11.113)"><path d="M46 650h247q112 0 169 -41q64 -49 64 -141q0 -91 -56.5 -143t-138.5 -62q-23 -3 -45 -1l-74 19v-157q0 -63 15.5 -77t83.5 -19v-28h-271v28q60 5 74 20t14 77v400q0 65 -12.5 79t-69.5 18v28zM212 559v-245q31 -12 73 -12q59 0 102 38t43 124q0 77 -42.5 115.5
t-108.5 38.5q-44 0 -56 -12q-11 -8 -11 -47z" id="x50"></path></g><g transform="matrix(.017,-0,0,-.017,20.104,11.113)"><path d="M517 162l29 -4q-18 -100 -35 -158h-482v28q68 5 83.5 20t15.5 77v398q0 66 -13 80t-73 19v28h446q0 -8 8 -144l-29 -4q-11 51 -32 85q-11 17 -35 23t-78 6h-76q-23 0 -28 -5t-5 -27v-228h94q67 0 83 12t25 63h29v-191h-29q-10 54 -26 65.5t-81 11.5h-95v-189
q0 -65 24 -81q19 -13 104 -13q57 0 83 7t43 27q28 33 50 94z" id="x45"></path></g>
</svg> versus phase angle for the Crank-Nicolson scheme.</table>

Journal of Applied Mathematics

fig4

Figure 4

Figure 4: Numerical Solution of the 1D Advection-Diffusion Equation Using Standard and Nonstandard Finite Difference Schemes 