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.

2D � � �<br />

x<br />

��<br />

b 2D � � �<br />

y<br />

��<br />

b 2D � � �<br />

z<br />

��<br />

�<br />

��<br />

��� � � �<br />

V ��� ��<br />

V<br />

t<br />

�<br />

�x � �� �u��<br />

�� t dV<br />

�x ��<br />

�<br />

�� ��<br />

�� � ���<br />

�<br />

�� t<br />

cf ��<br />

�� �u��<br />

� �<br />

�� t e<br />

�x �� x � dS<br />

S<br />

� ��t �e �� t�s<br />

�� e<br />

�<br />

��<br />

��<br />

��<br />

�� �� ��<br />

����u�� ����u�� Se,x ��� t�<br />

��<br />

w<br />

�x<br />

���x<br />

�� s<br />

�<br />

– 4.21 –<br />

�<br />

�u��<br />

�� ��<br />

e<br />

�x �� x ��<br />

�� �� dV<br />

�<br />

�u��<br />

��<br />

S<br />

�x �� x��<br />

��<br />

�<br />

�� cf<br />

� � n<br />

Sw, x � �t ��<br />

w<br />

����u�� ����u�� ����u�� Ss,x ��� t�<br />

S<br />

��<br />

t<br />

t,x � ��t �<br />

�x<br />

��<br />

b<br />

�x<br />

���x<br />

�<br />

�y � �� �v��<br />

�� t<br />

�x ��<br />

�� t<br />

�<br />

��<br />

��� dV � � �<br />

V ��� ��<br />

V<br />

t<br />

�<br />

�� ��<br />

�� � ���<br />

�<br />

�� t<br />

cf ��<br />

�� �v��<br />

� �<br />

�� t e<br />

�x �� y � dS<br />

S<br />

� ��t �n ��t�b �� n<br />

�<br />

��<br />

��<br />

����v�� ����v�� Sn,y � �� t�<br />

��<br />

s<br />

�x<br />

���x<br />

�� b<br />

�<br />

����v�� ����v�� Sb,y � �� t�<br />

��<br />

e<br />

�x<br />

���x<br />

�<br />

�z � �� �w��<br />

�� t<br />

�x ��<br />

�<br />

��<br />

�� �� ��<br />

�<br />

�v��<br />

�� ��<br />

e<br />

�x �� y��<br />

�� �� dV<br />

�<br />

�v��<br />

��<br />

S<br />

�x��<br />

y��<br />

��<br />

�<br />

�� cf<br />

� � t<br />

Ss, y � �t ��<br />

s<br />

�<br />

� � s<br />

Se, y � �t ��<br />

e<br />

��<br />

�� �� ��<br />

�<br />

��<br />

��� dV � � �<br />

V<br />

��� ��<br />

V<br />

t<br />

�<br />

�� ��<br />

�� � ���<br />

�<br />

�� t<br />

cf ��<br />

�� �w��<br />

� �<br />

�� t e<br />

�x �� z � dS<br />

S<br />

� ��t �t �� t�w<br />

�� t<br />

�<br />

��<br />

��<br />

����w�� ����w�� St,z � ��t �<br />

��<br />

b<br />

�x<br />

���x<br />

� ���w<br />

��<br />

���x<br />

�<br />

� � n<br />

Sw, z � �t ��<br />

w<br />

����u�� Sn,x �<br />

���x<br />

�� n<br />

�<br />

�<br />

Sb,x ��<br />

b<br />

����v�� St,y �<br />

���x<br />

�� t<br />

�<br />

����v ��<br />

���x<br />

�<br />

�w��<br />

�� ��<br />

e<br />

�x �� z ��<br />

�� �� dV<br />

�<br />

�w��<br />

��<br />

S<br />

�x �� z ��<br />

��<br />

�<br />

�� cf<br />

� � e<br />

Sb, z � �t ��<br />

b<br />

�� n<br />

�<br />

�<br />

Sw,y ��<br />

w<br />

����w�� Se,z �<br />

���x<br />

�� e<br />

�<br />

����w ��<br />

����w ��<br />

Sn,z � ��t �<br />

��<br />

s<br />

�x<br />

���x<br />

�<br />

Ss,z ��<br />

s<br />

(4.49a)<br />

(4.49b)<br />

(4.49c)<br />

The last four terms on the right-hand side are due to the grid non-orthogonality; they<br />

vanish for orthogonal cells. The cell face values of the velocity gradients are obtained by<br />

linear interpolation (see Sect. 4.3.3 and Eq. 4.21).<br />

4.3.9 Assembly of the coefficients<br />

After evaluating all terms of the convection, diffusion, and sources over the entire<br />

computational domain and rearranging the coefficients, the discretized transport equation

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

Saved successfully!

Ooh no, something went wrong!