Abstract:
The research presented in this thesis is concerned with the development of
numerical techniques and mathematical models for non-Newtonian
uids
and two-phase
ows in pipes and channels.
Single phase, turbulent
ow calculations of non-Newtonian
uids were performed
initially. Based on the literature a revised approach to wall modelling
is proposed and implemented. The approach uses analytical and
experimental analyses of the turbulent boundary layer structure. A comparison
with the standard approach is presented.
The interaction between turbulence and non-Newtonian behaviour is studied
by examining the rate of strain induced by
uctuating components of
velocity. The statistical analysis of published DNS data is performed. Finally,
a model is proposed where the turbulent rate of strain is determined
from turbulence quantities used by the Reynolds-averaged Navier{Stokes
model and used in the calculation of molecular viscosity.
For two-phase
ow, the solution procedure using periodic boundary conditions
was developed under an assumption of a
at interface. The numerical
technique was veri ed by comparing to an analytical result obtained for
laminar
ow in a channel. An extension to three dimensional
ow is performed.
With periodic boundary conditions standard turbulence models are applied
to two-phase strati ed
ow. Several models and their corrections for twophase
ow are assessed and a new model is proposed. The numerical studies
were carried out primiarily in the open-source code OpenFOAM, but initial
attempts were made in commercial packages such as STAR-CD and FLUENT.
Experimental data collected from the literature are used to verify the
results showing good agreement in pressure drops and phase fractions.