旋风分离器计算
旋风分离器计算程序

备注
6A.1示例中的计算值,HOFFMAN等2001年的实验数据,使用时将此 SHEET整体COPY成一个副本,进行计算
尺寸参数
(m)
计算结果
50切割粒径
8.07
微米
D
外径
D
0.075
总效率
26%
内径
Dx
0.030
粒径
效率 进口组成
出口
放料口直径
Dd
0.060
微米
eff
%
1
0%
13.50%
Vin Theta
Co alpha Rm Vzw
Ar Vow Vocs Vom ReR Vx
53.333 m/s 0.533333333 -
0.004 0.576 0.024 m 13.588 m/s 0.029 m^2 67.86 m/s 18.53 m/s 35.46 m/s 5695.8 56.62 m/s
Din50 D50 CC0>0=L0. C0<10.1
8.000 MU 8.07 mu
0.00578 kg/kg 0.01541 kg/kg 0.00578 kg/kg
迭代计算程序 相对粗糙度 插值得桶体摩擦系数 预计总效率 升气管弗劳德数 总摩擦系数
Ks/R fair eff Frx f
内旋涡旋转速度计算值 #NAME? #NAME?
0.000
b
进料口宽度
b
0.020
2
0%
16.50%
0.000
进料口高度
a
0.038
3
1%
14.00%
0.001
4
3%
9.00%
多管旋风分离器的设计计算公式

多管旋风分离器的设计计算公式多管旋风分离器的设计计算公式是根据气体和固体颗粒的流动特性和分离原理进行推导的。
该分离器通过产生旋流在固体颗粒与气体之间产生离心力,使得固体颗粒被扔到分离器的外墙,而纯净的气体则从分离器的上部排出。
以下是多管旋风分离器的设计计算公式:1.设计分离器尺寸:-内径(D):根据气体流量和分离效果要求来确定,通常选择在100mm到2000mm之间。
-高度(H):根据气体流速和旋流的惯性力要求来确定,通常选择在2到4倍D之间。
2.分离器的旋流衰减公式:- Vc = K * (Q / A) ^ (2/3)其中,Vc是旋流速度(m/s),K是校正系数(通常在0.35到0.55之间),Q是气体流量(m^3/s),A是旋流器断面积(m^2)。
3.分离器的分离效率公式:- η = 1 - exp(-0.35 * B * (Vc / U) ^ (0.35 - 0.159 * log10(Vc / U)))其中,η是分离效率,B是分离器高度与内径的比值(H/D),U是分离器的总进气速度(m/s)。
需要注意的是,以上公式是基于经验公式和试验结果得出的,并具有一定的应用范围和适用条件。
在实际设计中,还需要考虑分离器的材质、结构和运行参数等因素,以确保设计的有效性和可靠性。
另外,关于多管旋风分离器的设计拓展,可以考虑以下方面:-分离器的材质选择:根据分离介质的性质和工况条件,选择合适的耐磨、耐腐蚀材料,如不锈钢、钛合金等。
-分离器的结构改进:优化旋流器的结构和尺寸,增加分离效率和处理能力,如采用多级分离器、多出口设计等。
-分离器的控制和优化:结合自动化控制和流体力学模拟技术,优化分离器的运行参数和分离效果,提高分离器的稳定性和可调节性。
-分离器的节能降耗:采用节能措施,如热回收和余热利用,减少分离器的能耗和环境影响。
-分离器的应用领域拓展:除了气固分离外,还可以应用于气液分离、液固分离等领域,如石油化工、环保工程等。
旋风分离器计算程序--Muschelknautz模型方法

0.00543 kg/kg 0.00543 kg/kg 0.00523 kg/kg
迭代计算程序 ReR初值 相对粗糙度 查表得桶体摩擦系统 预计总效率 升气管弗劳德数 总摩擦系数 内旋涡旋转速度计算值
ReR0 Ks/R
fair eff Frx f
V0cs
压降计算 分离器中损失 旋转涡核与升气管损失 加速压力损失 总压力损失
Vin Theta
Co alpha Rm Vzw
Ar Vow Vocs Vom ReR Vx
Din50 D50 CC0>0=L0. C0<10.1
24.482 m/s 0.3 -
0.160 0.741 0.150 m 3.105 m/s 2.506 m^2 14.000 m/s 18.000 m/s 15.875 m/s 3570.324 17.038 m/s
Rho_p
1800 kg/m^3
90 100%
0.37%
0.004
堆积密度度
Rho_b
900 kg/m^3
100 100%
0.09%
0.001
>100
0
-
1.00
0.92
Muschelknautz Modeld Rin
入口速度 T因子 颗粒/气体 (m/m) 入口收缩系数 几何平均半径 器壁表面轴向速度 摩擦阻力总面积 器壁表面切向速度 内旋涡旋转速度 气体平均旋转速度 旋风分离器ReR 升气管中气流速度 分离效率 进口中位径 切割粒径 极限浓度
40 100%
14.65%
0.147
操作参数
50 100%
8.24%
0.082
旋风分离器计算结果

旋风分离器计算结果标准化工作室编码[XX968T-XX89628-XJ668-XT689N]旋风除尘器性能的模拟计算一、下图为旋风除尘器几何形状及尺寸,如图1所示,图中D、L及入口截面的长宽比在数值模拟中将进行变化与调整,其余参数保持不变。
图1 旋风分离器几何形状及尺寸(正视图)旋风分离器的空间视图如图2所示。
图2 旋风分离器空间视图二、旋风分离器数值仿真中的网格划分仿真计算时,首先对旋风除尘器进行网格划分处理,计算网格采用非结构化正交网格,如图3所示。
图3 数值仿真时旋风分离器的网格划分(空间)图4为从空间不同角度所观测到的旋风分离器空间网格。
图4 旋风分离器空间网格空间视图本数值仿真生成的非结构化空间网格数大约为125万,当几何尺寸(如D、L及长宽比)改变时,网格数会略有变化。
三、对旋风分离器的数值模拟仿真采用混合模型,应用Eulerian(欧拉)模型,欧拉方法,对每种工况条件下进行旋风分离器流场与浓度场的计算,计算残差<10-5,每种工况迭代约50000步,采用惠普工作站计算,CPU耗时约12h。
以下是计算结果的后处理显示结果。
由于计算算例较多,此处仅列出了两种工况条件下的计算后处理结果。
图5是L=1.3m,D=1.05m 入口长宽比1:3,入口速度10m/s时,在y=0截面(旋风分离器中心截面)上粒径为88微米烟尘的体积百分数含量分布图。
可以明显看出由于旋风除尘器的离心作用,灰尘被甩到外壁附近,而在靠近中心排烟筒下方筒壁四周,烟尘的体积浓度最大。
粒径88微米烟尘的空间浓度分布(空间)粒径88微米烟尘的浓度分布(旋风分离器中心截面)粒径200微米烟尘的空间浓度分布(空间)粒径200微米烟尘的浓度分布(旋风分离器中心截面)图5 L=1.3m、D=1.05m、长宽比1:3,入口速度10m/s时烟尘空间分布粒径88微米烟尘的空间浓度分布(空间)粒径88微米烟尘的浓度分布(旋风分离器中心截面)粒径200微米烟尘的空间浓度分布(空间)粒径200微米烟尘的浓度分布(旋风分离器中心截面)图6 L=2.3m、D=1.5m、长宽比1:1,入口速度15m/s时烟尘空间分布四、计算结果计算中,首先确定几何尺寸L,按照给定的两种烟尘颗粒,分别对L=2.3m、L=1.8m、L=1.3m、L=0.8m四种情况进行对比计算,对比计算结果为L=2.3m、L=1.3m时除尘效率较高。
旋风分离器计算公式

旋风分离器计算公式
汽量:kg/h178047水量:kg/h1148636汽密度:kg/m316.52水密度:kg/m3814旋分数量:个60旋分直径:mm290锅炉给水量kg/h1351336补充水kg/h165200入口引入管直径mm100
过程量:
旋分的蒸汽负荷t/h 2.96745判断旋分个数是否合适
每个旋风汽量:m3/h179.628每个旋风水量:m3/h23.5183循环倍率1.17647入口引入管的折算汽速度Wo'm/s 6.35626入口引入管的折算水速度Wo"m/s0.02276总速度w1m/s 6.37902判断汽水混合物的速度
循环速度W01m/s0.15176
结果:
旋分的阻力系数4阻力损失:Pa160.824最大阻力损失:Pa13932.1
布袋的截面积x过滤风速x60=每小时的处理风量(m3/h)布袋除尘器的理风量是指除尘设备在单位时间内所能净化气体的体积量。
布袋除尘器的理风量,是袋式除尘器设计中最重要的因素之一。
设计原理:根据风量设计袋式除尘器时,一般不能使除尘器在超过规定风量的情况下运行,否则,除尘器滤袋容易堵塞,寿命缩短,压力损失大幅度上升,除尘效率和除尘效果也要大大降低。
但也不能将粉尘处理的风量选的过大,否则增加除尘器除尘设备投资和占地面
积,而且浪费资源,不节能。
合理的选择处理风量常常是根据工艺情况和经验来决定的。
由此可以看出不是单纯的公式或者是有固定数值就可以算出风时,得多方面考虑。
旋风分离器的工艺计算

旋风分离器的工艺计算》:*目录一.前言 ............................................................................................................. 错误!未定义书签。
应用范围及特点....................................................................................... 错误!未定义书签。
分离原理................................................................................................... 错误!未定义书签。
分离方法................................................................................................... 错误!未定义书签。
)性能指标 ................................................................................................. 错误!未定义书签。
二.旋风分离器的工艺计算.............................................................................. 错误!未定义书签。
旋风分离器直径的计算........................................................................... 错误!未定义书签。
由已知求出的直径做验算....................................................................... 错误!未定义书签。
旋风分离器的设计计算

旋风分离器的设计计算
R.K.Sinnott
【期刊名称】《医药工程设计》
【年(卷),期】1989(000)004
【摘要】旋风分离器是气固分离器的主要型式,靠离心力进行分离,应用广泛。
旋风分离器结构简单,可用各种材料制作,能用于高温高压。
旋风分离器适用于分离粒径为5μm 以上的颗粒,对于粒径在0.5μm 以下的较小微粒,在旋风分离器中会发生附聚作用。
旋风分离器最通用的结构型式是反流
【总页数】4页(P11-14)
【作者】R.K.Sinnott
【作者单位】
【正文语种】中文
【中图分类】R
【相关文献】
1.旋风分离器旋风长度的分析计算 [J], 高翠芝;孙国刚;董瑞倩
2.一种新型旋风分离器--均匀料流旋风分离器 [J], 谭岳云
3.旋风分离器简捷设计计算方法 [J], 刘忠文
4.旋风分离器自然旋风长的影响因素 [J], 高助威; 王娟; 王江云; 毛羽
5.旋风分离器设计计算的研究 [J], 蔡安江
因版权原因,仅展示原文概要,查看原文内容请购买。
旋风分离器计算

作成作成::时间时间::2009.5.14一、問題提出PHLIPS FC9262/01這款吸塵器不是旋風除塵式的,現在要用這款吸塵器測參數選擇旋風分離裝置。
二、計算過程1.選擇工作狀況選擇工作狀況::根據空氣曲線選擇吸入效率最高點的真空度和流量作為旋風分離器的工作狀態。
吸塵器旋風分離器選擇Bryan_Wang已知最大真空度h和最大流量Q,則H-Q曲線的兩個軸截距已知,可確H-Q直線的方程。
再在這個直線上求得吸入功率H*Q最高點(求導數得)。
求解過程不再詳述。
求得最大吸入功率時真空度H=16.5kPa;流量Q=18.5L/s;吸入功率P2=305.25w發現計算得到的吸入功率最大值與產品標稱值375W相差一些,可能是由于測量誤差存在以及壓力損失的原因。
現將真空度及流量按照吸入功率計算值與實際值的比例放大,得真空度H=18.3kPa;流量Q=20.5L/s;2.選擇旋風分離器為使旋風分離裝置體積最小,選擇允許的最小旋風分離器尺寸。
一般旋風分離器筒體直徑不小于50mm,故選擇筒體直徑為50mm。
按照標準旋風分離器的尺寸比例,確定旋風除塵器的結構尺寸。
D0=50mmb=12.5mma=25mmde=25mmh0=20mmh=75mmH-h=100mmD2=12.5mm計算α約為11度一般要求旋風分離器進氣速度不超過25m/s,這里取旋風分離器進氣速度為22m/s.計算入口面積為S=3.125e-4平方米。
則單個旋風除塵器流量為Q=6.9e-3平方米/秒則所需旋風除塵器個數為3個計算分級效率根據GB/T 20291-2006吸塵器標準,這里使用標準礦物灰塵,為大理石沙。
进气粒径分布1030581001903757501500201010102016113顆粒密度ρp=2700kg/m3進口含塵濃度取為10g/Nm3,大致選取空氣粘度μ=1.8e-6Pa*s按照以下公式計算顆粒分級效率:平均粒徑(μm)比重(%)計算結果為d(m)1E-053E-056E-051E-042E-044E-048E-040.0023E-071E-075E-08ηi 111111110.91140.6750.5校核分割粒徑校核分割粒徑x x 5050::按照以下公式計算:計算得知在所有平均粒徑計算得到的分級效率都為100%,而分級效率為50%的粒徑為0.05微米。
旋风分离器的工艺计算

旋风分离器的工艺计算目录一.前言 (3)1.1应用范围及特点 (3)1.2分离原理 (3)1.3分离方法 (4)1.4性能指标 (4)二.旋风分离器的工艺计算 (4)2.1旋风分离器直径的计算 (5)2.2由已知求出的直径做验算 (5)2.2.1计算气体流速 (5)2.2.2计算旋风分离器的压力损失 (5)2.2.3旋风分离器的工作范围 (6)2.3进出气管径计算 (6)三.旋风分离器的性能参数 (6)3.1分离性能 (6)3.1.1临界粒径d pc (7)3.1.2分离效率 (8)3.2旋风分离器的压强降 (8)四.旋风分离器的形状设计 (9)五.入口管道设计 (10)六.尘粒排出设计 (10)七.算例(以天然气作为需要分离气体) (11)7.1工作原理 (11)7.2基本计算公式 (12)7.3算例 (13)八.影响旋风分离器效率的因素 (15)8.1气体进口速度 (15)8.2气液密度差 (15)8.3旋转半径 (15)参考文献 (15)旋风分离器的工艺计算摘要:分离器已经使用十分广泛无论在家庭生活中还是工业生产,而且种类繁多每种都有各自的优缺点。
现阶段旋风分离器运用比较广泛,它的性能的好坏主要决定于旋风分离器性能的强弱。
这篇文章主要是讨论旋风分离器工艺计算。
旋风分离器是利用离心力作用净制气体,主要功能是尽可能除去输送介质气体中携带的固体颗粒杂质和液滴,以达到气固液分离,以保证管道及设备的正常运行。
在本篇文章中,主要是对旋风分离器进行工艺计算。
关键字:旋风分离器、工艺计算一.前言旋风分离器设备的主要功能是尽可能除去输送介质气体中携带的固体颗粒杂质和液滴,达到气固液分离,以保证管道及设备的正常运行。
它是利用旋转气流产生的离心力将尘粒从气流中分离出来。
旋风分离器结构简单,没有转动部分制造方便、分离效率高,并可用于高温含尘气体的分离,而得到广泛运用。
旋风分离器采用立式圆筒结构,内部沿轴向分为集液区、旋风分离区、净化室区等。
旋风分离器计算计算

数量:6
3×160m^3补糖罐、1×120m^3补料罐 旋风分离器工艺设计
大罐排气空气流量 旋风进口设计流速
高度系数 宽度系数
D H b L1 L2 L0 D1 D2
2.5 m^3/S 18 m/s
0.5 圆整 0.22 1.26263 1.2 m
0.6 600mm 0.264 2600m
1.92 1900mm 3 3000mm
120m^3发酵大罐阻沫器参数
大罐排气空气流量 1.4 旋风进口设计流/S m/s
圆整
进口管道直径 0.273 0.273
L1
1.2 1200mm
L2
0.8 800mm
L0
0.6 600mm
D1
0.25 DN250
D2
0.125 DN125
0.66 600mm 0.6 DN600 0.3 DN300
数量:2
3×20m^3补油罐、3×20m^3补醇罐、1 ×20m^3补小料罐旋风分离器工艺设计
大罐排气空气流量 旋风进口设计流速
高度系数 宽度系数
D H b L1 L2 L0 D1 D2
0.83333 20 0.5 0.22 0.37879 0.31 0.1364 0.992 1.55 0.341 0.31 0.155
m^3/S m/s
圆整
0.62 m 310mm 140mm 1000mm 1600mm 340mm DN300 DN150
数量:2
二、阻 沫器工 艺计算
500m^3发酵大罐阻沫器工艺设计
L0
D1
L1
D
L
大罐排气空气流量 旋风进口设计流速
高度系数 宽度系数
D
进口管道直径
分离器设计计算

分离器设计题:气量110万方/日,压力为8MPa (表),温度270C ,相对密度为0.6,要求尽可能分离掉大量的固体颗粒(300um )和少量液滴,经查图确定颗粒的沉降速度为0.38m/s 。
试分别设计立式、卧式及旋风分离器。
解:已知:Q g =110万m 3/d P=8MPa T=27℃ S=0.6 ω=0.38m /s 由s=0.6<0.7可得,P pc =4.778-0.248s=4.778-0.248×0.6=4.629MPaT pc =106.1+152.21×0.6=197.4KP pr =P/ P pc =8/4.629=1.7T pr =T/ T pc =300/197.4=1.5根据P pr 和T pr 的值查表2-1可得Z=0.83取η =0.77则ν=ηω=0.77×0.38=0.29m/s实用流量:Q 1=86400g Q ×P 101325.0×293TZ =864001100000×8101325.0×293300×0.83=0.137m 3/s1、立式分离器筒体直径:D=(νπ14Q )1/2=0.78m高度:H=3.3D=3.3×0.78=2.57m进口直径:D 1=(15785.0137.0⨯)1/2=0.11m出口直径:D 2=(10785.0137.0⨯)1/2=0.13m 2、卧式分离器取A=3.3,77.0=η筒体直径:D=(AQ πηω14)1/2=0.43m 长度L=3.3×0.43=1.42m进口直径:D 1=(15785.0137.0⨯)1/2=0.11m 出口直径:D 2=(10785.0137.0⨯)1/2=0.13m 3、旋风分离器取k=1,筒体直径:D=3.39×10-5×(P TZQ g )1/2×K=3.39×10-5×(8110000083.0300⨯⨯)1/2 ×1=0.20m 进口直径:D 1=0.47D=0.47×0.20=0.09m 出口直径:D 2=0.67D=0.67×0.20=0.13m 进口流速:V 1=5.096×10-9×21g PD TZQ =5.096×10-9×2)09.0(8110000083.0300⨯⨯⨯=21.5m/s 出口流速:V 2=5.096×10-9×22g PD TZQ =5.096×10-9×2)13.0(8110000083.0300⨯⨯⨯=10.3m/s筒体流速:V=5.096×10-9×2g PD TZQ =5.096×10-9×2)20.0(8110000083.0300⨯⨯⨯=4.36m/s 因为15m/s<V 1<25m/s,5m/s<V 2<15m/s,2.45m/s<V<4.43m/s,所以K 的取值符合要求,就按照上面计算出来的D ,D 1,D 2设计旋风分离器。
旋风分离器的工艺计算

旋风分离器的工艺计算目录一.前言 (3)1.1应用范围及特点 (3)1.2分离原理 (3)1.3分离方法 (4)1.4性能指标 (4)二.旋风分离器的工艺计算 (4)2.1旋风分离器直径的计算 (5)2.2由已知求出的直径做验算 (5)2.2.1计算气体流速 (5)2.2.2计算旋风分离器的压力损失 (5)2.2.3旋风分离器的工作范围 (6)2.3进出气管径计算 (6)三.旋风分离器的性能参数 (6)3.1分离性能 (6)3.1.1临界粒径d pc (7)3.1.2分离效率 (8)3.2旋风分离器的压强降 (8)四.旋风分离器的形状设计 (9)五.入口管道设计 (10)六.尘粒排出设计 (10)七.算例(以天然气作为需要分离气体) (11)7.1工作原理 (11)7.2基本计算公式 (12)7.3算例 (13)八.影响旋风分离器效率的因素 (15)8.1气体进口速度 (15)8.2气液密度差 (15)8.3旋转半径 (15)参考文献 (15)旋风分离器的工艺计算摘要:分离器已经使用十分广泛无论在家庭生活中还是工业生产,而且种类繁多每种都有各自的优缺点。
现阶段旋风分离器运用比较广泛,它的性能的好坏主要决定于旋风分离器性能的强弱。
这篇文章主要是讨论旋风分离器工艺计算。
旋风分离器是利用离心力作用净制气体,主要功能是尽可能除去输送介质气体中携带的固体颗粒杂质和液滴,以达到气固液分离,以保证管道及设备的正常运行。
在本篇文章中,主要是对旋风分离器进行工艺计算。
关键字:旋风分离器、工艺计算一.前言旋风分离器设备的主要功能是尽可能除去输送介质气体中携带的固体颗粒杂质和液滴,达到气固液分离,以保证管道及设备的正常运行。
它是利用旋转气流产生的离心力将尘粒从气流中分离出来。
旋风分离器结构简单,没有转动部分制造方便、分离效率高,并可用于高温含尘气体的分离,而得到广泛运用。
旋风分离器采用立式圆筒结构,内部沿轴向分为集液区、旋风分离区、净化室区等。
旋风分离器计算模型

⎤ )⎥ ⎥ ⎥ ⎦
器壁表面轴向速度 0.9 ∙ 245000 = 5.57 m/s 5 3600 ∙ π ∙ 4 − 1.658 出风管(内筒)弗劳德数: 245000 F = = 3.86 π ∙ 2.2 2.2 3600 ∙ 4 ∙ 2 ∙ 9.8 ∙ 2 V = 雷诺准数初始值: Re 5 0.60 ∙ 5.57 ∙ 2 1.658 = ∙ = 30797 2.865 ∙ 10 15.7
—— —— —— —— —— ——
旋风筒总高度,m 旋风筒锥体高度,m 出料口直径,m 旋风筒直径,m 内筒直径,m 旋风筒出料口至内筒底部距离 m
1.2、 入口收缩系数 Kin
K 1 = T ⎡ T T ⎢ ∙ ⎢1 − 1 + 4 ∙ − 4 2 ⎢ ⎣ ∙ 1− 1−T ∙ 2T − T 1+C ⎤ ⎥ ⎥ ⎥ ⎦
固气比: C o =
C i 0 .5 = = 0.833 ρ a 0 .6
245000 = 17.62 3600 × 1.5 × 2.575 m/s
旋风筒进口风速: Vin = 入口收缩系数:
⎡ (1 − 0.6 ) ∙ (2 ∙ 0.6 − 0.6 1 ⎢ 0.6 0.6 K = ∙ ⎢1 − 1 + 4 ∙ − ∙ 1− 0.6 ⎢ 4 2 1 + 0.833 ⎣ = 0.763 器壁表面切向速度: 2 × 17.62 × 2.8 = 25.86 m/s V = 0.763 × 5 旋风筒筒体与内筒之间的平均半径:
器壁相对粗糙度: K SR = 查表得 fa=0.0058。 总摩擦系数:
KS D 1000 ⋅ 2
旋风分离器的选型计算

图 2 标 准 旋 风 分 离器 结构 图
摘 要 :旋风分离器是 指采用离心沉降的原理对气流 中所含 固体颗粒或液滴进行分离的设备 ,工业中一般 采用标准 旋风分离器的结构 比例参数进行旋 风分离器的设计 ,通过实 际项 目中选 用的旋风分离器的具体数据 , 总结 出进行旋风分离器选 型计算所需要 的参数 、所 采用 的具体计算方 法和步骤 。 关 键 词 :离心沉 降 ;旋 风分离器 ;人 口气 速 ;允许 压降 ;阻力系数 中 图分 类 号 :TE 624.4 1 文献 标 识 码 : A 文 章 编 号 : 1671—0460 (2018)02—0415—03
1 基础资料
1.1 工 艺流 程 有关 此旋 风 分离 器 的工艺 系统 流程 图如 图 1。
图 1 旋 风 分 离 器 工 艺 流 程 图
Fig.1 Process draw ing of cyclone separator
在流 程 图 中 ,D一002即 为新 增 的旋 风分 离 器 。 1.2 相关 数据
M P
p 而
, 、
(2)
式 中 : 气体 密度 ,kg/m ; _一。空气 摩 -7质 量 ,g/mol;
P一 压 强 , kPa: 一 理 想气体 常 数 ,8.3 14J/(g·K】。
则 有 p p =面 面MP =
=0. 33~ 0.83 kg/ m
根据 (1)有 :
Type Selection Calculation of Cyclone Separator
LIU Hng (CNPC Northeast Refining& Chem ical Engineering Co.,Ltd.,Liaoning Dalian 116085,China)