Search Menu
DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information
Search Scholarly Publications

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

Abstract

We introduce a novel Newton–Krylov (NK)-based fully implicit algorithm for solving fluid flows in a wide range of flow conditions—from variable density nearly incompressible to supersonic shock dynamics. The key enabling feature of our all-speed solver is the ability to efficiently solve conservation laws by choosing a set of independent variables that produce a well-conditioned Jacobian matrix for the linear iterations of the global nonlinear iterative solver. In particular, instead of choosing to discretize the conservative variables (density, momentum, total energy), which is traditionally used in Eulerian high-speed compressible fluid dynamics, we demonstrate superior performance by discretizing the primitive variables—pressure–velocity–temperature in the very low-Mach flow limits or density–velocity–temperature/entropy in the shock dynamics range. Moreover, our method allows us to avoid direct inversion of the mass matrix in discrete time derivatives, which is usually an additional source for stiffness, especially pronounced when going to very high-order schemes with non-orthogonal basis functions. Here, we show robust solutions obtained for discontinuous finite element discretization up to seventh-order accuracy. Another important aspect of the solution algorithm is the Advection Upstream Splitting Method (AUSM), adopted to compute numerical fluxes within our reconstructed discontinuous Galerkin (rDG) spatial discretization scheme. The use of the low-Mach modification ofmore » the hyperbolic flux operator is found to be necessary for enabling robust simulations of very stiff liquids and metals for Mach numbers below M=10-5, which is well known to be very computationally challenging for compressible solvers. We demonstrate that our fully implicit rDG-NK solver with the AUSM+ -up flux treatment produces efficient and high-resolution numerical solutions at all speeds, ranging from vanishing Mach numbers to transonic and supersonic, without substantial modifications of the solution procedures. (At high speed, we add limiting and use a simpler preconditioning of the Krylov solver.) Numerical examples include nearly incompressible constant-property flow past a backward-facing step with heat transfer, low-Mach variable-property channel flow of water at supercritical state, phase change and melt pool dynamics for laser spot welding and selective laser melting in additive manufacturing, and Mach 3 flow in a wind tunnel with a step.« less

Authors:
ORCiD logo [1]; ORCiD logo [1];  [1];  [1];  [1];  [1]; ORCiD logo [2]
+ Show Author Affiliations
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States). Design Physics Dept.
  2. Univ. of California, Davis, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1562804
Report Number(s):
LLNL-JRNL-745590
Journal ID: ISSN 0938-1287; 900894; TRN: US2000762
Grant/Contract Number:  
AC52-07NA27344
Resource Type:
Accepted Manuscript
Journal Name:
Shock Waves
Additional Journal Information:
Journal Volume: 29; Journal Issue: 5; Journal ID: ISSN 0938-1287
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; reconstructed/recovery discontinuous Galerkin method; all-speed AUSM+-up scheme; fully implicit; Jacobian-free Newton-Krylov; preconditioning

Citation Formats

Nourgaliev, R., Greene, P., Weston, B., Barney, R., Anderson, A., Khairallah, S., and Delplanque, J. -P. High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics. United States: N. p., 2018. Web. doi:10.1007/s00193-018-0871-8.
Copy to clipboard
Nourgaliev, R., Greene, P., Weston, B., Barney, R., Anderson, A., Khairallah, S., & Delplanque, J. -P. High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics. United States. https://doi.org/10.1007/s00193-018-0871-8
Copy to clipboard
Nourgaliev, R., Greene, P., Weston, B., Barney, R., Anderson, A., Khairallah, S., and Delplanque, J. -P. Wed . "High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics". United States. https://doi.org/10.1007/s00193-018-0871-8. https://www.osti.gov/servlets/purl/1562804.
Copy to clipboard
@article{osti_1562804,
title = {High-order fully implicit solver for all-speed fluid dynamics: AUSM ride from nearly incompressible variable-density flows to shock dynamics},
author = {Nourgaliev, R. and Greene, P. and Weston, B. and Barney, R. and Anderson, A. and Khairallah, S. and Delplanque, J. -P.},
abstractNote = {We introduce a novel Newton–Krylov (NK)-based fully implicit algorithm for solving fluid flows in a wide range of flow conditions—from variable density nearly incompressible to supersonic shock dynamics. The key enabling feature of our all-speed solver is the ability to efficiently solve conservation laws by choosing a set of independent variables that produce a well-conditioned Jacobian matrix for the linear iterations of the global nonlinear iterative solver. In particular, instead of choosing to discretize the conservative variables (density, momentum, total energy), which is traditionally used in Eulerian high-speed compressible fluid dynamics, we demonstrate superior performance by discretizing the primitive variables—pressure–velocity–temperature in the very low-Mach flow limits or density–velocity–temperature/entropy in the shock dynamics range. Moreover, our method allows us to avoid direct inversion of the mass matrix in discrete time derivatives, which is usually an additional source for stiffness, especially pronounced when going to very high-order schemes with non-orthogonal basis functions. Here, we show robust solutions obtained for discontinuous finite element discretization up to seventh-order accuracy. Another important aspect of the solution algorithm is the Advection Upstream Splitting Method (AUSM), adopted to compute numerical fluxes within our reconstructed discontinuous Galerkin (rDG) spatial discretization scheme. The use of the low-Mach modification of the hyperbolic flux operator is found to be necessary for enabling robust simulations of very stiff liquids and metals for Mach numbers below M=10-5, which is well known to be very computationally challenging for compressible solvers. We demonstrate that our fully implicit rDG-NK solver with the AUSM+ -up flux treatment produces efficient and high-resolution numerical solutions at all speeds, ranging from vanishing Mach numbers to transonic and supersonic, without substantial modifications of the solution procedures. (At high speed, we add limiting and use a simpler preconditioning of the Krylov solver.) Numerical examples include nearly incompressible constant-property flow past a backward-facing step with heat transfer, low-Mach variable-property channel flow of water at supercritical state, phase change and melt pool dynamics for laser spot welding and selective laser melting in additive manufacturing, and Mach 3 flow in a wind tunnel with a step.},
doi = {10.1007/s00193-018-0871-8},
journal = {Shock Waves},
number = 5,
volume = 29,
place = {United States},
year = {Wed Nov 28 00:00:00 EST 2018},
month = {Wed Nov 28 00:00:00 EST 2018}
}
Copy to clipboard
Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.1007/s00193-018-0871-8
Copyright Statement
Other availability
Search WorldCat Search WorldCat to find libraries that may hold this journal
Citation Metrics:
Cited by: 2 works
Citation information provided by
Web of Science
Save / Share:
Export Metadata
Save to My Library
You must Sign In or Create an Account in order to save documents to your library.

Works referenced in this record:

A pseudocompressibility method for the numerical simulation of incompressible multifluid flows
journal, July 2004

Preconditioning applied to variable and constant density time-accurate flows on unstructured meshes
conference, February 2013

A Parallel Reconstructed Discontinuous Galerkin Method for Compressible Flows on Arbitrary Grids
conference, June 2012

  • Luo, Hong; Ali, Amjad; Nourgaliev, Robert
  • 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition
  • DOI: 10.2514/6.2010-366

Very high order PNPM schemes on unstructured meshes for the resistive relativistic MHD equations
journal, October 2009

A Reconstructed Discontinuous Galerkin Method for the Compressible Navier-Stokes Equations on Arbitrary Grids
conference, June 2012

  • Luo, Hong; Luo, Luqing; Nourgaliev, Robert
  • 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition
  • DOI: 10.2514/6.2010-364

A Parallel, Reconstructed Discontinuous Galerkin Method for the Compressible Flows on Arbitrary Grids
journal, February 2011

A Projection Method for Low Speed Flows
journal, March 1999

  • Colella, Phillip; Pao, Karen
  • Journal of Computational Physics, Vol. 149, Issue 2
  • DOI: 10.1006/jcph.1998.6152

Direct Numerical Simulation of Interfacial Flows: Implicit Sharp-Interface Method (I-SIM)
conference, June 2012

  • Nourgaliev, Robert; Theofanous, Theo; Park, HyeongKae
  • 46th AIAA Aerospace Sciences Meeting and Exhibit
  • DOI: 10.2514/6.2008-1453

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

Effect of stepheight on the separated flow past a backward facing step
journal, March 1989

  • Thangam, S.; Knight, Doyle D.
  • Physics of Fluids A: Fluid Dynamics, Vol. 1, Issue 3
  • DOI: 10.1063/1.857430

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

Surface tension of binary metal—surface active solute systems under conditions relevant to welding metallurgy
journal, June 1988

  • Sahoo, P.; Debroy, T.; McNallan, M. J.
  • Metallurgical Transactions B, Vol. 19, Issue 3
  • DOI: 10.1007/BF02657748

Review of preconditioning methods for fluid dynamics
journal, May 1993

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

Comparison of the AUSM+ and H-CUSP Schemes for Turbomachinery Applications
conference, June 2012

  • Chima, Rodrick; Liou, Meng-Sing
  • 16th AIAA Computational Fluid Dynamics Conference
  • DOI: 10.2514/6.2003-4120

A Numerical Scheme for the Compressible Low-Mach Number Regime of Ideal Fluid Dynamics
journal, January 2017

  • Barsukow, Wasilij; Edelmann, Philipp V. F.; Klingenberg, Christian
  • Journal of Scientific Computing, Vol. 72, Issue 2
  • DOI: 10.1007/s10915-017-0372-4

A Reconstructed Discontinuous Galerkin Method for the Euler Equations on Arbitrary Grids
journal, November 2012

  • Luo, Hong; Luo, Luqing; Nourgaliev, Robert
  • Communications in Computational Physics, Vol. 12, Issue 5
  • DOI: 10.4208/cicp.250911.030212a

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

Numerical solution of the Navier-Stokes equations
journal, January 1968

The numerical simulation of two-dimensional fluid flow with strong shocks
journal, April 1984

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

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

Discontinuous Galerkin for Diffusion
conference, June 2012

  • van Leer, Bram; Nomura, Shohei
  • 17th AIAA Computational Fluid Dynamics Conference
  • DOI: 10.2514/6.2005-5108

Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics
journal, November 2019

  • Weston, Brian; Nourgaliev, Robert; Delplanque, Jean-Pierre
  • Journal of Computational Physics, Vol. 397
  • DOI: 10.1016/j.jcp.2019.07.045

The influence of surfactants on thermocapillary flow instabilities in low Prandtl melting pools
journal, June 2016

  • Kidess, Anton; Kenjereš, Saša; Kleijn, Chris R.
  • Physics of Fluids, Vol. 28, Issue 6
  • DOI: 10.1063/1.4953797

A Reconstructed Discontinuous Galerkin Method for the Compressible Euler Equations on Arbitrary Grids
conference, June 2012

  • Luo, Hong; Luo, Luqing; Nourgaliev, Robert
  • 19th AIAA Computational Fluid Dynamics
  • DOI: 10.2514/6.2009-3788

An implicit Hermite WENO reconstruction-based discontinuous Galerkin method on tetrahedral grids
journal, June 2014

A low-Mach number fix for Roe’s approximate Riemann solver
journal, June 2011

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

Numerical solutions of pulsating flow and heat transfer characteristics in a channel with a backward-facing step
journal, February 1997

A numerical method for solving incompressible viscous flow problems
journal, August 1967

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

Modified multi-dimensional limiting process with enhanced shock stability on unstructured grids
journal, January 2018

A robust and accurate approach to computing compressible multiphase flow: Stratified flow model and AUSM+-up scheme
journal, July 2007

Computational study of laminar heat transfer downstream of a backward-facing step
journal, February 1993

The lattice Boltzmann equation method: theoretical interpretation, numerics and implications
journal, January 2003

  • Nourgaliev, R. R.; Dinh, T. N.; Theofanous, T. G.
  • International Journal of Multiphase Flow, Vol. 29, Issue 1
  • DOI: 10.1016/S0301-9322(02)00108-8

A reconstructed discontinuous Galerkin method for the compressible Navier–Stokes equations on arbitrary grids
journal, September 2010

  • Luo, Hong; Luo, Luqing; Nourgaliev, Robert
  • Journal of Computational Physics, Vol. 229, Issue 19
  • DOI: 10.1016/j.jcp.2010.05.033

Contribution of the recent AUSM schemes to the OVERFLOW code - Implementation and validation
conference, February 2013

  • Liou, Meng-Sing; Buning, Pieter
  • 18th Applied Aerodynamics Conference
  • DOI: 10.2514/6.2000-4404

A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes on unstructured meshes
journal, September 2008

  • Dumbser, Michael; Balsara, Dinshaw S.; Toro, Eleuterio F.
  • Journal of Computational Physics, Vol. 227, Issue 18
  • DOI: 10.1016/j.jcp.2008.05.025

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

Accurate Projection Methods for the Incompressible Navier–Stokes Equations
journal, April 2001

  • Brown, David L.; Cortez, Ricardo; Minion, Michael L.
  • Journal of Computational Physics, Vol. 168, Issue 2
  • DOI: 10.1006/jcph.2001.6715

L 2 Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers : L
journal, September 2015

  • Oßwald, K.; Siegmund, A.; Birken, P.
  • International Journal for Numerical Methods in Fluids, Vol. 81, Issue 2
  • DOI: 10.1002/fld.4175

Low-Diffusion Flux-Splitting Methods for Flows at All Speeds
journal, September 1998

  • Edwards, Jack R.; Liou, Meng-Sing
  • AIAA Journal, Vol. 36, Issue 9
  • DOI: 10.2514/2.587

An algorithm for unsteady viscous flows at all speeds
journal, January 2000

A numerical fluid dynamics calculation method for all flow speeds
journal, October 1971

An evaluation of several differencing methods for inviscid fluid flow problems
journal, February 1968

Hierarchical slope limiting in explicit and implicit discontinuous Galerkin methods
journal, January 2014

Characteristic time-stepping or local preconditioning of the Euler equations
conference, February 2013

  • Va, Bram; Lee, Wen-Tzong; Roe, Philip
  • 10th Computational Fluid Dynamics Conference
  • DOI: 10.2514/6.1991-1552

Numerical speed of sound and its application to schemes for all speeds
conference, February 2013

  • Liou, Meng-Sing; Edwards, Jack
  • 14th Computational Fluid Dynamics Conference
  • DOI: 10.2514/6.1999-3268

The design and application of upwind schemes on unstructured meshes
conference, February 2013

  • Barth, Timothy; Jespersen, Dennis
  • 27th Aerospace Sciences Meeting
  • DOI: 10.2514/6.1989-366

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

Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries
book, January 1997

  • Balay, Satish; Gropp, William D.; McInnes, Lois Curfman
  • Modern Software Tools for Scientific Computing
  • DOI: 10.1007/978-1-4612-1986-6_8

The Application of Preconditioning in Viscous Flows
journal, April 1993

Marangoni driven turbulence in high energy surface melting processes
journal, June 2016

Experimental and theoretical investigation of backward-facing step flow
journal, February 1983

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

On the accuracy of limiters and convergence to steady state solutions
conference, February 2013

ALE3D: An Arbitrary Lagrangian-Eulerian Multi-Physics Code
text, January 2017

Riemann Solvers and Numerical Methods for Fluid Dynamics
book, January 2009

A Hermit WENO Reconstruction-based Discontinuous Galerkin Method for the Euler Equations on Tetrahedral Grids
conference, November 2012

  • Luo, Hong; Li, Shujie; Xia, Yidong
  • 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition
  • DOI: 10.2514/6.2012-461

Low-diffusion flux-splitting methods for flows at all speeds
journal, January 1998

  • Edwards, Jack R.; Liou, Meng-Sing
  • AIAA Journal, Vol. 36
  • DOI: 10.2514/3.14013

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

Marangoni driven turbulence in high energy surface melting processes
text, January 2016

The influence of surfactants on thermocapillary flow instabilities in low Prandtl melting pools
text, January 2016

    Sort by:
    [ × clear filter / sort ]

    Similar Records in DOE PAGES and OSTI.GOV collections: