An accuracy progressive sixth-order finite difference scheme
Chu, Peter C.
MetadataShow full item record
How to reduce the computational error is a key issue in numerical modeling and simulation. The higher the order of the difference scheme, the less the truncation error and the more complicated the computation. For compromise, a simple, three-point accuracy progressive (AP) finite-difference scheme for numerical calculation is proposed. The major features of the AP scheme are three-point, high-order accuracy, and accuracy progressive. The lower-order scheme acts as a ‘‘source’’ term in the higher-order scheme. This treatment keeps three-point schemes with high accuracy. The analytical error estimation shows the sixth-order accuracy that the AP scheme can reach. The Fourier analysis of errors indicates the accuracy improvement from lower-order to higher-order AP schemes. The Princeton Ocean Model (POM) implemented for the Japan/East Sea (JES) is used to evaluate the AP scheme. Consider a horizontally homogeneous and stably stratified JES with realistic topography.Without any forcing, initially motionless ocean will keep motionless forever; that is to say, there is a known solution (V 5 0). Any nonzero model velocity can be treated as an error. The stability and accuracy are compared with those of the second-order scheme in a series of calculations of unforced flow in the JES. The three-point sixthorder AP scheme is shown to have error reductions by factors of 10–20 compared to the second-order difference scheme. Due to their three-point grid structure, the AP schemes can be easily applied to current ocean and atmospheric models.
Journal of Atmospheric and Oceanic Technology, American Meteorological Society, 18, 1245-1257.
RightsThis publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.
Showing items related by title, author, creator and subject.
Fan, C.W.; Chu, Peter C. (1997);How to reduce the horizontal pressure gradient error is a key issue of using s-coordinate ocean models, especially of using primitive equation models for coastal regions. The error is caused by the splitting of the horizontal ...
Combining spectral and spatial information into hidden Markov models for unsupervised image classifciation Tso, B.; Olsen, R.C. (2005-05);Unsupervised classification methodology applied to remote sensing image processing can provide benefits in automatically converting the raw image data into useful information so long as higher classification accuracy is ...
Fan, C.W.; Chu, Peter C. (2010);With a finite volume approach, a flux-form semi-Lagrangian (TFSL) scheme with space-time transformation was developed to provide stable and accurate algorithm in solving the advection-diffusion equation. Different from the ...