Life-Cycle Performance and Cost Analysis of Sand Mitigation Measures: Toward a Hybrid Experimental-Computational Approach
Publication: Journal of Structural Engineering
Volume 148, Issue 7
Abstract
Windblown sand action affects civil structures and infrastructures in sandy environments, such as deserts and coasts. The wind interacts with human built structures of any kind leading to harmful effects, endangering their serviceability and users’ safety. To counter it, a number of sand mitigation measures (SMM) have been proposed, primarily through the trial-and-error empirical approach. As such, innovative approaches to properly quantify windblown sand action and to design SMM are needed in the current state-of-art and practice. In this study, the authors propose a novel hybrid approach to derive the life-cycle performance of SMM based on the combination of reliable wind-sand tunnel tests and innovative wind-sand computational simulations. Wind-sand tunnel tests are carried out to characterize the incoming sand flux in open field conditions. In a hybrid approach perspective, wind-sand tunnel measurements allow to properly tuned cheaper wind-sand computational simulations of the full-scale SMM performance. A probabilistic approach for determining windblown sand action and frequencies of sand removal maintenance is applied to a case study on a desert railway. Finally, a life-cycle cost analysis is carried out to assess extra-costs and savings derived from the implementation of the SMM. The proposed approach paves the way toward a comprehensive hybrid approach to the performance assessment of SMM.
Introduction
Windblown sand hazard affects several civil structures and infrastructures in sandy desert and coastal environments, such as single buildings, farms, towns, solar plants, pipelines, industrial facilities, roads, and railways (Bruno et al. 2018b). On one hand, coastal regions are experiencing the increased frequency of windstorms induced by climate change, giving rise to sand transport events from sandy coasts to built-up areas. On the other hand, desert regions are increasingly hosting human activities and built structures. Within this framework, railway infrastructures are particularly affected given their specific sensitivity to windblown sand and the increasing number of projects currently ongoing and planned in sandy regions across North Africa, Middle East, and Southeast Asia. This allows civil and structural engineering to familiarize with emerging challenging design issues, analogously to other design problems born from other application fields but resulting in key issues in structural engineering (e.g., Ribeiro et al. 2011; Harris et al. 2015).
Windblown sand interacts with surface-mounted human built obstacles of any kind inducing sand erosion and sedimentation around them. From a structural design perspective, windblown sand effects have been categorized into sand limit states (SLS) (Raffaele and Bruno 2019). Sand ultimate limit states (SULS) are defined as the threshold performance level beyond which structures are no longer safe, while attaining sand serviceability limit states implies their loss of functionality. Some examples are shown in Fig. 1: (1) passive lateral sand pressure on a wall and vertical load on a gable roof undermining safety; (2) indoor sand infiltration precluding serviceability; (3) railway full sand coverage compromising train passengers’ safety; and (4) permanent rail deformation affecting serviceability, induced by increased stiffness and decreased damping of contaminated ballast bed.
To cope with the effects above, the demand for the development of innovative approaches to properly quantify the so-called windblown sand action (Raffaele and Bruno 2020) and for the design and performance assessment of design solutions to reduce windblown sand induced effects, i.e., sand mitigation measures (SMMs) (Bruno et al. 2018b), has gained momentum in the structural and wind engineering literature.
Windblown sand action has been recently defined as an environmental action in analogy to snow (Raffaele and Bruno 2019). Such an analogy results from their phenomenological and modeling resemblance, but also from the comparable detrimental effects they induce on engineering structures and infrastructures (O’Rourke et al. 2005; Tominaga 2018). On one hand, windblown sand action primarily translates into sand accumulation, giving rise to variable, fixed, static loads directly applied to the affected structure. Remarkably, the Algerian snow and wind code is the sole standard defining global and local distributed vertical sand loads for flat and multispan roofs and providing the sand zone mapping of the country (D.T.R. C 2-4.7 2013). On the other hand, windblown sand translates into indirect actions undermining the performance of the built structure/infrastructure and resulting in periodic maintenance operations.
Once more in analogy to windblown snow, several design solutions to mitigate windblown sand effects have been proposed so far (see e.g., Alghamdi and Al-Kahtani 2005; Sañudo-Fontaneda et al. 2011; Basnet et al. 2014). SMMs aim to prevent sand from reaching the protected structure/infrastructure. Most of them are located between the sand source and the protected structure, and they are intended to trap incoming sand by promoting sedimentation [Fig. 2(a)]. As a result, the rigorous design shall be performed by considering both the SMM sand-trapping performance and the environmental loads induced by wind and sand. Such a kind of SMMs usually translates into berms and ditches realized through earthworks [Fig. 2(b)], reinforced concrete porous barriers [Fig. 2(c)], solid barriers [Fig. 2(d)], or a combination of them (Bruno et al. 2018b).
However, with some remarkable exceptions, the rigorous assessment of windblown sand action, and the design and performance assessment of SMMs remain at their early stage in the structural engineering literature, whereas they are mostly based on trial-and-error approaches in the technical practice. According to the authors, this is due to the multidisciplinary and multiphysic nature of the phenomenon coupling fluid dynamics and aeolian processes. As a result, to fill this gap of knowledge, research should benefit from studies in disciplines adjacent and partially overlapping structural engineering (such as wind engineering, applied mathematics, and aeolian geomorphology), and from experimental and numerical approaches.
Physical experiments usually translate into wind tunnel scale tests. wind-sand tunnel (WST) tests have been carried out both on flat ground conditions and around surface-mounted obstacles since the nineteen sixties (White 1996). WST tests allow to reproduce and measure with high accuracy and in a controlled setup the spatial and temporal evolution of wind and sand state variables. Nevertheless, WST tests show some deficiencies related to: (1) the technical difficulty of measuring wind shear stress and sand flux close to the sand bed, and (2) the experimental distortion arising from similarity mismatching related to the impossibility of jointly satisfying all multiphase/multiscale geometric and kinematic similarity requirements when scale models are tested (Raffaele et al. 2021).
Numerical simulations of multiphase flows of relevance to structural engineering applications have dramatically increased in the last several decades, particularly as regards wind-driven rain and windblown snow (see e.g., Kubilay et al. 2013; Tominaga 2018). The numerical simulation of windblown sand flow, herein called erosion-transport-deposition (ETD) simulation, is primarily carried out through the resolution of Eulerian–Lagrangian or fully Eulerian models coupling wind-flow aerodynamics and aeolian processes accounting for the morphodynamic evolution of the sand bed, as reviewed in Lo Giudice et al. (2019). Among them, Eulerian ETD simulations are emerging because they adapt well to the engineering needs of modeling large-scale processes and cutting costs with respect to WST tests. However, they shall be carefully adopted only after their calibration on physical experiments.
In this study, the authors pave the way towards a novel hybrid approach (Meroney 2016) to derive the life-cycle performance (LCP) of SMMs based on the combination of highly reliable WST tests on flat ground conditions, and innovative ETD simulations of the full-scale SMM behavior to overcome the limitations of the standalone methodologies. WST tests are carried out on a flat sand bed to characterize the incoming sand flux in open field conditions. Within a hybrid approach perspective, WST measurements allow to tune cheaper ETD simulations. ETD simulations are carried out by adopting an Eulerian multiphase first order computational fluid dynamics (CFD) model coupling wind flow aerodynamics and sand erosion, transport, sedimentation, and avalanching (Lo Giudice and Preziosi 2020). LCP is assessed through ETD simulations by taking into account the progressive loss of performance of the SMM caused by the gradual accumulation of sand around it. Then, the probabilistic approach to assess windblown sand action and plan sand removal maintenance operations proposed in Raffaele and Bruno (2019) is applied. This allows accounting for multiple uncertainties in environmental in-field conditions, i.e., wind speed, wind direction and threshold velocity for sand erosion, and for their propagation to the resulting windblown sand action. Finally, a life-cycle cost analysis (LCCA) is carried out. The LCCA is increasingly being adopted in the structural engineering domain to evaluate the monetary impact in a performance-based engineering framework (see e.g., Ierimonti et al. 2018; Le and Caracoglia 2019). LCCA allows to assess extra-costs and savings derived from the adoption of the SMM with respect to the unmitigated design scenario. The technical feasibility of the approach is demonstrated by discussing its application to a topical case study dealing with an endangered desert railway.
The paper is organized into four further sections. First, the experimental-computational techniques adopted within the hybrid approach and the modeling framework to assess the life-cycle performance and cost are introduced. The tested case study is introduced, and WST and ETD setups are outlined. Finally, results are discussed and conclusions and perspectives are outlined.
Hybrid Approach
The adopted wind tunnel experimental facility, multiphase computational fluid dynamics model, and probabilistic framework to assess the life-cycle performance and cost analysis are introduced in the following.
Wind Tunnel Facility
The wind tunnel test is carried out in the wind tunnel L-1B of von Karman Institute for Fluid Dynamics. The facility is a closed-circuit wind tunnel with a test section length of about 20 m, and a cross section of height and width [Fig. 3(a)].
A low-roughness boundary layer is reproduced to simulate open terrain conditions typical for sand desert. To characterize the clean wind boundary layer, avoiding interference of scattered light from flying sand particles with measuring equipment, a flat wooden board of length , width , and thickness with sand grains glued on it was set-up in the rectangular test section [Fig. 3(b)]. A ramp with gentle slope approximately equal to 3° is installed to smooth the transition between wind tunnel floor and the wooden board. The downwind edge is located at the distance from the inlet of the test section. The reference wind velocity is measured just upwind of the ramp through Prandtl pitot tube at 0.57 m from the wind tunnel floor. Some initial exploratory tests with a sand bed have been performed to ascertain the threshold velocity , defined as the minimum value of the wind speed at which quasisteady sand transport occurs at position . Three increasing wind speeds, respectively, equal to , are tested.
The mean velocity profile and mean turbulence intensity profile at the distances and from the upwind edge of the wooden board are measured along the test section centerline through particle image velocimetry (PIV) technique adopting a smoke generator to seed the flow with oil particles ranging from . Each measurement is taken along the test section centerline to avoid the influence of the boundary layer developed on the lateral sides of the wind tunnel. A 200 mJ Nd:YAG laser source pulsating at 10 Hz is located far downwind of the measuring section. A laser sheet is generated along the test section centerline and perpendicular to the floor. Its width inevitably varies along the fetch length resulting equal to 8.6, 7.2, and 5.4 mm, respectively, at positions , , and . Two CMOS cameras with resolution pixels and a 50 mm objective are located outside the test section to acquire the wind flow field with field of views (FoVs) equal to . The measured profiles at and results are almost unchanged. Fig. 3(c) shows the measured profiles at . The wind speed measurements are fitted with the log-law , being the wind shear velocity, the von Karman constant, and the aerodynamic roughness, leading to and . The corresponding profile of the streamwise turbulence intensity are included in Fig. 3(c), whereas the related streamwise integral length scales measured at pitot location are .
Computational Model
The wind flow is modeled as an unsteady incompressible turbulent flow through Unsteady Reynolds Average Navier Stokes (URANS) equations. URANS equations are chosen because we are primarily interested in the long-term behavior of the sand transport which induces sand erosion and accumulation around SMMs and which takes place on a much larger time scale than turbulence. The SST turbulence model is adopted because of its proven accuracy in simulating wind flow separation around bluff bodies (Menter et al. 2003), such as SMMs. The same CFD model has been validated in Bruno and Fransos (2015) and adopted in Bruno et al. (2018a), Horvat et al. (2020), and Horvat et al. (2021) on the same class of problems, i.e., nominal 2D bluff bodies immersed in a turbulent atmospheric boundary layer. The set of governing equations reads:where = time; = air density; = average pressure; = air kinematic viscosity; = so-called turbulent kinematic viscosity; = turbulent kinetic energy; and = its specific dissipation rate. The kinetic energy production term is modeled by introducing a production limiter to prevent the build-up of turbulence in stagnation regions, i.e., , where and = strain-rate tensor. The standard blending function and model constants , , , , and are obtained from Menter et al. (2003).
(1)
Sand-grain roughness wall functions are complemented by SST turbulence model because of their wide use in computational wind engineering (Blocken et al. 2007) and their proven adequacy from past simulations on the same class of problem (e.g., Liu et al. 2011; Bruno and Fransos 2015). Standard wall functions (Launder and Spalding 1974) with roughness modification (Cebeci and Bradshaw 1977) are applied. The equivalent sand-grain roughness height is determined equal to , where is the roughness constant.
The sand phase is considered as a passive scalar and modeled through the conservation equation of sand volume fraction :where the sand flux is given by the combination of advection by wind, sedimentation effects due to gravity, and diffusive flux. In particular:where = transport velocity of the sand particles by wind taken proportional to , = vertical sedimentation velocity, and takes into account the mixing-diffusive contribution resulting from the combination of and the viscous effect due to random collisions at the sand surface , being and model parameters to calibrate the concentration profile, and the second invariant of (Preziosi et al. 2015).
(2)
(3)
The morphodynamic evolution of the sand surface is accounted for by imposing the continuity of sand flux through the wind-sand interface and the triggering of sand sliding when the sand slope exceeds the critical angle of repose , i.e., the steepest slope angle sand grains can pile up. This leads to the modified Exner equation (Lo Giudice and Preziosi 2020):where = height of the wind-sand interface defined on sandy patches only; = diffusion coefficient controlling the sliding speed; stands for positive part, = velocity of the wind-sand interface due to erosion or sedimentation; = sand close-packing volume ratio; and = direction normal to the surface.
(4)
Within the framework of this study: (1) steady-state simulations have been carried out on the WST setup to tune the model free parameters, and (2) unsteady simulations have been carried out on the full-scale SMM to assess its LCP. The adopted boundary conditions (b.c.) and 2D computational domains are schematically shown in Fig. 4. Null initial conditions are imposed in the whole domains. No-slip b.c. are imposed at the solid walls. At the inlet, Neumann zero-gradient b.c. is imposed for and , whereas Dirichlet b.c. is imposed on , , . At the outlet, Neumann zero-gradient b.c. is imposed for all flow variables, except for Dirichlet b.c. for . Concerning the WST domain, the inlet profile of is prescribed using the power law for turbulent boundary layer, whereas the inlet profile of and are set constant and equal to and . Concerning the full-scale SMM domain, a log-law inlet profile is set, whereas the profiles of and are set in accordance to Richards and Norris (2011) to replicate a neutral atmospheric boundary layer, and symmetry b.c. is imposed at the top. Finally, a Dirichlet sand erosion b.c. is set to properly model erosion on sandy surfaces, whereas a Neumann b.c. is imposed on nonerodible surfaces. In particular, erosion occurs when the wind shear velocity is higher than the threshold one (Kok et al. 2012). According to Ho et al. (2011):where = model parameter depending on the physical properties of the sand; = mean sand diameter; = acceleration due to gravity, and directly follows from and takes into account the local effect due to inclined slopes; and being the ground slope angle (Iversen and Rasmussen 1994).
(5)
Space discretization follows a predominantly structured grid of hexahedral control volumes. The mesh is refined close to the ground, so that the height of the wall-adjacent cell: (1) provides a sufficiently high mesh resolution along the normal direction to the surface to adequately resolve the gradients of wind-sand state variables, and (2) complies with the wall function requirement on dimensionless wall unit , being the cell center height. The total number of cells for WST and full-scale SMM domains is equal to about and , respectively.
The cell center values of the state variables are interpolated at face locations using the second-order central difference scheme for the diffusive terms. The convection terms are discretized by means of the so-called limited linear scheme. SIMPLE and PIMPLE algorithms are used for pressure-velocity coupling for steady-state and unsteady cases, respectively. The finite volume open source code OpenFOAM© is adopted. On average, a CPU time of about 1.5 h is required for each steady-state simulation and a CPU time from 72 to 96 h is required for each unsteady simulation on Intel(R) Core(TM) dual-processor Sandy Bridge server @ 2.6 GHz with 16 cores employed.
Life-Cycle Performance and Cost Analysis
The modeling framework to assess windblown sand action on a generic surface-mounted obstacle is shown in Fig. 5(a). The sand transport under open field conditions is commonly quantified by the incoming sand transport rate , defined as the integral over the vertical direction of the sand flux . In the present study, is quantified by taking into account the probabilistic description of and (Raffaele and Bruno 2019) through the semiempirical expression:where is a reference sand grain diameter. Whenever encounters an obstacle, it splits into sedimentation rate and outgoing transport rate . In turn, the windblown sand action results from the time-cumulated and translates into the accumulated sand volume .
(6)
Within this framework, the life-cycle performance of the th obstacle is expressed by the sand trapping efficiency . which depends on the sand transport direction and on the overall obstacle shape, that depends in turn on the initial geometry and on the volume of sand sedimented around it, . Therefore, the resulting accumulated sand volume after time is given by the random sum:where = independent and identically distributed copies of ; = sampling interval of the wind speed; and = number of sand transport occurrences along the direction over time . The basic condition for a satisfactory state is given by , being the resistant sand volume. From this, the probability of failure results equal to:where and = probability density and cumulative distribution functions of , respectively. As a result, the characteristic time of failure is given by:where = characteristic probability of failure. Fig. 5(b) sketches the generic trend of the time-variant sand action through its density , the mean value , and the increasing trend of the probability of failure .
(7)
(8)
(9)
Finally, LCCA (Fabrycky and Blanchard 1991) can be carried out based on , considering the cost of SMM design and construction , and sand maintenance related costs associated with the th obstacle :where = cumulated life-cycle cost at time ; = time required for SMM design and construction; = the number of infrastructure components inducing ; and = the discount rate.
(10)
Setup of the Study
Study Layout
The proposed hybrid approach is applied to a case study railway segment located near Al Ain along the Ethiad Rail line in the United Arab Emirates. The chosen site is threatened by the sand of Rub’ al-Khali desert. The examined railway segment develops along the NE-SW direction. The local sand granulometry is composed of fine grained, moderately well-sorted sand with mean diameter (Edgell 2006). The probability density function of is derived from Raffaele et al. (2016) as a function of . Conversely, the probability density functions of is derived from an anemometric station in proximity of Al Ain. The employed dataset of the 10 min average wind velocity at 10 m from the ground refers to 10 years, from January 2008 to December 2017. The time sampling corresponds to 1 h, whereas the yaw angle discretization is equal to 10°. The aerodynamic roughness length is set equal to , according to the recommendations given in EN 1991-1-4 (CEN 2005).
Fig. 6 shows the windblown sand action modeling framework and the related state variables for a single side of the railway infrastructure. The same scheme is mirrored for the opposite side. Such a scheme directly results from the general one in Fig. 5, by putting in series three successive obstacles, i.e., the SMM (), the embankment (), and track (). The berm and ditch is a well-known SMM for desert railway applications [see e.g., Phillips (2011)]. Here, we investigate a generalization of it. The tested SMM results from the combination of a high long berm, followed by a deep long ditch. The berm follows the sinusoidal profile . The ditch side walls are inclined with 1/3 slope gradient. The height of the embankment is set equal to . A double-track railway is considered, with a deep ballast bed and a 5 m wide line-side access track.
The LCP of the SMM is obtained through the proposed hybrid approach. The LCP of the embankment is derived from WST tests in Hotta and Horikawa (1990), whereas the LCP of the railway track is conjectured maximum and constant up to the filling of ballast voids and then linearly decreasing until it reaches its full capacity , in analogy to Raffaele and Bruno (2020). and are set equal to the volume corresponding to 80% and 50% of the initial value of , respectively. Conversely, is set equal to the volume giving rise to the SULS at full rail coverage.
The LCCA spans the whole service life of the infrastructure, set equal to . The railway construction time is set equal to . The costs retained in the LCCA comprehend: (1) SMM design and construction , (2) removal and disposal of accumulated sand around the SMM and embankment , (3) railway ballast cleaning . For the sake of generality, the discount rate is set equal to , whereas costs induced by land expropriation due to SMM construction and by loss of capacity due to sand maintenance operations are not considered in the present study. Indeed, they go beyond the single analyzed railway segment and highly depend on the specific features of the whole railway network and country economic system.
Wind-Sand Tunnel Test Setup
A flat sand bed of length , width , and thickness has been set up in the wind tunnel test section replacing the wooden board [Fig. 7(a)]. The sand bed is confined by wooden slats and fillet to the wind tunnel floor through the upwind ramp in analogy with the setup proposed by Tominaga et al. (2018).
The particle size distribution of the tested sand has been obtained through microscopic imaging technique and is plotted in Fig. 7(b) through its cumulative distribution . The mean grain diameter is equal to , very close to the one on site. The threshold shear velocity obtained through preliminary wind tunnel tests results equals to (corresponding to ). Such a value of agrees with past wind tunnel measurements on sand samples with same (see e.g., Raffaele et al. 2016).
Three tests have been performed by progressively increasing the incoming wind speed . Each test started from the same initial condition, a perfectly flat uniform sand bed, and lasted to ensure reaching a quasisteady sand transport. To ensure repeatability, each test has been performed twice. The similarity of the saltation layer with respect to open field condition is taken into account by: (1) assuring that the wind flow is fully rough, i.e., the friction Reynolds number satisfies the criterion (Anno 1984), (2) avoiding disturbance in streamwise pressure by satisfying the Froude number criterion, i.e., (White 1996), and (3) adopting a long sand fetch equal to to ease sand transport saturation and let the wind flow adjust to the sand (Kok et al. 2012).
The installed measuring equipment allows for the detection of the instantaneous sand particles concentration and velocity over the flat sand bed through particle tracking velocimetry (PTV). The laser source pulsates at 2 Hz. The two CMOS cameras are moved along the sand fetch to acquire PTV pictures in , , and . The acquired FoVs are sketched in Fig. 7(c). The velocity and concentration fields of sand particles are detected with FoVs for each position, with a 35 mm objective.
Wind-Sand Computational Simulations: Incoming Wind Flow and Sand Features
The incoming wind flow features differ as a function of the computational domain. On the one hand, the incoming wind flow of the WST scale simulation reflects the experimental conditions. The ground aerodynamic roughness length is set equal to , according to the boundary layer PIV measurements. The free parameter of the inlet wind speed power law is set to reproduce the wind speed profiles measured at and , resulting in . On the other hand, the incoming wind flow of the SMM full-scale simulations reflects actual desert conditions. The ground aerodynamic roughness length is set equal to . The incoming wind shear velocity is set equal to to exceed and induce windblown sand transport upwind from the SMM. During sand transport events, the lower bound () of the reference wind speed at the top of the berm is equal to , and the corresponding Reynolds number is . In the case of nominally 2D sinusoidal berms, such a value suggests that the flow triggering windblown sand transport is predominantly within the Reynolds supercritical regime (Ferreira et al. 1995). As such, significant Reynolds effects are not expected to take place for .
For both computational domains, the adopted sand diameter is set equal to the mean value of the tested sand, i.e., . The corresponding mean values of threshold shear velocity and sedimentation velocity are respectively set equal to , in accordance to wind tunnel PIV measurements, and to , according to the statistical characterization proposed in Raffaele et al. (2020). Finally, the remaining model constants are set equal to , , and .
The adopted time step in unsteady ETD simulations is set equal to , giving rise to a maximum Courant number . Given the time demanding unsteady multiphase simulations, the whole problem is decomposed into shorter simulations to assess the discrete piece-wise LCP. For each tested sand level: (1) an unsteady ETD simulation is carried out up to reaching of quasisteady sand transport conditions by assuring the convergence of incoming and outgoing sand transport rates, (2) LCP is assessed, (3) a new geometry corresponding to a new extrapolated sand level is obtained by shifting the wind-sand interface proportionally to the local value of the erosion-deposition velocity , and (4) a new mesh is built on the obtained geometry. The above steps are repeated systematically to describe the trend of LCP, up to reaching null SMM performance, i.e., sand flux at inflow equal to sand flux at outflow. The shifting of the wind-sand interface is done manually relying on , being an arbitrarily chosen time interval. As such, the convergence of sand transport shall be assured for each extrapolated sand level.
Results
In the following, the tuned steady ETD simulations are compared with WST measurements, then SMM life-cycle performance is assessed through unsteady ETD simulations. Finally, the time-varying windblown sand action is quantified on the considered case study and life-cycle cost analysis is applied. The proposed hybrid approach can be applied to any kind of structure/infrastructure as long as the LCP of any affected obstacle is available. At the time being, scarce laboratory measurements and even more scarce computational simulations only allow for the quantification of the LCP in the railway infrastructure domain.
Comparison between WST Measurements and ETD Results
Steady-state ETD simulations are tuned on WST measurements by setting the value of free parameters in Eqs. (3) and (5). In this section, the comparison between time-averaged WST measurements and ETD results is provided. In particular, Fig. 8 collects the vertical profiles of the sand volume fraction , streamwise component of the sand velocity , and flux at positions (a,d,g), (b,e,h), and (c,f,i), providing the mean value and standard deviation .
The sand volume fraction follows the typical exponentially decreasing law (Ho et al. 2011). For a given position, the rate of decay of is almost the same for and , whereas it sensibly varies for . Conversely, for a given wind speed, increases streamwise. According to the authors, the lower than expected and high in for are caused by the experienced intermittent unsteady sand transport induced by still close to the threshold and the short fetch length. Indeed, the adopted first order ETD model is not able to capture such out-of-equilibrium erosion conditions. Discrepancies in between WST and ETD fade out with the increasing of sand fetch.
The sand velocity profiles can be divided into two layers resulting from strong variation of (Valance et al. 2015). Within the near-wall region where is high (), is weakly sensitive to and linearly increasing with . Conversely, where is low (), is highly sensitive to and follows a logarithmic-like increasing trend. Despite the overall promising matching of profiles for a first order model, ETD simulations are not able to reproduce at the ground, and particularly overestimates close to the ground () for . The ETD first order model suits well highly diluted flows (i.e., far from the sand bed), whereas it shows some deficiency close to the wall where sand concentration is high (Lo Giudice et al. 2019).
The longitudinal sand flux directly results from , being the sand density. Whenever measurements were missing close to the wall due to nonuniform laser reflections, is extrapolated from fitted exponentially decreasing functions [see blank circles in Figs. 8(g–i)]. For , the measured profiles follow a typical decreasing exponential trend. Conversely, for the trend is no more monotonic and a maximum appears at about , analogously to e.g., Raffaele et al. (2021). The discrepancies in and between WST and ETD propagate to giving rise to: (1) the strong mismatching in for , (2) the nonnegligible overestimation in correspondence of the maximum value of for , and (3) the nonphysical null values of at the wall.
Finally, sand transport rate is evaluated to ascertain discrepancies between WST and ETD in bulk terms. The convergence in time is quantitatively checked by means of the weighted residual of the mean value of sand transport rate . The weighted residual is defined for growing dimensionless time as . The obtained results are plotted in Fig. 9(a) for , for the sake of conciseness. For each position, the trend of convergence gives rise to at about , consistently with Raffaele et al. (2021). Pointwise measurements and ETD results are compared in Fig. 9(b). For , ETD results overall quantitatively agree well with WST measurements (apart from ). For , the mismatching between WST and ETD is significant. In particular, ETD overestimates WST results. WST discrepancies may be related to the longer fetch required to reach equilibrium conditions due to the higher , as highlighted in Dong et al. (2004).
SMM Life-Cycle Performance
The wind flow features around the SMM are characterized though a steady RANS CFD simulation. Fig. 10(a) shows the wind flow topology around the clean SMM by means of streamlines coupled with a vorticity magnitude field. The vorticity magnitude is high close to the wall due to the atmospheric boundary layer, whereas the outer free flow is quasiirrotational. The boundary layer separates on the berm downwind slope, inducing a large clockwise vortex. The recirculation area is characterized by the coalescence of two smaller vortices, respectively, induced by the berm and by the ditch. The wind flow then reattaches within the ditch and accelerates along its downwind slope. The vertical dimensionless wind speed profiles upwind (p1), across (p2, p3, p4), and downwind (p5) from the SMM are plotted in Fig. 10(b) versus the dimensionless height , being the ground height, to further examine the wind flow pattern. p1 shows the upwards flow deflection induced by the berm, that results in a decrease of and . p2 shows a significant flow acceleration induced by the berm upwind slope resulting in an increment of close to the ground and . p3 and p4 are located within the recirculation region. p3 is located just downwind the berm, showing a large negative . p4 is located within the ditch and closer to the reattachment point, showing an overall lower wind speed magnitude and compared with p3. p5 is located downwind the ditch showing progressive wind speed recovery and a weak .
To assess the full-scale SMM life-cycle performance , unsteady ETD simulations are carried out starting from the above-mentioned wind flow initial condition. Nine consecutive ETD simulations have been carried out to describe the morphodynamic evolution of the sand bed, starting from the clean SMM (fill 0) up to its complete filling (fill 8). Fig. 10(c) shows the multiphase flow topology around the SMM by means of wind flow streamlines coupled with a dimensionless sand concentration field, for the increasing filling levels 0, 4, 6, and 7. Sand is eroded from the upwind side of the SMM and gradually transported downwind. Sand starts depositing along the berm upwind foot and upwind slope due to the upwards flow deflection highlighted in p1 [Fig. 10(b)]. After that, sand climbs over the berm due to wind flow acceleration and progressively deposits along the berm downwind slope and along the ditch upwind cut because of wind flow separation [see p2, p3, and p4 in Fig. 10(b)]. The gradual accumulation of sand within the SMM body progressively moves the separation and reattachment points downwind and upwind, respectively. This results in the shrinking of the clockwise vortex up to its vanishing. The adjustment of the wind flow features around the SMM modifies in turn and, therefore, erosion and deposition zones along the SMM profile. Accordingly, sand is transported up to the ditch downwind slope, slightly depositing on it, and eventually escaping from the SMM.
To assess the reaching of equilibrium conditions, the convergence of incoming and outgoing is visually checked for growing dimensionless time units . The dimensionless ratio sampled at every second is plotted in Fig. 10(d) for each performed ETD simulation. In particular, is evaluated in open field, and is evaluated just downwind the ditch. Overall, monotonically increases before reaching a plateau. The dimensionless time at which convergence is reached decreases with increasing filling level. This is primarily due to sedimented sand closer to the position at which is evaluated for high filling levels. Once quasisteady sand transport is reached, LCP can be evaluated. The resulting life-cycle performance is plotted in Fig. 10(e) for increasing filling ratios . The trend of reflects the wind flow pattern variation induced by sand morphodynamics. is initially constant and maximum for , i.e., the majority of sand is accumulated on the upwind side of the berm. After that, starts decreasing when sand consistently reaches the berm downwind slope and the vortex size decreases. As soon as the recirculation region disappears and the sand reaches the ditch downwind slope, drops to zero. The SMM resistant sand volume is then identified as the accumulated volume for which the efficiency drops to 80% of the initial value, i.e., .
Windblown Sand Action and SMM Life-Cycle Cost
Windblown sand transport attacks the railway line from the SE and NW sides. Figs. 11(a–d) provide the detailed description of incoming wind through the wind rose of : (1) the wind speed time histories, (2) together with their probability density functions , (3) direction frequencies , and (4) by referring to SE and NW sides. The Weibull-shaped are almost equivalent on both sides of the alignment, whereas wind occurrences are higher on the NW side. Because of that, the NW side is expected to be the most sensitive to windblown sand action. Fig. 11(e) shows the probability density function of the threshold shear velocity for sand transport .
The time growing windblown sand action on each infrastructure component is assessed by means of Monte Carlo simulations based on the bootstrapping technique, accounting for about numerical realizations. The resulting probabilistic windblown sand action on the NW side of the track is plotted in Figs. 11(f–h) through the mean value , the 5th percentile , and 95th percentile moved, and some probability density functions . The time-constant is plotted as well to extrapolate the time increasing failure probability. For each component, the trend of is nonlinear monotonic increasing and tends to the horizontal asymptote defined by the obstacle nominal capacity . Once realizations reach , the probability of failure monotonically increases [Figs. 11(i–k)]. The characteristic times of failure of SMMs, embankment, and track correspond to the time at which the th percentile curve crosses . In Figs. 11(i–k), and are shown. For any considered , failure occurs in chronological order on the SMM, then on the embankment, and finally on the track. In particular, on the NW side of the track, sand maintenance shall be planned every 7–8 months on the SMM, every 11–13 months on the embankment, and every 2.2–2.4 years on the track. It is worth stressing that exceptional sandstorm events are not taken into account by the present modeling framework.
To assess the economic impact of the SMM implementation along the railway service life, the LCCA is carried out. The unmitigated design scenario is retained for reference to assess the extra costs and savings.
Fig. 12 plots the cumulated extra costs and savings per kilometer of alignment together with the normalized costs considering both sides of the railway line, being the cumulated sand maintenance costs resulting from the unmitigated scenario. Environmental uncertainties, embodied by , , and variability, are transposed into the LCCA by adopting resulting from different windblown sand action statistics, i.e., , , and . This in turn leads to the probabilistic description of and . The design and construction of the SMM involve initial additional costs. However, they are low if compared with sand maintenance costs and the break-even time (i.e., time at which ) takes place after about 3–5 years from completion. Depending on the chosen performance level, costs variation is significant. The savings with respect to the unmitigated scenario after 20 years amount to about , i.e., from 40% to 50% the cost of the unmitigated railway. Finally, they increment up to about after 100 years, i.e., 50%–60% of the cost of the unmitigated scenario. The average cost of civil works, track works, and construction contract for Ethiad railway stage two package D amounts to about (Railway Technology 2019). As a result, sand maintenance costs and savings are relevant even if compared with both the railway infrastructure worth and the SMM design and construction costs.
In a (SULS) performance level perspective, high percentiles of windblown sand action shall be considered. However, the gap between LCCA curves reflects the uncertainty in windblown sand action and depicts the windblown sand induced costs and SMM induced savings in a conservative perspective. Furthermore, it demonstrates that a careful description of randomness in windblown sand action is fundamental to infer windblown sand induced costs.
Conclusion
The present study introduces a novel hybrid approach to assess the life-cycle performance and cost analysis of SMMs for the design of infrastructures in sandy regions and the retrofitting of existing ones. Within an engineering design perspective, the technical applicability of this novel approach has been demonstrated by referring to a case study dealing with a desert railway protected by a design solution against windblown sand action.
Open field windblown sand transport is characterized through wind-sand tunnel tests on a flat sand bed. Multiphase wind-sand computational simulations are tuned on WST measurements and put into practice to assess the full-scale SMM life-cycle performance for the first time. Windblown sand action is then quantified by considering uncertainties related to wind and sand subfields in the specific building site. A time-variant reliability analysis is performed to assess characteristic times of failure and plan life-cycle sand removal maintenance on any infrastructure component, i.e., the SMM, the embankment, or the track. Finally, a life-cycle cost analysis is adopted as a design tool to quantify the economic impact of the SMM design.
From the obtained results, some general comments can be drawn in an SMM hybrid performance assessment and design perspective. First, the characterization of windblown sand transport and quantification windblown sand action is of paramount importance for civil structures and infrastructures in sandy environments, such as SMMs. In this context, the combination of WST tests and ETD simulation allow to: (1) overcome WST shortcomings resulting from the testing of scale SMMs (Raffaele et al. 2021), and (2) assess SMM life-cycle performance in virtual full-scale conditions. In turn, this would allow toially replace and/or assist expensive field tests, usually incompatible with infrastructure designer/stakeholders time requirements. Second, the robust design and assessment of SMMs involve the accurate estimation of windblown sand action and SMM life-cycle cost. This implies in turn the correct estimation of the SMM life-cycle performance and wind-sand field variability. This is in analogy with the well-known wind engineering analysis chain first introduced by Davenport (1961), in which each chain link is crucial for the correct assessment of action and effects on structures.
Given the wide research field, several experimental-computational research perspectives arise. In particular: (1) the ETD model can be improved by including momentum balance of the sand phase and turbulence corrections. This would allow to model the two-way interaction between sand and wind phases so to correctly simulate wind and sand velocity close to the ground. However, two-way coupling models are difficult to be set up due to the number of modeling parameters and high computational cost. As a result, they are not commonly used for engineering problems (Lo Giudice et al. 2019); (2) despite the overall good SMM life-cycle performance, the railway land width should always be minimized to minimize in turn the soil consumption and the related expected costs of land expropriation. The adoption of the investigated SMM implies a large railway land footprint, equal to about 130 m. As such, berm and ditch geometry variations but also other SMM solutions meeting such a requirement are worth to be investigated; (3) in a fully probabilistic approach perspective, it is worth considering the randomness of uncertain LCCA setup parameters, such as discount rate and unit costs, to investigate the uncertainty propagation to the life cycle cost; (4) from Bruno et al. (2018a) and Horvat et al. (2021), we can conjecture that clean SMMs performs better at high wind speed and under orthogonal winds. As such, further 3D ETD simulations are envisaged to evaluate the SMM life-cycle performance under yawed winds and for different wind speeds. The thorough description of LCP according to wind directionality and magnitude will allow in turn to refine the estimation of the life-cycle cost; (5) in an even more general perspective, this study paves the way towards a full hybrid approach to SMM performance assessment. In this framework, ETD simulations shall be always fully validated against scale WST tests. As a result, highly reliable WST tests on surface-mounted SMMs are required to characterize the wind flow, sand transport, and sand morphodynamics around them. In a fully hybrid performance assessment perspective, this will allow in turn to quantify the experimental distortion of life-cycle performance by recurring to full-scale constraints-free ETD simulations (Raffaele et al. 2021).
Data Availability Statement
Some or all data, models, or code generated or used during the study are proprietary or confidential in nature and may only be provided with restrictions. Available items: experimental wind-sand tunnel data. Restricted items: erosion-transport-deposition model code is restricted by a nondisclosure agreement.
Acknowledgments
The study has been developed in the framework of the MSCA-IF-2019 research project Hybrid Performance Assessment of Sand Mitigation Measures (HyPer SMM, hypersmm.vki.ac.be). This project has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant Agreement No. 885985. In particular, we acknowledge J. van Beeck, S. Buckingham from the von Karman Institute for Fluid Dynamics, and S. Khris, A. Lo Giudice, R. Nuca from Optiflow Company and Politecnico di Torino for the helpful discussions about the topics of the study. The study has been jointly developed in the framework of the Windblown Sand Modeling and Mitigation (WSMM, www.polito.it/wsmm) research, development, and consulting group established between Politecnico di Torino and Optiflow Company. Computational resources were provided by Optiflow Company.
References
Alghamdi, A. A., and N. S. Al-Kahtani. 2005. “Sand control measures and sand drift fences.” J. Perform. Constr. Facil. 19 (4): 295–299. https://doi.org/10.1061/(ASCE)0887-3828(2005)19:4(295).
Anno, Y. 1984. “Requirements for modeling of a snowdrift.” Cold Reg. Sci. Technol. 8 (3): 241–252. https://doi.org/10.1016/0165-232X(84)90055-7.
Basnet, K., G. Constantinescu, M. Muste, and H. Ho. 2014. “Method to assess efficiency and improve design of snow fences.” J. Eng. Mech. 141 (3): 04014136. https://doi.org/10.1061/(ASCE)EM.1943-7889.0000871.
Blocken, B., T. Stathopoulos, and J. Carmeliet. 2007. “CFD simulation of the atmospheric boundary layer: Wall function problems.” Atmos. Environ. 41 (2): 238–252. https://doi.org/10.1016/j.atmosenv.2006.08.019.
Bruno, L., and D. Fransos. 2015. “Sand transverse dune aerodynamic: 3D coherent flow structures from a computational study.” J. Wind Eng. Ind. Aerodyn. 147 (Dec): 291–301. https://doi.org/10.1016/j.jweia.2015.07.014.
Bruno, L., D. Fransos, and A. Lo Giudice. 2018a. “Solid barriers for windblown sand mitigation: Aerodynamic behavior and conceptual design guidelines.” J. Wind Eng. Ind. Aerodyn. 173 (Feb): 79–90. https://doi.org/10.1016/j.jweia.2017.12.005.
Bruno, L., M. Horvat, and L. Raffaele. 2018b. “Windblown sand along railway infrastructures: A review of challenges and mitigation measures.” J. Wind Eng. Ind. Aerodyn. 177 (Jun): 340–365. https://doi.org/10.1016/j.jweia.2018.04.021.
Cebeci, T., and P. Bradshaw. 1977. Momentum transfer in boundary layers. Washington, DC: Hemisphere.
CEN (European Committee for Standardization). 2005. Eurocode 1: Actions on structures—Part 1-4: General actions—Wind actions. EN 1991-1-4. Brussels, Belgium: CEN.
Centre National d’Etudes et de Recherches Integrees du Batiment. 2013. Document Technique Réglementaire C 2-4.7: Règlement Neige & Vent. DTR C 2-4.7. Algiers, Algeria: Centre National d’Etudes et de Recherches Integrees du Batiment.
Davenport, A. G. 1961. “The application of statistical concepts to the wind loading of structures.” Proc. Inst. Civ. Eng. 19: 449–471.
Dong, Z., H. Wang, X. Liu, and X. Wang. 2004. “The blown sand flux over a sandy surface: A wind tunnel investigation on the fetch effect.” Geomorphology 57 (1–2): 117–127. https://doi.org/10.1016/S0169-555X(03)00087-4.
Edgell, H. S. 2006. Arabian desert. Berlin: Springer.
Fabrycky, W., and B. Blanchard. 1991. Life cycle cost and economic analysis. Englewood Cliffs, NJ: Prentice Hall.
Ferreira, A. D., A. M. G. Lopes, D. X. Viegas, and A. C. M. Sousa. 1995. “Experimental and numerical simulation of flow around two-dimensional hills.” J. Wind Eng. Ind. Aerodyn. 54–55 (Feb): 173–181. https://doi.org/10.1016/0167-6105(94)00040-K.
Guo, Z., N. Huang, Z. Dong, R. Van Pelt, and T. Zobeck. 2014. “Wind erosion induced soil degradation in northern China: Status, measures and perspective.” Sustainability 6 (12): 8951–8966. https://doi.org/10.3390/su6128951.
Harris, A. D., W. J. Grogg, A. Akoma, B. J. Hayes, R. F. Reidy, E. F. Imhoff, and P. C. Collins. 2015. “Revisiting (some of) the lasting impacts of the liberty ships via a metallurgical analysis of rivets from the SS ‘John W. Brown’.” Jom 67 (12): 2965–2975. https://doi.org/10.1007/s11837-015-1668-1.
Ho, T. D., A. Valance, P. Dupont, and A. O. El Moctar. 2011. “Scaling laws in aeolian sand transport.” Phys. Rev. Lett. 106 (9): 094501. https://doi.org/10.1103/PhysRevLett.106.094501.
Horvat, M., L. Bruno, and S. Khris. 2021. “CWE study of wind flow around railways: Effects of embankment and track system on sand sedimentation.” J. Wind Eng. Ind. Aerodyn. 208 (Jan): 104476. https://doi.org/10.1016/j.jweia.2020.104476.
Horvat, M., L. Bruno, S. Khris, and L. Raffaele. 2020. “Aerodynamic shape optimization of barriers for windblown sand mitigation using CFD analysis.” J. Wind Eng. Ind. Aerodyn. 197 (Feb): 104058. https://doi.org/10.1016/j.jweia.2019.104058.
Hotta, S., and K. Horikawa. 1990. “Function of sand fence placed in front of embankment.” In Proc., 22nd Int. Conf. on Coastal Engineering, 2754–2767. Reston, VA: ASCE.
Ierimonti, L., L. Venanzi, and L. Caracoglia. 2018. “Life-cycle damage-based cost analysis of tall buildings equipped with tuned mass dampers.” J. Wind Eng. Ind. Aerodyn. 176 (May): 54–64. https://doi.org/10.1016/j.jweia.2018.03.009.
Iversen, J. D., and K. D. Rasmussen. 1994. “The effect of surface slope on saltation threshold.” Sedimentology 41 (4): 721–728. https://doi.org/10.1111/j.1365-3091.1994.tb01419.x.
Kok, J. F., E. J. R. Parteli, T. I. Michaels, and D. B. Karam. 2012. “The physics of wind-blown sand and dust.” Rep. Prog. Phys. 75 (10): 106901. https://doi.org/10.1088/0034-4885/75/10/106901.
Kubilay, A., D. Derome, B. Blocken, and J. Carmeliet. 2013. “CFD simulation and validation of wind-driven rain on a building facade with an Eulerian multiphase model.” Build. Environ. 61 (Mar): 69–81. https://doi.org/10.1016/j.buildenv.2012.12.005.
Launder, B. E., and D. B. Spalding. 1974. “The numerical computation of turbulent flows.” Comput. Methods Appl. Mech. Eng. 3 (2): 269–289. https://doi.org/10.1016/0045-7825(74)90029-2.
Le, V., and L. Caracoglia. 2019. “Life-cycle cost analysis of a point-like structure subjected to tornadic wind loads.” J. Struct. Eng. 146 (2): 04019194. https://doi.org/10.1061/(ASCE)ST.1943-541X.0002480.
Liu, B., J. Qu, W. Zhang, and G. Qian. 2011. “Numerical simulation of wind flow over transverse and pyramid dunes.” J. Wind Eng. Ind. Aerodyn. 99 (8): 879. https://doi.org/10.1016/j.jweia.2011.06.007.
Lo Giudice, A., R. Nuca, L. Preziosi, and N. Coste. 2019. “Wind-blown particulate transport: A review of computational fluid dynamics models.” Math. Eng. 1 (3): 508–547. https://doi.org/10.3934/mine.2019.3.508.
Lo Giudice, A., and L. Preziosi. 2020. “A fully Eulerian multiphase model of windblown sand coupled with morphodynamic evolution: Erosion, transport, deposition, and avalanching.” Appl. Math. Model. 79 (Mar): 68–84. https://doi.org/10.1016/j.apm.2019.07.060.
Menter, F. R., M. Kuntz, and R. Langtry. 2003. “Ten years of industrial experience with the SST turbulence model.” In Vol. 4 of Proc., 4th Int. Symp. on Turbulence, Heat and Mass Transfer: Turbulence Heat and Mass Transfer Series, edited by K. Hanjalić, Y. Nagano, and J. Tummers, 1208. Antalya, Turkey: Begell House.
Meroney, R. N. 2016. “Ten questions concerning hybrid computational/physical model simulation of wind flow in the built environment.” Build. Environ. 96 (Feb): 12–21. https://doi.org/10.1016/j.buildenv.2015.11.005.
O’Rourke, M., A. DeGaetano, and J. D. Tokarczyk. 2005. “Analytical simulation of snow drift loading.” J. Struct. Eng. 131 (4): 660–667. https://doi.org/10.1061/(ASCE)0733-9445(2005)131:4(660).
Phillips, D. A. 2011. “Analysis of potential sand dune impacts on railway tracks and methods of mitigation.” Accessed July 10, 2021. http://www.iktissadevents.com/files/events/gtrc/1/presentations/d2-s8-duncan-phillips.pdf.
Plaza, J. S., M. Lombardero Barceló, and P. Rodríguez de Lema Tapetado. 2012. “Sand and wind: An outline of the study of aeolian action on infrastructure with reference to Haramain High Speed Railway, Makkah–Al-Madinah.” Revista de Obras Públicas 159 (3537): 7–36.
Preziosi, L., D. Fransos, and L. Bruno. 2015. “A multiphase first order model for non-equilibrium sand erosion, transport and sedimentation.” Appl. Math. Lett. 45 (Jul): 69–75. https://doi.org/10.1016/j.aml.2015.01.011.
Raffaele, L., and L. Bruno. 2019. “Windblown sand action on civil structures: Definition and probabilistic modelling.” Eng. Struct. 178 (Jan): 88–101. https://doi.org/10.1016/j.engstruct.2018.10.017.
Raffaele, L., and L. Bruno. 2020. “Windblown sand mitigation along railway megaprojects: A comparative study.” Struct. Eng. Int. 30 (3): 355–364. https://doi.org/10.1080/10168664.2020.1714530.
Raffaele, L., L. Bruno, F. Pellerey, and L. Preziosi. 2016. “Windblown sand saltation: A statistical approach to fluid threshold shear velocity.” Aeolian Res. 23 (Dec): 79–91. https://doi.org/10.1016/j.aeolia.2016.10.002.
Raffaele, L., L. Bruno, and D. J. Sherman. 2020. “Statistical characterization of sedimentation velocity of natural particles.” Aeolian Res. 44 (Jun): 100593. https://doi.org/10.1016/j.aeolia.2020.100593.
Raffaele, L., J. van Beeck, and L. Bruno. 2021. “Wind-sand tunnel testing of surface-mounted obstacles: Similarity requirements and a case study on a sand mitigation measure.” J. Wind Eng. Ind. Aerodyn. 214 (Jul): 104653. https://doi.org/10.1016/j.jweia.2021.104653.
Railway Technology. 2019. “Etihad Rail awards $1.25B contract for Package D work.” Accessed July 10, 2021. https://www.railway-technology.com/news/etihad-rail-contract-package-d/.
Ribeiro, A. S., A. L. L. Silva, and A. M. P. de Jesus. 2011. “Evolution of fatigue history.” In Proc., Congresso Brasileiro de Engenharia Mecânica. Hampton, IA: ABCM.
Richards, P., and S. Norris. 2011. “Appropriate boundary conditions for computational wind engineering models revisited.” J. Wind Eng. Ind. Aerodyn. 99 (4): 257–266. https://doi.org/10.1016/j.jweia.2010.12.008.
Sañudo-Fontaneda, L. A., D. Castro-Fresno, J. J. Coz-Díaz, and J. Rodriquex-Hernandez. 2011. “Classification and comparison of snow fences for the protection of transport infrastructures.” J. Cold Reg. Eng. 25 (4): 162–181. https://doi.org/10.1061/(ASCE)CR.1943-5495.0000031.
Tominaga, Y. 2018. “Computational fluid dynamics simulation of snowdrift around buildings: Past achievements and future perspectives.” Cold Reg. Sci. Technol. 150 (Jun): 2–14. https://doi.org/10.1016/j.coldregions.2017.05.004.
Tominaga, Y., T. Okaze, and A. Mochida. 2018. “Wind tunnel experiment and CFD analysis of sand erosion/deposition due to wind around an obstacle.” J. Wind Eng. Ind. Aerodyn. 182 (Nov): 262–271. https://doi.org/10.1016/j.jweia.2018.09.008.
Valance, A., K. R. Rasmussen, A. O. El Moctar, and P. Dupont. 2015. “The physics of aeolian sand transport.” C.R. Phys. 16 (1): 105–117. https://doi.org/10.1016/j.crhy.2015.01.006.
White, B. R. 1996. “Laboratory simulation of aeolian sand transport and physical modeling of flow around dunes.” Ann. Arid Zone 35 (3): 187–213.
Xin, G., N. Huang, J. Zhang, and H. Dun. 2021. “Investigations into the design of sand control fence for Gobi buildings.” Aeolian Res. 49 (Feb): 100662. https://doi.org/10.1016/j.aeolia.2020.100662.
Information & Authors
Information
Published In
Copyright
This work is made available under the terms of the Creative Commons Attribution 4.0 International license, https://creativecommons.org/licenses/by/4.0/.
History
Received: Jul 30, 2021
Accepted: Jan 25, 2022
Published online: May 3, 2022
Published in print: Jul 1, 2022
Discussion open until: Oct 3, 2022
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
- Luca Caracoglia, Jamie E. Padgett, Ilaria Venanzi, Risk-Informed and Life-Cycle Analyses of Structures and Infrastructures, Journal of Structural Engineering, 10.1061/(ASCE)ST.1943-541X.0003495, 148, 12, (2022).