27.11.2012 Views

FLOW AROUND A CYLINDER - istiarto

FLOW AROUND A CYLINDER - istiarto

FLOW AROUND A CYLINDER - istiarto

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

– 4.16 –<br />

formulae which combine the convective and diffusive transports in a special way.<br />

Depending on the grid Peclet number Pe, being the ratio of the convective and diffusive<br />

conductance, Pe � q L �S where L is the nodal distance, either the upwind-scheme<br />

convection, central-difference diffusion, or combination of the two, is considered to<br />

transport any scalar quantity � across a cell face.<br />

The hybrid scheme (Spalding, 1972) uses the upwind scheme for large Peclet numbers<br />

(|Pe| ≥ 2) and central difference for small Peclet numbers (|Pe| < 2). According to this<br />

C D<br />

scheme, the total flux across the east face, Fe � Fe � Fe , is defined as follows:<br />

� for Pe e � q eL PE � eS e � �2 , only the convective transport is taken into account:<br />

n�n�1 � ��1<br />

��<br />

Fe � qe �E<br />

(4.35)<br />

� for �2 � Pe e � q eL PE � eS e � 0, a part of the diffusive transport is also taken into<br />

consideration:<br />

n�n�1 � ��1<br />

Fe � qe �E � 1 � 0.5Pee<br />

��<br />

��<br />

��<br />

��<br />

��<br />

��<br />

� � � e<br />

�<br />

Se ��1 ��1 � �� E � �P �� �e Se<br />

LPE �<br />

����<br />

��<br />

�<br />

���n��<br />

e<br />

��<br />

��<br />

�<br />

�� ��<br />

�� ��<br />

��<br />

�� ��<br />

��<br />

��<br />

��<br />

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

��<br />

��<br />

e��<br />

��<br />

��<br />

(4.36)<br />

� for 0 � Pe e � q eL PE � eS e � 2, a part of the diffusive transport is also taken into<br />

consideration:<br />

n�n�1 � ��1<br />

Fe � qe �P � 1 � 0.5Pee<br />

��<br />

��<br />

��<br />

��<br />

��<br />

��<br />

� � � e<br />

�<br />

Se ��1 ��1 � �� E � �P �� �e Se<br />

LPE �<br />

����<br />

��<br />

�<br />

���n<br />

��<br />

e<br />

��<br />

��<br />

�<br />

�� ��<br />

����<br />

��<br />

�� ��<br />

��<br />

��<br />

��<br />

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

��<br />

��<br />

e��<br />

��<br />

��<br />

� for Pe e � q eL PE � eS e � 2 , only the convective transport is taken into account:<br />

(4.37)<br />

n�n�1 � ��1<br />

��<br />

Fe � qe �P<br />

(4.38)<br />

The power-law scheme (Patankar, 1980, p. 90-91) sets the limiting value of Pe where the<br />

diffusion no longer affects the transport at Pe = 10, instead of Pe = 2 used in the hybrid<br />

scheme.<br />

� for Pe e � q eL PE � eS e � �10 :<br />

n�n�1 � ��1<br />

��<br />

Fe � qe �E<br />

(4.39)<br />

� for �10 � Pe e � q eL PE � eS e � 0:<br />

n�n�1 � ��1<br />

Fe � qe �E � 1 � 0.1Pee<br />

��<br />

��<br />

��<br />

��<br />

��<br />

��<br />

� � 5 � e<br />

�<br />

Se ��1 ��1 � �� E � �P �� �e Se<br />

LPE �<br />

������<br />

�<br />

���n��<br />

e<br />

��<br />

��<br />

�<br />

�� ��<br />

�� ��<br />

��<br />

�� ��<br />

��<br />

��<br />

�� (4.40)<br />

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

��<br />

��<br />

e��<br />

��<br />

��

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

Saved successfully!

Ooh no, something went wrong!