Numerical Simulation of Snow Deposition Around living Snow Fences

In this study, computational fluid dynamics (CFD) was used to investigate the air flow around porous snow fences to gain insight into snow transport and deposition in the vicinity of fences. Numerical simulations were performed to validate the CFD approach using experimental data from a wind tunnel study. Subsequent simulations were used to test the use of a porosity model to represent fence geometry and determine the effect of fence spacing for fences comprised of multiple rows. The results demonstrate that CFD simulations can reproduce the aerodynamics around porous fences. Additionally, the flow field generated with a porosity model is in close agreement with that from a model with explicit representation of fence porosity. Simulations of fences comprised of two rows spaced at various distances demonstrate that when the row spacing is small the fence behaves as a single row

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

Wind Shielding Analysis for Cold Regions Using Experimental and Numerical Techniques

The thesis presents a systematic experimental and numerical study on the interactions among porous fence, airflow, and windblown snowdrifts, a knowledge that will contribute to optimize the performance of porous wind shielding system in Cold Regions. A comprehensive review of the concepts, theories, techniques, and key findings associated with the research work has been undertaken. The key technical parameters influencing fence performance have been systematically studied by means of wind tunnel experimental investigations and Computational Fluid Dynamics (CFD) simulations. The study has found that porosity is the most influential structural parameter affecting the performance of porous fences in many aspects. Fence height stands a significant positive position in terms of its performance. It was found that fence performance is not sensitive to the changes of approaching atmospheric airflow velocity. Nevertheless, a bottom gap can improve snow fence trap efficiency. All of those findings agree with most of the findings of other researchers, which affirms that the research methodology adopted in this research is sound. Physical experimental work was performed to assess the reliability and credibility of the numerical models. Those models have been intentionally simplified, which made them easier to construct and quicker to obtain numerical solutions at a lower computational cost. Furthermore, the numerical models demonstrate the level of competence acquired through this research that is implemented in the optimisation of fence design. Special attention has been paid to the issues where elaborate research work has not been systematically reached in the open literature, this includes areas such as the effects of arrangement of porous holes, fence surface shear, and directions of wind load with respect to the fence, etc. Correlation between the reattachment length, the shelter distance, and the creation and distribution of fence surface shear is reported, to the author’s knowledge, for the first time in the open literature. General guidelines for the design of shelters based on porous fences have been established through this study. For example, the desirable size of hole range should be identified beforehand, and porous holes with sharp angular corners should usually be avoided in the fence design. It is recommended to place the fence within an angle of 30° to the wind load, where the effective shelter distance can be estimated in a linearized equation, and the normal drag coefficient can be described as a function of cos2θ. Optimal design of the arrangement of porous holes will maximize the fence performance, especially when the close fence environment is of concern. Although the definition of fence effective zone is still vague in the research field, the key factors influencing the fence effective zone have been investigated by evaluating the reduction of wind velocity leeward of the fence in this thesis. It is found that the fence effective zone is not sensitive to the change of approaching airflow velocity, and that increasing fence height will increase the physical size of the fence effective zone, but not in a proportional manner. It is also concluded that fence effective zone will be significantly reduced when the non-normal wind load is inclined at an angle greater than 30° to the fence. The effective zone increases effectively when the fence porosity is optimal. In contrast to the majority of published research work, the transient snow transport model presented in this work considers the snow transport rate as a whole without distinguishing the rate in saltation and suspension layer. The numerical study indicated that the position of the snow crest is mainly determined by the fence height, while porosity and bottom gap mainly affect the downwind deposition length. The optimal porosity for snow fences is in the range of 0.4 to 0.5, which is greater than the one for wind fences, which lies in the range from 0.25 to 0.35. Two snow crests have been observed leeward the fence at the onset of snow deposition, when the fence was placed without a bottom gap to the snow ground. This finding has not been encountered in any of the reported research work. Wind tunnel simulations of snowdrift around the fences have marginally under-predicted the sizes of snow deposition. The numerical predictions were quantitatively and qualitatively in good agreement with the field observations. This incompetence of wind tunnel experiments on porous fences implies that numerical modelling can play a more important role in snow fence research

Shield for Sand: An Innovative Barrier for Windblown Sand Mitigation

Abstract: Background: Windblown sand mitigation for civil structures in arid environment is crucial. Indeed, the number of railways crossing deserts and arid lands is increasing. A number of sand mitigation measures already exist. Among them, sand barriers are particularly intended for line-like infrastructures. We reviewed patented sand barriers on the basis of their shape and porosity. Objective: A new solid barrier for windblown sand mitigation called Shield for Sand is presented. Shield for Sand has been designed with the aim of maximizing the sand trapping efficiency through an upper windward deflector and simplifying its maintenance by complaining to sand removal ma-chines. The development of Shield for Sand follows the path traced by the Technology Readiness Level scale. Methods: The preliminary design of Shield for Sand has been supported by computational simulations of the wind flow around the barrier. Then, Shield for Sand has been tested in a wind tunnel with drifting sand in order to assess its efficiency. Both computational and experimental approaches allow an increase of the Technology Readiness Level. Results: The reversed flow induced by Shield for Sand increases its sand accumulation potential with respect to similar existing sand mitigation measures, such as the straight vertical wall. The efficiency of Shield for Sand resulting from the wind tunnel test is very high and almost constant with increasing sand accumulation level. Conclusion: The Shield for Sand working principles and performances are confirmed excellent. Final full-scale in-situ experiments are necessary to test the barrier under real environmental operational conditions

Performance analysis of wind fence models when used for truck protection under crosswind through numerical modeling

This paper is focused on truck aerodynamic analysis under crosswind conditions by means of numerical modeling. The truck was located on the crest of an embankment during the study. In order to analyze the performance of three wind fence models, the truck's aerodynamic coefficients were obtained and compared in two different situations either with or without the wind fences installed. In addition, the effect of both height and porosity of wind fence models on the aerodynamic coefficients acting on truck with respect to separation distance between the truck and the wind fence, was analyzed. A finite volume (or computational fluid dynamic) code was used to carry out the numerical modeling. The Reynolds-averaged Navier?Stokes (RANS) equations along with the k?? SST turbulence model were used to predict the behavior of turbulent flow. With respect to the results, the influence of the distance on the rollover coefficient is soft for all height values studied except for the lowest value (1 m of fence height), where the maximum value of rollover coefficient was obtained for the truck position closer to the fence. Regarding fence porosity, its effect on rollover coefficient is stronger for truck positions on road closer to the wind fence model.This work was supported by the OASIS Research Project that was co-financed by CDTI (Spanish Science and Innovation Ministry) and developed with the Spanish companies: Iridium, OHL Concesiones, Abertis, Sice, Indra, Dragados, OHL, Geocisa, GMV, Asfaltos Augusta, Hidrofersa, Eipsa, PyG, CPS, AEC and Torre de Comares Arquitectos S.L. and 16 research centres. The authors would also like to thank the GICONSIME research group of the University of Oviedo (Spain) for their collaboration in this research

Wind-sand tunnel testing of surface-mounted obstacles: Similarity requirements and a case study on a Sand Mitigation Measure

Windblown sand flow interacts with a number of surface-mounted human-built obstacles. The wind-sand flow perturbation and resulting morphodynamic response of the sand bed cannot be assessed in analytical terms. Therefore, wind-sand tunnel studies around scale physical models are often carried out. They should be driven by physical similarity theory based on dimensionless numbers referred to the whole multiphase and multiscale flow. However, similarity requirements cannot be fully satisfied under typical testing conditions and attention should be paid on the extent of the similarity relaxation. In this study, the background of wind-sand tunnel testing of surface-mounted obstacles is recalled by reviewing wind tunnel setups and similarity requirements. Then, a wind-sand tunnel campaign on a Sand Mitigation Measure is described and critically discussed. The setup dimensionless numbers are compared with statistics on those of past studies. The inescapable relaxation of similarity requirements is motivated by the test goals. The time evolution towards in-equilibrium conditions of both sand bed morphodynamics and sand transport is discussed. Finally, the results of engineering interest are described: the Sand Mitigation Measure sand trapping performance is assessed in dimensionless terms through the measurements of the incoming and outgoing sand concentration in air

Numerical Simulation of Snow Deposition Around Living Snow Fences

In this study, computational fluid dynamics (CFD) was used to investigate the air flow around porous snow fences to gain insight into snow transport and deposition in the vicinity of fences. Numerical simulations were performed to validate the CFD approach using experimental data from a wind tunnel study. Subsequent simulations were used to test the use of a porosity model to represent fence geometry and determine the effect of fence spacing for fences comprised of multiple rows. The results demonstrate that CFD simulations can reproduce the aerodynamics around porous fences. Additionally, the flow field generated with a porosity model is in close agreement with that from a model with explicit representation of fence porosity. Simulations of fences comprised of two rows spaced at various distances demonstrate that when the row spacing is small the fence behaves as a single row

Aerodynamics study of the flowfield at a shelterbelt

Shelterbelts are used world-wide for such purposes as reduction of soil errosion, control of snow drift, and provision of an effective agrometeorological method of field microclimate management and yield enhancement. Whether performing a wind tunnel test, conducting a field observation, or implementing a numerical simulation to investigate shelterbelt effects, researchers are more interested in an optimum reduction in a thin air lasier near the ground on the leeside of the shelterbelt rather than total wind-speed reduction in the whole flowfield. The purpose of this study is to formulate a Navier-Stokes based scheme to simulate the turbulent aerodynamic characteristics of a shelterbelt. Qualitative results from field observation of a living-tree shelterbelt under real atmospheric flow conditions and a wind-tunnel flow visualization of scale-model fences were used to explore the fundamental phenomena of the shelterbelt flow to help in the numerical modeling. A modified higher-order numerical scheme using the Lagrange interpolation to represent the interface convection terms is developed and applied to better simulate the turbulent shelterbelt flowfield. It is shown that this new scheme not only can enhance accuracy during computation but also is capable of retaining the numerical stability and good convergence characteristics which are lost in most higher-order numerical schemes. The flow retardation and porosity of shelterbelts are modelled via momentum sources with the help of the aerodynamic parameters, normal pressure drag and skin friction drag. The results obtained from this newly developed numerical scheme show satisfactory agreement with both field experiments and other numerical simulations. In addition, this procedure offers a generalized technique for simulating more complicated shelterbelt configurations

Optimization of Snow Drifting Mitigation and Control Methods for Iowa Conditions

TR-626 1 80224 00Blowing and drifting of snow is a major concern for transportation efficiency and road safety in regions where their development is common. One common way to mitigate snow drift on roadways is to install plastic snow fences. Correct design of snow fences is critical for road safety and maintaining the roads open during winter in the US Midwest and other states affected by large snow events during the winter season and to maintain costs related to accumulation of snow on the roads and repair of roads to minimum levels. Of critical importance for road safety is the protection against snow drifting in regions with narrow rights of way, where standard fences cannot be deployed at the recommended distance from the road. Designing snow fences requires sound engineering judgment and a thorough evaluation of the potential for snow blowing and drifting at the construction site. The evaluation includes site-specific design parameters typically obtained with semi-empirical relations characterizing the local transport conditions. Among the critical parameters involved in fence design and assessment of their post-construction efficiency is the quantification of the snow accumulation at fence sites. The present study proposes a joint experimental and numerical approach to monitor snow deposits around snow fences, quantitatively estimate snow deposits in the field, asses the efficiency and improve the design of snow fences. Snow deposit profiles were mapped using GPS based real-time kinematic surveys (RTK) conducted at the monitored field site during and after snow storms. The monitored site allowed testing different snow fence designs under close to identical conditions over four winter seasons. The study also discusses the detailed monitoring system and analysis of weather forecast and meteorological conditions at the monitored sites. A main goal of the present study was to assess the performance of lightweight plastic snow fences with a lower porosity than the typical 50% porosity used in standard designs of such fences. The field data collected during the first winter was used to identify the best design for snow fences with a porosity of 50%. Flow fields obtained from numerical simulations showed that the fence design that worked the best during the first winter induced the formation of an elongated area of small velocity magnitude close to the ground. This information was used to identify other candidates for optimum design of fences with a lower porosity. Two of the designs with a fence porosity of 30% that were found to perform well based on results of numerical simulations were tested in the field during the second winter along with the best performing design for fences with a porosity of 50%. Field data showed that the length of the snow deposit away from the fence was reduced by about 30% for the two proposed lower-porosity (30%) fence designs compared to the best design identified for fences with a porosity of 50%. Moreover, one of the lower-porosity designs tested in the field showed no significant snow deposition within the bottom gap region beneath the fence. Thus, a major outcome of this study is to recommend using plastic snow fences with a porosity of 30%. It is expected that this lower-porosity design will continue to work well for even more severe snow events or for successive snow events occurring during the same winter. The approach advocated in the present study allowed making general recommendations for optimizing the design of lower-porosity plastic snow fences. This approach can be extended to improve the design of other types of snow fences. Some preliminary work for living snow fences is also discussed. Another major contribution of this study is to propose, develop protocols and test a novel technique based on close range photogrammetry (CRP) to quantify the snow deposits trapped snow fences. As image data can be acquired continuously, the time evolution of the volume of snow retained by a snow fence during a storm or during a whole winter season can, in principle, be obtained. Moreover, CRP is a non-intrusive method that eliminates the need to perform man-made measurements during the storms, which are difficult and sometimes dangerous to perform. Presently, there is lots of empiricism in the design of snow fences due to lack of data on fence storage capacity on how snow deposits change with the fence design and snow storm characteristics and in the estimation of the main parameters used by the state DOTs to design snow fences at a given site. The availability of such information from CRP measurements should provide critical data for the evaluation of the performance of a certain snow fence design that is tested by the IDOT. As part of the present study, the novel CRP method is tested at several sites. The present study also discusses some attempts and preliminary work to determine the snow relocation coefficient which is one of the main variables that has to be estimated by IDOT engineers when using the standard snow fence design software (Snow Drift Profiler, Tabler, 2006). Our analysis showed that standard empirical formulas did not produce reasonable values when applied at the Iowa test sites monitored as part of the present study and that simple methods to estimate this variable are not reliable. The present study makes recommendations for the development of a new methodology based on Large Scale Particle Image Velocimetry that can directly measure the snow drift fluxes and the amount of snow relocated by the fence