CN106021820A - 一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法 - Google Patents
一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法 Download PDFInfo
- Publication number
- CN106021820A CN106021820A CN201610490783.5A CN201610490783A CN106021820A CN 106021820 A CN106021820 A CN 106021820A CN 201610490783 A CN201610490783 A CN 201610490783A CN 106021820 A CN106021820 A CN 106021820A
- Authority
- CN
- China
- Prior art keywords
- granule
- fly ash
- collision
- monte carlo
- furnace
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法,属于火力发电厂炉内积灰预测领域。该方法包括如下步骤:1)采用计算流体动力学方法计算积灰采样点附近流场;2)将计算得到的流场条件导入颗粒运动方程;3)采用直接模拟蒙特卡洛方法进行飞灰颗粒间碰撞的计算;4)采样枪提供壁面热流,然后求解飞灰粘附率方程计算飞灰颗粒在采样面上的沉积过程。计算结果表明,本发明构建的炉内积灰模型能准确预测炉内飞灰颗粒的沉积速率,并给出详细的积灰内层飞灰颗粒粒径分布,成功实现炉内积灰过程的动态预测。
Description
技术领域
本发明涉及一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法,属于火力发电厂炉内积灰预测领域。
背景技术
炉膛积灰是电站锅炉运行中常见也是较难解决的问题,它极大地影响了锅炉运行的安全性和经济性。目前电站锅炉设计中针对积灰问题主要采用经验估计的方法,如何准确的预测炉内积灰过程,是实现燃煤火电机组节能减排的关键。
影响积灰的因素有很多。发生炉内积灰的两个关键条件是:第一,飞灰颗粒能够碰撞到壁面上;第二,这些飞灰颗粒具有适当的粘度能够附着在壁面上。显然经验参数模型很难准确的描述积灰沉积过程,现有的离散动力学方法可以求解单颗粒的运动方程,从而跟踪单颗粒的运动轨迹。但当模拟较小颗粒直径,且颗粒数较大的高浓度颗粒物运动时,确定性方法就会因为过高的计算负荷和计算机储存需求而无能为力。
发明内容
本发明的目的是提出一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法,准确的描述炉内积灰沉积过程。直接模拟蒙特卡洛方法使用较少的模拟颗粒代表大量的真实颗粒,因而有效地节省了计算时间。为了准确的描述积灰沉积过程,获得积灰沉积速率以及积灰内层的颗粒粒径分布,现提出一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法。
本发明的技术方案如下:
一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法,其特征在于所述方法包括如下步骤:
1)作为直接模拟蒙特卡洛计算的输入条件,需要先对采样点附近流场进行计算,采用计算流体动力学方法,其主要控制方程是稳态不可压缩流动纳维叶-斯托克斯方程:
式中,ρ为当地流体的密度,U为当地流体的速度,p为当地流体的压力,μ为动力学粘度,S为应变率张量;
2)飞灰颗粒在外力作用下运动满足牛顿第二运动定律,仅考虑飞灰颗粒的重力、气体施加的曳力以及热泳力,颗粒运动方程可表示为:
Fd=-3πdμ(v-u)f (3)
式中,Fd为流体施加给颗粒的曳力,Fth为热泳力,m为颗粒质量,g为重力加速度,d为颗粒直径,μ为空气动力学粘度系数,v为颗粒运动速度,u为气体速度,f为曳力系数,κtr为气体导热系数,r为颗粒半径,为气体分子的平均运动速度,T为温度,kB为玻尔兹曼常数,mg为气体分子质量;
3)采用直接模拟蒙特卡洛方法计算飞灰颗粒的碰撞时,在时间步长Δt内求解步骤2)中的颗粒运动方程,首先确定飞灰颗粒的位移且不处理碰撞,接着采用蒙特卡洛方法进行飞灰颗粒碰撞的计算,颗粒i的碰撞概率为Pi:
式中,Pij为颗粒i与颗粒j的碰撞概率,N为计算网格内的模拟颗粒数,n为网格内真实颗粒数,g0为径向分布函数,vij为颗粒i与颗粒j的相对速度;
通过时间计数器法搜索出飞灰颗粒碰撞对,利用颗粒碰撞动力学方程处理碰撞后颗粒的速度:
mi(vi,1-vi,0)=J (8)
mj(vj,1-vj,0)=-J (9)
式中,mi为颗粒i的质量,mj为颗粒j的质量,vi,0和vj,0表示飞灰颗粒i和颗粒j碰撞前的速度,vi,1和vj,1表示飞灰颗粒i和颗粒j碰撞后的速度,J表示碰撞冲量;
4)根据颗粒位移判断颗粒是否与壁面发生碰撞,记录下发生了壁面碰撞的飞灰颗粒,求解飞灰粘附率模型,计算飞灰颗粒在采样面上的沉积;
飞灰颗粒的捕集率fdep:
式中,pi(Tps)是来流飞灰颗粒粘附效率,μref是临界粘度,ps(Ts)是积灰表面已沉积颗粒的粘附效率;
根据飞灰颗粒的捕集率确定颗粒是粘附在采样面上后,利用采样枪提供的壁面热流q和传热方程就可获得更新后的壁面温度Ts:
Ts=T2+q·(δ2/λ+δs/λs) (12)
式中,T1和T2分别是采样枪的内外壁温,δ1和δ2是厚度,λ是采样枪的导热系数,λs是积灰的导热系数,δs是积灰厚度;
5)经过步骤4)后判断累计时间步长与设定时间大小,如果小于设定时间,则更新壁面温度进行下一次时间步长的计算,如果大于设定时间,则流程终止,输出结果,给出炉内飞灰颗粒的沉积速率,并能得到详细的积灰内层飞灰颗粒粒径分布。
本发明具有以下优点及突出性的技术效果:该方法构建的炉内积灰模型能准确预测炉内飞灰颗粒的沉积速率,并给出详细的积灰内层飞灰颗粒粒径分布,成功实现炉内积灰过程的动态预测。
附图说明
图1方法总体流程图。
图2积灰采样枪结构图,图中:1一号热电偶;2-二号热电偶;3-积灰采样面。
图3热电偶读数随时间变化。
图4a、图4b分别为积灰厚度随时间的变化以及粘附颗粒粒径分布随时间的变化。
图5积灰采样面实验结果。
具体实施方式
下面结合附图和实施例对本发明做进一步的说明。本发明提供的一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法,其具体的实现方法如下:
1)作为直接模拟蒙特卡洛计算的输入条件,需要先对采样点附近流场进行计算,采用计算流体动力学方法,其主要控制方程是稳态不可压缩流动纳维叶-斯托克斯方程:
式中,ρ为当地流体的密度,U为当地流体的速度,p为当地流体的压力,μ为动力学粘度,S为应变率张量;
2)飞灰颗粒在外力作用下运动满足牛顿第二运动定律,仅考虑飞灰颗粒的重力、气体施加的曳力以及热泳力,颗粒运动方程可表示为:
Fd=-3πdμ(v-u)f (3)
式中,Fd为流体施加给颗粒的曳力,Fth为热泳力,m为颗粒质量,g为重力加速度,d为颗粒直径,μ为空气动力学粘度系数,v为颗粒运动速度,u为气体速度,f为曳力系数,κtr为气体导热系数,r为颗粒半径,为气体分子的平均运动速度,T为温度,kB为玻尔兹曼常数,mg为气体分子质量;
3)采用直接模拟蒙特卡洛方法计算飞灰颗粒的碰撞时,在时间步长Δt内求解步骤2)中的颗粒运动方程,首先确定飞灰颗粒的位移且不处理碰撞,接着采用蒙特卡洛方法进行飞灰颗粒碰撞的计算,颗粒i的碰撞概率为Pi:
式中,Pij为颗粒i与颗粒j的碰撞概率,N为计算网格内的模拟颗粒数,n为网格内真实颗粒数,g0为径向分布函数,vij为颗粒i与颗粒j的相对速度。
通过时间计数器法搜索出飞灰颗粒碰撞对,利用颗粒碰撞动力学方程处理碰撞后颗粒的速度:
mi(vi,1-vi,0)=J (8)
mj(vj,1-vj,0)=-J (9)
式中,mi为颗粒i的质量,mj为颗粒j的质量,vi,0和vj,0表示飞灰颗粒i和颗粒j碰撞前的速度,vi,1和vj,1表示飞灰颗粒i和颗粒j碰撞后的速度,J表示碰撞冲量;
4)根据颗粒位移判断颗粒是否与壁面发生碰撞,记录下发生了壁面碰撞的飞灰颗粒,求解飞灰粘附率模型,计算飞灰颗粒在采样面上的沉积。
飞灰颗粒的捕集率fdep:
式中,pi(Tps)是来流飞灰颗粒粘附效率,μref是临界粘度,ps(Ts)是积灰表面已沉积颗粒的粘附效率;
根据飞灰颗粒的捕集率确定颗粒是粘附在采样面上后,利用采样枪提供的壁面热流q和传热方程就可获得更新后的壁面温度Ts:
Ts=T2+q·(δ2/λ+δs/λs) (12)
式中,T1和T2分别是采样枪的内外壁温,δ1和δ2是厚度,λ是采样枪的导热系数,λs是积灰的导热系数,δs是积灰厚度;
5)经过步骤4)后判断累计时间步长与设定时间大小,如果小于设定时间,则更新壁面温度进行以一次时间步长的计算,如果大于设定时间,则流程终止,输出结果,给出炉内飞灰颗粒的沉积速率,并能得到详细的积灰内层飞灰颗粒粒径分布;
实施例:
本发明具体的实现方法如下:
以某660MW超临界煤粉锅炉为研究对象,将积灰采样枪伸入炉膛某一层看火孔内,如图2所示。其中采样枪由3个部分组成,枪头是一个12CrMoV材料的圆柱形合金钢,直径D1=90mm,D2=3mm厚度L1=30mm,两个K型热电偶分别埋在厚度L2=28mm以及外表面处,通过两个温度显示器分别记录热电偶读数。实验持续30分钟,图3为两个热电偶记录下的读数,再根据公式(13)获得热流。
采用计算流体动力学方法对采样点附近流场进行计算,获得采样点附近的速度场,然后将图3的热流数据作为直接模拟蒙特卡洛计算的边界条件输入模型,由步骤2)和步骤3)计算飞灰颗粒的碰撞和运动。图4a为模拟计算的积灰厚度随时间变化过程,图4b为粘附在采样面上的飞灰颗粒粒径分布。模拟计算时间为0.01s,积灰形状呈现中间高两边低,积灰厚度最大值4.84×10-8mm出现在长度27.78mm处,平均积灰速率为2.16×10-6m/s,以实验时间30分钟为例,则模拟计算的积灰厚度应该为3.89mm。
图5为积灰采样枪在实验结束后表面所沉积的灰样,通过投影积灰轮廓线可以计算出平均积灰厚度为3.69mm,模拟计算的结果与实验测量值误差为5.14%。实验证明了这种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法的准确性。
Claims (1)
1.一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法,其特征在于所述方法包括如下步骤:
1)采用计算流体动力学方法计算采样点附近流场,获得采样点附近的速度场;
2)采用颗粒运动方程计算飞灰颗粒的重力、气体施加的曳力以及热泳力,其表达式为:
Fd=-3πdμ(v-u)f (2)
式中,Fd为流体施加给颗粒的曳力,Fth为热泳力,m为颗粒质量,g为重力加速度,d为颗粒直径,μ为空气动力学粘度系数,v为颗粒运动速度,u为气体速度,f为曳力系数,κtr为气体导热系数,r为颗粒半径,为气体分子的平均运动速度,T为温度,kB为玻尔兹曼常数,mg为气体分子质量;
3)在时间步长Δt内求解步骤2)中的颗粒运动方程获得颗粒位移,然后采用蒙特卡洛方法进行飞灰颗粒碰撞的计算,颗粒i的碰撞概率为Pi:
式中,Pij为颗粒i与颗粒j的碰撞概率,N为计算网格内的模拟颗粒数,n为网格内真实颗粒数,g0为径向分布函数,vij为颗粒i与颗粒j的相对速度;
利用颗粒碰撞动力学方程计算碰撞后颗粒的速度:
mi(vi,1-vi,0)=J (7)
mj(vj,1-vj,0)=-J (8)
式中,mi为颗粒i的质量,mj为颗粒j的质量,vi,0和vj,0表示飞灰颗粒i和颗粒j碰撞前的速度,vi,1和vj,1表示飞灰颗粒i和颗粒j碰撞后的速度,J表示碰撞冲量;
4)根据颗粒位移判断颗粒是否与壁面发生碰撞,记录下发生了壁面碰撞的飞灰颗粒,求解飞灰粘附率模型,计算飞灰颗粒在采样面上的沉积;
飞灰颗粒的捕集率fdep:
式中,pi(Tps)是来流飞灰颗粒粘附效率,μref是临界粘度,ps(Ts)是积灰表面已沉积颗粒的粘附效率;
利用采样枪提供的壁面热流q和传热方程获得更新后的壁面温度Ts:
Ts=T2+q·(δ2/λ+δs/λs) (11)
式中,T1和T2分别是采样枪的内外壁温,δ1和δ2是厚度,λ是采样枪的导热系数,λs是积灰的导热系数,δs是积灰厚度;
5)经过步骤4)后判断累计时间步长与设定时间大小,如果小于设定时间,则更新壁面温度进行下一次时间步长的计算,如果大于设定时间,则流程终止,输出结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610490783.5A CN106021820B (zh) | 2016-06-28 | 2016-06-28 | 一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610490783.5A CN106021820B (zh) | 2016-06-28 | 2016-06-28 | 一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106021820A true CN106021820A (zh) | 2016-10-12 |
CN106021820B CN106021820B (zh) | 2019-02-22 |
Family
ID=57084076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610490783.5A Expired - Fee Related CN106021820B (zh) | 2016-06-28 | 2016-06-28 | 一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106021820B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885947A (zh) * | 2019-02-27 | 2019-06-14 | 重庆大学 | 基于cfd技术的高温活化粉煤灰显热回收模拟分析方法 |
CN113221423A (zh) * | 2021-05-21 | 2021-08-06 | 华北电力大学(保定) | 确定污秽颗粒是否粘附在电力设备外表绝缘面上的方法 |
CN114814076A (zh) * | 2022-04-12 | 2022-07-29 | 华中科技大学 | 一种飞灰颗粒沾污平均界面能的计算方法及其测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1715339A2 (de) * | 2005-03-30 | 2006-10-25 | Technische Universität Wien | Verfahren zur Ermittlung der Anteile biogener und fossiler Energieträger sowie biogener und fossiler Kohlendioxidemissionen beim Betrieb von Verbrennungsanlagen |
CN102354336A (zh) * | 2011-10-12 | 2012-02-15 | 西安交通大学 | 一种对多晶硅cvd反应器辐射热损失的预估方法 |
-
2016
- 2016-06-28 CN CN201610490783.5A patent/CN106021820B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1715339A2 (de) * | 2005-03-30 | 2006-10-25 | Technische Universität Wien | Verfahren zur Ermittlung der Anteile biogener und fossiler Energieträger sowie biogener und fossiler Kohlendioxidemissionen beim Betrieb von Verbrennungsanlagen |
CN102354336A (zh) * | 2011-10-12 | 2012-02-15 | 西安交通大学 | 一种对多晶硅cvd反应器辐射热损失的预估方法 |
Non-Patent Citations (1)
Title |
---|
J.R. FAN 等: "Simulation of ash deposit in a pulverized coal-fired boiler", 《FUEL》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109885947A (zh) * | 2019-02-27 | 2019-06-14 | 重庆大学 | 基于cfd技术的高温活化粉煤灰显热回收模拟分析方法 |
CN109885947B (zh) * | 2019-02-27 | 2020-10-20 | 重庆大学 | 基于cfd技术的高温活化粉煤灰显热回收模拟分析方法 |
CN113221423A (zh) * | 2021-05-21 | 2021-08-06 | 华北电力大学(保定) | 确定污秽颗粒是否粘附在电力设备外表绝缘面上的方法 |
CN114814076A (zh) * | 2022-04-12 | 2022-07-29 | 华中科技大学 | 一种飞灰颗粒沾污平均界面能的计算方法及其测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106021820B (zh) | 2019-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106021820A (zh) | 一种基于直接模拟蒙特卡洛的炉内积灰模型的建模方法 | |
Rajneesh et al. | Prediction of raceway size in blast furnace from two dimensional experimental correlations | |
Kasper et al. | Numerical modeling and simulation of particulate fouling of structured heat transfer surfaces using a multiphase Euler-Lagrange approach | |
Bu et al. | Numerical analysis of ABS deposition and corrosion on a rotary air preheater | |
CN105468844A (zh) | 管道内水-气耦合瞬变流的模拟方法 | |
CN103400035B (zh) | 一种高可信度快速预测飞行器滚转动导数的方法 | |
CN103473396B (zh) | 一种运载火箭低温贮箱换热过程仿真方法 | |
Abd-Elhady et al. | Modelling the impaction of a micron particle with a powdery layer | |
CN115587551B (zh) | 一种高超声速飞行器防热结构烧蚀行为多尺度预测方法 | |
CN104881588A (zh) | 铸锭宏观偏析数值模拟方法 | |
CN104634708A (zh) | 基于计算流体力学流化床内颗粒密度和粒径分布预测方法 | |
Usov et al. | SAFR/V1 (EVKLID/V2 integral code module) aided simulation of melt movement along the surface of a fuel element in a fast reactor during a serious accident | |
CN108197070A (zh) | 天然气不完全堵塞管道数值模拟方法 | |
CN107122516A (zh) | 一种灭火系统自段管道的沿程压力损失确定方法 | |
Peng et al. | A novel nonsmooth dynamics method for multibody systems with friction and impact based on the symplectic discrete format | |
Ji et al. | Critical sinking of hydrophobic micron particles | |
CN103279645A (zh) | 基于gpu并行运算的碳纳米管分子动力学仿真方法 | |
Chen et al. | Investigation of particle distribution after the energetic module broken in the ignition process of gun | |
Florio | Estimation of particle impact based erosion using a coupled direct particle—Compressible gas computational fluid dynamics model | |
CN104036084B (zh) | 塔式太阳能热电站管状接收器的分布参数建模方法 | |
CN113591345A (zh) | 一种基于广义黎曼解法器的爆炸反应流高精度预测方法 | |
CN109753679B (zh) | 一种空预器积灰堵塞可视化监控方法及系统 | |
Liu et al. | Numerical simulation and mathematical models of ash deposition behavior considering particle properties and operating conditions | |
Chilamkurti et al. | Characterizing particle-wall contact behavior and fluctuations in gravity-driven dense granular flows in cylindrical tubes using DEM | |
Lutro | The effect of thermophoresis on the particle deposition on a cylinder |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190222 Termination date: 20190628 |