DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes

Abstract

In this study, an optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.

Authors:
 [1];  [2];  [3]
  1. Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States); The Univ. of Chicago, Chicago, IL (United States)
  2. The Univ. of Chicago, Chicago, IL (United States)
  3. The Univ. of Chicago, Chicago, IL (United States); Fermi National Accelerator Lab. (FNAL), Batavia, IL (United States)
Publication Date:
Research Org.:
Fermi National Accelerator Laboratory (FNAL), Batavia, IL (United States)
Sponsoring Org.:
USDOE Office of Science (SC), High Energy Physics (HEP)
OSTI Identifier:
1420907
Alternate Identifier(s):
OSTI ID: 1548848
Report Number(s):
arXiv:1801.03108; FERMILAB-PUB-18-025-A
Journal ID: ISSN 0021-9991; 1647346; TRN: US1801505
Grant/Contract Number:  
AC02-07CH11359
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 359; Journal Issue: C; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS; Numerical methods; Computational fluid dynamics; Partial differential equations

Citation Formats

Gnedin, Nickolay Y., Semenov, Vadim A., and Kravtsov, Andrey V. Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes. United States: N. p., 2018. Web. doi:10.1016/j.jcp.2018.01.008.
Gnedin, Nickolay Y., Semenov, Vadim A., & Kravtsov, Andrey V. Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes. United States. https://doi.org/10.1016/j.jcp.2018.01.008
Gnedin, Nickolay Y., Semenov, Vadim A., and Kravtsov, Andrey V. Tue . "Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes". United States. https://doi.org/10.1016/j.jcp.2018.01.008. https://www.osti.gov/servlets/purl/1420907.
@article{osti_1420907,
title = {Enforcing the Courant–Friedrichs–Lewy condition in explicitly conservative local time stepping schemes},
author = {Gnedin, Nickolay Y. and Semenov, Vadim A. and Kravtsov, Andrey V.},
abstractNote = {In this study, an optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally constrained time steps. However, such a scheme can result in violation of the Courant-Friedrichs-Lewy (CFL) condition, which is manifestly non-local. Although the violations can be considered to be "weak" in a certain sense and the corresponding numerical solution may be stable, such calculation does not guarantee the correct propagation speed for arbitrary waves. We use an experimental fluid dynamics code that allows cubic "patches" of grid cells to step with independent, locally constrained time steps to demonstrate how the CFL condition can be enforced by imposing a condition on the time steps of neighboring patches. We perform several numerical tests that illustrate errors introduced in the numerical solutions by weak CFL condition violations and show how strict enforcement of the CFL condition eliminates these errors. In all our tests the strict enforcement of the CFL condition does not impose a significant performance penalty.},
doi = {10.1016/j.jcp.2018.01.008},
journal = {Journal of Computational Physics},
number = C,
volume = 359,
place = {United States},
year = {Tue Jan 30 00:00:00 EST 2018},
month = {Tue Jan 30 00:00:00 EST 2018}
}

Journal Article:

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

Save / Share:

Works referenced in this record:

�ber die partiellen Differenzengleichungen der mathematischen Physik
journal, December 1928

  • Courant, R.; Friedrichs, K.; Lewy, H.
  • Mathematische Annalen, Vol. 100, Issue 1
  • DOI: 10.1007/BF01448839

Adaptive mesh refinement for hyperbolic partial differential equations
journal, March 1984


Local adaptive mesh refinement for shock hydrodynamics
journal, May 1989


Numerical approximations to nonlinear conservation laws with locally varying time and space grids
journal, January 1983


Time-accurate local time stepping method based on flux updating
journal, September 1994

  • Zhang, X. D.; Trepanier, J. -Y.; Reggio, M.
  • AIAA Journal, Vol. 32, Issue 9
  • DOI: 10.2514/3.12195

High resolution upwind-mixed finite element methods for advection-diffusion equations with variable time-stepping
journal, September 1995

  • Dawson, Clint
  • Numerical Methods for Partial Differential Equations, Vol. 11, Issue 5
  • DOI: 10.1002/num.1690110508

Interpolation finite difference schemes on grids locally refined in time
journal, November 2000

  • Shishkin, G. I.; Vabishchevich, P. N.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 190, Issue 8-10
  • DOI: 10.1016/S0045-7825(99)00451-X

Discontinuous Galerkin Time-Domain Methods for Multiscale Electromagnetic Simulations: A Review
journal, February 2013


A Discontinuous Galerkin Scheme Based on a Space–Time Expansion. I. Inviscid Compressible Flow in One Space Dimension
journal, February 2007


An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes - V. Local time stepping and p -adaptivity
journal, November 2007


An explicit discontinuous Galerkin scheme with local time-stepping for general unsteady diffusion equations
journal, May 2008

  • Lörcher, Frieder; Gassner, Gregor; Munz, Claus-Dieter
  • Journal of Computational Physics, Vol. 227, Issue 11
  • DOI: 10.1016/j.jcp.2008.02.015

Arbitrary-Lagrangian–Eulerian ADER–WENO finite volume schemes with time-accurate local time stepping for hyperbolic conservation laws
journal, October 2014


A conservative finite volume scheme with time-accurate local time stepping for scalar transport on unstructured grids
journal, December 2015


Multirate Timestepping Methods for Hyperbolic Conservation Laws
journal, September 2007


Multirate Explicit Adams Methods for Time Integration of Conservation Laws
journal, August 2008


Cosmological hydrodynamics with adaptive mesh refinement: A new high resolution code called RAMSES
journal, April 2002


Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence
journal, July 2011


Arbitrary-Lagrangian–Eulerian ADER–WENO finite volume schemes with time-accurate local time stepping for hyperbolic conservation laws
journal, October 2014


Works referencing / citing this record:

Application of a Projection Method for Simulating Flow of a Shear-Thinning Fluid
journal, July 2019