Consider a nonlinear advection-diffusion equation of the form

$$ \frac{\partial u}{\partial t} = \nabla \cdot (a(u) \nabla b(u) - \vec{c}(u)u) \tag{1} $$

on a rectangular domain with Dirichlet boundary conditions. We have $a(u) > 0$ for all $u \in \mathbb{R}$ and $b(u)$ is a monotonically increasing function of $u$. I wish to solve (1) by finite differences.

Can (1) always be discretized in such a way such that all matrix solves involve only symmetric positive definite matrices? If not, are there conditions on $a$, $b$, and $\vec{c}$ such that this is the case?

  • $\begingroup$ Can you specify your desired temporal discretization? $\endgroup$ – Charles Jun 11 '17 at 2:22
  • $\begingroup$ @Charlie Either Crank Nicholson or Backward Euler $\endgroup$ – eepperly16 Jun 11 '17 at 2:24

One strategy may be to use a semi-implicit temporal discretion, treating the diffusive term implicitly and the advective term explicitly. Doing so alleviates the extreme timestep restriction typically imposed by explicit diffusion, and it keeps your linear system symmetric. This is a common trick with (eg) the Navier Stokes equations to avoid precisely the same issues.

| cite | improve this answer | |

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.