POD/DEIM nonlinear model order reduction of an ADI implicit shallow water equations model

被引:112
作者
Stefanescu, R. [1 ]
Navon, I. M. [1 ]
机构
[1] Florida State Univ, Dept Comp Sci, Tallahassee, FL 32306 USA
基金
美国国家科学基金会;
关键词
Shallow water equations; Proper orthogonal decomposition; Reduced-order models (ROMs); Finite difference methods; Discrete empirical interpolation method (DEIM); INTERPOLATION METHOD; POD; FLOW; DEIM;
D O I
10.1016/j.jcp.2012.11.035
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In the present paper we consider a 2-D shallow-water equations (SWE) model on a beta-plane solved using an alternating direction fully implicit (ADI) finite-difference scheme on a rectangular domain. The scheme was shown to be unconditionally stable for the linearized equations. The discretization yields a number of nonlinear systems of algebraic equations. We then use a proper orthogonal decomposition (POD) to reduce the dimension of the SWE model. Due to the model nonlinearities, the computational complexity of the reduced model still depends on the number of variables of the full shallow - water equations model. By employing the discrete empirical interpolation method (DEIM) we reduce the computational complexity of the reduced order model due to its depending on the nonlinear full dimension model and regain the full model reduction expected from the POD model. To emphasize the CPU gain in performance due to use of POD/DEIM, we also propose testing an explicit Euler finite difference scheme (EE) as an alternative to the ADI implicit scheme for solving the swallow water equations model. We then proceed to assess the efficiency of POD/DEIM as a function of number of spatial discretization points, time steps, and POD basis functions. As was expected, our numerical experiments showed that the CPU time performances of POD/DEIM schemes are proportional to the number of mesh points. Once the number of spatial discretization points exceeded 10000 and for 90 DEIM interpolation points, the CPU time decreased by a factor of 10 in case of POD/DEIM implicit SWE scheme and by a factor of 15 for the POD/DEIM explicit SWE scheme in comparison with the corresponding POD SWE schemes. Moreover, our numerical tests revealed that if the number of points selected by DEIM algorithm reached 50, the approximation errors due to POD/DEIM and POD reduced systems have the same orders of magnitude, thus supporting the theoretical results existing in the literature. (c) 2012 Elsevier Inc. All rights reserved.
引用
收藏
页码:95 / 114
页数:20
相关论文
共 30 条
[1]   Inverse Shallow-Water Flow Modeling Using Model Reduction [J].
Altaf, Muhammad Umer ;
Heemink, Arnold W. ;
Verlaan, Martin .
INTERNATIONAL JOURNAL FOR MULTISCALE COMPUTATIONAL ENGINEERING, 2009, 7 (06) :577-594
[3]   An 'empirical interpolation' method: application to efficient reduced-basis discretization of partial differential equations [J].
Barrault, M ;
Maday, Y ;
Nguyen, NC ;
Patera, AT .
COMPTES RENDUS MATHEMATIQUE, 2004, 339 (09) :667-672
[4]   Reduced-order modeling of the upper tropical Pacific Ocean model using proper orthogonal decomposition [J].
Cao, Yanhua ;
Zhu, Jiang ;
Luo, Zhendong ;
Navon, I. M. .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2006, 52 (8-9) :1373-1386
[5]  
Chaturantabut S., 2009, Dimension reduction for unsteady nonlinear partial differential equations via empirical interpolation methods
[6]   A STATE SPACE ERROR ESTIMATE FOR POD-DEIM NONLINEAR MODEL REDUCTION [J].
Chaturantabut, Saifon ;
Sorensen, Danny C. .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2012, 50 (01) :46-63
[7]   Application of POD and DEIM on dimension reduction of non-linear miscible viscous fingering in porous media [J].
Chaturantabut, Saifon ;
Sorensen, Danny C. .
MATHEMATICAL AND COMPUTER MODELLING OF DYNAMICAL SYSTEMS, 2011, 17 (04) :337-353
[8]   NONLINEAR MODEL REDUCTION VIA DISCRETE EMPIRICAL INTERPOLATION [J].
Chaturantabut, Saifon ;
Sorensen, Danny C. .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2010, 32 (05) :2737-2764
[9]   A dual-weighted approach to order reduction in 4DVAR data assimilation [J].
Daescu, D. N. ;
Navon, I. M. .
MONTHLY WEATHER REVIEW, 2008, 136 (03) :1026-1041
[10]   Comparative study with data assimilation experiments using proper orthogonal decomposition method [J].
Dimitriu, Gabriel ;
Apreutesei, Narcisa .
LARGE-SCALE SCIENTIFIC COMPUTING, 2008, 4818 :393-+