how to define meshwalls in both in.liggghts_init and in.liggghts_run

Submitted by yirenxindudu on Fri, 01/18/2019 - 04:04

Dear all,

recently, I faced a problem with CFD-DEM coupling.
both in.liggghts_init and in.liggghts_run need to define "mesh/surface" . I defined wall using STL files located in folder "meshes", the commends are as follows:
fix cv all mesh/surface file meshes/wall.stl type 1 #surface_vel -4.5 0. 0.
fix bx all mesh/surface file meshes/inlet.stl type 1
fix inface all mesh/surface file meshes/outlet.stl type 1
fix wall all wall/gran model hertz tangential history mesh n_meshes 2 meshes bx cv

However, it went wrong as follows:
fix cv all mesh/surface file meshes/wall.stl type 1 #surface_vel -4.5 0. 0.
ERROR on proc 0: Cannot open mesh file meshes/wall.stl (../input_mesh_tri.cpp:105)
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 DUP FROM 0
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------

I think it is also possible to define walls using "wall/gran" method, and the above problem can be solved, but I don't know how to define cone wall.

Thank you!

AttachmentSize
Plain text icon log file of my case4.74 KB

Daniel Queteschiner | Fri, 01/18/2019 - 09:25

A pure LIGGGHTS simulation (init) is usually executed in the DEM folder. CFDEMcoupling simulations on the other hand are controlled by the CFD side, i.e. the execution typically takes place in the CFD folder of the coupled case, not in the DEM folder. Hence the path to the mesh files must be adapted accordingly. A solution that works for executing the LIGGGHTS part of the simulation from the CFD or DEM folder may look like this:
fix cv all mesh/surface file ../DEM/meshes/wall.stl type 1 #surface_vel -4.5 0. 0.
fix bx all mesh/surface file ../DEM/meshes/inlet.stl type 1
fix inface all mesh/surface file ../DEM/meshes/outlet.stl type 1

yirenxindudu | Mon, 01/21/2019 - 15:41

Thank you very much! It works!
In addition, I found that the quality of mesh (STL) file can effect the calculation speed, it should be further meshed using netgen etc after being exported from solidworks.
Best regards