Numerical Methods for Layered Models of Ocean Circulation

该项目的目标是为一般海洋环流的分层（等密度线）模型开发数值方法。 在这种模型中，垂直坐标不是线性距离，而是潜在的密度或其他一些相关的量。 当等密度线模式在垂直方向上离散化时，其效果是将海洋表示为近似不混溶的层的堆叠。 这与海洋的上部和下部区域大致相互隔离的观察结果相对应。 然而，这种隔离并不完全，在等密度模型中，确实发生的交换是在建模者的控制下。 这种类型的模式表现出很大的承诺，代表海洋动力学的长期行为。 在本项目中，首席研究员将开发和纳入用于求解通量形式的动量方程的方法，其中动量而不是速度是因变量。 看来，动量公式可以给出更可靠的结果，例如在层厚度趋于零的情况下，因为层界面与流体域的顶部或底部相交。 这种方法在模拟海洋中的垂直运输方面也可能是有价值的。 将使用非振荡或几乎非振荡的平流算法来求解动量方程。 这种方法将被分析，并纳入一个两级的时间步进算法，研究人员最近开发的制剂的控制方程中的快速和慢速的时间尺度被分成单独的子系统。 由此产生的算法将进行测试的模型问题的解决方案的行为是已知的。 如果时间允许，PI将在现有的海洋环流操作模型中实施和测试这些方法。 在本项目期间，PI将与俄勒冈州州立大学和洛斯阿拉莫斯国家实验室的海洋模型专家保持定期联系，这项工作的长期目标是促进对全球气候系统的了解。 理解这一系统的一个重要工具是计算机模拟，包括大气、海洋、海冰和陆地效应的耦合模型。 通过这种模拟，人们可以测试地球气候对外部强迫变化的反应，例如温室气体排放量的增加。 世界上的海洋储存和运输大量的热能，因此海洋在我们的气候系统中发挥着至关重要的作用。 当在计算机上模拟海洋环流时，一个是使用计算机获得描述流体流动演变的数学方程的近似解。 本项目的目标是改进用于求解这些方程的方法。

The goal of this project is to develop numerical methods for layered (isopycnal) models of the general ocean circulation. In such models, the vertical coordinate is not linear distance, but instead is potential density or some other related quantity. When an isopycnal model is discretized in the vertical, the effect is to represent the ocean as a stack of layers that are approximately immiscible. This corresponds to the observation that the upper and lower regions of the ocean are approximately isolated thermodynamically from each other. However, this isolation is not complete, and in an isopycnal model the exchanges that do occur are under the control of the modeler. This type of model shows great promise for representing the long-term behavior of ocean dynamics. In the present project, the principal investigator will develop and incorporate methods for solving the momentum equation in a flux form for which momentum, not velocity, is the dependent variable. It appears that the momentum formulation may give more reliable results, such as in situations where layer thicknesses tend to zero because of layer interfaces intersecting the top or bottom of the fluid domain. This approach may also be valuable in modeling vertical transport in the ocean. Non-oscillatory or nearly non-oscillatory advection algorithms will be used to solve the momentum equation. Such methods will be analyzed and incorporated into a two-level time-stepping algorithm that the investigator has recently developed for formulations of the governing equations in which the fast and slow time scales are split into separate subsystems. The resulting algorithms will be tested on model problems for which the solution behavior is known. As time permits, the PI will then pursue the implementation and testing of these methods in existing, operational models of ocean circulation. During this project the PI will be in regular contact with ocean modelers at Oregon State University and at Los Alamos National Laboratory.The long-term objective of this work is to contribute to an understanding of the global climate system. An important tool towards understanding this system consists of computer simulations involving coupled models of the atmosphere, ocean, sea ice, and terrestrial effects. With such simulations one can test the response of the earth's climate to changes in external forcing, such as increased emissions of greenhouse gases. The world's oceans store and transport tremendous amounts of heat energy, so the oceans play a crucial role in our climate system. When the circulation of the ocean is simulated on a computer, one is using the computer to obtain approximate solutions to mathematical equations that describe the evolution of fluid flows. The goal of the present project is to improve the methods that are used to solve those equations.