
FLUENT's controlvolume formulation requires that convection and diffusion fluxes through the control volume faces be computed and balanced with source terms within the control volume itself. There are four schemes in FLUENT for the calculation of face fluxes for the VOF model: geometric reconstruction, donoracceptor, Euler explicit, and implicit.
In the geometric reconstruction and donoracceptor schemes, FLUENT applies a special interpolation treatment to the cells that lie near the interface between two phases. Figure 22.2.1 shows an actual interface shape along with the interfaces assumed during computation by these two methods.
The Euler explicit scheme and the implicit scheme treat these cells with the same interpolation as the cells that are completely filled with one phase or the other (i.e., using the standard upwind, secondorder, or QUICK scheme), rather than applying a special treatment.
The Geometric Reconstruction Scheme
In the geometric reconstruction approach, the standard interpolation schemes that are used in FLUENT are used to obtain the face fluxes whenever a cell is completely filled with one phase or another. When the cell is near the interface between two phases, the geometric reconstruction scheme is used.
The geometric reconstruction scheme represents the interface between fluids using a piecewiselinear approach. In FLUENT this scheme is the most accurate and is applicable for general unstructured meshes. The geometric reconstruction scheme is generalized for unstructured meshes from the work of Youngs [ 300]. It assumes that the interface between two fluids has a linear slope within each cell, and uses this linear shape for calculation of the advection of fluid through the cell faces. (See Figure 22.2.1.)
The first step in this reconstruction scheme is calculating the position of the linear interface relative to the center of each partiallyfilled cell, based on information about the volume fraction and its derivatives in the cell. The second step is calculating the advecting amount of fluid through each face using the computed linear interface representation and information about the normal and tangential velocity distribution on the face. The third step is calculating the volume fraction in each cell using the balance of fluxes calculated during the previous step.
!! When the geometric reconstruction scheme is used, a timedependent solution must be computed. Also, if you are using a conformal grid (i.e., if the grid node locations are identical at the boundaries where two subdomains meet), you must ensure that there are no twosided (zerothickness) walls within the domain. If there are, you will need to slit them, as described in Section 5.7.8.
The DonorAcceptor Scheme
In the donoracceptor approach, the standard interpolation schemes that are used in FLUENT are used to obtain the face fluxes whenever a cell is completely filled with one phase or another. When the cell is near the interface between two phases, a ``donoracceptor'' scheme is used to determine the amount of fluid advected through the face [ 104]. This scheme identifies one cell as a donor of an amount of fluid from one phase and another (neighbor) cell as the acceptor of that same amount of fluid, and is used to prevent numerical diffusion at the interface. The amount of fluid from one phase that can be convected across a cell boundary is limited by the minimum of two values: the filled volume in the donor cell or the free volume in the acceptor cell.
The orientation of the interface is also used in determining the face fluxes. The interface orientation is either horizontal or vertical, depending on the direction of the volume fraction gradient of the th phase within the cell, and that of the neighbor cell that shares the face in question. Depending on the interface's orientation as well as its motion, flux values are obtained by pure upwinding, pure downwinding, or some combination of the two.
!! When the donoracceptor scheme is used, a timedependent solution must be computed. Also, the donoracceptor scheme can be used only with quadrilateral or hexahedral meshes. In addition, if you are using a conformal grid (i.e., if the grid node locations are identical at the boundaries where two subdomains meet), you must ensure that there are no twosided (zerothickness) walls within the domain. If there are, you will need to slit them, as described in Section 5.7.8.
The Euler Explicit Scheme
In the Euler explicit approach, FLUENT's standard finitedifference interpolation schemes are applied to the volume fraction values that were computed at the previous time step.
(22.28) 
where  =  index for new (current) time step  
=  index for previous time step  
=  face value of the th volume fraction, computed from the first  
or secondorder upwind or QUICK scheme  
=  volume of cell  
=  volume flux through the face, based on normal velocity 
This formulation does not require iterative solution of the transport equation during each time step, as is needed for the implicit scheme.
!! When the Euler explicit scheme is used, a timedependent solution must be computed.
The Implicit Scheme
In the implicit interpolation method, FLUENT's standard finitedifference interpolation schemes are used to obtain the face fluxes for all cells, including those near the interface.
(22.29) 
Since this equation requires the volume fraction values at the current time step (rather than at the previous step, as for the Euler explicit scheme), a standard scalar transport equation is solved iteratively for each of the secondaryphase volume fractions at each time step.
The implicit scheme can be used for both timedependent and steadystate calculations. See Section 22.6.4 for details.