Skip to main content
OSTI.GOVU.S. Department of Energy Office of Scientific and Technical Information
U.S. Department of Energy
Office of Scientific and Technical Information
 

Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics

Journal Article · · Journal of Computational Physics
 [1];  [1];  [2];  [1]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Univ. of California, Davis, CA (United States)
+ Show Author Affiliations
n this paper, we introduce a multigrid block-based preconditioner for solving linear systems arising from a Discontinuous Galerkin discretization of the all-speed Navier-Stokes equations with phase change. The equations are discretized in conservative form with a reconstructed Discontinuous Galerkin (rDG) method and integrated with fully-implicit time discretization schemes. To robustly converge the numerically stiff systems, we use the Newton-Krylov framework with a primitive-variable formulation (pressure, velocity, and temperature), which is better conditioned than the conservative-variable form at low-Mach number. In the limit of large acoustic CFL number and viscous Fourier number, there is a strong coupling between the velocity-pressure system and the linear systems become non-diagonally dominant. To effectively solve these ill-conditioned systems, an approximate block factorization preconditioner is developed, which uses the Schur complement to reduce a 3 x 3 block system into a sequence of two 2 x 2 block systems: velocity-pressure,vP, and velocity-temperature, vT. We compare the performance of the vP-vT Schur complement preconditioner to classic preconditioning strategies: monolithic algebraic multigrid (AMG), element-block SOR, and primitive variable block Gauss-Seidel. The performance of the preconditioned solver is investigated in the limit of large CFL and Fourier numbers for low-Mach lid-driven cavity flow, Rayleigh-Bénard melt convection, compressible internally heated convection, and 3D laser-induced melt pool flow. Here, numerical results demonstrate that the vP-vT Schur complement preconditioned solver scales well both algorithmically and in parallel, and is robust for highly ill-conditioned systems, for all tested rDG discretization schemes (up to 4th-order).
Research Organization:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Organization:
USDOE; USDOE National Nuclear Security Administration (NNSA)
Grant/Contract Number:
AC52-07NA27344
OSTI ID:
1577941
Alternate ID(s):
OSTI ID: 1701812
Report Number(s):
LLNL-JRNL--745515; 900481
Journal Information:
Journal of Computational Physics, Journal Name: Journal of Computational Physics Journal Issue: C Vol. 397; ISSN 0021-9991
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English

References (56)

On a multilevel preconditioning module for unstructured mesh Krylov solvers: two-level Schwarz journal March 2002
Modelling liquid-solid phase changes with melt convection journal August 1989
Performance of fully coupled algebraic multilevel domain decomposition preconditioners for incompressible flow and transport journal January 2006
Numerical simulations of a transient injection flow at low Mach number regime journal October 2008
The Application of Preconditioning in Viscous Flows journal April 1993
A New Flux Splitting Scheme journal July 1993
A Sequel to AUSM: AUSM+ journal December 1996
A Numerical Method for Solving Incompressible Viscous Flow Problems journal August 1997
Physics-Based Preconditioning and the Newton–Krylov Method for Non-equilibrium Radiation Diffusion journal May 2000
Preconditioning Techniques for Large Linear Systems: A Survey journal November 2002
High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics journal November 2018
Jacobian–Free Newton–Krylov Methods for the Accurate Time Integration of Stiff Wave Systems journal October 2005
Simulation of time-dependent pool shape during laser spot welding: Transient effects journal December 2003
A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows journal October 1972
A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems journal August 1987
High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method journal December 1982
Preconditioned methods for solving the incompressible and low speed compressible equations journal October 1987
High-resolution viscous flow simulations at arbitrary Mach number journal January 2003
A parallel block multi-level preconditioner for the 3D incompressible Navier–Stokes equations journal May 2003
The extension of incompressible flow solvers to the weakly compressible regime journal February 2003
On the behaviour of upwind schemes in the low Mach number limit journal January 1999
BoomerAMG: A parallel algebraic multigrid solver and preconditioner journal April 2002
Laser powder-bed fusion additive manufacturing: Physics of complex melt flow and formation mechanisms of pores, spatter, and denudation zones journal April 2016
A finite element method for low-speed compressible flows journal May 2004
A set of parallel, implicit methods for a reconstructed discontinuous Galerkin method for compressible flows on 3D hybrid grids journal July 2014
Jacobian-free Newton–Krylov methods: a survey of approaches and applications journal January 2004
Performance of fully coupled domain decomposition preconditioners for finite element transport/reaction simulations journal May 2005
A sequel to AUSM, Part II: AUSM+-up for all speeds journal May 2006
Development of a 2-D algorithm to simulate convection and phase transition efficiently journal November 2006
A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier–Stokes equations journal January 2008
On physics-based preconditioning of the Navier–Stokes equations journal December 2009
Stabilization and scalable block preconditioning for the Navier–Stokes equations journal January 2012
A Hermite WENO reconstruction-based discontinuous Galerkin method for the Euler equations on tetrahedral grids journal June 2012
A reconstructed discontinuous Galerkin method based on a Hierarchical WENO reconstruction for compressible flows on tetrahedral grids journal March 2013
A Newton method with adaptive finite elements for solving phase-change problems with natural convection journal October 2014
Fully-implicit orthogonal reconstructed Discontinuous Galerkin method for fluid dynamics with phase change journal January 2016
A mathematical and numerical framework for the analysis of compressible thermal convection in gases at very high temperatures journal May 2016
A scalable, fully implicit algorithm for the reduced two-field low-β extended MHD model journal December 2016
Mesoscopic simulation model of selective laser melting of stainless steel powder journal November 2014
Modelling of Natural Convection Flows with Large Temperature Differences: A Benchmark Problem for Low Mach Number Solvers. Part 1. Reference Solutions journal May 2005
Solid velocity correction schemes for a temperature transforming model for convection phase change journal February 2006
Reducing Complexity in Parallel Algebraic Multigrid Preconditioners journal January 2006
Newton-GMRES Preconditioning for Discontinuous Galerkin Discretizations of the Navier–Stokes Equations journal January 2008
Estimation of Sparse Jacobian Matrices and Graph Coloring Blems journal February 1983
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems journal July 1986
Choosing the Forcing Terms in an Inexact Newton Method journal January 1996
Physics-Based Preconditioners for Ocean Simulation journal January 2013
A New Approximate Block Factorization Preconditioner for Two-Dimensional Incompressible (Reduced) Resistive MHD journal January 2013
A Multigrid-Preconditioned Newton--Krylov Method for the Incompressible Navier--Stokes Equations journal January 2001
Preconditioning Strategies for Fully Implicit Radiation Diffusion with Material-Energy Transfer journal January 2001
On Preconditioning Newton--Krylov Methods in Solidifying Flow Applications journal January 2001
Robustness and Scalability of Algebraic Multigrid journal January 2000
Preconditioning Techniques in Computational Fluid Dynamics journal January 1999
Truchas – a multi-physics tool for casting simulation journal August 2009
Preconditioning applied to variable and constant density flows journal November 1995
Newton-Krylov methods for low-Mach-number compressible combustion journal May 1996

Cited By (1)

High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics journal November 2018

Figures / Tables (27)

Figure 1(p. 8)figure Figure 1
Table 1(p. 9)table Table 1
Table 2(p. 11)table Table 2
Table 3(p. 11)table Table 3
Figure 2(p. 18)figure Figure 2
Figure 3(p. 19)figure Figure 3
Table 4(p. 19)table Table 4
Figure 4(p. 20)figure Figure 4
Figure 5(p. 21)figure Figure 5
Figure 6(p. 22)figure Figure 6
Figure 7(p. 22)figure Figure 7
Figure 8(p. 23)figure Figure 8
Figure 9(p. 23)figure Figure 9
Figure 10(p. 24)figure Figure 10
Table 5(p. 24)table Table 5
Figure 11(p. 25)figure Figure 11
Table 6(p. 25)table Table 6
Figure 12(p. 26)figure Figure 12
Figure 13(p. 27)figure Figure 13
Figure 14(p. 28)figure Figure 14
Table 7(p. 28)table Table 7
Figure 15(p. 29)figure Figure 15
Figure 16(p. 29)figure Figure 16
Figure 17(p. 30)figure Figure 17
Table 8(p. 30)table Table 8
Figure 18(p. 31)figure Figure 18
Table 9(p. 31)table Table 9

Similar Records

On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing
Thesis/Dissertation · Wed May 17 00:00:00 EDT 2017 · OSTI ID:1361587
A Scalable Preconditioner for a Primal Discontinuous Petrov--Galerkin Method
Journal Article · Mon Apr 23 20:00:00 EDT 2018 · SIAM Journal on Scientific Computing · OSTI ID:1734606
A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints
Journal Article · Sat Nov 29 19:00:00 EST 2025 · Numerical Linear Algebra with Applications · OSTI ID:3014007

Related Subjects

97 MATHEMATICS AND COMPUTING
All speed fluid dynamics
Block preconditioning
Fully implicit
Newton Krylov
Physics-based preconditioning
Reconstructed discontinuous Galerkin method