A communication-avoiding implicit-explicit method for a free-surface ocean model

详细信息    查看全文

• 作者：Christopher Newman^a ; ^{cnewman@lanl.gov" class="auth_mail" title="E-mail the corresponding author} ; Geoffrey Womeldorff^b ; ^{womeld@lanl.gov" class="auth_mail" title="E-mail the corresponding author} ; Dana A. Knoll^c ; ^{nol@lanl.gov" class="auth_mail" title="E-mail the corresponding author} ; Luis Chacó ; n^d ; ^{chacon@lanl.gov" class="auth_mail" title="E-mail the corresponding author}
• 关键词：Ocean modeling ; IMEX method ; Jacobian-free Newton&ndash ; Krylov methods ; Nonlinear elimination ; Communication-avoiding algorithm
• 刊名：Journal of Computational Physics
• 出版年：2016
• 出版时间：15 January 2016
• 年：2016
• 卷：305
• 期：Complete
• 页码：877-894
• 全文大小：953 K

文摘

We examine a nonlinear elimination method for the free-surface ocean equations based on barotropic–baroclinic decomposition. The two dimensional scalar continuity equation is treated implicitly with a preconditioned Jacobian-free Newton–Krylov method (JFNK). The remaining three dimensional equations are subcycled explicitly within the JFNK residual evaluation with a method known as nonlinear elimination. In this approach, the memory footprint of the underlying Krylov vector is greatly reduced over that required by fully coupled implicit methods. The method is second-order accurate and scales algorithmically, with allowed timesteps much larger than fully explicit methods. Moreover, the hierarchical nature of the algorithm lends itself readily to emerging architectures. In particular, we introduce a communication staging strategy for the three dimensional explicit system that greatly reduces the communication costs of the algorithm and provides a key advantage as communication costs continue to dominate relative to floating point costs in emerging architectures.