Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

I am using a box-scheme for solving partial differential equation. The function is approximated with: $$ f_p=\psi\cdot\left(\theta\cdot f_{j}^{k+1}+(1-\theta)\cdot f_{j}^{k} \right) + (1-\psi)\cdot \left(\theta\cdot f_{j+1}^{k+1}+(1-\theta)\cdot f_{j+1}^{k} \right) $$ Time derivative is approximated with: $$ \psi\frac{f_j^{k+1}-f_j^k}{\Delta t}+(1-\psi)\frac{f_{j+1}^{k+1}-f_{j+1}^k}{\Delta t} $$ Time spatial derivative is approximated with: $$ (1-\theta)\frac{f_{j+1}^{k}-f_j^k}{\Delta x}+\theta\frac{f_{j+1}^{k+1}-f_{j}^{k+1}}{\Delta x} $$

My question is how to approximate with box-scheme following PDE: $$ \frac{\partial f}{\partial t} + \frac{\partial}{\partial x} \left(\sqrt{c- \frac{\partial f}{\partial x}} \right) $$

I am stuck at that moment at which I approximated the outer spatial derivative: $$ \frac{\partial}{\partial x} \left(\sqrt{c- \frac{\partial f}{\partial x}} \right)\approx (1-\theta)\frac{\left(\sqrt{c- \frac{\partial f}{\partial x}} \right)_{j+1}^{k}-\left(\sqrt{c- \frac{\partial f}{\partial x}} \right)_j^k}{\Delta x}+\theta\frac{\left(\sqrt{c- \frac{\partial f}{\partial x}} \right)_{j+1}^{k+1}-\left(\sqrt{c- \frac{\partial f}{\partial x}} \right)_{j}^{k+1}}{\Delta x} $$ Where:

$c$ - constant, $t$ - time, $x$ - spatial variable, $j$ - space node index, $k$ - time node index.

What is the best method to approximate the derivatives inside the square root?

share|improve this question

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.