3.5 Matrix construction
The construction of each matrix equation involves building the matrix and source coeﬃcients from the terms in the equation being solved, with further adjustments for boundary conditions.
The ﬁgure below illustrates the process of building a matrix equation for a ﬁeld from an equation including advection, diﬀusion and source .
The overall coeﬃcients are calculated as the sum of coeﬃcients for each term in the the equation. Most terms contribute to both the matrix and source coeﬃcients, although this depends on the choice of discretisation scheme.
Finally, boundary conditions are incorporated into the equation through further adjustments to coeﬃcients in and as shown below. The adjustments, principally from the advection and diﬀusion terms, are applied to coeﬃcients corresponding to cells at the domain boundary.
Implicit and explicit
The equation for a ﬁeld on the previous page is discretised to form the matrix equation . The discretisation of a term is implicit when it contributes to coeﬃcients in by treating as the solved ﬁeld .
Explicit discretisation calculates coeﬃcients in only, using current values of ﬁelds. When solving an equation for , derivatives without must be explicit. Terms with could be treated explicitly by using current values of , but they generally are not, since explicit solutions are unstable beyond a limiting time step as described in Sec. 3.17 . A notable exception to this are the terms discussed in Sec. 3.20 .
The curl derivative, e.g. , includes terms in and in the decoupled matrix equation for . These terms must be treated explicitly since they do not include itself. The situation for and is the same, so the curl derivative can be explicit only.
This leaves the following terms which are generally treated implicitly:3
This chapter details: the discretisation of these terms, which are generally treated implicitly; and, other terms such as divergence and gradient which can only be discretised explicitly within a segregated solution.