Windblown sand transport, accumulation and erosion around structures and infrastructures in arid regions poses significant threats about their safety and serviceability. The mathematical modeling of the above phenomena is a challenging task, because of the interaction between the multiple phases of wind and sand, the time-dependent free interface among them, the multiple scales from km-long infrastructures to mm-sized sand grains through the turbulent wind energy cascade. The accurate and efficient computational simulation of the above phenomena is of paramount scientific and industrial importance, because of the lack of similarity in scale wind tunnel tests and of the costs and uncertainties in full scale field tests.