Algorithm
A static meniscus shape is initially specified and then grid generation is conducted using the static meniscus. In subroutine INI, initial value for each variables are given. Derivative values are calculated in MAT. Momentum equations are calculated by being followed by pressure correction. Finally, energy equation is solved, based on calculated velocity values from the momentum equations. Outside governing equation subroutines, a new meniscus surface is calculated by using calculated velocity and temperature profiles. Revised grid generation is conducted to accommodate the newly calculated meniscus shape. This iterative calculations updating the surface shape is repeated until a specified convergence criterion is met.
