Distribute drag forces cfd-dem simulation

[yashuuzi]

Description

In current code of cfd-dem solver, Each quadrature point have the same value, and only consider the particle inside the cell. To solve this problem, I propose to evaluate the drag force at each quadrature point using reference sphere.

1. When calculate the particle-fluid-interaction, stored the beta_drag as each particle's propertied.
2. When assemble the drag force, evaluate drag force at each quadrature point using reference sphere. To do so, we need to consider the neighbor cell like QCM method.

Documentation

I add the parameter "distribute drag force" to enable and disable this future easily.

Miscellaneous (will be removed when merged)

I modified only Rong drag model

Checklist (will be removed when merged)

See this page for more information about the pull request process. Code related list:

• [ ] All in-code documentation related to this PR is up to date (Doxygen format)
• [ ] Lethe documentation is up to date
• [ ] New feature has unit test(s) (preferred) or application test(s)
• [ ] The branch is rebased onto master
• [ ] Changelog (CHANGELOG.md) is up to date
• [ ] Code is indented with indent-all and .prm files (examples and tests) with prm-indent

Pull request related list:

• [ ] Labels are applied
• [ ] There are at least 2 reviewers (or 1 if small feature) excluding the responsible for the merge
• [ ] If this PR closes an issue or is related to a project, it is linked in the "Projects" or "Development" section
• [ ] If the fix is temporary, an issue is opened
• [ ] The PR description is cleaned and ready for merge

[yashuuzi]

When it comes to calculating drag force in "calculate_particle_fluid_interaction", I realized I don't need new set of function. But I think I need new function for assemble_matrix and assemble_rhs because these function need to have an access to the cell, particlehandler, and data members in GLSVANSSolver. That's why I included gls_vans.h in vans_assembler.cc.

I think my explanation is not enough, so you can take a look at my document about implementation