CN110472270A - 一种钠冷堆非能动余热排出***的建模优化方法 - Google Patents
一种钠冷堆非能动余热排出***的建模优化方法 Download PDFInfo
- Publication number
- CN110472270A CN110472270A CN201910579525.8A CN201910579525A CN110472270A CN 110472270 A CN110472270 A CN 110472270A CN 201910579525 A CN201910579525 A CN 201910579525A CN 110472270 A CN110472270 A CN 110472270A
- Authority
- CN
- China
- Prior art keywords
- air
- sodium
- heat exchanger
- heat
- tube
- 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
Landscapes
- Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
Abstract
本发明涉及一种钠冷堆非能动余热排出***的建模优化方法,属于钠冷堆余热排出***的技术领域,初步设定热力参数、空气热交换器结构参数、拔风烟囱结构参数等;预设空气侧出口温度,设定传热量等相关计算判断空气侧流动平衡;预设钠侧进口温度与钠侧出口温度设定传热量等相关计算判断钠侧流动平衡;空气侧流动平衡且钠侧流动平衡后经相关计算判断***能量是否平衡,并输出相关设计参数。本发明模型便于电脑语言编程并查找编写错误,可准确高效的完成设备及工艺流程设计工作,并可根据设计需要添加或删减相关模块。本发明的建模优化方法为建立准确的非能动余热排出***模型创造了前提,并显著提高参数估计的精度和速度。
Description
技术领域
本发明涉及钠冷堆余热排出***的技术领域,具体是涉及一种钠冷堆非能动余热排出***的建模优化方法。
背景技术
日本福岛核泄漏事故后,核电***的安全性和可靠性越来越受到人们的重视。非能动安全作为***核技术的关键考虑因素,可有效降低堆芯融毁的概率。钠冷堆相对于传统的压水堆可将铀资源的利用率最大化,同时热响应时间长,冷却剂难发生沸腾现象,反应堆有着更好的安全特性。
目前,对非能动余热排出***进行高效设计并优化尚缺乏完善的解决方法。
发明内容
为了解决非能动余热排出***的高效设计和技术优化的技术问题,本发明提供一种钠冷堆非能动余热排出***的建模优化方法。
为了实现本发明的目的,本发明采用了以下技术方案:一种钠冷堆非能动余热排出***的建模优化方法,钠冷堆非能动余热排出***包括水平高度依次升高的堆芯、事故热交换器、空气热交换器;所述空气热交换器设于拔风烟囱的内部;所述事故热交换器的壳程为堆芯内待冷却的热介质,管程为冷却介质;所述事故热交换器的冷却介质与空气热交换器的加热介质为形成中间自然循环回路的液钠;所述拔风烟囱设有入风口,所述入风口设于空气热交换器的底部;所述空气热交换器的冷却介质为空气,方法包括以下步骤:
S1、特定物性参数为特定条件下的固有参数;初步设定热力参数、空气热交换器结构参数、拔风烟囱结构参数,及事故热交换器液钠出口与空气热交换器液钠出口之间的高度差H2;
所述热力参数为空气热交换器中的各参数,分别包括空气侧进口温度Te、设定传热量Q、钠侧进口温度Ts,i与钠侧出口温度Ts,o;
所述空气热交换器为液钠走管程,空气走壳程的翅片管式换热器;翅片管式换热器的换热管为翅片管,所述翅片管包括基管及设置在基管上的翅片;空气热交换器结构参数包括基管结构参数和翅片结构参数;所述基管结构参数包括基管总管数ns、基管外直径do、基管内直径di、基管的总管长Ls、基管的排列方式与排列的间距,及根据以上数据得到的沿空气流通方向基管的管排数Na、单位长度基管外表面积Ao*,基管内流通截面积Ss,空气热交换器翅片管间最小流通面积Sa,min;所述翅片结构参数包括翅片的翅高Hf、节距y,及根据以上数据得到的单位长度翅片表面积Af与翅化比β;
所述液钠的中间自然循环回路中,钠侧弯头总数为Ns,
所述拔风烟囱结构参数包括空气换热器的壳程空气出口与拔风烟囱出口高度差H1、烟道内直径Di、烟道流通截面积Sc、烟道弯头数Nc;
所述特定物性参数包括钠侧比热cps与空气侧比热cpa;不同温度对应下的空气密度ρa;不同温度对应下的液钠密度ρs;热耗散q、直管阻力系数f、局部阻力系数ξ、进出口阻力系数ε、污垢热阻r、钠侧沿管程的雷诺数Re、导热系数λ、对数平均温差ΔTm、粘度μ、贝克来数Pes与对数平均温差修正系数φ;
S2、预设空气热交换器的空气侧出口温度Ta,o1,根据设定传热量Q确定空气质量流量ma;计算得到空气侧浮升力Fa和空气侧总阻力降ΔPa后;计算空气侧浮升力Fa与空气侧总阻力降ΔPa的相对误差当相对误差大于1%时,则返回调整预设空气侧出口温度Ta,o1;当相对误差小于等于1%时,则判断空气侧流动平衡;
预设钠侧进口温度Ts,i与钠侧出口温度Ts,o,由设定传热量Q确定液钠的质量流量ms,计算得到钠侧提升力Fs与钠侧总阻力降ΔPs;计算钠侧提升力Fs与钠侧总阻力降ΔPs的相对误差当相对误差大于1%时,则返回调整钠侧进口温度Ts,i与钠侧出口温度Ts,o;当相对误差小于等于1%时,则判断钠侧流动平衡;
S3、待S2中的空气侧流动平衡且钠侧流动平衡后,计算得到空气侧传热系数ha和钠侧传热系数hs;进而再求得空气热交换器总传热系数K、对数平均温差ΔTm与总传热量Qs;再计算总传热量Qs与设定传热量Q的相对误差当相对误差大于1%时,则判断***能量不平衡,并返回步骤S2重新计算;当相对误差小于等于1%时,则判断***能量平衡,并输出空气质量流量ma、钠侧质量流量ms、空气侧进口温度Te、空气侧出口温度Ta,o1、钠侧进口温度Ts,i、钠侧出口温度Ts,o、空气热交换器总传热系数K、总传热量Qs、总传热量Qs与设定传热量Q的相对误差
钠侧比热cps是指管路翅片管中液钠的比热容;与空气侧比热cpa是指管路翅片管中空气的比热容。
空气热换热器壳程空气的入口与拔风烟囱的入风口可以是在同一水平高度,也可以是空气热换热器壳程空气的入口的水平高度高于拔风烟囱的入风口。
进一步地,步骤S2所述空气侧浮升力Fa和空气侧总阻力降ΔPa的计算步骤如下:
空气热交换器中换热的能量平衡方程如下式(1),钠侧能量损失量与空气侧能量增加量均分别等于设定传热量Q:
mscps(Ts,i-Ts,o)=macpa(Ta,o1-Te)=Q (1)
其中,Ts,i、Ts,o分别为空气热交换器的钠侧进口温度、钠侧出口温度;Te、Ta,o1分别为空气热交换器的空气侧进口温度、空气侧出口温度;ms、ma分别为空气热交换器的钠侧质量流量、空气质量流量;cps、cpa分别为空气热交换器中的钠侧比热、空气侧比热;
拔风烟囱内的空气侧浮升力Fa的计算公式如下式(2):
其中,H1为空气换热器的壳程空气出口与拔风烟囱出口高度差;ρe为空气热交换器的空气侧进口温度Te下的空气密度;ρa,om为对环境散热后拔风烟囱内平均温度Ta,om下的空气密度,式(2)中的ρa,om与式(3)中的Ta,om对应,并满足公式Ta,om由式(3)求得,其中,Ta,o1为空气热交换器的空气侧出口温度,qc为拔风烟囱与外部环境之间的热耗散,ma为空气热交换器的空气质量流量,cpa为空气侧比热;g为重力加速度;dh为高度h的微分;
空气侧总阻力降ΔPa的计算公式如下式(4):
其中,fa为翅片管外阻力系数,Na为沿空气流通方向基管的管排数,ma为空气质量流量,Sa,min为空气热交换器翅片管间最小流通面积,ρa,m为空气热交换器中的空气平均密度;fc为烟囱内沿程阻力系数,Sc为烟道流通截面积,ρa,om为对环境散热后拔风烟囱内平均温度Ta,om下的空气密度,Di为拔风烟囱烟道内直径;Nc为拔风烟囱烟道弯头数,ξc为拔风烟囱局部阻力系数,εc为拔风烟囱进出口阻力系数;
步骤S2所述的钠侧提升力Fs与钠侧总阻力降ΔPs计算步骤如下:忽略事故热交换器液钠出口与空气热交换器液钠入口之间管路的热量损耗,事故热交换器液钠出口温度等于空气热交换器液钠入口温度;事故热交换器液钠出口与经过换热后的空气热交换器液钠出口之间存在温度差,并存在密度差,两者之间存在高度差H2;事故热交换器液钠出口与空气热交换器液钠出口之间钠侧提升力Fs的计算公式如下式(5):
其中,H2为事故热交换器液钠出口与空气热交换器液钠出口之间的高度差;ρs,o、ρs,i分别为空气热交换器的出口液钠密度、进口液钠密度;g为重力加速度;dh为高度h的微分;
钠侧总阻力降ΔPs的计算公式如下式(6):
其中,式(6)中的fs为钠侧沿程阻力系数,钠侧沿程阻力系数fs与雷诺数Res的关系如式(7);ms为钠侧质量流量,Ss为基管内流通截面积;ρs,m为自然循环回路中液钠的平均密度;Lsur为自然循环回路中除基管外的剩余段管长;Ls为基管的总管长;di为基管内直径;Ns为钠侧弯头总数;ξs为钠侧局部阻力系数。
进一步地,空气热交换器中的翅片管为三角形错列排布;所述壳程空气流动方向与管程液钠流动方向垂直;空气为竖向流通,液钠为水平方向流动;
步骤S3所述空气对流传热系数ha的计算公式如下式(8-1):
其中:λa为空气导热系数;do为基管外直径;ma为空气质量流量;Sa,min为空气热交换器翅片管间最小流通面积;μa为空气粘度;cpa为空气侧比热;Ao为单位长度翅片管总外表面积;为单位长度基管外表面积;
液钠冷却过程,于均匀热流边界条件下,液钠流动的对流传热系数hs的计算公式如下式(9):
其中:λs为液钠的导热系数;di为基管内直径;Pes为贝克来数;
根据式(8-1)和式(9),空气热交换器总传热系数K的计算公式如下式(10):
其中:β为翅化比,ha为空气侧对流传热系数;ra为空气侧污垢热阻;Ao为单位长度翅片管总外表面积;λw为管材导热系数;do为基管外直径;di基管内直径;rs为钠侧污垢热阻;η为翅片壁面总效率;hs为钠侧对流传热系数;管材导热系数是指翅片管选用材料的传导效率,基管和翅片为相同材质,具有相同的传导效率。
***总传热量Qs的计算公式如下式(11):
Qs=φKnsLsAoΔTm (11)
其中,φ为对数平均温差修正系数;K为空气热交换器总传热系数;ns为基管总管数;Ls为基管的总管长;Ao为单位长度翅片管总外表面积;ΔTm为对数平均温差。
进一步地,所述翅片管为正方形直列排布;所述壳程空气流动方向与管程液钠流动方向垂直;空气为竖向流通,液钠为水平方向流动;
步骤S3所述空气对流传热系数ha的计算公式如下式(8-2):
其中:λa为空气导热系数;do为基管外直径;ma为空气质量流量;Sa,min为空气热交换器翅片管间最小流通面积;μa为空气粘度;cpa为空气侧比热;Af单位长度翅片表面积;为单位长度基管外表面积;
液钠冷却过程,于均匀热流边界条件下,液钠流动的对流传热系数hs的计算公式如下式(9):
其中:λs为液钠的导热系数;di为基管内直径;Pes为贝克来数;
根据式(8-2)和式(9),空气热交换器总传热系数K的计算公式如下式(10):
其中:β为翅化比,ha为空气侧对流传热系数;ra为空气侧污垢热阻;Ao为单位长度翅片管总外表面积;λw为管材导热系数;do为基管外直径;di基管内直径;rs为钠侧污垢热阻;η为翅片壁面总效率;hs为钠侧对流传热系数;
***总传热量Qs的计算公式如下式(11):
Qs=φKnsLsAoΔTm (11)
其中,φ为对数平均温差修正系数;K为空气热交换器总传热系数;ns为基管总管数;Ls为基管的总管长;Ao为单位长度翅片管总外表面积;ΔTm为对数平均温差。
进一步地,当翅片管为三角形错列排布时,公式(4)中的翅片管外阻力系数fa为:
当翅片管为正方形直列排布时,公式(4)中的翅片管外阻力系数fa为:
其中,ma为空气质量流量,Sa,min为空气热交换器翅片管间最小流通面积;do为基管外直径;μa为空气粘度;Pt为管束横向管间距;Pl为管束纵向管间距;y为翅片管节距;Hf为翅高。其中,以翅片管为正方形直列排布的情况考虑;空气是自低处向高处流动,流通方向为竖直方向,其管束纵向管间距即是为管束沿空气流通方向上的相邻两管束自间的距离;液钠的流通在空气热交换器的管程内,其是水平流动的,管束横向管间距是指同一水平高度上的相邻两管束之间的距离。
进一步的,构建物性参数计算模块、空气侧流动平衡计算模块、液态钠侧流动平衡计算模块与***能量平衡计算模块;所述特定物性参数存于物性参数计算模块内;步骤S2从预设空气侧出口温度至判断空气侧流动平衡的过程属于空气侧流动平衡计算模块;步骤2从预设钠侧进口温度与钠侧出口温度至判断钠侧流动平衡的过程属于液态钠侧流动平衡计算模块;步骤S3的过程属于***能量平衡计算模块;所述物性参数计算模块分别与空气侧流动平衡计算模块、液态钠侧流动平衡计算模块连接;所述空气侧流动平衡计算模块判断空气侧流动平衡后进入***能量平衡计算模块;液态钠侧流动平衡计算模块判断钠侧流动平衡后进入***能量平衡计算模块。
本发明的有益效果在于:本发明钠冷堆非能动余热排出***主要是由空气热交换器、事故热交换器、拔风烟囱组成,拔风烟囱底部设有风门(风门为拔风烟囱的入风口)。***包括三个自然循环部分,热钠池中事故热交换器和堆芯构成了一回路自然循环,事故热交换器和空气热交换器之间的液钠循环构成了中间回路自然循环,位于拔风烟囱中的空气热交换器和外部环境构成了空气冷却自然循环。事故发生时,一回路自然循环将堆芯热量转移至事故热交换器中,热量通过事故热交换器经中间回路传递至空气热交换器,最终依靠拔风烟囱中空气冷却循环将热量转移至大气中。本发明主要基于钠冷堆余热排出***的工艺过程进行模块化分析,根据设计传热量与传热元件的物理特性参量、自然动力源高度和内直径等参数,即可获得***各状态点热力参数。计算结果为钠冷反应堆余热排出***的安全设计提供理论依据和实用价值。
数学模型的模块化分析方法依据实现功能的差异将设计模型分为不同模块,便于电脑语言编程并查找编写错误,可准确高效的完成设备及工艺流程设计工作,并可根据设计需要添加或删减相关模块。选择合理的模块化结构是建立准确的非能动余热排出***模型的重要前提;而恰当地确定模型参数的范围,可以显著提高参数估计的精度和速度。
并且,本发明在构建的排出***模型的基础上还验证了传热元件不同的物理特性参量、自然动力源高度和内直径参数对***模型数据选取的影响;本发明为非能动自然循环***,总驱动压等于公式(2)求得的空气侧浮升力Fa;结果表明,1)、***的总传热系数和总驱动压随着环境温度升高而逐渐减小,且环境温度对驱动压力的影响更为明显。2)、随着拔风烟囱高度增加,***总传热系数和总驱动压均增大,且增大趋势不断变缓,存在设计最优值;3)、翅片管翅高减小,***总传热系数及单位压降传热系数大幅增加,对***的传热性能影响明显;4)、拔风烟囱内直径增加,***总传热系数增大,总驱动压减小,且变化趋势不断变缓,内直径增加对***传热和流动有积极影响。5)、翅片厚度增加,***总传热系数及单位压降传热系数大幅增加,对***的传热性能影响明显。以上趋势的总结为之后实际建筑的参数设置有着指导意义。
附图说明
图1为本发明钠冷堆余热排出***模型的结构示意图;
图2为钠冷堆余热排出***模型计算流程图;
图3(a)-3(c)为不同环境温度下的影响实验线形图;
图4(a)-4(d)为不同烟囱高度及翅片高度的影响实验线形图;
图5(a)-5(d)为不同烟囱内直径及翅片厚度的影响实验线形图。
具体实施方式
以下结合附图和具体实施方式对本发明技术方案作进一步详细描述。
图1为钠冷堆余热排出***模型示意图,所述钠冷堆非能动余热排出***包括空气热交换器、拔风烟囱、事故热交换器,及堆芯与事故热交换器的作用形成的对液钠升温提升的热钠池。
表1各公式中字母含义
表2空气物性与空气侧温度Ta(K)的关系
表3液钠物性与钠侧温度Ts(K)的关系
实施例1
非能动余热排出***初始设定传热量与传热元件的物理特性参量、自然动力源烟囱的高度和内直径参数如下表4所示;
表4
翅片管结构
设备结构
如图2钠冷堆余热排出***模型计算流程图所示的,经以下步骤进行计算:一种钠冷堆非能动余热排出***的建模优化方法,钠冷堆非能动余热排出***包括水平高度依次升高的堆芯、事故热交换器、空气热交换器;所述空气热交换器设于拔风烟囱的内部;所述事故热交换器的壳程为堆芯内待冷却的热介质,管程为冷却介质;所述事故热交换器的冷却介质与空气热交换器的加热介质为形成中间自然循环回路的液钠;所述拔风烟囱设有入风口,所述入风口设于空气热交换器的底部;所述空气热交换器的冷却介质为空气,方法包括以下步骤:
S1、特定物性参数为特定条件下的固有参数;初步设定热力参数、空气热交换器结构参数、拔风烟囱结构参数,及事故热交换器液钠出口与空气热交换器液钠出口之间的高度差H2;
所述热力参数为空气热交换器中的各参数,分别包括空气侧进口温度Te、设定传热量Q、钠侧进口温度Ts,i与钠侧出口温度Ts,o;
所述空气热交换器为液钠走管程,空气走壳程的翅片管式换热器;翅片管式换热器的换热管为翅片管,所述翅片管包括基管及设置在基管上的翅片;空气热交换器结构参数包括基管结构参数和翅片结构参数;所述基管结构参数包括基管总管数ns、基管外直径do、基管内直径di、基管的总管长Ls、基管的排列方式与排列的间距,及根据以上数据得到的沿空气流通方向基管的管排数Na、单位长度基管外表面积Ao *,基管内流通截面积Ss,空气热交换器翅片管间最小流通面积Sa,min;所述翅片结构参数包括翅片的翅高Hf、节距y,及根据以上数据得到的单位长度翅片表面积Af与翅化比β;
所述液钠的中间自然循环回路中,钠侧弯头总数为Ns,
所述拔风烟囱结构参数包括空气换热器的壳程空气出口与拔风烟囱出口高度差H1、烟道内直径Di、烟道流通截面积Sc、烟道弯头数Nc;
所述特定物性参数包括钠侧比热cps与空气侧比热cpa;不同温度对应下的空气密度ρa;不同温度对应下的液钠密度ρs;热耗散q、直管阻力系数f、局部阻力系数ξ、进出口阻力系数ε、污垢热阻r、钠侧沿管程的雷诺数Re、导热系数λ、对数平均温差ΔTm、粘度μ、贝克来数Pes与对数平均温差修正系数φ;
S2、预设空气热交换器的空气侧出口温度Ta,o1,根据设定传热量Q确定空气质量流量ma;计算得到空气侧浮升力Fa和空气侧总阻力降ΔPa后;计算空气侧浮升力Fa与空气侧总阻力降ΔPa的相对误差当相对误差大于1%时,则返回调整预设空气侧出口温度Ta,o1;当相对误差小于等于1%时,则判断空气侧流动平衡;
预设钠侧进口温度Ts,i与钠侧出口温度Ts,o,由设定传热量Q确定液钠的质量流量ms,计算得到钠侧提升力Fs与钠侧总阻力降ΔPs;计算钠侧提升力Fs与钠侧总阻力降ΔPs的相对误差当相对误差大于1%时,则返回调整钠侧进口温度Ts,i与钠侧出口温度Ts,o;当相对误差小于等于1%时,则判断钠侧流动平衡;
S3、待S2中的空气侧流动平衡且钠侧流动平衡后,计算得到空气侧传热系数ha和钠侧传热系数hs;进而再求得空气热交换器总传热系数K、对数平均温差ΔTm与总传热量Qs;再计算总传热量Qs与设定传热量Q的相对误差当相对误差大于1%时,则判断***能量不平衡,并返回步骤S2重新计算;当相对误差小于等于1%时,则判断***能量平衡,并输出空气质量流量ma、钠侧质量流量ms、空气侧进口温度Te、空气侧出口温度Ta,o1、钠侧进口温度Ts,i、钠侧出口温度Ts,o、空气热交换器总传热系数K、总传热量Qs、总传热量Qs与设定传热量Q的相对误差
构建物性参数计算模块、空气侧流动平衡计算模块、液态钠侧流动平衡计算模块与***能量平衡计算模块;所述特定物性参数存于物性参数计算模块内;步骤S2从预设空气侧出口温度至判断空气侧流动平衡的过程属于空气侧流动平衡计算模块;步骤2从预设钠侧进口温度与钠侧出口温度至判断钠侧流动平衡的过程属于液态钠侧流动平衡计算模块;步骤S3的过程属于***能量平衡计算模块。
步骤S2所述空气侧浮升力Fa和空气侧总阻力降ΔPa的计算步骤如下:
空气热交换器中换热的能量平衡方程如下式(1),钠侧能量损失量与空气侧能量增加量均分别等于设定传热量Q:
mscps(Ts,i-Ts,o)=macpa(Ta,o1-Te)=Q (1)
其中,Ts,i、Ts,o分别为空气热交换器的钠侧进口温度、钠侧出口温度;Te、Ta,o1分别为空气热交换器的空气侧进口温度、空气侧出口温度;ms、ma分别为空气热交换器的钠侧质量流量、空气质量流量;cps、cpa分别为空气热交换器中的钠侧比热、空气侧比热;
拔风烟囱内的空气侧浮升力Fa的计算公式如下式(2):
其中,H1为空气换热器的壳程空气出口与拔风烟囱出口高度差;ρe为空气热交换器的空气侧进口温度Te下的空气密度;ρa,om为对环境散热后拔风烟囱内平均温度Ta,om下的空气密度,式(2)中的ρa,om与式(3)中的Ta,om对应,并满足公式Ta,om由式(3)求得,其中,Ta,o1为空气热交换器的空气侧出口温度,qc为拔风烟囱与外部环境之间的热耗散,ma为空气热交换器的空气质量流量,cpa为空气侧比热;g为重力加速度;dh为高度h的微分;
空气侧总阻力降ΔPa的计算公式如下式(4):
其中,fa为翅片管外阻力系数,Na为沿空气流通方向基管的管排数,ma为空气质量流量,Sa,min为空气热交换器翅片管间最小流通面积,ρa,m为空气热交换器中的空气平均密度;fc为烟囱内沿程阻力系数,Sc为烟道流通截面积,ρa,om为对环境散热后拔风烟囱内平均温度Ta,om下的空气密度,Di为拔风烟囱烟道内直径;Nc为拔风烟囱烟道弯头数,ξc为拔风烟囱局部阻力系数,εc为拔风烟囱进出口阻力系数;
当翅片管为三角形错列排布时,公式(4)中的翅片管外阻力系数fa为:
当翅片管为正方形直列排布时,公式(4)中的翅片管外阻力系数fa为:
其中,ma为空气质量流量,Sa,min为空气热交换器翅片管间最小流通面积;do为基管外直径;μa为空气粘度;Pt为管束横向管间距;Pl为管束纵向管间距;y为翅片管节距;Hf为翅高。
步骤S2所述的钠侧提升力Fs与钠侧总阻力降ΔPs计算步骤如下:忽略事故热交换器液钠出口与空气热交换器液钠入口之间管路的热量损耗,事故热交换器液钠出口温度等于空气热交换器液钠入口温度;事故热交换器液钠出口与经过换热后的空气热交换器液钠出口之间存在温度差,并存在密度差,两者之间存在高度差H2;事故热交换器液钠出口与空气热交换器液钠出口之间钠侧提升力Fs的计算公式如下式(5):
其中,H2为事故热交换器液钠出口与空气热交换器液钠出口之间的高度差;ρs,o、ρs,i分别为空气热交换器的出口液钠密度、进口液钠密度;g为重力加速度;dh为高度h的微分;
钠侧总阻力降ΔPs的计算公式如下式(6):
其中,式(6)中的fs为钠侧沿程阻力系数,钠侧沿程阻力系数fs与雷诺数Res的关系如式(7);ms为钠侧质量流量,Ss为基管内流通截面积;ρs,m为自然循环回路中液钠的平均密度;Lsur为自然循环回路中除基管外的剩余段管长;Ls为基管的总管长;di为基管内直径;Ns为钠侧弯头总数;ξs为钠侧局部阻力系数。
进一步地,所述壳程空气流动方向与管程液钠流动方向垂直;空气为竖向流通,液钠为水平方向流动;
所述壳程空气流动方向与管程液钠流动方向垂直;空气为竖向流通,液钠为水平方向流动;
空气热交换器中的翅片管为三角形错列排布;步骤S3所述空气对流传热系数ha的计算公式如下式(8-1):
所述翅片管为正方形直列排布;步骤S3所述空气对流传热系数ha的计算公式如下式(8-2):
其中:λa为空气导热系数;do为基管外直径;ma为空气质量流量;Sa,min为空气热交换器翅片管间最小流通面积;μa为空气粘度;cpa为空气侧比热;Af单位长度翅片表面积;Ao为单位长度翅片管总外表面积;为单位长度基管外表面积;
液钠冷却过程,于均匀热流边界条件下,液钠流动的对流传热系数hs的计算公式如下式(9):
其中:λs为液钠的导热系数;di为基管内直径;Pes为贝克来数;
根据式(8-1)和式(9),或是式(8-2)和式(9),空气热交换器总传热系数K的计算公式如下式(10):
其中:β为翅化比,ha为空气侧对流传热系数;ra为空气侧污垢热阻;Ao为单位长度翅片管总外表面积;λw为管材导热系数;do为基管外直径;di基管内直径;rs为钠侧污垢热阻;η为翅片壁面总效率;hs为钠侧对流传热系数;
***总传热量的计算公式如下式(11):
Qs=φKnsLsAoΔTm (11)
其中,φ为对数平均温差修正系数;K为空气热交换器总传热系数;ns为基管总管数;Ls为基管的总管长;Ao为单位长度翅片管总外表面积;ΔTm为对数平均温差。
对本发明所构建***的合理性与有效性进行验证,以下分别对不同的空气侧进口温度、换热器的翅片参数和烟囱参数进行测试;表4是对应的当翅片管为三角形错列排布时,余热排出***中的设计参数、换热管结构设计参数以及设备机构参数,其本验证实验为单一因素考察实验,除变量外的其他参数均同表4;表4中拔风烟囱高度近似等于空气换热器的壳程空气出口与拔风烟囱出口高度差H1;事故热交换器与空气热交换器高度差近似等于事故热交换器液钠出口与空气热交换器液钠出口之间的高度差H2。
环境温度的影响:环境温度即是空气侧进口温度;预设不同空气侧进口温度,并分别进行***热力参数计算结果,不同环境温度下***热力参数计算结果如下表5和图3(a)-3(c)所示。本发明为非能动自然循环***,总驱动压等于公式(2)求得的空气侧浮升力Fa;结果表明,***的总传热系数和总驱动压随着环境温度升高逐渐减小,且环境温度的变化对驱动压力的影响更为明显。如图3所示的,***的总驱动压和总传热系数随着环境温度升高逐渐减小,且环境温度对驱动压力的影响更为明显。
表5不同环境温度下***热力参数计算结果
拔风烟囱高度及翅片高度的影响:于相同翅片高度,不同烟囱高度下计算***热力参数,计算结果如表6-1、表6-2、表6-3与图4(a)-4(d)所示,由表中数据可知,翅片高度固定时,随着拔风烟囱高度增加,***总驱动压和总传热系数均增大,且增大趋势不断变缓,则可预料其总传热系数必然会存在设计最优值;另外,随着翅片高度减小,***总传热系数及单位压降传热系数大幅增加,对***的传热性能也有较高的影响。结果表明,拔风烟囱高度增加,***总传热系数和总驱动压均增大,且增大趋势不断变缓,存在设计最优值;翅片高度减小,***总传热系数及单位压降传热系数大幅增加,对***的传热性能影响明显。
表6-1翅片高度9mm时不同烟囱高度***热力参数
表6-2翅片高度11mm时不同烟囱高度***热力参数
表6-3翅片高度13mm时不同烟囱高度***热力参数
烟囱内直径及翅片厚度的影响:于相同翅片厚度时,不同烟囱内直径下***热力参数的测试,测得结果如下表7-1、表7-2、表7-3与图5(a)-5(d)所示。由表中数据可以明确得出,随着拔风烟囱内直径增加,***总传热系数增大,总驱动压减小,且变化趋势不断变缓,内直径增加对***传热和流动有积极影响。随着翅片厚度增加,***总传热系数及单位压降传热系数大幅增加,对***的传热性能影响明显。
表7-1翅片厚度0.4mm时不同烟囱内直径***热力参数
表7-2翅片厚度0.6mm时不同烟囱内直径***热力参数
表7-3翅片厚度0.8mm时不同烟囱内直径***热力参数
由上测试结果可以看出,相关参数的变化具有规律可寻,本发明的模拟***也具有较高的实用性,可以为实际的核电站施工建筑的各种结构参数进行指导,可以提高工程效率,减少风险,并且本发明的模拟***具有操作简便、计算速度快的优点。
以上实施方式仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施方式对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施方式所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施方式技术方案的精神和范围。
Claims (6)
1.一种钠冷堆非能动余热排出***的建模优化方法,其特征在于,钠冷堆非能动余热排出***包括水平高度依次升高的堆芯、事故热交换器、空气热交换器;所述空气热交换器设于拔风烟囱的内部;所述事故热交换器的壳程为堆芯内待冷却的热介质,管程为冷却介质;所述事故热交换器的冷却介质与空气热交换器的加热介质为形成中间自然循环回路的液钠;所述拔风烟囱设有入风口,所述入风口设于空气热交换器的底部;所述空气热交换器的冷却介质为空气;该方法包括以下步骤:
S1、特定物性参数为特定条件下的固有参数;初步设定热力参数、空气热交换器结构参数、拔风烟囱结构参数,及事故热交换器液钠出口与空气热交换器液钠出口之间的高度差H2;
所述热力参数为空气热交换器中的各参数,分别包括空气侧进口温度Te、设定传热量Q、钠侧进口温度Ts,i与钠侧出口温度Ts,o;
所述空气热交换器为液钠走管程,空气走壳程的翅片管式换热器;翅片管式换热器的换热管为翅片管,所述翅片管包括基管及设置在基管上的翅片;空气热交换器结构参数包括基管结构参数和翅片结构参数;所述基管结构参数包括基管总管数ns、基管外直径do、基管内直径di、基管的总管长Ls、基管的排列方式与排列的间距,及根据以上数据得到的沿空气流通方向基管的管排数Na、单位长度基管外表面积Ao *,基管内流通截面积Ss,空气热交换器翅片管间最小流通面积Sa,min;所述翅片结构参数包括翅片的翅高Hf、节距y,及根据以上数据得到的单位长度翅片表面积Af与翅化比β;
所述液钠的中间自然循环回路中,钠侧弯头总数为Ns,
所述拔风烟囱结构参数包括空气换热器的壳程空气出口与拔风烟囱出口高度差H1、烟道内直径Di、烟道流通截面积Sc、烟道弯头数Nc;
所述特定物性参数包括钠侧比热cps与空气侧比热cpa;不同温度对应下的空气密度ρa;不同温度对应下的液钠密度ρs;热耗散q、直管阻力系数f、局部阻力系数ξ、进出口阻力系数ε、污垢热阻r、钠侧沿管程的雷诺数Re、导热系数λ、对数平均温差ΔTm、粘度μ、贝克来数Pes与对数平均温差修正系数φ;
S2、预设空气热交换器的空气侧出口温度Ta,o1,根据设定传热量Q确定空气质量流量ma;计算得到空气侧浮升力Fa和空气侧总阻力降ΔPa后;计算空气侧浮升力Fa与空气侧总阻力降ΔPa的相对误差当相对误差大于1%时,则返回调整预设空气侧出口温度Ta,o1;当相对误差小于等于1%时,则判断空气侧流动平衡;
预设钠侧进口温度Ts,i与钠侧出口温度Ts,o,由设定传热量Q确定液钠的质量流量ms,计算得到钠侧提升力Fs与钠侧总阻力降ΔPs;计算钠侧提升力Fs与钠侧总阻力降ΔPs的相对误差当相对误差大于1%时,则返回调整钠侧进口温度Ts,i与钠侧出口温度Ts,o;当相对误差小于等于1%时,则判断钠侧流动平衡;
S3、待S2中的空气侧流动平衡且钠侧流动平衡后,计算得到空气侧传热系数ha和钠侧传热系数hs;进而再求得空气热交换器总传热系数K、对数平均温差ΔTm与总传热量Qs;再计算总传热量Qs与设定传热量Q的相对误差当相对误差大于1%时,则判断***能量不平衡,并返回步骤S2重新计算;当相对误差小于等于1%时,则判断***能量平衡,并输出空气质量流量ma、钠侧质量流量ms、空气侧进口温度Te、空气侧出口温度Ta,o1、钠侧进口温度Ts,i、钠侧出口温度Ts,o、空气热交换器总传热系数K、总传热量Qs、总传热量Qs与设定传热量Q的相对误差
2.根据权利要求1所述的钠冷堆非能动余热排出***的建模优化方法,其特征在于:步骤S2所述空气侧浮升力Fa和空气侧总阻力降ΔPa的计算步骤如下:
空气热交换器中换热的能量平衡方程如下式(1),钠侧能量损失量与空气侧能量增加量均分别等于设定传热量Q:
mscps(Ts,i-Ts,o)=macpa(Ta,o1-Te)=Q (1)
其中,Ts,i、Ts,o分别为空气热交换器的钠侧进口温度、钠侧出口温度;Te、Ta,o1分别为空气热交换器的空气侧进口温度、空气侧出口温度;ms、ma分别为空气热交换器的钠侧质量流量、空气质量流量;cps、cpa分别为空气热交换器中的钠侧比热、空气侧比热;
拔风烟囱内的空气侧浮升力Fa的计算公式如下式(2):
其中,H1为空气换热器的壳程空气出口与拔风烟囱出口高度差;ρe为空气热交换器的空气侧进口温度Te下的空气密度;ρa,om为对环境散热后拔风烟囱内平均温度Ta,om下的空气密度,式(2)中的ρa,om与式(3)中的Ta,om对应,并满足公式Ta,om由式(3)求得,其中,Ta,o1为空气热交换器的空气侧出口温度,qc为拔风烟囱与外部环境之间的热耗散,ma为空气热交换器的空气质量流量,cpa为空气侧比热;g为重力加速度;dh为高度h的微分;
空气侧总阻力降ΔPa的计算公式如下式(4):
其中,fa为翅片管外阻力系数,Na为沿空气流通方向基管的管排数,ma为空气质量流量,Sa,min为空气热交换器翅片管间最小流通面积,ρa,m为空气热交换器中的空气平均密度;fc为烟囱内沿程阻力系数,Sc为烟道流通截面积,ρa,om为对环境散热后拔风烟囱内平均温度Ta,om下的空气密度,Di为拔风烟囱烟道内直径;Nc为拔风烟囱烟道弯头数,ξc为拔风烟囱局部阻力系数,εc为拔风烟囱进出口阻力系数;
步骤S2所述的钠侧提升力Fs与钠侧总阻力降ΔPs计算步骤如下:事故热交换器液钠出口与空气热交换器液钠出口之间钠侧提升力Fs的计算公式如下式(5):
其中,H2为事故热交换器液钠出口与空气热交换器液钠出口之间的高度差;ρs,o、ρs,i分别为空气热交换器的出口液钠密度、进口液钠密度;g为重力加速度;dh为高度h的微分;
钠侧总阻力降ΔPs的计算公式如下式(6):
其中,式(6)中的fs为钠侧沿程阻力系数,钠侧沿程阻力系数fs与雷诺数Res的关系如式(7);ms为钠侧质量流量,Ss为基管内流通截面积;ρs,m为自然循环回路中液钠的平均密度;Lsur为自然循环回路中除基管外的剩余段管长;Ls为基管的总管长;di为基管内直径;Ns为钠侧弯头总数;ξs为钠侧局部阻力系数。
3.根据权利要求1所述的钠冷堆非能动余热排出***的建模优化方法,其特征在于:空气热交换器中的翅片管为三角形错列排布;所述壳程空气流动方向与管程液钠流动方向垂直;空气为竖向流通,液钠为水平方向流动;
步骤S3所述空气对流传热系数ha的计算公式如下式(8-1):
其中:λa为空气导热系数;do为基管外直径;ma为空气质量流量;Sa,min为空气热交换器翅片管间最小流通面积;μa为空气粘度;cpa为空气侧比热;Ao为单位长度翅片管总外表面积;为单位长度基管外表面积;
液钠冷却过程,于均匀热流边界条件下,液钠流动的对流传热系数hs的计算公式如下式(9):
其中:λs为液钠的导热系数;di为基管内直径;Pes为贝克来数;
根据式(8-1)和式(9),空气热交换器总传热系数K的计算公式如下式(10):
其中:β为翅化比,ha为空气侧对流传热系数;ra为空气侧污垢热阻;Ao为单位长度翅片管总外表面积;λw为管材导热系数;do为基管外直径;di基管内直径;rs为钠侧污垢热阻;η为翅片壁面总效率;hs为钠侧对流传热系数;
***总传热量Qs的计算公式如下式(11):
Qs=φKnsLsAoΔTm (11)
其中,φ为对数平均温差修正系数;K为空气热交换器总传热系数;ns为基管总管数;Ls为基管的总管长;Ao为单位长度翅片管总外表面积;ΔTm为对数平均温差。
4.根据权利要求1所述的钠冷堆非能动余热排出***的建模优化方法,其特征在于:所述翅片管为正方形直列排布;所述壳程空气流动方向与管程液钠流动方向垂直;空气为竖向流通,液钠为水平方向流动;
步骤S3所述空气对流传热系数ha的计算公式如下式(8-2):
其中:λa为空气导热系数;do为基管外直径;ma为空气质量流量;Sa,min为空气热交换器翅片管间最小流通面积;μa为空气粘度;cpa为空气侧比热;Af单位长度翅片表面积;为单位长度基管外表面积;
液钠冷却过程,于均匀热流边界条件下,液钠流动的对流传热系数hs的计算公式如下式(9):
其中:λs为液钠的导热系数;di为基管内直径;Pes为贝克来数;
根据式(8-2)和式(9),空气热交换器总传热系数K的计算公式如下式(10):
其中:β为翅化比,ha为空气侧对流传热系数;ra为空气侧污垢热阻;Ao为单位长度翅片管总外表面积;λw为管材导热系数;do为基管外直径;di基管内直径;rs为钠侧污垢热阻;η为翅片壁面总效率;hs为钠侧对流传热系数;
***总传热量Qs的计算公式如下式(11):
Qs=φKnsLsAoΔTm (11)
其中,φ为对数平均温差修正系数;K为空气热交换器总传热系数;ns为基管总管数;Ls为基管的总管长;Ao为单位长度翅片管总外表面积;ΔTm为对数平均温差。
5.根据权利要求2所述钠冷堆非能动余热排出***的建模优化方法,其特征在于:
当翅片管为三角形错列排布时,公式(4)中的翅片管外阻力系数fa为:
当翅片管为正方形直列排布时,公式(4)中的翅片管外阻力系数fa为:
其中,ma为空气质量流量,Sa,min为空气热交换器翅片管间最小流通面积;do为基管外直径;μa为空气粘度;Pt为管束横向管间距;Pl为管束纵向管间距;y为翅片管节距;Hf为翅高。
6.根据权利要求1所述钠冷堆非能动余热排出***的建模优化方法,其特征在于:构建物性参数计算模块、空气侧流动平衡计算模块、液态钠侧流动平衡计算模块与***能量平衡计算模块;所述特定物性参数存于物性参数计算模块内;步骤S2从预设空气侧出口温度至判断空气侧流动平衡的过程属于空气侧流动平衡计算模块;步骤2从预设钠侧进口温度与钠侧出口温度至判断钠侧流动平衡的过程属于液态钠侧流动平衡计算模块;步骤S3的过程属于***能量平衡计算模块;所述物性参数计算模块分别与空气侧流动平衡计算模块、液态钠侧流动平衡计算模块连接;所述空气侧流动平衡计算模块判断空气侧流动平衡后进入***能量平衡计算模块;液态钠侧流动平衡计算模块判断钠侧流动平衡后进入***能量平衡计算模块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910579525.8A CN110472270B (zh) | 2019-06-28 | 2019-06-28 | 一种钠冷堆非能动余热排出***的建模优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910579525.8A CN110472270B (zh) | 2019-06-28 | 2019-06-28 | 一种钠冷堆非能动余热排出***的建模优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110472270A true CN110472270A (zh) | 2019-11-19 |
CN110472270B CN110472270B (zh) | 2022-09-27 |
Family
ID=68507288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910579525.8A Active CN110472270B (zh) | 2019-06-28 | 2019-06-28 | 一种钠冷堆非能动余热排出***的建模优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110472270B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110991092A (zh) * | 2020-01-03 | 2020-04-10 | 中冶赛迪技术研究中心有限公司 | 一种余能一体化利用计算方法及*** |
CN112100815A (zh) * | 2020-08-14 | 2020-12-18 | 西安工程大学 | 一种用于空压机进气预处理的表冷器系列化设计方法 |
CN114893715A (zh) * | 2022-04-02 | 2022-08-12 | 安徽宇航派蒙健康科技股份有限公司 | 加热控制方法及其装置、***、计算机设备和存储介质 |
CN116204975A (zh) * | 2023-01-30 | 2023-06-02 | 哈尔滨工程大学 | 一种船用非能动余热排出***的设计优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050135544A1 (en) * | 2003-05-21 | 2005-06-23 | Jae-Hyuk Eoh | Passive safety-grade decay-heat removal method and decay-heat removal system for LMR with pool direct heat cooling process |
CN105551551A (zh) * | 2015-12-17 | 2016-05-04 | 中国原子能科学研究院 | 一种不需要常规岛投入的池式钠冷快堆低功率运行方法 |
US20180365356A1 (en) * | 2016-11-09 | 2018-12-20 | China University Of Petroleum (East China) | Design method for creep-fatigue strength of plate-fin heat exchanger |
CN109446623A (zh) * | 2018-10-19 | 2019-03-08 | 沈阳工业大学 | 一种基于传热速率平衡的固态蓄热加热特性匹配设计方法 |
-
2019
- 2019-06-28 CN CN201910579525.8A patent/CN110472270B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050135544A1 (en) * | 2003-05-21 | 2005-06-23 | Jae-Hyuk Eoh | Passive safety-grade decay-heat removal method and decay-heat removal system for LMR with pool direct heat cooling process |
CN105551551A (zh) * | 2015-12-17 | 2016-05-04 | 中国原子能科学研究院 | 一种不需要常规岛投入的池式钠冷快堆低功率运行方法 |
US20180365356A1 (en) * | 2016-11-09 | 2018-12-20 | China University Of Petroleum (East China) | Design method for creep-fatigue strength of plate-fin heat exchanger |
CN109446623A (zh) * | 2018-10-19 | 2019-03-08 | 沈阳工业大学 | 一种基于传热速率平衡的固态蓄热加热特性匹配设计方法 |
Non-Patent Citations (3)
Title |
---|
张国强等: "钠冷快堆非能动余热排出***实时仿真", 《兵器装备工程学报》 * |
李晓伟等: "模块式高温气冷堆非能动余热排出***分析与研究", 《原子能科学技术》 * |
王严冬等: "钠冷堆非能动余热排出***建模研究", 《热能动力工程》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110991092A (zh) * | 2020-01-03 | 2020-04-10 | 中冶赛迪技术研究中心有限公司 | 一种余能一体化利用计算方法及*** |
CN112100815A (zh) * | 2020-08-14 | 2020-12-18 | 西安工程大学 | 一种用于空压机进气预处理的表冷器系列化设计方法 |
CN112100815B (zh) * | 2020-08-14 | 2024-04-16 | 西安工程大学 | 一种用于空压机进气预处理的表冷器系列化设计方法 |
CN114893715A (zh) * | 2022-04-02 | 2022-08-12 | 安徽宇航派蒙健康科技股份有限公司 | 加热控制方法及其装置、***、计算机设备和存储介质 |
CN114893715B (zh) * | 2022-04-02 | 2023-11-21 | 安徽宇航派蒙健康科技股份有限公司 | 加热控制方法及其装置、***、计算机设备和存储介质 |
CN116204975A (zh) * | 2023-01-30 | 2023-06-02 | 哈尔滨工程大学 | 一种船用非能动余热排出***的设计优化方法 |
CN116204975B (zh) * | 2023-01-30 | 2023-08-15 | 哈尔滨工程大学 | 一种船用非能动余热排出***的设计优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110472270B (zh) | 2022-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110472270A (zh) | 一种钠冷堆非能动余热排出***的建模优化方法 | |
Li et al. | Performance enhancement for the natural draft dry cooling tower under crosswind condition by optimizing the water distribution | |
Xu et al. | The optimization for the straight-channel PCHE size for supercritical CO2 Brayton cycle | |
CN111707010B (zh) | 一种太阳能***均温管的设计方法 | |
Huang et al. | Performance analysis and design optimization of heat pipe sink with a variable height fin array under natural convection | |
CN107367078A (zh) | 一种利用旁路管路智能控制的太阳能蓄热*** | |
CN107166772A (zh) | 一种智能控制的太阳能蓄热*** | |
CN110220405A (zh) | 基于傅里叶数一致性的固体蓄热传热控制方法 | |
CN113689966B (zh) | 一种小型氟盐冷却高温堆非能动余排***综合实验装置 | |
Zeng et al. | Numerical study on the enhanced heat transfer characteristics of steam generator with axial economizer | |
CN112580272A (zh) | 一种基于数值模拟的lng空温式气化器的优化设计方法 | |
CN105546850B (zh) | 一种根据进口温度智能流量控制的太阳能蓄热利用*** | |
CN110186291A (zh) | 一种混流型闭式冷却塔校核计算方法 | |
CN105444440B (zh) | 一种多参数智能控制太阳能***水温的方法 | |
CN111238066B (zh) | 根据蓄热器出口温度控制泵的太阳能*** | |
CN204598554U (zh) | 一种水冷机柜及含该机柜的地下水冷*** | |
CN110398076B (zh) | 根据集热器出水温度智能控制阀门开度的太阳能*** | |
CN110285592B (zh) | 智能控制阀门开度的太阳能*** | |
CN105627597A (zh) | 一种智能计算热损失的太阳能蓄热*** | |
CN104748996B (zh) | 电厂排烟余热利用***动态响应仿真测试实验装置及方法 | |
CN112652414B (zh) | 反应堆蒸汽发生器c型管束 | |
CN209134358U (zh) | 光伏实验装置 | |
CN105466046B (zh) | 一种智能控制水位高度的太阳能集热器*** | |
CN105509347B (zh) | 一种太阳能***的智能控制方法 | |
CN209151101U (zh) | 光伏实验装置 |
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 |