From: Manabendra Saharia on
Hello sir/Ma'am,

As a part of a research project, I have this Virus transport equation that I have to solve first in 1D and then in 2D using finite difference scheme. I am using Crank Nicholson Method. Please help me as I am stuck. I know I can implement in finite element using PDEPE but I need to do it in finite difference.

R dC/Dt = D(d^2 C/dx^2) - v(dC/dx) - Lambda*C - ( (Lambda*) * (pk/o)*C

where
R=retardation Coefficient= 1 + pk/o
C=Concentration,
p=density,
Lambda,lambda*= inactivation coeeficient.
All these values are known.

Boundary Conditions: C(x,0)=0; C(0,t)=C; dC(infinity,t)/dx=0.

I'll be very obliged if you help. Ask me if more data is required.