A new oceanic general circulation model in pressure coordinates is formulated. Since the bottom pressure changes with time, the vertical coordinate is actually a pressure-sigma coordinate. The numerical solution of th...A new oceanic general circulation model in pressure coordinates is formulated. Since the bottom pressure changes with time, the vertical coordinate is actually a pressure-sigma coordinate. The numerical solution of the model is based on an energy-conservation scheme of finite difference. The most important new feature of the model is that it is a truly compressible ocean model and it is free of the Boussinesq approximations. Thus, the new model is quite different from many existing models in the following ways: 1) the exact form of mass conservation, 2) the in-situ instantaneous pressure and the UNESCO equation of state to calculate density, 3) the in-situ density in the momentum. equations, 4) finite difference schemes that conserve the total energy. Initial tests showed that the model code runs smoothly, and it is quite stable. The quasi-steady circulation patterns generated by the new model compare well with existing models, but the time evolution of the new model seems different from some existing models. Thus, the non-Boussinesq models may provide more accurate information for climate study and satellite observations.展开更多
The mass transport velocity in a thin layer of muddy fluid is studied theoretically. The mud motion is driven by a periodic pressure load on the free surface, and the mud is described by a power-law model. Based on th...The mass transport velocity in a thin layer of muddy fluid is studied theoretically. The mud motion is driven by a periodic pressure load on the free surface, and the mud is described by a power-law model. Based on the key assumptions of the shallowness and the small deformation, a perturbation analysis is conducted up to the second order to find the mean Eulerian velocity in an Eulerian coordinate system. The numerical iteration method is adopted to solve these non-linear equations of the leading order. From the numerical results, both the first-order flow fields and the second-order mass transport velocities are examined. The verifications are made by comparing the numerical results with experimental results in the literature, and a good agreement is confirmed.展开更多
文摘A new oceanic general circulation model in pressure coordinates is formulated. Since the bottom pressure changes with time, the vertical coordinate is actually a pressure-sigma coordinate. The numerical solution of the model is based on an energy-conservation scheme of finite difference. The most important new feature of the model is that it is a truly compressible ocean model and it is free of the Boussinesq approximations. Thus, the new model is quite different from many existing models in the following ways: 1) the exact form of mass conservation, 2) the in-situ instantaneous pressure and the UNESCO equation of state to calculate density, 3) the in-situ density in the momentum. equations, 4) finite difference schemes that conserve the total energy. Initial tests showed that the model code runs smoothly, and it is quite stable. The quasi-steady circulation patterns generated by the new model compare well with existing models, but the time evolution of the new model seems different from some existing models. Thus, the non-Boussinesq models may provide more accurate information for climate study and satellite observations.
基金supported by the National Natural Science Foun-dation of China(Grant No.40376028)the Application and Basic research of Tianjin(Grant No.11JCYBJC03200)
文摘The mass transport velocity in a thin layer of muddy fluid is studied theoretically. The mud motion is driven by a periodic pressure load on the free surface, and the mud is described by a power-law model. Based on the key assumptions of the shallowness and the small deformation, a perturbation analysis is conducted up to the second order to find the mean Eulerian velocity in an Eulerian coordinate system. The numerical iteration method is adopted to solve these non-linear equations of the leading order. From the numerical results, both the first-order flow fields and the second-order mass transport velocities are examined. The verifications are made by comparing the numerical results with experimental results in the literature, and a good agreement is confirmed.
