Numerical Simulation of Barchan-barchan repulsion

Published: 13 March 2024| Version 1 | DOI: 10.17632/ypkgwjfr4r.1
Nicolao Lima, Willian Righi Assis, Erick Franklin


The dataset shared within these files consists of the full setup to perform a numerical simulation of a subaquous barchan-barchan dune interaction. The simulation was carried out using the open-source free CFDEM software, which couples the open-source CFD code OpenFOAM with the open-source DEM code LIGGGHTS. The simulations were previously tested and validated (considering only one single barchan dune) against experimental results and results from the literature. The dataset contains a complete example case. We also share the full Octave code for post-processing those results. To understand and run the case, the user should be noticed that previously experience with both OpenFOAM and LIGGGHTS are recommended.


Steps to reproduce

The full setup of the case is stored inside the Base_Case folder. It works basically like this: the turbulence is initially developed by running only the numerical simulation of the fluid with OpenFOAM; then, the initial piles of grains are set within a stationary fluid; finally, both of these results are put together to develop the barchans and watch their dynamics. 1) Before starting the numerical simulations, the user should run a LES simulation using the same mesh and turbulence model provided in this setup with only OpenFOAM (without particles) which will be used as startup conditions for the numerical simulations on the next steps. It is recommended to use the perturbU utility of OpenFOAM to reach faster turbulence development. The dictionary used for this purpose, "perturbUDict", is provided inside the CFD/constant folder; 2) Then, the initial piles will be set by running the "allrunMesh", "allrunDEM", "allrunDEM_restart" and "allrunCFDEM" executable files in that order. They will automatically place and settle the grains inside the domain; 3) Once the grains have settled, the user needs to copy all the last OpenFOAM results (initial conditions) from the turbulence development calculation to the 0 "zero" folder of this case; 4) Finally, the user can run the "allrun_restart" executable to run the final simulation. To analyze the results, an Ocatve code is provided and located inside the Octave_Code folder. There, the user will find two different folders: the "Double_Dune_Octave_Traj" contains the codes to analyze the results from the LIGGGHTS software, and will extract informations such as morphodynamics, trajectories forces and so on. This code needs to be executed before the codes located in the "Double_Dune_Octave_Flow_L" folder. The latter, is responsible for analyzing the fluid flow results, and will provide the results for the turbulence flow over the dunes. All those codes are fully commented and the specific steps to execute them are also highlighted inside the "inputData" files.


Universidade Estadual de Campinas


Dune, Computational Fluid Dynamics, Discrete Element Method