Laser Interferometer Space Antenna double black holes: dynamics in gaseous nuclear discs

We study the inspiral of double black holes, with masses in the Laser Interferometer Space Antenna ( LISA ) window of detectability, orbiting inside a massive circumnuclear, rotationally supported gaseous disc. Using high-resolution smoothed particle hydrodynamics simulations, we follow the black hole dynamics in the early phase when gas-dynamical friction acts on the black holes individually, and continue our simulation until they form a close binary. We find that in the early sinking the black holes lose memory of their initial orbital eccentricity if they corotate with the gaseous disc. As a consequence, the massive black holes bind forming a binary with a low eccentricity, consistent with zero within our numerical resolution limit. The cause of circularization resides in the rotation present in the gaseous background where dynamical friction operates. Circularization may hinder gravitational waves from taking over and leading the binary to coalescence. In the case of counter-rotating orbits, the initial eccentricity (if present) does not decrease, and the black holes may bind forming an eccentric binary. When dynamical friction has subsided, for equal mass black holes and regardless their initial eccentricity, angular momentum loss, driven by the gravitational torque exerted on the binary by surrounding gas, is nevertheless observable down to the smallest scale probed (≃1 pc). In the case of unequal masses, dynamical friction remains efficient down to our resolution limit, and there is no sign of formation of any ellipsoidal gas distribution that may further harden the binary. During inspiral, gravitational capture of gas by the black holes occurs mainly along circular orbits; eccentric orbits imply high relative velocities and weak gravitational focusing. Thus, the active galactic nucleus activity may be excited during the black hole pairing process and double active nuclei may form when circularization is completed, on distance scales of tens of parsecs.     

冰封期水动力水质模型在松花江水污染事件中的应用

根据研究流域特点建立了冰封期水动力水质模型,结合冰封期水位测量的特点,对水位模拟结果进行了修正,分别利用2006和2007年水文监测数据对模型进行了率定和验证,利用模型对2005年松花江硝基苯水污染事件进行模拟。计算结果表明:建立的冰封期水动力水质模型能够很好地适用于该地区,部分断面硝基苯质量浓度峰值的相对误差为5%~13%,峰值出现时间的绝对误差为2~3 h。利用模型预测了上游边界流量增加20%及减少5%对水质的影响,结果表明:流量增加导致污染团峰值提前22 h到达摆渡河断面,硝基苯质量浓度峰值降低了15%,流量减少导致峰值延迟6 h到达摆渡河断面,质量浓度峰值升高了4.6%。     

冻融循环作用下黄土孔隙率变化规律

冻融循环作用通过破坏黄土颗粒的大小和土体的骨架及组构影响黄土的孔隙率. 通过颗粒分析实验发现冻融循环后土体颗粒有变小趋势,且集中分布在0.01～0.05 mm,在冻融10次后颗粒大小趋于稳定. 对经历冻融循环后的黄土进行压汞实验,发现随着冻融次数的增加,土体孔隙中的大孔径先减少,后增多;小孔径先增多,后减少,最后向5～10 μm范围内集中. 土的孔隙率也随冻融次数的增加而发生变化. 原状黄土的孔隙率在冻融5次时达到最小值,在冻融10次时达到最大值,冻融10次之后原状黄土的孔隙率在40%上下波动,达到稳定状态. 重塑黄土在冻融3次时孔隙率减少了6%并达到最小值,在冻融5次时孔隙率增加3%并达到最大值,最后在40%上下波动,达到稳定状态. 渗透系数的变化曲线与通过压汞实验获得的孔隙率变化曲线相似,印证了随着冻融次数的增加黄土的孔隙率呈先减小,后增大,然后趋于稳定的变化规律.     

北部湾防城港东湾白骨壤潮滩近底层动力响应台风作用的耗散过程

揭示极端天气影响的红树林潮滩潮流能量耗散过程是理解生物海岸演变及滨海湿地生态修复工程的核心内容。以北部湾防城港东湾白骨壤红树林潮滩为例,基于声学多普勒流速仪（ADV）获取区域2020年8月连续6天的水动力和白骨壤红树林下垫面植被实测数据,分析白骨壤潮滩近底层动力响应台风“森拉克”的耗散过程。结果表明：1）正常天气涨、落潮期间,自光滩到白骨壤红树林林内近底层湍流动能变化具有潮汐不对称性特征。涨潮期间光滩—白骨壤红树林边缘、白骨壤红树林边缘—林内潮流挟沙能量分别通过泥沙沉降以及搬运泥沙两种方式沿程消耗,落潮期间潮流挟沙能量主要以泥沙净沉降的方式消耗。2）白骨壤通过其潮滩表层向上发育长约10 cm的呼吸根以降低湍流垂向紊动;红树林边缘的枝、叶通过影响水平方向上的水流动力,致使潮流挟沙能耗降低。3）与正常天气比较,台风“森拉克”期间东湾自光滩到白骨壤红树林林内近底层水体流速无明显变化,但流向偏转幅度明显变大。同时湍流动能的耗散率和用以搬运泥沙为主的潮流挟沙能耗亦均增大。     

海拉尔盆地贝尔凹陷兴安岭群储层粘土矿物组成及成因研究

油田注水开发是保持地层能量,提高采收率的有效手段.海拉尔地区由于粘土矿物含量高,在开发过程中储层遇水易泥化、坍塌,严重影响油气的有效开采.通过X射线衍射、偏光显微镜观察和扫描电镜分析,确定了海拉尔盆地贝16断块兴安岭群储层自上而下存在三种粘土矿物组合:蒙脱石 高岭石 (绿泥石)→高岭石 伊/蒙间层 伊利石 (绿泥石)→伊利石 绿泥石 绿/蒙间层 伊/蒙间层.粘土总量在9.1%～69.1%之间,由地层浅部到深部呈现由高到低的变化规律.详细研究了水敏性极强的钠蒙脱石及其间层矿物,认为凝灰物质的水解蚀变以及长石和中、酸性岩屑在成岩过程中的次生变化是造成海拉尔地区粘土含量高的主要因素.     

锦屏大理岩蠕变损伤演化细观力学特征的数值模拟研究

岩石的蠕变损伤和断裂是岩石流变效应的重要表现形式,但其损伤演化过程往往难以直观观测,为此,采用二维颗粒流数值模拟方法(PFC2D)对岩石的蠕变损伤和断裂的细观力学机制进行了分析。在锦屏大理岩室内试验基础上,利用颗粒流应力腐蚀模型(PSC),建立了能反映其短期和长期强度特征的柱状岩样数值模型,并开展了大量数值试验。结果表明,蠕变损伤的演化过程与暂态的损伤演化过程具有明显的差异。在岩样蠕变损伤过程中,其内部微裂纹多沿加载方向开裂且分布均匀,先快速增加再稳定扩展,最后则发生快速断裂。当荷载较小时,岩样宏观上呈现劈裂破坏特征,当荷载较大时,岩样呈现剪切破坏特征。在岩样蠕变损伤初始和稳定演化阶段的前期,荷载大小对岩样的损伤演化过程影响不大;在稳定演化阶段的后期至断裂过程中,低荷载下岩样的损伤增速比高荷载下快。     

Large-scale galactic turbulence: can self-gravity drive the observed H i velocity dispersions?

Observations of turbulent velocity dispersions in the H  i component of galactic discs show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high-resolution adaptive mesh simulations of gaseous discs embedded within dark matter haloes to explore the roles of cooling, star formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disc slowly cools until gravitational and thermal instabilities give rise to a multiphase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large-scale driving of global non-axisymmetric modes as well as cloud–cloud tidal interactions and merging. At low star formation rates these processes alone can explain the observed H  i velocity dispersion profiles and the characteristic value of  ∼10 km s⁻¹  observed within a wide range of disc galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star formation rate per unit area  ≳10⁻³ M_⊙ yr⁻¹ kpc⁻²  .     

The biaxial compression mechanical properties of crushed rock

Crushed rock subgrade, as one of the roadbed-cooling methods, has been widely used in the Qinghai-Tibet Railway. Much attention has been paid on the cooling effect of crushed rock; however, the mechanical properties of crushed rock are somehow neglected. Based on the discrete element method, biaxial compression test condition for crushed rock is compiled in FISH language in PFC^2D, and the natural shape of crushed rock is simulated with super particle "cluster". The effect of particle size, crushed rock strength and confining pressure level on overall mechanical properties of the crushed rock aggregate are respectively analyzed. Results show that crushed rock of large particle size plays an essential framework role, which is mainly responsible for the deformation of crushed rock aggregate. The strength of gravel has a great influence on overall mechanical properties which means that strength attenuation caused by the freeze thaw cycles cannot be ignored. The stress-strain curves can be divided into two stages including shear contraction and shear expansion at different confining pressures.     

粗颗粒间液桥毛细力演化规律的动态计算方法

水力特性是非饱和土力学理论与工程的重要课题之一,土样水力特性的变化过程本质上可反映为土颗粒间毛细力的演化规律。为此,以粗颗粒为研究对象,将其简化为一对不等径球体颗粒,而其间的水分形态可视为形如圆环的液桥,不考虑颗粒重力和浮力的影响。以Young-Laplace方程为基础,先推得计算液桥毛细力的控制方程组,再结合液桥的无量纲体积最大、最小值及其外曲率半径的割线迭代算法提出求解毛细力的动态计算方法,进而研究毛细力与颗粒间距、颗粒半径比和液桥体积的无量纲关系,结果表明：当颗粒间距一定时,液桥的毛细力随其体积和颗粒半径比增大均呈递增趋势;当液桥体积一定时,其抗拉刚度随颗粒间距增大呈递减趋势。最后,利用已有文献中液桥毛细力与颗粒间距的实测关系,验证了该动态计算方法在表征液桥从形成至断裂时毛细力演化规律的有效性。