Omonov shahzod shoimovich

Ming-Fan Li, Ji-Rong Ren, and Tao Zhu., Fractional Vector Calculus and Fractional Special Function // Mathematical Physics.2010. arXiv:1001.2889 [math-ph].

Nigmatullin R. R. The realization of the generalized transfer equation in a medium with fractal geometry. Phys. Stat. Sol. (b). 133. 1986. 425–430.

Nkedi-Kizza P., Biggar J.W., Selim H.M., van Genuchen M.Th., Wierenga P.J., Davidson J.M., Nielsen D.R. On the equivalence of two concentual models for describing ion exchange during transport through an aggregated oxisol // Water Resour. Res. 20:1123-1130, 1984.

O’shaughnessy B., Procaccia I. Diffusion on fractals // Phys. Rev. A 32. 1985. Pp. 3073–3083.

Pachepsky Y., Benson D. and Rawls W. Simulating scale - dependent solute transport in soils with the fractional advective-dispersive equation // Soil Sci. Soc. Am. J. 4. 2000. Pp. 1234–1243.

Podlubny I. Fractional Differential Equations. Academic Press, San Diego. 1999.

Rahman M.M., Liedl R., Grathwohl P. Sorption kinetics during macropore transport of organic contaminants in soils: Laboratory experiments and analytical modeling // Water Resour. Res. 40. 2004. doi:10.1029/2002WR001946.

Rao P.S.C., Davidson J.M., Jessup R.E., Selim H.M. Evaluation of conceptual models for describing nonequilibrium adsorption-desorption of pesticides during steady-flow in soils // Soil Sci. Soc. Am. J., 43(1), 22-28, 1979.

Rao P.S.C., Jessup J. E., Rolston D.E., Davidson J.M., Kilgrease D.P. Experimental and mathematical description of nonadsorbed solute transfer by diffusion in spherical aggregates // Soil Sci. Soc. Am. J., Vol. 44. 1980a. Pp. 684 – 688.

Rao P.S.C., Jessup R.E., Addison T.M. Experimental and theoretical aspects of solute diffusion in spherical and non-spherical aggregates // Soil Sci. Soc. Am. J., 133. 1982. Pp. 342-349

Rao P.S.C., Rolston D.E., Jessup R. E., Davidson J.M. Solute transport in aggregated porous media: Theoretical and experimental evaluation. // Soil Sci. Soc. Am. J., Vol. 44. 1980b. Pp. 1139 – 1146.

Rasmuson A., Neretnieks I., Exact solution for diffusion in particles and longitudinal dispersion in packed beds. // AIChE J., 26. 1980. Pp. 686-690.

Rezanezhad F., Jonathan S.P., James R.G. The effects of dual porosity on transport and retardation inpeat: A laboratory experiment // Can. J. Soil Sci. 92. 2012. 723-732 p.

Ruiz-Medina M. D., Anh V. V., Angulo J. M. Fractional generalized random fields of variable order // Stoch. Anal. Appl. 22, 2004, pp. 775-799.

Samko S.G., Kilbas A.A., Marichev O.I. Fractional integrals and derivatives: Theory and applications.Gordon and breach. London. 1993.

Sanskrityayn A., Bharati V.K., Kumar N. Analytical solution of advection dispersion equation with spatiotemporal dependence of dispersion coefficient and velocity using greens function method // Journal of Groundwater Research. Vol. 5. No.1. 2016. Pp. 24-31.

Schumer R., Benson D.A., Meerschaert M.M., Baeumer B. Fractal mobile/immobile solute transport // Water Resources Research, Vol. 39, No. 10, 1296, doi:10.1029/2003WR002141, 2003.

Scotter D.R. Preferential solute movement through large soil voids, I. Some computations using simple theory // Australian Journal of Soil Research Vol. 16. 1978. Pp. 257-267.

Selim H.M. Transport & Fate of Chemicals in Soils: Principles & Applications. CRC Press Taylor & Francis Group. Boca Raton London New York. 2015.

Sharma P.K., Shukla S.K., Choudhary R., Swami D. Modeling for solute transport in mobile–immobilesoil column experiment // ISH Journal of Hydraulic Engineering, 2016.

Shen Ch., Phanikumar M.S. An efficient space-fractional dispersion approximation for stream solute transport modeling // Advances in Water Resources 32. 2009. Pp 1482–1494.

Singh H., Sahoo M.R.,·Singh O.P. Numerical Method Based on Galerkin Approximation for the Fractional Advection-Dispersion Equation // Int. J. Appl. Comput. Math. 3. 2017. Pp 2171–2187.

Skopp J., Gardner W.R., Tyler E.J. Solute movement in structured soils: Two-region model with small interaction // Soil Sci. Soc. Am. J., 45, 837-842, 1981.

Skopp J., Warrick A.W. A two-phase model for the miscible displacement of reactive solutes in soil // Soil Sci. Soc. Am. J. 38(4). 1974. Pp. 545-550.

Sousa E. How to approximate the fractional derivative of order . Proceedings of the 4th IFAC Workshop on Fractional Differentiation and its Applications. Badajoz. Spain. 2010.

Sudicky E.A., Frind E.O. Contaminant transport in fractured porous media: Analytical solutions for a system of parallel fractures // Water Resour. Res. 18(6). 1982. Pp.1634-1642.

Swami D., Sharma P.K. and Ojha C.S.P. Simulation of experimental breakthrough curves using multiprocess non-equilibrium model for reactive solute transport in stratified porous media // Indian Academy of Sciences. Sadhana Vol. 39, Part 6, 2014. Pp. 1425-1446.

Tang D.H., Frind E.O., Sudicky E.A. Contaminant transport in fractured porous media: Analytical solution for a single fracture // Water Re sour. Res. 17(3). 1981. Pp. 555-564.

Toride N., Inoue M., Leij, F. Hydrodynamic dispersion in an unsaturated dune sand // J. Soil Sci. Soc. Am., 67(3), 2003. 703–712p.

Van Genuchten M.Th., Tang D. H. Some Exact Solutions for Solute Transport Through Soils Containing Large Cylindrical Macropores // Water Resources Research. Vol. 20. No. 3. 1984. Pp. 335-346.

Van Genuchten M.Th., Wagenet R.J. Two-Site/Two-Region Models for Pesticide Transport and Degradation: Theoretical Development and Analytical Solutions // Soil Sci. Soc. American J.,Vol.53, 1989.1303-1310 p.

Van Genuchten M.Th., Werenga P.J. Mass Transfer Studies in Sorbing Porous Media: II. Experimental Evaluation with Tritium (H₂O) // Soil Sci. Soc. American J., Vol. 41, 1977. 272-278 p.

Van Genuchten M.Th., Wierenga P.J. Mass transfer studies in sorbing porous media. 1. Analytical solutions // Soil Sci. Soc. Am. J., 40(4), 1976, 473-480.

Van Genuchten M.Th., Wierenga P.J., O’Conner G.A. Mass Transfer Studies in Sorbing Porous Media: III. Experimental Evaluation with 2,4,5-T // Soil Sci. Soc. American J., Vol. 41, 1977. 278-285 p.

Villermaux J., and W. P. M. van Swaajj, Modele repre"sentatif de la distribution des temps de sejour dans un re"acteur semi-infmi a dispersion axiale avec zones stagnarites. Application a I'fecoule-ment ruisselant des colonnes d'anneaux Raschig // Chem. Eng. Sci.,24, 1097-1111, 1969.

Weihua D. Numerical algorithm for the time fractional Fokker-Planck equation // Journal of Computational Physics. 227. 2007. Pp. 1510-1522.

Xia Y. Wu J. C. Numerical Solutions of Fractional Advection–Dispersion Equations // Journal of Nanjing University (Natural Sciences). 43(4). 2007. Pp. 441-446.

Xia Yuan, Wu Jichun, Zhou Luying Numerical solutions of time-space fractional advection–dispersion equations // ICCES, Vol.9, no.2. 2009. Pp. 117−126.

Yadav R.R., Kumar L.K. Two-Dimensional Conservative Solute Transport with Temporal and Scale-Dependent Dispersion: Analytical Solution. // International Journal of Advance in Mathematics.Vol 2018(2). 2018. Pp 90-111.

Young D.F., Ball W.P. Estimating diffusion coefficients in low-permeability porous media using a macropore column. // Environ. Sci. Technol. Vol. 32. 1998. Pp. 2578 – 2584.

Young D.F., William P.B. Effects of column conditions on the first-order rate modeling of nonequilibrium solute breakthrough: Cylindrical macropores versus spherical media // Water Resour.Res. Vol. 33. 1997. 1149 –1156.

Young D.F., William P.B. Two-region linear/nonlinear sorption modeling: Batch and column experiments // Environ. Toxicol. Chem.Vol.18. 1999. Pp. 1686 – 1693.

Zhang Y., Benson D.A., Reeves D.M. Time and space nonlocalities underlying fractional-derivative models: Distinction and literature review of field applications // Advances in Water Resources. Vol. 32. 2009. Pp. 561–581.

Zhou L., Selim H.M. Application of the fractional advection-dispersion equation in porous media // Soil Sci. Soc. Am. J. 67. 2003. Pp. 1079–1084.

Бейбалаев В.Д., Шабанова М.Р. Численный метод решения начально-граничной задачи для двумерного уравнения теплопроводности с производными дробного порядка // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2010. Выпуск 5(21). С. 244–251.

ILOVALAR.

function Aziz;

%o'zgarmaslarni qiymatlari

C0=0.01; Da=1e-6; Vm=1e-4; Dm=1e-5;

n=50; m=40;
aa=0.1; bb=0.5;
tetf=1.0; teta=0.6;
tau=1;
h1=0.035;

tmax=3600;

bet1=1;
bet=2;

%boshlang'ich va chegaraviy shart

Ca(1:n,1:m,1:tmax)=0;

Cim(1:n,1:tmax)=0;
Cm(1:n,1:tmax)=0;
Cm(1,1:tmax)=C0;
Ca(1:n,1,1)=Cm(1:n,1);
Vf=(aa*aa)/(bb*bb);
Va=1-Vf;
tetm=Vf*tetf
tetim=Va*teta
h2=(bb-aa)/m;

for k=1:tmax-1

for i=1:n
Ca(i,1,k)=Cm(i,k);
for j=2:m-1
if (2*bet1)==2
s0=Ca(i,j+1,k)-2*Ca(i,j,k)+Ca(i,j-1,k);
else
s0=0;
for l=1:j-1
s0=s0+((l-1+1)^(2-2*bet1)-(l-1)^(2-2*bet1))*(Ca(i,j+2-l,k)-2*Ca(i,j+1-l,k)+Ca(i,j-l,k));
end;
end;
Ca(i,j,k+1)=Ca(i,j,k)+tau*Da/((j*h2)^bet1)*(((j+1/2)^bet1-bet1*(j-1/2)^bet1)/gamma(2-bet1)*...
(Ca(i,j,k)-bet1*Ca(i,j-1,k))/(gamma(2-bet1)*h2^bet1)+...
(j*h2)^bet1/(gamma(3-2*bet1)*h2^(2*bet1))*s0);
end;
Ca(i,m,k+1)=Ca(i,m-1,k+1);
Ca(i,1,k+1)=Cm(i,k);

end;
q1(k)=-Da*((Ca(1,2,k)-bet1*Ca(1,1,k)))/(gamma(2-bet1)*h2^bet1);

q2(k)=-Da*((Ca(4,2,k)-bet1*Ca(4,1,k)))/(gamma(2-bet1)*h2^bet1);
q3(k)=-Da*((Ca(10,2,k)-bet1*Ca(10,1,k)))/(gamma(2-bet1)*h2^bet1);
A1(k)=Cm(1,k);
A2(k)=Cm(4,k);
A3(k)=Cm(6,k);
A4(k)=Cm(11,k);
B1(k)=Cim(1,k);
B2(k)=Cim(4,k);
B3(k)=Cim(6,k);
B4(k)=Cim(11,k);
for i=1:n
ff=h2/2*(0*Ca(i,1,k+1)+(bb-aa)*Ca(i,m,k+1));
for j=1:m-1
ff=ff+h2*(0+j*h2)*Ca(i,j,k+1);
end;
Cim(i,k+1)=2*ff/((bb-aa)^2-0);
end;
tet=tetim/tetm;
% G0=gamma(2-gam)*tau^gam;
D=tau*Dm/(gamma(3-bet)*h1^bet);
Cm(1,k)=C0;
for i=2:n-1
if bet==2
s0=Cm(i+1,k)-2*Cm(i,k)+Cm(i-1,k);
else
s0=0;
for l=1:i-1
s0=s0+((l-1+1)^(2-bet)-(l-1)^(2-bet))*(Cm(i+2-l,k)-2*Cm(i+1-l,k)+Cm(i-l,k));
end;
end;
Cm(i,k+1)=D*s0-tau*Vm/h1*(Cm(i,k)-Cm(i-1,k))-tet*(Cim(i,k+1)-Cim(i,k))+Cm(i,k);
end;
Cm(n,k+1)=Cm(n-1,k+1);
k
end;

%plot((0:l)*h2,Ca(1:l1,tmax,1:l1));

%Ca(1:n,tmax-1,1:l1);
mesh((0:m-1)*h1,(0:m-1)*h2,Ca(1:m,1:m,tmax-1)'),
set(gca,'FontSize',12);
% Cim(1:n,tmax/2);
%Cim(1:n,900)
%Cim(1:n,1800)
%Cim(1:n,3600)
filename1 = 'q1.xlsx';
sheet = 1;
xlRange = 'A1';
xlswrite(filename1,q1',sheet,xlRange);
xlRange = 'B1';
xlswrite(filename1,q2',sheet,xlRange);
xlRange = 'C1';
xlswrite(filename1,q3',sheet,xlRange);
filename2 = 'Caa.xlsx';
xlswrite(filename2,Ca(1:n,1:m,tmax-1));
filename3 = 'A.xlsx';
sheet = 1;
xlRange = 'A1';
xlswrite(filename3,A1',sheet,xlRange);
xlRange = 'B1';
xlswrite(filename3,A2',sheet,xlRange);
xlRange = 'C1';
xlswrite(filename3,A3',sheet,xlRange);
xlRange = 'D1';
xlswrite(filename3,A4',sheet,xlRange);
filename4 = 'B.xlsx';
sheet = 1;
xlRange = 'A1';
xlswrite(filename4,B1',sheet,xlRange);
xlRange = 'B1';
xlswrite(filename4,B2',sheet,xlRange);
xlRange = 'C1';
xlswrite(filename4,B3',sheet,xlRange);
xlRange = 'D1';
xlswrite(filename4,B4',sheet,xlRange);
%Cm(1:n,900)
%Cm(1:n,1800)
Cm(1:n,3600)
hold on;