Need MATLAB program for 2DFATMIC code for 1D reactant solute transport with distance dépend flow and dispersion and a pulse injection

4 次查看（过去 30 天）

Thomas TJOCK-MBAGA 2022-5-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1723070-need-matlab-program-for-2dfatmic-code-for-1d-reactant-solute-transport-with-distance-depend-flow-and

评论： Torsten 2022-5-22

I would like to applied the 2DFATMIC code to one-dimensional advection-dispersion transport with variable velocity flow and dispersion coefficient. In my problem (single-well push-test applying to PCE reactant solute) the inlet boundary condition is assumed to be a pulse type (first type), the contamination is assumed to be injected until t = 5 days. I want to obtain the concentration distribution during the injection phase and after. But i don't have the 2DFATMIC code. And i wanted to use this code un MATLAB. Thanks !!! I have the following transport équation:

R*dC/dt = d/dx*(D(x)*dC/dx-v(x)*C)-mu*C The velocity is express like this v(x) = 30/(2*pi*10*0.1*x); thé dispersion coefficient is equal to D(x)= 0.1*v(x)

mu is the first order reaction coefficient (mu = 0.075 day^(-1)) Initial condition C(x,t=0) = 0

Boundaries conditions

C(x=0,t) = 1 , 0=<t<t0 and 0 , t> t0 where t0 = 5 days dC(x=L, t) = 0

L is the length of the domain L = 10 m

If it is not possible to obtain the 2DFMTMIC code i wanted as well as possible the cranck Nicholson code or other schème code.