Development of a 2-D algorithm to simulate convection and phase transition efficiently

Katherine J. Evans; D. A. Knoll; Michael Pernice

doi:10.1016/j.jcp.2006.03.025

Development of a 2-D algorithm to simulate convection and phase transition efficiently

Katherine J. Evans, D. A. Knoll, Michael Pernice

Research output: Contribution to journal › Article › peer-review

22 Scopus citations

Abstract

We develop a Jacobian-Free Newton-Krylov (JFNK) method for the solution of a two-dimensional convection phase change model using the incompressible Navier-Stokes equation set and enthalpy as the energy conservation variable. The SIMPLE algorithm acts as a physics-based preconditioner to JFNK. This combined algorithm is compared to solutions using SIMPLE as the main solver. Algorithm performance is assessed for two benchmark problems of phase change convection of a pure material, one melting and one freezing. The JFNK-SIMPLE method is shown to be more efficient per time step and more robust at larger time steps. Overall CPU savings of more than an order of magnitude are realized.

Original language	English
Pages (from-to)	404-417
Number of pages	14
Journal	Journal of Computational Physics
Volume	219
Issue number	1
DOIs	https://doi.org/10.1016/j.jcp.2006.03.025
State	Published - Nov 20 2006
Externally published	Yes

Funding

This work was carried out under the auspices of the National Nuclear Security Administration of the US Department of Energy at Los Alamos National Laboratory under Contract No. W-7405-ENG-36 (LA-UR-06-0051).

Funders	Funder number
U.S. Department of Energy
National Nuclear Security Administration
Los Alamos National Laboratory	LA-UR-06-0051, W-7405-ENG-36

Keywords

Newton-Krylov
Phase change convection
SIMPLE preconditioner

Access to Document

10.1016/j.jcp.2006.03.025

Cite this

@article{e9beeb5848f741b8acc090a965305786,

title = "Development of a 2-D algorithm to simulate convection and phase transition efficiently",

abstract = "We develop a Jacobian-Free Newton-Krylov (JFNK) method for the solution of a two-dimensional convection phase change model using the incompressible Navier-Stokes equation set and enthalpy as the energy conservation variable. The SIMPLE algorithm acts as a physics-based preconditioner to JFNK. This combined algorithm is compared to solutions using SIMPLE as the main solver. Algorithm performance is assessed for two benchmark problems of phase change convection of a pure material, one melting and one freezing. The JFNK-SIMPLE method is shown to be more efficient per time step and more robust at larger time steps. Overall CPU savings of more than an order of magnitude are realized.",

keywords = "Newton-Krylov, Phase change convection, SIMPLE preconditioner",

author = "Evans, {Katherine J.} and Knoll, {D. A.} and Michael Pernice",

year = "2006",

month = nov,

day = "20",

doi = "10.1016/j.jcp.2006.03.025",

language = "English",

volume = "219",

pages = "404--417",

journal = "Journal of Computational Physics",

issn = "0021-9991",

publisher = "Elsevier",

number = "1",

}

TY - JOUR

T1 - Development of a 2-D algorithm to simulate convection and phase transition efficiently

AU - Evans, Katherine J.

AU - Knoll, D. A.

AU - Pernice, Michael

PY - 2006/11/20

Y1 - 2006/11/20

N2 - We develop a Jacobian-Free Newton-Krylov (JFNK) method for the solution of a two-dimensional convection phase change model using the incompressible Navier-Stokes equation set and enthalpy as the energy conservation variable. The SIMPLE algorithm acts as a physics-based preconditioner to JFNK. This combined algorithm is compared to solutions using SIMPLE as the main solver. Algorithm performance is assessed for two benchmark problems of phase change convection of a pure material, one melting and one freezing. The JFNK-SIMPLE method is shown to be more efficient per time step and more robust at larger time steps. Overall CPU savings of more than an order of magnitude are realized.

AB - We develop a Jacobian-Free Newton-Krylov (JFNK) method for the solution of a two-dimensional convection phase change model using the incompressible Navier-Stokes equation set and enthalpy as the energy conservation variable. The SIMPLE algorithm acts as a physics-based preconditioner to JFNK. This combined algorithm is compared to solutions using SIMPLE as the main solver. Algorithm performance is assessed for two benchmark problems of phase change convection of a pure material, one melting and one freezing. The JFNK-SIMPLE method is shown to be more efficient per time step and more robust at larger time steps. Overall CPU savings of more than an order of magnitude are realized.

KW - Newton-Krylov

KW - Phase change convection

KW - SIMPLE preconditioner

UR - http://www.scopus.com/inward/record.url?scp=33750304690&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2006.03.025

DO - 10.1016/j.jcp.2006.03.025

M3 - Article

AN - SCOPUS:33750304690

SN - 0021-9991

VL - 219

SP - 404

EP - 417

JO - Journal of Computational Physics

JF - Journal of Computational Physics

IS - 1

ER -

Development of a 2-D algorithm to simulate convection and phase transition efficiently

Abstract

Funding

Keywords

Access to Document

Other files and links

Fingerprint

Cite this