Muroran-IT Academic Resources Archive Title Author(s) Citation Issue Date URL 円管内における粒体の水力輸送(第7報) : 流動様式と付加圧力損失内容の変化について奥田, 教海; 山岸, 英明室蘭工業大学研究報告.理工編 Vol.10 No.3, pp.419-425, 1981 1981-11-30 http://hdl.handle.net/10258/3727 Rights Type Journal Article See also Muroran-IT Academic Resources Archive Copyright Policy Muroran Institute of Technology 円管内における粒体の水力輸送第 7報流動様式と付加圧力損失内容の変化について奥田教海・山岸英明 HydraulicTransportofSolidsinPipes Part7 .FlowPatternsofSphere-WaterMixtureFlowsandVarious FactorsofAdditionalPressureLossesaboveThoseof ClearWateronly KyokaiOkudaandHideakiYamagishi A b s t r a c t Thea n a l y s i si sp r e s e n t e d，w h i c hi sonv a r i a t i o no ft h ef a c t o r sf o ra d d i t i o n a lp r e s s u r e1 0 鈴 e si ns p h e r e w a t e rm i x t u r ef l o wi np i p e sa b o v et h o s eo fc 1e a rw a t e ro n l y . Thef a c t o r sa r ea sf o l l o w sfromt h ea n a l y s i s: ( a )R o l l i n gF r i c t i o no ft h es p h e r e sont h ep i p ew a l l ( b ) D r agV a r i a t i o no ft h es p h e r e ss u b j e c tt of l u c t u a t i n gp i p ef l o w ( c ) MeanDr agd u et ov e l o c i t yd i f f e r e n c e sb e t w e e nw a t e ra n ds p h e r e si nt h ep i p e t h er e s u l t so ft h ea n a l y s i sa r e R e g a r d i n gac e r t a i na s s u m p t i o n sont h ed a t ao ft h et h r e ef l o wp a t t e r n s， o b t a i n e da sf o l l o w s: ( a ) Ther o l l i n gf r i c t i o ni sl a r g e s ti nmovingb e df l o wp a t t e r n s ， andl a r g e ri nh e t e r o g e n e o u sf l o wp a t t e r n s r 1yhomQgeneousf l o wp a t t e r n s . t h a ni nn e a ( b ) T hemoret h ec o n c e n t r a t i o no fs p h e r e si n c r e a s e s， t h eg r e a t e r出 ee f f e c to fr o l l i n gf r i c t i o na r i s e si ne v e r y f l o wp a t t e r n s . 1.まえカずき問題名の第 6報 1)で報告した実験の結果を更に検討し，観測された流動様式ごとの諸値を輸送管内の全圧力損失に関連させて考え，水のみの場合の圧力損失に上積みされる付加圧力損失が，粒子の転動摩擦，粒子の衝突による速度変動に基づく抗力，流体と粒子の平均速度差から求められる抗力などから構成きれること，およびそれらが全圧カ損失を占有する割合が流動様式ごとにどの様に変化するかということについて，実験値にある推定を加えて解析したので，それらについて報告する。この解析の基礎には，第 5報 2) で報告した粒子のランダム運動過程を考慮に入れている。上述の観測された諸値というのは，管内濃度，管内濃度比分布，水と粒子の速 ( 3 1 ) 4 2 0 奥田教海・山岸英明度比，粒子と管の直径比，粒子ブルード数，粒子レイノルズ数などを指す。 I I . 実験装置と実験方法詳細は第 6報を参照されたい。図 -1には実験装置全体の平面図を示す。①④は観測点で，各点では図 -2に示す管断面内の水平な層を粒子が通過する頻度を求め，その値より各層の濃度比分布を計算した。 ④ 四ーーー・・ ‘ ーー四ー 4 図 -1 実験装置 - -5層一ー 4層一ー 3層・一回 2層 1 ーー 1層図 -2 断面内側定位置 I I I . 全力損失内容の解析水力輸送の全圧力損失は，次式に示すような諸項から構成されるものと考える。 (1.全圧力損失項)= ( 2 . 粒子転動摩擦項)+ ( 3 . 粒子速度変動抗力項) + (4. 全粒子抗力項)+ (5. 水のみの圧力損失項) (1) 管内濃度が高い場合，管下部における粒子が相互に摩擦し，そのための圧力損失が考えられるが，後に述べる計算法に従えば，本報告の実験例では微小で、あることが解ったので省略する。観測きれた諸値と，他の研究者の実験式から，妥当な仮定を置いて式(1)の各項を数値的に見積り，更に(1)の 2， 3， 4項については粒子のランダ、ム運動過程を考慮、して，それらの事象が確率的にある割合しか起らないと仮定して，確率的係数 P 2，P 3，P 4をそれぞれ掛けることとした。 ( 3 2 ) 円管内における粒体の水カ輸送 4 2 1 第 7報流動様式ごとの実験例を選び出し，式(1)を次のような計算式により数値化し， 3 個以上の測定方程式を立てて最小自乗法により P 2，P 3，P 4を求める。それらが求められたならば更めて 2，3， 4項の値を求め， 1項内を占有する 2， 3，4，5項の割合が計算される。次の式は一般化するために各項をそれぞれの実験の場合の粒子の重量で割って無次元化しである。 j b (与 ) ' .; .. ( 宏 ) プ 2. 長y { ( 与一 1)会+1}企詐づ+ρ3・告す 4 3 C.( ム・2: AU jー仙 2十 ] _ 会竺ム U-U.)2 +2 f__.(~)~ -J-.( l 'n nj¥U j - U.jI ，ρ j J. 4; r a r . g，.( U -U .I ， n.. ¥ a ) .Y s . . ¥的 dx)， j) w (2) C:管内粒子濃度， CDS1 流速が変動するときの粒子抗力係数，CDS2 :管内定常流の粒記号子抗力係数， Cv :吐出粒子濃度，D:管径， d:粒子(球)直径， !:転がり摩擦の腕(図- 4 )，g:重力の加速度， k:係数， L:管長，n 管長 L=D内の球数，均:図 -6の ( D / 6 区)xLの中の球数， ρ2，t3，ρ 4 .それぞれ 2， 3， 4項の確率的係数，d P /命:圧力損失，U 管内平均流速，t も:平均粒子速度，Uj図 -6の D/6区内の流速，的:向上内の粒子速度， γ:比重量， θ:図 -5の粒子堆積角， ρ:密度，。 i:全圧力勾配，加水のみの圧力勾配， s 粒子， w:水添字ー凶HM ド " 、p h41-4 式 (2)の各項の誘導は次の通りである。 nHVFhd 唱ZA (1) 全圧力損失項: ( π/4)D2L(φ/ 批);である。こ、で(ゆ/批)i は図 -3の D u r a n da n dC o n d o l i o sの値を修正したZ a n d ia n dG o v a t o sの値3) より求める。 2 1 2 1 0 0 の Cvは， )去+1} C v=1/ { (~ -1 た。 (2) 粒子転動摩擦項: 図 -4に示す転がり摩擦の考え方より F mg!/(d/2) 二とすれば，粒子転動摩擦項は k .( π/6).d3n.(ρsρ )g'2!/d 出 ( 3 3 ) φの中 4 2 2 奥田教海・山岸英明 F 寸ト f 図 -4. 転がり摩擦図ー 5 . 管底に堆積している球群となる。浮遊粒子が管底に堆積したと仮定したとき管壁に接触する粒子数を図 -5により計算 C・(D/d)Sを考慮すれば，浮遊状態で接触粒子数の nにし，L=D内の粒子数 n= (3/2)・対する割合は {2 ¥ {d， ¥( {D 8 ， k=ρ 2 .~ 3~Cハ万 J l\ d- 1 Jz一 +q 1 ¥ 1 ) (3) 粒子速度変動抗力項: /7乗則より求め，粒子速度図 -6に示す D/6区ごとに流速と粒子速度を考える。流速は 1 は，平均速度比 u s / uの値より各区ごとに想定し，半数の粒子は衝突によりその速度の 95%に減速されると仮定する。その速度とその区の流速の差より粒子速度変動抗力項を求めると， ~+nj" PS. CDS!. ( 7 r/4)d2.( P w / 2 ) . (Uj-USj)2 となる。流速分布 (4) 全粒子抗力項: 濃度比分布管断面の区分前項の外に，全区の粒子に対して，流れと粒子の平均速度差に関する粒子の抗力を考えると， ρ 4・CDs2・n (π / 4 )d2.(ρ / 2 )・(U-Us)2 却となる。 (5) 水のみの圧力損失項: 国 -6 . 粒子速度変動抗力項の計算 π (/4) D2L.( φ/命)加である。以上を式(1)の形にまとめ，両辺を n ( 1 C /6) d S y sで割って無次元化すると，式(2)が得られる。Co s ! = 2 . 0， CDS2=1 .5 ，/=d/4 と仮定すれば，実験結果より，式(2)の各項の値があ， / J s ， P 4を除き求められる。 ( 3 4 ) 円管内における粒体の水力輸送 4 2 3 第 7報 I V .計算結果採り上げた実験例は， (1)低濃度管底流 (C=1 .5-4.8%) 18例 (2)高濃度管底流表 ~1 分類 (C=8.2-10.9%) 6例 (3)低濃度部分流 (C=2.1-7.7%) 12例 i j l [ 5 ] 確率的係数 P 2 P 3 0 . 5 5 8 0 . 7 8 8 0 . 4 7 0 0 . 9 3 8 0 . 7 4 3 0 . 0 6 3 0 . 0 3 3 0 . 0 1 7 P 4 O . ~0.033 ~0.084 (4)低濃度均一流 (C=2.8-6.8%) 1 0例 (5)高濃度均一流。 OIl-A A T - t o FBEEトEEE﹄各項の割合% 各項の割合の見方 1 0 0 ?よ ∞ o (C=12.5-14.2%J 4例、、、、 4 。 r " ， 5 デ4 1 5 1 0 F r o 図一 7 .低濃度管底流 c =1. 5~4.8% I D O 各項の割合% 5 0 。 1 0 1 5 丹d 2 . 0 図 - 8低濃度部分流， C=2.1~7.7% ( 3 5 ) 。 0 . 0 8 6 。 1I l 1 0 . 1 1 3 4 2 4 奥田教海・山岸英明 1 0 0 各項の割合% $ 0 。 20 阪 F r d 25 図 -9 低濃度均一流， 30 C=2.8-6.8% -管底流高濃度均一流 0部分流高濃度管底流。均一流割合% 低濃度管底流。低濃度部分流 5 F r d 低濃度均一流。 3 国一 1 0 粒子転動摩擦項の割合である。最小自乗法による計算の結果，確率的係数は表 -1のように求められた。この九九九の値を用いて更めて各項の値を出した。各項の全圧力損失項内に占める割合を求め図示したも 0である。のは，図 7ー図 1 0 )，また転動摩擦項の割合は，管底流，部分流，均一流へと遷移するに従い小さくなり(図 1 管内濃度の高い場合は大きくなることが解る。また (2J ， (3J ， (5Jの場合， 3， 4項に負の値が出ていることは，粒子がランダム運動を行っていることを考慮すれば，粒子と流体との聞に運動量の交換が行われているものと考えられるが，これについては，なお今後検討を要する。 ( 3 6 ) 円管内における粒体の水力輸送 v .まと 4 2 5 第 7報め粒子浮遊状態の観測から付加圧力損失の内容の変化を，いくつかの妥当な仮定をおくことにより，求めることができた。水力輸送管内の粒子流動機構を探る一知見としたい。 6年 5月 2 0日受理) (昭和 5 文献 1) 奥田教海，山岸英明ー室蘭工業大学研究報告， 9，( 2 )5 3 5( 1 9 7 7 ) 2) 奥田教海室蘭工業大学研究報告， 8，( 3 )6 7 1( 1 9 7 6 ) I . :Advancesi nS o l i d L i q u i dFlowi nP i p e sandi t sA p p l i c a t i o n， p . 1 2 9(PergamonP r e s s， 1 9 7l ) 3 ) Zandi， ( 3 7 )