3.4 Equation discretisation
Equation discretisation converts partial diﬀerential equations for continuous ﬁelds, e.g. pressure , into sets of linear equations for discrete ﬁelds.
The values of the principal ﬁelds, e.g. , are associated with cells. A ﬁeld is then represented by an array of values, , for cell indices . is the total number of cells.
Equation discretisation creates a linear equation for each cell. For cell 43 above, the equation might have the following form:
where and are coeﬃcients corresponding to cell indices , (diagonal coeﬃcient in bold). The set of linear equations for all cells can be written as a matrix equation of the form:
The matrix equation can be represented as , where: are the matrix coeﬃcients; are the source coeﬃcients ; and, is the discretised pressure ﬁeld. It may also be illustrated as follows, where ‘’ indicates non-zero coeﬃcients.
A CFD simulation generally solves a set of physical equations, e.g. for mass, momentum and energy conservation. The ﬁnite volume method traditionally discretises each physical equation separately to form individual matrix equations for single solution variables, e.g. , , , rather than creating a single matrix equation that represents all the physical equations.
The segregated matrix equations are solved one variable at a time, e.g. solving for , and in separate steps. Where the solution variable is a vector or tensor, e.g. , it is decoupled into individual matrix equations for each component, e.g. , , .
The matrix equations are solved in an iterative sequence, in which the equation for one variable, e.g. , incorporates current values of other variables, e.g. , and , into the source vector, as shown below.