大沽河地下咸淡水过渡带的水-岩作用模拟

大沽河是青岛市最大的地下水源地,过量开采导致严重的地下咸水入侵,研究地下咸-淡水的地球化学作用对于地下咸水的治理具有重要的指导作用。本文选取大沽河地下咸淡水过渡带作为研究区,在地下水动态监测的基础上,通过地下水水质和含水介质矿物分析,借助当前成熟的水文地球化学模拟软件—PHREEQC 2.0,计算地下水中主要离子和络合离子的含量,定量确定矿物的饱和指数和地下水流向上矿物相与水之间的物质交换量。研究结果表明:(1)地下水中钠、钾、镁、钙主要以单一离子形态存在,硫、碳主要以络阴离子形态存在,而硅和铝主要以络合物形态存在;(2)在地下水流动路径上,方解石、白云石、伊利石、高岭土处于过饱和状态;而石膏、CO2、氯化钠、萤石、玉髓处于不饱和状态;(3)咸淡水过渡带地下水化学类型由HCO3-Ca·Na型水向Cl-Na型水转变的趋势。     

迎风有限元法在三维潮流数值模拟中的应用

为研究和开发迎风有限元法在河口海岸三维波数值模拟中的应用，并为建立三维物质扩散模型，最终为解决实际工程问题打下基础，以非线性浅水波方程为基础，采用一咱平面迎风有奶元与垂向工有限差分相结合地数值计算方法建立三维潮流数学模型，基地实际的物理过程，在计算中将三维流动分成外重力波和内重力波，耦联求解潮位和流速的空间分布。应用本模式计算了天津新港附近渤海海域的三维潮流运动，以较高的分辩率揭示了潮流空间结构特     

深海环境中溢油输移扩散的初步数值模拟

通过对溢油在深海环境中的输移过程及行为特点的分析,初步建立基于拉格朗日积分法的深海溢油模型.该模型除了能够模拟油气混合物在真实深海环境中的共同输移与分离输移扩散过程,还考虑了石油溶解、气体溶解、天然气水合物形成与分解等行为变化对溢油运动轨迹的影响.应用该模型初步数值模拟了一次实际深海溢油试验,结果表明溢油在水下的空间分...     

辐射应力对台风风暴潮预报的影响和数值研究

台风过程期间,风暴潮和海浪是相伴相生的,相互作用的.波致辐射应力对于近岸风暴增、减水起着十分重要的作用,传统的海浪模式计算辐射应力耗时较多,不能满足业务化预报的要求.根据已有波浪辐射应力的理论表达式,经过严密的数学推导,适当的简化处理,提出了一个较为简单的波浪辐射应力表达式,并将其应用到业务化风暴潮数值预报模式中去,通...     

满足给定数字特征和谱的随机信号的产生

本文给出了满足给定数字特征和谱的随机序列的产生方法，即用无记忆非线性变换实现非正态随机信号的模拟．该法与Ｇｕｊａｒ法相比较，放宽了变量概率分布的限制，只要求变量矩的存在。同时考虑无记忆非线性变换对谱的影响，用异于Ｌｉｕ的方法在频域上求解输入正态分布随机序列的谱．     

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.     

Green water overtopping analyzed with a SPH model

Wave overtopping on the decks of offshore platforms and ships can cause severe damage due to the high forces generated by the water. This phenomenon is analyzed within the framework of the Smoothed Particle Hydrodynamics (SPH) method. The presence of a fixed horizontal deck above the mean water level modifies strongly the wave kinematics. In particular, the flow in the wave crest is split into two, showing a different behavior above and below the deck. Numerical results generated by the SPH method are compared to laboratory experiments. The formation of a jet in the rear of the deck after overtopping is observed under extreme conditions.     

New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation. A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.     

一种更有效的测量短波的方法

使用传统的在海底（或海底以上几米处）固定压力测波仪的方法观测波浪,因水层的滤波作用而无法有效地观测到高频短波。本文提出一种新的观测方法－使用悬挂压力测波仪的方法观测波浪。通过数值方法模拟与实际比较吻合的波面。然后模拟固定压力测波仪和悬挂坟力测波仪两种不同的方法获取此波面的水下压力记录曲线。应用传统谱分析的方法得到各自的压力频谱,再根据不同的情况分别换算成表面谱。通过数值的方法可以证明使用悬挂压力测     

随机激励下非线性海洋结构物响应分析方法研究进展

对在随机载荷作用下非线性海洋结构物的运动响应分析的各种方法进行了综述,对有些方法的基本原理做了阐述,并指出一些最新的研究进展和今后进一步研究的方向。