海洋模式中潮致混合调和分析参数化方案的建立

潮致混合对海洋环流的调整起着重要作用。陆架环流的数值模拟中如果不考虑潮汐作用,往往不能得到与观测相符的垂向温盐结构。本文基于调和分析方法,建立了一套潮致混合参数化方案。该方案通过对垂向混合系数进行调和分析,从而得到随时间变化的潮致混合系数。将该方案用于黄海冷水团数值模拟的结果显示,其能够得到与在数值模式开边界直接加入潮汐强迫相当的冷水团温盐结构。和直接引入潮汐强迫相比,这一潮致混合参数化方案的优势在于,它能够大大节省数值模拟计算机时,因此有望显著提高大规模高分辨率的海洋环流及气候模式的模拟能力和计算效率。     

Langmuir环流在雷诺平均型模式中的参数化研究进展

Langmuir环流影响着海洋上层的能量输入,对海洋上混合层的形成和加深起着重要作用,对于海洋上混合层具有重要意义。近年来许多学者采用大涡模拟(LES)方法对Langmuir环流进行机制研究,并通过在雷诺平均模型中参数化Langmuir环流效应,将Langmuir环流过程引入到三维海洋环流或海洋耦合模式中,提出了一系列混合参数化方案。本文回顾了Langmuir环流在雷诺平均模式参数化中的研究进展,主要可分为以下几种方案:一种方法是用Langmuir数在KPP垂直混合参数化方案中引入湍流特征速度增强因子,并不断发展Langmuir数的定义;一种是在Mellor-Yamada2.5湍流闭合模型中增加斯托克斯漂流剪切效应项,此外还有通过修改模式中混合长方程来加入Langmuir效应等。通过在雷诺平均模式中应用的结果来看,现有的参数化方案在一定程度上改善了混合层深度和SST的模拟,肯定了Langmuir环流在加深混合等方面的作用,但仍存在一些问题需要在今后的研究中进一步改进。     

适应极地快速变化海冰模式的研发与挑战

极地海冰是地球气候系统的重要组成部分,也是气候环境变化的指示器和放大器.极地海冰复杂的多尺度物理过程和极地观测资料的匮乏,给海冰模式的研发带来了巨大的挑战.在过去的半个多世纪中,大气-海冰-海洋的复杂相互作用和冰内物理过程在海冰模式中的数学描述取得了重大的进展,但海冰模式对一些重要物理过程的描述仍很不完善,尤其是近年来...     

重力波非共振激发的数值研究

应用二阶算子分裂格式,采用全非线性模式,对可压大气中重力波波包非共振激发进行数值模拟和分析.研究表明：两个向上传播的重力波波包,通过非共振相互作用,也能激发出一个向下传播的重力波波包.生成波的能量几乎正比于主波的能量,也随次波的能量增大而增大.能量传递方向主要表现为从高频的主波向次波和生成波转移,这与共振相互作用不完全相同.非共振作用的特征时间由主次波的初始位置和波长频率决定,与主次波的初始振幅无关.生成波的波矢量并不严格满足弱非线性相互作用理论预言的波矢量匹配条件,而且生成波的波长和频率在传播的过程中会随时间发生变化,这表明了弱非线性相互作用理论的局限性.     

长清大学园区高架桥关键工序的施工方案

长清大学园区高架桥全长3.533 km,是济菏高速公路工程14个合同段中唯一一个全部是桥梁结构的标段。按照本工程的地理环境、地势条件、设计特点以及业主对工期和质量的要求,编制了一份以灌注桩和现浇连续箱梁为主要工序的施工方案。     

胶州湾多分潮漫滩数值模拟研究

在普林斯顿海洋模式(POM)数值模型的基础上加入了漫滩格式,通过对胶州湾1个月实测潮位资料的调和分析,选取了5个较大的分潮,首次对胶州湾进行了多分潮的潮位和潮流数值模拟,模拟结果与实测资料符合较好。     

风能输入参数对珠江口波浪模拟的影响

珠江口是粤港澳大湾区的核心区域,台风、风暴潮和巨浪等海洋灾害频发,对沿岸人民的生命财产安全构成严重威胁。准确的高分辨率波浪模拟/预报对区域经济建设和防灾减灾具有重要意义。波浪预报/模拟的质量很大程度上取决于风能输入的误差。本研究基于WAVEWATCH III(WW3)的南海-珠江口双重嵌套精细化海浪模式,探讨不同的风场产品与风能输入参数方案组合对珠江口波浪动力过程模拟的影响,确定最优的风场和参数方案组合。ERA5风场更适合珠江口海域的风浪模拟,其模式结果略优于采用CFSR风场的模式结果。ERA5风场+T500方案的组合对珠江口波高变化过程的模拟效果最好,ERA5风场+T471方案的组合次之,ERA5风场+ST6方案再次之。CFSR风场与T471f参数方案最为适配,其结果稍差于ERA5风场+ST6参数方案。T500方案调整高风速下的风能输入和涌浪对风能输入的反馈作用,并考虑水深引起的波浪破碎效应,更适合水深限制的珠江口浅水区域。另外,WW3模式开关ST4的参数方案的表现优于开关ST6的参数方案。     

Modeling piecewise planar fault discontinuities without element-partitioning in 3D reservoir-geomechanical models

Using the extended finite element method, faults can be introduced into a three-dimensional reservoir-geomechanical model without meshing to assess the potential for fault reactivation associated with industrial activities such as disposal of CO₂ or fluid extraction. The residuals to the governing equations include basis functions formed from the enrichment functions for strong and weak discontinuities. The traditional spatial integration scheme is based on an algorithm that partitions each enriched element into a collection of subtetrahedra, which is time and storage intensive. To avoid element partitioning, we adopt the homogeneous numerical integration scheme (referred as HNI hereafter) where the integration of homogeneous monomials over each polyhedron is converted into the integration of the same monomials over the one-dimensional edges of the polyhedron by using Stokes's theorem and Euler's homogeneous function theorem. The integrands for the strong and weak discontinuities are derived to implement the integration scheme in the three-dimensional reservoir-geomechanical model. An example of a jagged fault is presented to reveal the advantages of HNI scheme over the standard element partitioning approach. Several other examples that involve fluid flow, fault sliding, and fault sliding triggered by injection pressure are also presented to demonstrate that accurate and efficient computations are realized by the new integration scheme.     

Numerical stability for modelling of dynamic two‐phase interaction

Dynamic two‐phase interaction of soil can be modelled by a displacement‐based, two‐phase formulation. The finite element method together with a semi‐implicit Euler–Cromer time‐stepping scheme renders a discrete equation that can be solved by recursion. By experience, it is found that the CFL stability condition for undrained wave propagation is not sufficient for the considered two‐phase formulation to be numerically stable at low values of permeability. Because the stability analysis of the two‐phase formulation is onerous, an analysis is performed on a simplified two‐phase formulation that is derived by assuming an incompressible pore fluid. The deformation of saturated porous media is now captured in a single, second‐order partial differential equation, where the energy dissipation associated with the flow of the fluid relative to the soil skeleton is represented by a damping term. The paper focuses on the different options to discretize the damping term and its effect on the stability criterion. Based on the eigenvalue analyses of a single element, it is observed that in addition to the CFL stability condition, the influence of the permeability must be included. This paper introduces a permeability‐dependent stability criterion. The findings are illustrated and validated with an example for the dynamic response of a sand deposit. Copyright © 2015 John Wiley & Sons, Ltd.     

沙氏指数计算方案探讨

从沙氏指数的物理定义出发,给出了计算沙氏指数的4种方案并进行了误差分析。A方案没有考虑水汽凝结潜热,计算结果误差大;B方案虽然考虑水汽凝结潜热,但计算方法没有考虑湿绝热线的曲率等问题,计算结果虽优于A方案,但误差仍较大;C方案不仅考虑水汽凝结潜热,并从计算方法上作了改进,不仅考虑了湿绝热线的曲率等问题,同时避免了由气压直接求取气温造成误差较大的问题,计算结果优于A、B方案,比较理想,可在实际工作中使用;D方案由假绝热方程经过转化求得上升气块的温度,然后算出沙氏指数,虽然精确,但计算复杂,实际应用较困难。