礁坪上波浪传播的数值模拟

本文基于具备间断捕捉能力的二阶全非线性Boussinesq数值模型,对规则波和随机波在礁坪地形上的传播变形进行了数值模拟。该模型采用高阶有限体积法和有限差分方法求解守恒格式的控制方程,将波浪破碎视为间断,同时采用静态重构技术处理了海岸动边界问题。重点针对礁坪上波浪传播过程中的波高空间分布和沿程衰减,礁坪上的平均水位变化,以及波浪能量频谱的移动和空间差异等典型水动力现象开展数值计算。将数值结果与实验结果对比,两者吻合情况良好,验证了模型具有良好的稳定性,具备模拟破碎波浪和海-岸动边界的能力,能较为准确地模拟波浪在礁坪地形上的传播过程中发生的各种水动力现象。     

Modeling of iceberg drift in the marginal ice zone of the Barents Sea

The paper is focused on the analysis of the drift of tabular iceberg observed in 2009 in the marginal ice zone of the North-West Barents Sea. Momentum balance equations are derived from the Kirchhoff equations describing plane motion of solid body in an ideal fluid. Field works performed on the drift iceberg and on the drift icenear the iceberg are described. Results of the field works and numerical simulations of the iceberg drift and rotation are performed and discussed. It is shown that acceleration of water flow around the iceberg has visible influence on the iceberg drift. Kinetic energy balance of drift iceberg is used to estimate the forces applied to the iceberg by the drift ice.     

URANS investigation of the interaction between the free surface and a shallowly submerged underwater vehicle at steady drift

An axisymmetric underwater vehicle (UV) at a steady drift angle experiences the complex three-dimensional crossflow separation. This separation arises from the unfavorable circumferential pressure gradient developed from the windward side toward the leeward side. As is well known, the separated flow in the leeward side gives rise to the formation of a pair of vortices, which affects considerably the forces and moments acting on the UV. In this regard, the main purpose of the present study is to evaluate the role of the leeward vortical flow structure in the hydrodynamic behavior of a shallowly submerged UV at a moderate drift angle traveling beneath the free surface. Accordingly, the static drift tests are performed on the SUBOFF UV model using URANS equations coupled with a Reynolds stress turbulence model. The simulations are carried out in the commercial code STARCCM+ at a constant advance velocity based on Froude number equal to F_n = 0.512 over submergence depths and drift angles ranging from h = 1.1D to h = ∞ and from β = 0 to β = 18.11°, respectively. The validation of the numerical model is partially conducted by using the existing experimental data of the forces and moment acting on the totally submerged bare hull model. Significant interaction between the low-pressure region created by the leeward vortical flow structure and the free surface is observed. As a result of this interaction, the leeward vortical flow structure appears to be largely responsible for the behavior of the forces and moments exerted on a shallowly submerged UV at steady drift.     

基于GPU加速的Boussinesq类波浪传播变形数值模型

Boussinesq波浪模型是一类相位解析模型,在时域内求解需要较高的空间和时间分辨率以保证计算精度。为提高计算效率,有必要针对该类模型开展并行算法的研究。与传统的中央处理器(CPU)相比,图形处理器(GPU)有大量的运算器,可显著提高计算效率。基于统一计算设备架构CUDA C语言和图形处理器,实现了Boussinesq模型的并行运算。将本模型的计算结果同CPU数值模拟结果和解析解相比较,发现得到的结果基本一致。同时也比较了CPU端与GPU端的计算效率,结果表明,GPU数值模型的计算效率有明显提升,并且伴随数值网格的增多,提升效果更为明显。     

最小二乘定权方法对流动重力计算结果的影响

在处理流动重力观测数据时,采用最小二乘方法解算由权系数和观测值组成的超定方程组,得到最优权值,并利用最优权计算重力值。结果表明,等权方法计算得到的结果点值精度不高,但包含了较全的重力异常地震分析预报信息;不等权方法计算得到的结果点值精度较高,但削弱了有用的重力异常地震分析预报信息;最小二乘方法计算得到的结果不但有较高的点值精度,对M_S4.0及以上地震的分析预报效果也最优。     

扩展共线方程结合车载组合全景相机影像区域网平差

利用扩展共线方程实现组合相机单片影像的空中三角测量,通过该方法来提高全景影像的定位精度。由于组合全景相机中各子相机刚性固定,相机间的相对位姿参数可以精确标定,因此利用扩展共线方程模型,可以将组合全景相机的构像单元从单个子相机影像扩展为整个组合相机获取的多张影像。与已有的球模型方法,包括球理想模型和球严密模型方法相比,本文所用的扩展共线方程方法使用各子相机影像像点作为观测值,避免了球模型方法中由各子相机摄影中心不重合导致的球模型投影误差和摄站内重叠区域融合误差。试验结果表明:与球模型方法相比,本文方法无需模型投影和摄站内重叠区域融合等繁琐的处理过程,避免了由其引起的精度损失,保证了空三网稳定性的同时提高了定位精度。     

电力系统元件非线性微分-代数子系统模型的逆系统控制:一种新算法

针对电力系统元件非线性微分-代数子系统模型,本文提出一种新算法研究其逆系统控制问题.所提出的新算法不需要对控制输出及其高阶导数做复杂的变换,具有更好的应用性.本文的逆系统控制方法主要分为两步：第1步,利用所提出的新算法来判断被控元件的可逆性,若可逆,则基于状态反馈与动态补偿,构造出元件的α阶积分右逆系统,实现复合系统的线性化和解耦;第2步,利用线性控制的理论和方法设计闭环控制器,使得元件被控对象满足期望的性能指标.最后按照本文所提出的方法,研究了多机电力系统的分散非线性汽门控制问题.仿真结果验证了本文所提方法的有效性.     

斜压海洋动力学的一种三维数值模式──Ⅱ.温度、盐度和垂直涡动粘性系数的计算

根据海水温度和盐度平流扩散方程给出一种数值计算方案并采用混合长度理论给出垂直涡动粘性系数的计算方法。对于温度和盐度方程,其平流过程采用了Lax－Wendroff格式,水平扩散采用显格式,垂直扩散采用隐格式。时间步长主要受平流过程的Courant-Friedrichs-Lewy条件限制。垂直涡动粘性系数计算依据Prandtl混合长度理论,并考虑了海水层化的抑制作用,因而其数值与流场及密度场结构有关。温度、盐度及垂直涡动粘性系数的计算与动力方程中内模态的计算同步进行。应用本模式模拟渤、黄、东海由潮流、密度流和风海流迭加而成的综合海流,得出了良好的结果。     

官厅水库水质模型研究

根据官厅水库狭长、水浅和水文、水力、水质迁移扩散等特点与河流相似的事实,采用一维St.Venant方程组模拟其水流状况;采用完全混合反应器概念模拟其水质变化,两者的藕合就构成了官厅水库的水质模型。该模型在模拟官厅水库NH₃-N浓度时空变化时,精度令人满意,并已用于官厅水库的水质预测和水质管理。     

Refraction Effect in Geothermal Heat Flow Due to a 3-D Prism in a Two-Layered Earth

A mathematical model of the disturbance of the stationary geothermal field due to a three-dimensional perturbing body embedded near the surface of a two layered earth. The theoretical analysis is based on the generalized theory of the double-layer potential, similar to the boundary integral method likely to direct current geoelectricity problems. Special attention is paid to cases in which a pan of the body's surface is tangential to the planar boundaries of the substratum or of the surface of the earth The numerical calculations are performed for a three-dimensional prism using methods similar to the boundary element method (BEM). Numerous graphs are shown for the disturbance of the temperature field and of the heat flow on the surface of the earth and along vertical profiles (boreholes).