Browse ORBi

- What it is and what it isn't
- Green Road / Gold Road?
- Ready to Publish. Now What?
- How can I support the OA movement?
- Where can I learn more?

ORBi

Local Verlet buffer approach for broad-phase interaction detection in Discrete Element Method Mainassara Chekaraou, Abdoul Wahid ; Besseron, Xavier ; Rousset, Alban et al in Computer Physics Communications (in press) The Extended Discrete Element Method (XDEM) is a novel and innovative numerical simulation technique that extends the dynamics of granular materials or particles as described through the classical ... [more ▼] The Extended Discrete Element Method (XDEM) is a novel and innovative numerical simulation technique that extends the dynamics of granular materials or particles as described through the classical discrete element method (DEM) by additional properties such as the thermodynamic state, stress/strain for each particle. Such DEM simulations used by industries to set up their experimental processes are complexes and heavy in computation time. Those simulations perform at each time step a collision detection to generate a list of interacting particles that is one of the most expensive computation parts of a DEM simulation. The Verlet buffer method, which was first introduced in Molecular Dynamic (MD) (and is also used in DEM) allows to keep the interaction list for many time step by extending each particle neighborhood by a certain extension range, and thus broadening the interaction list. The method relies mainly on the stability of the DEM, which ensures that no particles move erratically or unpredictably from one time step to the next: this is called temporal coherency. In the classical and current approach, all the particles have their neighborhood extended by the same value which leads to suboptimal performances in simulations where different flow regimes coexist. Additionally, and unlike in MD (which remains very different from DEM on several aspects), there is no comprehensive study analyzing the different parameters that affect the performance of the Verlet buffer method in DEM. In this work, we apply a dynamic neighbor list update method that depends on the particles' individual displacement, and an extension range specific to each particle and based on their local flow regime for the generation of the neighbor list. The update of the interaction list is analyzed throughout the simulation based on the displacement of the particle allowing a flexible update according to the flow regime conditions. We evaluate the influence of the Verlet extension range on the performance of the execution time through different test cases and we empirically analyze and define the extension range value giving the minimum of the global simulation time. [less ▲] Detailed reference viewed: 77 (36 UL)Process analysis in thermal process engineering with high-performance computing using the example of grate firing Peters, Bernhard ; Rousset, Alban ; Besseron, Xavier et al in Process analysis in thermal process engineering with high- performance computing using the example of grate firing (in press) Biomass as a renewable energy source continues to grow in popularity to reduce fossil fuel consumption for environmental and economic benefits. In the present contribution, the combustion chamber of a 16 ... [more ▼] Biomass as a renewable energy source continues to grow in popularity to reduce fossil fuel consumption for environmental and economic benefits. In the present contribution, the combustion chamber of a 16 MW geothermal steam super-heater, which is part of the Enel Green Power "Cornia 2" power plant, is being investigated with high-performance computing methods. For this purpose, the extended discrete element method (XDEM) developed at the University of Luxembourg is used in a high-performance computing environment, which includes both the moving wooden bed and the combustion chamber above it. The XDEM simulation platform is based on a hybrid four-way coupling between the Discrete Element Method (DEM) and Computational Fluid Dynamics (CFD). In this approach, particles are treated as discrete elements that are coupled by heat, mass, and momentum transfer to the surrounding gas as a continuous phase. For individual wood particles, besides the equations of motion, the differential conservation equations for mass, heat, and momentum are solved, which describe the thermodynamic state during thermal conversion. The consistency of the numerical results with the actual system performance is discussed in this paper to determine the potentials and limitations of the approach. [less ▲] Detailed reference viewed: 65 (5 UL)Politiques de synchronisation dans les systèmes multi-agents distribués parallèles Rousset, Alban ; ; et al in Rousset, Alban (Ed.) Politiques de synchronisation dans les systèmes multi-agents distribués parallèles (2020, June 29) Among simulation or modelisation methods, multi-agent systems are an interesting candidate to simulate complex systems. As the size of the model's increases, the use of parallel multi-agent systems is ... [more ▼] Among simulation or modelisation methods, multi-agent systems are an interesting candidate to simulate complex systems. As the size of the model's increases, the use of parallel multi-agent systems is mandatory but comes with many issues. In this article, we are interested in the impact of synchronization on models design and on their execution. We highlight synchronization problems through model instances then we experimentally analyze the impact of synchronization on large scale test cases. [less ▲] Detailed reference viewed: 33 (0 UL)Detailed Numerical Three-dimensional and Transient Analysis of a Grate Firing Combustion Process by Innovative High Performance Computing Mainassara Chekaraou, Abdoul Wahid ; Rousset, Alban ; Besseron, Xavier et al Scientific Conference (2020, April 27) Detailed reference viewed: 51 (4 UL)6-way coupling of DEM+CFD+FEM with preCICE Besseron, Xavier ; Rousset, Alban ; et al Presentation (2020, February) In this work, we present our preliminary results on the 6-way coupling of 3 numerical solvers: XDEM for the Discrete Element Method (DEM), OpenFOAM for Computation Fluid Dynamics (CFD), and deal.II for ... [more ▼] In this work, we present our preliminary results on the 6-way coupling of 3 numerical solvers: XDEM for the Discrete Element Method (DEM), OpenFOAM for Computation Fluid Dynamics (CFD), and deal.II for Finite Element Method (FEM). We relied on the existing preCICE adapters for OpenFOAM and deal.II and we have implemented a new preCICE adapter for the eXtended Discrete Element Method (XDEM), an innovative DEM software developed at the University of Luxembourg. The XDEM adapter permits coupling of the particulate phase of DEM with CFD and FEM: - DEM+FEM is a surface coupling that performs the exchange of surface forces and displacement between the particles and a deformable solid; - DEM+CFD is a volume coupling that performs the exchange of porosity, momentum, drag force and buoyancy between the particles and the fluid. Put together with the pre-existing CFD+FEM coupling, we obtain a 6-way coupled multi-physics solver for particles, fluid and deformable solids. We have tested and evaluated our multi-physics solver on the tutorial case “Cylinder with a flap” derived from the benchmarking case of Turek and Hron, that we extended to include a particulate phase solved by XDEM. [less ▲] Detailed reference viewed: 138 (11 UL)High Performance Parallel Coupling of OpenFOAM+XDEM Besseron, Xavier ; ; Rousset, Alban et al Presentation (2019, June 21) Detailed reference viewed: 191 (19 UL)Co-located Partitioning Strategy and Dual-grid Multiscale Approach for Parallel Coupling of CFD-DEM Simulations Besseron, Xavier ; Pozzetti, Gabriele ; Rousset, Alban et al Presentation (2019, June 05) Detailed reference viewed: 56 (7 UL)A parallel dual-grid multiscale approach to CFD-DEM couplings Pozzetti, Gabriele ; ; Besseron, Xavier et al in Journal of Computational Physics (2019), 378 In this work, a new parallel dual-grid multiscale approach for CFD-DEM couplings is investigated. Dual- grid multiscale CFD-DEM couplings have been recently developed and successfully adopted in different ... [more ▼] In this work, a new parallel dual-grid multiscale approach for CFD-DEM couplings is investigated. Dual- grid multiscale CFD-DEM couplings have been recently developed and successfully adopted in different applications still, an efficient parallelization for such a numerical method represents an open issue. Despite its ability to provide grid convergent solutions and more accurate results than standard CFD-DEM couplings, this young numerical method requires good parallel performances in order to be applied to large-scale problems and, therefore, extend its range of application. The parallelization strategy here proposed aims to take advantage of the enhanced complexity of a dual-grid coupling to gain more flexibility in the domain partitioning while keeping a low inter-process communication cost. In particular, it allows avoiding inter- process communication between CFD and DEM software and still allows adopting complex partitioning strategies thanks to an optimized grid-based communication. It is shown how the parallelized multiscale coupling holds all its natural advantages over a mono-scale coupling and can also have better parallel performance. Three benchmark cases are presented to assess the accuracy and performance of the strategy. It is shown how the proposed method allows maintaining good parallel performance when operated over 1000 processes. [less ▲] Detailed reference viewed: 117 (27 UL)Verlet buffer for broad phase interaction detection in Discrete Element Method Mainassara Chekaraou, Abdoul Wahid ; Rousset, Alban ; Besseron, Xavier et al Poster (2018, September 24) The Extended Discrete Element Method (XDEM) is a novel and innovative numerical simulation technique that extends the dynamics of granular materials or particles as described through the classical ... [more ▼] The Extended Discrete Element Method (XDEM) is a novel and innovative numerical simulation technique that extends the dynamics of granular materials or particles as described through the classical discrete element method (DEM) by additional properties such as the thermodynamic state, stress/strain for each particle. Such DEM simulations used by industries to set up their experimental processes are complexes and heavy in computation time. Therefore, simulations have to be precise, efficient and fast in order to be able to process hundreds of millions of particles. To tackle this issue, such DEM simulations are usually parallelized with MPI. One of the most expensive computation parts of a DEM simulation is the collision detection of particles. It is classically divided into two steps: the broad phase and the narrow phase. The broad phase uses simplified bounding volumes to perform an approximated but fast collision detection. It returns a list of particle pairs that could interact. The narrow phase is applied to the result of the broad phase and returns the exact list of colliding particles. The goal of this research is to apply a Verlet buffer method to (X)DEM simulations regardless of which broad phase algorithm is used. We rely on the fact that such DEM simulations are temporal coherent: the neighborhood only changes slightly from the last time-step to the current time-step. We use the Verlet buffer method to extend the list of pairs returned by the broad phase by stretching the particles bounding volume with an extension range. This allows re-using the result of the broad phase for several time-steps before an update is required once again and thereby its reduce the number of times the broad phase is executed. We have implemented a condition based on particles displacements to ensure the validity of the broad phase: a new one is executed to update the list of colliding particles only when necessary. This guarantees identical results because approximations introduced in the broad phase by our approach are corrected in the narrow phase which is executed at every time-steps anyway. We perform an extensive study to evaluate the influence of the Verlet extension range on the performance of the execution in terms of computation time and memory consumption. We consider different test-cases, partitioners (ORB, Zoltan, METIS, SCOTCH, ...), broad phase algorithms (Link cell, Sweep and prune, ...) and grid configurations (fine, coarse), sequential and parallel (up to 280 cores). While a larger Verlet buffer increases the cost of the broad phase and narrow phase, it also allows skipping a significant number of broad phase execution (> 99 \%). As a consequence, our first results show that this approach can speeds up the total .execution time up to a factor of 5 for sequential executions, and up to a factor of 3 parallel executions on 280 cores while maintaining a reasonable memory consumption. [less ▲] Detailed reference viewed: 96 (23 UL)A co-located partitions strategy for parallel CFD-DEM couplings Pozzetti, Gabriele ; Besseron, Xavier ; Rousset, Alban et al in Advanced Powder Technology (2018) In this work, a new partition-collocation strategy for the parallel execution of CFD–DEM couplings is investigated. Having a good parallel performance is a key issue for an Eulerian-Lagrangian software ... [more ▼] In this work, a new partition-collocation strategy for the parallel execution of CFD–DEM couplings is investigated. Having a good parallel performance is a key issue for an Eulerian-Lagrangian software that aims to be applied to solve industrially significant problems, as the computational cost of these couplings is one of their main drawback. The approach presented here consists in co-locating the overlapping parts of the simulation domain of each software on the same MPI process, in order to reduce the cost of the data exchanges. It is shown how this strategy allows reducing memory consumption and inter-process communication between CFD and DEM to a minimum and therefore to overcome an important parallelization bottleneck identified in the literature. Three benchmarks are proposed to assess the consistency and scalability of this approach. A coupled execution on 280 cores shows that less than 0.1% of the time is used to perform inter-physics data exchange. [less ▲] Detailed reference viewed: 145 (29 UL)Hybrid MPI+OpenMP Implementation of eXtended Discrete Element Method Mainassara Chekaraou, Abdoul Wahid ; Rousset, Alban ; Besseron, Xavier et al in Proc. of the 9th Workshop on Applications for Multi-Core Architectures (WAMCA'18), part of 30th Intl. Symp. on Computer Architecture and High Performance Computing (SBAC-PAD 2018) (2018, September) The Extended Discrete Element Method (XDEM) is a novel and innovative numerical simulation technique that ex- tends classical Discrete Element Method (DEM) (which simulates the motion of granular material ... [more ▼] The Extended Discrete Element Method (XDEM) is a novel and innovative numerical simulation technique that ex- tends classical Discrete Element Method (DEM) (which simulates the motion of granular material), by additional properties such as the chemical composition, thermodynamic state, stress/strain for each particle. It has been applied successfully to numerous industries involving the processing of granular materials such as sand, rock, wood or coke [16], [17]. In this context, computational simulation with (X)DEM has become a more and more essential tool for researchers and scientific engineers to set up and explore their experimental processes. However, increasing the size or the accuracy of a model requires the use of High Performance Computing (HPC) platforms over a parallelized implementation to accommodate the growing needs in terms of memory and computation time. In practice, such a parallelization is traditionally obtained using either MPI (distributed memory computing), OpenMP (shared memory computing) or hybrid approaches combining both of them. In this paper, we present the results of our effort to implement an OpenMP version of XDEM allowing hybrid MPI+OpenMP simulations (XDEM being already parallelized with MPI). Far from the basic OpenMP paradigm and recommendations (which simply summarizes by decorating the main computation loops with a set of OpenMP pragma), the OpenMP parallelization of XDEM required a fundamental code re-factoring and careful tuning in order to reach good performance. There are two main reasons for those difficulties. Firstly, XDEM is a legacy code devel- oped for more than 10 years, initially focused on accuracy rather than performance. Secondly, the particles in a DEM simulation are highly dynamic: they can be added, deleted and interaction relations can change at any timestep of the simulation. Thus this article details the multiple layers of optimization applied, such as a deep data structure profiling and reorganization, the usage of fast multithreaded memory allocators and of advanced process/thread-to-core pinning techniques. Experimental results evaluate the benefit of each optimization individually and validate the implementation using a real-world application executed on the HPC platform of the University of Luxembourg. Finally, we present our Hybrid MPI+OpenMP results with a 15%-20% performance gain and how it overcomes scalability limits (by increasing the number of compute cores without dropping of performances) of XDEM-based pure MPI simulations. [less ▲] Detailed reference viewed: 176 (35 UL)Parallel Coupling of CFD-DEM simulations Besseron, Xavier ; Pozzetti, Gabriele ; Rousset, Alban et al Presentation (2018, August 20) Detailed reference viewed: 65 (2 UL)Parallel Coupling of CFD-DEM simulations Pozzetti, Gabriele ; Besseron, Xavier ; Rousset, Alban et al Presentation (2018, August) Eulerian-Lagrangian couplings are nowadays widely used to address engineering and technical problems. In particular, CFD-DEM couplings have been successfully applied to study several configurations ... [more ▼] Eulerian-Lagrangian couplings are nowadays widely used to address engineering and technical problems. In particular, CFD-DEM couplings have been successfully applied to study several configurations ranging from mechanical, to chemical and environmental engineering. However, such simulations are normally very computationally intensive, and the execution time represents a major issue for the applicability of this numerical approach to complex scenarios. With this work, we introduce a novel coupling approach aiming at improving the performance of the parallel CFD-DEM simulations. This strategy relies on two points. First, we propose a new partition-collocation strategy for the parallel execution of CFD–DEM couplings, which can considerably reduce the amount of inter-process communication between the CFD and DEM parts. However, this strategy imposes some alignment constraints on the CFD mesh. Secondly, we adopt a dual-grid multiscale scheme for the CFD-DEM coupling, that is known to offer better numerical properties, and that allows us to obtain more flexibility on the domain partitioning overcoming the alignment constraints. We assess the correctness and performance of our approach on elementary benchmarks and at a large scale with a realistic test-case. The results show a significant performance improvement compared to other state-of-art CFD-DEM couplings presented in the literature. [less ▲] Detailed reference viewed: 121 (13 UL)A Parallel Multiscale DEM-VOF Method For Large-Scale Simulations Of Three-Phase Flows Pozzetti, Gabriele ; Besseron, Xavier ; Rousset, Alban et al in Proceedings of ECCM-ECFD 2018 (2018) A parallel dual-grid multiscale DEM-VOF coupling is here investigated. Dual- grid multiscale couplings have been recently used to address different engineering problems involving the interaction between ... [more ▼] A parallel dual-grid multiscale DEM-VOF coupling is here investigated. Dual- grid multiscale couplings have been recently used to address different engineering problems involving the interaction between granular phases and complex fluid flows. Nevertheless, previous studies did not focus on the parallel performance of such a coupling and were, therefore, limited to relatively small applications. In this contribution, we propose an insight into the performance of the dual-grid multiscale DEM-VOF method for three- phase flows when operated in parallel. In particular,we focus on a famous benchmark case for three-phase flows and assess the influence of the partitioning algorithm on the scalability of the dual-grid algorithm. [less ▲] Detailed reference viewed: 91 (25 UL)Nested Graphs: a model to efficiently distribute multi-agent systems on HPC clusters Rousset, Alban ; ; et al in Concurrency & Computation : Practice & Experience (2017), 30(7), 22 Computational simulation is becoming increasingly important in numerous research fields. Depending on the modeled system, several methods such as differential equations or Monte-Carlo simulations may be ... [more ▼] Computational simulation is becoming increasingly important in numerous research fields. Depending on the modeled system, several methods such as differential equations or Monte-Carlo simulations may be used to represent the system behavior. The amount of computation and memory needed to run a simulation depends on its size and precision and large simulations usually lead to long runs thus requiring to adapt the model to a parallel system. Complex systems are often simulated using Multi-agent systems (MAS). While linear system based models benefit from a large set of tools to take advantage of parallel resources, multi-agent systems suffer from a lack of platforms that ease the use of such resources. In this paper, we propose the use of Nested Graphs for a new modeling approach that allows the design of large, complex and multi-scale multi-agent models which can efficiently be distributed on parallel resources. Nested Graphs are formally defined and are illustrated on the well-known predator-prey model. We also introduce PDMAS (Parallel and Distributed Multi-Agent System) a platform that implements the Nested Graph modeling approach to ease the distribution of multi-agent models on High Performance Computing clusters. Performance results are presented to validate the efficiency of the resulting models. [less ▲] Detailed reference viewed: 62 (2 UL)On the performance of an overlapping-domain parallelization strategy for Eulerian-Lagrangian Multiphysics software Pozzetti, Gabriele ; Besseron, Xavier ; Rousset, Alban et al in AIP Conference Proceedings ICNAAM 2017 (2017, September) In this work, a strategy for the parallelization of a two-way CFD-DEM coupling is investigated. It consists on adopting balanced overlapping partitions for the CFD and the DEM domains, that aims to reduce ... [more ▼] In this work, a strategy for the parallelization of a two-way CFD-DEM coupling is investigated. It consists on adopting balanced overlapping partitions for the CFD and the DEM domains, that aims to reduce the memory consumption and inter-process communication between CFD and DEM. Two benchmarks are proposed to assess the consistency and scalability of this approach, coupled execution on 252 cores shows that less than 1\% of time is used to perform inter-physics data exchange. [less ▲] Detailed reference viewed: 213 (61 UL)Parallelizing XDEM: Load-balancing policies and efficiency, a study Rousset, Alban ; Besseron, Xavier ; Peters, Bernhard Scientific Conference (2017, September) In XDEM, the simulation domain is geometrically decomposed in regular fixed-size cells that are used to distribute the workload between the processes. The role of the partitioning algorithm is to ... [more ▼] In XDEM, the simulation domain is geometrically decomposed in regular fixed-size cells that are used to distribute the workload between the processes. The role of the partitioning algorithm is to distribute the cells among all the processes in order to balance the workload. To accomplish this task, the partitioning algorithm relies on a computing/communication cost that has been estimated for each cell. A proper estimation of these costs is fundamental to obtain pertinent results during this phase. The study in the work is twofold. First, we integrate five partitioning algorithms (ORB, RCB, RIB, kway and PhG) in the XDEM framework [1]. Most of these algorithms are implemented within the Zoltan library [2], a parallel framework for partitioning and ordering problems. Secondly, we propose different policies to estimate the computing cost and communication cost of the different cells composing the simulation domain. Then, we present an experimental evaluation and a performance comparison of these partitioning algorithms and cost-estimation policies on a large scale parallel execution of XDEM running on the HPC platform of the University of Luxembourg. Finally, after explaining the pros and cons of each partitioning algorithms and cost-estimation policies, we discuss on the best choices to adopt depending on the simulation case. [less ▲] Detailed reference viewed: 108 (21 UL)Comparing Broad-Phase Interaction Detection Algorithms for Multiphysics DEM Applications Rousset, Alban ; Mainassara Chekaraou, Abdoul Wahid ; Liao, Yu-Chung et al in AIP Conference Proceedings ICNAAM 2017 (2017, September) Collision detection is an ongoing source of research and optimization in many fields including video-games and numerical simulations [6, 7, 8]. The goal of collision detection is to report a geometric ... [more ▼] Collision detection is an ongoing source of research and optimization in many fields including video-games and numerical simulations [6, 7, 8]. The goal of collision detection is to report a geometric contact when it is about to occur or has actually occurred. Unfortunately, detailed and exact collision detection for large amounts of objects represent an immense amount of computations, naively n 2 operation with n being the number of objects [9]. To avoid and reduce these expensive computations, the collision detection is decomposed in two phases as it shown on Figure 1: the Broad-Phase and the Narrow-Phase. In this paper, we focus on Broad-Phase algorithm in a large dynamic three-dimensional environment. We studied two kinds of Broad-Phase algorithms: spatial partitioning and spatial sorting. Spatial partitioning techniques operate by dividing space into a number of regions that can be quickly tested against each object. Two types of spatial partitioning will be considered: grids and trees. The grid-based algorithms consist of a spatial partitioning processing by dividing space into regions and testing if objects overlap the same region of space. And this reduces the number of pairwise to test. The tree-based algorithms use a tree structure where each node spans a particular space area. This reduces the pairwise checking cost because only tree leaves are checked. The spatial sorting based algorithm consists of a sorted spatial ordering of objects. Axis-Aligned Bounding Boxes (AABBs) are projected onto x, y and z axes and put into sorted lists. By sorting projection onto axes, two objects collide if and only if they collide on the three axes. This axis sorting reduces the number of pairwise to tested by reducing the number of tests to perform to only pairs which collide on at least one axis. For this study, ten different Broad-Phase collision detection algorithms or framework have been considered. The Bullet [6], CGAL [10, 11] frameworks have been used. Concerning the implemented algorithms most of them come from papers or given implementation. [less ▲] Detailed reference viewed: 252 (54 UL) |
||