Parallel simulations for a 2D x/z two-phase flow fluid-solid particle model

•A parallel algorithm for a 2D x/z two-phase Eulerian approach is presented.•The numerical scheme is based on a projection technique and a finite volume method.•The parallelization is designed using block domain decomposition.•Convergence of the SOR method and the Bi-CGSTAB method are compared.•The...

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Computers & fluids Ročník 173; s. 103 - 110
Hlavní autori: Uh Zapata, M., Pham Van Bang, D., Nguyen, K.D.
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Amsterdam Elsevier Ltd 15.09.2018
Elsevier BV
Predmet:
ISSN:0045-7930, 1879-0747
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí:•A parallel algorithm for a 2D x/z two-phase Eulerian approach is presented.•The numerical scheme is based on a projection technique and a finite volume method.•The parallelization is designed using block domain decomposition.•Convergence of the SOR method and the Bi-CGSTAB method are compared.•The erosion of granular beds is adopted to verify the algorithm performance. The importance of the two-phase flow model relies upon the correct simulation of specific problems in which the passive tracer model fails; however, the mathematical and numerical models need to be improved in order to reproduce fluid-solid particle interactions of greater complexity. Such is the case with regard to simulating erosion patterns suffered upon horizontal granular beds by means of a vertical water jet as shown here. Moreover, sequential platforms have proven to be insufficient in providing the required computational power needed to obtain fast and detailed simulations. In this paper, a fully parallel algorithm for a two-dimensional x/z two-phase Eulerian approach is presented and applied for the numerical solution of the erosion of sediment beds. The parallelization is designed by a row and column block domain decomposition technique using a distributed memory platform with Message Passing Interface (MPI). Arising from the numerical method, a Poisson problem for the pressure is solved at each time step. The discretization results in a non-symmetric variable-coefficient linear system which is solved using several parallel Successive Over-Relaxation (SOR) algorithms, including partitioning and coloring methods. The specification of the optimal relaxation parameter to achieve efficiency is found numerically. Results show that SOR methods achieve faster convergence rates and the simulation time achieves an order similar to that required for the typical, widely-used Bi-Conjugate Gradient Stabilized (Bi-CGSTAB) method. The performances of the algorithms are evaluated in terms of speedup and efficiency. The results indicate that the parallel code significantly improves the results of the sequential calculation in general.
Bibliografia:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:0045-7930
1879-0747
DOI:10.1016/j.compfluid.2018.03.019