Fully-implicit orthogonal reconstructed Discontinuous Galerkin method for fluid dynamics with phase change

详细信息    查看全文

• 作者：R. Nourgaliev^a ; ^{nourgaliev1@llnl.gov" class="auth_mail" title="E-mail the corresponding author} ; H. Luo^b ; B. Weston^c ; A. Anderson^a ; S. Schofield^a ; T. Dunn^a ; J.-P. Delplanque^c
• 关键词：Discontinuous Galerkin method ; High-order space&ndash ; time discretization ; Fully-implicit solvers ; Newton&ndash ; Krylov algorithm ; Fluid dynamics ; Phase-change
• 刊名：Journal of Computational Physics
• 出版年：2016
• 出版时间：15 January 2016
• 年：2016
• 卷：305
• 期：Complete
• 页码：964-996
• 全文大小：17852 K

文摘

A new reconstructed Discontinuous Galerkin (rDG) method, based on orthogonal basis/test functions, is developed for fluid flows on unstructured meshes. Orthogonality of basis functions is essential for enabling robust and efficient fully-implicit Newton–Krylov based time integration. The method is designed for generic partial differential equations, including transient, hyperbolic, parabolic or elliptic operators, which are attributed to many multiphysics problems. We demonstrate the method's capabilities for solving compressible fluid–solid systems (in the low Mach number limit), with phase change (melting/solidification), as motivated by applications in Additive Manufacturing (AM). We focus on the method's accuracy (in both space and time), as well as robustness and solvability of the system of linear equations involved in the linearization steps of Newton-based methods. The performance of the developed method is investigated for highly-stiff problems with melting/solidification, emphasizing the advantages from tight coupling of mass, momentum and energy conservation equations, as well as orthogonality of basis functions, which leads to better conditioning of the underlying (approximate) Jacobian matrices, and rapid convergence of the Krylov-based linear solver.