The overview of the F-DEMPack2 program is the following:
General Application Data
These are the contents of the General Application Data menu:
It has the following sections:
The user can fill here several fields with information related to the current simulation.
Here the user can choose the Debug option in order to see the complete menu despite a possible loss of robustness.
It consists of a few settings concerning the whole simulation, in particular the simulation duration, the output delta time and the number of threads to use in the simulation. At the moment the code runs only in OpenMP, MPI is currently under development.
It deals with general physical settings. At this moment it consists only of the gravity value.
It refers to those parameters concerning the inter-component coupling.
It has 4 possible values: No interactions, One way, Soft two-way and full two-way, which correspond, respectively, to no-interactions, only fluid to DEM, to unmodified fluid equations or to averaged Navier-Stokes.
Time averaging mode
- Take newest values: Take the present value of the hydrodynamic forces from DEM for the fluid integration.
- Take uniform average over DEM substeps: Average values of the hydrodynamic force over all substeps within one fluid step.
Interaction start time
The user can set here the time to start the interaction between DEM and fluid.
Hydrodynamic Force Model
- Pick each individual force: Pick specific expressions for the usual terms, one at a time, even though compatibility may be not guaranteed.
- Buoyancy Force type: Here the desired buoyancy law for the DEM particles can be chosen.
- Maxey-Riley Equations: To use these equations for the simulation, even though they are addressed to small Reynolds numbers.
- Drag Modifier Type: The user can choose between different drag laws available.
- Include Faxen terms: Including them may improve the accuracy of the results. Nevertheless, they require the calculation of the local laplacian and its derivative, which involves a certain overhead.
Those are options concerning variables that are calculated by post-processing the fluid-field solution. Here the user can choose the Gradient and the Velocity Laplacian Calculation Modes. Both are computed as a post-process from the fluid solution.
This section contains the options concerning the results to be post-processed. It allows the printing of many nodal values for both DEM elements and fluid mesh nodes. Discretion is advised when printing results in postprocess, because much memory may be needed to be allocated, harming performance in many cases.
These are the contents of the Fluid menu:
It has the following sections:
The user can choose between Monolithic and Fractional Step. The monolithic solver is slower, but remarkably more accurate than the fraction step solver.
The user can decide to use no turbulence model or to choose between Smagorinsky-Lilly or Spalart-Allmaras. In general, Smagorinsky-Lilly gets better results for separated boundary layers but needs finer meshers.
It has many options to adjust the behaviour of the solver. The user should not bother to change those default values.
The user can enter here the delta time to use in the computations of the fluid part as well as the step for the divergence clearance.
In this section, as many Fluid Properties as desired can be created. The standard Propertry, for example, is that related to the Water Material. The user, when later creating the different Fluid Elements, will have to assign every Element to its corresponding Property.
You must assign here the properties of the existing fluid groups. First of all, the Element Type that will be used in the computations, in this case the Tetrahedra element being the only available for 3D, and the corresponding Property for the group. Check the D-DEMPack Tutorial 4: Swimming DEM to know how this section works.
The user can assign here an initial velocity field for a given mass of fluid.
Here an initial pressure field can be assigned to the desired fluid group.
Here a fluid inlet entity can be assigned to a given group. The user must assign the group an inlet velocity field for the group as well as its fixity depending on the three axis.
The user can assign here an outlet pressure for a given group representing an outlet of fluid.
When assigning this condition to a solid boundary group, the fluid will have in general a non-zero velocity relative to the boundary.
On the contrary, when assigning this condition to a boundary, the fluid will have a null velocity in relation to the boundary.
Please, refer to D-DEMPack2 manual for a full description of the menus in this section.