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; TRN: US2102668
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. doi: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 = {Thu Jul 25 00:00:00 EDT 2019},
month = {Thu Jul 25 00:00:00 EDT 2019}
}

Journal Article:

Citation Metrics:
Cited by: 3 works
Citation information provided by
Web of Science

Figures / Tables:

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

Save / Share:

Works referenced in this record:

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

On the behaviour of upwind schemes in the low Mach number limit
journal, January 1999


A Hermite WENO reconstruction-based discontinuous Galerkin method for the Euler equations on tetrahedral grids
journal, June 2012


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

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


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

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


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

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


Preconditioning Techniques in Computational Fluid Dynamics
journal, January 1999


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

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

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

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


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

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

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 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

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

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

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

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


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

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

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

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

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


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


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


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

A reconstructed discontinuous Galerkin method based on a Hierarchical WENO reconstruction for compressible flows on tetrahedral grids
journal, March 2013


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

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


BoomerAMG: A parallel algebraic multigrid solver and preconditioner
journal, April 2002


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

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

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


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


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

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

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


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


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


The Application of Preconditioning in Viscous Flows
journal, April 1993


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


A Sequel to AUSM: AUSM+
journal, December 1996


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

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


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


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

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


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

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

Works referencing / citing this record: