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

<table>Plot of <svg height="11.425" id="M377" style="vertical-align:-0.0pt" version="1.1" viewbox="0 0 36.012501 11.425" width="36.012501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.363)"><path d="M673 0h-245v28q47 6 58 16t1 39l-55 154h-213q-31 -87 -36 -129q-16 -47 -5 -61t68 -19v-28h-226v28q50 7 68.5 23.5t43.5 80.5l214 524l24 9l193 -536q23 -63 41.5 -79.5t68.5 -21.5v-28zM418 280l-81 255h-2q-74 -197 -101 -255h184z" id="x41"></path></g><g transform="matrix(.017,-0,0,-.017,11.809,11.363)"><path d="M493 503l-29 -4q-13 56 -28 80q-13 22 -35.5 29.5t-78.5 7.5h-76q-22 0 -27.5 -5t-5.5 -27v-238h103q66 0 82.5 12t25.5 62h29v-193h-29q-9 54 -25.5 67t-82.5 13h-103v-181q0 -64 15 -78.5t84 -19.5v-28h-272v28q61 5 74.5 19.5t13.5 78.5v399q0 64 -13 78t-72 19v28
h440q1 -40 10 -147z" id="x46"></path></g><g transform="matrix(.017,-0,0,-.017,20.801,11.363)"><path d="M861 0h-263v28q61 5 73 22t12 97l-4 408h-4l-244 -548h-20l-224 534h-2l-15 -285q-5 -98 -3 -166q1 -34 19.5 -45t72.5 -17v-28h-233v28q45 5 62.5 17t22.5 43q11 65 19 183l14 208q3 36 3.5 55t-2.5 37.5t-6.5 25t-16.5 13t-25.5 8t-39.5 4.5v28h164l228 -501l232 501
h170v-28q-68 -6 -81 -21.5t-11 -91.5l9 -362q2 -80 15 -97t78 -22v-28z" id="x4D"></path></g>
</svg> and <svg height="11.3125" id="M378" 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 NSFD scheme.</table>

Journal of Applied Mathematics

fig6

Figure 6

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