上部バナーに（↑）広告している方へ：著作権に抵触しています。本サイト（http://homepage3.nifty.com/HTU/binran/index/）に無断で広告することを禁じます。

充填物ミラックス　恩田モデルソフトホームへ

     Web　PACKING HANDBOOK
             ウエブ充填物便覧
NTU

NTUgraph


　　　メールによるお問い合わせ先
　
　　　　　　　　　　　　　　　　　　　　　　　恩田モデルソフトの使い方

HTU 計算ソフト、充填物、物質移動、恩田モデル、へンリー定数、物性　
ミラックス327(at＝125ｍ2/ｍ3）
HTU Calculation soft, Mass-Transfer, Onda-model, Henry-constant

                                                                    DEK   下井洋一　Yoichi Shimoi 　 2007/9/1                       
               DEK（ミラックス:Mirax）Home       物性推算（大江修造）     NIST    
  post    連絡先　メール

目次　（Index)   §１ Onda model　＆　HTU              1-5)NTUOL            1-6)Pressure drop         
                          2-1)Property of water & Air       2-2)Diffusion coefficient in Gas        2-3)Diffusion coefficient in Water     
                    §3 Henry's constant             §§ Reference


蒸留には多くの計算モデルが提案されていますが、ガス吸収、ガス放散などの充填物を用いた定温物質移動操作には恩田モデルが適しています。
本HTU計算ソフトは目次のような内容で構成されていますが、HTUに関しては恩田モデルを用いて、aw（濡れ面積）→ｋL→ｋGの順に計算し、
最終的にＨＴＵOG、HTUOLを求められます。
計算例題としてマグネシウム水溶液による脱硫を赤字数値で、またベンゼンの放散を緑字数値で記載しましたので、
これらを順次空白に代入して、最後にクリックすれば、夫々のＨＴＵを計算できます。（数値をDrag & Pasteして貼り付けてください）
awに関しては初めに適当な数値 (like 60 m2/m3)を入れてクリックし、正しいawが出たら、これを改めて代入、もう一度クリックします。
有機化合物の放散は対象が広いが、本ソフトでは８０種類の有害な有機化合物放散のHTUOLを計算できます。
粘度、表面張力、液側拡散係数、ガス側拡散係数、へンリー定数などの物性およびＮＴＵＯＬや圧損の計算ソフトを掲載したので、広範な応用が可能です。

For distillation many good models have been already proposed. But Onda model is suitable for the iso-thermal operation like gas absorption or stripping.
This HTU calculation software has the contents as shown in the index.
HTU is calculated using Onda-model according to the procedure; aw(wetted area) → KL → KG and finally HTUOG or HTUOL can be gotten.
As the calculation example De-sulfurization is loaded in the red figures and stripping of organic matter in the green figure.
So you can get each HTU filling the blanks with figures in turn and clicking.(It is convenient to drag the figures and paste in the blank)
First input contemporal "aw" (like 60 m2/m3). Next click again after inputting the correct "aw" gotten by first click.
This soft is corresponded to 80 kinds of organic matters stripping.
Viscosity, Surface tension, Diffusion coefficient in liquid phase, Diffusion coefficient in gas phase, Henry constant and other properties are also loaded.
So you can calculate the broad absorption and stripping problems.

   §１恩田の式、Onda model　＆HTU & NTU & ΔP   (Top)

           １-１）充填物濡れ面積の計算/ To get wetted area/aw

PROJECT NAME, Job Number,Date and so on          (For your MEMO)

σC (Critical surface tention of packing material)  A1= 　 (Example: 386208 Kg/hr2 for PP)
σ (Surface tention of liquid (Water))   B2= (Example: 933120  761983Kg/hr2 for water See 2-1) )
L (Liquid mass velosity)   C3= (Example: 40000  100000Kg/m2/hr)
at (Geometric surface area of packing)　　D4= 　(Example: 125 m2/m3/hr for Mirax327)
μL (Viscosity of liquid)   E5=　 (Example: 3.6  2.39Kg/m/hr for water25℃  100℃See 2-1))
ρL (Density of liquid)　F6= 　 (Example: 997 Kg/m3 　958Kg/m3 or water; reference：997kg/m3 25℃： 958kg/m3 100℃)
g (gravity coefficient)  G7= 　 (Example: 127010000 m/hr2 on the Earth)
                 

        1-２）液側物質移動係数の計算/ To get Liquid phase mass-transfer coefficient/ kL (m/hr) &KLa(1/hr)

aw (Wetted area of packing)  H8= 　 (First input 60 m2/m3 approximately, then input again correct "aw" gotten by click)
DL (Diffusion coefficient of transfered gas in liquid (Water))   I9= (Example: 0.000005364   m2/hr for SO2 in water 25℃  0.0000109 Benzene in water 100℃See 2-3))
atDpL (Constant for Packing)   J10= (Example: 5.8 [-] for Mirax 327)


       １-３）ガス側物質移動係数の計算/ To get kG (kmol/m2/hr/atm) & KGa(kmol/m3/hr/atm)

R (Gas constant)  K11= 　 (Example: input 0.08205 ｍ3*atm/kmol/K)
T (Tempetarute)   L12= (Example: 298.2   373.2°K)
DG (Diffusion coefficient of transfered gas in inert gas (air))   M13= (Example: 0.0432 ｍ2/ｈｒ 25℃ ;  in case of Benzene gas phase is negligible 0.0432See 2-2) )
atDpG (Constant for Packing)   N14= (Example: 3.4 [-] for Mirax 327)
G (Gas mass velosity)   O15= (Example: 10000 100Kg/m2/hr  "G" should be decided for λ to be 2 to 3 See 1-5))
μG(Viscosity of gas    P16= (Example: 0.0648 Kg/m/hr for air 25℃ ;  in case of Benzene gas phase is negligible 0.0648See 2-1)
ρG(Gas density) Q17= (Example: 1.2 Kg/m3 for air 25℃, 0.5977 Kg/m3 for steam 100℃)


         １-４）充填物のＨＴＵ計算/ To get ＨＴＵ(Height of Transfer Unit)(m)

β(Enhancement factor)  R18= 　 (Example: β=6 for Desulfurization by (Mg(OH)2 aqueouse solution in case of Benzene β=1)
                  β=1 for physical absorption/desorption in water, β=∞=1000000 for instanteneous chemical reaction like SO2-NaOH aq. sol.)
H (Henry's constant   S19= (Example: 40 (atm/molfr) for SO2 in water 25℃;  4530 (atm/molfr) for Benzen in wate 100℃See 3 )
π (Pressure)   T20= (Example: 1 atm　1 atm)
Gas(Vapor)Molecular weight MWG   U21= (Example: 29 　　(Air);　18 (Steam)

                                                Primary Calculation   
                                                

         １-５）NTUOL計算/To get NTUOL  (Top)
This calculation is for gas stripping. (Waste water treatment, for example)
In the case of steam stripping, the assumption is that the liquid is suplied at boiling point
and heat loss is neglected.(In gas absorption NTU can be calculated as LN(y1/y2) in many case.
Because gas equilibrium pressure is Zero for chemical absortion)

Water rate  L A1=　 50000(Kg/hr)
Inlet concentration  x top B2=　 (ppm/ppb)(Inlet Methyl Mercaptan=1000ppb)
Outlet concentration  x bottom   C3= (ppm/ppb)(Outlet Methyl Mercaptan=50ppb)
Henry constant H  D4= (atm/molfr) (H=200 atm/molfr for Methyl Mercaptan in waterSee 3)
Total pressure π E5= (atm)(1 atm)
Lambdaλ(=m/(LM/GM)=(H/π)/(LM/GM) F6= (-)λ> 1 for stripping (λ=2.16) (λ=0.955 at Minimum required gas)
                                                                                      
         １-６）充填物圧損計算/To get Pressure Drop (Top)
L (Liquid mass velosity)  A1= 　 (Example: 40000  100000Kg/m2/hr)
G (Gas mass velosity)   B2= (Example: 10000 100Kg/m2/hr  "G" should be decided
for λ to be 2 to 3)
ρG(Gas density)  C3= (Example: 1.2 Kg/m3 for air 25℃, 0.5977 Kg/m3 for steam 100℃)
                                                                                      



   §２　物性、Matter property 

         ２-１）水、水蒸気、空気の物性　 (The property of Water & Water Vapor & Air;
        Viscosity, Surface tention, etc)   (Top)

                    Temperature   A1= (℃(centigrade or Celsius)    



  Preliminary Calculation         
                        

      
  ２-２）ガス側拡散係数；　藤田の式　DG (Gas phase diffusion coefficient, Fujita's Equation) (Top)
Inert gas(1) 　

Transfered gas(2) 　


        Inert gas(1)                               

                             Molecular　Weight(MW1)  A1=     
                            Critical Temperature(Tc1)  B2= (℃)         
                            Critical Pressure(Pc1)   C3=  (atm)  


       Transfered gas(2) 
                           Molecular Weight(MW2)  D4=      
                           Critical Temperature(Tc2)  E5=    (℃)      
                           Critical Pressure(Pc2)  F6=   (atm) 


                          Temperature   G7= (℃(centigrade or Celsius)    


                      

 
  ２-３）液側拡散係数　　DL (Liquid phase diffusioncoefficient)   (Top)
        水中の拡散係数の概算式(Wilkeの式）（１は溶質(solute)、２は溶媒(solvent)）

           DL=7.4*10^(-8)*(x*M2)^0.5*(T)/(μ2*V1^0.6)     (cm2/sec)  

                　Ｍ２＝水の分子量(MW)＝18、　Ｔ＝温度(Temp)（Ｋ）、　　μ２〔水の粘度(Vis. of water)　cp）は温度の関数としてプログラムに内臓。
              　 V1（沸点分子容：　ml/mol）は分子量と比重から計算する。（例えばトルエンの比重は常温で0.87、沸点（111℃）で0.78で体積は10％膨張する。
                 有機化合物の沸点分子容はデータがないので、分子量と常温の比重から常温分子容を計算し、1.2倍して沸点分子容とする）
                 pulldownメニューに比重と分子量を掲載する。これを枠に代入すればプログラム内部でV1を計算する。 
                  （比重のデータが見つからない時は、　0.75　を代入）
                 水の会合度（ｘ）の温度データが見当たらないので、20℃の値（ｘ＝2.6）を温度に関係なく用いた。
                                                   

         ２-３-１）(有害)有機物質 ((Toxic) Organic　matter)の比重と分子量

Organic Material and Coefficient.　

 If you can not find data　（MW & Specific gravity) in pull-down menu, please see the hand book like "Perry's Chemical Handbook" or 「化学便覧」                         

   MW  A1=                   (From the pulldown menu) 
  Specific Gravity   B2=      (From the pulldown menu)

  Temperature   C3= (℃(centigrade or Celsius)


                          

 


    §３へンリー定数/ Henry's Constant(Top)

The Coefficients are by the latest literature from the web site of NIST.
If you can not find the name of your serching matte,r please serch NIST
The unit of Henry's Constant is [atm/molfr]
(The number of the Material will be newly added on)


         ３-１）(有害)有機物質 ((Toxic) Organic matter)


Organic Material and Coefficient.　

         ３-２）無機ガス (Inorganic Gas)


Various Gases and Coefficient.　

         ３-３）計算 (Calculation Click!)

Coefficient 1( k°H )  A1=                   (From the pulldown menu) 
Coefficient 2(d(ln(kH))/d(1/T))   B2=      (From the pulldown menu)

Temperature   C3= (℃(centigrade or Celsius)


                          

§§　Reference  (Top)


DATA　1(Onda model aw/kL)

Onda model aw/kL

DATA　２ (Onda model kG)
Onda model kG
DATA　３ (HTU)
HTU
DATA　４ (DG Fujita Equation)
DG Fujita Equation

DATA　５ 充填物　Mirax327
Mirax327

あなたはcounter番目の訪問者です  (Top)