27.11.2012 Views

FLOW AROUND A CYLINDER - istiarto

FLOW AROUND A CYLINDER - istiarto

FLOW AROUND A CYLINDER - istiarto

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

– 4.27 –<br />

Next, we need to define the equation of the pressure correction, which is carried out by<br />

using the continuity equation. The discretized continuity equation can be obtained by<br />

writing the governing equation, Eq. 4.15, with � = 1, � = 0, and R = 0 that gives:<br />

� � cf<br />

� � � qcf � un,cf Scf cf� ewnstb cf�ewnstb<br />

V �d<br />

��<br />

��<br />

�� S � V �d<br />

S<br />

��<br />

S<br />

cf�ewnstb<br />

� �<br />

� � 0 (4.66)<br />

Using Eq. 4.65 with the development of the implicit term in Eq. 4.60, and inserting the<br />

result to the above continuity equation, one gets the pressure correction equation of the<br />

form:<br />

� 1 ��V<br />

��<br />

� ��a<br />

˜<br />

� 1 ��V<br />

��<br />

� ��a<br />

˜<br />

�<br />

Se c c<br />

�� �pE � pP��<br />

e LPE 1<br />

�<br />

Ss c c<br />

�� �p S � pP��<br />

s LPS 1<br />

�<br />

� q �<br />

cf�ewnstb<br />

where:<br />

q � � � � � u<br />

cf cf<br />

RC 1 ��V<br />

��<br />

qcf �<br />

� ��a<br />

˜<br />

q cf<br />

��V<br />

��<br />

��˜<br />

a<br />

��V<br />

��<br />

��a<br />

˜<br />

RC no � � � q<br />

cf cf � qcf � �<br />

Sw c c<br />

�� �pW � pP��<br />

w LPW 1<br />

�<br />

St c c<br />

�� �pT � pP��<br />

t LPT 1<br />

�<br />

� 0<br />

�<br />

�<br />

� �Scf,x ��v cf�Scf,y<br />

��w cf�Scf,z<br />

�� cf<br />

no 1 ��V<br />

��<br />

� �<br />

���˜<br />

a<br />

�� cf<br />

�p � ����<br />

��<br />

����<br />

��<br />

���n<br />

��<br />

��<br />

��<br />

�� ��<br />

�� ��<br />

��<br />

��<br />

��<br />

cf<br />

��<br />

�p c<br />

��<br />

�n ��<br />

cf<br />

� p nb<br />

� �<br />

� pP<br />

L nbP<br />

� �pc �� ��<br />

��<br />

����<br />

��<br />

�� cf<br />

After arranging the terms, one has:<br />

� �<br />

��<br />

��<br />

��<br />

�� Scf old<br />

��<br />

��<br />

��<br />

��<br />

S cf<br />

Sn c c<br />

�� �pN � pP�<br />

n L PN<br />

��V<br />

��<br />

��a<br />

˜<br />

��V<br />

��<br />

��˜<br />

a<br />

Sb c c<br />

�� �p B � pP�<br />

b LPB interpolated discharge<br />

correction due to cell - face interpolation<br />

correction due to non - orthogonal terms<br />

(4.67)<br />

p c p c<br />

aP pP � � anb pnb<br />

nb<br />

� b p (4.68)<br />

where:<br />

p 1<br />

� �<br />

��V<br />

���˜<br />

a<br />

��<br />

a E<br />

p 1<br />

aS � �<br />

��V<br />

��<br />

� ��a<br />

˜<br />

a P<br />

S e<br />

��<br />

e LPE S s<br />

��<br />

s LPS p 1<br />

, aW � �<br />

��V<br />

���˜<br />

a<br />

��<br />

p 1<br />

, aT � �<br />

��V<br />

��<br />

���˜<br />

a<br />

� , b<br />

nb<br />

p � � q �<br />

�<br />

cf<br />

p p<br />

� � anb<br />

S w<br />

��<br />

w L PW<br />

S t<br />

��<br />

t L PT<br />

p 1<br />

, aN � �<br />

��V<br />

� ��a ˜<br />

��<br />

p 1<br />

, aB � �<br />

��V<br />

��<br />

� ��a<br />

˜<br />

RC no � � � q<br />

cf cf � qcf � �.<br />

S n<br />

��<br />

n LPN S b<br />

��<br />

w LPB ,<br />

,

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!