基于BP神经网络的天然气水合物相平衡计算及预测

针对天然气水合物相平衡问题,文中提出用基于带动量因子的BP神经网络进行计算和预测。首先用遗传算法优化确定BP神经网络的结构和参数,得到最优化结构的神经网络;其次结合Levenberg-Marquart优化算法,建立天然气水合物相平衡计算及预测的神经网络模型;最后以实验测定的(CH4 CO2 H2S)三元酸性天然气水合物体系的平衡数据为训练和预测样本进行了计算。计算表明,预测结果与实验数据有良好的一致性,而且由于BP神经网络作为所谓的“纯粹”的算法不需要热力学模型,这对于相平衡计算是非常方便的,所以是研究天然气水合物相平衡计算及预测的一种新的有效方法。     

拖缆动力学离散化计算的后处理方法研究

拖缆的计算一般都是将拖缆离散化,建立非线性的动力学平衡方程组,通过龙格库塔法、Newton-Raphson迭代法等方法对非线性方程组进行求解.本方法是先将拖缆离散成若干段,列出离散段节点处的水动力学方程,并进一步建立节点处关于角加速度的线性平衡方程组,通过对角加速度的求解计算来解决拖揽的动力学问题.在计算过程中,将拖揽节点处关于角度的非线性问题转化为角加速度的线性问题,简化了求解非线性方程组的计算过程,避免了非线性方程组迭代求解中的初始点问题.     

青海一里平盐湖卤水Li-Na-K-Mg-Cl-SO_4-H_2O体系25℃时渗透系数、活度系数和饱和度的预测

将Ritrer离子作用模型的应用扩展到中国富含锂的一里平盐湖卤水Li－Na-K-Mg－Cl-SO4－H2O体系,综合1990－1997年作者对该体系的含锂三元次级体系25℃时的活度系数和渗透系数的研究结果和文献数据,给出了该体系的Pitzer离子作用模型参数;计算了溶解平衡常数;预测了青海一里平盐湖晶间卤水中离子和盐的活度系数、卤水的渗透系数、水活度、蒸汽压以及盐类的饱和度;揭示了卤水在蒸发和稀释过程中上述热力学性质随离子强度变化的规律,用离子强度与盐类饱和度的关联图（I-ln（K’／K）图）分析了该卤水的折盐顺序。计算结果表明,一旦平盐湖晶间因水在25℃时只有石盐饱和,与野外观测结果相一致;用参数化的Pitzer离子作用模型预测该卤水多组分体系的热力学性质和溶解平衡是可行的。     

近岸波浪在刚性植被区域传播的数值模型

基于扩展型Boussinesq水波方程,建立了波浪在刚性植被覆盖的近岸海域传播的数值模型。通过在动量方程源项中引入拖曳阻力项考虑植被对波浪的衰减作用。控制方程采用有限差分和有限体积混合格式求解,模型稳定性强,具备间断捕捉能力,能有效模拟近岸区域波浪的传播变形、破碎和处理海岸动边界问题。利用所建立模型对典型物理模型实验进行模拟,计算结果与实验结果吻合良好,表明模型可用于波浪在刚性植被覆盖海域的数值计算。     

简易海洋平台疲劳可靠性优化设计的进化算法

分别用进化计算中的遗传算法和进化规划进行了简易海洋平台基于非线性疲劳损伤模型的疲劳可靠性优化设计.计算结果表明了计算模型的可靠性以及遗传算法和进化规划解决结构疲劳可靠性优化设计的有效性,并对两种优化算法进行了比较和讨论。     

二维透声窗遮挡损失计算方法研究

为了水下探测和声对抗,舰艇艏部舱室通常安装主动声呐设备,声呐发射声波经由透声窗向外传播,因此,开展透声窗声学特性的相关研究具有重大现实意义。为解决实际工程中透声窗透声性能计算较困难这一问题,一种用于快速预报透声窗的遮挡损失的解析–数值混合方法被提出。该方法用球谐波展开法和叠加原理求解阵列活塞声源声场,得到声激励下透声窗模型表面的声压和法向振速,进而利用COMSOL的Kirchhoff-Helmholtz模块通过再辐射的方法求解模型的遮挡损失。最后,建立了一个二维透声窗模型,分别采用解析–数值混合方法和纯有限元法计算了其全方位的遮挡损失,验证了该方法的准确性、有效性。     

对南海主要出流通道的诊断及讨论

通过一个基于质量、热量和盐量守恒的最优化盒子模型,探讨南海的主要出流通道,并详细阐述了该盒子模型中最优化的物理含义,指出这是一种平流排热、排盐最高效的系统.计算结果表明,卡里马塔海峡是南海最主要的出流通道;南海海表有巨大的热通量和淡水通量输入,作为一个热盐守恒的系统,将高温低盐的海水从宽阔的卡里马塔海峡排出是维持整个海区热盐平衡最快速有效的途径.应用更为复杂的、考虑了热盐结构和动力学因素的逆方法进一步地验证了以上结果,同时逆方法结果还揭示了吕宋海峡断面流速分布的一些特征.     

港池中波浪折射、绕射的数值模拟方法

从流体运动方程和动量方程出发,引入海底摩擦力和港池不完全反射作用,推导出了Berkhoff折射、绕射方程。在深海部分解析方法用精确解(含待定系数)、在复杂地形用数值离散方法求解,中间过度段用的光滑匹配将离散数值(有限元)和解析解(精确解)一同求解。通过极值问题建立泛函,利用泛函的驻定性将海岸(港湾)问题进行数值离散,建立了可行的数值模拟模型。     

KCl+KHCO_3水溶液中氯化钾的平均活度系数和KHCO_3~0的缔合常数

用钾离子选择电极和银—氯化银电极构成的无液界电池测量了KCl+KHCO_3混合水溶液中KCl的平均活度系数,发现其活度系数的对数与混合溶液中KHCO_3的浓度m_2服从哈奈特关系式。用最小二乘法计算机处理求得了各离子强度条件下相应的哈奈特系数。还用适用于混合电解质溶液的皮茨公式计算了该体系电解质的平均活度系数和溶液的超额自由能,计算出活度系数与实验值之间的平均标准偏差为±0.006。溶液的超额自由能随溶液组成变化表明,对于有相同阳离子(K~+)的该混合溶液中离子间短程作用力主要受Cl~-作用的影响。用累代逼近法用微机计算,求得了KHCO_3~0的化学计量缔合常数经验公式及其热力学缔合常数。     

基于工程反问题的圆管湍流流体热物性参数反演

给出一种利用测定圆管对流换热温度场求解湍流流体热物性参数的方法。根据边界层理论,对管内湍流强制对流换热控制方程组进行简化,并通过考虑湍流运输中的密度脉动,对传统的混合长度湍流模型进行改进,基于该计算模型求解圆管湍流模型的温度场。采用Levenberg-Marquardt方法通过测量温度场反演了湍流对流换热控制方程中的黏度系数。结果表明,利用反问题方法通过测定温度场反演湍流流动流体的热物性参数具有可行性,精确度较高。     

Optimization Design of Fairings for VIV Suppression Based on Data-Driven Models and Genetic Algorithm

Vortex induced vibration(VIV)is a challenge in ocean engineering.Several devices including fairings have been designed to suppress VIV.However,how to optimize the design of suppression devices is still a problem to be solved.In this paper,an optimization design methodology is presented based on data-driven models and genetic algorithm(GA).Data-driven models are introduced to substitute complex physics-based equations.GA is used to rapidly search for the optimal suppression device from all possible solutions.Taking fairings as example,VIV response database for different fairings is established based on parameterized models in which model sections of fairings are controlled by several control points and Bezier curves.Then a data-driven model,which can predict the VIV response of fairings with different sections accurately and efficiently,is trained through BP neural network.Finally,a comprehensive optimization method and process is proposed based on GA and the data-driven model.The proposed method is demonstrated by its application to a case.It turns out that the proposed method can perform the optimization design of fairings effectively.VIV can be reduced obviously through the optimization design.     

不同形态海湾的水底平衡剖面研究

研究海湾平衡剖面对理解海湾地貌演变具有重要意义。本文给出了收缩型、扩张型和矩形3种典型海湾平面形态对平衡剖面的影响。建立了海湾长度远小于潮汐波长的短尺度海湾的平衡剖面和对应的时均悬沙浓度的解析解。采用水深平均的水动力方程、泥沙输移方程和地形演变方程的耦合模型对以上3种类型海湾的平衡剖面进行了数值模拟,得到了这3种类型海湾的水面、流速、时均悬沙浓度和平衡剖面的计算结果,并利用水面数值结果确定了海湾水面解析解所含的一个待定常数。研究结果给出了3种不同海湾平面形态所对应的平衡剖面形态:矩形海湾对应斜坡型;收缩型海湾对应下凹型;扩张型海湾对应上凸型。所得海湾平衡剖面和时均悬沙浓度的解析解与数值解一致。     

潮汐汊道均衡断面面积计算与稳定性分析中的问题 ——以海南陵水新村为例

通过对潮汐汊道稳定性与均衡断面面积定义的探讨,明确了O Brien等、Kreeke与Bruun对潮汐汊道稳定性定义的差异。通过不同方法求解一维水力学与连续方程,计算了海南陵水新村潮汐汊道的封闭曲线,并采用华南沿海A-P关系得出的均衡流速曲线及Soulsby临界流速曲线,计算了新村潮汐汊道的均衡断面面积,根据O Brien等与Kreeke的稳定性定义对新村潮汐汊道的稳定性进行了判别。结果表明,在大、小潮期用不同方法及不同的断面变化方式得出的稳定均衡断面面积(A2)在同样的均衡流速曲线下基本一致,而小潮期的稳定均衡断面面积要小于大潮期,最大流速处的断面面积(Acc)要大于大潮期。A2、Acc及饱满系数为0.8时对应的断面面积三者之间具有较大的差异。用封闭曲线求取潮汐汊道均衡断面面积并进行稳定性判断,在概念上有一定意义,可以作为潮汐汊道稳定性判断的一个依据,但通过求解一维水力学与连续方程以获取均衡断面面积则具有较大的局限性与不确定性。应通过中长期的野外观测、数值模拟与物理模型,以确定潮汐汊道的水动力特征以及地形演变规律,得出A-P关系进而求出均衡断面面积,并结合Bruun方法判定潮汐汊道的稳定性。     

砂质底质潮汐水道均衡态模拟初探

以过程模拟法探讨了发育在砂质堆积体上的潮汐水道的均衡态问题.在潮汐作用为主导、初始状况为平坦海底、底质粒度均一等条件下,利用修正的Bagnold推移质输沙率公式和沉积物连续方程,建立砂质潮汐水道均衡态的一维模型.数值实验结果表明,潮汐水道达到均衡态的时间尺度为10a以上,水道形成的早期冲刷速度较快,随着时间的推移,冲刷减缓.从水道冲刷出来的物质在水道前方形成堆积体.潮汐水道的演变受控于水动力变化和沉积物冲淤的负反馈机制,水道均衡态主要与沉积物粒径、涨落潮流速、历时、流速沿水道轴向分布的组合特征有关.流场特征决定了水道均衡态的形态;在沉积物组分均一的情况下水道均衡深度与潮流流速呈正相关;达到均衡态的时间尺度与沉积物粒径、涨落潮历时、流速的变化以及水道规模有关.沉积物临界起动流速随深度的变化对水道均衡态有较大影响.     

An Engineering Application Tool for Visual Assessment of the Equilibrium Beach

An engineering application tool for prediction of the static equilibrium bay(Beach Mod)is established to describe two bay shape formulas by use of the programming software "MATLAB" with a graphic user interface(GUI).The tool is user-friendly for engineering students for the design of beach shapes.This tool was tested through application on three types of beaches in Taiwan and Australia.By implementing the concept of Headland Control,the Beach Mod program allows users to draw a structure and create an artificial headland.The results indicate that Beach Mod can efficiently forecast beach changes as well as MEPBAY,a competing software package,while boasting a better user interface.     

Towards a true 4‐dimensional data assimilation algorithm: application of a cycling representer algorithm to a simple transport problem

3‐dimensional variational algorithms are widely used for atmospheric data assimilation at the present time, particularly on the synoptic and global scales. However, mesoscale and convective scale phenomena are considerably more chaotic and intermittent and it is clear that true 4‐dimensional data assimilation algorithms will be required to properly analyze these phenomena. In its most general form, the data assimilation problem can be posed as the minimization of a 4‐dimensional cost function with the forecast model as a weak constraint. This is a much more difficult problem than the widely discussed 4DVAR algorithm where the model is a strong constraint. Bennett and collaborators have considered a method of solution to the weak constraint problem, based on representer theory. However, their method is not suitable for the numerical weather prediction problem, because it does not cycle in time. In this paper, the representer method is modified to permit cycling in time, in a manner which is entirely internally consistent. The method was applied to a simple 1‐dimensional constituent transport problem where the signal was sampled (perfectly and imperfectly) with various sparse observation network configurations. The cycling representer algorithm discussed here successfully extracted the signal from the noisy, sparse observations     

A hybrid simplex genetic algorithm for estimating geoacousticparameters using matched-field inversion

Matched-fieId inversion (MFI) undertakes to estimate the geometric and geoacoustic parameters in an ocean acoustic scenario by matching acoustic field data recorded at hydrophone array with numerical calculations of the field. The model which provides the best fit to the data is the estimate of the actual experimental scenario. MFI provides a comparatively inexpensive method for estimating ocean bottom parameters over an extensive area. The basic components of the inversion process are a sound propagation model and matching (minimization) algorithm. Since a typical MFI problem requires a large number of computationally intensive sound propagation calculations, both of these components have to be efficient. In this study, a hybrid inversion algorithm which uses a parabolic equation propagation model and combines the downhill simplex algorithm with genetic algorithms is introduced. The algorithm is demonstrated on synthetic range-dependent shallow-water data generated using the parabolic equation propagation model. The performance for estimating the model parameters is compared for realistic signal-to-noise ratios in the synthetic data     

Punishment in optional public goods games

In this work, the optional public goods games with punishment are studied. By adopting the approximate best response dynamics, a micro model is given to explain the evolutionary process. Simultaneously, the magnitude of rationality is also considered. Under the condition of bounded rationality which provides a light to interpret phenomena in human society, the model leads to two types of equilibriums. One is the equilibrium without punishers and the other is the equilibrium including only punishers and cooperators. In addition, the effects of rationality on equilibriums are briefly investigated.     

Numerical modeling of crenulate bay shapes

A numerical model is developed to compute the shoreline planform in a crenulate bay beach. The new model combines polar and Cartesian coordinates and can be used effectively to compute a hooked zone shoreline in the lee of upcoast headland. The model is calibrated using laboratory data with an incident wave angle ranging from 25° to 60°. The results of calibration and verification suggest that the ratio of the sediment transport parameters by wave and longshore current in this model is close to unity, and the computed shoreline planforms for the hooked and unhooked zones are in good agreement with the ones measured, especially when a bay is close to static equilibrium. In addition, the bay shape calculated by the present model is similar to that given by the well-known empirical parabolic equation for a bay in static equilibrium. The process of bay shape development from a straight beach to a static equilibrium bay is studied using laboratory experiments and the present numerical model. The temporal variations in the computed longshore sediment transport at different locations within a bay beach are analyzed. From this the decrease in the sediment transport becomes apparent while a bay beach changes its shape from straight toward a state of equilibrium. Based on this experience, it may be concluded that the present numerical model can produce a temporal change in the shoreline planform of a crenulate bay beach from a transition state to static equilibrium subject to seasonal wave action.