优化了五对角紧致差分格式.通过Fourier分析,将优化目标转化为一个求带有约束的多元非线性函数的最小值问题,利用序列二次规划(SQP)方法获得最佳系数.通过三种措施保证优化的格式具有高精度和分辨率:①直接对量化波数误差积分求其最小值;②采用绝对误差准则,使各种波长的波具有同一误差限;③优化的波数空间和精确求解区间一致.通过调整内部和边界格式的Taylor精度及误差限,保证整个格式的稳定性,并从理论上证明了优化格式具有渐进稳定性.一维和二维基准算例体现了优化格式的性能改进.
The optimized pentadiagonal compact finite difference schemes were presented in this paper. Through Fourier analysis, the optimization object was converted to find the minimum of a nonlinear multi-variable function with multi-constraint. The advanced sequential quadratic programming (SQP) method was employed to find the minimum. In order to obtain high accuracy and resolution, the following strategies were applied: (1) the minimum is computed directly from integral errors of the scaled wavenumbers; (2) absolute error criterion is used for evaluating the deviation between effective and exact wavenumbers, and the waves with different wavelengths have the same error limit under this criterion; (3) the wavenumber domain for optimization is identical with the well-resolved domain. The asymptotical stability of the schemes is satisfied by adjusting the Taylor accuracy and error limits, and proved through the theoretical analysis. The increased performances of the optimized schemes are demonstrated through application to benchmark problems.