混合禁忌搜索算法在湖泊富营养化评价中的应用

本文提出了幂函数加和型湖泊富营养化评价指数公式,同时将禁忌搜索算法与模拟退火、免疫进化算法有机地结合起来,建立了求解优化问题的混合禁忌搜索算法(HTS).将此算法用于评价指数公式的参数优化,从而建立起湖泊富营养化综合评价模型,为湖泊富营养化评价和预测提供了一条有效的途径.将优化好的湖泊富营养化模型用于东湖富营养化评价的实例研究,结果表明此评价模型能较好地反映湖泊富营养化水平,并具有较好的通用性和实用性.     

基于QTM的海平面上升分析与模拟

针对海平面上升影响范围分析与模拟这一国际前沿问题,通过研究基于球面四元三角网(QTM)的关键技术问题,包括复杂拓扑关系计算、LOD剖分、球面水淹分析、基于QTM的多分辨率的DEM数据组织方法和分析精度的相关性评定等,以.Net和Direct3D为开发工具,设计开发了基于QTM的海平面上升影响范围评估模型。该研究结果可为全球海平面上升影响的防灾减灾决策提供有效支撑,并推动了球面数据模型和球面格网拓扑分析的理论成果在全球变化预测相关领域的应用进展。     

区域Delaunay三角剖分法在全国平均降水量中的应用

介绍了Delaunay三角网的性质及其算法类型;根据1980～2009年全国2200个观测站的降水量资料,将观测点和采集的边界点共同进行普通的Delaunay三角剖分,通过删除边界点及其区域外的三角形以实现区域Delaunay三角剖分,得到了较理想全国陆地的Delaunay三角网;随后对球面上的三角片进行面积计算,在已知站点的经纬度情况下,将大地坐标系转换到空间直角坐标系中,应用平面三角余弦定理获得球面三角内角,从而求得三角片面积,并以面积大小确定各个站点降水量的权重系数,得到全国平均降水量值。对比分析了30年的全国不同时间尺度(日、月、年)平均降水量,Delaunay三角法对应全国平均降水量均值和标准差都明显低于算术平均法,但是两种方法计算的降水量值的相关系数较高;通过Shapiro-Wilk方法进行正态性检验分析,两种计算方法求得的年平均降水量总体服从正态分布;在方差奇性的F检验中,两者的方差具有非奇性特点;使用t检验,在显著性α=0.05时,Delaunay三角剖分法计算的全国平均降水量总体均值偏小。最后,根据欧洲和日本数值模式2009年的降水预报,对于两方法计算结果进行了比较,分析表明在较大区域的平均降水量计算中,较之于传统的算术平均法,基于区域的Delaunay三角剖分法更为合理;区域平均降水量不仅和计算方法有关,还和区域气候特点有密切关系。     

GPU加速的改进PAM聚类算法研究与应用

空间聚类是空间数据挖掘的重要方法,而K-Medoids是一种常用的空间聚类算法。K-Medoids聚类算法存在初始点选择问题,而且计算复杂。为了提高算法的有效性和时间效率,本文结合模拟退火算法思想,改进了传统的K-Medoids算法PAM,提出一种基于GPU计算的并行模拟退火PAM算法。类比矩阵乘法运算,定义了一种新的矩阵计算方法,可以有效减少数据在GPU全局内存和共享内存之间的传输,提高了算法在GPU中的执行效率。利用模拟退火算法搜索聚类中心点,保证了聚类结果的全局最优性。基于不同的数据集,将串行和并行模拟退火PAM算法以及已有的遗传PAM算法进行比较,结果表明并行模拟退火PAM算法聚类结果正确,且时间效率高。最后,应用本文改进算法对贵州省安监系统的安全监管隐患数据进行聚类分析,发现了隐患聚集中心,相关结果对政府的决策具有一定的实际应用价值。     

Mesh grading approach for wave propagation in high velocity-contrast media

A mesh grading approach based on investigated lump method has been presented for simulating wave propagation in high velocity-contrast media. Unstructured fine grids are used to discretize the low wave-velocity medium in order to ensure the accuracy of numerical computation, and unstructured coarse grids are used for the high wave-velocity medium in order to substantially reduce the computational cost. On the interface, one coarse grid can match the fine grids of arbitrary odd number. The key feature of the proposed method is the constructions of investigated lumps on the interfaces of media. The transition zone, which is commonly used in the discontinuous grid scheme based on the staggered-grid finite-difference method, will not be used any more. Moreover, the computational instability that the discontinuous grid schemes frequently encountered does not arise in the proposed method. The comparisons with the analytical solutions and the application in studying the effects of sedimentary basin demonstrated that the mesh grading approach is a valid, accurate, convenient and flexible algorithm in simulating wave propagations in high velocity-contrast media with irregular interfaces.     

Turbulence characteristics of flows passing through a tetrahedron frame in a smooth open-channel

The turbulence characteristics of flows passing through a tetrahedron frame were investigated by using a 2-dimensional fiber-optic laser Doppler velocimeter (2-D FLDV). Experiments for uniform flows with different bed slopes under both submerged and un-submerged conditions were carried out in a re-circulating flume with glass side walls. The experimental bed was a smooth fixed bed. It was observed that with the tetrahedron frame the mean longitudinal velocity decrease in the retardation zone. However, both the longitudinal and the vertical turbulence intensities are larger than those for the undisturbed approach flow. The tetrahedron frame may reduce the probability of sediment entrainment by retarding the flow and reducing the boundary shear stress. In addition, it may induce sediment deposition in a sediment laden flow by changing the flow direction and increasing the energy dissipation.     

网格剖分及其精度和计算量分析

网格剖分程度直接影响着地震波正演数值模拟的计算精度及其计算量.以均匀介质模型为例,分析不同网格大小对波场模拟精度和计算量的影响,得出精细化网格剖分是实现高精度地震波正演模拟的有效方法,然而其计算量较大.以均匀倾斜介质模型为例,探讨倾斜地层网格剖分问题,数值实例分析不同震源频率对不同网格剖分方案引起的波场传播精度的影响,...     

On the determination of the mesh size for numerical simulations of shock wave propagation in near field underwater explosion

It is well known that the accuracy of mesh-based numerical simulations of underwater explosion strongly relies on the mesh size adopted in the analyses. Although a numerical analysis of underwater explosion can be performed with enough accuracy by using considerably fine meshes, such fine meshes may lead to substantially increase in the CPU time and the usage of computer memory. Thus, how to determine a suitable mesh size in numerical simulations is always a problem confronted when attempting to study the shock wave propagation resulting from underwater explosion and the subsequent response of structures. Considering that there is currently no universally accepted method for resolving this problem, this paper aims to propose a simple method to determine the mesh size for numerical simulations of near field underwater explosion. To this end, the mesh size effects on the shock wave propagation of underwater explosion are carefully investigated for different charge weights, through which the correlation between mesh sizes and charge weights is identified. Based on the numerical study, a dimensionless variable (λ), defined as the ratio of the radius of charge to the side length of element, is introduced to be the criterion for determining the mesh size in simulations. It is interesting to note that the presented method is suitable for various charge weights. By using the proposed meshing rule, adequate balance between solution accuracy and computational efficiency can be achieved for different blast scenarios in numerical simulations of underwater explosion.