Numerical Simulation of Barchans interacting with dune-size obstacles
Description
The dataset consists of the setup to reproduce a simulation of a subaqueous barchan interacting with an obstacle. 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 dataset contains a complete example of a by-pass case. Additionally, we share all Octave codes used for post-processing the results.
Files
Steps to reproduce
The full setup of the case is stored inside the Base_Case folder. The turbulence is initially developed by running only the numerical simulation of the fluid with OpenFOAM; then, the initial pile of grains are set within a stationary fluid in the domain; finally, both of these results are put together to develop the barchan and watch its dynamics with the obstacle. 1) 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. Note that the provided example is for a by-pass case (if the user wants to reproduce another case, there is the need to change the mesh for the desirable obstacle). The user should run the "allrun" executable file in the folder "Fluid_Initialization". 2) The initial pile will be set by running the "allrun" executable file in the folder "particleInitialization". 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 in the folder "run" of this case. Additionally, the user should copy the file restart from "particleInitialization" into the DEM/post/restart path of the "run" folder; 4) Finally, the user can run the "allrunMesh" and "allrun_CFDEM" executable files in the folder "run" to run the final simulation. 5) To analyze the results, the user can find within the folder "run" a folder named "Dune_Block_Octave_Traj" which contains the codes to analyze the results from the LIGGGHTS software, and will extract information such as trajectories, forces, velocity vectors and so on. 5.1) The file "inputData.m" needs to be adjusted with the initial conditions of the test; 5.2) The user should run mainCode.m file. 5.3) In the end, there is an option to adjust the plotData.m file with the desirable graphs format;