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.31 –<br />

The above relation holds as well for p � and p c values.<br />

The boundary node denoted by W (instead of w) upon which the inflow boundary values<br />

are defined allows the discretized equations formerly established for interior nodes to be<br />

applied to the boundary node W without the need to change the notation.<br />

Momentum and k-� equations<br />

Fig. 4.6 Inflow boundary.<br />

In forming the coefficients in Eq. 4.52, the following steps apply:<br />

� all variables at the inlet are given: � W � � in<br />

� evaluate the convective-diffusive terms as for normal interior cells:<br />

C D C D<br />

aW, a W,<br />

aP<br />

� � , �a W P�<br />

, b<br />

W D � � , b<br />

W 1D � � , b<br />

W 2D � �W � bring the contribution of node W to the source term: b � b � b D<br />

� set the coefficient at node W to zero: a W � 0<br />

Pressure and velocity corrections<br />

� � W � a W<br />

C D � � a W��<br />

W<br />

The contribution of the discharge across the west face, q w , to the continuity equation, Eq.<br />

4.66, is replaced by the imposed discharge, q in . In forming the coefficients in Eq. 4.68,<br />

the following steps apply:<br />

p<br />

� set the coefficient at node W to zero: aW � 0<br />

� set the contribution of the inflowing discharge to the source term: b p<br />

� � w � �q in<br />

For the velocity correction, Eq. 4.57, the pressure correction gradient, �p c �x i , is<br />

computed by the finite-volume technique, Eq. 4.20, which requires the value of<br />

p c � � � p<br />

w c<br />

� � W . This latter is obtained by Eq. 4.73.

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

Saved successfully!

Ooh no, something went wrong!