共查询到18条相似文献,搜索用时 109 毫秒
1.
任意形冒落非均质采空区流场流态数值模拟 总被引:8,自引:2,他引:8
为研究冒落非均质采空区流场流态及其随边界条件的变化,基于多孔介质渗流方程建立了采空区漏风流态的有限元数值模型,结合现场实例,从宏观上描绘了复杂几何形状采空区的漏风流态(风压分布等值线和流函数线),给出更准确的流态图形分布解.指出用流函数Laplace方程修正函数,对复杂形状和沿多边界方向非均质的采空区流场计算,仍能满足流网正交,其正交性可通过剖分网格精度来控制,得到采空区漏风强度沿固壁边界附近的冒落非压实区相对较大,工作面上下端附近处的漏风强度很大,并直通采空区深部。 相似文献
2.
前进式U形通风采空区风流与自燃的数值模拟 总被引:1,自引:0,他引:1
基于渗流方程、多组分气体渗流-扩散方程和多孔介质渗流传热方程,建立前进式开采U型
通风采空区自然发火模型. 通过迭代方法计算沿空巷道的侧漏风量,以此边界条件数值模拟
计算出前进式开采U型通风采空区的漏风规律;根据瓦斯-Ohttp://d.wanfangdata.com.cn/Periodical_lxysj200804005.aspx 前进式开采U型通风;采空区;风流规律;自燃温度场;有限元; 基于渗流方程、多组分气体渗流-扩散方程和多孔介质渗流传热方程,建立前进式开采U型
通风采空区自然发火模型. 通过迭代方法计算沿空巷道的侧漏风量,以此边界条件数值模拟
计算出前进式开采U型通风采空区的漏风规律;根据瓦斯-Ohttp://d.wanfangdata.com.cn/Periodical_lxysj200804005.aspx 前进式开采U型通风;采空区;风流规律;自燃温度场;有限元; 基于渗流方程、多组分气体渗流-扩散方程和多孔介质渗流传热方程,建立前进式开采U型
通风采空区自然发火模型. 通过迭代方法计算沿空巷道的侧漏风量,以此边界条件数值模拟
计算出前进式开采U型通风采空区的漏风规律;根据瓦斯-Ohttp://d.wanfangdata.com.cn/Periodical_lxysj200804005.aspx 前进式开采U型通风;采空区;风流规律;自燃温度场;有限元; 基于渗流方程、多组分气体渗流-扩散方程和多孔介质渗流传热方程,建立前进式开采U型通风采空区自然发火模型.通过迭代方法计算沿空巷道的侧漏风量,以此边界条件数值模拟计算出前进式开采U型通风采空区的漏风规律;根据瓦斯-O2-CO浓度和温度分布状态及其变化来描述自燃过程.在计算中,采空区非均质性按自然冒落考虑,兼顾瓦斯涌出对自燃的耦合作用;模拟得到,在前进式开采特有的3面漏风边界条件下,采空区自燃位置(最高温点)在开切眼冒落区进风一侧,自燃产生的有害气体和瓦斯随漏风流从回风边界渗出.在前进开采U型通风形式下,通风量在到达工作面时衰减很大,因此,前进开采布置一次推进距离不易过大;解决好沿空留巷的封堵材料和正确的堵漏措施,是保证前进式开采顺利进行的技术关键. 相似文献
3.
开区注氮采空区自燃温度场的数值模拟研究 总被引:4,自引:1,他引:4
为了研究开区注氮采空区遗煤自燃和被抑制的力学过程,用非均质多孔介质中的渗流连续性方程、气体弥散方程和综合传热方程的联立,建立了采空区注氮非定常数值模型.结合实例,用迎风格式有限元方法求解.给出了采空区的漏风流态、氮气流态,描绘了采空区自燃过程中氧、CO浓度和温度分布的变化过程.计算考虑了瓦斯涌出和工作面推进的影响.得到注氮条件下采空区高温区形状变窄;随着注氮量的提高,使自然发火期逐渐变长,直至不自燃.理论计算与实际情况吻合. 相似文献
4.
5.
饱和多孔介质中的混合有限元法和有限应变下应变局部化分析 总被引:1,自引:0,他引:1
对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元.固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值.基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式,导出了单元公式.进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法.对几何非线性分析,采用了共旋公式途径.数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能. 相似文献
6.
基于图像识别的裂隙煤层气非Darcy渗流模拟 总被引:1,自引:0,他引:1
煤层气在煤岩裂隙中流动时,由于其粘度系数较小,很容易产生非Darcy现象.因此,考虑煤层气的非Darcy渗流将会使得研究结果更加客观和准确,对煤层气开采具有实际意义.首先将实物裂隙高清拍照导入Matlab软件中进行图像二值化处理,再利用COMSOL软件对二值图像进行模型重构和数值模拟.模拟结果显示,煤层气渗流速度随着裂隙尺度的变化而不断改变,速度随时间不断增加,但区域渗流速度分布状态保持不变.与Darcy结果的对比分析可知,相同压力梯度下,考虑非Darcy因素的产气量要明显小于仅考虑Darcy因素的情况.同时,结果显示,在渗透率不同的煤裂隙中,煤层气的渗流速度不同.渗透率较大时,渗流速度相对较大,同样,渗透率较小时,渗透速度较小. 相似文献
7.
提出了一种不可压缩流体与弹性薄膜耦合问题的特征线分裂有限元解法. 首先, 给出了流场和结构的控制方程. 然后, 对流场、结构以及流固耦合的具体求解过程进行了描述. 其中, 流场求解采用改进特征线分裂方法和双时间步方法相结合的隐式求解方式, 并利用艾特肯加速法对每个时间步的迭代收敛过程进行了加速处理;结构部分的空间离散和时间积分分别采用伽辽金有限元方法和广义方法, 并通过牛顿迭代法对所得非线性代数方程组进行了求解;流场网格的更新采用弹簧近似法;流场、结构两求解模块之间采用松耦合方式.最后, 采用该方法对具有弹性底面的方腔顶盖驱动流问题进行了求解, 验证了算法的准确性和稳定性.此外, 计算结果表明艾特肯加速法可以显著地提高双时间步方法迭代求解过程的收敛速度. 相似文献
8.
水力压裂在页岩气开采中被广泛使用,采用数值方法研究压裂机理具有重要意义.基于连续-非连续单元法(CDEM) 和中心型有限体积法(FVM),提出解决水力压裂流固耦合问题的二维混合数值计算模型.该混合模型中,使用CDEM 求解应力场和裂缝扩展过程,使用FVM 求解裂隙渗流场.应力场裂缝扩展和渗流场均使用显式迭代求解, 并通过相互之间数据交换实现流固耦合.通过与KGD 理论模型进行对比, 验证数值模型的正确性.通过与颗粒离散元数值结果进行对比,验证数值模型的有效性.通过计算复杂缝网压裂模型,研究水力压裂机理,并说明该数值模型在水力压裂模拟中具有很好的前景. 相似文献
9.
10.
建立了一种基于压力梯度和渗流速度时间序列提取变质量破碎岩体渗透性参量(渗透率、非Darcy流β因子、加速度系数)的方法。首先,根据破碎岩体渗透性参量之间的幂指数关系和Forchheimer关系,建立了采样时刻渗透率的代数方程,并利用Newton切线法求得此代数方程的根;其次,对渗透性参量的参考值及幂指数进行优化。通过算例分析了该方法用于计算变质量破碎岩体渗透性参量的可行性及准确性。研究结果表明:随着细小颗粒的迁移流失,破碎岩体在恒定的压力梯度下,渗流速度和渗透率增加,非Darcy流β因子和加速度系数降低。 相似文献
11.
12.
基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合.验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明:非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。 相似文献
13.
It was derived that micro-scale amount level of average pore radius of clay changed from 0.01 to 0.1 micron by an equivalent concept of flow in porous media.There is good agreement between the derived results and test ones.Results of experiments show that flow in micro-scale pore of saturated clays follows law of nonlinear flow.Theoretical analyses demonstrate that an interaction of solid-liquid interfaces varies inversely with permeability or porous radius.The interaction is an important reason why nonlinear flow in saturated clays occurs.An exact mathematical model was presented for nonlinear flow in micro-scaie pore of saturated clays.Dimension and physical meanings of parameters of it are definite.A new law of nonlinear flow in saturated clays was established.It can describe characteristics of flow curve of the whole process of the nonlinear flow from low hydraulic gradient to high one.Darcy law is a special case of the new law.A math- ematical model was presented for consolidation of nonlinear flow in radius direction in saturated clays with constant rate based on the new law of nonlinear flow.Equations of average mass conservation and moving boundary,and formula of excess pore pressure distribution and average degree of consolidation for nonlinear flow in saturated clay were derived by using an idea of viscous boundary layer,a method of steady state in stead of transient state and a method of integral of an equation.Laws of excess pore pressure distribution and changes of average degree of consolidation with time were obtained.Re- suits show that velocity of moving boundary decreases because of the nonlinear flow in saturated clay.The results can provide geology engineering and geotechnical engineering of saturated clay with new scientific bases.Calculations of average degree of consolidation of the Darcy flow are a special case of that of the nonlinear flow. 相似文献
14.
V.A.F.Costa 《中国颗粒学报》2009,7(3)
This work reports on fluid flow in a fluid-saturated porous medium, accounting for the boundary and inertial effects in the momentum equation. The flow is simulated by Brinkman-Forchheimer-extended Darcy formulation (DFB), using MAC (Marker And Cell) and Chorin pressure iteration method. The method is validated by comparison with analytic results. The effect of Reynolds number, Darcy number, porosity and viscosity ratio on velocity is investigated. As a result, it is found that Darcy number has a decisive influence on pressure as well as velocity, and the effect of viscosity ratio on velocity is very strong given the Darcy number. Additional key findings include unreasonable choice of effective viscosity can involve loss of important physical information. 相似文献
15.
Effects of porous medium have been investigated on the steady flow of a third grade fluid between two stationary porous plates. The continuity and momentum equations along with modified Darcy??s law are used for the development of mathematical problem. The governing nonlinear problem is solved by a homotopy analysis method. The dimensionless velocity and shear stresses at the plates are analyzed. 相似文献
16.
The combined effects of the magnetic field, permeable walls, Darcy velocity, and slip parameter on the steady flow of a fluid in a channel of uniform width are studied. The fluid flowing in the channel is assumed to be homogeneous, incompressible,and Newtonian. Analytical solutions are constructed for the governing equations using Beavers-Joseph slip boundary conditions. Effects of the magnetic field, permeability,Darcy velocity, and slip parameter on the axial velocity, slip velocity, and shear stress are discussed in detail. It is shown that the Hartmann number, Darcy velocity, porous parameter, and slip parameter play a vital role in altering the flow and in turn the shear stress. 相似文献
17.
The effect of melting heat transfer on the two dimensional boundary layer flow of a micropolar fluid near a stagnation point embedded in a porous medium in the presence of internal heat generation/absorption is investigated. The governing non-linear partial differential equations describing the problem are reduced to a system of non-linear ordinary differential equations using similarity transformations solved numerically using the Chebyshev spectral method. Numerical results for velocity, angular velocity and temperature profiles are shown graphically and discussed for different values of the inverse Darcy number, the heat generation/absorption parameter, and the melting parameter. The effects of the pertinent parameters on the local skin-friction coefficient, the wall couple stress, and the local Nusselt number are tabulated and discussed. The results show that the inverse Darcy number has the effect of enhancing both velocity and temperature and suppressing angular velocity. It is also found that the local skin-friction coefficient decreases, while the local Nusselt number increases as the melting parameter increases. 相似文献
18.
PIV measurements for gas flow under gradient magnetic fields 总被引:1,自引:0,他引:1
Particle Image Velocimetry (PIV) techniques were developed to measure the convective N2-air flow under gradient magnetic fields. The velocity fields were calculated by the Minimum Quadratic Difference (MQD) algorithm and spurious vectors were eliminated by Delaunay Tessellation.The N2-air flow was measured as the magnetic flux density varying from 0~1.5T. A strengthened vortex flow of air was observed under the condition that the magnetic field was applied, and the velocity of N2 jet rose with the increase of the magnetic density. The experimental results show that the magnetic force will induce a vortex flow and cause a convection flow of the air mixture when both gradients of the O2 concentration and the magnetic field intensity exist. 相似文献