The UMIST-N Near-Wall Treatment Applied to Periodic Channel Flow
The UMIST-N Near-Wall Treatment Applied to Periodic Channel Flow
The UMIST-N Near-Wall Treatment Applied to Periodic Channel Flow
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 4. NUMERICAL IMPLEMENTATION 51<br />
Before performing any dimensional analysis on the above equations, one must<br />
actually integrate with respect <strong>to</strong> volume, applying limits of 0 <strong>to</strong> 1 on the in-<br />
tegrations in x and z. It is also notable that the PASSABLE code uses a grid<br />
that is unit length in z but not unit length in x (∆x = 1). This detail must<br />
be remembered when working with the code, particularly when specifying<br />
volume-integrated source terms. Furthermore, in the PASSABLE implemen-<br />
tation, Equations 4.1-4.4 are multiplied through by ρ. This convention is<br />
maintained in the subgrid.<br />
4.3 Discretisation<br />
<strong>The</strong> process of discretisation involves approximating the governing differen-<br />
tial equations by algebraic equations so that they can be solved by a com-<br />
puter.<br />
A fully implicit scheme is used in this work. This means that the previous<br />
time-step values of adjacent nodes are not used. In contrast, a fully explicit<br />
scheme would use the previous time-step values, but not the current time-<br />
step values of adjacent nodes. <strong>The</strong> fully implicit method ensures that the<br />
solution can converge at large time steps.<br />
Central differencing is used in this work <strong>to</strong> approximate derivatives. Other<br />
differencing schemes include PLDS 2 [51] and QUICK 3 [37]. <strong>The</strong>se alternative<br />
schemes involve higher-order approximations. <strong>The</strong> QUICK scheme is also<br />
upstream-biased, in the sense that it makes use of more information on the<br />
2 Power-Law Differencing Scheme<br />
3 Quadratic Upstream Interpolation for Convective Kinetics