边界条件对碎石层降温效果及机理的影响

对平均粒径为22.1 cm, 厚度1.3 m, 边界为开放和封闭的碎石层进行了一系列的实验. 从环境温度为最高和最低时刻封闭碎石层内的温度和速度分布可知: 当碎石层上表面温度低于下表面时,碎石层自下向上的散热由空气对流和碎石间的热传导来完成; 当上表面温度高于下表面时, 自上向下传递的热量由碎石接触面间的热传导完成, 此时, 由于其内部的空气几乎静止, 能阻隔热量的传入, 因此封闭边界的碎石层具有热半导体特性. 而开放边界的碎石层, 当平均温度为 0.5℃的空气从上表面吹过时, 碎石体内的热量传递主要是靠强迫对流来完成, 热半导体效果不明显, 不利于使其下面的冻土降温.     

局部和全局边界条件的组合

时域局部人工边界是以行波概念为基础的,原则上不能用于模拟零频分量（消减运动）。Ｗｏｌｆ提出了能同时模拟零频分量和行波的双渐近－多向透射边界。本文指出,在这一组合中的双渐近边界的高频渐近项是不合理的,还查能降低模拟精度。据此,本文提出一种概念清楚、简便实用的组合方案－多次射透－静力边界,并通过数值试验对此方案进行检验。     

地质雷达正演中的频散压制和吸收边界改进方法

从麦克斯韦方程组出发,建立了地质雷达的时域有限差分法（ＦＤＴＤ）数学模型,导出了理想 散关系和超级吸收边界条件。理想频散关系 考虑了ＦＤＴＤ法的收敛性和稳定性,也考虑了高频电磁波在Ｙｅｅ氏网格中的传播特点;超吸收边界条件则用磁场分量来提高电场分量精度。数值试验表明,理想频散关系能真实地反映雷达波在地下介质中的传播规律,超吸收边界条件能有效减小截断边界的伪反射,提高正演精度。将之应用于实际计算,取得     

冻融作用对系统与环境间能量交换的影响

通过求解古典Neumann问题,研究土壤在冻结和融化过程的热交换特征,从理论上分析了冻融作用对系统与环境间能量交换的影响.以草炭亚粘土和亚粘土为例,计算了广泛环境条件下的冻融作用对系统温度变化和与环境间热交换的影响.结果表明,冻融作用使相应过程表层温度梯度加大,热交换强度增强.     

构造分析方法探讨

构造不连续带是与构造区段相对应的构造单元,在构造变形中起着极为重要的作用。①控制其两侧的构造变形,成就不同的构造样式;②运动面的作用,两则地块相向、相背或旋转运动,错综复杂的构造图像,就是构造不连带控制、调整的结果,大尺度的构造不连续带围限同尺度或小尺度的构造区段,大尺度的构造区段包容小尺度的构造不连续带,这种围限与被围限、包容与被包容的关系,揭示了不同尺度构造的控制与序次的关系。提出构造不连续带     

非规则极坐标系弹性波场正演模拟（英文）

在研究隧道、井孔等特殊地质环境下的波场传播时,极坐标正演模拟比传统的笛卡尔坐标正演模拟精度更高。而且在实际情况下,极坐标环境大多是不规则的,为了解决这一问题,提出了一种非规则极坐标系正演模拟方法提高模拟精度。该方法首先将具有不规则表面的空间在非规则极坐标系下进行网格剖分,并通过变换将不规则表面映射为规则表面,然后在非规则极坐标系下计算波场。在求解非规则极坐标系下的波动方程时,利用全交错网格对其进行差分离散求解。此外,利用人工吸收边界、极坐标系自由边界和圆周边界条件对虚假边界反射进行吸收。选取了三个模型来验证该方法,结果显示,相比于常规极坐标系正演模拟,本论文非规则极坐标系弹性波场正演模拟方法在对具有不规则表面的极坐标空间进行正演模拟时,可以更获得准确、稳定的波场模拟结果。非规则极坐标系弹性波全波形反演结果进一步证明了本文提出的弹性波正演模拟方法可以准确地模拟不规则极空间中的波场。     

复杂3D瞬变电磁场FDTD模拟中需要解决的问题

为了实现缓变瞬态电磁场的时域有限差分(FDTD)模拟,引入了虚拟位移电流,使FDTD随时间步的递推能够进行。在以往较简单模型的计算中,利用时域瞬变场传播需要时间的性质,可将激励源化为初始条件代入,但是对于复杂3D模型,原有的方法不再适用,应当重新考虑源和边界条件的处理方法。由于实际源带有“斜坡”,故可采取将激励源模拟注入的方式;地下边界根据电磁波在导电媒质中快速衰减的特性设定,将改进的完全匹配层用于空中边界。分析表明,对于由空中和地下区域非均匀网格离散引起的时间稳定性问题,交替方向隐式FDTD是可行的解决办法。      

MT有限元模拟中截断边界的影响

利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。     

广义塑性梯度模型的变分原理和边界条件

在广义塑性力学的双屈服面模型的体积屈服面和剪切屈服面中引进应变梯度项,构造梯度依赖的双屈服面,可以建立广义塑性梯度模型的理论框架,此时的位移率 、剪切塑性乘子率 、体积塑性乘子率 成了各自独立的变量。为了建立相应的数值模型,构造了增量泛函,建立并证明了3类变量( , , )变分原理,得到应力边界条件,同时塑性乘子作为自然边界条件给出。     

VisualModflow模型在白涧铁矿区矿井涌水量预测中的应用

在分析确定白涧矿区水文地质边界条件的基础上,运用实测资料对Visu-almodflow模型进行识别和验证,同时修正了含水层的参数,修正后的模型较为可靠,可用于白涧矿区涌水量的预测。依据矿区奥灰含水层边界条件和渗透系数的不同,设计了-200m和-400m水平面不同边界及渗透系数三个方案对矿井的涌水量进行预测。结果表明:方案2的涌水量比较合理,能较真实的反映白涧铁矿矿体分布区水位疏降至-400m时的涌水量情况,预测其涌水量为7505m3/h。