skip to main content
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

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

Abstract

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 heatedmore » 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).« less

Authors:
 [1]; ORCiD logo [1];  [2];  [1]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Univ. of California, Davis, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1577941
Alternate Identifier(s):
OSTI ID: 1701812
Report Number(s):
LLNL-JRNL-745515
Journal ID: ISSN 0021-9991; 900481
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 397; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Block preconditioning; Physics-based preconditioning; Fully implicit; Newton Krylov; All speed fluid dynamics; Reconstructed discontinuous Galerkin method

Citation Formats

Weston, Brian, Nourgaliev, Robert, Delplanque, Jean -Pierre, and Barker, Andrew T. Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics. United States: N. p., 2019. Web. https://doi.org/10.1016/j.jcp.2019.07.045.
Weston, Brian, Nourgaliev, Robert, Delplanque, Jean -Pierre, & Barker, Andrew T. Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics. United States. https://doi.org/10.1016/j.jcp.2019.07.045
Weston, Brian, Nourgaliev, Robert, Delplanque, Jean -Pierre, and Barker, Andrew T. Thu . "Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics". United States. https://doi.org/10.1016/j.jcp.2019.07.045. https://www.osti.gov/servlets/purl/1577941.
@article{osti_1577941,
title = {Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics},
author = {Weston, Brian and Nourgaliev, Robert and Delplanque, Jean -Pierre and Barker, Andrew T.},
abstractNote = {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).},
doi = {10.1016/j.jcp.2019.07.045},
journal = {Journal of Computational Physics},
number = C,
volume = 397,
place = {United States},
year = {2019},
month = {7}
}

Journal Article:

Figures / Tables:

Figure 1 Figure 1: Thermal model, u(T), for equilibrium phase change.

Save / Share:

Works referenced in this record:

Numerical simulations of a transient injection flow at low Mach number regime
journal, October 2008

  • Beccantini, A.; Studer, E.; Gounand, S.
  • International Journal for Numerical Methods in Engineering, Vol. 76, Issue 5
  • DOI: 10.1002/nme.2331

Preconditioning Techniques for Large Linear Systems: A Survey
journal, November 2002


High-resolution viscous flow simulations at arbitrary Mach number
journal, January 2003


Preconditioning Strategies for Fully Implicit Radiation Diffusion with Material-Energy Transfer
journal, January 2001


A scalable, fully implicit algorithm for the reduced two-field low-β extended MHD model
journal, December 2016


The Application of Preconditioning in Viscous Flows
journal, April 1993


A Numerical Method for Solving Incompressible Viscous Flow Problems
journal, August 1997


Robustness and Scalability of Algebraic Multigrid
journal, January 2000

  • Cleary, Andrew J.; Falgout, Robert D.; Henson, Van Emden
  • SIAM Journal on Scientific Computing, Vol. 21, Issue 5
  • DOI: 10.1137/S1064827598339402

Estimation of Sparse Jacobian Matrices and Graph Coloring Blems
journal, February 1983

  • Coleman, Thomas F.; Moré, Jorge J.
  • SIAM Journal on Numerical Analysis, Vol. 20, Issue 1
  • DOI: 10.1137/0720013

Stabilization and scalable block preconditioning for the Navier–Stokes equations
journal, January 2012

  • Cyr, Eric C.; Shadid, John N.; Tuminaro, Raymond S.
  • Journal of Computational Physics, Vol. 231, Issue 2
  • DOI: 10.1016/j.jcp.2011.09.001

A New Approximate Block Factorization Preconditioner for Two-Dimensional Incompressible (Reduced) Resistive MHD
journal, January 2013

  • Cyr, Eric C.; Shadid, John N.; Tuminaro, Raymond S.
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 3
  • DOI: 10.1137/12088879X

A Newton method with adaptive finite elements for solving phase-change problems with natural convection
journal, October 2014

  • Danaila, Ionut; Moglan, Raluca; Hecht, Frédéric
  • Journal of Computational Physics, Vol. 274
  • DOI: 10.1016/j.jcp.2014.06.036

Modelling liquid-solid phase changes with melt convection
journal, August 1989

  • Dantzig, Jonathan A.
  • International Journal for Numerical Methods in Engineering, Vol. 28, Issue 8
  • DOI: 10.1002/nme.1620280805

Reducing Complexity in Parallel Algebraic Multigrid Preconditioners
journal, January 2006

  • De Sterck, Hans; Yang, Ulrike Meier; Heys, Jeffrey J.
  • SIAM Journal on Matrix Analysis and Applications, Vol. 27, Issue 4
  • DOI: 10.1137/040615729

Simulation of time-dependent pool shape during laser spot welding: Transient effects
journal, December 2003

  • Ehlen, Georg; Ludwig, Andreas; Sahm, Peter R.
  • Metallurgical and Materials Transactions A, Vol. 34, Issue 12
  • DOI: 10.1007/s11661-003-0194-x

Choosing the Forcing Terms in an Inexact Newton Method
journal, January 1996

  • Eisenstat, Stanley C.; Walker, Homer F.
  • SIAM Journal on Scientific Computing, Vol. 17, Issue 1
  • DOI: 10.1137/0917003

A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier–Stokes equations
journal, January 2008

  • Elman, Howard; Howle, V. E.; Shadid, John
  • Journal of Computational Physics, Vol. 227, Issue 3
  • DOI: 10.1016/j.jcp.2007.09.026

A parallel block multi-level preconditioner for the 3D incompressible Navier–Stokes equations
journal, May 2003


Development of a 2-D algorithm to simulate convection and phase transition efficiently
journal, November 2006

  • Evans, Katherine J.; Knoll, D. A.; Pernice, Michael
  • Journal of Computational Physics, Vol. 219, Issue 1
  • DOI: 10.1016/j.jcp.2006.03.025

High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method
journal, December 1982


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


Mesoscopic simulation model of selective laser melting of stainless steel powder
journal, November 2014


Laser powder-bed fusion additive manufacturing: Physics of complex melt flow and formation mechanisms of pores, spatter, and denudation zones
journal, April 2016


Jacobian-free Newton–Krylov methods: a survey of approaches and applications
journal, January 2004


Newton-Krylov methods for low-Mach-number compressible combustion
journal, May 1996

  • Knoll, D. A.; McHugh, P. R.; Keyes, D. E.
  • AIAA Journal, Vol. 34, Issue 5
  • DOI: 10.2514/3.13174

On Preconditioning Newton--Krylov Methods in Solidifying Flow Applications
journal, January 2001

  • Knoll, D. A.; VanderHeyden, W. B.; Mousseau, V. A.
  • SIAM Journal on Scientific Computing, Vol. 23, Issue 2
  • DOI: 10.1137/S1064827500374303

Jacobian–Free Newton–Krylov Methods for the Accurate Time Integration of Stiff Wave Systems
journal, October 2005

  • Knoll, D. A.; Mousseau, V. A.; Chacón, L.
  • Journal of Scientific Computing, Vol. 25, Issue 1
  • DOI: 10.1007/s10915-004-4640-8

Truchas – a multi-physics tool for casting simulation
journal, August 2009


Modelling of Natural Convection Flows with Large Temperature Differences: A Benchmark Problem for Low Mach Number Solvers. Part 1. Reference Solutions
journal, May 2005

  • Le Quéré, Patrick; Weisman, Catherine; Paillère, Henri
  • ESAIM: Mathematical Modelling and Numerical Analysis, Vol. 39, Issue 3
  • DOI: 10.1051/m2an:2005027

Performance of fully coupled algebraic multilevel domain decomposition preconditioners for incompressible flow and transport
journal, January 2006

  • Lin, Paul T.; Sala, Marzio; Shadid, John N.
  • International Journal for Numerical Methods in Engineering, Vol. 67, Issue 2
  • DOI: 10.1002/nme.1624

A Sequel to AUSM: AUSM+
journal, December 1996


A sequel to AUSM, Part II: AUSM+-up for all speeds
journal, May 2006


A New Flux Splitting Scheme
journal, July 1993

  • Liou, Meng-Sing; Steffen, Christopher J.
  • Journal of Computational Physics, Vol. 107, Issue 1
  • DOI: 10.1006/jcph.1993.1122

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


Solid velocity correction schemes for a temperature transforming model for convection phase change
journal, February 2006

  • Ma, Zhanhua; Zhang, Yuwen
  • International Journal of Numerical Methods for Heat & Fluid Flow, Vol. 16, Issue 2
  • DOI: 10.1108/09615530610644271

A finite element method for low-speed compressible flows
journal, May 2004

  • Martinez, Mario J.; Gartling, David K.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 193, Issue 21-22
  • DOI: 10.1016/j.cma.2003.12.049

Physics-Based Preconditioning and the Newton–Krylov Method for Non-equilibrium Radiation Diffusion
journal, May 2000

  • Mousseau, V. A.; Knoll, D. A.; Rider, W. J.
  • Journal of Computational Physics, Vol. 160, Issue 2
  • DOI: 10.1006/jcph.2000.6488

The extension of incompressible flow solvers to the weakly compressible regime
journal, February 2003


Physics-Based Preconditioners for Ocean Simulation
journal, January 2013

  • Newman, C.; Knoll, D. A.
  • SIAM Journal on Scientific Computing, Vol. 35, Issue 5
  • DOI: 10.1137/120881397

Fully-implicit orthogonal reconstructed Discontinuous Galerkin method for fluid dynamics with phase change
journal, January 2016


On physics-based preconditioning of the Navier–Stokes equations
journal, December 2009

  • Park, HyeongKae; Nourgaliev, Robert R.; Martineau, Richard C.
  • Journal of Computational Physics, Vol. 228, Issue 24
  • DOI: 10.1016/j.jcp.2009.09.015

A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows
journal, October 1972


A Multigrid-Preconditioned Newton--Krylov Method for the Incompressible Navier--Stokes Equations
journal, January 2001


Newton-GMRES Preconditioning for Discontinuous Galerkin Discretizations of the Navier–Stokes Equations
journal, January 2008

  • Persson, P. -O.; Peraire, J.
  • SIAM Journal on Scientific Computing, Vol. 30, Issue 6
  • DOI: 10.1137/070692108

GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

Performance of fully coupled domain decomposition preconditioners for finite element transport/reaction simulations
journal, May 2005

  • Shadid, J. N.; Tuminaro, R. S.; Devine, K. D.
  • Journal of Computational Physics, Vol. 205, Issue 1
  • DOI: 10.1016/j.jcp.2004.10.038

On a multilevel preconditioning module for unstructured mesh Krylov solvers: two-level Schwarz
journal, March 2002

  • Tuminaro, R. S.; Tong, C. H.; Shadid, J. N.
  • Communications in Numerical Methods in Engineering, Vol. 18, Issue 6
  • DOI: 10.1002/cnm.478

Preconditioned methods for solving the incompressible and low speed compressible equations
journal, October 1987


Preconditioning Techniques in Computational Fluid Dynamics
journal, January 1999


A fixed grid numerical modelling methodology for convection-diffusion mushy region phase-change problems
journal, August 1987


Preconditioning applied to variable and constant density flows
journal, November 1995

  • Weiss, Jonathan M.; Smith, Wayne A.
  • AIAA Journal, Vol. 33, Issue 11
  • DOI: 10.2514/3.12946

    Works referencing / citing this record: