A modeling study is undertaken to better understand the physics of katabatic flows. This study is divided into three topics; a comparison between a large eddy simulation (LES) and a mesoscale model of katabatic flows, a sensitivity study of katabatic flows to various physical parameters, and an investigation into the effect of subgrid scale terrain features on katabatic flow models. In the first topic, a comparison between LES, and a mesoscale model, ARPS, of katabatic flows is made to better quantify the accuracy of subgrid parametenzation in ARPS. It is shown that, although the modeled flows agree on a number of parameters, the LES model produces a lower and faster jet than that of ARPS, and also cools more near the surface. The momentum budgets of the two models agree well with each other. The ARPS model has a higher amount of TKE than the LES model, due to an overproduction by shear in the ARPS subgrid parameterizations. The second portion of this thesis represents a sensitivity study of katabatic flows to various physical parameters. The depth and strength of katabatic flows are shown to vary with surface heat fluxes, slope angle, and ambient stratification. Katabatic flows are shown to grow in depth and magnitude as slope angle increases, due to an increase in entrainment of overlying ambient air. The ratio of advection to mixing is shown to collapse to a near universal value regardless of surface heat fluxes. With increasing ambient stratification, entrainment in katabatic flows becomes small and the momentum equation is reduced to a two-way balance between buoyancy and drag. In this case, the heat flux of entrained air into the katabatic flow approaches that of the surface cooling, and the flow ceases to grow in the down-slope direction. Finally, predictions for bulk velocity and buoyancy strength scales are developed as a function of slope angle and surface heat fluxes. The last portion of this study focuses on the effect of subgrid scale terrain features on katabatic flows. It is shown that in areas of inadequate terrain resolution, the effect of the terrain smoothing routine in ARPS is to increase the slope height in areas of concave mountains. The concept of energy conversion in katabatic flows is introduced, and it is shown that the effect of raising terrain is to assign parcels more buoyant potential energy than they would otherwise have, and thus over-predict the magnitude of katabatic flows. Finally, an investigation into the effect of changing upper slope angle on katabatic flows over combined slopes is made. It is concluded that a combined slope cannot be predicted using a linear combination of simple slopes, since the transition portion of the slope results in a turbulent hydraulic jump with enhanced mixing. The magnitude of mixing in the turbulent hydraulic jump in combined slopes is shown to depend on the difference between upper and lower slope angle. / Graduation date: 2004

No description available.

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.