where ~u and u* might be viewed as intermediate velocities occurring
between the time steps n and n þ 1, although they are, in fact,
non-divergence-free approximations to un þ 1. The discrete operators
Vh and Dh represent the usual centered-difference approximations
for the gradient and Laplace operators, respectively.
2.3. Forcing strategies
The forces applied to the computational node due to presence of
a solid body must be determined at the advanced time level before
the solution procedure presented in Eqs. (4)e(8) can be started.
These forcing terms of any node near the solid body can be