An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling
In this paper， a new numerical scheme was proposed to solve the advection equation in a multi-moment nonhydrostatic dynamical core.To guarantee the shape-preserving property， the limiting operations are devised for a hybrid discretization framework adopted by the multi-moment dynamical core， consist...
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | zho |
| Published: |
Science Press, PR China
2024-04-01
|
| Series: | Gaoyuan qixiang |
| Subjects: | |
| Online Access: | http://www.gyqx.ac.cn/EN/10.7522/j.issn.1000-0534.2023.00068 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1850200716071665664 |
|---|---|
| author | Zhaoyang SUN Chungang CHEN Xingliang LI Xueshun SHEN |
| author_facet | Zhaoyang SUN Chungang CHEN Xingliang LI Xueshun SHEN |
| author_sort | Zhaoyang SUN |
| collection | DOAJ |
| description | In this paper， a new numerical scheme was proposed to solve the advection equation in a multi-moment nonhydrostatic dynamical core.To guarantee the shape-preserving property， the limiting operations are devised for a hybrid discretization framework adopted by the multi-moment dynamical core， consisting of the multi-moment finite-volume and the conservative finite-difference schemes for the horizontal and vertical discretizations respectively.In the horizontal direction， a nonoscillaory scheme is accomplished by adjusting the slope of the multi-moment reconstruction polynomial at the cell center with the application of a WENO （weighted essentially non-oscillatory） algorithm.The resulting multi-moment scheme can achieve the fourth-order accuracy in the convergence test.In the vertical direction， a TVD （total variation diminishing） slope limiter is applied in the finite-difference discretization to remove the non-physical oscillations around the discontinuities.To accomplish the time marching in the proposed advection model， a second-order paired explicit Runge-Kutta scheme is adopted， which is expected to be an efficient and practical method for the advection solvers in the atmospheric models with very high spatial resolutions.The explicit time marching， without the dimension splitting， is useful to avoid the divergence errors in the advection transport calculations.Two Runge-Kutta schemes， requiring different times of conducting the spatial discretization within a time step， are combined， and used for the time marching in the different directions.The finite-difference discretization is called for six times within a time step in order to increase the maximum available CFL （Courant-Friedrichs-Lewy） number in the vertical direction， while the horizontal multi-moment spatial discretization is conducted for two times as the regular second-order schemes.As a result， the difference between the maximum time steps determined by the horizontal and vertical discretizations， due to the very large aspect ratio of the computational cells in atmospheric modeling， can be diminished.The non-negativity property of the proposed advection scheme is assured by devising a new flux-correction algorithm.It improves the existing positivity-preserving algorithm through further considering the mass flowing into the computational cell in an iterative procedure during the flux-correction operations.The proposed flux-correction algorithm can approach the necessary and sufficient condition for assuring the non-negative solutions and is more accurate for the advection calculations with CFL numbers larger than one.The widely used two-dimensional benchmark tests were checked in this study and the numerical results verified the performance the proposed advection scheme， which has the practical potential to build an accurate and efficient advection equation solver for the scalable high-resolution nonhydrostatic atmospheric models. |
| format | Article |
| id | doaj-art-492cbc85cfd346b7b3fc72c9e5d7a31a |
| institution | OA Journals |
| issn | 1000-0534 |
| language | zho |
| publishDate | 2024-04-01 |
| publisher | Science Press, PR China |
| record_format | Article |
| series | Gaoyuan qixiang |
| spelling | doaj-art-492cbc85cfd346b7b3fc72c9e5d7a31a2025-08-20T02:12:15ZzhoScience Press, PR ChinaGaoyuan qixiang1000-05342024-04-0143252052810.7522/j.issn.1000-0534.2023.000681000-0534（2024）02-0520-09An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric ModelingZhaoyang SUN0Chungang CHEN1Xingliang LI2Xueshun SHEN3Department of Mechanics， Xi’an Jiaotong University， Xi’an 710049， Shaanxi， ChinaDepartment of Mechanics， Xi’an Jiaotong University， Xi’an 710049， Shaanxi， ChinaCenter for Earth System Modelling and Predication， China Meteorological Administration， Beijing 100081， ChinaCenter for Earth System Modelling and Predication， China Meteorological Administration， Beijing 100081， ChinaIn this paper， a new numerical scheme was proposed to solve the advection equation in a multi-moment nonhydrostatic dynamical core.To guarantee the shape-preserving property， the limiting operations are devised for a hybrid discretization framework adopted by the multi-moment dynamical core， consisting of the multi-moment finite-volume and the conservative finite-difference schemes for the horizontal and vertical discretizations respectively.In the horizontal direction， a nonoscillaory scheme is accomplished by adjusting the slope of the multi-moment reconstruction polynomial at the cell center with the application of a WENO （weighted essentially non-oscillatory） algorithm.The resulting multi-moment scheme can achieve the fourth-order accuracy in the convergence test.In the vertical direction， a TVD （total variation diminishing） slope limiter is applied in the finite-difference discretization to remove the non-physical oscillations around the discontinuities.To accomplish the time marching in the proposed advection model， a second-order paired explicit Runge-Kutta scheme is adopted， which is expected to be an efficient and practical method for the advection solvers in the atmospheric models with very high spatial resolutions.The explicit time marching， without the dimension splitting， is useful to avoid the divergence errors in the advection transport calculations.Two Runge-Kutta schemes， requiring different times of conducting the spatial discretization within a time step， are combined， and used for the time marching in the different directions.The finite-difference discretization is called for six times within a time step in order to increase the maximum available CFL （Courant-Friedrichs-Lewy） number in the vertical direction， while the horizontal multi-moment spatial discretization is conducted for two times as the regular second-order schemes.As a result， the difference between the maximum time steps determined by the horizontal and vertical discretizations， due to the very large aspect ratio of the computational cells in atmospheric modeling， can be diminished.The non-negativity property of the proposed advection scheme is assured by devising a new flux-correction algorithm.It improves the existing positivity-preserving algorithm through further considering the mass flowing into the computational cell in an iterative procedure during the flux-correction operations.The proposed flux-correction algorithm can approach the necessary and sufficient condition for assuring the non-negative solutions and is more accurate for the advection calculations with CFL numbers larger than one.The widely used two-dimensional benchmark tests were checked in this study and the numerical results verified the performance the proposed advection scheme， which has the practical potential to build an accurate and efficient advection equation solver for the scalable high-resolution nonhydrostatic atmospheric models.http://www.gyqx.ac.cn/EN/10.7522/j.issn.1000-0534.2023.00068advection equationpaired explicit runge-kutta schemeslope limiternon-negative schememulti-moment method |
| spellingShingle | Zhaoyang SUN Chungang CHEN Xingliang LI Xueshun SHEN An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling Gaoyuan qixiang advection equation paired explicit runge-kutta scheme slope limiter non-negative scheme multi-moment method |
| title | An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling |
| title_full | An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling |
| title_fullStr | An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling |
| title_full_unstemmed | An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling |
| title_short | An Advection Scheme using Paired Explicit Runge-Kutta Time Integration for Atmospheric Modeling |
| title_sort | advection scheme using paired explicit runge kutta time integration for atmospheric modeling |
| topic | advection equation paired explicit runge-kutta scheme slope limiter non-negative scheme multi-moment method |
| url | http://www.gyqx.ac.cn/EN/10.7522/j.issn.1000-0534.2023.00068 |
| work_keys_str_mv | AT zhaoyangsun anadvectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT chungangchen anadvectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT xingliangli anadvectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT xueshunshen anadvectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT zhaoyangsun advectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT chungangchen advectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT xingliangli advectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling AT xueshunshen advectionschemeusingpairedexplicitrungekuttatimeintegrationforatmosphericmodeling |