基于随机生长四参数法的多孔介质模型.pdf
《基于随机生长四参数法的多孔介质模型.pdf》由会员分享,可在线阅读,更多相关《基于随机生长四参数法的多孔介质模型.pdf(7页珍藏版)》请在咨信网上搜索。
1、能源研究与管理2023,15(2)能源研究与管理2023,15(2)收稿日期:2023-02-20基金项目:核技术应用教育部工程研究中心开放基金(HJSYB2022-1)第一作者:王召(1997),女,硕士研究生,主要从事物理与数值模拟研究。E-mail:*通信作者:杜俊杰(1987),男,讲师,博士,主要从事超导磁体及常规磁体技术研究。E-mail:引用格式:王召,杜俊杰,刘俊雄.基于随机生长四参数法的多孔介质模型J.能源研究与管理,2023,15(2):88-94.摘要:为了研究所列因素的改变对生成孔隙空间以及程序运行时间的影响,使用向量乘法建立四参数随机生长方法(QSGS)模型,并讨论模
2、型参数对生成孔隙结构及计算机生成孔隙模型耗时的影响,使用Matlab软件针对不同孔隙度、固相核分布和生长概率参数建立二维和三维孔隙模型。结果表明:孔隙率与固相核占比呈正相关,孔隙率选取过大会生成无效孔隙,过小会导致固相核孤立。在二维和三维模型中,模孔隙率参数越大生成孔隙模型时间越长,而生长概率与孔隙分布参数改变对模型生成时间影响较小。在三维模型情况下由于消耗更多计算资源,因此研究计算区域网格数对孔隙模型生成时间影响,且随网格数增加模型生成时间显著增加。关键词:随机生长四参数法;多孔介质;土壤;孔隙模型中图分类号:X502文献标志码:A文章编号:20967705(2023)02008807WAN
3、G Zhao1,2,DU Junjie1,*,LIU Junxiong2(1.Engineering Research Center of Ministry of Education for Nuclear Technology Application,East China Institute of Technology,Nanchang 330013,China;2.School of Nuclear Science andEngineering,East China Institute of Technology,Nanchang 330013,China)To investigate t
4、he impact of changes in the listed factors on the generation of pore space and program runtime.Multiplication of vectors is used to the establish of four-parameter stochastic growth method(QSGS)model.The influence ofmodel parameters on the generation of pore structure and the time consumption of com
5、puter generated pore model isdiscussed.Using Matlab software,two-dimensional and three-dimensional pore models are established for different porosity,solid core distribution and growth probability parameters.Numerical simulation and analysis are carried out for differentporosity,distribution probabi
6、lity of solid phase nuclei and growth probability,and the changes of program running timecaused by parameter changes are counted.In two dimensions,the porosity is positively correlated with the proportion of thesolid phase.High porosity generates invalid pores and low porosity leads to the isolation
7、 of solid phase nuclei.The smaller isthe distribution probability,the fewer is the number of initial solid phase nuclei.The growth probability affects the poredistribution in precision.In three dimensions,the effect of different parameters on the pore distribution is similar to that of 2D,but the ru
8、nning time is significantly increased.In addition,it is found that in three dimensions the running time increasesDOI:10.16056/j.2096-7705.2023.02.013基于随机生长四参数法的多孔介质模型王召1,2,杜俊杰1,*,刘俊雄2(1.东华理工大学核技术应用教育部工程研究中心,南昌330013;2.东华理工大学核科学与工程学院,南昌330013)研究与探讨88窑窑能源研究与管理2023,15(2)能源研究与管理2023,15(2)significantly w
9、hen the number of grids increases.quartet structure generation set(QSGS);porous medium;soil;heavy metal speciation引言土壤孔隙指的是土壤颗粒之间或聚团物之间存在的间隙。土壤的孔隙作为容纳水、空气和溶质的空间,对于提高作物产量以及治理环境污染等问题起到非常重要的作用。例如在土壤孔隙中含水、空气和溶质容量和种类不同会直接影响到土壤里的动物以及微生物的数量以及质量,对农作物产量产生不同程度影响。同时在我国重金属污染土壤的环境问题研究中1-4,张荣辉等5通过高氧官能团的多孔碳材料为铁氧化物
10、结合土壤中的重金属离子提供了有效的表面配体,可以协同铁氧化物修复重金属污染土。因此孔隙结构对于相关问题的研究具有重要的意义。但在理论计算孔隙流,水迁移和污染物扩散等问题中,由于孔隙结构复杂往往不能使用传统方法构造孔隙通道,并且随着计算机技术不断更新发展,因此使用计算机生成孔隙结构成为具有成本更低、效率更高且更直观方法之一。传统的测量土壤孔隙率的方法有直接测出密度和比重来计算及压汞法6,但这2种方法对于研究孔隙的具体形态效果较差。随着图像处理技术的快速发展,土壤切片技术7、X射线断层扫描8-10以及核磁共振6等技术能更为直观地观测孔隙,但也有技术复杂或有成本高等缺陷。简化孔隙结构建立模型,通过计
11、算机进行模拟分析,既能解决成本问题,又能解决其他方法样本量小代表性不强等问题是很好的选择。选择更优的建模方法,建立更准确实用的模型成为了新课题。首先,由WANG等11提出了可人为控制的多孔介质微观结构重构方法,即四参数随机生长方法(quartet structure generation set,QSGS)。与其他随机生长算法相比,QSGS融入更多物理意义明确的参数,使得多相之间的接触连接问题得到了较好解决,同时颗粒分布更加随机,统计特性能够得到更效提取。然而该模型中颗粒尺寸单一,与自然多孔介质中颗粒尺寸多变特性相比仍存在一定差别,不能够很好地表征多孔介质的尺度不变属性。张季如等12对QSGS
12、方法的微观孔隙结构的分析特征进行了研究,为参数的合理化取值提供参考依据。罗齐彬等13及张朔等14利用QSGS方法建立二维微观孔隙模型,分别使用有限元和格子玻尔兹曼方法(lattice Boltzmann method,LBM)对孔隙内流体进行模拟,得出了空隙度、缩放覆盖率等变化对孔隙分布的影响结果。李滔等15使用QSGS方法建立的三维孔隙模型利用LBM模型对三维孔隙结构模拟并对孔隙流体进行研究。基于前人的研究以及深入的讨论,本文利用Matlab来实现QSGS方法,选取孔隙度、固相核分布概率s、生长概率g与网格数作为变量,建立二维和三维孔隙模型,进一步利用LBM模型对孔隙流体的数值模拟。1模型建
13、立1.1二维随机生长模拟为得到最终的多孔介质模型,QSGS算法将土壤基质分为固相以及非固相(孔隙)。在固定区域内首先随机构建固相生长内核,再在选择方向上让此内核按照给定生长概率连续随机生长。土的团聚体由细小的矿物颗粒经凝聚和胶结作用而形成3,16-17,是一个渐进的过程,而QSGS算法的固相生长法则类似于土的团聚体形成过程18-19。以土粒和孔隙二相构成的二维土体孔隙结构为例,若令土粒为生长相,孔隙为非生长相,构造区域的初始相全为孔隙,则构造过程如下:1)在构造区内的每个网格节点上,以分布概率d随机分布固相生长核,在最终构造出的孔隙结构中,d不能大于土粒面积的百分比,即d臆1;2)网格节点上的
14、各生长核按照给定的方向生长概率(代表生长方向,=1,2,8),二维生长方向示意图如图1所示,从图中所示的四面八方向周围相邻网格节点生长。若8个方向的生长概率18相同,则生成各向同性的孔隙结构;3)重复步骤2)直至孔隙面积达到给定的孔隙度。应该指出的是,QSGS算法通过、d、(,)4个参数来影响和控制多相多孔介质微观结构的生成,其中(,)为概率密度,表示在方向上第相在第相上的生长概率,反映了各相之间的相互影响。对于土粒和孔隙二相构成的土体,只需采用、d和3个参数,即可获得表征土体孔隙结构的QSGS模型20-22。研究与探讨89窑窑能源研究与管理2023,15(2)能源研究与管理2023,15(2
15、)由图1可知,在二维情况下只有8个生长方向,从正上方开始沿顺时针方向的标号分别为18,为了计算方便,将生成随机数小于固体核生成概率的点定义为(,),沿18方向对应的位置坐标行向量为0 1 1 1 0111、列向量1 1 01110 1。由此将每次生存随机数小于生长方向概率向量坐标定位1,其他为0,这样通过简单的对应相乘便可得到生长方向的坐标。可以发现整个过程类似于随机行走,模型在这种移动模型中,移动节点随机选择一个方向和速度来从当前位置移动到新的位置。而随机行走模型也被称为布朗运动。使用以上方法另外一个难点是在固体核生长到边界时由于边界外未定义的原因使得代码不能继续进行,因此根据4边、4顶点的
16、限制条件设定了生长边界,并且假设生长到边界时不允许“逆”生长。1.2三维随机生长模型随着维度的增加生长方向成指数式增长,正方体的共有26个生长方向如图2所示,左图为其中14个记为1,右图为其余12个记为2。其中设定了6面、8顶点和12边的限制生长条件。通过点、线、面及逆时针和由上至下的命名规则对26个生长方分别向进行标号其中18为顶点、920为边、2126为面,其对应向量分别为行向量()1 1111 1 1111 1 1 1 01 0 0 1 011 0 1 0 00、列向量()1 1 1111 1 1 111 1 0 1011 0 1 0 01 0 1 0 0、纵向量()1 1 1 1111
17、1 0 0 0 0 1 1 1 11111 0 0 0 0 11。同样假设在固相核生长到边界时不允许“逆”生长。图 1二维生长方向示意图Fig.12D growth direction diagram初始生成的随机数分布,通过与分布概率的大小比较来确定固相核的位置,然后对应每个固相核生成26个生长方向的向量,则构成了初始位置矩阵,同样与三维向量、简单的对应相乘随后便可得到随机生长方向的位置坐标。整个生长过程较为简单,依然与布朗运动情况类似,而在处理边界固相核生长时较为复杂。循环以上过程直到孔隙率达到要求为止,便可构建所需参数下的孔隙分布模型。2模拟结果与讨论2.1二维随机生长法模拟结果选用Ma
18、tlab采用本文1.1部分的具体方法建立循环,成立模型。在这一模型下影响最终成像结果的变量为以下所列,通过控制变量法分别观测不同变量的变化对模拟结果的影响,以及导致的计算时间变化。孔隙度由孔隙体积占基质总体积的百分比来表示。指基质中通气孔隙与持水孔隙的总和。孔隙度可分为有效孔隙度和无效孔隙度。有效孔隙度是指那些互相连通的,在一般压力条件下,可允许流体在其中流动的孔隙体积之和与基质总体积的比值,与之相对的则是无效孔隙度。在考虑土壤实际应用价值以及污染治理时,被考虑较多的为有效孔隙度。生长概率即固相内核在选定方向生长的概率,在实际中可以选择各向同性也可选择各向异性。固相核分布概率会影响所选区域内固
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 随机 生长 参数 多孔 介质 模型
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。