Rotational test spaces for a fully-implicit FVM and FEM for the DNS of fluid-particle interaction
KAUST DepartmentComputer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division
Applied Mathematics and Computational Science Program
Extreme Computing Research Center
Online Publication Date2019-05-23
Print Publication Date2019-09
Permanent link to this recordhttp://hdl.handle.net/10754/656472
MetadataShow full item record
AbstractThe paper presents a fully-implicit and stable finite element and finite volume scheme for the simulation of freely moving particles in a fluid. The developed method is based on the Petrov-Galerkin formulation of a vertex-centered finite volume method (PG-FVM) on unstructured grids. Appropriate extension of the ansatz and test spaces lead to a formulation comparable to a fictitious domain formulation. The purpose of this work is to introduce a new concept of numerical modeling reducing the mathematical overhead which many other methods require. It exploits the identification of the PG-FVM with a corresponding finite element bilinear form. The surface integrals of the finite volume scheme enable a natural incorporation of the interface forces purely based on the original bilinear operator for the fluid. As a result, there is no need to expand the system of equations to a saddle-point problem. Like for fictitious domain methods the extended scheme treats the particles as rigid parts of the fluid. The distinguishing feature compared to most existing fictitious domain methods is that there is no need for an additional Lagrange multiplier or other artificial external forces for the fluid-solid coupling. Consequently, only one single solve for the derived linear system for the fluid together with the particles is necessary and the proposed method does not require any fractional time stepping scheme to balance the interaction forces between fluid and particles. For the linear Stokes problem we will prove the stability of both schemes. Moreover, for the stationary case the conservation of mass and momentum is not violated by the extended scheme, i.e. conservativity is accomplished within the range of the underlying, unconstrained discretization scheme. The scheme is applicable for problems in two and three dimensions.
CitationHöllbacher, S., & Wittum, G. (2019). Rotational test spaces for a fully-implicit FVM and FEM for the DNS of fluid-particle interaction. Journal of Computational Physics, 393, 186–213. doi:10.1016/j.jcp.2019.05.004
SponsorsThe authors are grateful to D. Logaschenko for the useful discussions and thank to S. Reiter and A. Vogel for the support with the implementation in UG4.
JournalJournal of Computational Physics
Except where otherwise noted, this item's license is described as This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).