K = 4 



K = 3 



K-2 



K = l 



j=l j=2 j=3 j=4 j=5 



Figure 14. Position of variables. 



and g s the components of gravitational acceleration, and u* and w* 

 convective terms. The terms on the right side of equations (129) and 

 (130) are taken at time t. Then, the pressures are recomputed from 

 the equation 



J 3,k 



l/ Pj+l,fc + Pj-1, 

 B \ Ax 2 



k Pj,fe+i + Pj> 



Az' 



1 * R i,k) 



(131) 



where 



and 



B = 2 



\Ax 2 Az 2 / 



1 A J+l/2,fc ,7-1/2, fe J,k+l/2 3 ,k- 



J,k 



!* 



Ax 



Az 



,fe-l/2\ 



(132) 



(133) 



The convective terms are defined as shown by Fromm (1968) . The 



equation for u*. , ,„ , is 

 n 3+l/2,k 



j+i/2, k j+i/2, k-i 

 (3 - D 2 



- 1 r ~ i 



-r (U . , - U . , ) 



2 J+l/2,k-2 J+l/2,k J 



C«J + i/ 2f fc-2 " 2 Vl/ 2 ,fe-l + Vl/ 2 ,fe } U34) 



62 



