F-DEMPack Tutorial 1: Curved pipe
Before starting with this tutorial, the user is strongly encouraged to follow the G-DEMPack Tutorial 1: Conveyor belt 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.[Fluid-DEMPack_Tutorial http://kratos-wiki.cimne.upc.edu/images/0/0b/D_DEMPack2_Tutorial_4.gid.zip]
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:
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.
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.
The FEM2DEM mesher was used to obtain this initial mesh of DEM spheres. See G-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
The parameters in this section are identical to those in the G-DEMPack problem type, so no extra information is necessary in this case.
Inlet DEM group
See G-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
This section is also identical to its counterpart in G-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:
The DEM part is identical to G-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:
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:
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:
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:
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:
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.
The fluid solver needs some initial conditions in the fluid to solve the problem. The next figure shows the section:
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.
- Important remarks
- Make sure that the whole fluid boundary has some fluid condition applied. Each surface of the fluid external skin must be either Inlet, Outlet, Slip or No-slip.
- Make sure that entities not belonging to the external skin are NOT selected when applying the fluid 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:
A vertical constant flow of 1m/s entering the bottom opening of the tube was chosen.
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:
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.
- Outlet pressure
The upper surface has been marked as an outlet. The value of the pressure can be chosen. Typically 0.0.
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:
On the other hand, the initial DEM spheres mesh should look like the following:
While the inlet surface mesh should be very similar to this:
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:
The next figure shows the stationary pressure field in the fluid:
The user can also make a cut in the mass of fluid and obtained a 2D view of the velocity field in the fluid: