065三维非线性温度渗流耦合
- 格式:pdf
- 大小:904.16 KB
- 文档页数:8
中文摘要摘要随着国家经济建设的发展,越来越多的岩石工程涉及到多场耦合问题,裂隙岩体温度场-渗流场-应力场的耦合问题已经成为当前岩石工程的研究热点和研究难点。
由于实际岩石工程中裂隙岩体多场耦合作用所处地质环境的复杂性,以及室内试验方法的局限性,数值模拟方法是目前研究裂隙岩体多场耦合作用最有效的手段之一。
近场动力学理论是一种非局部理论,它采用空间积分法描述物质力学行为,在求解不连续问题时能够有效的避免解微分方程而产生的奇异性问题,对于处理材料的不连续问题具有较大的优势,同时由于近场动力学基于非局部理论,能很好的模拟热传导与地下水的渗流问题。
本文根据近场动力学的基本原理,建立了裂隙岩体破裂过程中温度场-渗流场-应力场耦合的数值计算模型,并编制相应的计算程序对该模型进行了验证和分析,论文主要的研究工作如下:①通过在基于键作用的近场动力学理论中引入切向键以模拟材料的剪切变形,从微观机理上完善了基于键作用的近场动力学本构模型,建立了近场动力学微观力学参数与宏观弹性常数之间的关系。
根据近场动力学柯西应力张量,建立了基于非普通状态的近场动力学理论损伤破坏模型,将物质点上的应力转化为键上的应力,并运用最大拉应力强度准则、莫尔-库仑强度准则双剪强度准则来判断键的破坏与否,再将每个物质点上断裂的键的数量与该物质点上包含键的总数的比值作为该物质点的损伤函数。
该模型成功的运用于模拟岩石三维裂纹的起裂、扩展和连接,并且得到了岩石破裂过程的应力应变曲线。
②根据热传导理论,并运用欧拉-拉格朗日方程推导了基于非局部理论的近场动力学热传导方程,建立了近场动力学微导热系数与材料宏观导热系数之间的关系;运用材料的热膨胀特性,将根据近场动力学热传导方程求解出的温度场转换为近场动力学物质点的变形梯度张量,再将变形梯度张量代入非普通状态近场动力学的力状态函数中,从而实现了岩体温度场与应力场的耦合。
③根据达西定律,推导了基于非局部理论的近场动力学渗流基本方程,运用质量守恒原理建立了一维和二维情况下宏观渗透系数与微观近场动力学渗透系数之间的关系。
低温冻融条件下岩体温度-渗流-应力-损伤(thmd)耦合模型研究及其在寒区隧道中的应用低温冻融条件下岩体温度-渗流-应力-损伤(thmd)耦合模型的研究,对于寒区隧道的建设和运营具有重要意义。
以下将从模型建立、数值模拟及实例应用等方面对其进行论述。
一、模型建立1.1 温度模型考虑到低温环境下的岩体温度变化,建立了温度耦合方程式,包括岩体内部热平衡方程和表面边界条件等,得到岩体内部温度分布。
1.2 渗流模型在岩体中加入水分后可导致岩体内的渗透性增加,因此需要建立渗流耦合方程式,包括岩体内部水分扩散方程和表面边界条件等,得到岩体内部渗流水压分布。
1.3 应力模型由于冻融环境下温度和水分变化,岩体中的应力状态也会发生变化,考虑到应力涉及到岩体的力学性质,建立了应力耦合方程式,包括岩体内部力学平衡方程和表面边界条件等,得到岩体内部应力分布。
1.4 损伤模型根据岩体断裂的损伤特性,建立了损伤耦合方程式,包括耗能损伤方程和裂纹扩展方程等,得到岩体的损伤特性及破坏状态。
二、数值模拟在建立了岩体的低温冻融条件下的温度-渗流-应力-损伤耦合模型后,进行了数值模拟。
通过数值模拟得到了不同参数情况下的模型响应情况和岩体的损伤状态等,为后面的实例应用提供了重要参考。
三、实例应用本研究将所建立的模型应用于寒区隧道的建设和运营,对隧道的环境适应性、寿命等方面进行了研究。
3.1 隧道建设阶段在隧道建设阶段,岩体的温度、水分等参数的变化会影响到隧道的施工和安全,采用所建立的模型可以对岩体环境进行预测,实现对施工和隧道的安全保障。
3.2 隧道运营阶段在隧道运营阶段,隧道内部的温度、水分等参数的变化也会影响到隧道的运营和安全。
采用所建立的模型,可以对隧道内部的环境进行监测和分析,及时调整运营管理策略,保障隧道的安全和可靠性。
综上所述,低温冻融条件下岩体温度-渗流-应力-损伤(thmd)耦合模型的研究,对于寒区隧道的建设和运营具有重要意义。
FLAC3D流固耦合(手册翻译)1.1简介FLAC3D通过具有渗流性的实体(比如土)来模拟流体的流动。
流动模型的建立可以独立于力学计算而自动完成,或者说可以与力学模型同时建立,这样就可以考虑流体与土体之间的相互作用。
流固耦合的一种类型是“固结”,即:空隙水压力逐渐消散而导致土体的沉降。
这个过程包括两种力学反映:一,空隙水压的改变导致有效应力的变化,这将影响到土体的力学反映(如:有效应力的减小可能导致塑性区的产生);二,力学实体中某一区域的流动会随着空隙水压的改变而改变。
该程序可以计算完全饱和情况下的流动,也可以模拟具有自由水面的流动。
模拟具有自由水面的流动时,自由水面以上的部分空隙水压等于0,气相将不参与计算。
对于不考虑毛细水压力颗粒较粗的材料可以采用这种模拟方法。
流体计算就有以下特点:1 根据各项同性和各项异性的渗流计算,相应采用两种流体运动定律。
流动中的null材料用来模拟流动范围内的非渗流材料。
2 不同区域可以拥有不同的流动模型(isotropic, anisotropic or null)和模型参数。
3 可以事先指定流体的压力、流量、非渗流区边界条件。
4 流体源可以以电源,也可以以体源的形式插入到材料中,这些源对应于流体的流入或流出,可以随着时间而变化。
5 对于完全饱和流动,可以采用显式和隐式两种算法,但对于非饱和流动则只能采用显示计算。
6 任何力学和温度计算模型都可以与流体模型一起使用,在耦合计算中,可以考虑饱和体的压缩性和热膨胀性。
7.流体与力学计算的耦合通过提供比奥系数来实现。
和不排水温度系数β8.与温度的耦合计算可以通过提供线性热膨胀系数αt(undrained thermal coefficient,可能翻译的不对)来实现。
9.热-流体计算以线性理论为基础,假定材料参数为常数,不考虑对流。
流体与实体的温度保持局部平衡。
非线性行为可以采用FISH语言改变孔隙压力、材料特性来实现。
ANSYS--热力耦合分析单元简介挑选了部分常用的,希望能方便大家的使用,其中自己翻译了一部分,不准确之处还望见谅,大家还可以继续补充哦!:SOLID5-三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。
本单元由8个节点定义,每个节点有6个自由度。
在静态磁场分析中,可以使用标量势公式(对于简化的RSP,微分的DSP,通用的GSP)。
在结构和压电分析中,具有大变形的应力钢化功能。
与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。
INFIN9-二维无限边界用于模拟一个二维无界问题的开放边界。
具有两个节点,每个节点上带有磁向量势或温度自由度。
所依附的单元类型可以为PLANE13和PLANE53磁单元,或PLANE55和PLANE77和PLANE35热单元。
使用磁自由度(AZ)时,分析可以是线性的也可以是非线性的,静态的或动态的。
使用热自由度时,只能进行线性稳态分析。
PLANE13-二维耦合场实体具有二维磁场、温度场、电场和结构场之间有限耦合的功能。
由4个节点定义,每个节点可达到4个自由度。
具有非线性磁场功能,可用于模拟B-H曲线和永久磁铁去磁曲线。
具有大变形和应力钢化功能。
当用于纯结构分析时,具有大变形功能,相似的耦合场单元有SOLID5、SOLID98和SOLID62。
LINK31-辐射线单元用于模拟空间两点间辐射热流率的单轴单元。
每个节点有一个自由度。
可用于二维(平面或轴对称)或三维的、稳态的或瞬态的热分析问题。
允许形状因子和面积分别乘以温度的经验公式是有效的。
发射率可与温度相关。
如果包含热辐射单元的模型还需要进行结构分析,辐射单元应当被一个等效的或(空)结构单元所代替。
LINK32-二维传导杆用于两节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
可用于二维(平面或轴对称)稳态或瞬态的热分析问题。
如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。
ansys热力耦合分析单元简介SOLID5-三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。
本单元由8个节点定义,每个节点有6个自由度。
在静态磁场分析中,可以使用标量势公式(对于简化的RSP,微分的DSP,通用的GSP)。
在结构和压电分析中,具有大变形的应力钢化功能。
与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。
INFIN9-二维无限边界用于模拟一个二维无界问题的开放边界。
具有两个节点,每个节点上带有磁向量势或温度自由度。
所依附的单元类型可以为PLANE13和PLANE53磁单元,或PLANE55和PLANE77和PLANE35热单元。
使用磁自由度(AZ)时,分析可以是线性的也可以是非线性的,静态的或动态的。
使用热自由度时,只能进行线性稳态分析。
PLANE13-二维耦合场实体具有二维磁场、温度场、电场和结构场之间有限耦合的功能。
由4个节点定义,每个节点可达到4个自由度。
具有非线性磁场功能,可用于模拟B-H曲线和永久磁铁去磁曲线。
具有大变形和应力钢化功能。
当用于纯结构分析时,具有大变形功能,相似的耦合场单元有SOLID5、SOLID98和SOLID62。
LINK31-辐射线单元用于模拟空间两点间辐射热流率的单轴单元。
每个节点有一个自由度。
可用于二维(平面或轴对称)或三维的、稳态的或瞬态的热分析问题。
允许形状因子和面积分别乘以温度的经验公式是有效的。
发射率可与温度相关。
如果包含热辐射单元的模型还需要进行结构分析,辐射单元应当被一个等效的或(空)结构单元所代替。
LINK32-二维传导杆用于两节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
可用于二维(平面或轴对称)稳态或瞬态的热分析问题。
如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。
LINK33-三维传导杆用于节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
《低渗透非线性渗流规律研究》篇一一、引言随着能源需求持续增长,低渗透储层因其巨大的开发潜力逐渐成为油气开采的重要领域。
低渗透非线性渗流规律研究是当前石油工程领域的重要课题,对提高采收率、优化开发策略具有重大意义。
本文旨在深入探讨低渗透非线性渗流的基本原理、影响因素及研究方法,以期为相关领域的研究和实践提供理论支持。
二、低渗透非线性渗流基本原理低渗透非线性渗流是指低渗透储层中流体在多孔介质中的流动过程,其流动规律与常规线性渗流存在显著差异。
低渗透储层具有孔隙度小、渗透率低、非均质性强等特点,导致流体在其中的流动表现出明显的非线性特征。
这种非线性渗流规律主要受储层物理性质、流体性质及外部条件等因素的影响。
三、影响低渗透非线性渗流的主要因素1. 储层物理性质:储层的孔隙度、渗透率、孔喉比等物理性质对非线性渗流具有重要影响。
其中,孔隙度和渗透率是决定流体流动能力的主要因素,孔喉比则影响流体的传输效率。
2. 流体性质:流体的粘度、密度、表面张力等性质也会对非线性渗流产生影响。
粘度较大的流体在低渗透储层中流动时,更容易表现出非线性特征。
3. 外部条件:温度、压力等外部条件的变化也会对非线性渗流产生影响。
例如,温度升高可能导致流体粘度降低,从而改变渗流规律。
四、低渗透非线性渗流规律研究方法1. 实验研究:通过室内实验,模拟低渗透储层中的流体流动过程,观察并记录非线性渗流现象,分析影响因素及作用机制。
2. 数值模拟:利用数值模拟软件,建立低渗透储层的地质模型和流体流动模型,通过计算分析非线性渗流规律及影响因素。
3. 理论分析:结合储层物理性质、流体性质及外部条件等因素,建立非线性渗流的数学模型和物理模型,进行理论分析和预测。
五、研究现状及展望目前,低渗透非线性渗流规律研究已取得一定成果,但仍存在诸多挑战和问题。
未来研究方向包括:深入探讨低渗透储层的物理性质和流体性质对非线性渗流的影响机制;优化实验和数值模拟方法,提高研究精度和可靠性;开发适用于低渗透储层的开采技术和策略,提高采收率等。
断层突水温度-非线性渗流-应力耦合模型研究刘伟韬;赵吉园;丁希阳【摘要】为了研究深部矿井高地应力、高承压水压及高地温条件下断层突水过程中的底板承压水运移规律和地温变化规律,运用COMSOL建立温度场(T)-非线性渗流场(H)-应力场(M)耦合模型,其中非线性渗流方程耦合Darcy定律,Brinkman方程和Navier-Stokes方程.结果表明:模型能够很好地描述温度场、渗流场和采动应力场之间的相互作用关系以及底板承压水从含水层经过断层破碎带到达巷道的非线性运移过程;采动应力场和渗流速度对断层带温度场的影响很大,通过监测断层带的温度变化,可以预测断层带渗流场的变化,从而预测断层突水.【期刊名称】《煤矿安全》【年(卷),期】2018(049)011【总页数】5页(P59-62,68)【关键词】断层突水;THM耦合;数值模型;非线性流动;温度【作者】刘伟韬;赵吉园;丁希阳【作者单位】山东科技大学矿业与安全工程学院,山东青岛 266590;山东科技大学矿业与安全工程学院,山东青岛 266590;山东济矿鲁能煤电股份有限公司阳城煤矿,山东济宁 272000【正文语种】中文【中图分类】TD745+.2在深部岩体高地应力、高承压水压及高地温的条件下,深部煤炭资源开采诱导的断层突水问题是我国煤矿安全生产的重大隐患。
目前,在煤层采动引起的断层突水预测上,许多学者从理论、数值模拟、实验和实测等方面进行了研究,得到了断层突水相关控制准则[1],建立了断层突水数值模型[2],进行了断层突水相关试验[3-4],对水压和温度进行实测,制定了预警等级及其判识准则[5]。
在THM耦合研究方面,国外学者Tsang C F基于Biot固结理论和质量、能量、动量守恒原理以及其他不同假设,提出了各种不同的THM耦合力学模型[6]。
国内学者在多孔材料小变形和热弹性线性理论的基础上推导出了1组较为完全THM耦合方程组[7],建立了热-流-固耦合数学模型,并进行了数值求解[8]。
《非线性渗流方程解析方法研究及应用》篇一一、引言渗流现象在自然界和工程领域中广泛存在,如地下水流动、油藏开发等。
非线性渗流方程是描述渗流现象的重要数学模型,其解析方法的研究对于理解渗流机制、预测渗流行为以及优化工程实践具有重要意义。
本文旨在探讨非线性渗流方程的解析方法及其应用,以期为相关领域的研究和实践提供理论支持。
二、非线性渗流方程的基本概念非线性渗流方程是一类描述多孔介质中流体流动的偏微分方程,其形式复杂,包含了流体性质、介质特性以及边界条件等多种因素。
根据不同的流体和介质特性,非线性渗流方程的具体形式也会有所不同。
在研究非线性渗流方程时,需要首先明确所研究的问题类型和边界条件,然后根据实际情况选择合适的方程形式。
三、非线性渗流方程的解析方法针对非线性渗流方程的解析方法,本文主要介绍以下几种:1. 线性化方法:将非线性渗流方程进行线性化处理,使其变为易于求解的线性方程。
这种方法适用于某些特定情况下的非线性渗流问题,如某些简单的多孔介质模型。
2. 有限元法:将求解区域划分为有限个单元,通过求解每个单元的近似解来得到整个区域的解。
有限元法具有较高的求解精度和灵活性,适用于复杂的多孔介质模型和边界条件。
3. 数值模拟方法:通过计算机模拟来求解非线性渗流方程。
这种方法可以处理复杂的非线性问题和多因素交互问题,但需要较高的计算资源和计算时间。
4. 反演分析法:根据实测数据反演求解非线性渗流方程的参数,以获得更准确的模型预测结果。
这种方法需要大量的实测数据和较高的数据处理能力。
四、非线性渗流方程的应用非线性渗流方程在地下水资源评价、油藏开发、环境污染控制等领域具有广泛的应用。
例如,在地下水资源评价中,可以通过求解非线性渗流方程来预测地下水的流动路径和流量,为合理开发利用地下水资源提供依据;在油藏开发中,可以通过求解非线性渗流方程来优化油藏开采方案,提高采收率;在环境污染控制中,可以通过求解非线性渗流方程来预测污染物的扩散范围和速度,为制定污染控制措施提供依据。
《非线性渗流方程解析方法研究及应用》篇一一、引言渗流现象广泛存在于自然界和工程领域中,如地下水流动、油藏开发等。
非线性渗流方程是描述这些现象的重要数学工具,其解析方法的研究对于理解渗流机制、优化工程设计和提高资源利用效率具有重要意义。
本文将重点研究非线性渗流方程的解析方法,并探讨其在实际应用中的价值。
二、非线性渗流方程概述非线性渗流方程是一类描述多孔介质中流体流动的偏微分方程,其形式复杂且具有高度非线性。
这类方程通常涉及到流体在多孔介质中的压力、饱和度、渗透率等参数,以及复杂的边界条件和初始条件。
由于非线性渗流方程的复杂性,传统的解析方法往往难以求解,因此需要发展新的解析方法。
三、非线性渗流方程的解析方法(一)渐近法渐近法是一种常用的非线性渗流方程解析方法。
该方法通过引入适当的近似和假设,将原问题简化为更易于处理的子问题。
在渐近法中,可以采用摄动法、分离变量法、尺度变换等方法来求解非线性渗流方程。
这些方法可以有效地求解一些具有特定边界条件和初始条件的问题,但在处理复杂问题时仍具有一定的局限性。
(二)数值解析法数值解析法是一种基于数值计算的非线性渗流方程解析方法。
该方法通过将原问题离散化,将连续的偏微分方程转化为离散的代数方程组,然后采用数值计算方法求解。
常见的数值解析法包括有限差分法、有限元法、边界元法等。
这些方法可以有效地求解复杂的非线性渗流方程,具有较高的精度和灵活性。
(三)其他解析方法除了渐近法和数值解析法外,还有一些其他的非线性渗流方程解析方法,如反演法、多尺度法等。
这些方法可以根据具体问题选择使用,以获得更准确的解。
四、非线性渗流方程的应用非线性渗流方程在地下水流动、油藏开发、污染物运移等领域具有广泛的应用。
通过解析方法求解非线性渗流方程,可以了解流体在多孔介质中的流动机制和分布规律,为工程设计和优化提供重要的依据。
例如,在油藏开发中,通过求解非线性渗流方程可以了解油气的流动规律和分布情况,为制定合理的开采方案提供依据;在地下水流动中,通过解析方法可以预测地下水的运动轨迹和水位变化情况,为地下水资源的保护和利用提供依据。
渗流-应力耦合作用对储层裂缝发育的影响研究宋子怡;王昊;李静;孙鲁宁;张加太;刘晨【摘要】裂缝是影响储层高产、稳产的重要因素,而储层处在复杂的地质环境中,裂缝的形成和发育受众多因素的影响,研究各因素间的耦合作用对储层裂缝发育的影响,对指导油气勘探开发具有重要意义.为此,针对任丘油田任11井区雾迷山组碳酸盐岩储层进行了渗流-应力耦合作用对储层裂缝发育的影响研究.研究结果表明:未考虑渗流-应力耦合作用时,研究区最大水平主应力范围为82~100 MPa,从西南到东北逐渐增大;最小水平主应力范围为72~88 MPa,从研究区中心向西南、东北两侧逐渐递增;考虑耦合作用后,研究区最大水平主应力范围为84~102 MPa,最小水平主应力范围为76~91 MPa,最大及最小水平主应力增加.渗流-应力耦合作用后,研究区裂缝发育指数分布在0.027~1.156之间,山头顶部和近东西向断层的内部区域裂缝发育指数在0.7左右,为裂缝较发育区域;而研究区西南和东北边缘区域裂缝发育指数在0.2以下,为裂缝欠发育区域.随着耦合作用时间的增长,储层裂缝发育指数逐渐增大,在注入井和产油井附近区域的裂缝发育指数增大幅度尤为明显;储层裂缝线密度也呈增大趋势,仅产油井周围的裂缝线密度呈现为先减小后增大的趋势;未考虑耦合作用时的储层裂缝参数小于考虑耦合作用后的裂缝参数,说明仅考虑应力场进行储层裂缝预测所得结果偏小.【期刊名称】《地质力学学报》【年(卷),期】2019(025)004【总页数】9页(P483-491)【关键词】储层裂缝;渗流;应力;耦合作用【作者】宋子怡;王昊;李静;孙鲁宁;张加太;刘晨【作者单位】中海油研究总院有限责任公司, 北京100028;中国石油大学 (华东) 地质力学与工程研究所, 山东青岛266580;中国石油大学 (华东) 地质力学与工程研究所, 山东青岛266580;中国石油大学 (华东) 地质力学与工程研究所, 山东青岛266580;中国石油大学 (华东) 地质力学与工程研究所, 山东青岛266580;中国石油大学 (华东) 地质力学与工程研究所, 山东青岛266580【正文语种】中文【中图分类】TE3190 引言储层裂缝是影响储层获得高产和稳产的重要因素,储层裂缝评价是油气勘探、开采的重要工作,也是油气地质学研究的重点和难点。
渗流力学渗流力学研究的内容流体通过多孔介质的流动称为渗流。
多孔介质是指由固体骨架和相互连通的孔隙、裂缝或各种类型毛细管所组成的材料。
渗流力学就是研究流体在多孔介质中运动规律的科学。
它是流体力学的一个重要分支,是流体力学与岩石力学、多孔介质理论、表面物理、物理化学以及生物学交叉渗透而形成的。
渗流现象普遍存在于自然界和人造材料中。
如地下水、热水和盐水的渗流;石油、天然气和煤层气的渗流;动物体内的血液微循环和微细支气管的渗流;植物体内水分、气体和糖分的输送;陶瓷、砖石、砂模、填充床等人造多孔材料中气体的渗流等。
渗流力学在很多应用科学和工程技术领域有着广泛的应用。
如土壤力学、地下水水文学、石油工程、地热工程、给水工程、环境工程、化工和微机械等等。
此外,在国防工业中,如航空航天工业中的发汗冷却、核废料的处理以及诸如防毒面罩的研制等都涉及渗流力学问题。
渗流的特点在于:(1)多孔介质单位体积孔隙的表面积比较大,表面作用明显。
任何时候都必须考虑粘性作用;(2)在地下渗流中往往压力较大,因而通常要考虑流体的压缩性;(3)孔道形状复杂、阻力大、毛管力作用较普遍,有时还要考虑分子力;(4)往往伴随有复杂的物理化学过程。
渗流力学是一门既有较长历史又年轻活跃的科学。
从Darcy定律的出现已过去一个半世纪。
20世纪石油工业的崛起极大地推动了渗流力学的发展。
随着相关科学技术的发展,如高性能计算机的出现,核磁共振、CT扫描成像以及其它先进试验方法用于渗流,又将渗流力学大大推进了一步。
近年来,随着非线性力学的发展,将分叉、混沌以及分形理论用于渗流,其它诸如格气模型的建立等等,更使渗流力学的发展进入一个全新的阶段。
渗流力学的应用范围越来越广,日益成为多种工程技术的理论基础。
由于多孔介质广泛存在于自然界、工程材料和人体与动植物体内,因而就渗流力学的应用范围而言,大致可划分为地下渗流、工程渗流和生物渗流3个方面。
地下渗流是指土壤、岩石和地表堆积物中流体的渗流。
变形多孔介质流固耦合非线性渗流模型及解析摘要:本文探讨了变形多孔介质流固耦合非线性渗流模型的特点,并提出了一种新的解析方法来解决该模型所表示的渗流问题。
文中首先介绍了变形多孔介质流固耦合非线性渗流模型的数学模型和算法,然后提出了一种用于解析该模型的解析方法,并分析了该方法的有效性。
最后,通过实际问题的计算验证,研究表明,该解析方法能够有效解决变形多孔介质流固耦合非线性渗流模型的计算问题。
变形多孔介质流固耦合非线性渗流是一种常见的多物理场耦合现象,它具有空间、时间多尺度和非线性等特点。
这种类型的渗流进程不仅会导致渗流发生变化,而且还可能影响到包括多物理场在内的其他特性。
因此,研究变形多孔介质流固耦合非线性渗流模型是多物理场耦合问题的重要研究内容。
变形多孔介质流固耦合非线性渗流模型的数学模型由多个多物理场的偏微分方程组构成,其中包括流体流动的变形Navier-Stokes 方程、热传导的传热方程、物质渗透的物质输运方程等。
此外,模型还具有多个参量(如流体非线性参量、渗流系数、温度变化、介质密度、介质体积等)。
进一步地,由于变形多孔介质流固耦合非线性渗流模型是一个多尺度耦合系统,因此其计算难度较高。
在传统的基于数值方法的研究中,科学家们使用离散化方法解决了变形多孔介质流固耦合非线性渗流模型的计算问题,但是存在系数稳定性的问题以及计算成本的问题。
为了解决这些问题,本文提出了一种新的解析方法来解决变形多孔介质流固耦合非线性渗流模型。
该解析方法可以有效解决变形多孔介质流固耦合非线性渗流模型的计算问题。
该解析方法首先利用蒙特卡洛方法获取流体、热传导和物质输运的局部解,然后将局部解用于多物理场耦合动力学系统的本征分析,最终得到全局的解。
该方法不仅具有计算精度高、计算成本低的优点,而且还可以有效提取变形多孔介质流固耦合非线性渗流模型的宏观行为特征,进一步分析和推理问题的性质。
本文中,通过对实际问题数据进行计算验证,研究表明,该解析方法能够有效解决变形多孔介质流固耦合非线性渗流模型的计算问题。
FLAC3DFLAC3D(Fast Lagrangian Analysis of Continua快速拉格朗日分析)由美国Itasca 公司开发的。
目前,FLAC有二维和三维计算程序两个版本,二维计算程序V3.0以前的为DOS 版本,V2.5版本仅仅能够使用计算机的基本内存(64K),所以,程序求解的最大结点数仅限于2000个以内。
1995年,FLAC2D已升级为V3.3的版本,其程序能够使用护展内存。
因此,大大发护展了计算规模。
FLAC3D是一个三维有限差分程序,目前已发展到V2。
1版本。
FLAC3D是美国Itasca Consulting Group Inc开发的三维显式有限差分法程序,它可以模拟岩土或其它材料的三维力学行为。
FLAC软件的基本原理是拉格朗日差分法(源于流体力学,主要研究流体质点在任一时段内的运动轨迹、速度、压力等特征,将其移植到固体力学中,把所研究的区域划分成网格,其结点相当于流体质点,然后按时步来用拉格朗日法来研究网格结点的运动),它是一种利用拖带坐标系分析大变形问题的数值方法,并利用差分格式按时步积分求解。
随着构形的不断变化,不断更新坐标,允许介质有较大的变形。
模型经过网格剖分,物理网格映射成数学网格,数学网格上的某个结点就与物理网格上相应的结点坐标相对应。
对于某一个结点而言,在每一时刻它受到来自其周围区域的合力的影响。
如果合力不等于零,结点具有了失稳力(Unbalanced Force) ,就要产生运动。
假定结点上集中有连接该结点的质量,于是在失稳力作用下,根据牛顿定律,结点就要产生加速度,进而可以在一个时步中求得速度和位移的增量。
对于每一个区域而言,可以根据其周围结点的运动速度求得它的应变率,然后根据材料的本构关系求得应力的增量。
由应力增量求出t和t+Δt时刻各个结点的不平衡力和各个节点在t+Δt时的加速度。
对加速度进行积分,即可得出结点的新的位移值,由此可求得各结点新的坐标值。
1求解步骤
1.1模块挂载
(1)启动SciFEA,选择“工具”—>“模块挂载”菜单弹出如图1-1所
示的对话框,在对话框当中选择“三维非线性温度渗流耦合”。
图1-1 选择项目类型对话框
(2)点击“SciFEA模块挂载”栏中的“三维非线性温度渗流耦合”选项。
如图1-1所示。
(3)点击“确定”按钮完成项目类型的选择。
(4)单击工具条中的按钮弹出如图1-2所示,找到SciFEA下的tempudm 目录,并且选择打开udm.pro文件。
图1-2
(5)点击“打开”按钮之后,弹出如图1-3:
图1-3
(6)最后点击“OK”按钮。
1.2设置材料参数和边界条件
(1)单击工具条中的按钮弹出如图1-4所示材料参数数据输入表格。
图1-4 材料参数输入
(2)单击工具条中的按钮弹数据输入表格。
按照问题描述中的参数依次填入边界条件数据表格。
填写完成后如图1-5所示。
图1-5填写完成扩散边界
2建模、设置材料属性和施加边界条件
(1)启动GID以创建模型。
单击工具条中的按钮弹出如图1-6所示前
处理初始化窗口。
图1-6前处理初始化窗口
2.1建模。
a. 点击菜单【Geometry】-【Create】-【Line】,在命令栏依次输入点坐标,最后得到模型图如1-7所示:
图1-7建成的模型图
b.定义材料特性、施加载荷约束。
点击菜单【Data】-【Conditions】弹出“Conditions”对话框。
①施加载荷。
点击线单元按钮,选中下拉菜单中的“surface-htcn3da”,最后得
到如图1-8所示。
图1-8 施加载荷
②施加载荷。
点击线单元按钮,选中下拉菜单中的“surface-htcn3db”,最后得
到如图1-9所示。
图1-9
③施加载荷。
点击线单元按钮,选中下拉菜单中的“volume-D1htcn3db”,最后得到如图1-10所示。
图1-10
④施加载荷。
点击线单元按钮,选中下拉菜单中的“volume-D1htcn3dc”,最后得到如图1-11所示。
图1-11
⑤施加载荷。
点击线单元按钮,选中下拉菜单中的“volume-aec8”,最后得到
如图1-12所示。
图1-12
需要说明的是,这几步定义材料,施加载荷的操作没有先后顺序要求,但必须保证不能遗漏某个条件的定义。
2.2划分网格和导出数据
(1).划分网格。
点击菜单【Meshing】-【Generate】,要求定义单元尺寸大小,
点击“OK”按钮,此时弹出“Progress in meshing”进度条,网格划分完毕后弹出消息框显示总的单元数和节点数,如图1-13所示。
图1-13 划分网格得到的单元节点信息
点击“OK”按钮后即可得到划分完的网格模型,如图1-14
图1-14网格图
(1)把前处理数据转化为SciFEA所需计算格式。
点击【Files】-【Save】
菜单,保存为GID文件。
点击菜单【Calculate】-【Calculate】,弹出
“process info”消息框,如图1-15所示,点击“ok”按钮,然后退出
Gid,至此前处理工作结束。
图1-15 转化数据消息框
3有限元计算
(1) 单击工具条中的按钮弹出如图1-16所示计算模拟窗口。
图1-16 计算模拟窗口
4结果分析
对计算结果进行分析属于后处理,GID提供了丰富的后处理操作,可以从不同角度,不同方式来输出计算结果。
(1)点击工具菜单中的进入后处理,如图1-17所示。
图1-17 进入后处理结果分析
(2)点击菜单【Files】-【Postprocess】进入后处理程序。
(3)点击菜单【View results】-【Contour fill】-【unoda0】显示的位移云图,如图1-18所示。
图1-18。