Advancements for the Numerical Simulation of Free Fall Penetrometers and the Analysis of Wind Erosion of Sands

The coastal population is growing, putting extra stress on coastal sediments and protection features, such as beach dunes. Moreover, global warming will increase the frequency of storms, and coastal dunes and other defense infrastructure will be subjected to increased erosion and scouring, endangering the people they are meant to protect. Understanding soil dynamics and fluid interaction is crucial to predict the effects of sand erosion. In particular, the study of wind erosion of sands in coastal dunes is essential due to the protective role these earthen structures have during storm events.

One of the challenges about predicting wind erosion in coastal dunes is its extended spatial scale and the associated economic and logistics costs of sampling and characterizing the sediments. Because of this, in-situ testing for sediment characterization is essential. In particular, the usage of free-fall penetrometers (FFP) is appealing due to their portability and robustness. The sediment properties obtained with this type of testing can later be used to assess wind erosion susceptibility by determining, for example, the wind velocity to initiate the erosion process.

FFP testing involves dropping an instrumented probe that impacts the soil and measures the kinematics or kinetics during the penetration process. For example, deceleration measurements are used to compute an equivalent quasi-static failure, which is not in line with the dynamic process characteristic of FFP testing. This preassumed failure mechanism is used to back-calculate the sand's geomechanical properties. However, soil behavior is highly complex under rapid loading, and incorporating this behavior into FFP sediment characterization models is challenging. Advanced numerical modeling can improve the understanding of the physics behind FFP testing.

This thesis presents various advancements in numerical modeling and erosion models to bridge FFP in-situ testing with predicting the initiation of wind erosion of sands. First, improvements oriented to the Material Point Method (MPM) for modeling in-situ FFP testing are proposed. The numerical results show that the simulation of FFP deployment in sands is affected by strain localization and highlight the importance of considering constitutive models sensitive to different loading rates. Because of the importance of rate effects in soil behavior, the second aspect of this thesis proposes a novel consistency framework. Two constitutive models are adapted to study strain-rate sensitive non-cohesive materials: i) a strain-softening Mohr-Coulomb, and ii) a NorSand model. In addition to increased strength, the proposed framework captures increased dilatation, an early peak deviatoric stress, and relaxation.

Finally, a novel sand erosion model is derived using a continuum approximation and limit equilibrium analysis. The erosion law considers geotechnical parameters, the effects of slope, and moisture suction, in a combined manner. The proposed model is theoretically consistent with existing expressions in the literature. It covers a wide range of environmental and geometrical conditions and helps to reconcile the results from FFP testing with the prediction of the initiation of wind erosion. The model was validated in a wind tunnel and is demonstrated to be a viable alternative for predicting sand erosion initiation.

This thesis opens up new research prospects, such as improving the soil characterization models or the direct prediction of sand erosion using rapid, reliable, and efficient in-situ testing methods.