首页 理论教育热流固耦合模型探讨

热流固耦合模型探讨

【摘要】:为研究传热、水流和力学等参数的横观各向同性耦合效应,采用FLAC3D软件对该地下岩石实验室中的黏土岩的热流固耦合效应进行了数值仿真模拟。图12.3热-流-固耦合关系图温度改变引起的有效应力通过改变计算模型的应变量来实现;孔隙压力对有效应力的影响由太沙基有效应力原理来计算。计算采用的黏土岩热-流-固耦合关系见图12.3。

为研究传热、水流和力学等参数的横观各向同性耦合效应,采用FLAC3D软件对该地下岩石实验室中的黏土岩的热流固耦合效应进行了数值仿真模拟。

在FLAC3D中,热流固耦合分析通过求解以下平衡方程加以实现:

能量平衡方程

流体质量平衡方程

Cauchy运动方程

式中:qi,i为热流量向量,W/m2;qv为体积热源,W/m3;ρ为孔隙介质干密度,kg/m3;cv为体积比热,J/(kg·℃);T为温度,℃;qi,i为单位流量向量,m/s;qv为体积流体源强度,1/s;ζ为孔隙介质体内单位体积的流体体积改变量;t为时间,s;bi为单位体积的体力分量;dvi/dt为速度对时间的导数

孔隙压力增量与温度和体积应变增量等的关系式为

式中:M为比奥模量,N/m2;α为比奥系数;β为材料的不排水体积热膨胀系数,1/℃;s为饱和度;φ为孔隙率

图12.3 热-流-固耦合关系图

温度改变引起的有效应力通过改变计算模型的应变量来实现;孔隙压力对有效应力的影响由太沙基有效应力原理来计算。

式中:σij为总应力;为有效应力;α为比奥系数;Δp为孔隙压力增量;Kb为饱和岩体的体积模量;βTb为饱和岩体的线性膨胀系数;ΔT为温度变化量;δij为克罗内克函数。

大量研究表明孔隙水运动和岩体骨架变形对饱和岩体介质的温度影响程度较弱,因此耦合计算模型中忽略了孔隙流动和岩体骨架变形引起的温度改变。计算采用的黏土岩热-流-固耦合关系见图12.3。