Spurious Dianeutral Mixing in a Global Ocean Model Using Spherical Centroidal Voronoi Tessellations

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 2＂, from the perspectives of resolution, viscosity, Horizontal Grid Reynolds Number （HGRN）, Rea, 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., MITocM, 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 MIT6c 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 hexa- gon 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.

Stability-Controllable Second-Order DifferenceScheme for Convection Term

A new finite difference scheme-SCSD scheme has been proposed based on CD (Central Difference)scheme and SUD (Secondr-order Upwind Difference) scheme. Its basic feature is controIlable convectivesfability and always second-order accuracy (Stability-Controllable Second-order Difference ). It hasbeen proven that this scheme is convective-stable if the grid Peclet number .The advanage of this new scheme has been discussed based on the modified wavenumber analysis byusing Fourier transform. This scheme has been applied to the 2-D incompressible convective-diffusiveequation and 2-D compressible Euler equation, and corresponding finite difference equations have beenderived. Numerical examples of 1-D Burgers equation and 2-D transport equation have been presentedto show its good accuracy and controllable convective stability

Numerical simulation of flow and heat transfer in a random packed bed

Random packed beds have more complex interior structure than structured beds and are widely used in industry and engineering. CFD simulation was carried out to investigate and analyze the local flow and heat transfer in a 120-sphere random packed bed. 3D Navier-Stokes equation was solved with a finite volume formulation based on the Chimera meshing technique. Investigation was focused on low Reynolds number flow （Re=4.6-56.2）, which typically occurs in packed bed reactors in bio-chemical fields. Detailed temperature field information was obtained. Inhomogeneity of flow and heat transfer due to the non-uniform distribution of void fraction was discussed and analyzed.