I am trying to discretise the PDE: $$\phi \frac{\partial c}{\partial t}+\frac{\partial j}{\partial x}=0$$ where $c$ is a function of $x$ and $t$, and $j=qc-D(q)\frac{\partial c}{\partial x}$, $q$ is the Darcy flux (function of $x$ and $t$) and $D(q)$ is the coefficient of diffusion/dispersion, and $\phi$ is a constant.

over the interval $(x_{i-1/2}, x_{i+1/2})$

My Progress

$$\begin{align*} &\phi \frac{\partial c}{\partial t}+\frac{\partial j}{\partial x}=0 \\ &\Rightarrow \frac{\partial c}{\partial t}+\frac{\partial}{\partial x}\left(\frac{qc}{\phi}-D(q)\frac{\partial c}{\partial x}\right)=0 \\&\Rightarrow \frac{\partial}{\partial t}\int_{x_{i-1/2}}^{x_{i+1/2}}c \ dx \ +\left|\frac{qc}{\phi}-D(q)\frac{\partial c}{\partial x}\right|_{x_{i-1/2}}^{x_{i+1/2}}=0 \end{align*}$$

What do I do to simplify the expression in terms of the values of the $c_i, q_i$ terms. (I am trying to get matlab to run this iteratively).


$x_i$ is the midpoint of the $i$th interval.

$\delta x$ is the length of the interval.

$$f_i=\frac1{\delta x}\int_{x_{i-1/2}}^{x_{i+1/2}}f \ dx$$



Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Browse other questions tagged or ask your own question.