近断层应力变化快速计算方法及对强余震空间分布的指示意义:以M_w7.9汶川地震为例

构造地震一般由断层摩擦失稳所致.断层内部及周边所累积的剪切形变则通过同震滑动位移得到局部释放.因此,震后断层内部及近断层周边的静态剪切应力变化量的空间分布可通过断层面上的滑动位移分布计算得到.本文采用傅氏变换方法(FTM)计算单一有限断层同震滑移场所形成的静态剪切应力变化分布,近断层两侧的应力变化计算可由波数域内应力近似算法获得.结果表明,FTM快速有效、易于实现,有效地避免了常规应力计算中奇异值的出现.以2008年Mw7.9中国汶川大地震为例,采用前人所得有限断层滑动位移模型,得到了断层面和近断层周边准3D剪切应力分布解,并同主震后中强余震的空间分布特征作了比较.结果表明,大部分的中强余震震源位置处于剪切应力变化值为正的区域,由同震位移所产生的局部应力降峰值和均值大小同应力变化的正值大致相当,从而表明了快速且有效地计算断层内部及近断层附近的应力变化分布可以为主震后强余震发生的潜在区域提供指示意义.需要强调的是,应力变化空间分布特征的获取强烈地依赖于有限断层滑移模型解.有关滑动位移反演解的多解性对应力变化计算结果的影响,本文作了必要的讨论.

管内液膜流动失稳现象的直接数值模拟

为分析液膜流动时垂直液膜的重力加速度分量和表面张力对扰动波振幅的影响,通过FTM模拟了液膜在重力作用下流动失稳过程。结果表明:重力分量和表面张力系数越大,扰动振幅衰减得越快;并且表面张力系数越小,扰动振幅波动的频率越高。此外还分析了促进液膜流动失稳的两个因素,即密度比和雷诺数Re。当增大密度比时,液膜会出现两种极端的变化:一种是液膜不再完全覆盖底面,会形成锲状液膜,随着速度变大和振幅脉动,当锲状液膜前锋出现新的扰动波时,新的扰动会与初始扰动发生干涉产生二次波动;另一种是增加密度比时,液膜将会在剪切力的作用下产生飞脱的现象,在Re数较大时,液膜流动失稳会在剪切力的作用下产生K-H(开尔文-赫姆霍兹)不稳定性,且初始扰动对形成开尔文-赫姆霍兹波浪具有诱导作用。