(57)【要約】【課題】撹拌反応槽における気泡の実際の挙動及び撹拌翼の撹拌動力の低下を考慮して、従来法よりも精度の高い気液界面積の算出方法を提供すること。【解決手段】本発明の気液界面積の算出方法は、撹拌翼を備える気液撹拌槽において、液相と、液相に分散する気泡との気液界面積の算出方法であり、相間交換係数を算出する工程と、補正関数を算出する工程と、補正相間交換係数を算出する工程と、補正相間交換係数より算出される撹拌動力計算値と撹拌動力実測値とを対比し、その結果に基づいて調整補正関数を算出する工程と、調整補正関数に基づいて気液界面積を算出する工程と、を含むものであり、補正関数としてη＝ｍａｘ［ｍｉｎ｛（１＋ｑ−ｒｄ）ｐ，１｝，ｍ］（ηは補正関数、ｒｄは気泡の体積分率である。また、０．１≦ｑ≦０．５、０．１≦ｍ≦０．３、ｐ＝−１，−０．５，０，２，４のいずれかである。）の式を用いる。【選択図】図４
2019141797.tif 000002
【特許請求の範囲】
【請求項１】
撹拌翼を備える気液撹拌槽において、液相と、該液相に分散する気泡との気液界面積の算出方法であって、
下記（１）式〜（３）式より、相間交換係数Ｋを算出する相間交換係数算出工程と、
下記（４）式に示す補正関数の式より、補正関数ηを算出する補正関数算出工程と、
下記（５）式に示す補正相間交換係数の式より、補正相間交換係数Ｋ’を算出する補正相間交換係数算出工程と、
前記補正相間交換係数Ｋ’より算出される撹拌動力計算値と前記撹拌動力実測値とを対比し、該対比の結果に基づいて調整補正関数η’を算出する調整補正関数算出工程と、
前記調整補正関数η’に基づいて気液界面積を算出する気液界面積算出工程と、
を含む
気液界面積の算出方法。
【数１】
2019141797.tif 000018
（上式において、ρは液相の密度、ｄは気泡径、Ａは気泡表面積、ｆは抵抗関数、τは緩和時間である。）
【数２】
2019141797.tif 000019
（上式において、Ｃｄは抵抗係数、Ｒｅはレイノルズ数である。）
【数３】
2019141797.tif 000020
（上式において、ρは液相の密度、ｄは気泡径、μは液相の粘度である。）
【数４】
2019141797.tif 000021
（上式において、ηは補正関数、ｒｄは気泡の体積分率である。また、０．１≦ｑ≦０．５、０．１≦ｍ≦０．３、ｐ＝−１，−０．５，０，２，４のいずれかである。）
【数５】
2019141797.tif 000022
（上式において、Ｋ’は補正相間交換係数、Ｋは相間交換係数、ηは補正関数である。）
【請求項２】
撹拌翼及び液相へのガス吹込み口を備える気液撹拌槽における前記ガス吹込み口の位置設計方法であって、
請求項１に記載の気液界面積の算出方法により算出される前記気液界面積が所定の値となるように、前記ガス吹込み口の位置を設計する
ガス吹込み口の位置設計方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、気液界面積の算出方法及びガス吹込み口の位置設計方法に関する。
【背景技術】
【０００２】
液中に気泡が分散する気液反応系においては、気液界面積がその反応性に大きく影響を及ぼすため、気液界面積を正確に算出する方法が求められている。この気液界面積を正確に算出して最適化するように反応槽を設計することができれば、反応効率が高まることが期待される。
【０００３】
しかしながら、特に、非鉄金属を分散させたスラリーに空気を吹き込んで撹拌し、加圧下で酸素と反応させて非鉄金属を液中に浸出させる、いわゆる加圧浸出等に用いられる、密閉型の加圧型反応槽では、槽内の気泡の分散状態を直接観察することが困難である。したがって、このような系においては特に気液界面積の算出方法が重要であるが、その方法は充分に確立されていない。
【０００４】
液中に気泡が分散する系における計算には、相間抗力、揚力、壁面潤滑力等を考慮するが、このうち相間抗力のモデルとしては、例えばＳｃｈｉｌｌｅｒ−Ｎａｕｍａｎｎの式（非特許文献１）、Ｇｒａｃｅの式（非特許文献２）等が用いられている。このうち、Ｓｃｈｉｌｌｅｒ−Ｎａｕｍａｎｎの式は、流体粒子を剛体球とするモデルである。また、Ｇｒａｃｅの式は、流体粒子の楕円変形を考慮するモデルである。しかしながら、これらのモデルは一様な流れの中の１つの気泡を対象としてモデル化されたものであるため、撹拌され、且つ気泡が多数存在する条件において、このようなモデルでは気液撹拌槽中の気液界面積を求めるには充分でない。
【０００５】
また、液中に気泡が分散していると、撹拌翼の回転に伴う遠心力によって、密度の低い気泡には中心へ向かう力が働き、通気量が多いほど撹拌翼近傍のガスの割合が大きくなる。このようにして撹拌翼近傍のガスの割合が大きくなると、気泡が撹拌翼にまとわりつく、いわゆるフラッディング現象が起こり、撹拌動力が大きく低下し、それによって気泡の挙動も変化するが、このような撹拌動力の低下を考慮したモデルはこれまでない。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】Ｌ．Ｓｃｈｉｌｌｅｒ  ａｎｄ  Ａ．Ｎａｕｍａｎｎ，ＶＤＩ  Ｚｅｉｔｓ，７７，３１８（１９３３）．
【非特許文献２】Ｒ．Ｃｌｉｆｔ，Ｊ．Ｒ．Ｇｒａｃｅ，Ｍ．Ｅ．Ｗｅｂｅｒ，“Ｂｕｂｂｌｅｓ，Ｄｒｏｐｓ  ａｎｄ  Ｐａｒｔｉｃｌｅｓ”，ＡｃａｄｅｍｉｃＰｒｅｓｓ，Ｎｅｗ  Ｙｏｒｋ，Ｕ．Ｓ．Ａ．（１９７８）．
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、以上のような実情に鑑みてなされたものであり、撹拌反応槽における気泡の実際の挙動及び撹拌翼の撹拌動力の低下を考慮して、従来法よりも精度の高い気液界面積の算出方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本発明者らは、上述した課題を解決するために鋭意検討を重ねた。その結果、相間抗力として気泡の分裂・合体を考慮した補正関数を用いて撹拌動力計算値を求め、撹拌動力実測値と対比して調整した補正関数に基づき気液界面積を算出することにより、従来よりも精度の高い気液界面積の算出方法を提供することができることを見出し、本発明を完成するに至った。具体的に、本発明は、以下のものを提供する。
【０００９】
（１）本発明の第１の発明は、撹拌翼を備える気液撹拌槽において、液相と、該液相に分散する気泡との気液界面積の算出方法であって、下記（１）式〜（３）式より、相間交換係数Ｋを算出する相間交換係数算出工程と、下記（４）式に示す補正関数の式より、補正関数ηを算出する補正関数算出工程と、下記（５）式に示す補正相間交換係数の式より、補正相間交換係数Ｋ’を算出する補正相間交換係数算出工程と、前記補正相間交換係数Ｋ’より算出される撹拌動力計算値と前記撹拌動力実測値とを対比し、該対比の結果に基づいて調整補正関数η’を算出する調整補正関数算出工程と、前記調整補正関数η’に基づいて気液界面積を算出する気液界面積算出工程と、を含む気液界面積の算出方法である。
【００１０】
【数１】
2019141797.tif 000003
（上式において、ρは液相の密度、ｄは気泡径、Ａは気泡表面積、ｆは抵抗関数、τは緩和時間である。）
【００１１】
【数２】
2019141797.tif 000004
（上式において、Ｃｄは抵抗係数、Ｒｅはレイノルズ数である。）
【００１２】
【数３】
2019141797.tif 000005
（上式において、ρは液相の密度、ｄは気泡径、μは液相の粘度である。）
【００１３】
【数４】
2019141797.tif 000006
（上式において、ηは補正関数、ｒｄは気泡の体積分率である。また、０．１≦ｑ≦０．５、０．１≦ｍ≦０．３、ｐ＝−１，−０．５，０，２，４のいずれかである。）
【００１４】
【数５】
2019141797.tif 000007
（上式において、Ｋ’は補正相間交換係数、Ｋは相間交換係数、ηは補正関数である。）
【００１５】
（２）本発明の第２の発明は、撹拌翼及び液相へのガス吹込み口を備える気液撹拌槽における前記ガス吹込み口の位置設計方法であって、第１の発明に係る気液界面積の算出方法により算出される前記気液界面積が所定の値となるように、前記ガス吹込み口の位置を設計する、ガス吹込み口の位置設計方法である。
【発明の効果】
【００１６】
本発明によれば、撹拌反応槽における気泡の実際の挙動及び撹拌翼の撹拌動力の低下を考慮して、従来よりも精度の高い気液界面積の算出方法を提供することができる。
【図面の簡単な説明】
【００１７】
【図１】本実施の形態に係る気液界面積の算出方法を示す図である。
【図２】撹拌動力計算値及び撹拌動力の実測値の撹拌動力値対撹拌回転数のプロットである。
【図３】本実施の形態に係る気液撹拌槽の断面模式図である。
【図４】実施例１、比較例１及び２の気液撹拌層における補正関数対気泡の体積分率のプロットである。
【図５】実施例１及び比較例１の撹拌動力対撹拌回転数のプロットである。
【発明を実施するための形態】
【００１８】
以下、本発明の具体的な実施形態（以下「本実施の形態」という）について詳細に説明するが、本発明は以下の実施形態に何ら限定されるものではなく、本発明の要旨を変更しない範囲内において、適宜変更を加えて実施することができる。
【００１９】
なお、本明細書において、「ｍａｘ（Ｘ，Ｙ）」、「ｍｉｎ（Ｘ，Ｙ）」とは、それぞれＸ及びＹのうち、数値の大きいもの、小さいものとする。
【００２０】
≪１．気液界面積の算出方法≫
気液界面積の算出方法は、撹拌翼を備える気液撹拌槽において、液相と、その液相に分散する気泡との気液界面積の算出する方法である。
【００２１】
図１は、本実施の形態に係る気液界面積の算出方法を示す図である。図１に示すように、本発明に係る気液界面積の算出方法は、相間交換係数算出工程Ｓ１と、補正関数算出工程Ｓ２と、補正相間交換係数算出工程Ｓ３と、調整補正関数算出工程Ｓ４と、気液界面積算出工程Ｓ５とを含むものである。以下、気液界面積の算出方法における各工程について、図１を参照して順を追って説明する。
【００２２】
（相間交換係数算出工程）
相間交換係数算出工程Ｓ１は、下記（１）式〜（３）式より、相間交換係数Ｋを算出する工程である。
【００２３】
【数６】
2019141797.tif 000008
（上式において、ρは液相の密度、ｄは気泡径、Ａは気泡表面積、ｆは抵抗関数、τは緩和時間である。）
【００２４】
【数７】
2019141797.tif 000009
（上式において、Ｃｄは抵抗係数、Ｒｅはレイノルズ数である。）
【００２５】
【数８】
2019141797.tif 000010
（上式において、ρは液相の密度、ｄは気泡径、μは液相の粘度である。）
【００２６】
ここで、気泡径ｄとしては、任意の値を選択する。例えば、ガス吹込み口を用いて液相にガスを供給するような系では、ガス吹込み口の径の大きさやガス吹込みの速度を考慮して算出することができる。また、予備実験をして気泡径ｄを測定してもよい。
【００２７】
また、気泡表面積Ａとしては、気泡を剛体球と仮定して、気泡径ｄから比表面積を求める。
【００２８】
（補正関数算出工程）
補正関数算出工程Ｓ２は、下記（４）式に示す補正関数の式より、補正関数ηを算出する工程である。具体的に、補正関数算出工程Ｓ２においては、ｐ、ｑ及びｍとして、任意の値を選択して補正関数ηを算出する。
【００２９】
【数９】
2019141797.tif 000011
（上式において、ηは補正関数、ｒｄは気泡の体積分率である。また、０．１≦ｑ≦０．５、０．１≦ｍ≦０．３、ｐ＝−１，−０．５，０，２，４のいずれかである。）
【００３０】
これに対し、従来のＧｒａｃｅモデルにおける補正関数は下記（４’）式のように示される。
【００３１】
【数１０】
2019141797.tif 000012
【００３２】
上記（４）式の技術的意義について、図４を参照して、従来の（４’）式のＧｒａｃｅモデルと対比しながら説明する。図４は気液撹拌層における補正関数対気泡の体積分率のプロットである。図４において、実施例１が本発明のモデル、比較例１がＧｒａｃｅモデルである。比較例１のＧｒａｃｅモデルにおいては、気泡の体積分率ｒｄが増加するごとに、補正関数が減少し、これに伴い補正関数を係数とする抵抗係数が減少するものとしている。
【００３３】
これに対し、実施例１の本発明のモデルにおいて、補正関数対気泡の体積分率のプロットは、図４に示すようにシグモイド型のような形状を有している。なお、図４においては、ｑ＝０．３、ｍ＝０．３である。本発明のモデルでは、上記（４）式のうち、ｑの項を設けることにより、ｒｄの値が所定の値ｑ（図４において０．３）を超えるまでは補正関数が低下しないものとした。これは、気泡の密集度が低いときすなわちｒｄが小さいときは気泡同士の影響が無いとして、液相中の気泡の合体や分裂は、液相や気相（気泡）の粘度や密度に影響を受けるため、この影響を反映させることを目的としている。例えば、液相の粘度や密度が大きい場合や気泡が少ない場合、気泡の合体はあまり起こらず、これによって補正関数及び抵抗係数が低下しないとしたものである。また、気泡が増加しｒｄの値がさらに大きくなった場合には、補正関数の最小値をｍとすることで、補正関数ηが所定の値ｍ（図４において０．３）よりも低下しないものとした。気泡が多い場合、気泡の合体が起こり液相からの影響は小さくなるものの、液相からの影響を受け続けるため、補正関数及び抵抗係数がゼロにはならないものとした。
【００３４】
このようにして、上記（４）式において、Ｇｒａｃｅの式と同様に気泡の楕円変形を考慮するとともに、さらにｑ、ｍの変数を設け、これらの変数によって気泡の分離及び合体を考慮する。このようにして、相間交換係数Ｋに補正を加えることで、後段の撹拌動力計算値算出工程で気泡の分離及び合体が考慮された撹拌動力計算値が算出され、これに基づき気液界面積を精度良く算出することができる。
【００３５】
（補正相間交換係数算出工程）
補正相間交換係数算出工程Ｓ３は、下記（５）式に示す補正相間交換係数の式より、補正相間交換係数Ｋ’を算出する工程である。具体的に補正相間交換係数Ｋ’は、相間交換係数算出工程において算出した相間交換係数Ｋと、補正関数算出工程において算出した補正関数ηの積である。
【００３６】
【数１１】
2019141797.tif 000013
（上式において、Ｋ’は補正相間交換係数、Ｋは相間交換係数、ηは補正関数である。）
【００３７】
（調整補正関数算出工程）
調整補正関数算出工程Ｓ４は、補正相間交換係数Ｋ’より算出される撹拌動力計算値と撹拌動力実測値とを対比し、その対比の結果に基づいて調整補正関数η’を算出する工程である。
【００３８】
ここで、「撹拌動力計算値」及び「撹拌動力実測値」とは、それぞれ気液撹拌槽の撹拌翼の撹拌動力の計算値及び実測値のことをいう。
【００３９】
撹拌動力計算値は、液相及び気泡それぞれにおいて、下記（６）式に示す連続の方程式、下記（７）式に示す運動量の方程式及び下記（８）式に示す相間抗力の式を解き、撹拌翼に作用する力を求め、これに基づき算出されたトルクから算出する。具体的には、トルク（Ｎｍ）に回転数（ｒｐｍ）を掛けて、９５４９で割ることによって、撹拌機の仕事率すなわち撹拌動力（ｋＷ）を算出する。
【００４０】
【数１２】
2019141797.tif 000014
（上式において、αは体積分率、ρは密度、ｕは流速である。また、ｋはＧ又はＬであり、ＧとＬはそれぞれ気泡、液相における変数を意味するものとする。）
【００４１】
【数１３】
2019141797.tif 000015
（上式において、αは体積分率、ρは密度、ｕは流速、Ｐは圧力、τは応力ひずみテンソル、ｇは重力加速度、Ｍは相間抗力、Ｆは揚力、ｋはＧ又はＬであり、ＧとＬはそれぞれ気泡、液相における変数を意味するものとする。）
【００４２】
【数１４】
2019141797.tif 000016
（上式において、Ｍは相間抗力、Ｋは相間交換係数、ｕは流速である。）
【００４３】
また、撹拌動力実測値は、撹拌機のモーターに接続された電力計によって測定する。なお、「撹拌動力」とは、撹拌機モーターが単位時間あたりの流体に与えるエネルギーと定義される。
【００４４】
攪拌動力計算値と攪拌動力実測値との対比方法としては、特に限定されるものではない。例えば、撹拌回転数２点以上におけるそれぞれの、撹拌動力計算値と撹拌動力実測値を求める。撹拌動力実測値を、撹拌回転数の３乗に比例する関数として近似する。図２は、撹拌動力計算値及び撹拌動力の実測値の撹拌動力値対撹拌回転数のプロットの一例である。この図２のように、撹拌動力計算値と撹拌動力の実測値が概ね重なるようであれば、両者は相関があると判断できる。そして、両者に相関があると判断した場合における補正関数ηを調整補正関数η’とする。
【００４５】
その相関の判断の基準としては特に限定されないが、例えば任意の撹拌回転数において、撹拌動力実測値に対する、撹拌動力実測値と撹拌動力計算値の差の絶対値が所定の閾値（例えば、１０％以下、５％以下）であるか否かにより判断する。そして、閾値以下である場合、撹拌動力計算値に相関があるものとする。
【００４６】
一方で、決定係数が所定値未満である場合、上記（４）式におけるｐ、ｑ及びｍの少なくともいずれかを変更し、再度、補正関数算出工程からこの調整補正関数算出工程までの工程を、決定係数が所定値以上となるまで繰り返す。なお、具体的には、撹拌動力計算値が撹拌動力実測値よりも大きい場合、抗力を小さくすべくｑを小さく、ｍを小さくする。一方で、撹拌動力計算値が撹拌動力実測値よりも小さい場合、抗力を大きくするようにｑを大きく、ｍを大きくする。そして、両者に相関があると判断した場合の補正関数ηを調整補正関数η’とする。
【００４７】
このように、撹拌動力計算値と撹拌動力実測値を対応させるように、上記（４）式の変数ｑ及びｍを変更する。これにより、ｑ及びｍが実際の実験系における撹拌動力実測値を反映したものとなり、また、気泡の存在による撹拌動力の低下（フラッディング）の影響を反映したものとなり、精度良く気液界面積を算出することができる。
【００４８】
（気液界面積算出工程）
気液界面積算出工程Ｓ４は、調整補正関数η’に基づいて気液界面積を算出する工程である。具体的には、上記（６）式〜（８）式の計算結果より、撹拌動力と同時に気泡の量、気泡径が求められ、気泡量を気泡径で割ることによって気液界面積を求める。
【００４９】
≪２．ガス吹込み口の位置設計方法≫
本実施の形態に係るガス吹込み口の位置設計方法は、撹拌翼及び液相へのガス吹込み口を備える気液撹拌槽におけるものであり、上述した気液界面積の算出方法により算出される気液界面積が所定の値となるように、ガス吹込み口の位置を設計するものである。
【００５０】
図３は、本実施の形態に係る気液撹拌槽の断面模式図である。気液撹拌槽１は、反応槽１１と撹拌翼１２とガス吹込み口１３とを備え、その反応槽１１の内部に液相Ｌが充填されている。また、液相Ｌ中にはガス吹込み口１３から気泡Ｇが供給され、これにより気泡Ｇが分散している。
【００５１】
なお、図３においては、反応槽１１として球状のものを用いているが、その形状は液相Ｌを格納し得るものであれば特に限定されず通常の円筒状のもの等を用いることができる。また、反応槽１１の素材としても、液相Ｌの水圧及び反応槽内に圧力を印加する場合にはその圧力に耐えるものであれば特に限定されず、例えばステンレス等、各種金属、合金等を用いることができる。さらに、撹拌翼１２の形状、位置及び素材、ガス吹込み口１３の形状、素材についても、特に限定されるものではなく、適宜選択することができる。
【００５２】
以上のような気液撹拌槽において、液相中に分散される気泡の挙動は、ガス吹込み口の位置に大きく影響を受ける。より具体的に、ガス吹込み口の位置が変化すると、気液撹拌槽内の気泡の分布が変わる。したがって、ガス吹込み口の位置の設計は非常に重要である。一方で、上述した気液界面積の算出方法によれば、精度良く気液界面積を算出することができる。例えば、気液撹拌槽のモデルで吹込み位置を、例えば撹拌機から距離を少しずつ大きくしたり、吹き込む角度を変えて計算したりすることで、槽内の吹込み位置の座標に対し気液界面積をプロットして、所定の値になる位置を特定する。そこで、ガス吹込み口を変化させて、それぞれの場合において上述した気液界面積の算出方法により気液界面積を算出する。このようにして算出した気液界面積のうち、所定の値を示す場合のガス吹込み口の位置に設計する。
【００５３】
なお、気液間の反応性を高くするならば気液界面積を最大とすればよく、また、反応性を低くするならば気液界面積を最小とすればよい。また、気液間の反応性に最適値があるならば、そのように気液界面積を決定すればよい。
【００５４】
このようなガス吹込み口の位置設計方法によれば、気液界面積が例えば最大となるガス吹込みを精度良く特定することができ、その気液撹拌槽における化学反応の効率を高めることができる。
【実施例】
【００５５】
以下、本発明の実施例を示してさらに詳細に説明するが、本発明はこれらの実施例に何ら限定されるものではない。
【００５６】
図１の模式図に示すような撹拌翼を備える球状の気液撹拌槽内に、液相として密度１４５０ｋｇ／ｍ３、粘度０．０１Ｐａ・ｓのスラリーを収容して、撹拌翼の側方に設置されたガス吹込み口から気相として密度：１１ｋｇ／ｍ３の空気を吹き込み、１８ＭＰａに加圧した気液撹拌槽において、撹拌翼の回転数を１３５ｒｐｍ、１４０ｒｐｍ、１４６ｒｐｍとした時の撹拌動力の実測値を測定した。回転数：１３５ｒｐｍを基準としてそれぞれの相対値を表１に示す。
【００５７】
【表１】
2019141797.tif 000017
【００５８】
（実施例１）
図１に示す気液撹拌槽について、オイラー・オイラー法で分散相における合体、分裂を考慮したポピュレーションバランスモデルを用いた連続相−分散相の混相流でモデル化を行い、液相（連続相）中の気相（分散相）の体積分率を解析した。具体的には、解析ソフトＡＮＳＹＳ  Ｆｌｕｅｎｔを用いて、図３に示すフローチャートに従って、液相と気相の連続の方程式と運動方程式を立て、撹拌翼の回転数１３５ｒｐｍ、１４０ｒｐｍ、１４６ｒｐｍにおける槽内の気相の体積分率の分布を解析した。ここで相間交換係数は上記（５）式を用い、Ｋの計算は（１）式〜（３）式に基づいて行った。なお、ガス吹込み口から供給された空気は直径５ｍｍの気泡として液相に導入されると仮定し、抵抗関数ｆは非特許文献２であるＧｒａｃｅの文献で開示された数式を用いた。
【００５９】
求められた気相の体積分率の分布における撹拌翼面の結果に基づき撹拌翼に及ぼすトルクを計算し、各回転数の撹拌動力を計算した。この計算結果と表１に示す実測値とを比較し、乖離が大きい場合には、（４）式に示す相間交換係数Ｋの補正関数ηのｑ及びｍを調整し、再度槽内の気相の体積分率の分布を解析し直し、その結果に基づいて各回転数における撹拌動力を再計算した。なお、本反応系は（４）式のｐの値は後述する比較例１において、ｐ＝２が最適であったのでこの値に固定し、ｑ及びｍの値のみを調整することとした。この操作を繰り返し、補正関数ηの変数を絞り込んだ結果、ｑ＝０．３、ｍ＝０．３が最適値であると判断し、その値に基づき調整補正関数η’を算出した。
【００６０】
（比較例１）
相間交換係数Ｋの算出にあたって上記（４）式の代わりに、Ｇｒａｃｅモデルの（４’）式を用い、且つ各回転数の撹拌動力の計算値と実測値の比較も行わない以外は、実施例１と同様にして各回転数における撹拌動力を計算した。なお、これは、（４）式においてｑ＝０、ｍ＝０とした場合に相当する。
【００６１】
（比較例２）
補正関数ηが気泡の体積増加に伴い減少せずη＝１で一定とし、且つ各回転数の撹拌動力の計算値と実測値の比較も行わない以外は、実施例１と同様にして各回転数における撹拌動力を計算した。なお、これは、（４）式においてｐ＝０、ｑ＝１、ｍ＝１とした場合に相当する。
【００６２】
図４は、実施例１、比較例１及び２の気液撹拌槽における補正関数対気泡の体積分率のプロットである。図４から、比較例１の補正関数ηは、気泡の体積分率ｒｄが増加するにしたがい、減少することが分かる。これに対し、実施例１の補正関数ηは、体積分率ｒｄが増加しても、ｒｄが０．３を超えるまで減少せず、０．３を超えて減少が開始する。また、さらにｒｄが増加してηが減少しても、０．３以下となることがない。なお、比較例２は、気泡の体積分率ｒｄに関わらず１である。
【００６３】
図５は、実施例１及び比較例１の撹拌動力対撹拌回転数のプロットである。実施例１では、計算モデルから得られた撹拌動力の相対値は実測値を精度良く再現できているのに対し、従来の計算モデルとした比較例１や比較例２では再現の精度が低い。
【符号の説明】
【００６４】
１    気液撹拌槽
１１    反応槽
１２    撹拌翼
１３    ガス吹込み口
【図１】
2019141797.tif 000023
【図２】
2019141797.tif 000024
【図３】
2019141797.tif 000025
【図４】
2019141797.tif 000026
【図５】
2019141797.tif 000027