{"title":"Monolithic Newton-Multigrid Solver for Multiphase Flow Problems with Surface Tension","authors":"M. Afaq, S. Turek, A. Ouazzi, Arooj Fatima","doi":"10.4995/yic2021.2021.12390","DOIUrl":null,"url":null,"abstract":"We have developed a monolithic Newton-multigrid solver for multiphase flow problems which solves velocity, pressure and interface position simultaneously. The main idea of our work is based on the formulations discussed in [1], where it points out the feasibility of a fully implicit monolithic solver for multiphase flow problems via two formulations, a curvature-free level set approach and a curvature-free cutoff material function approach. Both formulations are fully implicit and have the advantages of requiring less regularity, since neither normals nor curvature are explicitly calculated, and no capillary time restriction. Furthermore, standard Navier-Stokes solvers might be used, which do not have to take into account inhomogeneous force terms. The reinitialization issue is integrated with a nonlinear terms within the formulations.The nonlinearity is treated with a Newton-type solver with divided difference evaluation of the Jacobian matrices. The resulting linearized system inside of the outer Newton solver is a typical saddle point problem which is solved using the geometrical multigrid with Vanka-like smoother using higher order stable FEM pair $Q_2/P^{\\text{disc}}_1$ for velocity and pressure and $Q_2$ for all other variables. The method is implemented into an existing software packages for the numerical simulation of multiphase flows (FeatFlow). The robustness and accuracy of this solver is tested for two different test cases, i.e. static bubble and oscillating bubble, respectively [2].REFERENCES[1] Ouazzi, A., Turek, S. and Damanik, H. A curvature-free multiphase flow solver via surface stress-based formulation. Int. J. Num. Meth. Fluids., Vol. 88, pp. 18–31, (2018).[2] Afaq, M. A., Turek, S., Ouazzi, A. and Fatima, A. Monolithic Newton-Multigrid Solver for Multiphase Flow Problems with Surface Tension. Ergebnisberichte des Instituts fuer Angewandte Mathematik Nummer 636, Fakultaet fuer Mathematik, TU Dortmund University, 636, 2021.","PeriodicalId":406819,"journal":{"name":"Proceedings of the YIC 2021 - VI ECCOMAS Young Investigators Conference","volume":"12 1","pages":"0"},"PeriodicalIF":0.0000,"publicationDate":"2021-07-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Proceedings of the YIC 2021 - VI ECCOMAS Young Investigators Conference","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.4995/yic2021.2021.12390","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
We have developed a monolithic Newton-multigrid solver for multiphase flow problems which solves velocity, pressure and interface position simultaneously. The main idea of our work is based on the formulations discussed in [1], where it points out the feasibility of a fully implicit monolithic solver for multiphase flow problems via two formulations, a curvature-free level set approach and a curvature-free cutoff material function approach. Both formulations are fully implicit and have the advantages of requiring less regularity, since neither normals nor curvature are explicitly calculated, and no capillary time restriction. Furthermore, standard Navier-Stokes solvers might be used, which do not have to take into account inhomogeneous force terms. The reinitialization issue is integrated with a nonlinear terms within the formulations.The nonlinearity is treated with a Newton-type solver with divided difference evaluation of the Jacobian matrices. The resulting linearized system inside of the outer Newton solver is a typical saddle point problem which is solved using the geometrical multigrid with Vanka-like smoother using higher order stable FEM pair $Q_2/P^{\text{disc}}_1$ for velocity and pressure and $Q_2$ for all other variables. The method is implemented into an existing software packages for the numerical simulation of multiphase flows (FeatFlow). The robustness and accuracy of this solver is tested for two different test cases, i.e. static bubble and oscillating bubble, respectively [2].REFERENCES[1] Ouazzi, A., Turek, S. and Damanik, H. A curvature-free multiphase flow solver via surface stress-based formulation. Int. J. Num. Meth. Fluids., Vol. 88, pp. 18–31, (2018).[2] Afaq, M. A., Turek, S., Ouazzi, A. and Fatima, A. Monolithic Newton-Multigrid Solver for Multiphase Flow Problems with Surface Tension. Ergebnisberichte des Instituts fuer Angewandte Mathematik Nummer 636, Fakultaet fuer Mathematik, TU Dortmund University, 636, 2021.
我们开发了一种求解多相流问题的单片牛顿-多网格求解器,可以同时求解速度、压力和界面位置。我们工作的主要思想是基于[1]中讨论的公式,其中指出了通过两种公式(无曲率水平集方法和无曲率截止材料函数方法)求解多相流问题的全隐式单片求解器的可行性。这两个公式都是完全隐式的,并且具有要求较少规则性的优点,因为法线和曲率都没有显式计算,也没有毛细时间限制。此外,可以使用标准的Navier-Stokes解,它不必考虑非齐次力项。重新初始化问题与公式中的非线性项相结合。用雅可比矩阵的分差求值的牛顿型求解器处理非线性问题。所得到的外牛顿求解器内部线性化系统是一个典型的鞍点问题,该问题采用具有Vanka-like光滑的几何多重网格,采用高阶稳定有限元对$Q_2/P^{\text{disc}}_1$表示速度和压力,$Q_2$表示其他变量。该方法已在现有的多相流数值模拟软件包(FeatFlow)中实现。通过静态气泡和振荡气泡两种不同的测试用例对求解器的鲁棒性和准确性进行了测试[2]。[1]刘建军,刘建军,刘建军,等。一种基于表面应力的无曲率多相流求解方法。Int。J. Num.冰毒。液体。, Vol. 88, pp. 18-31, (2018).[2]Afaq, M. A, Turek, S., Ouazzi, A.和Fatima, A.考虑表面张力的多相流问题的单片牛顿-多网格求解器。德国多特蒙德大学数学研究所(第636期),德国多特蒙德大学数学研究所(第636,2021)。