A time-dependent numerical model of the mild-slope equation

On the basis of the previous studies, the simplest hyperbolic mild-slope equation has been gained and the linear time-dependent numerical model for the water wave propagation has been established combined with different boundary conditions. Through computing the effective surface displacement and transforming into the real transient wave motion, related wave factors will be calculated. Compared with Lin’s model, analysis shows that calculation stability of the present model is enhanced efficiently, because the truncation errors of this model are only contributed by the dissipation terms, but those of Lin’s model are induced by the convection terms, dissipation terms and source terms. The tests show that the present model succeeds the merit in Lin’s model and the computational program is simpler, the computational time is shorter, and the computational stability is enhanced efficiently. The present model has the capability of simulating transient wave motion by correctly predicting at the speed of wave propagation, which is important for the real-time forecast of the arrival time of surface waves generated in the deep sea. The model is validated against analytical solution for wave diffraction and experimental data for combined wave refraction and diffraction over a submerged elliptic shoal on a slope. Good agreements are obtained. The model can be applied to the theory research an d engineering applications about the wave propagation in a biggish area.     

Model simulations of carbon sequestration in the northwest Pacific by patch fertilization

Iron fertilization of nutrient-rich surface waters of the ocean is one possible way to help slow the rising levels of atmospheric CO₂ by sequestering it in the oceans via biological carbon export. Here, I use an ocean general circulation model to simulate a patch of nutrient depletion in the subpolar northwest Pacific under various scenarios. Model results confirm that surface fertilization is an inefficient way to sequester carbon from the atmosphere (Gnanadesikan et al., 2003), since only about 20% of the exported carbon comes initially from the atmosphere. Fertilization reduces future production and thus CO₂ uptake by utilizing nutrients that would otherwise be available later. Effectively, this can be considered as leakage when compared to a control run. This “effective” leakage and the actual leakage of sequestered CO₂ cause a significant, rapid decrease in carbon retention (only 30–45% retained after 10 years and less than 20% after 50 years). This contrasts markedly with the almost 100% retention efficiency for the same duration using the same model, when carbon is disposed directly into the northwest Pacific (Matsumoto and Mignone, 2005). As a consequence, the economic effectiveness of patch fertilization is poor in two limiting cases of the future price path of carbon. Sequestered carbon in patch fertilization is lost to the atmosphere at increasingly remote places as time passes, which would make monitoring exceedingly difficult. If all organic carbon from one-time fertilization reached the ocean bottom and remineralized there, acidification would be about −0.05 pH unit with O₂ depletion about −20 μmol kg⁻¹. These anomalies are probably too small to seriously threaten deep sea biota, but they are underestimated in the model because of its large grid size. The results from this study offer little to advocate purposeful surface fertilization as a serious means to address the anthropogenic carbon problem.     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.     

波峰特征量的统计分布

本文在作者另文数值模拟得到的三维海浪基础上 ,进一步分析给出了三维波峰长度、高度及方向角度等特征量的统计分布。发现考虑海浪的方向性质后波峰分布不再是均匀连续的 ;波峰的长度分布受方向函数影响 ,方向分布越宽 ,波峰的平均长度越短 ,波峰的方向角分布越宽 ;波峰高度和长度在波峰高度较小时有很大的相关性 ,而在波峰高度很大时无关     

海南东寨港红树林港湾潮汐动力研究

在海南东寨港林市村附近红树林港湾潮汐动力测量的基础上，建立二维可移动边界有限元数值模型，研究红树林港湾浅水湖滩潮流场特征。实测与计算结果表明，潮沟区流速最大，具驻波特征，呈明显时间－速度不对称，东潮沟涨潮优势，百潮沟落潮优势，整个港湾水体呈顺时针循环。白滩与林地周期性淹没，白滩流速较湖沟略小。林地曼宁系数为白滩的１０倍，流速小于１０ｃｍ·ｓ（－１），约为相应潮沟流速的１／１０。林地如去掉林木，流速可增加３倍以上，流向也有变化。     

洋浦近岸海域污染控制规划方案优化——Ⅰ.排污口位置优选

在洋浦近岸海域计算潮流场基础上,计算预选排污口附近海水质点运动轨迹及预测污染物浓度分布,最后,从环保角度出发,推荐了洋浦地区各开发区的排污口位置,为洋浦近岸海域污染控制规划方案的制定及其优化提供了依据。     

油田高含水期基于微观渗流机理的宏观岩心水驱油数值模拟研究

建立地层孔隙网络模型描述微观孔喉特征,利用逾渗方法计算基于微观渗流机理的相对渗透率曲线,并将之应用到数值模拟计算中,给出考虑微观渗流特征的方程,并进行数值求解,得到地层中流体分布特征。实例分析表明该方法可以得到基于微观特征的剩余油分布数量和形式,为调剖堵水、化学驱提供技术参数,实现微观渗流机理与宏观计算的结合。     

环台湾岛海域半日潮波特征的三维模拟

用1997版POM海洋模式,首次应用于环台湾岛海域的潮波数值研究.得到该海域的半日潮波主要为23°N以南西太平洋传来的胁振潮.影响台湾海峡的半日潮波分别由海峡南北口传入的两支潮波,且北支强于南支.福建沿岸湄州湾-兴化湾为最强潮区,其M₂分潮最大振幅可达240cm.最强潮流区位于澎湖水道,M₂分潮最大潮流达196cm/s.环台湾岛海域潮波潮流水平结构上除海峡北部原有一个圆流点外,还发现另外存在4个新的圆流点.潮流垂直结构上主要为右偏,接近底层处为左偏.     

川西盆地雾和能见度的气候特征及其对飞行的影响

利用川西盆地广汉机场气象台1986年1月~1995年12月共10 a地面观测资料,对广汉机场低能见度的生消时间、持续时间和能见度的年季节变化及其日变化做了统计分析,研究广汉机场全年能见度的飞行气候特征,特别是低能见度气候特征对飞行的影响。结果表明,能见度日变化特征表现为:白天早、晚时段能见度较低,午后能见度达到最高。能见度的季节变化特征表现为:冬(12,1,2月)、春季(3,4,5月)能见度较低,而夏(6,7,8月)、秋季(9,10,11月)能见度较高。冬季能见度的变化主要受辐射雾影响,夏季能见度的变化主要受降水的影响,随着降水强度的变化而变化,低能见度时间短暂。各标准低能见度逐时频率的日变化特征在各季基本相似,低能见度高频率区均出现在早晨到上午,低频率区均出现在下午。     

南海东北部台风潮的数值后报及极值增水估计

本研究对近40年来影响南海东北部陆架海区的28次台风引起的风暴潮进行了数值后报,其中8个过程的沿岸后报增水值与实测值进行了比较,表明后报值与实测值符合良好,90%以上的最大增水值偏差在30厘米以下。为了得出本海区多年一遇的台风增水极值,在后报台风路径密集处选择了9个不同水深点,对每点取出各次台风下的最大增水值,然后用Weibull分布进行拟合,得出了各点的极值分布。