37 research outputs found
Numerical modelling of the snow flow characteristics surrounding Sanae IV Research Station, Antarctica
Thesis (PhD)--University of Stellenbosch, 2004.ENGLISH ABSTRACT:This work is concerned with the numerical simulation of the aeolian snow transportation
process (drifting or wind blown snow) and especially the snow deposition and erosion
phenomenon (snow drift). The research work is interested in modelling the atmospheric
boundary layer wind flow and its associated snow drifting processes around threedimensional
obstacles by means of computational fluid dynamics (CFD).
A modelling method is required to predict and evaluate the snow drifting phenomenon
surrounding the SANAE IV research station in Antarctica. This station is of an elevated
design to ensure that wind blown snow may travel around the structure relatively undisturbed
and without deposition near the structure. This design is partly successful but localised drifts
are formed especially leeward of the interconnecting structures that join the main building
sections together.
The theoretical and numerical description to describe the turbulent transport of the two-phase
mixture of air and snow particles is investigated. This theory is subsequently employed to
describe the snow deposition and erosion process and two models are developed to determine
the deposition flux onto the snow surface. These models presented and discussed are a
threshold based approach and a conservative based approach. The first model is dependent on
a threshold shear velocity to determine the onset of either erosion or deposition. The second
model determines the deposition or erosion flux based on the conservation of the snow mass
transport in the near surface control volume. A numerical scheme that evaluates the snow
deposition flux at the surface and forces a temporal surface adaptation during the simulation is
established and implemented in a commercial CFD software code by means of user
subroutines.
Various test cases for which observed snow drift data are available are numerically modelled
to validate the snow drift schemes presented in this work. These tests include the wind driven
snow accumulation around a three-dimensional cube, around two adjacent three-dimensional
cubes and near a typical porous snow fence. The results indicate that both methods can predict
realistic snow drifts for a variety of wind flow conditions but also show that the conservative
approach is superior to the threshold based approach in describing the snow drift process
around obstacles. This model allows drifts to form not only in areas of low flow velocities but
also under high shear conditions. The theoretical investigation and the development and
validation of the conservatively based snow drift scheme shows that drift formation depends
strongly on the near surface flow divergence and secondary flow structures. To resolve the
snow drift formation under a variety of flow conditions a three-dimensional field solution is
required to determine velocity and snow concentration gradients and include the effects of
near surface convective and turbulent entrainment.
The model is applied to numerically simulate and predict snow drifting around the SANAE
IV base for a moderate as well as a high wind speed event. The predicted snow drift around
the base agrees favourably with the observed drifts at the station. Further numerical
simulations are carried out to evaluate the effects a few design modifications may have on the
snow deposition. These results suggest that a simple baffle plate installation near the bottom
of the interconnecting link structures may minimise the snow accumulation leeward of that
area.
This study shows that to achieve realistic numerical snow drift predictions around, on or near
obstacles, a conservative based snow drift scheme should be considered using some form of
temporal terrain adaptation strategy. Only then does one include a sufficient level of
important flow effects such as deposition along near surface boundaries of strong flow
divergence which plays as an important role as vertical settling and entrainment in
determining deposition rates.AFRIKAANSE OPSOMMING:Hierdie studie behels die numeriese simulasie van windgedrewe sneeubeweging asook die
daarmee gepaardgaande sneeu neerslag en erosie eienskappe. Die navorsing het verder
belang in die berekening van die atmosferiese grenslaag vloei en die simulasie van sneeu
neerslag naby drie-dimensionele strukture deur gebruik te maak van berekeningsvloeimeganika
(BVM).
‘n Berekeningsmetodiek is nodig om die eienskappe van die sneeu neerslag rondom die
SANAE IV navorsingsstasie in Antarktika te voorspel en te evalueer. Die bogrondse struktuur
is spesifiek so ontwerp om te verseker dat wind gedrewe sneeu hoofsaaklik onversteurd verby
die struktuur kan beweeg sonder neerslag teenaan die struktuur. Die ontwerp is grotendeels
suksesvol alhoewel sneeu neerslag wel lokaal plaasvind, wind af vanaf die aansluitings
strukture tussen die hoof geboue.
Die teoretiese en numeriese beskrywing van die twee-fase lug- en sneeumengsel beweging
word ondersoek en gebruik om die sneeu neerslag en erosie einskappe te beskryf. Twee
modelle wat hierdie verskynsel beskryf word beskryf en bespreek naamlik ‘n drumpel
gebaseerde benadering en ‘n konserwatief gebaseerde benadering. Die eerste model is
afhanklik van ‘n drumpel skuifsnelheid om die aanvang van of erosie of neerslag te bereken.
Die tweede model bereken die neerslag eerder gebaseer op die behoud van die sneeu massa
vloei in die kontrole volume naby aan die oppervlak. ‘n Numeriese metode is ontwikkel en
geimplementeer in ‘n kommersiële BVM sagteware pakket deur van gebruikerssubroetines
gebruik te maak. Die ontwikkelde kode evalueer die sneeu neerslag vloed by die oppervlak en
forseer ‘n tydafhanklike oppervlak aanpassing gedurende die simulasie.
Die sneeu neerslag metode wat beskryf word in hierdie studie word ge-evalueer teen verskeie
toetsgevalle waarvoor daar waargenome sneeu neerslag resultate beskikbaar is. Hierdie toetse
sluit in die wind gedrewe sneeu neerslag rondom ‘n drie-dimensionele kubus, rondom twee
naby geleë kubusse en naby ‘n tipiese poruese sneeu heining. Die resultate dui aan dat beide
die metodes realistiese sneeu neerslag voorspel vir verskeie wind toestande. Die studie wys
ook dat die konserwatief gebaseerde benadering vir die beskrywing van die sneeu neerslag
proses meer akkuraat is as die drumpel gebaseerde benadering aangesien die neerslagvoorspel kan word nie net alleenlik in gebiede met lae vloeisnelhede nie, maar ook in gebiede
waar hoë skuifsnelhede teenwoordig is. Die teoretiese ondersoek, ontwikkeling en toepassing
van die konserwatief gebaseerde model dui daarop dat die neerslag afhanklik is van die
divergensie van die vloeiveld asook van die sekondêre vloei patrone naby die oppervlak. Ten
einde die sneeu neerslag vir verskeie toestande op te los is dit nodig om snelheids- en
sneeukonsentrasie gradiënte te kan bereken in ‘n drie-dimensionele vloei veld om sodoende
die invloed van naby-oppervlak konveksie en turbulente verspreiding in ag te neem.
Die metode word toegepas deur die sneeu neerslag rondom die SANAE IV navorsingsstasie te
voorspel vir ‘n gematigde asook ‘n hoë wind snelheid toestand. Die sneeu neerslag
voorspelling stem gunstig ooreen met die waargenome neerslag by die struktuur. Verdere
numeriese simulasies is uitgevoer om die invloed van ontwerpsverandering op die neerslag te
evalueer. Uit hierdie resultate blyk dit dat ‘n eenvoudige plaat struktuur onder die
aansluitingsstrukture die sneeu neerslag wind af mag verminder.
Hierdie navorsingsstudie dui daarop dat ‘n tydafhanklike terrein aanpassing strategie saam
met die konserwatiewe neerslag model noodsaaklik is ten einde realistiese resultate te behaal
vir die sneeu opbou rondom of naby strukture. Sodoende word genoegsame vlakke van
belangrike vloei verskynsels, soos die invloed van vloei divergensie, in ag geneem wat net so
‘n belangrik rol in neerslag speel soos vertikale afsetting
Whole-island wind bifurcation and localized topographic steering: impacts on aeolian dune dynamics
Recommended from our members
Reversing transverse dunes: Modelling of airflow switching using 3D computational fluid dynamics
Airflow dynamics across dune surfaces are the primary agent of sediment transport and resulting dune migration movements. Using 3D computational fluid dynamic modelling, this study examined the behaviour of near surface airflow travelling over transverse (reversing) dunes on a beach system. Wind direction was modelled in two opposing directions (both perpendicular to dune crestline) to investigate surface alteration of flow on the dune topography. Surface shear stress, velocity streamlines and potential sediment flux were extracted from the modelling. The work shows that under SW winds the surface (under the configuration measured) underwent almost 10% more aeolian flux than with opposing NE winds of the same magnitude. Differences were also noted in the airflow behaviour with SW winds staying attached to the surface with less turbulence while NE winds had detached flow at dune crests with more localised turbulence. The work provides detailed insights into how 3D airflow behaviour is modified according to incident flow direction of reversing dune ridges and the resulting implications for their topographic modification. These dune types also provide interesting analogues for similarly scaled Transverse Aeolian Ridges found on Mars and the findings here provide important understanding of flow behaviour of such landforms and their potential movement
Seasonal variations in airflow over the Namib Dune, Gale Crater, Mars: Implications for dune dynamics
A finite volume method for the analysis of the thermo-flow field of a solar chimney collector
Thesis (MEng)--University of Stellenbosch, 2000.ENGLISH ABSTRACT: This study investigates the implementation of the finite volume numerical
method applicable to non-orthogonal control volumes and the application of
the method to calculate the thermo-flow field within the collector area of a
solar chimney power generating plant. The discretisation of the governing
equations for the transient, Newtonian, incompressible and turbulent fluid
flow, including heat transfer, is presented for a non-orthogonal coordinate
frame. The standard k - E turbulence model, modified to include rough
surfaces, is included and evaluated in the method.
An implicit solution procedure (SIP-semi implicit procedure) as an alternative
to a direct solution procedure for the calculation of the flow field on nonstaggered
grids is investigated, presented and evaluated in this study. The
Rhie and Chow interpolation practice was employed with the pressurecorrection
equation to eliminate the presence of pressure oscillations on nonstaggered
grids.
The computer code for the solution of the three-dimensional thermo-flow
fields is developed in FORTRAN 77. The code is evaluated against simple
test cases for which analytical and experimental results exist. It is also
applied to the analysis of the thermo-flow field of the air flow through a
radial solar collector.
KEYWORDS:
NUMERICAL METHOD, FINITE VOLUME, NON-ORTHOGONAL, k+-e
TURBULENCE MODEL, SIPAFRIKAANSE OPSOMMING: Die studie ondersoek die implementering van 'n eindige volume numeriese
metode van toepassing op nie-ortogonale kontrole volumes asook die
toepassing van die metode om die termo-vloei veld binne die kollekteerder
area van 'n sonskoorsteen krag aanleg te bereken. Die diskretisering van die
behoudsvergelykings vir die tyd-afhanlike, Newtonse, onsamedrukbare en
turbulente vloei, insluitende hitteoordrag, word beskryf vir 'n nie-ortogonale
koordinaatstelsel. Die standaard k - E turbulensiemodel, aangepas om growwe
oppervlakrandvoorwaardes te hanteer, is ingesluit en geevalueer in die studie.
'n Implisiete oplossings metode (SIP-semi implisiete prosedure) as alternatief
vir 'n direkte oplossingsmetode is ondersoek en geimplimenteer vir die
berekening van die vloeiveld met nie-verspringde roosters. 'n Rhie en Chow
interpolasie metode is gebruik tesame met die drukkorreksie-vergelyking ten
einde ossilasies in die drukveld in die nie-verspringde roosters te vermy.
Die rekenaarkode vir die oplossing van die drie dimensionele termo-vloeiveld
is ontwikkel in FORTRAN 77. Die kode is geevalueer teen eenvoudige
toetsprobleme waarvoor analitiese en eksperimentele resultate bestaan. Die
kode IS ook gebruik om die termo-vloeiveld binne 'n radiale son
kollekteerder te analiseer.
SLEUTELWOORDE:
NUMERIESE METODE, EINDIGE VOLUME, NIE-ORTOGONAAL, k - E
TURBULENSIE MODEL, SI