CN112560281B - 基于Fluent优化气流分离电工级氧化镁粉的方法 - Google Patents

基于Fluent优化气流分离电工级氧化镁粉的方法 Download PDF

Info

Publication number
CN112560281B
CN112560281B CN202011540321.2A CN202011540321A CN112560281B CN 112560281 B CN112560281 B CN 112560281B CN 202011540321 A CN202011540321 A CN 202011540321A CN 112560281 B CN112560281 B CN 112560281B
Authority
CN
China
Prior art keywords
model
classifier
grid
electrical grade
fluent
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
Application number
CN202011540321.2A
Other languages
English (en)
Other versions
CN112560281A (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.)
Shenyang Institute of Automation of CAS
Original Assignee
Shenyang Institute of Automation of CAS
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 Shenyang Institute of Automation of CAS filed Critical Shenyang Institute of Automation of CAS
Priority to CN202011540321.2A priority Critical patent/CN112560281B/zh
Publication of CN112560281A publication Critical patent/CN112560281A/zh
Application granted granted Critical
Publication of CN112560281B publication Critical patent/CN112560281B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205—Re-meshing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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/14—Force 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)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Combined Means For Separation Of Solids (AREA)

Abstract

本发明涉及基于Fluent优化气流分离电工级氧化镁粉的方法，包括：分区域建立气流分级机的几何模型；对气体流动过程采用三维稳态湍流模型进行模拟，建立数学模型，采用的控制方程包括：连续性方程、Navier‑Stokes方程和标准k‑ε模型；对气流分级机几何模型采用非结构性网格进行网格划分得到对应的气流分级机的网格文件；通过仿真计算进行两相流动模拟，得到一系列不同转速下气流分级机粗粉产品出口颗粒的粒度分布情况；获取的不同转速和风量下收集粗粉的粒度分布，并对比找出最佳运行参数，用于实际生产过程控制。

Description

基于Fluent优化气流分离电工级氧化镁粉的方法
技术领域
本发明涉及电工级氧化镁粉分离，尤其是涉及基于Fluent优化气流分离电工级氧化镁粉的方法。
背景技术
纵观镁化工的上下游产业，电工级氧化镁一直被誉为镁化工行业的明珠，也是镁化工行业的市场制高点。电工级氧化镁具有绝佳的电绝缘性、耐高温性和热传导性，是典型的科技含量高、附加值高的氧化镁产品，它是电加热(管状)元器件的生产中必不可少的绝缘填充材料，广泛用于核能、航空航天、家电等领域，具有极其重要的战略价值和市场空间。
电工级氧化镁粉是由菱镁矿石电熔成电熔镁砂，经破碎、振动筛筛分，再利用气流分级机的气流分离，将电工级氧化镁粉料进一步分离处理，得到符合粒度分布要求的电工级氧化镁粉。在制备过程中电工级氧化镁的主要形态为粉体颗粒，其目数分布是由用户所需密度及流速和制电热管大小而决定的。目数太大，其颗粒度小，单位重量上之表面积较大，容易吸潮；目数过小，其颗粒越大，压延后会损伤发热丝而影响电热管寿命。
目前，我国电工级氧化镁粉的加工设备陈旧、自动化程度低，导致产品质量波动大，产品一致性差。作为电工级氧化镁生产过程的关键工序，工人只凭经验调整气流分级机的参数，没有足够的理论依据支撑。因此，有必要提出一种优化气流分级机分离电工级氧化镁粉的方法，对实际生产提供理论支撑，使得电工级氧化镁粉的产品质量有所保证，从而提高成品率。
发明内容
为进一步提高电工级氧化镁粉的产品质量和一致性，本发明要解决的技术难点是找到一种优化气流分离电工级氧化镁粉的方法,对实际生产提供理论支撑，使得电工级氧化镁粉的产品质量有所保证，从而提高成品率。
本发明为实现上述目的所采用的技术方案是：
基于Fluent优化气流分离电工级氧化镁粉的方法，包括以下步骤：
步骤一、分区域建立气流分级机的几何模型，使得所建立的模型与生产电工级氧化镁粉的分级机设备尺寸一致；
步骤二、分级机内部为气－固两相流，气流作为连续相，对其流动过程采用三维稳态湍流模型进行模拟，并建立数学模型，采用的控制方程包括：连续性方程、Navier-Stokes方程和标准k-ε模型；
步骤三、对气流分级机几何模型采用非结构性网格进行网格划分，所建立模型需要对分级叶轮附近的网格划分进行细密处理，对于进料管、气流出口区域的网格划分采用较为疏松的处理方式，最后得到对应的气流分级机的网格文件；
步骤四、设置电工级氧化镁粉的物料性质、模拟计算的求解条件和边界条件，通过Fluent软件进行两相流动模拟，获得颗粒在气流分级机内部的运动情况，同时得到一系列不同转速下气流分级机粗粉产品出口颗粒的粒度分布情况；
步骤五、获取的不同转速和风量下收集粗粉的粒度分布，与实际生产过程中符合电工级氧化镁粉产品要求的粒度分布标准参数进行对比，找出对应的气流分级机的最佳运行参数，用于实际生产过程控制。
是将建模区域分两个部分：分级叶轮、分级机本体外壳分别建模。
所述网格划分是利用ICEM前处理软件实现的。
所述步骤三中的网格划分需要满足以下控制方程约束：连续性方程、Navier-Stokes方程和标准k-ε模型。
所述的步骤三中对气流分级机几何模型进行网格划分后，需要进行网格独立性验证，具体方法为：利用计算流体力学软件Fluent对模型中待划分区域进行网格加密处理，比对不同网格数量的计算结果的偏差，当观测参数发生变化小于5％，则认为该模型与网格无关，选择符合网格无关性要求的最小网格数量进行后续模拟计算。
所述通过仿真计算是通过将得到的气流分级机网格文件导入计算流体力学软件Fluent进行两相流动模拟实现的。
所述获取的不同转速和风量下收集粗粉的粒度分布为通过流体可视化后处理软件Tecplot实现的。
所述实际生产过程中符合电工级氧化镁粉产品要求的粒度分布标准参数为预先已知的。
本发明具有以下有益效果及优点：
1.利用Fluent优化气流分离电工级氧化镁粉的方法，能从理论上分析气流分级机转速与风量对电工级氧化镁粉的分离效果，得到较好的产品粒度分布；
2.该方法能够指导电工级氧化镁粉的实际生产，提高成品率，保证产品质量。
附图说明
图1为气流分级机示意图；
图2为基于Fluent优化气流分离电工级氧化镁粉的方法流程图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方法做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但本发明能够以很多不同于在此描述的其他方式来实施，本领域技术人员可以在不违背发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施的限制。
除非另有定义，本文所使用的所有技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。
如图1～2所示，基于Fluent优化气流分离电工级氧化镁粉的方法，包括以下步骤：
步骤一、建立气流分级机的几何模型，为计算方便，所述气流分级机几何模型进行简化处理，建模区域分两个部分，一部分为分级叶轮，一部分为分级机本体外壳。所建立的模型与生产电工级氧化镁粉的分级机设备尺寸一致。
步骤二、数学模型建立：分级机内部为气－固两相流，气流作为连续相，其流动过程采用三维稳态湍流模型进行模拟，控制方程主要包括：连续性方程、Navier-Stokes方程和标准k-ε模型。
连续性方程：
其中，ρ是密度，t是时间，u、v、w是速度矢量u在x、y和z方向的分量。
Navier-Stokes方程：
Su、Sv和Sw是单位体积流体受的外力，p是压力，常数μ是动力粘度。
标准k-ε模型：
其中t是时间，ρ是密度，ui是速度，k是湍动能，ε是耗散率，Gk是由于平均速度梯度引起的湍动能k的产生项，Gb是由于浮力引起的湍动能k的产生项，μt是湍流动力粘度，μ是层流动力粘度，Ym代表可压湍流中脉动扩张，C1ε、C2ε、σk和σε分别是经验常数。
模拟时忽略颗粒－颗粒之间的相互作用，忽略颗粒体积分数对作为连续相主体的气流流动产生的影响。模拟时考虑气体曳力、浮力、重力对颗粒运动产生的影响。
步骤三、利用ICEM前处理软件，对步骤一中的气流分级机几何模型采用非结构性网格进行网格划分，所建立模型需要对分级叶轮附近的网格划分进行细密处理，对于进料管、气流出口等区域的网格划分可以较为疏松，最后得到对应的气流分级机的网格文件。
步骤四、将步骤三得到的气流分级机网格文件导入计算流体力学软件Fluent进行两相流动模拟，同时设置电工级氧化镁粉的物料性质，模拟计算的求解条件和边界条件，当各个变量残差达到10-5即认为达到收敛，设计时间步长为10-2s，然后开始仿真计算，得到一系列不同转速下气流分级机粗粉(即电工级氧化镁合格产品)出口颗粒的粒度分布情况。所述物料性质为电工级氧化镁的密度、分离前电工级氧化镁粉各种粒径下颗粒所占百分比；所述求解条件为速度－压力耦合采用PISO算法，其他物理量采用二阶迎风格式。所述边界条件为分级机的进料管口和上端气流出口采用速度入口边界条件，下部出料口为标准压力出口，出口压力为常压。壁面为完全弹性壁面；所述粗粉产品出口颗粒的粒度分布情况为粗粉产品出口各种粒径的百分含量分布。
步骤五、将步骤四中得到的数据文件应用流体可视化后处理软件Tecplot获取不同转速和风量下收集粗粉(即电工级氧化镁合格产品)的粒度分布，同时与实际生产过程中符合电工级氧化镁粉产品要求的粒度分布(40-325目)进行对比，找出气流分级机的最佳运行参数。模拟结果表明，当分级机转速在1050r/min左右时，风流量为4030m3/h左右，能获得符合电工级氧化镁粉产品要求的粒度分布(40-325目)。转速较大，会使得比325目大的颗粒被细粉收集，成品粒度范围不合格；转速较小，则会使比40目大的颗粒进入粗粉中，同样会导致产品不合格。该方法可以指导电工级氧化镁粉的实际生产，提高电工级氧化镁粉的成品率，减少物料损失。运行参数为转速和风量。
进一步的，所述的步骤三中，对气流分级机几何模型进行网格划分后，进行网格独立性验证，具体方法为：利用计算流体力学软件Fluent对模型中分级叶轮附近的网格划分进行细密处理，比对不同网格数量的计算结果的偏差，当观测参数发生变化小于5％，则认为该模型与网格无关，选择符合网格无关性要求的最小网格数量进行后续模拟计算。所述计算结果为粗粉(即电工级氧化镁合格产品)的粒度分布结果。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰应视为本发明的保护范围。

Claims (5)

1.基于Fluent优化气流分离电工级氧化镁粉的方法，其特征在于，包括以下步骤：
步骤一、分区域建立气流分级机的几何模型，使得所建立的模型与生产电工级氧化镁粉的分级机设备尺寸一致；
步骤二、分级机内部为气－固两相流，气流作为连续相，对其流动过程采用三维稳态湍流模型进行模拟，并建立数学模型，采用的控制方程包括：连续性方程、Navier-Stokes方程和标准k-ε模型；
步骤三、对气流分级机几何模型采用非结构性网格进行网格划分，所建立模型需要对分级叶轮附近的网格划分进行细密处理，对于进料管、气流出口区域的网格划分采用较为疏松的处理方式，最后得到对应的气流分级机的网格文件；所述网格划分需要满足以下控制方程约束：连续性方程、Navier-Stokes方程和标准k-ε模型；所述对气流分级机几何模型进行网格划分后，需要进行网格独立性验证，具体方法为：利用计算流体力学软件Fluent对模型中待划分区域进行网格加密处理，比对不同网格数量的计算结果的偏差，当观测参数发生变化小于5％，则认为该模型与网格无关，选择符合网格无关性要求的最小网格数量进行后续模拟计算；
步骤四、设置电工级氧化镁粉的物料性质、模拟计算的求解条件和边界条件，将得到的气流分级机网格文件导入计算流体力学软件Fluent软件进行两相流动模拟，获得颗粒在气流分级机内部的运动情况，同时得到一系列不同转速下气流分级机粗粉产品出口颗粒的粒度分布情况；
步骤五、获取的不同转速和风量下收集粗粉的粒度分布，与实际生产过程中符合电工级氧化镁粉产品要求的粒度分布标准参数进行对比，找出对应的气流分级机的最佳运行参数，用于实际生产过程控制。
2.根据权利要求1所述的基于Fluent优化气流分离电工级氧化镁粉的方法，其特征在于，是将建模区域分两个部分：分级叶轮、分级机本体外壳分别建模。
3.根据权利要求1所述的基于Fluent优化气流分离电工级氧化镁粉的方法，其特征在于，所述网格划分是利用ICEM前处理软件实现的。
4.根据权利要求1所述的基于Fluent优化气流分离电工级氧化镁粉的方法，其特征在于，所述获取的不同转速和风量下收集粗粉的粒度分布为通过流体可视化后处理软件Tecplot实现的。
5.根据权利要求1所述的基于Fluent优化气流分离电工级氧化镁粉的方法，其特征在于，所述实际生产过程中符合电工级氧化镁粉产品要求的粒度分布标准参数为预先已知的。
CN202011540321.2A 2020-12-23 2020-12-23 基于Fluent优化气流分离电工级氧化镁粉的方法 Active CN112560281B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011540321.2A CN112560281B (zh) 2020-12-23 2020-12-23 基于Fluent优化气流分离电工级氧化镁粉的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011540321.2A CN112560281B (zh) 2020-12-23 2020-12-23 基于Fluent优化气流分离电工级氧化镁粉的方法

Publications (2)

Publication Number Publication Date
CN112560281A CN112560281A (zh) 2021-03-26
CN112560281B true CN112560281B (zh) 2023-08-01

Family

ID=75032319

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011540321.2A Active CN112560281B (zh) 2020-12-23 2020-12-23 基于Fluent优化气流分离电工级氧化镁粉的方法

Country Status (1)

Country Link
CN (1) CN112560281B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949227A (zh) * 2021-03-30 2021-06-11 中国华能集团清洁能源技术研究院有限公司 适用于复杂地形风场湍流强度确定方法、系统、设备及存储介质
CN113351472A (zh) * 2021-07-08 2021-09-07 大石桥市美尔镁制品有限公司 一种粉体除微粉装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1717529A (zh) * 2002-10-24 2006-01-04 国际壳牌研究有限公司 用于加热地下地岩层或者井孔的温控加热器
CN101107420A (zh) * 2004-04-23 2008-01-16 国际壳牌研究有限公司 用于对地下地层进行加热的温度受限加热器
CN101292115A (zh) * 2005-11-17 2008-10-22 美国莫博特克公司 反应物利用率得到改良的循环流化床锅炉
CN104384488A (zh) * 2014-11-04 2015-03-04 西宁特殊钢股份有限公司 一种防止钢水卷渣的钢包设计方法
CN105912745A (zh) * 2016-03-24 2016-08-31 安徽威达环保科技股份有限公司 袋式除尘器气流组织多参数优化方法
CN109985487A (zh) * 2019-04-03 2019-07-09 天津大学 一种基于Fluent的变压吸附碳捕集的优化方法
CN110795799A (zh) * 2019-10-21 2020-02-14 浙江理工大学 一种气力输送系统中阀门的磨损预测和优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100314788A1 (en) * 2006-08-18 2010-12-16 Cheng-Hung Hung Production of Ultrafine Particles in a Plasma System Having Controlled Pressure Zones

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1717529A (zh) * 2002-10-24 2006-01-04 国际壳牌研究有限公司 用于加热地下地岩层或者井孔的温控加热器
CN101107420A (zh) * 2004-04-23 2008-01-16 国际壳牌研究有限公司 用于对地下地层进行加热的温度受限加热器
CN101292115A (zh) * 2005-11-17 2008-10-22 美国莫博特克公司 反应物利用率得到改良的循环流化床锅炉
CN104384488A (zh) * 2014-11-04 2015-03-04 西宁特殊钢股份有限公司 一种防止钢水卷渣的钢包设计方法
CN105912745A (zh) * 2016-03-24 2016-08-31 安徽威达环保科技股份有限公司 袋式除尘器气流组织多参数优化方法
CN109985487A (zh) * 2019-04-03 2019-07-09 天津大学 一种基于Fluent的变压吸附碳捕集的优化方法
CN110795799A (zh) * 2019-10-21 2020-02-14 浙江理工大学 一种气力输送系统中阀门的磨损预测和优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于Fluent的岸边集装箱起重机外流场数值模拟;范磊;梅驷俊;宗晨;吕瑞军;;起重运输机械(第06期);全文 *
新型生物质颗粒回转燃烧器流动特性仿真分析;贾国海;李立君;高自成;李际平;陈喜龙;;江苏大学学报(自然科学版)(第06期);全文 *

Also Published As

Publication number Publication date
CN112560281A (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
CN112560281B (zh) 基于Fluent优化气流分离电工级氧化镁粉的方法
Vakamalla et al. Development of novel hydrocyclone designs for improved fines classification using multiphase CFD model
Yuan et al. Movement rules and screening characteristics of rice-threshed mixture separation through a cylinder sieve
Yu et al. Large eddy simulation of a planar jet flow with nanoparticle coagulation
CN110795778A (zh) 一种基于Fluent软件对旋流电解槽结构优化的方法
Sun et al. A new static cyclonic classifier: Flow characteristics, performance evaluation and industrial applications
CN110737935B (zh) 一种基于数字孪生的室内热环境建模方法
CN114912383A (zh) 一种基于fluent的旋转填料床流场分析方法
CN105045994B (zh) 一种优化带式焙烧机热工制度的方法
CN109285437A (zh) 一种可视化3d打印装备循环过滤系统模拟装置及方法
Huber et al. Evaluation of Surface Tension Models for SPH-Based Fluid Animations Using a Benchmark Test.
CN102567580A (zh) 基于计算流体力学的湿式摩擦片特性分析方法
Cao et al. A novel micro-spiral pneumatic selection system for the separation of fresh tea leaves
Chen et al. Rheological equations in asymptotic regimes of granular flow
Yang et al. Simulation of pore scale fluid flow of granular ore media in heap leaching based on realistic model
Peng et al. Experimental and numerical simulation study on the hydrodynamic characteristics of spherical and irregular-shaped particles in a 3D liquid-fluidized bed
Berdowski et al. 3D Lagrangian VPM: simulations of the near-wake of an actuator disc and horizontal axis wind turbine
CN113886964B (zh) 基于边界替代混合模型的飞行器气动分析方法
CN117113572A (zh) 一种基于数字孪生的高炉设计方法和系统
Chen et al. Numerical analysis of aerodynamic noise radiated from cross flow fan
Li et al. Study on Gas-solid Two-phase Flow Field of Ramon Mill Classifier
Liu et al. CFD simulation of the two-phase flow in a novel cooking tank for ethanol production
CN115935842A (zh) 一种模拟血管植入器械后血流量增幅的方法
Zhang et al. Numerical Simulation and Optimization Design of Particle Transport in 3D Printers Spiral Extrusion 3D Printers
Long Simulation analysis of bubble coalescence behavior characteristics in Newtonian fluids based on the phase field method

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