A Compact Explicit Difference Scheme of High Accuracy for Extended Boussinesq Equations

Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan.     

海岸工程中海浪数值模拟和环境参数选取的进展

主要介绍了我国近的来在海岸工程设计，建设中有关海浪数值模拟及环境参数选取方面的研究进展。     

水面皮层破坏-复原过程的实验室模拟研究

依据水面红外发射和红外遥感测温原理，采用HDG－高灵敏度红外测温仪和常规测量仪器相结合的方法，在实验室空气稳定条件下，模拟测得了水面皮层破坏－复原（重建）的热力过程和气－水温差对水面皮层复原过程的影响，获得了大量的测量数据。数据分析表明，当气－水温差从3．0℃变为11．5℃时，水面皮层破坏可导致皮温增量从气－水温差3．0℃时的0．15℃变到11．5℃时的0．45℃，并发现恢复时间与气－水温差呈负线性关系。     

闽江口水动力和污染扩散数值模拟

本文建立了闽江口二维水动力、污染扩散数值模型，模型采用ＡＤＩ双向隐式求解。在水动力模拟中，对岸界的丁坝有用线性边界进行拟合，并考虑到河口区侧向摩擦作用；在污染扩散模拟中，用动态拟合的方法对浓度水边界进行了较为合理的改进。计算结果表明，本模型是成功的。     

渤、黄、东海内潮的数值模拟

在全球的海洋中,中国东海和临近海域是最显著的内潮生成地之一。本文采用NODC(Levitus) World Ocean Atlas 1998提供的季平均温、盐资料,计算海水的密度,并计算垂向密度梯度的最大值点,得到一个较符合海水实际的密度分层。使用三维非线性数值模型(将海洋分为2层)研究了潮汐(M2,S2,K1,O1分潮)作用下渤黄东海的内潮,揭示了整个海区内潮起伏的空间分布,结果发现大振幅的波动均发生在台湾东北(冲绳海槽)海域和中国近海地形突变之处,其中前者更显著。对于各分潮模拟得到的表面潮与TOPEX/Poseidon高度计资料基本一致。研究结果表明上层海水的深度和厚度的梯度对内潮有一定的影响;冬季分布区域比夏季小,强度比夏季大。     

协方差阵二次型容许估计的一个必要条件

研究协方差阵Σ的二次型容许估计问题。设 y1,y2 ,… ,yniid,n≥ 2 ,y1与 p维正态分布N (β,Σ )有相同的前四阶矩。其中β =(β1,β2 ,… ,βp)′∈ Rp与Σ =(σij) p× p >0均未知。记 y =△ (y1,y2 ,… ,yn)′。在二次损失 L (d ,Σ ) =tr(d -Σ) 2下给出Σ的二次型估计 a S2 + nby-y-′是容许估计的必要条件为 :(n - 1) a + b + 2 max(a,b)≤ 1。此必要条件比张立振等协方差阵的二次型容许估计中的必要条件有了明显的加强     

湍流数值模拟中封闭模式应用的局限性

通过对琼州海峡的潮流场特征进行数值模拟，指出了选择不同的特征混合长度表达式对数值模拟结果的影响，表明了基于特征混合长度理论的流封闭模式在近海湍流数值模拟中应用的局限性。     

1994年发生在台湾海峡的一次地震海啸的数值模拟

建立了一个地震海啸数值模式,模式包含越洋海啸传播部分和近岸海啸变形部分,在越洋海啸传播部分中采用线性浅水方程,使用蛙跃格式求解,并且选择合适的空间步长与时间步长,使差分格式中产生的数值频散与包辛尼斯克方程中的物理频散一致,这样在不影响海啸数值计算精度的前提下,节省了计算机的机时与内存.在近岸海啸变形部分的计算中,考虑了非线性对流项与海底摩擦项.同时该模式采用了多重网格嵌套技术,提高了所关心地区的计算精度.利用这个地震海啸模式模拟了1994年发生在台湾海峡的一次地震海啸,结果与观测记录较吻合.这个模型已用于我国沿海核电站可能最大地震海啸的数值计算.