All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 43 11 1
PDF Downloads 7 3 1

A Comparative Study of Numerical Advection Schemes Featuring a One-Step Modified WKL Algorithm

Ching-Yuang HuangDepartment of Marine, Earth, and Atmospheric Sciences, North Carolina State University, Raleigh, North Carolina

Search for other papers by Ching-Yuang Huang in
Current site
Google Scholar
PubMed
Close
and
Sethu RamanDepartment of Marine, Earth, and Atmospheric Sciences, North Carolina State University, Raleigh, North Carolina

Search for other papers by Sethu Raman in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

A fourth-order Crowley-type advection scheme based on the multistep Warming-Kutler-Lomax (WKL) scheme is proposed in this study. This scheme utilizes a free parameter to minimize dispersion and dissipation and can be used to represent the advection of positive-definite scalars (such as moisture).

Linear Fourier component analyses indicate that the fourth-order Crowley-type scheme can reproduce the features of other modified Crowley-type schemes of third order, such as the scheme of Schlesinger and the quadratic upstream interpolation. Using the free parameter, the scheme may illustrate the limitation of the Crowley-type schemes for which diffusion is required for numerical stability of advective quantity. For these schemes, formulations that preserve amplitude are inevitably associated with smaller time steps. Adding the first cross-space term into these schemes could eliminate marginal instability or overshooting in linear advection.

Linear and nonlinear advection tests show that the performance of the proposed scheme is comparable to the fourth-order leapfrog scheme (which requires more computer memory) and the cubic upstream spline (which requires more computer time). This two-time-level advection scheme can thus be used in a numerical model to save computer resources.

Abstract

A fourth-order Crowley-type advection scheme based on the multistep Warming-Kutler-Lomax (WKL) scheme is proposed in this study. This scheme utilizes a free parameter to minimize dispersion and dissipation and can be used to represent the advection of positive-definite scalars (such as moisture).

Linear Fourier component analyses indicate that the fourth-order Crowley-type scheme can reproduce the features of other modified Crowley-type schemes of third order, such as the scheme of Schlesinger and the quadratic upstream interpolation. Using the free parameter, the scheme may illustrate the limitation of the Crowley-type schemes for which diffusion is required for numerical stability of advective quantity. For these schemes, formulations that preserve amplitude are inevitably associated with smaller time steps. Adding the first cross-space term into these schemes could eliminate marginal instability or overshooting in linear advection.

Linear and nonlinear advection tests show that the performance of the proposed scheme is comparable to the fourth-order leapfrog scheme (which requires more computer memory) and the cubic upstream spline (which requires more computer time). This two-time-level advection scheme can thus be used in a numerical model to save computer resources.

Save