Hello,

I would like to numerically solve the following system of PDEs using a stable computational scheme. The system involves three equations and three unknowns ($\displaystyle q_x$, $\displaystyle q_y$, and $\displaystyle p$). I've tried to solve this system using an explicit scheme, but I found that the solution could become unstable. An implicit scheme may work, but then how would I deal with the three unknowns in a computationally efficient fashion?

Could someone suggest a stable scheme to numerically solve this system?

$\displaystyle

\[

A\frac{{\partial p}}{{\partial t}} + Bp = \frac{{\partial q_x }}{{\partial x}} + \frac{{\partial q_y }}{{\partial y}}

\]

$

$\displaystyle

\[

\frac{{\partial q_x }}{{\partial t}} = \frac{{\partial p}}{{\partial x}}

\]

$

$\displaystyle

\[

\frac{{\partial q_y }}{{\partial t}} = \frac{{\partial p}}{{\partial y}}

\]

$