The baroclinic-barotropic mode splitting technique is commonly employed in numerical solutions of the primitive equations for ocean modeling to deal with the multiple time scales of ocean dynamics.In this paper,a seco...The baroclinic-barotropic mode splitting technique is commonly employed in numerical solutions of the primitive equations for ocean modeling to deal with the multiple time scales of ocean dynamics.In this paper,a second-order implicit-explicit(IMEX)scheme is proposed to advance the baroclinic-barotropic split system.Specifically,the baroclinic mode and the layer thickness of fluid are evolved explicitly via the second-order strong stability preserving Runge-Kutta scheme,while the barotropic mode is advanced implicitly using the linearized Crank-Nicolson scheme.At each time step,the baroclinic velocity is first computed using an intermediate barotropic velocity.The barotropic velocity is then corrected by re-advancing the barotropic mode with an improved barotropic forcing.Finally,the layer thickness is updated by coupling the baroclinic and barotropic velocities together.In addition,numerical inconsistencies on the discretized sea surface height caused by the mode splitting are alleviated via a reconciliation process with carefully calculated flux deficits.Temporal truncation error is also analyzed to validate the second-order accuracy of the scheme.Finally,two benchmark tests from the MPAS-Ocean platform are conducted to numerically demonstrate the performance of the proposed IMEX scheme.展开更多
基金partially supported by the U.S.Department of Energy,Office of Science,Office of Biological and Environmental Research through Earth and Environmental System Modeling and Scientific Discovery through Advanced Computing programs under university grants DE-SC0020270 and DE-SC0020418partially supported by Shandong Excellent Young Scientists Program(Overseas)under the grant 2023HWYQ-064OUC Youth Talents Project.
文摘The baroclinic-barotropic mode splitting technique is commonly employed in numerical solutions of the primitive equations for ocean modeling to deal with the multiple time scales of ocean dynamics.In this paper,a second-order implicit-explicit(IMEX)scheme is proposed to advance the baroclinic-barotropic split system.Specifically,the baroclinic mode and the layer thickness of fluid are evolved explicitly via the second-order strong stability preserving Runge-Kutta scheme,while the barotropic mode is advanced implicitly using the linearized Crank-Nicolson scheme.At each time step,the baroclinic velocity is first computed using an intermediate barotropic velocity.The barotropic velocity is then corrected by re-advancing the barotropic mode with an improved barotropic forcing.Finally,the layer thickness is updated by coupling the baroclinic and barotropic velocities together.In addition,numerical inconsistencies on the discretized sea surface height caused by the mode splitting are alleviated via a reconciliation process with carefully calculated flux deficits.Temporal truncation error is also analyzed to validate the second-order accuracy of the scheme.Finally,two benchmark tests from the MPAS-Ocean platform are conducted to numerically demonstrate the performance of the proposed IMEX scheme.