CN113609729B - 一种东亚地区地幔动力学的数值模拟方法 - Google Patents
一种东亚地区地幔动力学的数值模拟方法 Download PDFInfo
- Publication number
- CN113609729B CN113609729B CN202110868514.9A CN202110868514A CN113609729B CN 113609729 B CN113609729 B CN 113609729B CN 202110868514 A CN202110868514 A CN 202110868514A CN 113609729 B CN113609729 B CN 113609729B
- Authority
- CN
- China
- Prior art keywords
- phase change
- calculation module
- calculation formula
- viscosity
- component
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000004088 simulation Methods 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 46
- 230000008859 change Effects 0.000 claims abstract description 23
- 239000000463 material Substances 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 12
- 239000000126 substance Substances 0.000 claims abstract description 11
- 239000012530 fluid Substances 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims abstract description 4
- 238000009792 diffusion process Methods 0.000 claims description 10
- 230000004913 activation Effects 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 5
- 238000004134 energy conservation Methods 0.000 claims description 3
- 230000020169 heat generation Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 description 5
- 230000009189 diving Effects 0.000 description 4
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种东亚地区地幔动力学的数值模拟方法,包括以下步骤:化学组分计算模块调用不同组分场的参数,计算不同组分的影响以及其运动;相变计算模块调用相变参数,计算组分的相变过程以及其对流场的影响;并行算法模块根据需要启动并行计算节点、协调计算节点间的数据调度,并根据初边值条件设置全局或节点的信息;流体运动核心计算模块根据有限元节点信息建立矩阵系统,并根据配置文件确定是否调用化学组分计算模块和相变计算模块,最终计算物质速度与温度。本发明中的模型结合自适应网格方法,可包含与应变率相关的流变特性,对深部动力过程的模拟更加真实地球,结果更加合理。
Description
技术领域
本发明属于地球动力学数值模拟技术领域,特别涉及一种东亚地区地幔动力学的数值模拟方法。
背景技术
现有的东亚区域的俯冲过程及深部流场的数值模拟未包含与应变率相关的流变特性,这是由于算法的复杂性和计算时的巨大耗时造成的;
现有的地幔对流数值模拟程序模拟的介质模型多为二维平面;
现有的地幔对流数值模拟程序模拟多进行的是瞬态问题的研究;
现有的地幔对流数值模拟程序多只涉及单组分在地震波速不连续面的相变。
发明内容
为了至少解决或部分解决上述问题,提供一种东亚地区地幔动力学的数值模拟方法。
为了达到上述目的,本发明提供了如下的技术方案:
本发明一种东亚地区地幔动力学的数值模拟方法,包括以下步骤:
设置初边值条件,设置动力学参数;
化学组分计算模块调用不同组分场的参数,计算不同组分的影响以及其运动;相变计算模块调用相变参数,计算组分的相变过程以及其对流场的影响;
并行算法模块根据需要启动并行计算节点、协调计算节点间的数据调度,并根据初边值条件设置全局或节点的信息;
流体运动核心计算模块根据有限元节点信息建立矩阵系统,并根据配置文件确定是否调用化学组分计算模块和相变计算模块,最终计算物质速度与温度。
作为本发明的一种优选技术方案,所述流体运动核心计算模块的最终计算公式为
其中,代表应变率,p为压力,u为速度。
其中不可压缩条件下上述两个方程分别为:
能量守恒方程为
其中,Cp为热容,T为温度,k为热导率,方程右端为热源项,代表潜热生热,S为熵,X代表某材料的比例。
作为本发明的一种优选技术方案,所述流体运动核心计算模块的最终计算公式还包括组分平流方程,所述组分平流方程的计算公式为
其中,ci代表第i个组分值,qi为该组分的生成或减少值;
忽略化学组分的影响,由岩石的热膨胀性产生的密度异常的计算公式为
Δρ=-αρ0(T-T0)
其中,α,p0和T0分别为热膨胀系数、参考密度和参考温度;
全应力张量计算公式为
σij=τij-pδij
其中,τij,p和δij分别表示偏应力、动压力和Kroneckerδ函数。
偏应力张量计算公式为
其中μ为黏度。
作为本发明的一种优选技术方案,所述黏度的处理方法为:
计算扩散蠕变和位错蠕变;
计算屈服黏度;
根据屈服黏度、扩散蠕变和位错蠕变得到黏度。
作为本发明的一种优选技术方案,所述扩散蠕变和位错蠕变的计算公式为
其中,扩散蠕变μdiff中,取n=1,位错蠕变μdisl中,取n=3.5;A为参数因子,E为活化能,V为活化体积,P为压力,R为气体常数,为应变率张量的第二不变量,Tr为温度。
作为本发明的一种优选技术方案,所述屈服黏度的计算公式为
其中,τyield=min(τ0+fcP,τy,max),τ0为表面屈服强度,fc为摩擦系数,τy,max为最大屈服强度。
作为本发明的一种优选技术方案,所述黏度的计算公式为
与现有技术相比,本发明的有益效果如下:
本发明中的模型结合自适应网格方法,可包含与应变率相关的流变特性,对深部动力过程的模拟更加真实地球,结果更加合理。
本发明可以模拟二维和三维模型。
本发明既可以进行瞬态的模拟,也可以进行演化的模拟。
本发明既可以进行多组分的相变模拟,也可以进行某一组分在不同深度发生多次相变的模拟。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。在附图中:
图1是本发明的整体方法流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
此外,如果已知技术的详细描述对于示出本发明的特征是不必要的,则将其省略。
实施例1
如图1所示,本发明提供一种东亚地区地幔动力学的数值模拟方法,包括以下步骤:
设置初边值条件,设置动力学参数;
化学组分计算模块调用不同组分场的参数,计算不同组分的影响以及其运动;相变计算模块调用相变参数,计算组分的相变过程以及其对流场的影响;
并行算法模块根据需要启动并行计算节点、协调计算节点间的数据调度,并根据初边值条件设置全局或节点的信息;
流体运动核心计算模块根据有限元节点信息建立矩阵系统,并根据配置文件确定是否调用化学组分计算模块和相变计算模块,最终计算物质速度与温度。
对本发明中所使用的微分方程进行说明:
其中,代表应变率,p为压力,u为速度。
其中不可压缩条件下上述两个方程分别为:
能量守恒方程为
其中,Cp为热容,T为温度,k为热导率,方程右端为热源项,代表潜热生热,S为熵,X代表某材料的比例。
联立上述微分方程,可以得到速度场和温度场。除此之外,还包括组分平流方程
其中,ci代表第i个组分值,qi为该组分的生成/减少值。
忽略化学组分的影响,由岩石的热膨胀性产生的密度异常可以写为
Δρ=-αρ0(T-T0)
其中,α,ρ0和T0分别为热膨胀系数、参考密度和参考温度。
全应力张量可以表示为:
σij=τij-pδij
其中,τij,p和δij分别表示偏应力、动压力和Kroneckerδ函数。
偏应力张量为:
其中μ为黏度。
黏度的处理方法为:
(1)扩散蠕变和位错蠕变
其中,扩散蠕变μdiff中,取n=1,位错蠕变μdisl中,取n=3.5。A为参数因子,E为活化能,V为活化体积,P为压力,R为气体常数,为应变率张量的第二不变量,Tr为温度。
(2)屈服黏度
其中,τyield=min(τ0+fcP,τy,max),τ0为表面屈服强度,fc为摩擦系数,τy,max为最大屈服强度。
(3)黏度
下面是根据上述方法模拟的实施程序:该地幔对流数值模拟系统程序可以用来模拟全球和局部地区的地幔对流与板块俯冲过程。对于全球尺度的深部地幔动力学的研究,可以配合Gplates进行三维空间的模拟,得到不同深度上的流场模式。对于局部地区的俯冲动力学的研究,可针对研究区域的地质和地震学观测结果,构建区域模型的初条件。在原有程序基础上改进的材料模型插件可以对选用多种组分处理粘塑性相变问题。
(一)全球模型模拟实例
将需要建立的模型的各种条件写入.prm文件:
(1)三维全球几何模型,模型包括岩石圈和地幔部分:
(2)设置初边值条件:
模型设置中的动力学参数为:
(3)网格设置:
(4)规定计算时间及输出目录后,提交.prm文件。计算结果可由Paraview打开。
由于计算的是三维模型,根据研究需求,可以显示出不同计算时间下的不同深度,不同剖面,不同方向的速度场。
(二)俯冲相变模拟实例
本系统能够模拟:材料组分的速度、温度,并通过后处理程序得到应变率、密度、粘度场等数据;通过更改物性参数文件,可以模拟不同的材料介质。实验构建了如下二维实验模型,该模型长4000km,高1000Km。
该二维模型的介质物性参数取值入下表所示:
该模型的初始、边界条件设置如下:
(1)温度初始条件:岩石圈部分遵循半空间冷却模型,地幔部分遵循绝热剖面;温度边界条件:上下界面固定,上界面为273K,下界面温度遵循初始温度设置。
(2)速度边界条件:各边界均为自由滑动边界条件。
编译插件程序,生成相应.so文件。
在参数输入prm文件中写入相应的库文件名,并设置相应计算参数。
(3)设置并行计算节点数,以aspect*.prm命令的形式进行计算。
(4)后处理及可视化:将生成文件中的solution.pvd文件用paraview软件打开进行处理;可以显示出不同计算时间下的不同深度,不同剖面,不同方向的速度场。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种东亚地区地幔动力学的数值模拟方法,其特征在于,包括以下步骤:
设置初边值条件,设置动力学参数;
化学组分计算模块调用不同组分场的参数,计算不同组分的影响以及其运动;相变计算模块调用相变参数,计算组分的相变过程以及其对流场的影响;
并行算法模块根据需要启动并行计算节点、协调计算节点间的数据调度,并根据初边值条件设置全局或节点的信息;
流体运动核心计算模块根据有限元节点信息建立矩阵系统,并根据配置文件确定是否调用化学组分计算模块和相变计算模块,最终计算物质速度与温度;
所述流体运动核心计算模块的最终计算公式为
其中,代表应变率,p为压力,u为速度;
其中不可压缩条件下上述两个方程分别为:
能量守恒方程为
其中,Cp为热容,T为温度,k为热导率,方程右端为热源项,代表潜热生热,S为熵,X代表某材料的比例;
所述流体运动核心计算模块的最终计算公式还包括组分平流方程,所述组分平流方程的计算公式为
其中,ci代表第i个组分值,qi为该组分的生成或减少值;
忽略化学组分的影响,由岩石的热膨胀性产生的密度异常的计算公式为
Δρ=-αρ0(T-T0)
其中,α,ρ0和T0分别为热膨胀系数、参考密度和参考温度;
全应力张量计算公式为
σij=τij-pδij
其中,σij,p和δij分别表示偏应力、动压力和Kroneckerδ函数;
偏应力张量计算公式为
其中μ为黏度;
所述黏度的处理方法为:
计算扩散蠕变和位错蠕变;
计算屈服黏度;
根据屈服黏度、扩散蠕变和位错蠕变得到黏度;
所述扩散蠕变和位错蠕变的计算公式为
其中,扩散蠕变μdiff中,取n=1,位错蠕变μdisl中,取n=3.5;A为参数因子,E为活化能,V为活化体积,P为压力,R为气体常数,为应变率张量的第二不变量,Tr为温度;
所述屈服黏度的计算公式为
其中,τyield=min(τ0+fcP,τy,max),τ0为表面屈服强度,fc为摩擦系数,τy,max为最大屈服强度;
所述黏度的计算公式为
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110868514.9A CN113609729B (zh) | 2021-07-30 | 2021-07-30 | 一种东亚地区地幔动力学的数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110868514.9A CN113609729B (zh) | 2021-07-30 | 2021-07-30 | 一种东亚地区地幔动力学的数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113609729A CN113609729A (zh) | 2021-11-05 |
CN113609729B true CN113609729B (zh) | 2023-09-01 |
Family
ID=78306130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110868514.9A Active CN113609729B (zh) | 2021-07-30 | 2021-07-30 | 一种东亚地区地幔动力学的数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113609729B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005049477A1 (de) * | 2005-10-13 | 2007-04-19 | Freiberger Compound Materials Gmbh | Verfahren und Vorrichtung zur Kristallzüchtung mittels kombinierter Heizerregelung |
JP2015006650A (ja) * | 2013-06-26 | 2015-01-15 | 須知 晃一 | システム構成構造細胞複合諸物体の製造方法と構成材料 |
CN105092699A (zh) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | 岩石超声波测试系统及高温高压三分量超声探头制作方法 |
CN112364577A (zh) * | 2020-11-30 | 2021-02-12 | 季颖锋 | 一种板块俯冲热结构参数的数值模拟分析方法 |
CN113009559A (zh) * | 2021-03-02 | 2021-06-22 | 中国科学院地质与地球物理研究所 | 基于多类型地球物理数据的地震评价方法 |
-
2021
- 2021-07-30 CN CN202110868514.9A patent/CN113609729B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005049477A1 (de) * | 2005-10-13 | 2007-04-19 | Freiberger Compound Materials Gmbh | Verfahren und Vorrichtung zur Kristallzüchtung mittels kombinierter Heizerregelung |
JP2015006650A (ja) * | 2013-06-26 | 2015-01-15 | 須知 晃一 | システム構成構造細胞複合諸物体の製造方法と構成材料 |
CN105092699A (zh) * | 2014-05-20 | 2015-11-25 | 中国石油化工股份有限公司 | 岩石超声波测试系统及高温高压三分量超声探头制作方法 |
CN112364577A (zh) * | 2020-11-30 | 2021-02-12 | 季颖锋 | 一种板块俯冲热结构参数的数值模拟分析方法 |
CN113009559A (zh) * | 2021-03-02 | 2021-06-22 | 中国科学院地质与地球物理研究所 | 基于多类型地球物理数据的地震评价方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113609729A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | A three-dimensional thermo-hydro-mechanical coupled model for enhanced geothermal systems (EGS) embedded with discrete fracture networks | |
Mayer et al. | 3D fluid–structure-contact interaction based on a combined XFEM FSI and dual mortar contact approach | |
CN110298105A (zh) | 饱和多孔介质大变形分析的ccpdi-impm方法 | |
CN113609729B (zh) | 一种东亚地区地幔动力学的数值模拟方法 | |
Yang et al. | A simple gas kinetic scheme for simulation of 3D incompressible thermal flows | |
Yu et al. | A 3D Coupled Thermal-Hydraulic-Mechanical THM Model Using EDFM and XFEM for Hydraulic-Fracture-Dominated Geothermal Reservoirs | |
Ibrahimbegovic et al. | Fluid-structure interaction problems solution by operator split methods and efficient software development by code-coupling | |
Yu et al. | A lattice Boltzmann model for reaction dynamical systems with time delay | |
Jordan et al. | Parallel interactive visualization of 3d mantle convection | |
Ruhi et al. | A lattice Boltzmann relaxation scheme for incompressible fluid flows | |
Kim et al. | Visual simulation of turbulent fluids using mls interpolation profiles | |
Vollmer et al. | Revisiting the Boundary Conditions for Unsteady Flows Adjacent to Rigid and Dynamic Solid Walls | |
Das | An Introduction to Finite Element Analysis Using Matlab Tools | |
Profir | Mesh morphing techniques in CFD | |
Wang et al. | A prediction of the acoustical properties of induction cookers based on an FVM–LES-acoustic analogy method | |
Kawai | ADVENTURE AutoGL: A handy graphics and GUI library for researchers and developers of numerical simulations | |
Van Treeck et al. | Direct building energy simulation based on large eddy techniques and lattice Boltzmann methods | |
Talebian et al. | An extended finite element model for CO2 sequestration | |
Boylston | Quasi-One-Dimensional Flow for Use in Real-Time Facility Simulations | |
Ordonez Egas | Scalable linear solver for thermo-hydro-mechanics with a second gradient of dilation regularization problems | |
Gurugubelli et al. | Modeling of Fluid-Structure Interactions with Exact Interface Tracking Methods | |
Heidarinejad | Vortex simulation of the reacting shear layer | |
Hossain | Investigating the Nature of Compressible Boundary Layers at High Reynolds Number | |
Jia et al. | Sensitivity analysis of heat extraction performance and a data-mining based BP neural network prediction model of geothermal doublet system | |
Nguyen et al. | AMass Conservation AlgorithmFor Adaptive Unrefinement Meshes Used By Finite Element Methods |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |