申请/专利权人：浙江大学

申请日：2022-04-29

公开（公告）日：2024-06-07

公开（公告）号：CN114781285B

主分类号：G06F30/28

分类号：G06F30/28;G06F113/08;G06F119/08;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.06.07#授权;2022.08.09#实质审查的生效;2022.07.22#公开

摘要：本发明公开了一种基于球簇假设和Laguerre‑Voronoi结构的生物质大颗粒热解仿真方法。所述方法属于生物质热解领域，可用于探究生物质的热解行为和现象，包括初始球簇生成、时间步离散及热解迭代、结果输出三个主要步骤。本发明利用球簇假设表征不规则形状生物质颗粒；利用组成子球的体积减小和移动实现母颗粒的收缩；利用Laguerre‑Voronoi结构划分颗粒内部区域，不仅实现了颗粒内部热传输过程，还解析了颗粒内部孔隙结构；利用方向相关的热导率和收缩因子实现了颗粒各向异性的热传导和收缩过程。本方法有助于获取生物质热解过程中实验难以观察得到的内部细节，也可以扩展用于含多生物质颗粒的反应器尺度计算流体力学模型中，进而大幅度提高大尺度模拟的预测精度和准确度。

主权项：1.一种基于球簇假设和Laguerre-Voronoi结构的生物质大颗粒热解仿真方法，其特征在于，包括以下步骤：S1，初始球簇生成采用具有大小、密度、温度、物质组成属性的子球体模拟真实生物质材料，由若干子球体共同构成球簇，以球簇等效表示原始生物质大颗粒；S2，在时间尺度上将热解过程划分为若干时间区间，并在每一时间区间上进行颗粒热解迭代；颗粒热解迭代过程包括颗粒结构更新、传热计算、热解反应计算和颗粒收缩计算四个子步骤；所述的颗粒结构更新包括外包络面更新和内Laguerre-Voronoi结构更新；其中，外包络面更新具体方法为：采用若干与球簇表面子球体相切的平面将球簇整体包络在一凸多面体内，所形成的包络面集合构成所表征生物质颗粒的边界；内Laguerre-Voronoi结构更新具体方法为：利用Laguerre-Voronoi结构对任意子球体所属的空间区域进行划分：中心位置为xi，半径为ri的子球i所属的空间区域定义为Laguerre-Voronoi单元，其中任意一点x满足 对任意子球jj≠i恒成立，其中xj，rj分别为子球j的中心位置和半径，dx,xi和dx,xj表示点x与子球i和j中心的距离；所述的颗粒收缩计算包括子球体积更新和子球平移计算；其中，子球体积更新由热解过程中的质量变化和密度变化实现，其公式为 其中Vi，mi，ρi分别表示子球i的体积、质量和密度；子球平移计算包括主动平移子球计算和被动平移子球计算；所述的主动平移子球计算方法为：首先定义收缩阈值γthreshold，其值在[0,1]范围内选取；随后对子球体进行穷举，满足Ωi＞γthreshold·maxΩ关系式的子球均发生主动移动，相反则不具备主动移速，其中Ωi表示子球i的单位体积固体消耗速率，maxΩ则表示当前时刻最大的单位体积固体消耗速率；子球主动移速基于在时间步内对应Laguerre-Voronoi单元孔隙率不变的假设，其计算公式为： 其中ui为子球i的主动移速，Scelli和Vcelli则表示子球i所属的Laguerre-Voronoi单元的表面积和体积，xc,yc,zc和xi,yi,zi分别是生物质颗粒和其第i个组成子球的中心坐标，χx，χy和χz分别是x，y，z方向上的收缩因子，用于实现各向异性的生物质颗粒收缩；S3，完成迭代并输出计算结果，所述的计算结果包括颗粒内部温度分布及演化、密度及孔隙率分布及演化、颗粒形态演化、颗粒质量演化。

全文数据：

权利要求：

百度查询： 浙江大学 一种基于球簇假设和Laguerre-Voronoi结构的生物质大颗粒热解仿真方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD