16.08.2013 Views

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

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.

CHAPTER 5. RESULTS 80<br />

sudden increase in k + occurs for y/δ ≥ 0.5 (Figures 5.23 & 5.24). When the<br />

pressure gradient diminishes and becomes adverse, the effect on velocity is<br />

expected <strong>to</strong> manifest itself first near the wall and then propagate outwards.<br />

This happens because the laminar layer near the wall contains less kinetic<br />

energy and therefore less momentum. All of the modelled results exhibit<br />

this preferential slowing of the fluid near the wall when the pressure gradient<br />

becomes adverse, but they do so later in phase and <strong>to</strong> a lesser extent than<br />

what is seen in the DNS data. This can be seen in Figures 5.25 & 5.26.<br />

It is interesting <strong>to</strong> compare the snapshots for the phase angles of π/2, 3π/4,<br />

and π found in Figures 5.29 & 5.30. <strong>The</strong> excessive flatness of the 〈U〉 +<br />

profiles is apparent in all of these snapshots, but particularly at the phase<br />

angle of π, where the flattening occurs at a lower value of y + for all the models<br />

than at the phase angle of π/2. Between these two, the snapshot at phase<br />

angle 3π/4 shows the subgrid k-ω model already flattening at lower y + , while<br />

the low-Reynolds-number k-ε model persists in flattening only at higher y + .<br />

<strong>The</strong> subgrid k-ε modelled results follow the low-Reynolds-number solution<br />

for lower y + , but the 〈U〉 + profile flattens at y/δ ≈ 0.2 earlier <strong>to</strong> join the<br />

subgrid k-ω results for higher y + . <strong>The</strong> point at which the subgrid solution<br />

flattens corresponds roughly <strong>to</strong> the outer extent of the subgrid. Thus, at this<br />

transitional snapshot in time, when the 〈U〉 + profile is undergoing a change<br />

of shape, the solutions appears <strong>to</strong> be particularly sensitive <strong>to</strong> the averaging<br />

procedure used <strong>to</strong> calculate main grid boundary condition from the subgrid<br />

solution.<br />

An interesting feature of the <strong>UMIST</strong>-N subgrid approach is highlighted by<br />

the k-ω prediction of k + at phase angles of π/2 and 3π/4. Here, the k +<br />

profiles predicted by the subgrid k-ω model exhibit a double peak (two local

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

Saved successfully!

Ooh no, something went wrong!