In order to quantitatively evaluate the spurious dianeutral mixing in a global ocean model MPAS-Ocean(Model for Prediction Across Scales) using a spherical centroidal voronoi tessellations developed jointly by the Nat...In order to quantitatively evaluate the spurious dianeutral mixing in a global ocean model MPAS-Ocean(Model for Prediction Across Scales) using a spherical centroidal voronoi tessellations developed jointly by the National Center for Atmospheric Research and the Los Alamos National Laboratory in the United States, we choose z* vertical coordinate system in MPAS-Ocean, in which all physical mixing processes, such as convection adjustment and explicit diffusion parameter schemes, are omitted, using a linear equation of state. By calculating the Reference Potential Energy(RPE), front revolution position, time rate of RPE change, probability density function distribution and dimensionless parameter χ, from the perspectives of resolution, viscosity, Horizontal Grid Reynolds Number(HGRN), Re?, and momentum transmission scheme, using two ideal cases, overflow and baroclinic eddy channel, we qualitatively analyze the simulation results by comparison with the three non-isopycnal models in Ilicak et al.(2012), i.e., MITGCM, MOM, and ROMS. The results show that the spurious dianeutral mixing in the MPAS-Ocean increases over time. The spurious dianeutral transport is proportional to the HGRN directly and is reduced by increasing the lateral viscosity or using a finer resolution to control HGRN. When the HGRN is less than 10, spurious transport is reduced significantly. When using the proper viscosity closure, MPAS-Ocean performs better than MITGCM and MOM, closely to ROMS, in the 2D case without rotation, and much better than the above-mentioned three ocean models under the condition of 3D space with rotation due to the cell area difference between the hexagon cell and the quadrilateral cell with the same resolution. Both the Zalesak(1979) flux corrected transport scheme and Leith closure in MPAS-Ocean play an excellent role in reducing spurious dianeutral mixing. The performance of Leith scheme is preferable to the condition of three-dimensional baroclinic eddy.展开更多
This paper puts forth a simplified dynamic modeling strategy for the eddy viscosity coefficient parameterized in space and time.The eddy viscosity coefficient is dynamically adjusted to the local structure of the flow...This paper puts forth a simplified dynamic modeling strategy for the eddy viscosity coefficient parameterized in space and time.The eddy viscosity coefficient is dynamically adjusted to the local structure of the flow using two different nonlinear eddy viscosity functional forms to capture anisotropic dissipation mechanism,namely,(i)the Smagorinsky model using the local strain rate field,and(ii)the Leith model using the gradient of the vorticity field.The proposed models are applied to the one-layer and two-layer wind-driven quasigeostrophic ocean circulation problems,which are standard prototypes of more realistic ocean dynamics.Results show that both models capture the quasi-stationary ocean dynamics and provide the physical level of eddy viscosity distribution without using any a priori estimation.However,it is found that slightly less dissipative results can be obtained by using the dynamic Leith model.Two-layer numerical experiments also reveal that the proposed dynamic models automatically parameterize the subgrid-scale stress terms in each active layer.Furthermore,the proposed scale-aware models dynamically provide higher values of the eddy viscosity for smaller resolutions taking into account the local resolved flow information,and addressing the intimate relationship between the eddy viscosity coefficients and the numerical resolution employed by the quasigeostrophic models.展开更多
基金supported by the National Natural Science Foundation of China (Grant No.41175089)
文摘In order to quantitatively evaluate the spurious dianeutral mixing in a global ocean model MPAS-Ocean(Model for Prediction Across Scales) using a spherical centroidal voronoi tessellations developed jointly by the National Center for Atmospheric Research and the Los Alamos National Laboratory in the United States, we choose z* vertical coordinate system in MPAS-Ocean, in which all physical mixing processes, such as convection adjustment and explicit diffusion parameter schemes, are omitted, using a linear equation of state. By calculating the Reference Potential Energy(RPE), front revolution position, time rate of RPE change, probability density function distribution and dimensionless parameter χ, from the perspectives of resolution, viscosity, Horizontal Grid Reynolds Number(HGRN), Re?, and momentum transmission scheme, using two ideal cases, overflow and baroclinic eddy channel, we qualitatively analyze the simulation results by comparison with the three non-isopycnal models in Ilicak et al.(2012), i.e., MITGCM, MOM, and ROMS. The results show that the spurious dianeutral mixing in the MPAS-Ocean increases over time. The spurious dianeutral transport is proportional to the HGRN directly and is reduced by increasing the lateral viscosity or using a finer resolution to control HGRN. When the HGRN is less than 10, spurious transport is reduced significantly. When using the proper viscosity closure, MPAS-Ocean performs better than MITGCM and MOM, closely to ROMS, in the 2D case without rotation, and much better than the above-mentioned three ocean models under the condition of 3D space with rotation due to the cell area difference between the hexagon cell and the quadrilateral cell with the same resolution. Both the Zalesak(1979) flux corrected transport scheme and Leith closure in MPAS-Ocean play an excellent role in reducing spurious dianeutral mixing. The performance of Leith scheme is preferable to the condition of three-dimensional baroclinic eddy.
文摘This paper puts forth a simplified dynamic modeling strategy for the eddy viscosity coefficient parameterized in space and time.The eddy viscosity coefficient is dynamically adjusted to the local structure of the flow using two different nonlinear eddy viscosity functional forms to capture anisotropic dissipation mechanism,namely,(i)the Smagorinsky model using the local strain rate field,and(ii)the Leith model using the gradient of the vorticity field.The proposed models are applied to the one-layer and two-layer wind-driven quasigeostrophic ocean circulation problems,which are standard prototypes of more realistic ocean dynamics.Results show that both models capture the quasi-stationary ocean dynamics and provide the physical level of eddy viscosity distribution without using any a priori estimation.However,it is found that slightly less dissipative results can be obtained by using the dynamic Leith model.Two-layer numerical experiments also reveal that the proposed dynamic models automatically parameterize the subgrid-scale stress terms in each active layer.Furthermore,the proposed scale-aware models dynamically provide higher values of the eddy viscosity for smaller resolutions taking into account the local resolved flow information,and addressing the intimate relationship between the eddy viscosity coefficients and the numerical resolution employed by the quasigeostrophic models.