The equations in Sec. 5.10 are used to compute the velocity , ﬂux and pressure based on conservation of mass and momentum. At boundaries, the ﬂux corrector Eq. (5.18 ) must calculate in a manner consistent with and and their respective boundary conditions.
For example, at an impermeable stationary wall the calculated ﬂux must be , consistent with the no-slip condition . At boundary faces, and must be compatible to evaluate the correct according to Eq. (5.18 ).
The ﬁgure shows the fundamental boundary conditions from Sec. 4.3 and corresponding ﬂux evaluations. At the inlet and wall boundaries, is directly assigned from the boundary velocity by within the ﬂux corrector Eq. (5.18 ).
At the outlet, is not prescribed since is not a ﬁxed value condition. Instead, is evaluated from Eq. (5.18 ) using taken from cells adjacent to the boundary and calculated on the boundary.
When a body force is present in the momentum equation, the gradient condition for , e.g. at an inlet or wall, in principle becomes , as discussed in Sec. 4.4 . The precise details of the boundary condition in fact depend on how is incorporated within the coupling algorithm, discussed below.
With established, the condition is calculated based on the known by inverting Eq. (5.18 ). This approach causes to include a contribution from viscous stresses, as in Eq. (4.5 ), which may cause instability.