|
|
Improving computational efficiency of contact solution in fully resolved CFD-DEM simulations with arbitrarily-shaped solids
Studeník, Ondřej ; Kotouč Šourek, M. ; Isoz, Martin
The abundance of industrial processes containing both solid and liquid phases generate demand for fully resolved models allowing for detailed analysis and optimization of these processes. An established approach providing such models is based using a variant of an immersed boundary method to couple the computational fluid dynamics (CFD) and discrete element method (DEM). In the talk, we will present our custom and monolithic implementation of a fully-resolved CFDDEM solver and concentrate on the intricacies of solving contact between two arbitrarily-shaped solids. We shall propose an efficient contact treatment based on the concept of a virtual mesh, which provides the mesh resolution required by DEM through dividing the space around the contact point in a finite volume fashion without any changes to the CFD mesh itself. A substantial part of the talk will devoted to the parallelization of the contact solution, especially in the context of the domain decomposition method imposed by the CFD solver.
|
|
|
Simulating particle-laden flows: from immersed boundaries towards model order reduction
Isoz, Martin ; Kubíčková, Lucie ; Kotouč Šourek, M. ; Studeník, Ondřej ; Kovárnová, A.
Particle-laden flow is prevalent both in nature and in industry. Its appearance ranges from the trans-port of riverbed sediments towards the magma flow, from the deposition of catalytic material inside particulate matter filters in automotive exhaust gas aftertreatment towards the slurry transport in dredging operations. In this contribution, we focus on the particle-resolved direct numerical simulation (PR-DNS) of the particle-laden flow. Such a simulation combines the standard Eulerian approach to computational fluid dynamics (CFD) with inclusion of particles via a variant of the immersed boundary method (IBM) and tracking of the particles movement using a discrete element method (DEM). Provided the used DEM allows for collisions of arbitrarily shaped particles, PR-DNS is based (almost) entirely on first principles, and as such it is a truly high-fidelity model. The downside of PR-DNS is its immense computational cost. In this work, we focus on three possibilities of alleviating the computational cost of PR-DNS: (i) replacing PR-DNS by PR-LES or PR-RANS, while the latter requires combining IBM with wall functions, (ii) improving efficiency of DEM contact solution via adaptively refined virtual mesh, and (iii) developing a method of model order reduction specifically tailored to PR-DNS of particle-laden flows.
|
|
|
Octree-generated virtual mesh for improved contact resolution in CFD-Dem coupling
Studeník, O. ; Kotouč Šourek, M. ; Isoz, Martin
The present work is focused on improving the efficiency of a computational fluid dynamics (CFD) – discrete element method (DEM) solver allowing for computations with non-spherical solids. In general, the combination of CFD and DEM allows for simulations of freely moving solid particles within a computational domain containing fluid. The standard approach of CFD-DEM solvers is to approximate solid bodies by spheres, the geometry of which can be fully defined via its radius and center position. Consequently, the standard DEM contact models are based on an overlap depth between particles, which can be easily evaluated for a sphere-sphere contact. However, for a contact between two non-spherical particles, the overlap depth cannot be used and has to be replaced by the more general overlap volume. The precision of the overlap volume computation is (i) crucial for the correct evaluation of contact forces, and (ii) directly dependent on the computational mesh resolution. Still, the contact volume evaluation in DEM for arbitrarily shaped bodies is usually by at least one order of magnitude more demanding on the mesh resolution than the CFD. In order to improve the computational efficiency of our CFD-DEM solver, we introduce the concept of an OCTREEbased virtual mesh, in which the DEM spatial discretization is adaptively refined while the CFD mesh remains unchanged.
|
host ::
RSS.