Development and Application of Turbulence Models in Numerical Geodynamo Simulations

湍流模型在数值地球发电机模拟中的发展与应用

• 批准号: 1045277
• 负责人: Bruce Buffett
• 金额: $ 26.66万
• 依托单位: University of California-Berkeley
• 依托单位国家: 美国
• 项目类别: Continuing Grant
• 财政年份: 2011
• 资助国家: 美国
• 起止时间: 2011-01-15 至 2012-12-31
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=1045277&HistoricalAwards=false
• 关键词: Development; Application; Turbulence; Models; Numerical

Fluid motion in the liquid iron core of the Earth generates a magnetic field, which protects the surface environment from the harmful effects of solar wind and cosmic rays. The low viscosity of liquid iron (comparable to that of liquid water) ensures that the fluid motion is highly turbulent and time dependent. Accurate numerical simulations of the fluid motion and field generation are beyond the limit of current computational power, so modelers are forced to use values for the fluid viscosity and other physical parameters that are very far from Earth-like values. As a consequence, the nature of the dynamics is altered and our ability to address important scientific questions is limited. Only modest improvements in the models are expected through brute-force improvements in resolution, so there is a pressing need to develop better models for dealing with turbulence in the core.We propose to develop a second generation of geodynamo model that advances the current models in two ways. First, the new model will employ a method of solution that is better suited to parallel computing. Second, the new model will include a more sophisticated treatment of small-scale turbulence. To achieve these goals we will use the finite-element method to solve the governing equations and implement an adaptive, scale-similarity model to account for the effects of turbulence. The finite-element method reduces the amount of communication on a parallel computer and allows us to locally refine the resolution of the calculation, but only where it is needed. The scale-similarity method is a proven technique that has been successfully implemented in a plane-layer dynamo model to account for complex interactions between the flow and the magnetic field at small scales. We conservatively estimate that the introduction of the scale-similarity model in the geodynamo model will allow us to lower the fluid viscosity by two orders of magnitude, relative to current capabilities. This may be sufficient to reach the dynamical regime (e.g. a Taylor state) where the Earth's geodynamo is thought to operate. The new model will also allow us to probe the geophysically relevant behavior of dynamos at low magnetic Prandtl number Pm, which represents the ratio of fluid viscosity to magnetic diffusivity. Current state-of-the-art numerical simulations assume Pm 1, whereas the expected value for the Earth's core is Pm ~ 1.0e-6 . The operation of dynamos at low Pm raises a number of important scientific questions and poses a substantial computational challenge.

地球液态铁核中的流体运动产生磁场，保护地表环境免受太阳风和宇宙射线的有害影响。液态铁的低粘度（与液态水的粘度相当）确保了流体运动是高度湍流的和时间依赖的。对流体运动和场生成的精确数值模拟超出了当前计算能力的极限，因此建模者被迫使用与地球值相差甚远的流体粘度和其他物理参数值。因此，动态的性质被改变，我们解决重要科学问题的能力受到限制。只有适度的改进，预计通过蛮力提高分辨率，所以有迫切需要开发更好的模型来处理湍流的核心。我们建议开发第二代地球发电机模型，在两个方面推进当前的模型。首先，新模型将采用更适合并行计算的求解方法。其次，新模型将包括对小尺度湍流的更复杂处理。为了实现这些目标，我们将使用有限元方法来解决控制方程，并实现一个自适应的，规模相似的模型来考虑湍流的影响。 有限元方法减少了并行计算机上的通信量，并允许我们局部地改进计算的分辨率，但仅限于需要的地方。 尺度相似性方法是一种经过验证的技术，已成功地在平面层发电机模型中实现，以解释小尺度下流动和磁场之间的复杂相互作用。我们保守地估计，在地球发电机模型中引入尺度相似性模型将使我们能够将流体粘度降低两个数量级，相对于目前的能力。这可能足以达到地球的地球发电机被认为运行的动力学状态（例如泰勒状态）。新的模型还将使我们能够探测发电机在低磁普朗特数Pm，这代表了流体粘度的磁扩散率的比率的相关行为。 目前最先进的数值模拟假设Pm为1，而地球核心的预期值为Pm ~ 1.0e-6。发电机在低PM下的运行提出了一些重要的科学问题，并提出了一个重大的计算挑战。