F-DEMPack Tutorial 1: Curved pipe

From KratosWiki
Revision as of 13:43, 23 January 2017 by Mlnunez (Talk | contribs)
Jump to: navigation, search



Before starting with this tutorial, the user is strongly encouraged to follow the D-DEMPack Tutorial 2: DemPack 2.0 to get a feeling of how the problem type works, and in particular the DEM section. This tutorial will focus mainly in the Fluid section and its particularities. The user must start by downloading the File:D DEMPack2 Tutorial 4.gid.zip. This file has already created the groups that will be used in the simulation. It has also already assigned the DEM groups, as well as the mesh sizes for all groups.


The geometry of study consists of a curved tube through which a flux of water passes. Additionally, an initial mass of DEM elements exists inside the mass of water, as well as an inlet creating DEM particles with time. This geometry can be observed in the picture that follows:

DEM 2 tutorial swimming geometry.png

The geometry has five groups, as seen in the next picture. They consist of a sphere made up of DEM elements, placed at the middle of the tube, a cuadrilaterial inlet of DEM spheres near the bottom of the tube, the mass of fluid inside the tube, the inlet of fluid in the opening at the bottom of the tube and the walls of the tube. The next picture shows the assignation of group by different entities that the user is expected to find in the file available for downloading.

DEM 2 tutorial swimming groups.png

DEM Entities

The file contains several entities and conditions in relation to the DEM part of the problem. As previouly said, they are already preassigned so the user does not have to bother and can concentrate on the fluid aspects and details of the simulation. Nevertheless, figures showing the details on the DEM parts will be added here for the sake of completion and as a reference should the user lose these settings or in the case of have any problem when loading the file.

DEM group

DEM 2 tutorial swimming DEM group.png

The FEM2DEM mesher was used to obtain this initial mesh of DEM spheres. See D-DEMPack Tutorial 3: DEM Meshers for a reference on the different meshers in the program. In this case, the chosen distribution of diameters for the elements is of 2mm with no variance.

DEM-FEM wall group

DEM 2 tutorial swimming DEM FEM group.png

The parameters in this section are identical to those in the D-DEMPack problem type, so no extra information is necessary in this case.

Inlet DEM group

DEM 2 tutorial swimming DEM inlet group.png

See D-DEMPack Tutorial 2: DemPack 2.0 for details. It is important to note that it is still not possible to create neither clusters nor nanoparticles from an inlet entity. This aspect of the code is still under development.

DEM Initial Conditions

DEM 2 tutorial swimming DEM Initial Conditions.png

This section is also identical to its counterpart in D-DEMPack, so no further explanations are needed.


The Materials section in F-DEMPack contains the material data for both DEM and fluid elements. See the next figure for an overview of that section:

DEM 2 tutorial swimming materials section.png

The DEM part is identical to D-DEMPack, and a full explanation of every aspect in it can be found in the corresponding links given above. In this case, though, an additional section for the Fluid part exists, where the user can set the values of some fundamental fluid properties, as for example density, viscosity, bulk modulus or rheological characteristics.

DEM-Fluid Interaction Settings

Most of the interaction parameters between the DEM spheres and fluid are inside the General Application Data section, whose overview is given next:

DEM 2 tutorial swimming general application data.png

Some of the parameters in this tree are very straightforward, as for example the duration time, the output delta time, the number of threads to use in the simulation or the gravity vector. Others, though, carry a higher difficulty and are mostly related to the way the two subdomains interact. A deeper explanation of those parameters can be found in F-DEMPack2 manual. The previous screenshot has been given to the user as a reference as well as a guide for choosing some default values that give good results in this particular case. This tree also includes the Results section, which is pretty straightforward.


This section contains the information about the properties of the fluids, the different existing fluid elements and, when necessary, their assigned conditions. It also has some parameters in relation to the settings of the fluid solver. The figure that follow shows an overview of this section:

DEM 2 tutorial swimming fluid menu.png

The screenshot shows the chosen parameters in this example. The user is encouraged to play a little bit to those values and see the results. In this case a monolithic solver was chosen to get more accurate results but no turbulence model was considered necessary. The linear solver parameters require a much deeper understanding of the underlying theory and the associated numerical methods and will not be discussed here, check for more information. The default parameters in the problem type gave good results in this particular case. Finally, the user can enter the desired computational time step, which does not have to be too small as long as stable simulations are obtained.


In this section, the user can create a Property related with each of the fluids in the problem. In order to do this, the corresponding fluid must have been previously created in the Materials section. The next picture shows the process:

DEM 2 tutorial swimming Fluid Properties.png

In this particular case, the Water default material was assigned to the Property1, which is good enough in this sample simulation.


The user must specify the fluids that will be present in the problem, so the next step is to assign the desired groups to the Fluid Elements. The process is the same as in other tutorials. A screenshot is given next:

DEM 2 tutorial swimming Elements Fluid.png

To finish the elements assignation, though, an additional step is necessary. The user must specify the FEM element type to be used in the fluid mesh and the corresponding Property. The next figure shows this:

DEM 2 tutorial swimming Elements Fluid Property.png

In this case, the only available element in 3D is the tetrahedra, while the chosen property was number 1, corresponding to the Water material.


We finish the process by assigning the necessary conditions to their corresponding groups.

Initial Conditions

The fluid solver needs some initial conditions in the fluid to solve the problem. The next figure shows the section:

DEM 2 tutorial swimming Initial Conditions.png

For this simulation, an initial vertical value of 1m/s was assigned to the mass of fluid inside the tube. In this case, no initial pressure was necessary.

Boundary Conditions

  • Inlet velocity

An inlet velocity is necessary in this simulation for a flow to exist in time along the interior of the tube. See the next picture:

DEM 2 tutorial swimming inlet group.png

A vertical constant flow of 1m/s entering the bottom opening of the tube was chosen. No Outlet Pressure was needed in the simulation.

  • No-Slip

The user must tell the program if there exists any relative velocity between the fluid and the solid boundary. The capture that follows shows the section:

DEM 2 tutorial swimming No Slip.png

For this example, a no-slip behaviour was given to the walls. If there had been additional DEM-FEM entities, an Is-Slip behaviour could have been given to them.

Meshing and Running

If the user has succesfully assigned to the corresponding groups all the previous Properties, Elements and Conditions, the only thing that is still missing is the computational mesh. As already said, all the messing characteristics has been previously set in this problem, so the user has been able to focus in the Fluid aspects of the problem type. So to mesh, hit Ctr-G and enter a value of 0.02 and press OK. A global view of the resulting mesh should be similar to this one:

DEM 2 tutorial swimming mesh.png

On the other hand, the initial DEM spheres mesh should look like the following:

DEM 2 tutorial swimming mesh DEM detail.png

While the inlet surface mesh should be very similar to this:

DEM 2 tutorial swimming mesh DEM inlet detail.png


Once a mesh is obtained and the file saved, the simulation can be run. To start the computations, the user must only hit F5 and the simulation will begin. The next four captures show the resulting simulation at different times:

DEM 2 tutorial swimming velocities1of4.png DEM 2 tutorial swimming velocities2of4.png DEM 2 tutorial swimming velocities3of4.png DEM 2 tutorial swimming velocities4of4.png

The next figure shows the stationary pressure field in the fluid:

DEM 2 tutorial swimming fluid pressure.png

The user can also make a cut in the mass of fluid and obtained a 2D view of the velocity field in the fluid:

DEM 2 tutorial swimming fluid velocities.png

Personal tools