Second-order Stokes wave solutions for intefacial waves in three-layer stratified fluid with background current

In this paper, interfacial waves in three-layer stratified fluid with background current are investigated using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory, and the Kelvin-Helmholtz instability of interfacial waves is studied. As expected, for three-layer stratified fluid with background current, the first-order asymptotic solutions (linear wave solutions), dispersion relation and the second-order asymptotic solutions derived depend on not only the depths and densities of the three-layer fluid but also the background current of the fluids, and the second-order Stokes wave solutions of the associated elevations of the interfacial waves describe not only the second-order nonlinear wave-wave interactions between the interfacial waves but also the second-order nonlinear interactions between the interfacial waves and currents. It is also noted that the solutions obtained from the present work include the theoretical results derived by Chen et al (2005) as a special case. It also shows that with the given wave number k (real number) the interfacial waves may show Kelvin-Helmholtz instability.

孟加拉湾表层盐度季节变化和ENSO循环的拓展伴随模态分析

本文采用旋转经验正交函数分析和拓展伴随模态分析方法研究孟加拉湾表层盐度的季节变化，结果表明其主要特征是西北淡水羽（位于湾的西北部）和东北淡水羽（位于湾的东北部）的季节性演变。西北淡水羽始于7月初，8月底达到峰值，10月底转变为湾形淡水羽；湾形淡水羽集中分布在孟加拉湾的西北和东北沿岸附近，并在第2年1月转变为东北淡水羽。 孟加拉湾北部的环流（地转流+Ekman漂流）对淡水羽的形成和演变起重要作用。Sverdrup环流，代表对风应力的正压响应，能部分解释除了秋季湾形淡水羽的东支之外淡水羽的形成机制和演变过程。海表淡水通量与淡水羽在空间上分布不一致：当孟加拉湾东北部的海表淡水通量急剧增加（减少）时，西北淡水羽（湾形淡水羽的东支）产生，所以海表淡水通量不是淡水羽的形成和演变的主要原因。此外，孟加拉湾周边河流的径流量可能只是淡水羽的主要淡水来源，而与淡水羽的形成和演变关系不大。 初步研究了孟加拉湾北部盐度对海面高度异常的影响，结果表明Ekman 抽吸对孟加拉湾的海面高度异常起着主要作用，热量通量和淡水通量起次要作用。7－12月，孟加拉湾北部的淡水通量决定北部的海面比容高度异常；1－3月，热量通量对海面比容高度异常起主要作用；4－6 月，热量通量和淡水通量在不同时空对孟加拉湾北部的海面比容高度异常起不同的作用。 同样，采用拓展伴随模态分析方法研究ENSO循环的暖位相——El Niño，表明必须将中低纬度印度洋和太平洋海气系统作为一个整体来研究El Niño；证实了正是直接来自南北中纬度太平洋的异常西风爆发和海面风辐聚所驱动的近赤道异常表层海水东移和经向辐聚，而不是Kelvin波，造就了El Niño现象的增温信号。El Niño事件中海洋和大气具有明显不同的性态：气候系统中并不存在一个相对独立的热带太平洋大气变异，但是的确存在一个相对独立的热带太平洋海洋变异。总体上而言：El Niño事件早期主要是大气驱动海洋，晚期主要是海洋驱动大气。 研究结果显示拓展伴随模态分析方法是研究ENSO循环和海洋水文特征形成机制的有力工具。

太平洋西部边界及边界流对赤道暖水库形成和演化的影响

A model of equatorial ocean is used to study the roles of the Pacific western boundary and the Mindanao Current (MC) in the evolution of the equatorial warm pool. The model consists of the single baroclinic mode of a two-layer ocean, with the parameterization of the anomalous increment of the interface representing the SST difference from its long-term-space-mean. The ocean is driven by a wind path in the middle ocean with a real or an artificial geometry assigned at the western and eastern boundaries. In order to test the role of the MC, the western boundary current is introduced into the model by a boundary condition at a position, real and unreal, respectively. The model experiments show that the warm pool, which is insensitive to the longitudinal width of the wind band in middle ocean, results mainly from the accumulation o the eastly-drifted warm water in the equatorial western Pacific. It is the dominant factor for the formation of the warm pool that, at a very low latitude, the Papua New Guinea coast intersects the longitudinally lined Philippine Islands at an obtuse angle. In contrast, the western Atlantic boundary, which inclines poleward from the equator at some 135 degrees, could guide the warm water there moving to a higher latitude. On the other hand, the equatorial warm pool in the western equatorial Pacific is very sensitive to the assignment of th Mindanao Current at 7.5°N and displaces southward, with a stronger southern branch than the northern one. We attribute this asymmetry to the combined effect of the western boundary and the MC upon the equatorial warm away from the equator. A by-product of our solutions is the possible mechanism of the "secondary warm pool" in the eastern Pacific north of the equator. It is suggested that, mainly or partly, the "secondary warm pool" results from the cooperation of the southeast monsoon in eastern Pacific and the eastern boundary hindering the propagation of the Kelvin wave poleward alongshore.

热带印度洋低频Rossby波及其对上层热结构的影响

Rossby波是地球物理流体动力学中非常重要的一种波动，海洋斜压长Rossby波在海洋动力过程中起着相当重要的作用。它维持并影响强西边界流，是海盆内能量传播的主要机制，它所携带的变异信号从大洋的东边界传播到内部，对海气耦合系统起到很重要的作用。热带印度洋是季风爆发的源地，对季风的年际变化具有重要影响。研究热带印度洋对理解季风变率和提高季风预测水平有重要的科学和应用价值。 本文利用TOPEX/Poseidon等高度计资料、美国国家海洋数据中心（NODC）的世界海洋图集（WOA05）长期气候态水文资料、美国Scripps海洋研究所的上层海温资料、中国Argo资料中心提供的Argo资料、美国国家环境预测中心（NCEP）的海表面温度、FSU（Florida State University）月平均风场和海气界面热通量等观测数据，全面分析了热带印度洋低频Rossby波的基本特征，并深入研究了低频Rossby波的生成机制及其对上层海洋热结构的影响。 采用相关分析等统计方法，结合1.5层约化重力模式，研究了热带南印度洋低频Rossby波的生成机制。结果表明： （1）热带南印度洋低频Rossby波分为东边界扰动产生的Rossby波和南印度洋中部风强迫Rossby波；东边界激发的为自由Rossby波，沿12°S波速大约13 cm/s，向西最远传播到80°E左右，之后被局地变量调整；强迫Rossby波在西传的过程中不断加强，波速较快，沿12°S能超过20 cm/s； （2）东边界扰动由印度尼西亚贯通流（ITF）导致的地转调整过程引起；内区风强迫Rossby波生成和加强的关键区为（70°E–95°E，15°S–5°S）；显著的西传Rossby波同太平洋上的厄尔尼诺/南方涛动（ENSO）事件紧密相连，ENSO通过大气的遥驱动机制激发热带南印度洋低频Rossby波； （3）作为东边界低频变量扰动的一个重要因子，ITF的变化与ENSO事件密切相关，总的来讲，El Niño年ITF偏弱，La Niña年ITF偏强，这与前人的研究结果一致；但它在ENSO的不同位相时期，存在一定差异，并具有夏季锁相特征：El Niño事件发生年的春季到秋季，ITF偏强，夏季最强；从El Niño盛期（冬季）到次年秋季，ITF持续偏弱，夏季最弱。上述夏季锁相特征与夏季风的强弱变化相对应。La Niña期间情况相反。 西南印度洋（SWIO）（50°E–75°E，15°S–5°S）以及苏门答腊-爪哇沿岸地区是海表面高度异常（SSHA）和海表面温度异常（SSTA）相关显著的区域，表明内部海洋动力过程在次表层和表层变量的相互关系中起重要作用。本文以2006–2008年期间三个连续的同El Niño或者La Niña同时发生的正印度洋偶极子（IOD）事件为背景，基于Argo观测资料研究了表层和次表层IOD的演变以及二者的区别和联系，并首次采用垂直模态分解方法探讨了Rossby波对上层海洋热力结构影响的动力学特征，得到如下主要结论： （1）在热带印度洋，海洋动力过程一般主要由第一和第二低阶垂直斜压模态控制，而第一斜压模态处于主导地位——在SWIO海区，第一斜压模态运动的方差解释率为第二模态的2–3倍，在赤道和东南印度洋也达到2倍左右；另外，赤道印度洋地区高阶斜压模态运动对该地区的海洋动力过程也具有一定的贡献； （2）低频斜压Rossby波能影响海洋的垂直层结，尤其是强暖Rossby波使同第二斜压模态运动紧密相连的海洋上层层结减弱，加强第二斜压模态的贡献量，导致上层各等压线向下垂直位移增大，最终通过垂直混合过程调整上层海洋的热力结构；而低频斜压冷Rossby波会加强上层垂直层结，抑制该层内变量变化，因此第二斜压模态的贡献依然很小； （3）表层IOD和次表层IOD分布形态不同：表层东部冷异常主要集中在东南印度洋Sumatra-Java沿岸，次表层冷异常基本关于赤道对称；表层西部暖异常基本关于赤道对称，而在次表层赤道以南海温扰动强度远远大于赤道以北； （4）正IOD事件中，东南印度洋冷SSTA首先出现于Java沿岸，沿岸东南风引起的潜热释放增加以及沿岸上升流是该初始冷异常建立的主要机制，与之相关的SSTA东西梯度加强大气环流变化，并进一步强迫随后的海洋运动；1–2个月后，SST冷异常中心北跳到Sumatra沿岸并向西扩展，同时不断增强，其中Sumatra沿岸上升流、来自赤道印度洋的冷Kelvin波及其反射的西传冷Rossby波是这一演变过程的动力机制，而沿岸上升流起决定作用。

西太平洋暖池纬向变异规律研究

西太平洋暖池汇集着世界上开阔海洋中温度最高的海水，它的维持和变化对全球天气和气候的变化起着关键性的作用。从海洋学角度全面而系统地研究西太平洋暖池的纬向变异特征、机制以及暖池纬向变异对ENSO的影响等科学问题，对深入了解发生在暖池区的海气耦合过程，全面认识暖池变异在ENSO循环及在全球气候变化中的作用是不可缺少的，该项研究具有重要的理论意义和实践意义。 本文根据长时间序列的海洋和大气资料，利用EOF、小波等分析方法较系统地研究了西太平洋暖池整层暖水纬向运移特征、暖池水体的时空振荡及暖池纬向的运移、时空振荡对ENSO的影响，分析了西太平洋暖池和东太平洋暖池对ENSO的联合影响，探讨了风场和流场对西太平洋暖池纬向变异的影响，并用一个简单的海洋模式，对暖池纬向变异的动力机制进行了诊断分析。主要研究内容及结果如下： 利用加权平均法，建立了一个表征整体暖池纬向运移的指标序列，分析了暖池纬向运移的时频特征以及暖池内部暖水在纬向运移上的差异，并探讨了暖池的纬向运移对赤道西太平洋海平面高度和温跃层深度的影响以及与ENSO的关系。结果显示：暖池的纬向运移具有显著的2-7a的年际变化和10-16a的年代际变化，并于1976年前后经历了一次气候跃变。暖池内部大致可以50m深度水层为界分为上、下两部分，暖池的上半部分纬向运移幅度非常大，而下半部分则相对较小。暖池的纬向运移不仅是赤道西太平洋海平面高度和温跃层深度异常变化的一个非常重要的直接因素，而且对ENSO的形成与发展有非常重要的作用。 为进一步了解暖池整体的时空振荡及其影响，分别对暖池厚度场、热带太平洋风应力场、海表面高度和热含量场进行了EOF分析。结果显示，西太平洋暖池水体具有大致以赤道为横轴的经向“跷跷板”式季节性反位相振荡和大致以 170o E为纵轴的纬向“跷跷板”式年际反位相振荡。纬向振荡是暖池厚度异常场年际变率的的主导模态。西太平洋暖池暖水的纬向运移是造成暖池水体纬向振荡的主要原因。热带太平洋海表面高度与热含量异常场有着非常相似的ENSO时间尺度的振荡特征：即大致以5o N为横轴的经向“跷跷板”式反位相振荡和大致以170o W为纵轴的纬向“跷跷板”式反位相振荡。而与ENSO密切相关的热带太平洋风应力异常场也存在着纬向和经向两种主要模态，分别表现为近赤道的纬向风异常和远赤道的经向风异常。西太平洋暖池的纬向振荡对风应力异常存在滞后1个月左右的响应和显著的相互作用关系。西太平洋暖池纬向振荡能够造成热带太平洋水体质量与热量的重新分配，并通过强烈的海气相互作用对ENSO循环形成与发展起着重要的作用。 西太平洋暖池和东太平洋暖池是热带太平洋两个重要的暖水区。在系统分析西太平洋暖池的纬向变异和东太平洋暖池的经向变异特征基础上，提出了两暖池对ENSO循环联合影响的新观点，并建立了一个联合影响指数。结果显示，当联合指数达1.6时，则预示着一次新El Nino的发生。这为ENSO的研究和预测提供了一个新线索。 数据分析显示，赤道中太平洋的纬向风应力及上层纬向流异常是暖池纬向运移的两个重要的动力因素，且两者对暖池的纬向运移有良好的预报意义。利用一个简单的海洋模式模拟结果，分析了风驱动下的上层纬向流异常、波动(Kelvin波和Rossby波)及其在边界的反射效应对暖池纬向运移的动力影响。模式结果显示：中太平洋上层纬向流异常是西太平洋暖池纬向运移的主要驱动机制。暖池的纬向运移对波模态流具有滞后约4个月的响应，而风漂流对暖池纬向运移的影响则比较适时。波动在太平洋东、西两边界的反射效应对暖池西缩的影响较大。

Effects of spark discharge on the anodic coatings on magnesium alloy

Spark discharge was the representative phenomenon of Micro-arc oxidation (MAO) method distinguished from other electrochemical oxidation methods. Under the spark discharge treatment, characteristics of the anodic layer were significantly changed. To investigate the influences of the spark discharge, a piece of magnesium alloy AZ91D specimen was partly treated by MAO method in alkaline silicate solution. And the microstructure, element distributions as well as the surface potential distributions of the specimen were studied by scanning electron microscopy (SEM), energy dispersive spectroscopy (EDS) and scanning Kelvin probe (SKP) technique. As a result of intensive spark discharge treatment, porous external layer with dense internal layer were formed on Mg alloy surface. At the same time, the depositions of OH- and SiO32- ions were accelerated, which resulted in the enrichment of element oxygen and silicon at the spark discharge region. Moreover, due to the compact internal layer, the intensive spark discharge region exhibited more positive potentials with respect to other regions, which meant this region could restrain the ejection of electron and provide effective protection to the substrate. In addition, it was found that oxygen played a vital role in determining the intensity and size of sparks, and abundant oxygen resulted in intensive and larger sparks. (c) 2005 Elsevier B.V. All rights reserved.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.