辅助微分方程完全匹配层在声波方程数值模拟中的应用

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。     

有限差分格式与物理量守恒性

为了成功地开展大气数值模拟,如何构造相应的离散格式是一项十分重要的工作.本文以大气球面浅水方程为例,从不同形式的微分方程组出发,分别构造了能量守恒格式、准辛格式和多守恒格式,并作相应的数值检验和比较,从中可以清楚地看到三类格式的优劣.这三类格式构造方法不同,所保持的物理守恒性不完全相同,数值计算时所需机时也不同,这些比较结果为实际工作的不同使用需求提供了选择依据.     

Numerical Modeling of Sediment Mass Sliding and Tsunami Generation: The Case of February 7, 1963, in Corinth Gulf,Greece

Numerical modeling of a landslide process at the continental slope by taking into account slope parameters is performed for the tsunami event of February 7, 1963, reported in the Corinth Gulf, central Greece. A layered sediment structure was considered, and an initial external dynamic action at the landslide process was introduced. The results obtained were quite consistent with the observational runup data as well as with the results arrived at by conventional rigid-body and viscous-fluid models. However, this approach permits to describe in detail the formation of several wave groups and particular tsunami characteristics that are strongly dependent on the landslide model.     

具有精确色散性的非线性波浪水平二维数学模型

建立具有色散性的水平二维非线性波浪方程,方程的非线性近似到了三阶。方程以波面升高和自由表面速度势表达的微分-积分型数学方程,给出方程的数值求解方法和算例,对方程积分项的处理给出了计算方法。计算结果与Boussinesq方程模型和缓坡方程模型的对应计算结果进行了对比。     

Wave Height Transformation and Set-up Between a A Submerged Permeable Breakwater and A Seawall

In this study, we investigated wave transformation and wave set-up between a submerged permeable breakwater and a seawall. Modified time-dependent mild-slope equations, which involve parameters of the porous medium, were used to calculate the wave height transformation and the mean water level change around a submerged breakwater. The numerical solution is verified with experimental data. The simulated results show that modulations of the wave profile and wave set-up are clearly observed between the submerged breakwater and the seawall. In contrast to cases without a seawall, the node or pseudo-node of wave height evolution can be found between the submerged breakwater and the seawall. Higher wave set-up occurs if the nodal or pseudo-nodal point appears near the submerged breakwater. We also examined the influence of the porosity and friction factor of the submerged permeable breakwater on wave transformation and set-up.     

一维大变形主、次固结耦合新模型

针对传统一维大变形固结理论不能考虑土体次固结变形的缺点,在连续介质力学有限变形理论和土力学固结理论基础上,提出了全新的一维大变形主、次固结耦合模型。该模型采用超静孔压为控制变量,可以分别表述为欧拉坐标、拉格朗日坐标和固相体积坐标等三种形式,其中均可以显式地反映出土体重度变化、固相速度、地表沉降速度、地下水位变化、荷载变化和次固结等6个因素的影响。在不同简化条件下,新模型可以分别退化为以往学者提出的相关模型。当不考虑次固结效应时,新模型与经典的Gibson大变形固结模型完全等价。从便于工程应用的角度而言,新模型中仅考虑次固结和外荷载变化的退化形式更便于求解,可作为一种供参考的实用模型。     

Blasius方程中临界值的分析估计

证明了流体力学边界层理论中著名的Blasius方程的一个新性质,它在二维平板流的研究中具有重要的作用,由此建立它的积分形式,并给出方程中的临界值βc的分析估计。     

等效水力隙宽和水力梯度对岩体裂隙网络非线性渗流特性的影响

等效水力隙宽和水力梯度是影响岩体裂隙网络渗流特性的重要因素。制作裂隙网络试验模型,建立高精度渗流试验系统;求解纳维-斯托克斯方程,模拟流体在裂隙网络内的流动状态,研究等效水力隙宽和水力梯度对非线性渗流特性的影响。结果表明,当水力梯度较小时,等效渗透系数保持恒定的常数,流体流动属于达西流动区域,流量与压力具有线性关系,可采用立方定律计算流体流动;当水力梯度较大时,等效渗透系数随着水力梯度的增加而急剧减少,流体流动进入强惯性效应流动区域,流量与压力具有强烈的非线性关系,可采用Forchheimer方程计算流体流动。随着等效水力隙宽的增加,区别线性和非线性流动区域的临界水力梯度呈幂函数关系递减。当水力梯度小于临界水力梯度时,控制方程可选立方定律;当水力梯度大于临界水力梯度时,控制方程可选Forchheimer方程,其参数A和B可根据经验公式计算得到。其研究结果可为临界水力梯度的确定及流体流动控制方程的选取提供借鉴意义。