Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning

doi:10.1137/15M1044369

Title: Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning

Article Details
Other Related Research

Here, this paper presents a characteristic-based flux partitioning for the semi-implicit time integration of atmospheric flows. Nonhydrostatic models require the solution of the compressible Euler equations. The acoustic time scale is significantly faster than the advective scale, yet it is typically not relevant to atmospheric and weather phenomena. The acoustic and advective components of the hyperbolic flux are separated in the characteristic space. High-order, conservative additive Runge-Kutta methods are applied to the partitioned equations so that the acoustic component is integrated in time implicitly with an unconditionally stable method, while the advective component is integrated explicitly. The time step of the overall algorithm is thus determined by the advective scale. Benchmark flow problems are used to demonstrate the accuracy, stability, and convergence of the proposed algorithm. The computational cost of the partitioned semi-implicit approach is compared with that of explicit time integration.

Authors:

Ghosh, Debojyoti ^[1] ; Constantinescu, Emil M. ^[1]

Argonne National Lab. (ANL), Lemont, IL (United States)

Publication Date:: 2016-06-23

Grant/Contract Number:: AC02-06CH11357

Type:: Accepted Manuscript

Journal Name:: SIAM Journal on Scientific Computing

Additional Journal Information:: Journal Volume: 38; Journal Issue: 3; Journal ID: ISSN 1064-8275

Publisher:: SIAM

Research Org:: Argonne National Lab. (ANL), Argonne, IL (United States)

Sponsoring Org:: USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21)

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; 54 ENVIRONMENTAL SCIENCES; Euler equations; atmospheric flows; characteristic-based splitting; compressible; implicit-explicit time integration; nonhydrostatic

OSTI Identifier:: 1362122


                    Ghosh, Debojyoti, and Constantinescu, Emil M. Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning.  United States: N. p.,  
        Web.  doi:10.1137/15M1044369.

Copy to clipboard


                    Ghosh, Debojyoti, & Constantinescu, Emil M. Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning.  United States.  doi:10.1137/15M1044369.

Copy to clipboard


                    Ghosh, Debojyoti, and Constantinescu, Emil M. 2016.  
        "Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning".  United States.  
        doi:10.1137/15M1044369.  https://www.osti.gov/servlets/purl/1362122.

Copy to clipboard


                    
@article{osti_1362122,

  title        = {Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning},

  author       = {Ghosh, Debojyoti and Constantinescu, Emil M.},

  abstractNote = {Here, this paper presents a characteristic-based flux partitioning for the semi-implicit time integration of atmospheric flows. Nonhydrostatic models require the solution of the compressible Euler equations. The acoustic time scale is significantly faster than the advective scale, yet it is typically not relevant to atmospheric and weather phenomena. The acoustic and advective components of the hyperbolic flux are separated in the characteristic space. High-order, conservative additive Runge-Kutta methods are applied to the partitioned equations so that the acoustic component is integrated in time implicitly with an unconditionally stable method, while the advective component is integrated explicitly. The time step of the overall algorithm is thus determined by the advective scale. Benchmark flow problems are used to demonstrate the accuracy, stability, and convergence of the proposed algorithm. The computational cost of the partitioned semi-implicit approach is compared with that of explicit time integration.},

  doi          = {10.1137/15M1044369},

  journal      = {SIAM Journal on Scientific Computing},

  number       = 3,

  volume       = 38,

  place        = {United States},

  year         = {2016},

  month        = {6}

}

Copy to clipboard

Free Publicly Accessible Full Text
Accepted Manuscript (DOE)
Publisher's Version of Record
10.1137/15M1044369
https://doi.org/10.1137/15M1044369

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

Save / Share this Record
Export Metadata
- Endnote
- RIS
- Excel
- CSV
- XML
Save to My Library
Send to Email

Send to Email

Email address:

Content:

Similar Records

Similar records in DOE PAGES and OSTI.GOV collections:

Kinetic Simulation of Collisional Magnetized Plasmas with Semi-implicit Time Integration

Journal Article Ghosh, Debojyoti ; Dorf, Mikhail A. ; Dorr, Milo R. ; ... May 2018 - Journal of Scientific Computing

Plasmas with varying collisionalities occur in many applications, such as tokamak edge regions, where the flows are characterized by significant variations in density and temperature. While a kinetic model is necessary for weakly-collisional high-temperature plasmas, high collisionality in colder regions render the equations numerically stiff due to disparate time scales. In this study, we propose an implicit–explicit algorithm for such cases, where the collisional term is integrated implicitly in time, while the advective term is integrated explicitly in time, thus allowing time step sizes that are comparable to the advective time scales. This partitioning results in a more efficient algorithmmore »« less
A computationally-efficient, semi-implicit, iterative method for the time-integration of reacting flows with stiff chemistry

Journal Article Savard, B. ; Xuan, Y. ; Bobbitt, B. ; ... August 2015 - Journal of Computational Physics

Cited by 9Full Text Available
A comparative study of Rosenbrock-type and implicit Runge-Kutta time integration for discontinuous Galerkin method for unsteady 3D compressible Navier-Stokes equations

Journal Article Liu, Xiaodong ; Xia, Yidong ; Luo, Hong ; ... October 2016 - Communications in Computational Physics

A comparative study of two classes of third-order implicit time integration schemes is presented for a third-order hierarchical WENO reconstructed discontinuous Galerkin (rDG) method to solve the 3D unsteady compressible Navier-Stokes equations: — 1) the explicit first stage, single diagonally implicit Runge-Kutta (ESDIRK3) scheme, and 2) the Rosenbrock-Wanner (ROW) schemes based on the differential algebraic equations (DAEs) of Index-2. Compared with the ESDIRK3 scheme, a remarkable feature of the ROW schemes is that, they only require one approximate Jacobian matrix calculation every time step, thus considerably reducing the overall computational cost. A variety of test cases, ranging from inviscid flowsmore »« less
Cited by 3Full Text Available
Trust-region based return mapping algorithm for implicit integration of elastic-plastic constitutive models

Journal Article Lester, Brian ; Scherzinger, William January 2017 - International Journal for Numerical Methods in Engineering

Here, a new method for the solution of the non-linear equations forming the core of constitutive model integration is proposed. Specifically, the trust-region method that has been developed in the numerical optimization community is successfully modified for use in implicit integration of elastic-plastic models. Although attention here is restricted to these rate-independent formulations, the proposed approach holds substantial promise for adoption with models incorporating complex physics, multiple inelastic mechanisms, and/or multiphysics. As a first step, the non-quadratic Hosford yield surface is used as a representative case to investigate computationally challenging constitutive models. The theory and implementation are presented, discussed, andmore »« less
Cited by 1Full Text Available
A third-order implicit discontinuous Galerkin method based on a Hermite WENO reconstruction for time-accurate solution of the compressible Navier-Stokes equations

Journal Article Xia, Yidong ; Liu, Xiaodong ; Luo, Hong ; ... June 2015 - International Journal for Numerical Methods in Fluids

Here, a space and time third-order discontinuous Galerkin method based on a Hermite weighted essentially non-oscillatory reconstruction is presented for the unsteady compressible Euler and Navier–Stokes equations. At each time step, a lower-upper symmetric Gauss–Seidel preconditioned generalized minimal residual solver is used to solve the systems of linear equations arising from an explicit first stage, single diagonal coefficient, diagonally implicit Runge–Kutta time integration scheme. The performance of the developed method is assessed through a variety of unsteady flow problems. Numerical results indicate that this method is able to deliver the designed third-order accuracy of convergence in both space and time,more »« less
Cited by 4Full Text Available

Access journal articles and accepted manuscripts resulting from DOE research funding

Title: Semi-implicit time integration of atmospheric flows with characteristic-based flux partitioning

Access journal articles and accepted manuscripts
resulting from DOE research funding