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.

F w<br />

��<br />

– 4.18 –<br />

n� n�1 ��<br />

�<br />

� D �w Sw �� ��1 ��<br />

�<br />

� D �w Sw �� ��1<br />

� ��max�qw<br />

,0��<br />

fw ���<br />

P � ���max�<br />

�qw,0��<br />

fw ���<br />

W<br />

��<br />

LPW �� ��<br />

LPW �� �������������������������� �� �������������������������� �� ��<br />

D �<br />

�fw �� �w Sw<br />

implicit<br />

�<br />

������<br />

�<br />

���n��<br />

w<br />

��<br />

��<br />

�<br />

��<br />

�� ��<br />

�� ��<br />

��<br />

��<br />

�� ��<br />

��<br />

��<br />

��<br />

�� ����<br />

�� ��<br />

w<br />

��<br />

��<br />

��<br />

�� ��<br />

��<br />

�� ���������� �� �� ���������� �� ��<br />

explicit<br />

(4.43b)<br />

Shorter notations is used to write the expression of the convective-diffusive flux, for<br />

n�n�1 C D ��1 C D ��1 D<br />

example Eq. 4.43a maybe rewritten as F<br />

�� e ��aP � aP��<br />

P ��a E � a E��<br />

E � bE , in<br />

such a way that by summing up the convective-diffusive fluxes across the six faces of<br />

cell P, one obtains:<br />

F n�n�1 � a E<br />

��<br />

C D ��1 C D ��1 C D ��1 C D<br />

� � aE��<br />

E ��aW � aW��<br />

W ��aN � aN��<br />

N ��aS � aS �� S<br />

C D ��1 C D ��1 C D ��1 D<br />

�a T � a T��<br />

T ��aB � aB��<br />

B ��aP � a P��<br />

P � b<br />

��1 �<br />

(4.44)<br />

where the coefficients indicate contribution of convective-diffusive terms from<br />

neighboring cells as presented in Table 4.4. The diffusive-correction terms, being<br />

evaluated explicitly, are known from the previous iteration and are included in the<br />

independent coefficient b D .<br />

Table 4.4 Coefficients of the discretized convective-diffusive equations.<br />

Convective terms Diffusive terms<br />

��<br />

a C �<br />

E � � max �qe ,0<br />

��<br />

a C<br />

W � �max �qw<br />

��<br />

a C<br />

N � �max �qn<br />

��<br />

a C<br />

S � � max �qs<br />

��<br />

a C<br />

T � � max �qt<br />

��<br />

a C<br />

B � � max �qb<br />

C C<br />

aP � �a nb<br />

�� nb�EWNSTB<br />

� � ��<br />

� �<br />

, 0�<br />

��<br />

� �<br />

,0�<br />

��<br />

� �<br />

,0�<br />

��<br />

� �<br />

,0�<br />

��<br />

� �<br />

,0�<br />

��<br />

a D D �<br />

E � �f �e Se LPE a D D �<br />

W � �f �w Sw L PW<br />

a D D �<br />

N � �f �n Sn LPN a D D �<br />

S � �f �s Ss LPS a D D �<br />

T � �f �t St LPT a D D �<br />

B � �f �b Sb LPB � � �<br />

�<br />

qcf D<br />

aP � �<br />

D<br />

�anb<br />

cf�ewnstb<br />

nb�EWNSTB<br />

—<br />

b D � � f D � � S<br />

cf�ewnstb<br />

��<br />

�<br />

� � cf<br />

��<br />

������<br />

��<br />

�����n<br />

��<br />

�<br />

�� cf<br />

� �� �� ��<br />

��<br />

��<br />

����<br />

�<br />

�� cf<br />

��<br />

��<br />

��<br />

��

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

Saved successfully!

Ooh no, something went wrong!