有机污染物在含土工膜复合隔离墙和含水层系统的运移半解析模型

为评估土工膜与土-膨润土组成的复合隔离墙对有机污染物的防污性能,建立了污染物在源区域发生降解扩散时透过复合隔离墙的一维瞬态扩散模型。通过拉普拉斯变换和Talbot数值反演获得了模型解。第3类入流边界条件下,考虑源区的扩散和降解时,第1类隔离墙(土工膜/土-膨润土)和第2类隔离墙(土-膨润土/土工膜/土-膨润土)100 a的击穿浓度分别较源浓度恒定时减少了59%和53%。由于假设土工膜渗透系数高于10^–12 m/s,污染物在复合隔离墙中的对流作用不能忽略,因此,采用第2类复合隔离墙使土-膨润土能更好地发挥主导隔离作用。将土工膜的渗透系数从10^–10 m/s减少到10^–16 m/s,第1类复合隔离墙的击穿时间从26 a增加到188 a,第2类复合隔离墙的击穿时间从32 a增加到81 a。利用抽水井调整墙体内外为负水头差及提高源区的污染物降解能力,可使隔离墙的防污性能得到较大提升。     

1986—2021年雅江-尼洋河交汇段辫状河道演变过程

雅鲁藏布江与尼洋河交汇段的辫状河道形态及变化主要受尼洋河入汇、河谷边界条件及滨河植被的影响,时空变化的差异性和复杂性值得研究。利用1986—2021年遥感数据提取河道-沙洲-植被主要形态参数(主河道迁移率、弯曲率、面积等)分析辫状河道形态的复杂特征与演变过程。结果表明：主河道迁移率最大为483 m/a,曲率减小3.43%,植被面积逐年增加到8.05 km²。米林—派镇段辫状河道受边界条件和水文条件的差异影响,主河道在尼洋河汇入前呈横向单侧摆动,汇入后雅鲁藏布江主河道同时发生横向与纵向摆动;2013年以来该河段的植被仅有11.8%的覆盖区域对河床产生抑冲促淤的正向反馈并促进沙洲的稳定。     

基于负表皮层影响的径向溶质运移模型构建与新求解方法

对于模拟示踪试验的溶质运移,传统方法采用细网格离散小尺度表皮层,存在网格数量多和计算时间长等问题。由于表皮层外缘的参数变化剧烈,即使采用细网格,仍然存在显著误差。因此发展基于完整井的新示踪试验溶质运移模型,提出新的瞬态Robin边界条件反映表皮层吸附/解吸溶质的影响,实现表皮层的无网格。应用拉普拉斯变换导出模型的解析解,通过有限元法建立基于非完整井的数值解。结果显示新瞬态Robin边界条件精准反映负表皮层的影响。对于求解反演问题,表皮层宽度的估计范围由0.45～0.54 m缩至0.47～0.48 m,纵向弥散度由0.6～10.0 m缩至6.4～7.7 m,参数估计的可靠度显著提升。表皮层的无网格节省97%的数值解计算时间。     

多年冻区管道通风路基温度边界条件及温度场实测研究

管道通风路基在多年冻土地区是一种良好的主动保护冻土工程措施, 但在应用数值方法研究其长期效果中, 对边界条件的选取存在着困难, 已有的计算结果也缺乏实体工程的验证. 基于青藏铁路北麓河实体试验工程, 对通风管中气温以及管壁温度进行了分析和拟合. 结果表明: 通风管中气温年平均值普遍高于环境气温, 平均高出1.6～1.8 ℃, 但在正温期管中气温与环境气温的最高值相差不大（＜1 ℃）, 而负温期相差较大（达到2 ℃）; 路基阳坡面下0.5 m深度地温高于阴坡面3.5～5.5 ℃. 普通路基填筑后在抬升多年上限的同时, 也升高了上限附近的地温, 且地温场在整个范围内表现出横向上显著的不对称性. 通风路基可以对下伏多年冻土起到良好的主动保护作用, 表现为冻土上限的抬升和地温的降低. 但通风管埋设位置较高的情况下, 路基地温场横向不对称范围涉及到原地面以下. 降低通风管的埋设高度的情况下, 边坡面换热导致的温度横向不对称范围被限制在通风管以上的土体中, 而下部土体温度场横向不对称性将得到极大改善.     

下边界条件为水位流量关系的扩散波洪水演算方法研究

经大量的理论分析和实际洪水演算表明,扩散波可较好地反映洪水波的运动特点。本文基于河道下断面的水位流量关系和连续方程的联立求解,导出了水位流量关系型的下边界条件。利用Laplace变换,得到了扩散波方程在该边界条件下的解析解。文中的求解过程表明,可利用该法进行洪水演算,关键在于如何与实际洪水结合,这个问题有待进一步地分析研究。     

地幔对流环型场的激发(Ⅰ)

讨论了地幔内部的粘滞度及施加在地表和CMB 的边界条件对地幔对流环型场的激发分析表明,当粘滞度侧向均匀时,环型场与极型场自然解耦,且环型场不影响重力位,当粘滞度侧向不均时,环型场与极型场耦合在一起.两者共同影响重力位.当引入板块运动速度时,边界条件非零,也能激发环型场;对侧向均匀粘滞度地幔,零边界条件不能激发环型场     

二维横向各向同性介质中人为边界反射的吸收——差分法弹性波场模拟

本文对有限差分法弹性波场三分量正演模拟中遇到的人为边界反射问题进行了研究,分别获得2-D横向各向同性介质中U_x,U_z,与U_y位移分量模拟时的边界吸收条件.从吸收边界上反射系数与入射角之间的关系可以说明,本文提出的边界吸收条件的吸收能力比较强.     

浅水流动计算中—阶有限体积法Osher格式的实现

近年来,一阶有限体积法Osher格式已在二维浅水明流的一批模型问题和应用实例中获得成功。本文首次讨论其算法实现的种种问题。核心是建立单元水力模型-阶梯流,在数学上可用一类特殊的黎曼问题来描述。将该问题化作气体动力学中的黎曼问题近似求解,然后对结果加以校正。还在理论分析和数值试验的基础上详细讨论了各种外部边界条件、内部边界和动边界的处理,构成完整的算法。     

Ｅ网格变量分布下差分格式的性质

本文用线性化的浅水波方程系统地分析、讨论了Ｅ网格变量分布下差分计算的性质；指出了目前有关讨论Ｅ网格差分性质在方法上存在的问题；给出了一种解决Ｅ网格中两个Ｃ网格分离的方法；提出了一种Ｅ网格用于海洋计算时，海洋侧边界条件的处理方法。初步的数值试验证明了这些方法对抑制Ｅ网格中的计算躁声是有效的。     

关于应力边界条件下塑性问题的应力场求法：答王元汉同志

讨论文(见本刊本期)指出文[1]有两个观点是不正确的,但表明第二个不正确的观点实际上乃是由第一个造成的,因而下面仅对第一个观点的正确性作二点说明。 1.讨论文未对新方法的基本方程提出异议,因而我觉得首先可以从基本方程的求解上再做一些说明,新方法的基本方程如下: