Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
– 4.28 –<br />
The pressure gradients encountered in the source terms are computed according to the<br />
following relations:<br />
��<br />
��<br />
��<br />
��<br />
�p �<br />
��<br />
�n ��<br />
e<br />
e<br />
� �p� �� �� ��<br />
�� ��<br />
�� �x �� ��<br />
�� e<br />
�p c �� ��<br />
�� �� �<br />
���n<br />
��<br />
�pc �� �� ��<br />
�� �� ��<br />
�� �x �� ��<br />
��<br />
��<br />
��<br />
�p c<br />
��<br />
�� ��<br />
��<br />
e<br />
e<br />
� �pc �� �� ��<br />
��<br />
��<br />
�x �� ��<br />
��<br />
�� e<br />
S e,x<br />
S e<br />
S e,x<br />
S e<br />
�xPE LPE ��<br />
���<br />
��<br />
�p� �� �� ��<br />
�� ��<br />
���y<br />
�� ��<br />
�� e<br />
��<br />
���<br />
��<br />
�pc �� �� ��<br />
�� �� ��<br />
�� �y �� ��<br />
e<br />
S e,y<br />
S e<br />
S e,y<br />
S e<br />
�� ���pc<br />
�� ��<br />
��<br />
��� ��<br />
��<br />
�y �� ��<br />
��<br />
�� e<br />
��<br />
���<br />
��<br />
�p� �� �� ��<br />
�� ��<br />
�� �z �� ��<br />
�yPE LPE �� e<br />
��<br />
���<br />
��<br />
�pc �� �� ��<br />
�� �� ��<br />
�� �z �� ��<br />
e<br />
S e,z<br />
S e<br />
S e,z<br />
S e<br />
�� ���pc<br />
�� ��<br />
��<br />
��� ��<br />
��<br />
�z �� ��<br />
��<br />
�� e<br />
��<br />
��<br />
��<br />
��<br />
��<br />
��<br />
�zPE LPE ��<br />
��<br />
��<br />
(4.69)<br />
in which the interpolation of the pressure gradients along the Cartesian coordinates is<br />
done using Eq. 4.21.<br />
4.4.2 Pressure correction procedure<br />
The source term, b p , has pressure correction terms contained in the discharge due to non-<br />
no<br />
orthogonality of the cells, qe . These terms are evaluated explicitly by a double-step<br />
pressure correction procedure as follows:<br />
� Solve Eq. 4.68 for pc no<br />
by neglecting the non-orthogonal terms, qcf � 0 , and correct the<br />
velocities and pressure according to Eqs. 4.55a,b.<br />
� Solve again Eq.4.4 with the non-orthogonal terms now available from the first step<br />
and correct once again the velocities and pressure.<br />
4.4.3 Under-relaxation factor and time step<br />
To avoid instability of the computation, it is a common practice to put an underrelaxation<br />
factor to the pressure correction, 0 � � p �1, in updating the pressure:<br />
p � p � � � p p c (4.70)<br />
As mentioned in Sect. 4.2, the time step plays also as an under-relaxation factor for<br />
steady flow cases. This type of application, that is using the transient equations to solve<br />
steady flows, is generally known as a pseudo-transient computation. In order to achieve<br />
the effects of under-relaxed iterative steady-state computations from a given initial field<br />
by means of a pseudo-transient computation starting from the same initial field, the timestep<br />
size is taken such that (Fletcher, 1997, p. 365):<br />
� p 1<br />
�<br />
1 � E�t with E�t � ˜ a P<br />
�t (4.71)<br />
V P