CN116976162A - 基于生物质锅炉有限元分析的脱硝设计仿真方法 - Google Patents

基于生物质锅炉有限元分析的脱硝设计仿真方法 Download PDF

Info

Publication number
CN116976162A
CN116976162A CN202310718598.7A CN202310718598A CN116976162A CN 116976162 A CN116976162 A CN 116976162A CN 202310718598 A CN202310718598 A CN 202310718598A CN 116976162 A CN116976162 A CN 116976162A
Authority
CN
China
Prior art keywords
wall
heat exchange
boiler
combustion
setting
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.)
Granted
Application number
CN202310718598.7A
Other languages
English (en)
Other versions
CN116976162B (zh
Inventor
熊建
高勇
任鲁军
李哲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan Optics Valley Lanyan New Energy Co ltd
Original Assignee
Wuhan Optics Valley Lanyan New Energy Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan Optics Valley Lanyan New Energy Co ltd filed Critical Wuhan Optics Valley Lanyan New Energy Co ltd
Priority to CN202310718598.7A priority Critical patent/CN116976162B/zh
Publication of CN116976162A publication Critical patent/CN116976162A/zh
Application granted granted Critical
Publication of CN116976162B publication Critical patent/CN116976162B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)

Abstract

一种基于生物质锅炉有限元分析的脱硝设计仿真方法,包括:步骤1:构建3D网格模型,步骤2:设定边界条件,步骤3:计算区域设置,步骤4:设置燃烧模拟求解器,步骤5:对模型进行流速分析、温度场分析、NOX生成分析,步骤6:根据分析结果设计脱硝设备,步骤7:根据标注的核心位置设计脱硝喷枪,步骤8:对步骤7的结果进行分析验证。本设计不仅能够对设计结果进行仿真论证方案的可行性,而且能够通过仿真结果设计脱硝装置,并模拟脱硝效果对脱硝方案进行优化设计。

Description

基于生物质锅炉有限元分析的脱硝设计仿真方法
技术领域
本发明涉及一种基于生物质锅炉有限元分析的脱硝设计仿真方法,具体适用于生物质锅炉的燃烧仿真分析以及脱硝设备设计。
背景技术
随着计算机及数值计算的发展,流体、燃烧、传热过程仿真及研究越来越重要,并借助计算机模拟技术在计算流体力学(Computational Fluid Dynamics,CFD)的突出表现来实现研究工作的可视化。
通过数值仿真技术,可根据湍流、组分的扩散及燃烧等的基本理论建立所需要的数学模型,得到压力、速度、组分反映场、以及污染物的浓度分布等。数值模拟方法能够有效的反映流动、燃烧、污染物的扩散过程,在工程中的应用越来越广,技术也愈来愈趋于成熟。
生物质燃烧锅炉内流体力学CFD模拟中,对锅炉炉排及炉膛整体建立三维几何模型和数学模型,用FLUENT软件对反应器内的流场及组分、燃烧过程进行模拟,分析反应器内的压力、流动、温度及组分的分布情况,分析速度场、压力场、燃烧及污染物(NOx)的分布。
针对降低污染物(NOx)问题,在锅炉顶部装枪,喷枪内口径为19mm,采用SNCR模型模拟脱氮过程。
发明内容
本发明的目的是克服现有技术中存在的锅炉设计研究凭借设计人经验难以进行可视化仿真分析的问题,提供了一种实现可视化仿真分析的基于生物质锅炉有限元分析的脱硝设计仿真方法。
为实现以上目的,本发明的技术解决方案是:
一种基于生物质锅炉有限元分析的脱硝设计仿真方法,所述仿真方法基于如下锅炉结构:
所述锅炉结构包括炉膛结构和燃烧送料结构,所述炉膛结构包括设置于锅炉壳体内的膜式壁炉膛、进水管道和蒸汽管道,所述膜式壁炉膛的底部连通进水管道,所述膜式壁炉膛的顶部连通蒸汽管道;
所述膜式壁管道包括从燃料进口端到燃料出口端方向依次设置的前拱换热壁、后拱换热壁、隔墙换热壁和后墙换热壁,所述膜式壁管道还包括对称设置的两组换热侧顶壁,所述前拱换热壁、隔墙换热壁、后墙换热壁和换热侧顶壁均为竖向管道的膜式壁结构,所述后拱换热壁的中、下部为膜式壁结构,所述后拱换热壁的顶部为类似栅栏的管道结构,所述前拱换热壁、后拱换热壁、隔墙换热壁、后墙换热壁和两组换热侧顶壁管路的底部与进水管道相连通,所述前拱换热壁、后拱换热壁、隔墙换热壁、后墙换热壁和两组换热侧顶壁管路的顶部与蒸汽管道相连通;
所述前拱换热壁的下方为燃烧送料结构的入口,后墙换热壁下方为燃烧送料结构的排灰通道,所述后拱换热壁的下端与后墙换热壁的下端一起将锅炉内部分隔为上方的炉膛空间和下方的燃烧空间;所述燃烧送料结构包括:倾斜向下的送料传送带以及多个位于送料传送带下方或者侧部的送风装置,所述送料传送带位于前拱换热壁和后拱换热壁的正下方;
所述前拱换热壁、隔墙换热壁和后墙换热壁的左右两侧分别与其对应侧的换热侧顶壁密封连接,所述后拱换热壁下端的膜式壁段的左右两侧分别与其对应侧的换热侧顶壁以及锅炉壳体的内壁密封连接,所述换热侧顶壁顶部与蒸汽管道密封连接,所述后墙换热壁上部设置有两个烟气出口,所述前拱换热壁与后拱换热壁之间形成有主燃室,所述后拱换热壁与隔墙换热壁之间形成有二燃室,所述隔墙换热壁与后墙换热壁之间形成有沉降室,所述主燃室通过二燃室与沉降室相连通形成炉膛的内部空间;
所述仿真方法包括如下步骤:
步骤1:构建3D网格模型,将锅炉炉膛的三维图纸导入到仿真软件中进行简化,然后生成网格化的有限元分析模型;
步骤1.1:建立锅炉几何模型,将整体结构设计完成的生物质锅炉的图纸导入到有限元模型构建的软件中,对导入的图纸进行分析,提取仿真分析的主要部分,即生物质燃料燃烧区域和炉膛区域;
步骤1.2:对构建的锅炉几何模型进行简化处理,将两组换热侧顶壁、前拱换热壁和部分后墙换热壁作为炉膛区域的外部轮廓边界,然后将送料传送带与锅炉的侧部壳体进行简化作为燃烧区域的外部轮廓边界,然后将后拱换热壁、隔墙换热壁和剩余后墙换热壁的膜式壁部分进行挖空简化处理,使其起到分隔作用,然后根据送料传送带底部的送风装置的布置对燃烧区域的底部轮廓进行划分作为一次风入口,最后将侧部送风口简化后还原到燃烧区域的侧壁上作为二次风入口,生成生物质锅炉的三维几何模型;
步骤1.3:对已建立的锅炉几何模型定义材料属性,根据实际设计情况对生物质锅炉几何模型的各壁面的材料属性进行定义;
步骤1.4:对已建立的锅炉几何模型进行网格划分,定义完成材料属性后,采用先生成四面体再转化为六面体网格的方式进行网格化处理,然后为了提高局部计算精度,将炉膛前壁面、进料口区域、几何体复杂区域和小尺寸风口区域进行加密处理;
步骤1.5:将网格锅炉几何模型生成锅炉整体有限元分析模型;
步骤2:设定边界条件:分析燃烧反应的物质成分、设置锅炉环境参数;入口边界条件、出口边界条件和壁面边界条件;
步骤2.1设置参与燃烧反应的生物质燃料参数和空气参数,首先对生物质燃料的的进行元素成分分析和工业分析,然后估算出生物质燃料中的挥发分份额,然后根据设计的燃料使用工况,设定锅炉的燃料输入情况;
然后对空气都参数进行设置,并设置锅炉使用的环境参数信息:温度、压力、相对湿度;
步骤2.2设定锅炉的入口边界条件,包括:空气温度参数、湍动能强度参数、湍动能粘度比参数、一次风量参数、二次风量参数;
步骤2.3设定锅炉的出口边界条件,包括:出口截面积参数、出口压力参数、出口温度参数、回流湍流强度参数、回流湍动能粘度比参数;
步骤2.4设定壁面边界条件:分别包括炉膛、水冷壁和炉排的温度参数、壁厚参数、热条件、内部发射率参数;
步骤3:计算区域设置:设定炉膛区域参数,然后设定炉膛区域的特性;
3.1设定炉膛区域参数:先设定多孔介质区域:水分、挥发分、固定碳的材料流阻;然后设定燃料层的粘性阻力、惯性阻力和流体孔隙度;
3.2设定炉膛区域的特性:将炉膛区域的流体域属性设定为烟气;
步骤4:设置燃烧模拟求解器,将生物质燃料模拟锅炉内的燃烧环境进行多次燃烧实验,在实验过程中收集相关的燃烧实验数据,然后对实验结果数据进行分析,根据分析结果在仿真软件中设定燃烧和NOX生成的化学反应数据、温度传递数据、压力数据和流体流动数据;模型选择上,涵盖生物质锅炉燃烧所涉及的模型包括能量方程、湍流方程、辐射模型、组分输运与反应模型质量守恒模型即连续性模型默认打开状态,其中组分输运与反应模型可以添加化学反应方程,材料设置中,打开组分输运与反应模型后,会出现Mixture多组分混合材料;需要添加烧炉燃烧过程中所涉及的所有组分,设定完成后将数据在软件中整合后形成求解器;
步骤5:对模型进行流速分析、温度场分析、NOX生成分析,利用步骤3中生成的求解器对有限元分析模型进行求解,查看其模型中气体流速、温度以及NOX生成情况;
步骤6:根据分析结果设计脱硝设备,由于NOX脱硝温度在800度到1250度,NOX最佳脱硝温度在900度到1000度,将有限元温度场分析的结果以及NOX生成浓度云图的结果进行对比分析,在合适脱硝温度的区域内找出NOX浓度开始极速增加点的核心位置进行标注;
步骤7:根据标注的核心位置设计脱硝喷枪,在有限元模型中在步骤6中标注的核心位置附近,寻找多个合适的点位,在这些点位两侧的炉壁上设置对称的脱硝喷枪点位,并设定脱硝喷枪的边界条件,以及求解器中的脱硝反应模型,即SNCR模型,通过模拟开启一组或多组喷枪,分别利用求解器进行求解,得出相对应的添加脱硝反应后的NOX浓度云图;
步骤8:对步骤7的结果进行分析验证,判定出最佳的脱硝喷枪设定组合方式,同时对上述组合方式进行判定,看判定输出NOX是否满足设备的污染物排放要求,如不满足则回到步骤7通过调整脱硝喷枪位置或增加脱硝喷枪数量的方式对脱硝设备进行进一步改进、优化设计,然后在步骤8中再次进行验证,直到仿真结果符合污染物排放要求,此时仿真设计完成。
与现有技术相比,本发明的有益效果为:
1、本发明一种基于生物质锅炉有限元分析的脱硝设计仿真方法中将锅炉的整体设计图导入到仿真软件中,然后根据锅炉设计的使用状态设计求解器,模拟锅炉内的流体场、压力场、温度场、组分分布和燃烧过程,根据分析结果能够找到设计缺陷、优化脱硝设计方案,进而得出可行性更高的脱硝设计方案。
2、本发明一种基于生物质锅炉有限元分析的脱硝设计仿真方法中网格在划分过程中为了匹配几何,先将其整体划分为四面体网格,然后将其大部分转化为六面体网格,在边界及交界面上产生了部分金字塔形网格和多面体网格。网格生成转化有利于在不变化网格尺寸的情况下大大减少网格的数量,即在不降低计算精度的情况下减少了求解器的计算量,这样能加快计算。
附图说明
图1是本发明分析目标的结构示意图。
图2是本发明步骤1中几何模型简化并网格化的流程示意图。
图3是本发明炉膛中心截面速度流线图。
图4是本发明炉膛中心截面温度云图。
图5是本发明炉膛中心截面NOx分布云图。
图6是本发明步骤6中的核心位置标注示意图。
图7是本发明添加SNCR后NOx分布仿真示意图。
图中:膜式壁炉膛1、主燃室11、二燃室12、沉降室13、前拱换热壁2、后拱换热壁3、隔墙换热壁4、后墙换热壁5、换热侧顶壁6、进水管道7、蒸汽管道8。
具体实施方式
以下结合附图说明和具体实施方式对本发明作进一步详细的说明。
参见图1至图7,一种基于生物质锅炉有限元分析的脱硝设计仿真方法,所述仿真方法基于如下锅炉结构:
所述锅炉结构包括炉膛结构和燃烧送料结构,所述炉膛结构包括设置于锅炉壳体内的膜式壁炉膛1、进水管道7和蒸汽管道8,所述膜式壁炉膛1的底部连通进水管道7,所述膜式壁炉膛1的顶部连通蒸汽管道8;
所述膜式壁管道1包括从燃料进口端到燃料出口端方向依次设置的前拱换热壁2、后拱换热壁3、隔墙换热壁4和后墙换热壁5,所述膜式壁管道1还包括对称设置的两组换热侧顶壁6,所述前拱换热壁2、隔墙换热壁4、后墙换热壁5和换热侧顶壁6均为竖向管道的膜式壁结构,所述后拱换热壁3的中、下部为膜式壁结构,所述后拱换热壁3的顶部为类似栅栏的管道结构,所述前拱换热壁2、后拱换热壁3、隔墙换热壁4、后墙换热壁5和两组换热侧顶壁6管路的底部与进水管道7相连通,所述前拱换热壁2、后拱换热壁3、隔墙换热壁4、后墙换热壁5和两组换热侧顶壁6管路的顶部与蒸汽管道8相连通;
所述前拱换热壁2的下方为燃烧送料结构的入口,后墙换热壁5下方为燃烧送料结构的排灰通道,所述后拱换热壁3的下端与后墙换热壁5的下端一起将锅炉内部分隔为上方的炉膛空间和下方的燃烧空间;所述燃烧送料结构包括:倾斜向下的送料传送带以及多个位于送料传送带下方或者侧部的送风装置,所述送料传送带位于前拱换热壁2和后拱换热壁3的正下方;
所述前拱换热壁2、隔墙换热壁4和后墙换热壁5的左右两侧分别与其对应侧的换热侧顶壁6密封连接,所述后拱换热壁3下端的膜式壁段的左右两侧分别与其对应侧的换热侧顶壁6以及锅炉壳体的内壁密封连接,所述换热侧顶壁6顶部与蒸汽管道8密封连接,所述后墙换热壁5上部设置有两个烟气出口9,所述前拱换热壁2与后拱换热壁3之间形成有主燃室11,所述后拱换热壁3与隔墙换热壁4之间形成有二燃室12,所述隔墙换热壁4与后墙换热壁5之间形成有沉降室13,所述主燃室11通过二燃室12与沉降室13相连通形成炉膛的内部空间;
所述仿真方法包括如下步骤:
步骤1:构建3D网格模型,将锅炉炉膛的三维图纸导入到仿真软件中进行简化,然后生成网格化的有限元分析模型;
步骤1.1:建立锅炉几何模型,将整体结构设计完成的生物质锅炉的图纸导入到有限元模型构建的软件中,对导入的图纸进行分析,提取仿真分析的主要部分,即生物质燃料燃烧区域和炉膛区域;
步骤1.2:对构建的锅炉几何模型进行简化处理,将两组换热侧顶壁6、前拱换热壁2和部分后墙换热壁5作为炉膛区域的外部轮廓边界,然后将送料传送带与锅炉的侧部壳体进行简化作为燃烧区域的外部轮廓边界,然后将后拱换热壁3、隔墙换热壁4和剩余后墙换热壁5的膜式壁部分进行挖空简化处理,使其起到分隔作用,然后根据送料传送带底部的送风装置的布置对燃烧区域的底部轮廓进行划分作为一次风入口,最后将侧部送风口简化后还原到燃烧区域的侧壁上作为二次风入口,生成生物质锅炉的三维几何模型;
步骤1.3:对已建立的锅炉几何模型定义材料属性,根据实际设计情况对生物质锅炉几何模型的各壁面的材料属性进行定义;
步骤1.4:对已建立的锅炉几何模型进行网格划分,定义完成材料属性后,采用先生成四面体再转化为六面体网格的方式进行网格化处理,然后为了提高局部计算精度,将炉膛前壁面、进料口区域、几何体复杂区域和小尺寸风口区域进行加密处理;
步骤1.5:将网格锅炉几何模型生成锅炉整体有限元分析模型;
步骤2:设定边界条件,分析燃烧反应的物质成分、设置锅炉环境参数;入口边界条件、出口边界条件和壁面边界条件;
步骤2.1设置参与燃烧反应的生物质燃料参数和空气参数,首先对生物质燃料的的进行元素成分分析和工业分析,然后估算出生物质燃料中的挥发分份额,然后根据设计的燃料使用工况,设定锅炉的燃料输入情况;
然后对空气都参数进行设置,并设置锅炉使用的环境参数信息:温度、压力、相对湿度;
步骤2.2设定锅炉的入口边界条件,包括:空气温度参数、湍动能强度参数、湍动能粘度比参数、一次风量参数、二次风量参数;
步骤2.3设定锅炉的出口边界条件,包括:出口截面积参数、出口压力参数、出口温度参数、回流湍流强度参数、回流湍动能粘度比参数;
步骤2.4设定壁面边界条件:分别包括炉膛、水冷壁和炉排的温度参数、壁厚参数、热条件、内部发射率参数;
步骤3:计算区域设置,设定炉膛区域参数,然后设定炉膛区域的特性
3.1设定炉膛区域参数:先设定多孔介质区域:水分、挥发分、固定碳的材料流阻;然后设定燃料层的粘性阻力、惯性阻力和流体孔隙度;
3.2设定炉膛区域的特性:将炉膛区域的流体域属性设定为烟气;
步骤4:设置燃烧模拟求解器,将生物质燃料模拟锅炉内的燃烧环境进行多次燃烧实验,在实验过程中收集相关的燃烧实验数据,然后对实验结果数据进行分析,根据分析结果在仿真软件中设定燃烧和NOX生成的化学反应数据、温度传递数据、压力数据和流体流动数据;模型选择上,涵盖生物质锅炉燃烧所涉及的模型包括能量方程、湍流方程、辐射模型、组分输运与反应模型质量守恒模型即连续性模型默认打开状态,其中组分输运与反应模型可以添加化学反应方程,材料设置中,打开组分输运与反应模型后,会出现Mixture多组分混合材料;需要添加烧炉燃烧过程中所涉及的所有组分,设定完成后将数据在软件中整合后形成求解器;
步骤5:对模型进行流速分析、温度场分析、NOX生成分析,利用步骤3中生成的求解器对有限元分析模型进行求解,查看其模型中气体流速、温度以及NOX生成情况;
步骤6:根据分析结果设计脱硝设备,由于NOX脱硝温度在800度到1250度,NOX最佳脱硝温度在900度到1000度,将有限元温度场分析的结果以及NOX生成浓度云图的结果进行对比分析,在合适脱硝温度的区域内找出NOX浓度开始极速增加点的核心位置进行标注;
步骤7:根据标注的核心位置设计脱硝喷枪,在有限元模型中在步骤6中标注的核心位置附近,寻找多个合适的点位,在这些点位两侧的炉壁上设置对称的脱硝喷枪点位,并设定脱硝喷枪的边界条件,以及求解器中的脱硝反应模型,即SNCR模型,通过模拟开启一组或多组喷枪,分别利用求解器进行求解,得出相对应的添加脱硝反应后的NOX浓度云图;
步骤8:对步骤7的结果进行分析验证,判定出最佳的脱硝喷枪设定组合方式,同时对上述组合方式进行判定,看判定输出NOX是否满足设备的污染物排放要求,如不满足则回到步骤7通过调整脱硝喷枪位置或增加脱硝喷枪数量的方式对脱硝设备进行进一步改进、优化设计,然后在步骤8中再次进行验证,直到仿真结果符合污染物排放要求,此时仿真设计完成。
本发明的原理说明如下:
本设计通过仿真认证的方法对设计成品的工作状态进行仿真分析,然后根据仿真结果优化设计,进而得出理论上的设计依据,有效提高了设计效率,降低了实验成本。
实施例1:
一种基于生物质锅炉有限元分析的脱硝设计仿真方法,其特征在于:
所述仿真方法基于如下锅炉结构:
所述锅炉结构包括炉膛结构和燃烧送料结构,所述炉膛结构包括设置于锅炉壳体内的膜式壁炉膛1、进水管道7和蒸汽管道8,所述膜式壁炉膛1的底部连通进水管道7,所述膜式壁炉膛1的顶部连通蒸汽管道8;
所述膜式壁管道1包括从燃料进口端到燃料出口端方向依次设置的前拱换热壁2、后拱换热壁3、隔墙换热壁4和后墙换热壁5,所述膜式壁管道1还包括对称设置的两组换热侧顶壁6,所述前拱换热壁2、隔墙换热壁4、后墙换热壁5和换热侧顶壁6均为竖向管道的膜式壁结构,所述后拱换热壁3的中、下部为膜式壁结构,所述后拱换热壁3的顶部为类似栅栏的管道结构,所述前拱换热壁2、后拱换热壁3、隔墙换热壁4、后墙换热壁5和两组换热侧顶壁6管路的底部与进水管道7相连通,所述前拱换热壁2、后拱换热壁3、隔墙换热壁4、后墙换热壁5和两组换热侧顶壁6管路的顶部与蒸汽管道8相连通;
所述前拱换热壁2的下方为燃烧送料结构的入口,后墙换热壁5下方为燃烧送料结构的排灰通道,所述后拱换热壁3的下端与后墙换热壁5的下端一起将锅炉内部分隔为上方的炉膛空间和下方的燃烧空间;所述燃烧送料结构包括:倾斜向下的送料传送带以及多个位于送料传送带下方或者侧部的送风装置,所述送料传送带位于前拱换热壁2和后拱换热壁3的正下方;
所述前拱换热壁2、隔墙换热壁4和后墙换热壁5的左右两侧分别与其对应侧的换热侧顶壁6密封连接,所述后拱换热壁3下端的膜式壁段的左右两侧分别与其对应侧的换热侧顶壁6以及锅炉壳体的内壁密封连接,所述换热侧顶壁6顶部与蒸汽管道8密封连接,所述后墙换热壁5上部设置有两个烟气出口9,所述前拱换热壁2与后拱换热壁3之间形成有主燃室11,所述后拱换热壁3与隔墙换热壁4之间形成有二燃室12,所述隔墙换热壁4与后墙换热壁5之间形成有沉降室13,所述主燃室11通过二燃室12与沉降室13相连通形成炉膛的内部空间;
所述仿真方法包括如下步骤:
步骤1:构建3D网格模型,将锅炉炉膛的三维图纸导入到仿真软件中进行简化,然后生成网格化的有限元分析模型;
步骤1.1:建立锅炉几何模型,将整体结构设计完成的生物质锅炉的图纸导入到有限元模型构建的软件中,对导入的图纸进行分析,提取仿真分析的主要部分,即生物质燃料燃烧区域和炉膛区域;
步骤1.2:对构建的锅炉几何模型进行简化处理,将两组换热侧顶壁6、前拱换热壁2和部分后墙换热壁5作为炉膛区域的外部轮廓边界,然后将送料传送带与锅炉的侧部壳体进行简化作为燃烧区域的外部轮廓边界,然后将后拱换热壁3、隔墙换热壁4和剩余后墙换热壁5的膜式壁部分进行挖空简化处理,使其起到分隔作用,然后根据送料传送带底部的送风装置的布置对燃烧区域的底部轮廓进行划分作为一次风入口,最后将侧部送风口简化后还原到燃烧区域的侧壁上作为二次风入口,生成生物质锅炉的三维几何模型;
步骤1.3:对已建立的锅炉几何模型定义材料属性,根据实际设计情况对生物质锅炉几何模型的各壁面的材料属性进行定义;
步骤1.4:对已建立的锅炉几何模型进行网格划分,定义完成材料属性后,采用先生成四面体再转化为六面体网格的方式进行网格化处理,然后为了提高局部计算精度,将炉膛前壁面、进料口区域、几何体复杂区域和小尺寸风口区域进行加密处理;
经过ICEM软件对网格进行光顺,所有网格质量都在0.21以上,能够满足燃烧所需的高精度计算。
网格报告信息显示整体网格数量为85.5万。网格在划分过程中为了匹配几何,先将其整体划分为四面体网格,然后将其大部分转化为六面体网格,在边界及交界面上产生了部分金字塔形网格和多面体网格。网格生成转化有利于在不变化网格尺寸的情况下大大减少网格的数量(在原来的基础上至少减少了2/3),即在不降低计算精度的情况下减少了求解器的计算量,这样能加快计算。
步骤1.5:将网格锅炉几何模型生成锅炉整体有限元分析模型;
步骤2:设定边界条件:分析燃烧反应的物质成分、设置锅炉环境参数;入口边界条件、出口边界条件和壁面边界条件;
步骤2.1设置参与燃烧反应的生物质燃料参数和空气参数,首先对生物质燃料的的进行元素成分分析和工业分析,然后估算出生物质燃料中的挥发分份额,然后根据设计的燃料使用工况,设定锅炉的燃料输入情况;
表1燃料特性
表2工业分析
表3挥发分份额(估算值)
表4燃料输入情况
然后。对空气都参数进行设置,并设置锅炉使用的环境参数信息:温度、压力、相对湿度;
表5空气参数
表6环境参数
步骤2.2设定锅炉的入口边界条件,包括:空气温度参数、湍动能强度参数、湍动能粘度比参数、一次风量参数、二次风量参数;
表7入口边界条件
步骤2.3设定锅炉的出口边界条件,包括:出口截面积参数、出口压力参数、出口温度参数、回流湍流强度参数、回流湍动能粘度比参数;
表4-8出口边界条件
步骤2.4设定壁面边界条件:分别包括炉膛、水冷壁和炉排的温度参数、壁厚参数、热条件、内部发射率参数;
表9壁面边界条件
壁面 温度(℃) 壁厚(m) 热条件 内部发射率
炉膛 264(537.15K) 0 绝热 0.7
水冷壁 264(537.15K) 0.06 给定温度 0.7
炉排 264(537.15K) 0 绝热 0.7
步骤3:计算区域设置:设定炉膛区域参数,然后设定炉膛区域的特性;
3.1设定炉膛区域参数:先设定多孔介质区域:水分、挥发分、固定碳的材料流阻;然后设定燃料层的粘性阻力、惯性阻力和流体孔隙度;
表4-10炉排区域
3.2设定炉膛区域的特性:将炉膛区域的流体域属性设定为烟气;
步骤4:设置燃烧模拟求解器,将生物质燃料模拟锅炉内的燃烧环境进行多次燃烧实验,在实验过程中收集相关的燃烧实验数据,然后对实验结果数据进行分析,根据分析结果在仿真软件中设定燃烧和NOX生成的化学反应数据、温度传递数据、压力数据和流体流动数据;模型选择上,涵盖生物质锅炉燃烧所涉及的模型包括能量方程、湍流方程、辐射模型、组分输运与反应模型质量守恒模型即连续性模型默认打开状态,其中组分输运与反应模型可以添加化学反应方程,材料设置中,打开组分输运与反应模型后,会出现Mixture多组分混合材料;需要添加烧炉燃烧过程中所涉及的所有组分,设定完成后将数据在软件中整合后形成求解器;
步骤5:对模型进行流速分析、温度场分析、NOX生成分析,利用步骤3中生成的求解器对有限元分析模型进行求解,查看其模型中气体流速、温度以及NOX生成情况;
步骤6:根据分析结果设计脱硝设备,由于NOX脱硝温度在800度到1250度,NOX最佳脱硝温度在900度到1000度,将有限元温度场分析的结果以及NOX生成浓度云图的结果进行对比分析,在合适脱硝温度的区域内找出NOX浓度开始极速增加点的核心位置进行标注;
步骤7:根据标注的核心位置设计脱硝喷枪,在有限元模型中在步骤6中标注的核心位置附近,寻找多个合适的点位,在这些点位两侧的炉壁上设置对称的脱硝喷枪点位,并设定脱硝喷枪的边界条件,以及求解器中的脱硝反应模型,即SNCR模型,通过模拟开启一组或多组喷枪,分别利用求解器进行求解,得出相对应的添加脱硝反应后的NOX浓度云图;
步骤8:对步骤7的结果进行分析验证,判定出最佳的脱硝喷枪设定组合方式,同时对上述组合方式进行判定,看判定输出NOX是否满足设备的污染物排放要求,如不满足则回到步骤7通过调整脱硝喷枪位置或增加脱硝喷枪数量的方式对脱硝设备进行进一步改进、优化设计,然后在步骤8中再次进行验证,直到仿真结果符合污染物排放要求,此时仿真设计完成。

Claims (4)

1.一种基于生物质锅炉有限元分析的脱硝设计仿真方法,其特征在于:
所述仿真方法基于如下锅炉结构:
所述锅炉结构包括炉膛结构和燃烧送料结构,所述炉膛结构包括设置于锅炉壳体内的膜式壁炉膛(1)、进水管道(7)和蒸汽管道(8),所述膜式壁炉膛(1)的底部连通进水管道(7),所述膜式壁炉膛(1)的顶部连通蒸汽管道(8);
所述膜式壁管道(1)包括从燃料进口端到燃料出口端方向依次设置的前拱换热壁(2)、后拱换热壁(3)、隔墙换热壁(4)和后墙换热壁(5),所述膜式壁管道(1)还包括对称设置的两组换热侧顶壁(6),所述前拱换热壁(2)、隔墙换热壁(4)、后墙换热壁(5)和换热侧顶壁(6)均为竖向管道的膜式壁结构,所述后拱换热壁(3)的中、下部为膜式壁结构,所述后拱换热壁(3)的顶部为类似栅栏的管道结构,所述前拱换热壁(2)、后拱换热壁(3)、隔墙换热壁(4)、后墙换热壁(5)和两组换热侧顶壁(6)管路的底部与进水管道(7)相连通,所述前拱换热壁(2)、后拱换热壁(3)、隔墙换热壁(4)、后墙换热壁(5)和两组换热侧顶壁(6)管路的顶部与蒸汽管道(8)相连通;
所述前拱换热壁(2)的下方为燃烧送料结构的入口,后墙换热壁(5)下方为燃烧送料结构的排灰通道,所述后拱换热壁(3)的下端与后墙换热壁(5)的下端一起将锅炉内部分隔为上方的炉膛空间和下方的燃烧空间;所述燃烧送料结构包括:倾斜向下的送料传送带以及多个位于送料传送带下方或者侧部的送风装置,所述送料传送带位于前拱换热壁(2)和后拱换热壁(3)的正下方;
所述前拱换热壁(2)、隔墙换热壁(4)和后墙换热壁(5)的左右两侧分别与其对应侧的换热侧顶壁(6)密封连接,所述后拱换热壁(3)下端的膜式壁段的左右两侧分别与其对应侧的换热侧顶壁(6)以及锅炉壳体的内壁密封连接,所述换热侧顶壁(6)顶部与蒸汽管道(8)密封连接,所述后墙换热壁(5)上部设置有两个烟气出口(9),所述前拱换热壁(2)与后拱换热壁(3)之间形成有主燃室(11),所述后拱换热壁(3)与隔墙换热壁(4)之间形成有二燃室(12),所述隔墙换热壁(4)与后墙换热壁(5)之间形成有沉降室(13),所述主燃室(11)通过二燃室(12)与沉降室(13)相连通形成炉膛的内部空间;
所述仿真方法包括如下步骤:
步骤1:构建3D网格模型,将锅炉炉膛的三维图纸导入到仿真软件中进行简化,然后生成网格化的有限元分析模型;
步骤2:设定边界条件:分析燃烧反应的物质成分、设置锅炉环境参数;入口边界条件、出口边界条件和壁面边界条件;
步骤3:计算区域设置:设定炉膛区域参数,然后设定炉膛区域的特性;
步骤4:设置燃烧模拟求解器,将生物质燃料模拟锅炉内的燃烧环境进行多次燃烧实验,在实验过程中收集相关的燃烧实验数据,然后对实验结果数据进行分析,根据分析结果在仿真软件中设定燃烧和NOX生成的化学反应数据、温度传递数据、压力数据和流体流动数据;模型选择上,涵盖生物质锅炉燃烧所涉及的模型包括能量方程、湍流方程、辐射模型、组分输运与反应模型质量守恒模型即连续性模型默认打开状态,其中组分输运与反应模型可以添加化学反应方程,材料设置中,打开组分输运与反应模型后,会出现Mixture多组分混合材料;需要添加烧炉燃烧过程中所涉及的所有组分,设定完成后将数据在软件中整合后形成求解器;
步骤5:对模型进行流速分析、温度场分析、NOX生成分析,利用步骤3中生成的求解器对有限元分析模型进行求解,查看其模型中气体流速、温度以及NOX生成情况;
步骤6:根据分析结果设计脱硝设备,由于NOX脱硝温度在800度到1250度,NOX最佳脱硝温度在900度到1000度,将有限元温度场分析的结果以及NOX生成浓度云图的结果进行对比分析,在合适脱硝温度的区域内找出NOX浓度开始极速增加点的核心位置进行标注;
步骤7:根据标注的核心位置设计脱硝喷枪,在有限元模型中在步骤6中标注的核心位置附近,寻找多个合适的点位,在这些点位两侧的炉壁上设置对称的脱硝喷枪点位,并设定脱硝喷枪的边界条件,以及求解器中的脱硝反应模型,即SNCR模型,通过模拟开启一组或多组喷枪,分别利用求解器进行求解,得出相对应的添加脱硝反应后的NOX浓度云图;
步骤8:对步骤7的结果进行分析,判定出最佳的脱硝喷枪设定组合方式,同时对上述组合方式进行判定,看判定输出NOX是否满足设备的污染物排放要求,如不满足则回到步骤7通过调整脱硝喷枪位置或增加脱硝喷枪数量的方式对脱硝设备进行进一步改进、优化设计,然后在步骤8中再次进行验证,直到仿真结果符合污染物排放要求,此时仿真设计完成。
2.根据权利要求1所述的一种基于生物质锅炉有限元分析的脱硝设计仿真方法,其特征在于:所述步骤1包括:
步骤 1.1 :建立锅炉几何模型,将整体结构设计完成的生物质锅炉的图纸导入到有限元模型构建的软件中,对导入的图纸进行分析,提取仿真分析的主要部分,即生物质燃料燃烧区域和炉膛区域;
步骤 1.2 :对构建的锅炉几何模型进行简化处理,将两组换热侧顶壁(6)、前拱换热壁(2)和部分后墙换热壁(5)作为炉膛区域的外部轮廓边界,然后将送料传送带与锅炉的侧部壳体进行简化作为燃烧区域的外部轮廓边界,然后将后拱换热壁(3)、隔墙换热壁(4)和剩余后墙换热壁(5)的膜式壁部分进行挖空简化处理,使其起到分隔作用,然后根据送料传送带底部的送风装置的布置对燃烧区域的底部轮廓进行划分作为一次风入口,最后将侧部送风口简化后还原到燃烧区域的侧壁上作为二次风入口,生成生物质锅炉的三维几何模型;
步骤 1.3:对已建立的锅炉几何模型定义材料属性,根据实际设计情况对生物质锅炉几何模型的各壁面的材料属性进行定义;
步骤 1.4 :对已建立的锅炉几何模型进行网格划分,定义完成材料属性后,采用先生成四面体再转化为六面体网格的方式进行网格化处理,然后为了提高局部计算精度,将炉膛前壁面、进料口区域、几何体复杂区域和小尺寸风口区域进行加密处理;
步骤 1.5 :将网格锅炉几何模型生成锅炉整体有限元分析模型。
3.根据权利要求2所述的一种基于生物质锅炉有限元分析的脱硝设计仿真方法,其特征在于:所述步骤2包括:
步骤2.1设置参与燃烧反应的生物质燃料参数和空气参数,首先对生物质燃料的的进行元素成分分析和工业分析,然后估算出生物质燃料中的挥发分份额,然后根据设计的燃料使用工况,设定锅炉的燃料输入情况;
然后对空气都参数进行设置,并设置锅炉使用的环境参数信息:温度、压力、相对湿度;
步骤2.2设定锅炉的入口边界条件,包括:空气温度参数、湍动能强度参数、湍动能粘度比参数、一次风量参数、二次风量参数;
步骤2.3设定锅炉的出口边界条件,包括:出口截面积参数、出口压力参数、出口温度参数、回流湍流强度参数、回流湍动能粘度比参数;
步骤2.4设定壁面边界条件:分别包括炉膛、水冷壁和炉排的温度参数、壁厚参数、热条件、内部发射率参数。
4.根据权利要求3所述的一种基于生物质锅炉有限元分析的脱硝设计仿真方法,其特征在于:所述步骤3包括:
步骤3.1设定炉膛区域参数:先设定多孔介质区域:水分、挥发分、固定碳的材料流阻;然后设定燃料层的粘性阻力、惯性阻力和流体孔隙度;
步骤3.2设定炉膛区域的特性:将炉膛区域的流体域属性设定为烟气。
CN202310718598.7A 2023-06-16 2023-06-16 基于生物质锅炉有限元分析的脱硝设计仿真方法 Active CN116976162B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310718598.7A CN116976162B (zh) 2023-06-16 2023-06-16 基于生物质锅炉有限元分析的脱硝设计仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310718598.7A CN116976162B (zh) 2023-06-16 2023-06-16 基于生物质锅炉有限元分析的脱硝设计仿真方法

Publications (2)

Publication Number Publication Date
CN116976162A true CN116976162A (zh) 2023-10-31
CN116976162B CN116976162B (zh) 2024-01-30

Family

ID=88473928

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310718598.7A Active CN116976162B (zh) 2023-06-16 2023-06-16 基于生物质锅炉有限元分析的脱硝设计仿真方法

Country Status (1)

Country Link
CN (1) CN116976162B (zh)

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009228918A (ja) * 2008-03-19 2009-10-08 Hitachi Ltd ボイラの制御装置、及びボイラの制御方法
CN102252779A (zh) * 2011-04-28 2011-11-23 上海交通大学 基于烟气能量平衡的炉膛出口烟温优化测量方法
CN104147916A (zh) * 2014-07-31 2014-11-19 浙江天蓝环保技术股份有限公司 基于fluent设置循环流化床锅炉的SNCR喷枪方法
CN106649917A (zh) * 2016-09-12 2017-05-10 广东电网有限责任公司电力科学研究院 一种电站锅炉煤粉在变氧量下的燃烧特性生成数值的模拟方法及装置
CN107748816A (zh) * 2017-10-19 2018-03-02 中冶华天工程技术有限公司 一种脱硝反应器及其支架一体化设计方法
CN108870384A (zh) * 2018-04-28 2018-11-23 东南大学 低氮燃烧循环流化床锅炉的燃烧与sncr脱硝协同优化方法
CN110263467A (zh) * 2019-06-26 2019-09-20 东北电力大学 一种基于场协同原理的四角切圆锅炉数值分析方法
CN111695249A (zh) * 2020-05-29 2020-09-22 广东省特种设备检测研究院顺德检测院 一种燃气锅炉热效率的预测方法
CN111841275A (zh) * 2019-10-08 2020-10-30 玖龙纸业(东莞)有限公司 一种适用于锅炉超洁净排放方法
CN113941251A (zh) * 2021-11-08 2022-01-18 西安热工研究院有限公司 一种烟气脱硝流场模拟校核方法和系统
CN114239430A (zh) * 2021-12-03 2022-03-25 东南大学 一种基于数值模拟的炉膛出口NOx预测方法及系统
CN114611436A (zh) * 2022-03-31 2022-06-10 合肥希普蓝达环境科技发展有限公司 一种烟气脱硝设备cfd流场优化方法
CN114811595A (zh) * 2022-04-26 2022-07-29 中国电建集团河南工程有限公司 一种垃圾焚烧电站锅炉排放控制调试方法
CN115659760A (zh) * 2022-11-16 2023-01-31 河北省机电一体化中试基地有限公司 一种基于Fluent的热熔胶胶枪分析和结构优化方法
CN115935837A (zh) * 2022-11-10 2023-04-07 天津大唐国际盘山发电有限责任公司 一种用于模拟锅炉尾部烟气流动过程的数值计算方法
WO2023065719A1 (zh) * 2021-10-21 2023-04-27 西安热工研究院有限公司 一种锅炉炉膛数值模拟与水动力协同集成计算方法

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009228918A (ja) * 2008-03-19 2009-10-08 Hitachi Ltd ボイラの制御装置、及びボイラの制御方法
CN102252779A (zh) * 2011-04-28 2011-11-23 上海交通大学 基于烟气能量平衡的炉膛出口烟温优化测量方法
CN104147916A (zh) * 2014-07-31 2014-11-19 浙江天蓝环保技术股份有限公司 基于fluent设置循环流化床锅炉的SNCR喷枪方法
CN106649917A (zh) * 2016-09-12 2017-05-10 广东电网有限责任公司电力科学研究院 一种电站锅炉煤粉在变氧量下的燃烧特性生成数值的模拟方法及装置
CN107748816A (zh) * 2017-10-19 2018-03-02 中冶华天工程技术有限公司 一种脱硝反应器及其支架一体化设计方法
CN108870384A (zh) * 2018-04-28 2018-11-23 东南大学 低氮燃烧循环流化床锅炉的燃烧与sncr脱硝协同优化方法
CN110263467A (zh) * 2019-06-26 2019-09-20 东北电力大学 一种基于场协同原理的四角切圆锅炉数值分析方法
CN111841275A (zh) * 2019-10-08 2020-10-30 玖龙纸业(东莞)有限公司 一种适用于锅炉超洁净排放方法
CN111695249A (zh) * 2020-05-29 2020-09-22 广东省特种设备检测研究院顺德检测院 一种燃气锅炉热效率的预测方法
WO2023065719A1 (zh) * 2021-10-21 2023-04-27 西安热工研究院有限公司 一种锅炉炉膛数值模拟与水动力协同集成计算方法
CN113941251A (zh) * 2021-11-08 2022-01-18 西安热工研究院有限公司 一种烟气脱硝流场模拟校核方法和系统
CN114239430A (zh) * 2021-12-03 2022-03-25 东南大学 一种基于数值模拟的炉膛出口NOx预测方法及系统
CN114611436A (zh) * 2022-03-31 2022-06-10 合肥希普蓝达环境科技发展有限公司 一种烟气脱硝设备cfd流场优化方法
CN114811595A (zh) * 2022-04-26 2022-07-29 中国电建集团河南工程有限公司 一种垃圾焚烧电站锅炉排放控制调试方法
CN115935837A (zh) * 2022-11-10 2023-04-07 天津大唐国际盘山发电有限责任公司 一种用于模拟锅炉尾部烟气流动过程的数值计算方法
CN115659760A (zh) * 2022-11-16 2023-01-31 河北省机电一体化中试基地有限公司 一种基于Fluent的热熔胶胶枪分析和结构优化方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
江建军等: "醋纤生产中锅炉废气综合治理技术研究进展", 化工时刊, vol. 35, no. 9, pages 12 - 16 *
王宁;郭科宏;: "蓄热式燃气锅炉SNCR脱硝试验及数值模拟研究", 工业加热, no. 03 *
陈黎等: "烟气再循环脱硝技术的创新应用", 工业加热, vol. 51, no. 1, pages 51 - 54 *
高宁;张向宇;张波;陆续;向小凤;徐宏杰;: "SNCR脱硝系统长喷枪应力有限元模拟分析", 热力发电, no. 04 *

Also Published As

Publication number Publication date
CN116976162B (zh) 2024-01-30

Similar Documents

Publication Publication Date Title
CN103364212A (zh) 煤粉锅炉改烧烟煤变负荷燃烧特性测量方法与系统
CN115034138A (zh) 以计算流体力学与深度学习结合的锅炉温度场预测方法
Klačková et al. Simulation of wood biomass combustion in hot water boiler
CN106548032A (zh) 基于锅炉烟道网格检测烟气co和o2浓度的分布场重建方法
CN104461688B (zh) 一种详细燃烧化学反应机理骨架简化的方法
Bianco et al. CFD analysis and design optimization of an air manifold for a biomass boiler
CN116976162B (zh) 基于生物质锅炉有限元分析的脱硝设计仿真方法
Mazaheri et al. Numerical optimization of laboratory combustor geometry for NO suppression
CN103631997A (zh) 锅炉燃烧器的建模方法
Durcansky Comparison of mathematical models for heat exchangers of unconventional CHP units
CN104729855A (zh) 一种实验气体污染对超燃冲压发动机性能影响的评估方法
CN103324991B (zh) 煤粉锅炉NOx排放预测的方法与系统
Janiga et al. Reducing the CO emissions in a laminar burner using different numerical optimization methods
Roslyakov et al. Experience gained with CFD-modeling of liquid and gaseous fuel combustion processes in power installations
CN103323573A (zh) 煤粉锅炉改烧烟煤变磨煤机组合燃烧特性测量方法与系统
Malalasekera et al. Calculation of radiative heat transfer in combustion systems
Kovalnogov et al. Software and information complex for the coupled numerical solution of the equations of heat-and-moisture transfer and the studying of heat-and-humidity kinetics of biofuel cells
Husin et al. Modification of perforated plate in fluidized-bed combustor chamber through computational fluid dynamics simulation
Besseron et al. Toward High-Performance Multi-Physics Coupled Simulations for the Industry with XDEM
Hernik Numerical calculations of the WR-40 boiler with a furnace jet boiler system
CN103335327A (zh) 煤粉锅炉改烧烟煤变燃尽风量燃烧特性测量方法与系统
Srichat et al. Comparison of the heat for rectangular and circle salt boiling stove in Ban Dung area, Udon Thani Province.
Akinshilo et al. Thermal analysis and cost evaluation of boilers considering different insulation materials and fuel oil
Szubel et al. Application of various' response surface'based algorithms in optimization of air manifolds for batch boilers
CN115374563A (zh) 一种利用COMSOL Multiphysics软件优化干熄焦循环风量的方法

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