生命科学関連特許情報
| タイトル： | 公開特許公報(A)_圧力損失推定方法、圧力損失推定装置、圧力損失推定プログラム、および該プログラムを記録した記録媒体 |
| 出願番号： | 2005084428 |
| 年次： | 2006 |
| IPC分類： | G01N 15/08 |
特許情報キャッシュ
遠藤 禎行 JP 2006266824 公開特許公報(A) 20061005 2005084428 20050323 圧力損失推定方法、圧力損失推定装置、圧力損失推定プログラム、および該プログラムを記録した記録媒体 住友化学株式会社 000002093 特許業務法人原謙三国際特許事務所 110000338 遠藤 禎行 G01N 15/08 20060101AFI20060908BHJP JPG01N15/08 D 11 1 ＯＬ 18 本発明は、単一粒子群を複数混合して成る粒子充填層である複数粒子群充填層に関して、該複数粒子群充填層を透過する流体の圧力損失を推定する圧力損失推定方法、圧力損失推定装置、圧力損失推定プログラム、および該プログラムを記録した記録媒体に関するものである。 粒子充填層（固定層）の間隙を、液体、気体などの流体が透過する現象は、工業的には種々の装置で発生する。これらの装置の例としては、固液を分離する濾過装置、集塵のためのバグフィルタ、触媒粒子層にガスを通して反応させる触媒反応装置、および、湿った粉体を通気乾燥させる粉体乾燥装置が挙げられる。 上記現象が発生する場合、流体の流れは、粒子充填層からの抵抗などにより損失する。この損失を圧力の次元で表して圧力損失という。上記現象が発生する装置は、粒子充填層の圧力損失を予測した上で設計する必要がある。さもなくば、例えば、通液や通気が低下したり、それによって、反応や乾燥が不十分になったりといった不具合が発生する。 なお、本願では、粒子の周りの流体の流れが層流である、すなわち、ストークス域（粒径基準のレイノルズ数が５以下である領域）に限定している。この場合の圧力損失の関係式としては、最も簡単なダルシー（Darcy）の法則が知られている。ダルシーの法則は次式で表される。ΔＰ／Ｈ＝（μ／ｋ）×ｕＳ ・・・（１）。ここで、μは粘性率（Ｐa・ｓ）であり、Ｈは粒子充填層の高さ（厚さ）（ｍ）であり、ｕＳは流体の空塔速度（見かけ速度）（ｍ／ｓ）である。また、ｋは粒子充填層の透過率であり、粒子充填層の物理的性質によって決まる無次元の定数である。 ダルシーの法則に対して、粒子充填層を１つの管とみなして管抵抗の概念を導入し、表面係数をφｃとして球形以外の充填物を対象とすることにより、コゼニー・カーマン（Kozeny-Carman）の式が得られる。コゼニー・カーマンの式は次式の通りである。ΔＰ／Ｈ＝１８０×（μ×ｕＳ／（φｃ２×ｄｐ２））×（（１−ε）２／ε） ・・・（２）。ここで、ｄｐは粒子の直径（粒径）（ｍ）であり、εは空隙率（無次元）である。通常はコゼニー・カーマンの式を用いて、圧力損失を理論的に予測している（例えば、非特許文献１・２を参照）。社団法人化学工学会編，「化学工学便覧」，改訂６版，丸善株式会社，１９９９年２月２５日，ｐ．３１０−３１２社団法人化学工学会編，「現代の化学工学ＩＩ」，初版，株式会社朝倉書店，１９８９年２月２０日，ｐ．３９−４１I. F. Macdonald, et al., Ind. Eng. Chem. Fundamentals, 18 (1979), 199-208M. J. Macdonald, et al., AIChE Journal, 37 (1991), 1583-1588Y. Li, et al., Ind. Eng. Chem. Res., 37 (1998), 2005-2011Y. Endo, et al., Filtration & Separation, March 1998, 191-195 しかしながら、コゼニー・カーマンの式は、粒径分布の狭い１種類の粒子からなる粒子充填層の場合にしか適用できないという問題点が存在する（例えば、非特許文献３を参照）。この問題点を解消するため、例えば非特許文献４〜６に記載のように、種々の圧力損失および透過率の関係式が提案されている。 非特許文献４には、ブレーク・コゼニー（Blake-Kozeny）の式を複数サイズの球形粒子に一般化することが開示されている。また、非特許文献４では、粒子形状が同じ球形であり、粒径分布が異なる３つの粒子群を混合したフリットガラスの粒子充填層について、空隙率および透過率を求める実験が行われている。また、非特許文献５には、種々の粒径分布を有する球形粒子から成る粒子充填層の透過率に関して開示されている。 このように、非特許文献４・５には、粒径分布の広い粒子から成る粒子充填層や、粒径分布の異なる複数の粒子群から成る粒子充填層に関する圧力損失および透過率の関係式が開示されている。しかしながら、非特許文献４・５では、粒子が球形である場合に限定されているため、粒子が砂などの非球形である場合には、上記圧力損失および透過率の関係式を満たさない虞がある。 一方、本発明者は、非特許文献６において、異なる粒径分布を有する２つの粒子群を混合した粒子充填層に関する圧力損失の関係式を提案している。非特許文献６では、粒子の形状は、球形粒子でなくてもよい。しかしながら、非特許文献６では、各粒子群の粒径分布が対数正規分布に従うと仮定している。 本発明は、上記の問題点に鑑みてなされたものであり、その目的は、単一粒子群を複数混合して成る複数粒子群充填層に関して、各単一粒子群の粒径分布および粒子形状に関係なく、圧力損失を推定できる圧力損失推定方法などを提供することにある。 上記の課題を解決するため、本発明に係る圧力損失推定方法は、単一粒子群を複数混合して成る粒子充填層である複数粒子群充填層に関して、該複数粒子群充填層を透過する流体の圧力損失ΔＰを推定する圧力損失推定方法であって、前記単一粒子群は、該単一粒子群から成る粒子充填層である単一粒子群充填層が、ΔＰ’／Ｈ’＝（μ’／ｋ’）×ｕＳ’ ・・・（１）’、（但し、Ｈ’：前記単一粒子群充填層の高さ、μ’：前記単一粒子群充填層を透過する流体の粘性率、ｋ’：前記単一粒子群充填層を透過する流体の透過率、ｕＳ’：前記単一粒子群充填層を透過する流体の空塔速度）を充足する単一粒子群であって、各単一粒子群充填層における前記比μ’／ｋ’を取得するステップと、前記複数粒子群充填層を透過する流体の体積流量Ｑを取得するステップと、前記複数粒子群充填層を所定寸法のセルに細分し、前記複数粒子群充填層を構成する複数の前記単一粒子群の構成比に従って、複数の前記単一粒子群に関する複数の前記単一粒子群充填層を、前記複数のセルにそれぞれ配置したセルモデルを構築するステップと、構築されたセルモデルの各セルに関して、隣接するセルから流入する流体、または隣接するセルに流出する流体が受ける抵抗Ｒを、前記セルの寸法と、前記セルに配置された前記単一粒子群充填層の前記比μ’／ｋ’とに基づいて算出するステップと、各セルに流体の合流点を有し、隣接する２つのセルの前記合流点間を流れる流体は、該２つのセルにそれぞれ配置された２つの前記単一粒子群充填層の前記抵抗Ｒを受ける流路モデルを構築するステップと、構築された流路モデルに対して、前記体積流量Ｑを利用することにより、前記流体が前記流路モデルを透過することにより生じる圧力損失ΔＰを算出するステップとを含むことを特徴としている。 上記の方法によると、単一粒子群を複数混合して成る複数粒子群充填層を、上式（１）’を充足する単一粒子群充填層が配置されるセルに細分化している。前記単一粒子群充填層では、比μ’／ｋ’が一定であるから、比μ’／ｋ’とセルの寸法とを用いて、該セルにて流体が受ける抵抗Ｒを求めることができる。これにより、抵抗Ｒが既知の流路モデルを構築することができ、流体の体積流量Ｑに基づいて圧力損失ΔＰを求めることができる。したがって、各単一粒子群充填層が上式（１）’を充足することを前提として、単一粒子群の種類や、粒径分布や、粒子形状に関係なく、複数粒子群充填層を透過する流体の圧力損失ΔＰを推定することができる。 なお、前記比μ’／ｋ’を取得するステップにおいて、前記比μ’／ｋ’は、前記式（１）’を利用して、前記流体の空塔速度ｕＳ’と、単位厚さ当りの圧力損失ΔＰ’／Ｈ’とから取得することが好ましい。この場合、単一粒子群充填層に対して、流体の透過特性試験を行って、前記流体の空塔速度ｕＳ’と、単位厚さ当りの圧力損失ΔＰ’／Ｈ’との比を求めることにより、容易に比μ’／ｋ’を取得することができる。 また、複数の単一粒子群が充分に混合された複数粒子群充填層は、単一粒子群充填層をランダムに配置した構成となっていることが考えられる。そこで、前記セルモデルを構築するステップにおいて、前記セルモデルは、複数の前記単一粒子群充填層を、前記複数のセルにそれぞれランダムに配置したセルモデルであることが好ましい。この場合、圧力損失ΔＰをより精度良く推定することができる。 また、前記複数粒子群充填層の寸法を取得するステップをさらに含んでおり、前記セルモデルを構築するステップにおいて、前記複数粒子群充填層の寸法に基づいて、前記複数粒子群充填層を所定寸法のセルに細分することが好ましい。この場合、種々の寸法の複数粒子群充填層に関して圧力損失ΔＰを推定することができる。 また、前記圧力損失ΔＰを算出するステップは、前記流路モデルに対して、前記流体の体積流量Ｑを電流値に対応付け、前記流体の合流点を接続点に対応付け、前記抵抗Ｒを電気抵抗に対応付け、前記圧力損失を電圧降下に対応付けて、キルヒホッフの法則を利用することにより、前記流路モデル全体の圧力損失を求め、これを前記複数粒子群充填層の圧力損失ΔＰとするステップであることが好ましい。 この場合、流路モデルを電気回路網に置き換えることができる。この電気回路網は、キルヒホッフの法則を利用することにより、複数の線形一次方程式に置き換えることができる。したがって、線形一次方程式の解法を含むコンピュータプログラムをコンピュータ上で実行することにより、電気回路網全体の電圧降下を容易に求めることができる。これにより、上記コンピュータプログラムを利用して、圧力損失ΔＰを容易に推定することができる。 また、前記セルは同一寸法の立方体であることが好ましい。この場合、隣接するセルから流入する流体、または隣接するセルに流出する流体が受ける抵抗Ｒは、隣接する方向に関わらず同一となるので、計算を簡略化することができる。 なお、圧力損失ΔＰを精度良く推定するには、セルの寸法を小さくして、セルの数を増やすことが好ましい。具体的には、セルの一辺の長さを５ｍｍ以下にすることが好ましく、約３．４ｍｍにすることがより好ましい。 本発明に係る圧力損失推定装置は、前記複数粒子群充填層を透過する流体の圧力損失ΔＰを推定する圧力損失推定装置であって、前記単一粒子群は、前記単一粒子群充填層が、上式（１）’を充足する単一粒子群であって、各単一粒子群充填層の前記比μ’／ｋ’を取得する単一粒子群充填層特性取得手段と、前記複数粒子群充填層を透過する流体の体積流量Ｑを取得する流体特性取得手段と、前記複数粒子群充填層を所定寸法のセルに細分し、前記複数粒子群充填層を構成する複数の前記単一粒子群の構成比に従って、複数の前記単一粒子群に関する複数の前記単一粒子群充填層を、前記複数のセルにそれぞれ配置したセルモデルを構築するセルモデル構築手段と、構築されたセルモデルの各セルに関して、隣接するセルから流入する流体、または隣接するセルに流出する流体が受ける抵抗Ｒを、前記セルの寸法と、前記セルに配置された前記単一粒子群充填層の前記比μ’／ｋ’とに基づいて算出する抵抗算出手段と、各セルに流体の合流点を有し、隣接する２つのセルの前記合流点間を流れる流体は、該２つのセルにそれぞれ配置された２つの前記単一粒子群充填層の前記抵抗Ｒを受ける流路モデルを構築する流路モデル構築手段と、構築された流路モデルに対して、前記体積流量Ｑを利用することにより、前記流体が前記流路モデルを透過することにより生じる圧力損失ΔＰを算出する圧力損失算出手段とを備えることを特徴としている。 上記の構成によると、単一粒子群を複数混合して成る複数粒子群充填層を、上式（１）’を充足する単一粒子群充填層が配置されるセルに細分化している。前記単一粒子群充填層では、比μ’／ｋ’が一定であるから、比μ’／ｋ’とセルの寸法とを用いて、該セルにて流体が受ける抵抗Ｒを求めることができる。これにより、抵抗Ｒが既知の流路モデルを構築することができ、流体の体積流量Ｑに基づいて圧力損失ΔＰを求めることができる。したがって、各単一粒子群充填層が上式（１）’を充足することを前提として、単一粒子群の種類や、粒径分布や、粒子形状に関係なく、複数粒子群充填層を透過する流体の圧力損失ΔＰを推定することができる。 なお、上記の圧力損失推定方法における各ステップを、圧力損失推定プログラムによりコンピュータ上で実行させることができる。さらに、前記圧力損失推定プログラムをコンピュータ読取り可能な記録媒体に記憶させることにより、任意のコンピュータ上で前記圧力損失推定プログラムを実行させることができる。 以上のように、本発明に係る圧力損失推定方法は、単一粒子群を複数混合して成る複数粒子群充填層を、上式（１）’を充足する単一粒子群充填層が配置されるセルに細分化することにより、単一粒子群の種類や、粒径分布や、粒子形状に関係なく、複数粒子群充填層を透過する流体の圧力損失ΔＰを推定できるという効果を奏する。 本発明の実施形態について説明する前に、本発明の圧力損失推定方法の概要について、図２〜図８を参照しつつ説明する。 上述のように、粒子充填層（固定層）を流体が透過するときの圧力損失に関する最も簡単な法則としては、次式で示されるダルシーの法則が知られている。ΔＰ／Ｈ＝（μ／ｋ）×ｕＳ ・・・（１）。また、ストークス近似が成立する条件下では、上述のコゼニー−カーマンの式（式（２））、およびルース（Ruth）の式は、例えば濾材、フィルタなど、粒子充填層を支持する物体の抵抗を無視すると、上記のダルシーの法則に従う。ダルシーの法則を利用すれば、後述する流体透過特性試験を一度行って、流体の空塔速度ｕＳと単位厚さ当りの圧力損失ΔＰ／Ｈとから比（μ／ｋ）を取得することにより、粒子充填層の任意の高さＨと流体の任意の空塔速度ｕＳとにおける圧力損失ΔＰを推定することができる。 空塔速度ｕＳは粒子充填層の断面領域における単位面積当りの体積流量に等しいから、式（１）に示されるダルシーの法則は次式で表すことができる。ΔＰ＝（μ／ｋ）×Ｈ×（Ｑ／Ａ） ・・・（３）。ここで、Ｑは流体の体積流量（ｍ３／ｓ）であり、Ａは粒子充填層の断面積（ｍ２）である。流体の空塔速度ｕＳの代わりに、流体の体積流量Ｑを用いる理由は後述する。 式（３）は変形すると次式となる。ΔＰ＝Ｒｐ×Ｑ （但し、Ｒｐ≡（μ／ｋ）×（Ｈ／Ａ）） ・・・（４）。この式（４）から、比（μ／ｋ）のみが、流体の粘性と、粒子充填層の特性（例えば、粒径分布、粒子の形状、および空隙率）とに依存することが理解できる。 即ち、上記比（μ／ｋ）は、流体が同じであり、かつ粒子充填層を構成する粒子群が同じである場合には一定であるが、粒子群の構成または構成比が変わる場合には変化することになる。この場合、圧力損失ΔＰを推定するには、流体透過特性試験を再度行って比（μ／ｋ）を取得する必要があり煩雑である。 この問題点に対し、本発明者は、粒子充填層を構成する複数の粒子群（以下「単一粒子群」と称する。）のそれぞれに関して、該単一粒子群から成る粒子充填層（以下「単一粒子群充填層」と称する。）の比（μ’／ｋ’）を流体透過特性試験などにより予め取得しておき、各単一粒子群充填層の比（μ’／ｋ’）に基づいて、単一粒子群を複数混合してなる粒子充填層（以下「複数粒子群充填層」と称する。）の比（μ／ｋ）を算出する方法を案出した。これにより、単一粒子群の構成または構成比が変わっても、単一粒子群充填層の比（μ’／ｋ’）が既知であれば、複数粒子群充填層の比（μ／ｋ）を算出できるので圧力損失ΔＰを容易に推定することができる。 以下、各単一粒子群充填層の比（μ’／ｋ’）と複数粒子群充填層の比（μ／ｋ）との関係について説明する。なお、以下の「粒子充填層」は、単一粒子群充填層および複数粒子群充填層の何れにも適用可能なものを示している。 まず、本発明者は、粒子充填層を電気回路に対応付けることを考えた。式（４）は、電気回路におけるオームの法則に類似している。オームの法則は、理想的な導体における電圧Ｖと電流Ｉとの関係を次式で示すものである。Ｖ＝Ｒｅ×Ｉ （但し、Ｒｅ≡ρ×（Ｌ／Ｓ）） ・・・（５）。ここで、Ｒｅは導体の電気抵抗（Ω）であり、ρは電気抵抗率（比抵抗）（Ω・ｍ）であり、Ｌは導体の長さ（ｍ）であり、Ｓは導体の断面積（ｍ２）である。比抵抗ρは、導体に固有のものである。 式（４）および（５）を比較すると、電圧降下Ｖは圧力損失ΔＰに対応し、電流Ｉは体積流量Ｑに対応し、電気抵抗Ｒｅは関数Ｒｐに対応し、比抵抗ρは比（μ／ｋ）に対応し、粒子充填層の高さＨは導体の長さＬに対応し、粒子充填層の断面積Ａは導体の断面積Ｓに対応することが理解できる。このように、ダルシーの法則をオームの法則に対応付けるため体積流量Ｑを使用している。なお、以下では、上記関数Ｒｐを「抵抗」と称する。 したがって、抵抗器を流れる電流は、対象となっている粒子充填層の断面領域を横切る全体の流れ、すなわち体積流量Ｑと等価である。また、上記の類似性から、比抵抗ρおよび比（μ／ｋ）は、材料の特性に関連するものであり、したがって一定であると考えられる。 図２（ａ）・（ｂ）は、上記対応付けを図にて示すものである。同図（ａ）は、流体が粒子充填層を透過するときの圧力損失ΔＰ、体積流量Ｑ、および抵抗Ｒｐを示しており、同図（ｂ）は、電流が抵抗器を流れるときの電圧降下Ｖ、電流Ｉ、および電気抵抗Ｒｅを示している。また、図３（ａ）は、流体が２つの粒子充填層を直列に透過する状態を示しており、同図（ｂ）は、流体が２つの粒子充填層を並列に透過する状態を示している。図２（ａ）・（ｂ）の関係を参照すると、直列接続した粒子充填層は、直列接続した抵抗器に対応付けることができ、並列接続した粒子充填層は、並列接続した抵抗器に対応付けることができると考えられる。 したがって、一方の粒子充填層における抵抗をＲｐ１とし、他方の粒子充填層における抵抗をＲｐ２とすると、直列接続した２つの粒子充填層（図３（ａ））の全抵抗ＲｐＴは、直列接続した２つの抵抗器の全抵抗の公式と同様に次式となる。ＲｐＴ＝Ｒｐ１＋Ｒｐ２ ・・・（６）。一方、並列接続した２つの粒子充填層（図３（ｂ））の全抵抗ＲｐＴは、並列接続した２つの抵抗器の全抵抗の公式と同様に次式となる。１／ＲｐＴ＝１／Ｒｐ１＋１／Ｒｐ２ ・・・（７）。なお、以下では、粒子充填層の抵抗Ｒｐと抵抗器の電気抵抗Ｒｅとを区別するために用いた添え字ｅ・ｐを省略する。 次に、本発明者は、複数粒子群充填層が小寸法の単一粒子群充填層を多数組み合わせたものであると考え、前記複数粒子群充填層を所定寸法のセルＣｉに細分し、前記複数粒子群充填層を構成する複数の前記単一粒子群の構成比に従って、複数の前記単一粒子群に関する複数の前記単一粒子群充填層を、前記複数のセルにそれぞれ配置したセルモデルを構築することを考えた。そして、単一粒子群充填層の比（μ’／ｋ’）と式（４）とから、各セルＣｉの上記抵抗Ｒｉを求め、求めた抵抗Ｒｉを用いて、複数粒子群充填層の全抵抗ＲＴを求め、求めた全抵抗ＲＴと体積流量Ｑとを式（４）に適用して圧力損失ΔＰを求めることを考えた。 図４は、粒子充填層を多数（図示の場合は８×８×８＝５１２個）のセルに分割したセルモデルの例を示している。また、図５は、各セルＣｉの流路モデルを回路記号で示している。各セルＣｉは、隣接するセルと相互に接続している。このため、各セルＣｉの流路モデルは、中心に合流点（ノード）ｎｏｄを有し、合流点ｎｏｄから隣接するセルに向けて抵抗Ｒｉが設けられた構成となっている。 また、図６は、隣接する複数のセルの流路モデルを回路記号で示しており、同図（ａ）は平面的に示したものであり、同図（ｂ）は立体的に示したものである。なお、同図（ａ）にて黒点で示される合流点と側壁Ｗとの間の抵抗器は除外されている。これは、粒子充填層から側壁への流れが存在しないとみなしているためである。 上記のように、各セルが、隣接するセルと相互に接続しているため、複数粒子群充填層の流路モデルは、抵抗が相互に接続した複雑な回路網となる。したがって、複数粒子群充填層の全抵抗ＲＴを求めることは困難である。そこで、抵抗器の複雑な電気回路網による電圧降下を求めるために利用されるオームの法則とキルヒホッフの第一法則とを、本願発明に適用している。 図７は、キルヒホッフの第一法則を説明するための簡単な例を示している。キルヒホッフの第一法則によると、或る接続点に流入する電流の総和は、上記接続点から流出する電流の総和に等しい。図示の例では、キルヒホッフの第一法則は次式に対応付けることができる。Ｑ１＝Ｑ２＋Ｑ３ ・・・（８）。また、接続点に流入する電流を正とし、接続点から流出する電流を負とすると、キルヒホッフの第一法則は、或る接続点における電流の総和はゼロであるといえる。したがって、図示の例では、キルヒホッフの第一法則は次式に対応付けることができる。Ｑ１−Ｑ２−Ｑ３＝０ ・・・（９）。上式（８）または上式（９）を各合流点ｎｏｄに適用することにより、セルごとの線形一次方程式が得られる。 図８は、２次元の流路モデルを回路記号で示している。図示の場合では、合流点ａにおいて、キルヒホッフの第一法則は次式に対応付けることができる。Ｑ’(ｂ→ａ)＋Ｑ’(ｃ→ａ)＋Ｑ’(ｄ→ａ)＋Ｑ’(ｅ→ａ)＝０ ・・・（１０）。ここで、Ｑ’(ｂ→ａ)は、合流点ｂから合流点ａに流れる流体の体積流量を表し、その他の合流点に関しても同様である。 上式（１０）にオームの法則を対応付けると次式となる。（ΔＰ’(ｂ→ａ)／Ｒ(ｂ→ａ)）＋（ΔＰ’(ｃ→ａ)／Ｒ(ｃ→ａ)）＋（ΔＰ’(ｄ→ａ)／Ｒ(ｄ→ａ)）＋（ΔＰ’(ｅ→ａ)／Ｒ(ｅ→ａ)）＝０ ・・・（１１）。ここで、ΔＰ’(ｂ→ａ)は、合流点ｂから合流点ａまでの圧力損失を表し、その他の合流点に関しても同様である。また、Ｒ(ｂ→ａ)は、合流点ｂから合流点ａまでの全抵抗を表し、その他の合流点に関しても同様である。 次に、合流点ａ〜ｅのそれぞれにおける圧力をＰ’ａ〜Ｐ’ｅとする。また、合流点ｂから合流点ａまでには、合流点ｂに対応するセルの抵抗Ｒｂと、合流点ａに対応するセルの抵抗Ｒａとが含まれるとする。この場合、上式（１１）は次式となる。（（Ｐ’ｂ−Ｐ’ａ）／（Ｒｂ＋Ｒａ））＋（（Ｐ’ｃ−Ｐ’ａ）／（Ｒｃ＋Ｒａ））＋（（Ｐ’ｄ−Ｐ’ａ）／（Ｒｄ＋Ｒａ））＋（（Ｐ’ｅ−Ｐ’ａ）／（Ｒｅ＋Ｒａ））＝０ ・・・（１２）。 なお、各セルは、構成が均一である最小の体積単位とすることが望ましい。また、抵抗Ｒは、比（μ’／ｋ’）だけでなく、断面積に対する高さの比（Ｈ’／Ａ’）にも依存するので、セルは、内部で垂直方向の抵抗と水平方向の抵抗とが同じとなるように、正方形または立方体であることが好ましい。この場合、上式（１２）のように、合流点ａの周囲の抵抗全てがＲａに等しくなる。また、上式（１２）は、３次元の流路モデルに容易に拡張できる。具体的には、図８の記載された図面の法線方向に隣り合う２つの合流点に関する２つの項を上式（１２）に追加すればよい。 回路網に対応付けられる流路モデル、すなわち粒子充填層の複雑さが増すにつれて、線形一次方程式および未知数の数が急激に増加する。したがって、これらの線形一次方程式の解法を含むコンピュータプログラム（例えば、フォートラン（Compaq Visual Fortran Version 6 (Hewlett-Packard Development Company, L.P.)）で書かれたプログラム）をコンピュータ上で実行することにより、上記線形一次方程式の解を求めることが望ましい。この場合、流路モデル全体の圧力損失ΔＰを容易に求めることができる。 ユーザは、例えば、複数粒子群充填層の寸法および体積流量Ｑと、各セルに配置される単一粒子群充填層の比（μ’／ｋ’）といったパラメータを入力すると、上記プログラムが上記線形一次方程式を生成して解くことを自動的に行う。上記プログラムの出力は、合流点間の圧力損失ΔＰ’および体積流量Ｑ’を含んでいる。したがって、流路モデルにおける上端部の合流点から下端部の合流点までの圧力損失を求めることにより、複数粒子群充填層における圧力損失ΔＰを求めることができる。なお、上記プログラムは、体積流量Ｑを変えて再実行することにより、圧力損失と体積流量との関係を容易に取得できる。 〔実施の形態〕 次に、本発明の一実施形態について図１および図９〜図１５を参照しつつ説明する。図９は、或る粒子充填層を流体が透過する時の圧力損失を推定する圧力損失推定システムの概要を示している。図示のように、圧力損失推定システム１は、流体透過特性試験機２および圧力損失推定装置３を備える構成である。これらの装置２・３は、一般的なコンピュータと同様の構成を有している。なお、本実施形態では、流体として圧縮乾燥空気を用いているが、液体などその他の流体に対しても同様に適用することができる。 流体透過特性試験機２は、混合される各単一粒子群に関して、該単一粒子群から成る単一粒子群充填層の流体透過特性を測定して出力するものである。また、圧力損失推定装置３は、流体透過特性試験機２が測定した各単一粒子群の流体透過特性を用いて、各単一粒子群を混合して成る複数粒子群充填層における圧力損失を推定して出力するものである。以下、上記装置２・３のそれぞれについて詳述する。 なお、上記の装置２・３おいて、外部に情報を出力する形態としては、任意の形態を利用することができ、例えば、印刷媒体にデータを印刷する形態、記録媒体にデータを記録する形態、および、通信媒体を介して外部の装置にデータを送信する形態が挙げられる。同様に、外部から情報を入力する形態としては、任意の形態を利用することができ、例えば、印刷媒体のスキャニングを行う形態、記録媒体からデータを読み出す形態、および、外部の装置から通信媒体を介してデータを受信する形態が挙げられる。 まず、流体透過特性試験機２の詳細について図１０を参照しつつ説明する。図１０は、流体透過特性試験機２の概略構成を示している。図示のように、流体透過特性試験機２は、充填塔１１、ポンプ１２、圧力計１３、制御装置１４、流量計１５、およびバグフィルタ１６を備える構成である。 充填塔１１は、断面が正方形の筒状であり、略鉛直方向に立設している。充填塔１１内の中央から稍下方には分散板２０が設けられ、分散板２０の上には粒子充填層２１が配置される。分散板２０は、粉体が下方へ落下することを防止する一方で、圧縮乾燥空気を下方に透過させるものである。分散板２０には、例えば、多孔板や網が利用される。 また、充填塔１１内における分散板２０の下方と粒子充填層２１の上方とには、空気圧を検知して電気信号に変換する感圧素子２２・２３がそれぞれ設けられている。感圧素子２２・２３は、変換した電気信号を外部の圧力計１３に出力する。 なお、充填塔１１としては、乾燥空気を粉体に通して粉体を乾燥させる乾燥装置や気相で粉体を触媒と反応させる反応装置、気液または固液を分離する濾過装置などを流用することができる。 ポンプ１２は、充填塔１１に乾燥した圧縮空気（以下、「乾燥圧縮空気」と略称する。）を供給するものである。ポンプからの乾燥圧縮空気は、流量計１５を介して、充填塔１１の上端から内部に供給される。流量計１５は、ポンプからの圧縮乾燥空気の流量を測定するものである。 圧力計１３は、充填塔１１内に設けられた感圧素子２２・２３からの電気信号に基づいて、粒子充填層２１の下方の空気圧に対する分散板２０の上方の空気圧の差圧を求めるものである。すなわち、この差圧は、粒子充填層２１を透過する前の空気と、粒子充填層２１を透過した後の空気との圧力損失ΔＰとなる。 制御装置１４は、流体透過特性試験機２内の各種構成を統括的に制御するものである。バグフィルタ１６は、充填塔１１から排出される空気の流路に設けられており、充填塔１１から排出される空気に含まれる粉体を、フィルタバッグを用いて濾過収集するものであり、粉体が外部に漏れることを防止するものである。 次に、圧力損失推定装置３の詳細について、図１および図１１を参照しつつ説明する。図１１は、圧力損失推定装置３の概略構成を示している。図示のように、圧力損失推定装置３は、入力部３０、制御部３１、記憶部３２、および出力部３３を備える構成である。 入力部３０は、入力デバイス、スキャンデバイス、再生デバイス、受信デバイスなどを用いて、流体透過特性試験機２からの情報の入力を受け付けるものである。また、入力部３０は、入力された情報を制御部３１に送信する。 なお、入力デバイスの例としては、キーボード、テンキー、カーソルキー、例えばマウスなどのポインティングデバイス、およびタッチパネルが挙げられる。また、再生デバイスは、フレキシブルディスク、光ディスク、光磁気ディスクなどの外部記録媒体に記録されたデータを再生するものであり、例えば、フレキシブルディスクドライブ、ＣＤ−ＲＯＭドライブ、ＤＶＤドライブ、ＭＯドライブなどが挙げられる。 制御部３１は、圧力損失推定装置３における各構成を統括的に制御するものである。なお、制御部３１の詳細な構成については後述する。 記憶部３２は、各種情報を記憶する機能を有するものである。記憶部３２の例としては、ＲＡＭ（Random Access Memory）、フラッシュメモリなどの半導体メモリと、ハードディスクとが挙げられる。本実施形態では、記憶部３２は、セルの寸法を記憶している。 出力部３３は、制御部３１から受信した圧力損失ΔＰなどのデータを、印刷出力デバイス、記録デバイス、送信デバイスなどを用いて出力するものである。なお、記録デバイスは、フレキシブルディスク、光ディスク、光磁気ディスクなどの外部記録媒体にデータを記録するものであり、例えば、フレキシブルディスクドライブ、ＣＤ−Ｒドライブ、ＭＯドライブなどが挙げられる。 次に、制御部３１の構成について説明する。図１１に示されるように、制御部３１は、数値取得部（単一粒子群充填層特性取得手段、流体特性取得手段）４０、セルモデル構築部（セルモデル構築手段）４１、抵抗決定部（抵抗算出手段）４２、流路モデル構築部（流路モデル構築手段）４３、および圧力損失算出部（圧力損失算出手段）４４を備える構成である。 数値取得部４０は、流体透過特性試験機２から入力部３０を介して、各単一粒子群に関する比μ’／ｋ’を取得するものである。また、数値取得部４０は、入力部３０から入力された複数粒子群充填層の寸法と、圧縮乾燥空気の体積流量Ｑ（ｍ３／ｓ）とを取得するものである。数値取得部４０は、複数粒子群充填層の寸法をセルモデル構築部４１に送信し、各単一粒子群に関する比μ’／ｋ’を抵抗決定部４２に送信し、かつ、体積流量Ｑを圧力損失算出部４４に送信する。 セルモデル構築部４１は、数値取得部４０からの複数粒子群充填層の寸法と、記憶部３２からのセルの寸法とに基づいて、セルの配列（数および配置）を決定するものである。また、セルモデル構築部４１は、決定したセルの配列情報を流路モデル構築部４３に送信する。例えば、複数粒子群充填層の寸法が縦６０×横６０×高さ６０（ｍｍ３）であり、セルの寸法が一辺３．３ｍｍの立方体である場合、セルを縦方向に１８個、横方向に１８個、高さ方向１８個それぞれ配列した、合計５８３２個のセルからなるセルモデルが構築されることになる。 抵抗決定部４２は、数値取得部４０からの各単一粒子群に関する比μ’／ｋ’と、記憶部３２からのセルの寸法とに基づいて、上述のように、或る単一粒子群からなるセルの抵抗Ｒを単一粒子群ごとに決定するものである。また、抵抗決定部４２は、決定した抵抗Ｒを流路モデル構築部４３に送信する。 流路モデル構築部４３は、セルモデル構築部４１からのセルの配列情報と、抵抗決定部４２からの抵抗Ｒとに基づいて、複数粒子群充填層の流路モデルを構築する。具体的には、まず、複数粒子群充填層を構成する複数の単一粒子群の構成比に従って、複数の単一粒子群に関する複数の単一粒子群充填層を、複数のセルにそれぞれランダムに配置する。そして、各セルに配置された単一粒子群充填層を、図５に示される流路モデルに変換し、隣接するセルの流路モデルどうしを接続することにより、複数粒子群充填層の流路モデルを構築する。流路モデル構築部４３は、構築された流路モデルに基づいて、上式（１２）のような各合流点における方程式を作成し、圧力損失算出部４４に送信する。 圧力損失算出部４４は、流路モデル構築部４３からの方程式に、数値取得部４０からの体積流量Ｑを代入して解くことにより、流路モデル全体の圧力損失ΔＰを算出し、これを複数粒子群充填層における圧力損失ΔＰとするものである。また、圧力損失算出部４４は、求めた圧力損失ΔＰを体積流量Ｑと共に出力部３３に出力する。 図１は、上記構成の圧力損失推定システム１が行う圧力損失ΔＰの推定処理を示している。図示のように、まず、流体透過特性試験機２が、個々の単一粒子群について流体（本実施形態では圧縮乾燥空気）の透過特性試験を行って、比μ’／ｋ’を取得する（ステップＳ１０（以下「Ｓ１０」と略称することがある。他のステップについても同様である。））。 次に、数値取得部４０が、複数粒子群充填層の寸法と、流体の体積流量Ｑと、流体透過特性試験機２が取得した各単一粒子群に関する比μ’／ｋ’とを取得する（Ｓ１１）。次に、セルモデル構築部４１が、複数粒子群充填層の寸法とセルの所定寸法とに基づいて、セルの数および配置を決定する（Ｓ１２）。また、抵抗決定部４２が、各単一粒子群に関する比μ’／ｋ’とセルの所定寸法とに基づいて、各セルの抵抗Ｒを決定する（Ｓ１３）。 次に、決定されたセルの数および配置と、決定された各セルの抵抗Ｒとに基づいて、流路モデル構築部４３が流路モデルを構築する（Ｓ１４）。次に、構築された流路モデルと、流体の体積流量Ｑとに基づいて、圧力損失算出部４４が複数粒子群充填層の圧力損失ΔＰを算出する（Ｓ１５）。そして、算出された圧力損失ΔＰを出力部３３が出力する（Ｓ１６）。その後、圧力損失ΔＰの推定処理が終了する。 〔実施例〕 上記実施形態の圧力損失推定システム１を用いて、各単一粒子群に対して圧縮乾燥空気の透過特性試験を行って、各単一粒子群に関する比μ’／ｋ’を計測した。図１２および図１３には、上記構成の流体透過特性試験機２を用いた実施例を示している。 図１２は、使用した単一粒子群と、各単一粒子群の平均粒径（μｍ）、空隙率ε、および単一粒子群充填層の高さＨ’（ｍｍ）とを表形式で示している。同図において、ＧＢ Ａ〜Ｃは、それぞれ平均粒径および空隙率の異なる３種類のガラスビーズである。また、ＰＰ Ａ〜Ｃは、平均粒径の異なる３種類のポリプロピレンである。 また、図１３は、図１２に示される単一粒子群のそれぞれに対して流体透過特性試験を行うことにより求められたグラフである。このグラフは、流体の空塔速度ｕＳ’（ｍ／ｓ）と、単一粒子群充填層の単位高さ当りの圧力損失ΔＰ’／Ｈ’（Ｐａ／ｍ）との対応関係を示している。ダルシーの法則（上式（１））によると、上記グラフの傾きが比μ’／ｋ’に相当する。したがって、図１３に示されるグラフから、図１２に示されるような、各単一粒子群に関する傾き（比μ’／ｋ’）（Ｐａ・ｓ／ｍ２）の測定値を取得した。この測定値を圧力損失推定装置３に出力した。 そして、各単一粒子群に関する比μ’／ｋ’の測定値を用いて流路モデルを構築し、圧力損失ΔＰを算出した。その結果、図１４に示される表と、図１５に示されるグラフとが得られた。図１４は、混合したガラスビーズと、混合したポリプロピレンとに関して、圧縮乾燥空気の体積流量Ｑを５×１０−５ｍ３／ｓで一定とし、粒子充填層の高さを種々に変えて実験を行ったものである。図示において、実験値は、流体透過特性試験機２にて複数粒子群充填層を形成して実際に計測された圧力損失ΔＰであり、計算値は、圧力損失推定装置３が算出した圧力損失ΔＰである。図１５は、これらの実験値と計算値との関係を示すものである。図１５を参照すると、本実施例で算出された圧力損失ΔＰが、実験値とほぼ一致することが理解できる。 なお、圧力損失ΔＰを精度良く求めるには、セルの寸法を小さくして、セルの数を増やすことが考えられる。そこで、セルの寸法を種々に変更して圧力損失ΔＰを推定したところ、前記セルの一辺の長さが５ｍｍ以下であれば、実験値と計算値との誤差が良好に収束した。一方、セルの数が増えると、方程式および未知数の数が急激に増加して計算処理に費やす時間が急激に増加することになる。そこで、セルの寸法を種々に変更して圧力損失ΔＰの計算に費やす時間を測定したところ、セルの一辺の長さが約３．４ｍｍであれば、圧力損失ΔＰを短時間かつ精度良く求めることができた。 また、本実施例では、３つの粒子群を混合しているが、２つの粒子群を混合してもよいし、４つ以上の粒子群を混合してもよい。 本発明は上述した実施形態に限定されるものではなく、請求項に示した範囲で種々の変更が可能である。すなわち、請求項に示した範囲で適宜変更した技術的手段を組み合わせて得られる実施形態についても本発明の技術的範囲に含まれる。 例えば、上記実施形態では、断面正方形の充填塔１１が用いられ、これにより粒子充填層２１は、断面正方形の直方体となっている。しかしながら、任意の断面形状の充填塔１１を用いることができる。例えば、断面円形の充填塔１１を用いる場合、粒子充填層２１は円柱となる。このとき、円周側壁付近で三角柱のセルを利用し、その他の領域で立方体のセルを利用してもよい。 また、圧力損失推定装置３の各構成は、ハードウェアロジックによって構成してもよいし、次のようにＣＰＵを用いてソフトウェアによって実現してもよい。 すなわち、圧力損失推定装置３は、各機能を実現する制御プログラムの命令を実行するＣＰＵ（central processing unit）、上記プログラムを格納したＲＯＭ（read only memory）、上記プログラムを展開するＲＡＭ（random access memory）、上記プログラムおよび各種データを格納するメモリ等の記憶装置（記録媒体）などを備えている。そして、本発明の目的は、上述した機能を実現するソフトウェアである圧力損失推定装置３の制御プログラムのプログラムコード（実行形式プログラム、中間コードプログラム、ソースプログラム）をコンピュータで読み取り可能に記録した記録媒体を、圧力損失推定装置３に供給し、そのコンピュータ（またはＣＰＵやＭＰＵ）が記録媒体に記録されているプログラムコードを読み出し実行することによっても、達成可能である。 上記記録媒体としては、例えば、磁気テープやカセットテープ等のテープ系、フレキシブルディスク／ハードディスク等の磁気ディスクやＣＤ−ＲＯＭ／ＭＯ／ＭＤ／ＤＶＤ／ＣＤ−Ｒ等の光ディスクを含むディスク系、ＩＣカード（メモリカードを含む）／光カード等のカード系、あるいはマスクＲＯＭ／ＥＰＲＯＭ／ＥＥＰＲＯＭ／フラッシュＲＯＭ等の半導体メモリ系などを用いることができる。 また、圧力損失推定装置３を通信ネットワークと接続可能に構成し、上記プログラムコードを通信ネットワークを介して供給してもよい。この通信ネットワークとしては、特に限定されず、例えば、インターネット、イントラネット、エキストラネット、ＬＡＮ、ＩＳＤＮ、ＶＡＮ、ＣＡＴＶ通信網、仮想専用網（virtual private network）、電話回線網、移動体通信網、衛星通信網等が利用可能である。また、通信ネットワークを構成する伝送媒体としては、特に限定されず、例えば、ＩＥＥＥ１３９４、ＵＳＢ、電力線搬送、ケーブルＴＶ回線、電話線、ＡＤＳＬ回線等の有線でも、ＩｒＤＡやリモコンのような赤外線、Ｂｌｕｅｔｏｏｔｈ（登録商標）、８０２．１１無線、ＨＤＲ、携帯電話網、衛星回線、地上波デジタル網等の無線でも利用可能である。なお、本発明は、上記プログラムコードが電子的な伝送で具現化された搬送波あるいはデータ信号列の形態でも実現され得る。 以上のように、単一粒子群を複数混合して成る複数粒子群充填層の圧力損失ΔＰを、各単一粒子群から成る単一粒子群充填層における比μ’／ｋ’と、流体の体積流量Ｑとから精度良く推定できるので、固液を分離する濾過装置、集塵のためのバグフィルタ、触媒粒子層にガスを通して反応させる触媒反応装置、および、湿った粉体を通気乾燥させる粉体乾燥装置などの設計に利用できる。本発明の一実施形態である圧力損失推定システムの処理動作を示すフローチャートである。同図（ａ）は、流体が粒子充填層を透過する状態を示す図であり、同図（ｂ）は、電流が抵抗器を流れる状態を示す回路図である。同図（ａ）は、流体が２つの粒子充填層を直列に透過する状態を示す図であり、同図（ｂ）は、流体が２つの粒子充填層を並列に透過する状態を示す図である。複数粒子群充填層をセルに細分化した状態を示す斜視図である。或るセルを流体が透過するときの流路モデルを回路記号で示す斜視図である。隣接する複数のセルの流路モデルを回路記号で示す図であり、同図（ａ）は側面図であり、同図（ｂ）は斜視図である。キルヒホッフの第一法則を説明するための簡単な例を示す図である。２次元の流路モデルを回路記号で示す図である。上記圧力損失推定システムの概略構成を示すブロック図である。上記圧力損失推定システムにおける流体透過特性試験機の概略構成を示すブロック図である。上記圧力損失推定装置の概略構成を示すブロック図である。本発明の一実施例において使用した単一粒子群と、各単一粒子群の平均粒径、空隙率、単一粒子群充填層の高さ、および傾きとを表形式で示す図である。上記実施例において、流体の空塔速度と、単一粒子群充填層の単位厚さ当りの圧力損失との関係を示すグラフである。上記実施例において、単一粒子群を複数混合して成る複数粒子群充填層における圧力損失の実験値および計算値を表形式で示す図である。上記実験値と上記計算値との関係を示すグラフである。符号の説明 ２ 流体透過特性試験機 ３ 圧力損失推定装置 ２１ 粒子充填層 ４０ 数値取得部（単一粒子群充填層特性取得手段、流体特性取得手段） ４１ セルモデル構築部（セルモデル構築手段） ４２ 抵抗決定部（抵抗算出手段） ４３ 流路モデル構築部（流路モデル構築手段） ４４ 圧力損失算出部（圧力損失算出手段） Ｃｉ セル ΔＰ 圧力損失 Ｑ 流体の体積流量 ｎｏｄ 合流点 単一粒子群を複数混合して成る粒子充填層である複数粒子群充填層に関して、該複数粒子群充填層を透過する流体の圧力損失ΔＰを推定する圧力損失推定方法であって、 前記単一粒子群は、該単一粒子群から成る粒子充填層である単一粒子群充填層が、ΔＰ’／Ｈ’＝（μ’／ｋ’）×ｕＳ’ ・・・（１）’、（但し、Ｈ’：前記単一粒子群充填層の高さ、μ’：前記単一粒子群充填層を透過する流体の粘性率、ｋ’：前記単一粒子群充填層を透過する流体の透過率、ｕＳ’：前記単一粒子群充填層を透過する流体の空塔速度）を充足する単一粒子群であって、 各単一粒子群充填層における前記比μ’／ｋ’を取得するステップと、 前記複数粒子群充填層を透過する流体の体積流量Ｑを取得するステップと、 前記複数粒子群充填層を所定寸法のセルに細分し、前記複数粒子群充填層を構成する複数の前記単一粒子群の構成比に従って、複数の前記単一粒子群に関する複数の前記単一粒子群充填層を、前記複数のセルにそれぞれ配置したセルモデルを構築するステップと、 構築されたセルモデルの各セルに関して、隣接するセルから流入する流体、または隣接するセルに流出する流体が受ける抵抗Ｒを、前記セルの寸法と、前記セルに配置された前記単一粒子群充填層の前記比μ’／ｋ’とに基づいて算出するステップと、 各セルに流体の合流点を有し、隣接する２つのセルの前記合流点間を流れる流体は、該２つのセルにそれぞれ配置された２つの前記単一粒子群充填層の前記抵抗Ｒを受ける流路モデルを構築するステップと、 構築された流路モデルに対して、前記体積流量Ｑを利用することにより、前記流体が前記流路モデルを透過することにより生じる圧力損失ΔＰを算出するステップとを含むことを特徴とする圧力損失推定方法。 前記比μ’／ｋ’を取得するステップにおいて、前記比μ’／ｋ’は、前記式（１）’を利用して、前記流体の空塔速度ｕＳ’と、単位厚さ当りの圧力損失ΔＰ’／Ｈ’とから取得することを特徴とする請求項１に記載の圧力損失推定方法。 前記セルモデルを構築するステップにおいて、前記セルモデルは、複数の前記単一粒子群充填層を、前記複数のセルにそれぞれランダムに配置したセルモデルであることを特徴とする請求項１または２に記載の圧力損失推定方法。 前記複数粒子群充填層の寸法を取得するステップをさらに含んでおり、 前記セルモデルを構築するステップにおいて、前記複数粒子群充填層の寸法に基づいて、前記複数粒子群充填層を所定寸法のセルに細分することを特徴とする請求項１ないし３の何れか１項に記載の圧力損失推定方法。 前記圧力損失ΔＰを算出するステップは、前記流路モデルに対して、前記流体の体積流量Ｑを電流値に対応付け、前記流体の合流点を接続点に対応付け、前記抵抗Ｒを電気抵抗に対応付け、前記圧力損失を電圧降下に対応付けて、キルヒホッフの法則を利用することにより、前記流路モデル全体の圧力損失を求め、これを前記複数粒子群充填層の圧力損失ΔＰとするステップであることを特徴とする請求項１ないし４の何れか１項に記載の圧力損失推定方法。 前記セルは同一寸法の立方体であることを特徴とする請求項１ないし５の何れか１項に記載の圧力損失推定方法。 前記セルの一辺の長さが５ｍｍ以下であることを特徴とする請求項６に記載の圧力損失推定方法。 前記セルの一辺の長さが約３．４ｍｍであることを特徴とする請求項７に記載の圧力損失推定方法。 単一粒子群を複数混合して成る粒子充填層である複数粒子群充填層に関して、該複数粒子群充填層を透過する流体の圧力損失ΔＰを推定する圧力損失推定装置であって、 前記単一粒子群は、該単一粒子群から成る粒子充填層である単一粒子群充填層が、ΔＰ’／Ｈ’＝（μ’／ｋ’）×ｕＳ’、（但し、Ｈ’：前記単一粒子群充填層の高さ、μ’：前記単一粒子群充填層を透過する流体の粘性率、ｋ’：前記単一粒子群充填層を透過する流体の透過率、ｕＳ’：前記単一粒子群充填層を透過する流体の空塔速度）を充足する単一粒子群であって、 各単一粒子群充填層の前記比μ’／ｋ’を取得する単一粒子群充填層特性取得手段と、 前記複数粒子群充填層を透過する流体の体積流量Ｑを取得する流体特性取得手段と、 前記複数粒子群充填層を所定寸法のセルに細分し、前記複数粒子群充填層を構成する複数の前記単一粒子群の構成比に従って、複数の前記単一粒子群に関する複数の前記単一粒子群充填層を、前記複数のセルにそれぞれ配置したセルモデルを構築するセルモデル構築手段と、 構築されたセルモデルの各セルに関して、隣接するセルから流入する流体、または隣接するセルに流出する流体が受ける抵抗Ｒを、前記セルの寸法と、前記セルに配置された前記単一粒子群充填層の前記比μ’／ｋ’とに基づいて算出する抵抗算出手段と、 各セルに流体の合流点を有し、隣接する２つのセルの前記合流点間を流れる流体は、該２つのセルにそれぞれ配置された２つの前記単一粒子群充填層の前記抵抗Ｒを受ける流路モデルを構築する流路モデル構築手段と、 構築された流路モデルに対して、前記体積流量Ｑを利用することにより、前記流体が前記流路モデルを透過することにより生じる圧力損失ΔＰを算出する圧力損失算出手段とを備えることを特徴とする圧力損失推定装置。 請求項１ないし８の何れか１項に記載の圧力損失推定方法における各ステップをコンピュータに実行させるための圧力損失推定プログラム。 請求項１０に記載の圧力損失推定プログラムが記録されたコンピュータ読取り可能な記録媒体。 【課題】 単一粒子群を複数混合してなる複数粒子群充填層に関して、各単一粒子群の粒径分布および粒子形状に関係なく圧力損失を推定する。【解決手段】 まず、複数粒子群充填層の寸法と、流体の体積流量Ｑと、各単一粒子群に関する比μ’／ｋ’とを取得する（Ｓ１１）。次に、複数粒子群充填層の寸法とセルの所定寸法とから、セルの数と配置とを決定する（Ｓ１２）。次に、各単一粒子群に関する比μ’／ｋ’と、セルの所定寸法とから、各セルの抵抗を決定する（Ｓ１３）。次に、決定したセルの数および配置と、決定した各セルの抵抗とに基づいて、流路モデルを構築する（Ｓ１４）。そして、構築された流路モデルと、流体の体積流量Ｑとに基づいて、複数粒子群充填層の圧力損失ΔＰを算出する（Ｓ１５）。【選択図】 図１