|
|
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.
|
|
|
Estimating rheological properties of suspensions formed of arbitrarily-shaped particles via CDF-Dem
Kotouč Šourek, M. ; Isoz, Martin
In recent years, new methods combining computational fluid dynamics (CFD) and discrete element method (DEM) have been intensively studied. Usually, these methods are focused on simulations of spherical particles. Nevertheless, this is inadequate for a simulation of a common suspension, the rheology of which is affected by particle shapes. In this work, we leverage the capabilities of an in-house developed CFD-DEM solver to simulate suspensions formed of arbitrarily-shaped particles. Specifically, we simulate a rheological measurement to estimate the suspension viscosity. The CFD-DEM estimates are in very good agreement with available experimental data and correlations proving the new solver capabilities regarding firstprinciples-based simulations of complex non-Newtonian suspension behaviour. The practical potential of suspension simulation is illustrated in a numerical study of the washcoating process in the preparation of a catalytic filter for automotive exhaust gas after-treatment.
|
|
|
Recent improvements in CFD solver for fully coupled particle-laden flows
Šourek, M. ; Isoz, Martin
While new methods combining the computational fluid dynamics (CFD) and the discrete element method (DEM) have been developed to simulate freely moving solid particles, they tend to be focused on simulations with spherical particles only. Here, we present a strongly coupled CFD-DEM solver capable of simulating movement of arbitrarily-shaped particles dispersed in a fluid. The particles are assumed to be large enough to affect the fluid flow and distributed densely enough to come into contact with both the boundaries of the computational domain and with each other. In this paper, we will focus on the recent improvements of our solver, specifically, in the areas of (i) inclusion of solid bodies into the computational domain, (ii) general CFD-DEM coupling algorithm, and (iii) code parallelization and practical usability.
|
|
|
Development of CFD solver for four-way coupled particle-laden flows
Šourek, M. ; Isoz, Martin
Computational uid dynamics (CFD) simulations containing freely moving bodies are still a challenging topic. More so, if the bodies are large enough to a_ect the uid ow and distributed\ndensely enough to come in contact both with the boundaries of the computational domain and with each other. In this work, we concentrate on the topic of simulation of (i) irregular bodies\nwith ow-induced movement and contact with computational domain boundaries taken into account, and (ii) bodies entrained by the uid and coming in contact not only with the domain\nboundaries but also with each other. The developed modeling approach is based on the hybrid _ctitious domain-immersed boundary method extended by the discrete element method. The\npresent contribution is focused on presentation of simulation principles and results of initial benchmark cases.
|
|
|
POD-DEIM-based model order reduction for four-way coupled fluid-solid flows
Isoz, Martin ; Šourek, M.
Proper orthogonal decomposition (POD) and discrete empirical interpolation method (DEIM) have become established tools for model order reduction in simulations of fluid flows. However, including moving solid bodies in the computational domain poses additional issues with respect to the fluid-solid coupling and to the solution of the movement of the solids. Still, it seems that if the hybrid ctitious domain-immersed boundary method is used to include the solids in the flow domain, POD-DEIM based approaches may be extended for four-way coupled particleladen flows. The present work focuses on the construction of POD-DEIM based reduced order models for the aforementioned flows.
|
|
|
Geometrically realistic macro-scale model for multi-scalesimulations of catalytic filters for automotive gasaftertreatment
Hlavatý, Tomáš ; Isoz, Martin ; Plachá, M. ; Šourek, M. ; Kočí, P.
This paper is part of a research focused on simulating (i) the catalytic conversion of environment endangering gases, and (ii) trapping of the particulate matter in automotive exhaust gas aftertreatment. Historically, the catalytic conversion and the filtration of soot particles were performed in independent devices. However, recent trend is to combine the catalytic converter and soot filter into a single device, the catalytic filter. Compared to the standard two-device system, the catalytic filter is more compact and has lower heat losses. Nevertheless, it is highly sensitive to the catalyst distribution. This study extends our recently developed methodology for pore-scale simulations of flow, diffusion and reaction in the coated catalytic filters. The extension consists of enabling data transfer from macro- to pore-scale models by preparing geometrically realistic macro-scale CFD simulations. The simulation geometry is based on XRT scans of real-life catalytic filters. The flow data from the newly developed macro-scale model are mapped as boundary conditions into the pore-scale simulations and used to improve the estimates of the catalytic filter filtration efficiency.
|
|
|
Logics transforming the urban waterfront in Bratislava unravelling the decision-making dynamics behind urban development
Machala, Branislav ; Sýkora, Luděk (vedoucí práce) ; Šuška, Pavel (oponent) ; Šourek, Michal (oponent)
Disertační práce se zaměřuje na globalizující se městské nábřeží v Bratislavě. Jejím cílem je odkrývat klíčové logiky, které ovládají procesy rozhodování o revitalizaci břehů řeky v hlavním městě Slovenska. Zaměřuje se na rozšiřování logiky kapitalismu do měnícího se institucionálního prostředí v kontextu postsocialistických transformací. Analytickým rámcem pro tento výzkum je rozlišení mezi kapitalistickými a územními logikami moci (Jessop, 1999; Harvey, 2005). Rozvoj a revitalizace nábřeží se staly globální odpovědí na úpadek vnitřních měst a zasílající konkurenci mezi městy v podmínkách současného neoliberálního kapitalismu. Počátky transformací nábřeží jsou chápány jako jeden z městských projevů geograficky nerovného rozvoje, který byl spuštěn krizí fordismu v Evropě a Severní Americe (Smith, 1990; Jessop, 2000). Kapitál, který se přelévá do zastavěného prostředí (Harvey, 1978; 2005), je v důsledku rozhodnutí uskutečňovaných na více měřítkových úrovních dočasně ukotven na městských nábřežích (Brenner, 2001). Nová budoucnost měst měla být vytvořena pomocí různých typů strategií zaměřených na sociální i fyzické prostředí, jako např. ekosystémový přístup (Laidley, 2007), výstavba ambiciózních megaprojektů a organizace tzv. mega- events (Orueta & Fainstein 2009), nebo lokalizace děl globálních...
|
|
|
DEM-CFD study of flow in a random packed bed
Šourek, M. ; Isoz, Martin
Most catalytic surface reactions as well as other industrial applications take advantage of fixed packed bed reactors. Designers of these reactors rely mostly on empirical formulas derived for various simplifying assumptions, e.g. uniformly distributed porosity. The made simplifications and especially the assumption of uniformly distributed porosity fail if the tube to particle diameter ratio goes under 10 and the „wall effect“ becomes more significant. In such a case, the complete three-dimensional structure of the packed bed has to be considered. Thanks to ongoing improvements in numerical mathematics and computational power, the methods of computational fluid dynamics (CFD) have become a great tool for comprehensive description of the packed beds with low tube to particle diameter ratio. Three-dimensional simulations of the flow through two fixed beds differing in the type of the used particle are presented and compared with available experimental and empirical results. To generate the random fixed beds, we propose a custom approach based on the discrete element method (DEM) code implemented in open-source software Blender. Thereafter, OpenFOAM tools (snappyHexMesh, simpleFoam) are used for creation of the computational mesh and solution of the governing equations describing a single-phase flow in the packed bed.
|
host ::
RSS.