A parameterization scheme of vertical mixing due to inertial internal wave breaking in the ocean general circulation model

Based on the theoretical spectral model of inertial internal wave breaking （fine structure） proposed previ- ously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior be- low the surface mixed layer in the ocean general circulation model （OGCM） is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes （including inertial internal wave breaking product） on oceanic interior mixing is emphasized. We suggest that adding the inertial inter- hal wave breaking mixing scheme （F-scheme for short） put forward in this paper to the turbulence mixing scheme of Canuto et al. （T-scheme for short） in the OGCM, except the region from 15°S to 15°N. The numeri- cal results ofF-scheme by usingWOA09 data and an OGCM （LICOM, LASG/IAP climate system ocean model） over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation （AMOC） by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer.

Numerical Prediction of Added Resistance and Vertical Ship Motions in Regular Head Waves

The numerical prediction of added resistance and vertical ship motions of one ITTC （Intemational Towing Tank Conference） S-175 containership in regular head waves by our own in-house unsteady RANS solver naoe-FOAM-SJTU is presented in this paper. The development of the solver naoe-FOAM-SJTU is based on the open source CFD tool, OpenFOAM. Numerical analysis is focused on the added resistance and vertical ship motions （heave and pitch motions） with four very different wavelengths （ 0.8Lpp 〈 2 〈 1.5L ） in regular head waves. Once the wavelength is near the length of the ship model, the responses of the resistance and ship motions become strongly influenced by nonlinear factors, as a result difficulties within simulations occur. In the paper, a comparison of the experimental results and the nonlinear strip theory was reviewed and based on the findings, the RANS simulations by the solver naoe-FOAM-SJTU were considered competent with the prediction of added resistance and vertical ship motions in a wide range of wave lengths.

Experimental Study on Wave Attenuation Performance of A New Type of Floating Breakwater with Twin Pontoons and Multi Porous Vertical Plates

A floating breakwater(FB)has extensive potential applications in the fields of coastal,offshore,and ocean engineering owing to its advantages such as eco-friendliness,low cost,easy and rapid construction,and quick dismantling and reinstallation.An FB composed of twin pontoons and multi-porous vertical plates is proposed to improve the wave attenuation performance.The wave attenuation performance is investigated for different FB structures and vertical plate types under different incident wave heights and periods using 2D wave physical model tests in a wave flume.The results demonstrate that the proposed FB has a better performance than that of the conventional single pontoon-type FB.It reduces the wave transmission due to its enhanced wave reflection and energy loss.The wave transmission coefficient of the proposed FB decreases with an increase in the number of layers and relative draft depth of the vertical plates.However,a further decrease in the wave transmission coefficient is not observed when the number of porous vertical plates is increased from 4 to 5 layers.An equation has been derived to predict the wave transmission of the proposed FB based on the experimental results.

Simplified Method for Calculating Standing Wave Pressure on Vertical Breakwater

Through numerical modeling, a kind of simplified calculating method for standing wave pressure on vertical face breakwater have been put forward. Not only the formulas proposed in this paper are simple in form and very easy in use, but also they possess continuity on the full range of standing wave. And more, the precision requiremennts of calculation can be satisfied to a certain extent in engineering practice.

Contribution of surface wave-induced vertical mixing to heat content in global upper ocean

Compared with observations,the simulated upper ocean heat content(OHC)determined from climate models shows an underestimation bias.The simulation bias of the average annual water temperature in the upper 300 m is 0.2℃lower than the observational results.The results from our two numerical experiments,using a CMIP5 model,show that the non-breaking surface wave-induced vertical mixing can reduce this bias.The enhanced vertical mixing increases the OHC in the global upper ocean(65°S–65°N).Using non-breaking surface wave-induced vertical mixing reduced the disparity by 30%to 0.14℃.The heat content increase is not directly induced by air-sea heat fluxes during the simulation period,but is the legacy of temperature increases in the first 150 years.During this period,additional vertical mixing was initially included in the climate model.The non-breaking surface wave-induced vertical mixing improves the OHC by increasing the air-sea heat fluxes in the first 150 years.This increase in air-sea heat fluxes warms the upper ocean by 0.05–0.06℃.The results show that the incorporation of vertical mixing induced by nonbreaking surface waves in our experiments can improve the simulation of OHC in the global upper ocean.

Vertical two-dimensional non-hydrostatic pressure model with single layer

The vertical two-dimensional non-hydrostatic pressure models with multiple layers can make prediction more accurate than those obtained by the hydrostatic pres- sure assumption. However, they are time-consuming and unstable, which makes them unsuitable for wider application. In this study, an efficient model with a single layer is developed. Decomposing the pressure into the hydrostatic and dynamic components and integrating the x-momentum equation from the bottom to the free surface can yield a horizontal momentum equation, in which the terms relevant to the dynamic pressure are discretized semi-implicitly. The convective terms in the vertical momentum equation are ignored, and the rest of the equation is approximated with the Keller-box scheme. The velocities expressed as the unknown dynamic pressure are substituted into the continuity equation, resulting in a tri-diagonal linear system solved by the Thomas algorithm. The validation of solitary and sinusoidal waves indicates that the present model can provide comparable results to the models with multiple layers but at much lower computation cost.

The effect of the wave-induced mixing on the upper ocean temperature in a climate model

The significant underestimation of sea surface temperature （SST） and the temperature in the upper ocean is one of common problems in present climate models. The influence of the wave-induced mixing on SST and the temperature in the upper ocean was examined based on a global climate model. The results from the model coupled with wave-induced mixing showed a significant improvement in the simulation of SST and the temperature in the upper ocean compared with those of the original model without wave effects. Although there has still a cold bias, the new simulation is much closer to the climatology, especially in the northern ocean and tropical ocean. This study indicates that some important physical processes in the accurate simulation of the ocean may be ignored in present climate models, and the wave-induced mixing is one of those factors. Thus, the wave-induced mixing （ or the effect of surface waves） should be incorporated properly into climate models in order to simulate or forecast the ocean, then climate system, more accurately.

Seismic wave modeling in viscoelastic VTI media using spectral element method

Spectral element method（SEM） for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries.It is an advanced choice for wave simulations.Due to anelasticity of earth media,SEM for elastic media is no longer appropriate.On fundamental of the second-order elastic SEM,this work takes the viscoelastic wave equations and the vertical transversely isotropic（VTI） media into consideration,and establishes the second-order SEM for wave modeling in viscoelastic VTI media.The second-order perfectly matched layer for viscoelastic VTI media is also introduced.The problem of handling the overlapped absorbed corners is solved.A comparison with the analytical solution in a twodimensional viscoelastic homogeneous medium shows that the method is accurate in the wave-field modeling.Furtherly,numerical validation also presents its great flexibility in solving wave propagation problems in complex heterogeneous media.This second-order SEM with perfectly matched layer for viscoelastic VTI media can be easily applied in wave modeling in a limited region.

波浪能发电装置垂荡浮子PTO系统的建模与计算

为了便捷求出垂荡式波浪能发电装置垂荡浮子的平均功率,对现有的设计计算方法进行了研究。基于波浪理论和Froude-Krylov假定法求解波浪力,应用单自由度有阻尼系统受迫振动理论,将波浪激励力作为对垂荡浮体的输入动力,将发电机的电磁阻力作为主要的阻尼力,建立了质量-弹簧-阻尼振动力学模型,确定了垂荡浮子响应的计算方法,并对实际波浪能发电装置垂荡浮子Power Take Off(PTO)系统列举出计算实例。在算例基础上,对假设条件进行了扩展,应用频率响应曲线、MATLAB软件做出不同参数条件下平均功率三维曲线和输出平均功率数值表,清晰地表明了系统响应状况。研究显示,在波高和频率比一定的条件下,PTO系统平均机械功率随着浮体质量和阻尼力矩的变化具有峰值,围绕峰值区域设计装置的各项参数,能够取得较理想的效果。该建模与计算方法可为实际工程优化设计提供参考。

不同层结条件下海浪生湍流混合特征分析及其对上层海洋的影响

上层海洋混合过程的动力学研究和数值模拟一直是海洋和气象学家关注的重点之一,其中上层海洋垂向混合参数化方案的不准确直接影响海洋模式的模拟效果。风生海浪对海气相互作用以及上层海洋的内部混合都起到极为重要的作用。本文基于包含垂向层结效应的海浪生湍流混合系数理论模型,开展不同层结条件下海浪生湍流混合系数的分布变化特征分析。构建了一维理想解析实验和准全球算例,结果表明,在100 m以浅的上层海洋中,混合层深度与上述混合系数呈现较为显著的负相关关系,当跃层或混合层深度较大时能够显著增强垂向混合过程,进而使表层暖水下降至更深的位置,加深混合层深度模拟结果。主要原因是由于深跃层条件下表层的层结效应较弱,因而对混合的抑制作用较小。此外,在准全球算例中海浪生湍流混合方案可显著提高海洋环流模式对上层海洋温度结构的模拟能力,模式对跃层处WOA13实测数据的模拟精度最高可提升60%以上。

A comparison of two global ocean-ice coupled models with different horizontal resolutions

A global eddy-permitting ocean-ice coupled model with a horizontal resolution of 0.25° by 0.25° is estab- lished on the basis of Modular Ocean Model version 4 （MOM4） and Sea Ice Simulator （SIS）. Simulation results are compared with those of an intermediate resolution ocean-ice coupled model with a horizontal resolution of about 1° by 1°. The results show that the simulated ocean temperature, ocean current and sea ice concentration from the eddy-permitting model are better than those from the intermediate resolu- tion model. However, both the two models have the common problem of ocean general circulation models （OGCMs） that the majority of the simulated summer sea surface temperature （SST） is too warm while the majority of the simulated subsurface summer temperature is too cold. Further numerical experiments show that this problem can be alleviated by incorporating the non-breaking surface wave-induced vertical mixing into the vertical mixing scheme for both eddy-permitting and intermediate resolution models.

塔中地区中晚奥陶世内波、内潮汐沉积

现代海底普遍发育由内波、内潮汐引起的深水牵引流 ,内波、内潮汐不仅可形成各种厘米级的深水牵引流沉积构造 ,还可建造千米级的大型沉积物波。这种大型沉积物波既可沿斜坡向下迁移 ,也可沿斜坡向上迁移。在古代地层记录中已发现各种厘米级的内波、内潮汐沉积单元 ,但尚未发现内波成因的大型沉积物波。塔中地区中上奥陶统碎屑岩段内发育内波及内潮汐形成的各种牵引流沉积构造单元。通过详细的岩心观察 ,识别出 4种内波、内潮汐沉积微相类型及 5种基本垂向沉积层序。此外 ,通过地震剖面分析 ,在研究区中上奥陶统陆坡相中识别出顺坡向上迁移的大型沉积物波 ,该沉积物波的特征与现代海底发育的沉积物波的特征类似 ,运用内波理论可对其成因进行合理的解释。

黏弹性地基中基于虚土桩模型的桩顶纵向振动阻抗研究

基于虚土桩模型,对均质黏弹性地基中桩土纵向耦合振动问题进行了研究。首先,假定桩侧土为各向同性的线性黏弹性材料,并考虑土体的竖向波动效应,结合Euler-Bernoulli杆件理论,建立了桩土纵向耦合振动的定解问题;其次,采用分离变量法求解桩侧土纵向振动的控制方程,得到了桩侧土与桩身接触面上的剪切动刚度,将所得的剪切动刚度代入到桩身振动控制方程,采用Laplace变换技术,进一步求得了任意荷载作用下桩顶纵向振动阻抗的解析解。基于所得解,详细讨论了不同桩身设计参数时桩端土厚度对桩顶纵向振动阻抗的影响。最后,将虚土桩模型与其他桩端土支承模型进行了对比研究,结果表明,对虚土桩模型选用合适的材料参数和桩端土厚度,其得到的桩端支承复刚度值介于现有多种模型的计算值之间。

改进的正演模拟定位原理

本文在作者以前提出的模拟定位原理及其相应算法的基础上，提出了改进的原理──正演模拟定位原理2。改进后的原理是通过认识共炮记录正演模拟的实质，抛开偏移成像的因素，就正演而作正演。因此，在同样用单程声波方程获得完全一样的共炮点记录的正演记录的情况下，使用基于定位原理2的计算公式，可成百倍地提高运算速度。

溴化锂溶液垂直降膜过程中流动与传热的数值模拟

采用流体体积(VOF)法捕捉自由界面,并应用连续表面力(CSF)模型,建立了垂直降膜波动模型,研究了溴化锂溶液垂直降膜过程中的波动过程,根据波动模型得到的界面信息,建立了以波速为动坐标系的固定界面传热模型,同时考察了波动对传热的影响,并建立了垂直降膜波动模型和固定界面传热模型。结果表明:一定入口频率的扰动发展为孤波和毛细波的组合波,孤波波峰处出现了循环流动,且在第一个毛细波波谷处出现回流。波动强化机理与Pr数有关,Pr数较大时,对流起主要作用,Pr数较小时,液膜减薄起主要作用。

垂直地震剖面的有限元正演模拟

应用有限元法和差分法对弹性波波动方程及其定解条件进行数值模拟，得到了无震源形式下的波动方程递推解．引入震源，对有限波场区域进行垂直地震剖面模拟调试，发现空间离散间隔、时间采样间隔与弹性介质速度、震源最大频率之间存在制约关系，满足该关系可使模拟结果收敛、稳定．对４个地层模型（简单模型、垂直断层模型、正断层模型和多层介质模型）进行模拟计算表明，在模拟剖面上清晰可见直达波、反射纵波、反射横波、转换波和多次波能量，与理论分析一致．

斜向随机波对直立堤作用的试验研究

通过斜向随机波对直立堤作用的试验,研究了作用在直立堤上的单位堤长波浪力随波浪入射角的变化.讨论了波浪力沿堤长的纵向分布,以及单元堤长对波浪力纵向折减的影响.给出了波浪力纵向分布系数和折减系数与相对堤长的关系式.

长波高低不平顺与车体垂向加速度关联模型研究

轨道不平顺是车辆振动的主要激励源,随着车辆运营速度的增加,引起车辆振动加剧的不利波长也在变化。本文利用谱分析方法对0号高速综合检测列车实测数据进行分析,得到高速条件下长波高低不平顺和车体垂向加速度在空间频域上的分布规律,分别采用非参数模型和ARX模型构建长波高低不平顺与车体垂向加速度之间的关联关系,并利用系统辨识技术对两种模型的传递函数进行估计。通过将两种关联模型分别与长波高低不平顺结合,实现车体垂向加速度的预测;并与实测车体垂向加速度进行相关性和相干性分析。结果表明,两种模型均能较好反映长波高低不平顺与车体垂向加速度之间的传递关系。通过比较,ARX模型比非参数模型更精确。

黄海海雾WRF数值模拟中垂直分辨率的敏感性研究

基于WRF模式探究不同垂直分辨率下模式对黄海海雾的模拟表现。将3种垂直分层方案(35η、44η与63η)和2种边界层方案(YSU、MYNN)组合,对10次海雾做模拟研究。统计分析了水平雾区与雾顶高度对垂直分辨率的敏感性,并利用1次典型个例剖析了雾顶长波辐射与雾体湍流的作用。统计结果显示,提高垂直分辨率能显著改进水平雾区的模拟效果,改善明显的个例在不同垂直分辨率试验下的雾顶高度差异也较大;YSU方案相比于MYNN方案更敏感,从35η层增至44η层其试验个例的击中率(POD)和公正预兆得分(ETS)的平均改进率分别提高了13.29%和10.22%。典型个例研究结果揭示,模式对雾顶长波辐射和雾体湍流作用的模拟程度强烈依赖垂直分辨率:粗垂直分辨率给出的湍流强度很弱,导致模拟失败;雾顶存在"云水含量增多→长波辐射增强→降温加大→云水含量增多"的正反馈过程,细垂直分辨率比粗垂直分辨率更容易维持与增强此反馈;只有细垂直分辨率才能模拟出雾顶长波辐射冷却产生的贯穿雾体直抵海面、强度不弱于近海面机械湍流的浮力湍流,它导致雾体降温而出现符合观测事实的海温高于气温的现象。

斜压基本气流对东亚夏季风区气旋扰动低频发展影响的数值模拟研究

本文将夏季气候平均的基本气流分解为正压和斜压分量,使用一个线性斜压模式,研究了不同斜压基本气流对热带西北太平洋地区初始气旋性环流扰动低频发展演变的重要作用.其中,控制试验较好地模拟出初始气旋扰动向西北方向传播、在西北太平洋季风槽附近停滞增强、在东亚地区出现经向波列和在南海到海洋大陆地区形成西北—东南向波列等特征.改变斜压分量的敏感性试验结果表明,正压基流不能为西传的初始扰动供给足够的能量;海陆热力差异引起东亚地区的纬向温度梯度和北风垂直切变,是东亚太平洋型经向波列形成和维持的重要因素;当基本气流中的斜压纬向偏差部分线性增大时,扰动的能量会呈e指数迅速增强,提示在气候变化的背景下,基本气流微小的改变可能带来天气或季节内扰动强度的剧烈响应.