红黑着色的相场两相流并行投影算法
作者:
作者单位:

重庆大学 航空航天学院

中图分类号:

U448.213???????

基金项目:

国家自然科学基金项目(面上项目,重点项目,重大项目);中央高校基本科研业务费


A novel Red-Black Coloring parallel projection algorithm for two-phase flow using the phase field method
Author:
Affiliation:

College of Aerospace Engineering,Chongqing University

Fund Project:

The National Natural Science Foundation of China (General Program, Key Program, Major Research Plan);Fundamental Research Funds for the Central Universities

  • 摘要
  • | |
  • 访问统计
  • |
  • 参考文献 [39]
  • | | | |
  • 文章评论
    摘要:

    本文提出一种交错并行的有限体积投影算法求解基于相场法的两相流控制方程组。该算法的实施主要依赖于压力泊松方程的显式推进设计,从而突破投影算法求解不可压Navier-Stokes方程的效率瓶颈。同时,我们提出了一种交错扫描策略来更新节点上的变量,以实现更紧凑的时空耦合。本算法与相场模型相结合,能够高效、准确地捕获相界面的动态拓扑变化。测试算例结果表明:网格量为131072,采用8线程CPU并行时,新提出并行算法的效率达到串行标准投影算法的80倍以上。

    Abstract:

    In this study, an interleaved parallel finite volume projection algorithm is proposed to solve the two-phase flow control equations based on the phase field method. This strategy relies on the explicit advancement of the pressure Poisson equation, thus breaking through the efficiency bottleneck of the projection algorithm for solving the incompressible Navier-Stokes equations. Moreover, we propose an interleaved scanning strategy to update variables on nodes for compact spatiotemporal coupling. The combination with the phase field method enables accurate capture of interface dynamics and topology at a lower cost. Test results show that with 131072 grids and 8-thread CPU parallelization, the efficiency of the proposed parallel algorithm reaches more than 80 times that of the serial standard projection algorithm.

    参考文献
    [1] BARROS A P, PRAT O P, TESTIK F Y. Size distribution of raindrops[J]. Nature Physics, 2010, 6(4): 232–232.
    [2] O’SULLIVAN J J, NORWOOD E-A, O’MAHONY J A, et al. Atomisation technologies used in spray drying in the dairy industry: a review[J]. Journal of Food Engineering, 2019, 243: 57–69.
    [3] ATENCIA J, BEEBE D J. Controlled microfluidic interfaces[J]. Nature, 2005, 437(7059): 648–655.
    [4] GUERMOND J L, MINEV P, SHEN J. An overview of projection methods for incompressible flows[J]. Computer Methods in Applied Mechanics and Engineering, 2006, 195(44–47): 6011–6045.
    [5] DONG S, SHEN J. A time-stepping scheme involving constant coefficient matrices for phase-field simulations of two-phase incompressible flows with large density ratios[J]. Journal of Computational Physics, 2012, 231(17): 5788–5804.
    [6] DODD M S, FERRANTE A. A fast pressure-correction method for incompressible two-fluid flows[J]. Journal of Computational Physics, 2014, 273: 416–434.
    [7] SHEN J, XU J, YANG J. The scalar auxiliary variable (sav) approach for gradient flows[J]. Journal of Computational Physics, 2018, 353: 407–416.
    [8] LEE C S, HAMON F P, CASTELLETTO N, et al. An aggregation-based nonlinear multigrid solver for two-phase flow and transport in porous media[J]. Computers & Mathematics with Applications, 2022, 113: 282–299.
    [9] LIU T. A nonlinear multigrid method for inverse problem in the multiphase porous media flow[J]. Applied Mathematics and Computation, 2018, 320: 271–281.
    [10] WEYMOUTH G D. Data-driven multi-grid solver for accelerated pressure projection[J]. Computers & Fluids, 2022, 246: 105620.
    [11] FUKA V. PoisFFT – a free parallel fast poisson solver[J]. Applied Mathematics and Computation, 2015, 267: 356–364.
    [12] DOLEAN V, JOLIVET P, NATAF F. An introduction to domain decomposition methods: algorithms, theory, and parallel implementation[M/OL].Philadelphia, PA: Society for Industrial and Applied Mathematics, 2015[2023–02–22].
    [13] DJANALI V, ARMFIELD S W, KIRKPATRICK M P, et al. Parallel speed-up of preconditioned fractional step navier-stokes solvers[J]. Applied Mechanics and Materials, 2014, 493: 215–220.
    [14] ADAM N, FRANKE F, ALAND S. A simple parallel solution method for the navier–stokes cahn–hilliard equations[J]. Mathematics, 2020, 8(8): 1224.
    [15] KUHN M B, DESKOS G, SPRAGUE M A. A mass-momentum consistent coupling for mesh-adaptive two-phase flow simulations[J]. Computers & Fluids, 2022: 105770.
    [16] DING L, SHU C, DING H, et al. Stencil adaptive diffuse interface method for simulation of two-dimensional incompressible multiphase flows[J]. Computers & Fluids, 2010, 39(6): 936–944.
    [17] SAKANE S, TAKAKI T, ROJAS R, et al. Multi-gpus parallel computation of dendrite growth in forced convection using the phase-field-lattice boltzmann model[J]. Journal of Crystal Growth, 2017, 474: 154–159.
    [18] HA S, PARK J, YOU D. A gpu-accelerated semi-implicit fractional-step method for numerical solutions of incompressible navier–stokes equations[J]. Journal of Computational Physics, 2018, 352: 246–264.
    [19] QIU R, HUANG R, YAO X, et al. Physics-informed neural networks for phase-field method in two-phase flow[J]. Physics of Fluids, 2022.
    [20] RAISSI M, PERDIKARIS P, KARNIADAKIS G E. Physics-informed neural networks: a deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations[J]. Journal of Computational Physics, 2019, 378: 686–707.
    [21] DING H, SPELT P D M, SHU C. Diffuse interface model for incompressible two-phase flows with large density ratios[J]. Journal of Computational Physics, 2007, 226(2): 2078–2095.
    [22] LIU H, VALOCCHI A J, ZHANG Y, et al. Lattice boltzmann phase-field modeling of thermocapillary flows in a confined microchannel[J]. Journal of Computational Physics, 2014, 256: 334–356.
    [23] KIM J. A generalized continuous surface tension force formulation for phase-field models for multi-component immiscible fluid flows[J]. Computer Methods in Applied Mechanics and Engineering, 2009, 198(37–40): 3105–3112.
    [24] GOTTLIEB S, SHU C-W. Total variation diminishing runge-kutta schemes[J]. Mathematics of Computation, 1998, 67(221): 73–85.
    [25] FRIEDRICH O. Weighted essentially non-oscillatory schemes for the interpolation of mean values on unstructured grids[R]//JOURNAL OF COMPUTATIONAL PHYSICS.
    [26] MAGNINI M. Characterization of the velocity fields generated by flow initialization in the cfd simulation of multiphase flows[J]. Applied Mathematical Modelling, 2016: 20.
    [27] MIRJALILI S, IVEY C B, MANI A. Comparison between the diffuse interface and volume of fluid methods for simulating two-phase flows[J]. International Journal of Multiphase Flow, 2019, 116: 221–238.
    [28] KIM J. A continuous surface tension force formulation for diffuse-interface models[J]. Journal of Computational Physics, 2005, 204(2): 784–804.
    [29] Mirjalili S, Khanwale M A, Mani A. Assessment of an energy-based surface tension model for simulation of two-phase flows using second-order phase field methods[J]. Journal of Computational Physics, 2023, 474: 111795.
    [30] HYSING S, TUREK S, KUZMIN D, et al. Quantitative benchmark computations of two-dimensional bubble dynamics[J]. International Journal for Numerical Methods in Fluids, 2009, 60(11): 1259–1288.
    [31] KLOSTERMANN J, SCHAAKE K, SCHWARZE R. Numerical simulation of a single rising bubble by vof with surface compression: numerical simulation of a single rising bubble by vof with surface compression[J]. International Journal for Numerical Methods in Fluids, 2013, 71(8): 960–982.
    [32] XIAO Y, ZENG Z, ZHANG L, et al. A spectral element-based phase field method for incompressible two-phase flows[J]. Physics of Fluids, 2022, 34(2): 022114.
    [33] ALAND S, VOIGT A. Benchmark computations of diffuse interface models for two-dimensional bubble dynamics: benchmark computations of diffuse interface models[J]. International Journal for Numerical Methods in Fluids, 2012, 69(3): 747–761.
    [34] MARTIN, MOYCE. Part iv. an experimental study of the collapse of liquid columns on a rigid horizontal plane[J]. Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences, 1952, 244(882): 312–324.
    [35] YU C H, SHEU T W H. Development of a coupled level set and immersed boundary method for predicting dam break flows[J]. Computer Physics Communications, 2017, 221: 1–18.
    [36] MENG W, LIAO L, CHEN M, et al. An enhanced clsvof method with an algebraic second-reconstruction step for simulating incompressible two-phase flows[J]. International Journal of Multiphase Flow, 2022, 154: 104151.
    [37] LI Y L, WAN L, WANG Y H, et al. Numerical investigation of interface capturing method by the rayleigh-taylor instability, dambreak and solitary wave problems[J]. Ocean Engineering, 2019, 194: 106583.
    [38] KAMRA M M, MOHD N, LIU C, et al. Numerical and experimental investigation of three-dimensionality in the dam-break flow against a vertical wall[J]. Journal of Hydrodynamics, 2018, 30(4): 682–693.
    [39] KAMRA M M, MOHD N, LIU C, et al. Numerical and experimental investigation of three-dimensionality in the dam-break flow against a vertical wall[J]. Journal of Hydrodynamics, 2018, 30(4): 682–693.
    相似文献
    引证文献
引用本文
分享
文章指标
  • 点击次数:284
  • 下载次数: 0
  • HTML阅读次数: 0
  • 引用次数: 0
历史
  • 收稿日期:2023-03-02
  • 最后修改日期:2023-03-13
  • 录用日期:2023-03-13
文章二维码