|
Basilisk CFD
Adaptive Cartesian mesh PDE framework
|
#include <henry.h>
Data Fields | |
| vector | D |
| vector | beta |
We first define the relaxation and residual functions needed to solve the implicit discrete system
\[ \frac{c^{n + 1} - c^n}{\Delta t} = \nabla\cdot (D \nabla c^{n + 1} + \beta c^{n + 1}) \]
see section 3.2 of Farsoiya et al., 2021.
Note that these functions are close to that in [poisson.h]() and [diffusion.h]() but with the additional term \(\nabla\cdot (\beta c^{n + 1})\).