(57)【特許請求の範囲】
【請求項１】
流体の粘度を測定する方法であって、
ａ）特性横断寸法（Ｄ）のチャネル（１４、２４）内において層流条件下で前記流体の流れを確立するステップであって、特性寸法（ｄ）を有し、前記比Ｄ／ｄが５よりも大きい細長要素（１２、２２）を、前記チャネル内に、実質的にその長手方向に沿って実質的にその中心に配置し、前記細長要素の長さの部分（ｌ）を前記チャネル内に液浸する、確立するステップと、
ｂ）前記流体によって前記細長要素の壁に加わる摩擦力（ｆ）を測定するステップと、
ｃ）以下の式：
【数１】
0005474076.tif 000018
（式中：
Ｕは流れの平均速度であり、
λは幾何学的因子である）
を用いて前記流体の動的粘度（η）を測定するステップと、
を含み、
前記チャネル内で前記特性横断寸法（Ｄ）に沿って前記細長要素の様々な位置に関してステップｂ）及びｃ）を連続的に実施する、方法。
【請求項２】
前記流体がニュートン流体である場合、単一の値（Ｕ）に関してステップａ）〜ｃ）を実施する、請求項１に記載の方法。
【請求項３】
前記流体が非ニュートン流体である場合、様々な値（Ｕ）を用いてステップａ）〜ｃ）を連続的に実施する、請求項１に記載の方法。
【請求項４】
連続的に実施する、請求項１〜３のいずれか１項に記載の方法。
【請求項５】
請求項１〜４のいずれか１項に記載の測定する方法を実施する粘度計であって、
特性開口（Ｄ）及び長さ（Ｌ）のチャネル（１４、２４）と、
特性寸法（ｄ）を有すると共に、長さ（ｌ）の液浸部分を有し、前記比Ｄ／ｄが５よりも大きい細長要素（１２、２２）であって、前記チャネル内に、実質的に該チャネルの中心に、実質的に該チャネルの長手方向に沿って延びるように配置される、細長要素（１２、２２）と、
流体流が前記チャネル内で確立されると前記細長要素の壁に加わる摩擦力（ｆ）を測定する、測定手段（１３、２３）と、
前記細長要素の前記横断寸法に沿った位置に伴う前記摩擦力の変化を測定するために前記チャネルの前記開口内で前記細長要素を並進移動させる手段と、
を備える、粘度計。
【請求項６】
強制流量ポンプを含む前記流体を循環させる手段と関連する、請求項５に記載の粘度計。
【請求項７】
前記チャネル（１４）は、２つの実質的に平行なプレート（１０、１１）間の空きスペースによって画定される、請求項５又は６に記載の粘度計。
【請求項８】
前記チャネル（２４）は円筒形容器（２０）によって画定される、請求項５〜７のいずれか１項に記載の粘度計。
【請求項９】
前記チャネル及び前記細長要素は使い捨て用の要素である、請求項５〜８のいずれか１項に記載の粘度計。 
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、流体の粘度を測定する方法に関し、また、対応する装置に関する。
【背景技術】
【０００２】
研究所又は実際に工業用途において使用するのに好適な多くの粘度計が既知である。
【０００３】
流体を収容している管内に球を落下させる落球粘度計が知られており、流体の粘度を推測するためにその限界速度が測定される。
【０００４】
その粘度計は比較的正確な測定値を得ることを可能にする。それにもかかわらず、この粘度計は水等の低粘度流体の場合はあまり正確ではない。さらに、実際には、その装置は研究所でしか使用することができない。最後に、広範な粘度をカバーするために複数の球及び管を使用する必要がある。
【０００５】
圧力を加えることによって毛細管に沿って確立される流体流に基づいて動作する毛細管粘度計も存在する。流体の粘度は、流体の粘度及び圧力の両方に応じて変わる体積流量から推測することができる。
【０００６】
それらの粘度計は、粘度が広い範囲にわたって変化する可能性がある物質に使用する場合には正確さに欠ける。
【０００７】
（例えば円錐形、平面若しくは円筒形の形状の）既知の領域にかかる応力を測定するか又は既知の領域に応力を与える働きをすると共に、与えられたか若しくは測定された回転速度に基づいて２つの表面間の速度勾配（一定であると仮定される勾配）を推定する働きをする回転粘度計にも言及することができる。
【０００８】
その粘度計は、回転容器又は可動体の使用に基づく。動くことが意図される機械的部品は、正確に作製かつ位置決めする必要があり、摩擦が非常に低い回転装置の使用を要し、高価である。
【０００９】
振動粘度計は、振動棒である活動部を有する。振動粘度計は、流体の粘度の関数として変化する、振動の振幅を測定することを可能にする。
【００１０】
さらに、概して、上記粘度計は本質的に研究所において使用されるものであり、工業用途、及び／又は連続的な測定を行うにはあまり適していない。上記粘度計はまた、それらの構成要素である機械的部品に必要とされる許容差のために、原価が高い装置である。したがって、上記粘度計は、特に、該粘度計に損傷を与える可能性が高い流体に対して粘度測定を行うべきである場合に、非常に徹底した保守を必要とする器具である。さらに、上記粘度計は概して、限られた粘度範囲にわたってしか十分な精度の測定値を得ることが可能ではない。
【００１１】
また、ガスが沿って流れる管の入口及び出口間の差圧測定を用いる粘度計を記載している特許文献１にも言及することができる。
【００１２】
管内のガスの流れは、管の入口部分では層流ではないため測定が妨げられ、したがって粘度計の精度が大幅に下がる。
【００１３】
特許文献２は、同様に内部を流体が流れる管の入口及び出口間の差圧測定を用いる粘度計を記載している。
【００１４】
この粘度計の構造を考えると、流体は一方向にしか流れることができない。さらに、上記文献は、レイノルズ数が２０００を超えると述べている。管内の流体の流れは、レイノルズ数が大きい場合には層流となり得ないため、ここでもまた粘度計の精度が下がる。
【００１５】
最後に、特許文献３は、媒体中に懸濁している粒子の見掛け粘度を測定する器具を記載している。この器具は、器具内における後退位置と展開位置との間で移動可能なブレードを有する。展開位置では、ブレードはパイプ内に配置され、該パイプ内には媒体が流れ、その流れ方向は一定である。この器具はブレードに加わるせん断力を測定する。
【００１６】
それにも関わらず、流れの特性は、パイプの壁からその中心にかけて異なる。したがって、せん断力はブレードに沿って変化し、器具は媒体の粘度に関して精度の低い値を示す。この器具は本質的に、高い粘度を求めるように設計されている。
【先行技術文献】
【特許文献】
【００１７】
【特許文献１】米国特許第４７５０３５１号
【特許文献２】米国特許第６７５５０７９号
【特許文献３】国際公開第９８／３８４７７号
【発明の概要】
【発明が解決しようとする課題】
【００１８】
本発明の目的は、粘度が非常に広範な範囲にある流体、特に粘度が低い流体に対して実施することができる、流体の粘度を測定する方法を提案することによって上記の欠点を軽減することであり、その測定は、約数パーセントである許容可能な精度で行われ、本方法は、ニュートン流体にも非ニュートン流体にも等しく用いることができる。
【００１９】
本発明はまた、上記方法を実施する粘度計を提供する。本粘度計は、特にその構成部品が非常に正確な許容差を必要としないため低コストであり、研究所又は工業環境においても等しく用いることができる。
【課題を解決するための手段】
【００２０】
したがって、本発明は、流体の粘度を測定する方法であって、
ａ）特性横断寸法（characteristic transverse dimensions）Ｄのチャネル内において層流条件下で前記流体の流れを確立するステップであって、特性寸法（characteristic dimension）ｄの細長要素を、前記チャネル内に、実質的にその長手方向に沿って実質的にその中心に配置し、前記細長要素の長さの部分ｌを前記チャネル内に液浸する、確立するステップと、
ｂ）前記流体によって前記細長要素の壁に加わる摩擦力ｆを測定するステップと、
ｃ）以下の式：
【００２１】
【数１】
0005474076.tif 000002
【００２２】
（式中：
Ｕは流れの平均速度であり、
λは幾何学的因子である）
を用いて前記流体の動的粘度ηを計算するステップと、
を含む、方法を提供する。
【００２３】
本明細書を通して、「細長要素」という用語は、長さｌ’及び直径ｄを呈すると共に、長さｌを呈する液浸部分を有する要素を指すのに用いられ、比ｌ／ｄは１０〜１００の範囲にある。
【００２４】
したがって、本方法は、限られた幾何学的形状内に層流を確立することにあり、この流れはしたがって流体の粘度によって影響を及ぼされるのであって、流体の慣性によっては影響を及ぼされず、これは広範囲の流量に当てはまる。
【００２５】
さらに、本方法は、回転運動するか又は振動運動するかにかかわらず、可動部品を一切必要としない。したがって本方法の実施は必然的に、粘度を測定する従来の方法と比べて簡略化されたものとなる。
【００２６】
粘度を測定する別の技法は、比Ｄ／ｄが十分に大きく、粘度計のチャネルが円筒形である場合に通常５よりも大きければ、前記チャネル内の前記細長要素の、特性横断寸法Ｄに沿った様々な位置に対してステップｂ）及びｃ）を連続的に実施することにある。
【００２７】
これによって、摩擦力ｆ及び動的粘度を、チャネルの横断寸法に沿った細長要素の位置の関数として、したがって流れの速度プロファイルの関数として測定することが可能となる。同じ平均速度の場合、異なるせん断速度に関して粘度を求めることができる。
【００２８】
流体がニュートン流体である場合、平均流速Ｕの値を１つだけ用いてステップａ）〜ｃ）を実施することができるが、それでもなお、適切な精度を表す粘度の測定値が得られる。
【００２９】
流体が非ニュートン流体である場合、平均流速Ｕの様々な値を用いてステップａ）〜ｃ）を連続的に実施する。
【００３０】
この場合、本方法は、平均流速の関数として粘度の曲線を提供することにあり、その後この曲線を有効に分析して流体の特性を求めることができる。
【００３１】
流体が非ニュートン流体であり、比Ｄ／ｄに関して十分に大きく、粘度計のチャネルが円筒形である場合に通常は５よりも大きければ、本方法は、チャネルの横断寸法に沿った細長要素の様々な位置に関して平均流速の関数として粘度の曲線を得るように有利に実施され、流量は一定であるか、又は流量が異なる場合には細長要素の位置が固定される。
【００３２】
本発明の方法は、有利には連続的に実施することができる。
【００３３】
本発明はまた、本発明の測定方法を実施する粘度計であって、
特性開口（characteristic opening）Ｄ及び長さＬのチャネルと、
特性寸法ｄを有すると共に、長さｌの液浸部分を有する細長要素であって、前記チャネル内に、実質的に該チャネルの中心に、実質的に該チャネルの長手方向に沿って延びるように配置される、細長要素と、
流体流が前記チャネル内で確立されると前記細長要素の壁に加わる摩擦力ｆを測定する、測定手段と、
を備える、粘度計を提供する。
【００３４】
本発明の粘度計は、前記横断寸法に沿った前記要素の変化する位置に伴う前記摩擦力の変化を測定するために前記チャネルの前記開口内で前記細長要素を並進移動させる手段も含むことができる。
【００３５】
前記細長要素は可撓性とすることができる。細長要素は、可撓性であることによって、流体流の影響下で自身を前記チャネルの長手方向と容易に位置合わせすることができる。
【００３６】
例として、ポリエステル、ゴム又はナイロンのフィラメントに言及することができる。
【００３７】
粘度計は、流体を流動させる強制（imposed）流量ポンプを含むことができる。
【００３８】
前記粘度計の第１の変形形態では、前記チャネルは、２つの実質的に平行なプレート間の空きスペースによって画定される。
【００３９】
第２の変形形態では、前記チャネルは円筒形容器によって画定される。
【００４０】
本発明の粘度計の製造コストが低いことを踏まえれば、チャネル及び細長要素は使い捨てとすることができる。これによって、測定作業を行った後で洗浄作業を省くことが可能となる。これはまた、医療分野において、例えば本発明の粘度計が血液試料の粘度測定に用いられる場合に極めて有用であり得る。
【００４１】
添付の図面を参照してなされる以下の説明を読めば、本発明をよりよく理解することができ、本発明の他の目的、利点及び特徴がより明らかとなる。
【図面の簡単な説明】
【００４２】
【図１】本発明の例示的な粘度計の斜視図である。
【図２】図１に示される粘度計において用いられる細長要素が１．５ミリメートル（ｍｍ）に等しい直径のガラスシリンダーであり、流体が水である場合の、該細長要素に加わる力（ニュートン（Ｎ））を、流体の流量（ミリメートル／分（ｍＬ／ｍｉｎ））の関数としてプロットする曲線である。
【図３】流体が（２５％／７５％）水／グリセロール混合物である場合の、図１に示される装置に配置された直径２ｍｍの銅シリンダーに加わる力（Ｎ）を、流量（ｍＬ／ｍｉｎ）の関数としてプロットする曲線である。
【図４】図１の装置に配置された直径１．５ｍｍのガラスシリンダーに加わる力（Ｎ）を、そのオフセット（すなわち装置のチャネルの対称面に対するその位置）の関数として示す曲線である。
【図５】Ｄ／ｄ＝１００である場合の、図１に示されるような装置内の円筒形要素のオフセットの関数として（すなわち実際には装置のチャネルの対称面に対する上記要素の位置の関数として）幾何学的因子λを示すデジタルモデリングによって求められる曲線であり、オフセットはＤの関数として表される。
【図６】本発明の別の例示的な粘度計の概略斜視図である。
【図７】図６に示される粘度計の細長要素が１ｍｍに等しい直径の金属シリンダーであり、流体が水である場合の、該細長要素に加わる力（マイクロニュートン（μＮ））の変化を、流速（メートル／秒（ｍ／ｓ））の関数としてプロットする曲線である。
【図８】せん断速度（毎秒（ｓ−１））の関数として、図６に示される粘度計を用いて得られた水の図７の結果に対応する粘度値（ミリパスカル秒（ｍＰａ・ｓ））（◆点）と、Contraves社製Ｌｏｗ−ｓｈｅａｒ３０回転粘度計を用いて得られた値（＋点）とを示す図である。
【図９】図６に示される粘度計の細長要素が１ｍｍに等しい直径の金属シリンダーであり、流体が（１５％／８５％）水／グリセロール混合物である場合の、該細長要素に加わる力（μｍ）の変化を、流速（ｍ／ｓ）の関数として示す曲線である。
【図１０】せん断速度（ｓ−１）の関数として、図６に示される粘度計を用いて得られた（１５％／８５％）水／グリセロール混合物の図９の結果に対応する粘度値（ｍＰａ・ｓ）（●点）と、Anton Paar社製ＭＣＲ５０１レオメーターを用いて得られた粘度値（◇点）とをプロットする図であり、破線の直線は文献からの値を表す。
【図１１】Ｄ／ｄ＝５である場合、幾何学的因子λを、チャネルが管によって画定されている粘度計における円筒形要素のオフセットの関数として（すなわち管の中心に対する上記要素の位置の関数として）示す曲線であり、オフセットはＤの関数として表される。
【図１２】ｌｏｗ−ｓｈｅａｒ粘度計によって得られるような、粘度（ｍＰａ・ｓ）の変化を非ニュートン流体の場合のせん断速度（ｓ−１）の関数として示す曲線である。
【図１３】力（μＮ）の変化を流量（ｍＬ／ｍｉｎ）の関数として示す３つの曲線Ｃ１、Ｃ２及びＣ３をプロットする図であり、ここで、Ｃ１は図１２に示される結果に対応し、粘度計を流れる流体のレオロジー曲線から得られる理論上等しい力を示し、Ｃ２及びＣ３は、図６に示される粘度計を用いて同じ非ニュートン流体に関して行った試験に対応する。
【発明を実施するための形態】
【００４３】
図１に図式的に示される粘度計１は、実質的に同一である２つのプレート１０及び１１を備える。これらのプレートは、両端部以外は、実質的に滑らかであると共に平行な向かい合う表面１００及び１１０を呈する。このゾーンでは、細長要素を定位置に置くのを容易にするためにＹ字形が選択されている。このゾーンは省かれるか又は形状を変えられれば有利であろう。
【００４４】
これらの２つの向かい合う面は、流体流を通過させるのに使用されるチャネル１４を画定する。これらの２つの向かい合う面は、チャネルの特性開口又は横断寸法と称される距離Ｄだけ離間している。
【００４５】
図１に示される粘度計はまた、本明細書ではワイヤ１２の形態である細長要素を含む。図１に示されるように、細長要素１２は長さｌ’を呈し、この長さの一部ｌが、プレート１０及び１１によって画定されているチャネル１４内に配置される。当然ながら、要素１２が液体中に完全に液浸される場合、ｌ＝ｌ’である。
【００４６】
要素１２は、実質的にチャネル１４の中心に配置されるため、該要素とチャネルの壁とは全く接触しない。
【００４７】
プレートは、流体の液位から測定され得る長さＬを呈する。チャネル入口の長さＬ’がチャネル１４の長さＬに対して十分に短く、かつＬ’が要素のうち液浸される部分ｌと比較して十分に短かければ、入口の影響は無視できるほどである。図１に示される粘度計の場合、特にＬ’が０．３ｌ未満である場合に第２の条件が満たされる。
【００４８】
細長要素１２は力センサー１３に接続されている。この接続は、測定される力が細長要素の壁に加わる実際の力であるようになされる。
【００４９】
この組み立て体を、細長要素を開口内でかつ横断寸法Ｄに沿って移動させる働きをする移動テーブル（図１には示さず）上に配置する。したがって、力を、チャネルの特性開口において様々な場所に関して測定することができる。
【００５０】
細長要素は、例えばガラス若しくは金属から作製される、種々の形状、特に円筒形の形状をとることができるか、又はさらには撚り合わせたストランドの形態とすることができる。
【００５１】
原則として、管内に流体が流れる必要はない。
【００５２】
この要素の断面は、要素が円筒体である場合にはその直径に相当するか、又は撚り合わせたストランドに当てはまるようにその断面が不規則的である場合には実際にはその横断面の最大寸法に相当する、特性横断寸法ｄを呈すると考えられる。
【００５３】
図１には示されていない手段が、上記チャネル内で図１に示される矢印Ｆ１及びＦ２に沿って流体を循環させるように働き、流体の流れは層流条件下にある。流れ方向は逆にしてもよい。
【００５４】
プレートによって画定されるチャネル１４の特性寸法Ｄ及び細長要素１２の特性寸法ｄは、比Ｄ／ｄが１００以下になるように選択される。
【００５５】
この場合、図１に示される装置はいわゆる「限られた」幾何学的形状を画定する。その結果、チャネル１４内の流体の層流は本質的に、流体の慣性ではなく流体の粘度に依存する。したがって、層流条件下にある流体の流れが限られた幾何学的形状において確立される場合、流体の粘度と、該流体の流れの中に配置される細長要素に加わる粘性摩擦力との間には単純な式が存在する。この式は以下の通りである：
【００５６】
【数２】
0005474076.tif 000003
【００５７】
式中、
ｆは細長要素に加わる力であり、
Ｕは流れの平均速度（寸法Ｄにわたる速度の空間平均）であり、
ηは粘度であり、
λは、数値シミュレーション又は校正によって正確に求めることができる幾何学的因子である。
【００５８】
端部の影響が無視できるものであり、比Ｄ／ｄが５未満である場合、幾何学的因子λは、比ｄ／Ｄに伴って概ねアフィン形式で変化する。要素がまたチャネルの中心、すなわちチャネルの対称面にあり、プレート１０及び１１に平行である場合、λは（数値的及び実験的に観察される）以下の式に従う：
【００５９】
【数３】
0005474076.tif 000004
【００６０】
式（１）は、ニュートン流体、すなわち一定粘度を有する流体に有効である。同じ式を用いて非ニュートン流体を分析することができる方法を以下で説明する。
【００６１】
さらに、流れが２次元であり、かつλがｌ／Ｄとは無関係であると仮定することができるようにするには、ｌ／Ｄは十分に大きい（例えば５超である）必要がある。ニュートン流体で、かつレイノルズ数が小さい（Ｒｅが０．１未満）場合、また装置の校正を行えば、この条件は必要ではない。
【００６２】
ここで、レイノルズ数Ｒｅは以下の式によって定義されることが想起される：
【００６３】
【数４】
0005474076.tif 000005
【００６４】
式中：
ρは流体の密度（単位体積当たりの質量）であり、
ηは流体の動的粘度であり、
Ｕは平均流速であり、
Ｄは、流体が流れるチャネルの開口の特性寸法である。
【００６５】
実験的に得られた結果を詳細に調べる前に、流体流に影響を及ぼすのはどの効果であるのか、すなわち慣性効果であるのか又は粘性効果であるのかを判断するのにレイノルズ数が役立つことを想起するのが適切である。
【００６６】
レイノルズ数が２０未満である場合、流れに見られる慣性効果は無視できると仮定することができる。したがって、流れは粘性効果によって影響を及ぼされ、流れの中に配置されている細長要素に加わる力は流体の平均流速に対して線形である挙動を呈する。
【００６７】
また、本発明の方法によって、レイノルズ数の値が５０未満である場合に興味深い情報を得ることが可能となると仮定することができる。概して、慣性効果を無視し続けることができる最大レイノルズ数（したがって最大流量）を実験的に求めることが可能である。
【００６８】
ここで図２を参照すると、図２は、５ｍｍに等しい特性寸法Ｄを呈すると共に、１.５ｍｍに等しい直径ｄのガラスシリンダーが内部に配置されている、図１に示される粘度計を用いて得られた実験結果を示す。ガラスシリンダーは、プレート１０及び１１間に画定されているチャネル１４の中心に配置した。用いた流体は水であった。
【００６９】
力センサー１３によって、細長要素に加わる流体力学的力を測定した。
【００７０】
図２は、シリンダーに加わる力（ニュートンで表される）の測定点（＋）を、流速に比例する流体の流速の関数としてプロットし、これは２０．０℃の温度で蒸留水を用いて行った試験に当てはまる。
【００７１】
点線としてプロットされる曲線は、最初の測定点（最低流量）に基づく線形回帰に相当する。
【００７２】
この曲線から、流量の関数としての力の変化が実際に線形であることが確認される。直線の傾きは粘度を求めるのに役立つ。
【００７３】
図２に示される曲線に関して、測定された粘度は１．１ｍＰａ・ｓである。
【００７４】
５００ｍＬ／ｍｉｎよりも大きい流量に対応する測定点は、線形回帰に相当する点線の曲線からずれている。
【００７５】
これによって、高流量では慣性効果がなくなって無視できるほどになることが確認される。参考として、３００ｍＬ／ｍｉｎの流量に関しては、用いた粘度計内ではレイノルズ数は既に約５０である。
【００７６】
当然ながら、レイノルズ数の定義を考えると、測定装置の特性寸法Ｄが減るとこれらの慣性効果も減少し得る。
【００７７】
図３は、２ｍｍに等しい直径ｄの銅シリンダーが内部に配置されている、図１に示される粘度計を用いて得られた実験結果を示す。銅シリンダーは、実質的にチャネル１４の中心に配置し、この場合、装置の特性寸法Ｄは５ｍｍに等しい。用いた流体は水とグリセロール（ニュートン流体）との混合物であり、グリセロール画分は約７５重量％であり、温度は２１．９℃であった。
【００７８】
測定点は十字（＋）で表し、黒い点線は曲線に対する直線の当てはめ（linear fit）に相当する。
【００７９】
ここでもまた、力が流量と共に線形変化することが分かり、直線の傾きは粘度を求めるのに役立つ。測定された粘度値は３５．２ｍＰａ・ｓである。
【００８０】
図３は、上記実験において、高流量であっても測定点が線形回帰曲線からずれなかったことを示す。
【００８１】
実際には、上記実験において、レイノルズ数は常に５未満であった。
【００８２】
図２及び図３に示される結果は、用いられるニュートン流体（例えば水又は普通の油）とは関係なく、管に加わる力の線形変化が流量の変化と共に実際に観察されることを示す。
【００８３】
したがって、そのような流体に関して、また図２及び図３に示されるように、流体の粘度値を得るには所与の流量における１回の測定で十分である。測定される力の値は、慣性効果と比較して粘性効果の方が影響を及ぼす限り、流量に比例する。
【００８４】
式（１）がニュートン流体に関して線形であることは、流量又は流速が経時にわたって一定ではなくても式（１）が当てはまることを意味する。その場合、式（１）の速度及び力をそれらの経時にわたる平均値と置き換えればよい。
【００８５】
したがって、本発明の装置によって、装置を通る流量が一定ではなくてもニュートン流体の粘度を測定することが可能となる。これは特に、パルス流、又は実際には蠕動ポンプを用いて生成される流れの場合に起こり得る。
【００８６】
ここで図４を参照すると、図４は、図１に示される装置のチャネル内の細長要素の位置の影響を示すために得られた他の実験結果に関する。
【００８７】
したがって、図４は、５ｍｍに等しい特性寸法Ｄを有する装置１のチャネル１４内に配置された、直径ｄが１．５ｍｍのガラス管に加わる力の測定点（グレーのバー（Ｉ）によって表される）を示し、ガラス管の位置は、チャネルの対称面、すなわち両方のプレート１０及び１１から等距離に位置付けられる面に対して測定される。管とチャネルの中心との間の距離、すなわち実際には管のオフセットは、−１．７５ｍｍ〜＋１．７５ｍｍの範囲にわたって変化し得る。
【００８８】
測定は、水とグリセロールとの混合物から構成される流体を用いて行い、グリセロールの濃度は７５．７重量％であり、温度は２１．４℃であった。
【００８９】
図４に示されるグレーの垂直バーは、行った測定の最大分散を示し、一方でグレーの曲線は実験による測定値に対応する。
【００９０】
図４には、数値シミュレーションによって計算した流れの粘度値に相当する黒い破線曲線も見られる。Ｄ、ｄ及び細長要素のオフセットから開始して、（数値計算によって）オフセットの関数としてλを計算することが可能である。その後、式（１）を用いる。２つの曲線は、３％以内で一致しているため、実質的に重なっている。図１に示される例示的な粘度計では、プレート又はチャネルの幅Ｌａは９センチメートル（ｃｍ）であり、幅の影響はデジタル計算では考慮されない。これは、幅Ｌａがチャネルの寸法Ｄに対して大きい限り、測定値に影響を及ぼすのは主に寸法Ｄである、すなわち、その場合、寸法Ｌａの影響は無視できることを示す。
【００９１】
さらに、１．６ｍｍの範囲（−０．８ｍｍ〜＋０．８ｍｍ）にわたって、測定される力は、チャネルの中心における値から５％未満だけずれる。これは、幾何学的因子λが、管がチャネルの中心の近傍にあるときにチャネル内の管の位置にはほとんど左右されず、また、チャネルの壁及び管の壁に存在し得るいかなる凹凸（irregularities）にもほとんど左右されないことを示す。
【００９２】
図５は、比Ｄ／ｄ＝１００である場合の、オフセットの関数としての幾何学的因子λの変化を示す。幾何学的因子λは、円筒形要素のＤ正規化オフセット（すなわち、要素の長手方向軸と、図１において符号１００及び１１０が与えられている壁の対称面との間の距離）の関数として、有限要素モデルを用いた数値モデリングによって得られる。この図では、このオフセットは、２つの壁間の距離Ｄを用いることにより無次元にされる。オフセットは、細長要素が壁間の中間にある場合はゼロであり、細長要素が壁の一方に触れている場合（その場合、該要素は、その全長にわたって壁と接触している）は最大となる。この例では、最大オフセットは以下の通りである：
【００９３】
【数５】
0005474076.tif 000006
【００９４】
対称性の理由から、また図４に示されるように、幾何学的因子（又は摩擦力）は、オフセットの正負符号に左右されることはない。したがって、正のオフセット値のみを図５に示す。
【００９５】
これらの結果は、オフセットが±０．１Ｄ（すなわち±１０ｄ）未満である限り、幾何学的因子がほとんど変化しない（中心におけるその値に対して５％未満である）ことをはっきりと示す。このことは、細長要素のセンタリングに関する不正確さにもかかわらず本方法のロバスト性を示す。測定は、位置決めが少し不正確でも信頼性の高い結果を与える。
【００９６】
比Ｄ／ｄが５よりも大きく、増加してゆく場合、測定は、チャネルの開口内の速度プロファイルに対してますます敏感になる。開口内の移動を用いる力測定における変化を用いてせん断速度の関数として粘度の変化を得ることができることが有利である。
【００９７】
これはいくつかの結果につながる。
【００９８】
まず、本発明の方法を、図１に示される種類の装置を用いて実施することによって、細長要素がチャネルの中心に正確に置かれることを必須とすることなく、正しい測定値を得ることが可能となる。
【００９９】
さらに、細長要素の表面は必ずしも厳密に規則的でなくてもよい。
【０１００】
したがって、実験により、１本のストランドを共に撚り合わせた複数のストランドから構成される細長要素を用いても有効な測定を行うことが可能であることが示されている。そのような状況下では、細長要素の断面は円形ではない。それにもかかわらず、要素の平均直径を用いることによって因子λを正確に計算することができる。
【０１０１】
概して、細長要素の表面状態及びチャネルの壁の表面状態は非常に正確な許容差は必要としない。これは当然ながら、本発明の粘度計の製造コストに関して有利な影響を有する。
【０１０２】
さらに、細長要素が作製される材料が問題にならないことを示す実験も行った。図２及び図３に示される結果では、細長要素はガラス又は銅から作製される。縫糸等の撚り合わせたポリエステルストランドから構成される細長要素を用いることを想定することも可能である。
【０１０３】
このことはまた、本発明の粘度計を製造するコストを低減することに寄与する。
【０１０４】
以下の表（表１）を参照する。この表は、本発明の方法及び粘度計によって得ることができる精度を示す。
【０１０５】
【表１】
0005474076.tif 000007
【０１０６】
用いた流体それぞれに関して、表１は、本発明の方法及び粘度計を使用することによって測定された粘度を、図２〜図４において見られる実験結果に従って示す。この表はまた、幾つかの表から得られる粘度値、及び用いた流体の密度及び温度の測定値を示す。参考として、そのような幾つかの表は以下のアドレスで利用可能である。
http://www.dow.com/glycerine/resources/table18.htm
これらの値の比較から、測定された粘度値が基準値に非常に近いことが示される。さらに、これらの値の不確かさは検査中の粘度の範囲にわたって変化しない。上記実験から生じる観察結果を確認する他の実験を行った。
【０１０７】
さらに、非常に広範囲にわたって力を測定することを可能にする力センサーが市場で入手可能である。例えば、１０−５グラム（ｇ）〜８０ｇの範囲内にある見掛け重量に相当する１０−７Ｎ〜８×１０−１Ｎの範囲にわたる力の測定を可能にするSartorius社製ＣＰ  ２２５Ｄ天秤に言及することができる。
【０１０８】
これは、粘度値を実質的に一定である精度で広い範囲にわたって測定することを可能にすることに寄与する。
【０１０９】
したがって、振動粘度計のような既知の粘度計と比較して、本発明の粘度計は、水の粘度のような非常に低い粘度を測定することが可能である。
【０１１０】
上記説明から、本発明の粘度計は、回転移動するのであろうと振動移動するのであろうと機械的可動部品を全く有しないことになる。また、このことは、静止部品は回転運動又は振動運動する部品よりも低い精度で機械加工することができるという点で粘度計の製造コストに関して影響も有する。さらに、それによって破損のいかなるリスクも大幅に制限される。
【０１１１】
図１に示される例では、粘度計は、チャネルを画定する２つの向かい合うプレートから作製される。チャネルが管によって画定される粘度計を作製することも想定可能である。
【０１１２】
したがって、図６に示される粘度計２は、流体流を通過させるチャネル２４を画定する円筒形容器２０を備える。チャネルの直径Ｄはチャネルの特性寸法に相当する。
【０１１３】
粘度計２も、この場合は直径又は特性寸法ｄのシリンダー２２によって構成される細長要素を有する。細長要素は、チャネル２４内に液浸される部分ｌを有する長さｌ’を呈する。チャネルの長さＬは常にｌよりも長い。
【０１１４】
細長要素は、実質的にチャネル２４の中心に配置されるため、チャネルの壁とは全く接触しない。
【０１１５】
要素２２は、この場合は精密天秤である力センサー２３に接続されており、そのため測定された力は要素の壁に加わる実際の力である。
【０１１６】
手段（図示せず）が、流体を、チャネル２４内で注入流体に対応する矢印Ｆ３又は吸い上げ流体に対応する矢印Ｆ４に沿って流す働きをする。いずれにしても、流体流は層流条件下で生じる。
【０１１７】
図１を参照して上述したように、比Ｄ／ｄは１００以下である。
【０１１８】
例として、チャネルの直径Ｄは１ｃｍであり、シリンダーの直径ｄは１ｍｍであり、シリンダーの長さｌ’は１１ｃｍであり、液浸部分ｌは７ｃｍである。
【０１１９】
上記で説明したように、流れの慣性効果が無視できる（Ｒｅが５０未満、又はさらにはＲｅが２０未満）限り、センサー２３によって測定される力ｆは式（１）によって定義され、したがって、幾何学的因子λが求められれば、流体の粘度を求めることが可能となる。
【０１２０】
このように求めることは、粘度が既知である流体を用いる校正によって行うことができるか、又は数値的にも行うことができる。
【０１２１】
数値計算の場合、ストークスの法則を解くことによって速度場が計算される（この場合、慣性効果は無視できる）：
【０１２２】
【数６】
0005474076.tif 000008
【０１２３】
流れの速度は、長手方向軸から測定した場合、半径方向距離ｒにのみ伴って変化する：
【０１２４】
【数７】
0005474076.tif 000009
【０１２５】
ニュートン流体の場合、粘度は一定の値を有するため、Ｖｚ（ｒ）を得るために解析的解法を用いることができる。
【０１２６】
非ニュートン流体の場合、速度Ｖｚ（ｒ）を得るために、この２Ｄ表現はＦｒｅｅｆｅｍ＋＋ソフトウェアを用いて繰り返し解かれる。この場合、粘度はせん断速度に依存する。この新たな制限事項は、繰り返し毎に粘度場を推定する必要があることを意味する。
【０１２７】
細長要素に作用する力は、２つの寄与因子（contributions）、すなわち圧力及び摩擦力の和である。力ｆの一般式は以下の形を有する：
【０１２８】
【数８】
0005474076.tif 000010
【０１２９】
圧力
摩擦力
積分はシリンダーの外形に関して計算される。
【０１３０】
流体の粘度が既知であり、速度Ｕが強制され、既知の長さｌの場合、力ｆを計算するだけで十分であり、次いで幾何学的因子λが比ｆ／ηＵｌを計算することによって得られる。
【０１３１】
ここで図７を参照すると、図７は、１ｃｍに等しい特性寸法Ｄを呈すると共に１ｍｍに等しい直径ｄを有する金属シリンダーが内部に配置されている、図６に示される粘度計を使用することによって得られた実験結果を示す。シリンダーはチャネル２４の中心に配置した。用いた流体は水であった。
【０１３２】
そのような条件下では、幾何学的因子は、上記の数値方法を適用すると９．９であると推定された。
【０１３３】
２つの異なる流体温度（２２．６℃及び２２．８℃）に関して、流れ方向を変えることによって（流体Ｆ３の注入：正の速度；流体Ｆ４の吸い上げによって：負の速度）複数回の試験を行った。
【０１３４】
天秤２３によって測定した重量は、要素２２に加わる流体力学的力を求めるのに役立つ。
【０１３５】
図７は、水の速度（ｍ・ｓ−１）の関数としてμＮで表される力の測定点（２２．６℃の温度に関しては＋、２２．８℃の温度に関しては□＋）を示す。
【０１３６】
この曲線（２２．８℃の温度に関しては連続線、２２．６℃の温度に関しては破線）は線形回帰に相当する。
【０１３７】
ここでもまた、これらの曲線から、速度の関数としての力の変化が実際に線形であることが確認される。直線の傾き及び式（１）は粘度の測定を可能にする。
【０１３８】
したがって、測定された粘度は、２２．６℃では０．９８２ｍＰａ・ｓであり、２２．８℃では０．９２３８ｍＰａ・ｓである。これらの測定値を、対応の表に記載されている値（２２．６℃では０．９４２ｍＰａ・ｓ、２２．８℃では０．９３７ｍＰａ・ｓ）と比較することができる。
【０１３９】
本発明の粘度計を使用して得られた値の精度をよりよく評価するために、Contraves社製Ｌｏｗ−ｓｈｅａｒ  ３０回転粘度計を用いて２２．６℃の温度で比較測定を行った。チャネル２４内のせん断速度は以下の式を用いて求めた：
【０１４０】
【数９】
0005474076.tif 000011
【０１４１】
式中：
【０１４２】
【数１０】
0005474076.tif 000012
【０１４３】
：せん断速度；
Ｑ：セルに加わる体積流量；
Ｒ：セル半径。
【０１４４】
図８は、せん断速度の関数として、水の粘度に関して本発明の粘度計を用いて得られた値（◆）と、Contraves社製Ｌｏｗ−ｓｈｅａｒ  ３０粘度計を用いて得られた値（＋）とを示す。
【０１４５】
連続的な曲線は、Contraves社製Ｌｏｗ−ｓｈｅａｒ  ３０粘度計に対応する点を繋げたものである。
【０１４６】
さらに、直線の破線は、２２．６℃の温度で０．９４２ｍＰａ・ｓである水の粘度に関して表に記載されている値を示す。
【０１４７】
図８は、水の粘度に関する測定誤差が、本発明の粘度計を用いた場合には、Ｌｏｗ−ｓｈｅａｒ粘度計（この粘度計は高価であり使用が比較的複雑であるにもかかわらず）を用いて得られた測定値と比較して小さい（１０％以下）ことを示す。この測定誤差は、温度をよりよく制御することによって低減することができる。
【０１４８】
図６に示される粘度計を用いて、図７において設定された条件下で他の試験を行った。用いた流体は８５重量％のグリセロールと水との混合物であった。
【０１４９】
測定温度（２３℃）における上記混合物の表に記載されている粘度は９８．４ｍＰａ・ｓである。
【０１５０】
上記混合物によってシリンダーにおいて得られた力を、図９に詳細に示す。
【０１５１】
ここでもまた、図９は、流れの速度の関数としての力の変化が実際に線形であることを示す。この直線の傾きによって粘度を求めることが可能であり、粘度は９６．９ｍＰａ・ｓであると分かる。
【０１５２】
粘度に関して表に記載されている値（９８．４ｍＰａ・ｓ）との比較によって、本発明の粘度計を用いて測定された粘度値が基準値と非常に近いことが示される。
【０１５３】
本発明の粘度計を用いて得られた値の精度を力の値に基づいてよりよく評価するために、式（１）を用いて動的粘度を計算した。これらの結果を、表に記載されている値、及びAnton Paar社製ＭＣＲ５０１レオメーターを用いて行った測定値と比較した。
【０１５４】
図１０は、検討中の混合物の粘度に関して、せん断の関数として本発明の粘度計を用いて得られた値（●点）と、ＭＣＲ５０１レオメーターを用いて得られた値（◇点）とを示す。
【０１５５】
連続線の曲線はＭＣＲ５０１レオメーターに対応する点を繋げたものである。
【０１５６】
さらに、直線の破線は、混合物の粘度に関して表に記載されている値（９８．４ｍＰａ・ｓ）を示す。
【０１５７】
図１０は、本発明の粘度計によって、表に記載されている粘度に近い値（平均で３％以内、最悪でも７％以内）を得ることが可能となることを示す。Anton Paar社製レオメーターを用いて行った測定値に関しては、これらの測定値は表に記載されている値よりも５％低かった。
【０１５８】
この比較試験によって、本発明の粘度計は、Anton Paar社製レオメーターの精度よりも優れた精度で流体の粘度を求めるのに役立つことが示される。この観察結果は、粘度が低い流体に関してはさらにより顕著であるはずである。
【０１５９】
本発明の方法の重要な結果は、本発明によって粘度計のサイズをかなり小さくすることが可能となることである。これによって、測定された力の値の精度に対して何ら影響することなく、小さい体積の流体に関して測定を行うことが可能となる。
【０１６０】
次元解析によって（かつ実験的に検証したところ）、幾何学的因子λは、先験的にｌ／Ｄ、Ｄ／ｄ、及びオフセットとＤとの比に依存することを示すことができる。実際には、比ｌ／Ｄが非常に大きく、例えば５０よりも大きい場合、幾何学的因子λはそれに依存することはほとんどない。細長要素が（図１に示されるような）２つの平行なプレート間に配置される場合、細長要素がプレートの中心近くに配置されると、オフセットによるλの変化は特に小さい。この点が図４に示されており、図４は、Ｄ／ｄ＝３の場合、λはオフセットの関数として連続的に変化し、中心で最大となり、変化は小さい（５％未満）ことを示す。
【０１６１】
この点は、Ｄ／ｄ＝１００の場合の、幾何学的因子λの変化をオフセットの関数として示す図５においてさらに実証される。この因子は、−０．１Ｄから０．１Ｄに及ぶ広範な範囲にわたって５％未満だけ変化する。この幾何学的形状に関して、オフセットは測定値にはほとんど影響せず、λはオフセットに関係なく主に比Ｄ／ｄの関数である。その結果、物体の位置決めが不十分であっても測定値にはほとんど影響しない。
【０１６２】
中心の近傍において測定される力が細長要素を位置決めする精度にほとんど左右されないことが望ましい場合、シリンダーを管内に配置するときに比Ｄ／ｄが５よりも大きいのが適切である。
【０１６３】
この状況が図１１に示されており、図１１は、Ｄ／ｄ＝５の場合のλの変化をオフセットの関数として示す。
【０１６４】
対称性の理由から、また図５におけるように、正のオフセット値のみを図１１に示す。幾何学的因子λは、有限要素法を用いた数値モデリングによって、円筒形要素のＤ正規化オフセットの関数として得られる。ここで、要素は直径Ｄを有する管又はチャネル内に配置され、２つの軸は平行であり、直径の比はＤ／ｄ＝５に等しく選択される。この状況では、オフセットは、円筒形要素の長手方向軸と管の対称軸との間の距離に相当する。したがって、オフセットは、２つの軸が一致する場合はゼロであり、細長要素が管の内壁に触れる場合にその最大値に達する。そのような状況下では、直径Ｄに対して正規化された最大オフセットは：
【０１６５】
【数１１】
0005474076.tif 000013
【０１６６】
である。
【０１６７】
よって、幾何学的因子、及びしたがって（これに比例する）力は、細長要素のオフセットが０．１Ｄである場合には１０％未満だけ変化する。
【０１６８】
チャネルが平行なプレートによって画定されている、図１に示されるタイプの粘度計の場合、Ｄ／ｄが５よりも大きくなる（この状況では開口内の力の変化が５％よりも大きくなる）と、開口内の力の変化が相当なものとなる。この場合、チャネルの開口内の、その横断寸法に沿った要素の位置の関数として力の変化を測定することが、開口にわたって変化するせん断速度の関数として流体のレオロジー特性を決定する上で有利となり得る。
【０１６９】
細長要素がセンタリングされている場合、重要な結果は、平均流速Ｕ、細長要素の液浸部分の長さｌ及び特性寸法の比ｄ／Ｄが全て同じままであれば、また比ｌ／Ｄが５０よりも大きければ、装置を小型化することができることである。そのような小型化は測定される力の値に実際に影響することはない。
【０１７０】
例えば、長さｌが２ｃｍに等しく、所与の速度Ｕの場合、測定される力ｆは、直径Ｄ＝２ｍｍの管内の直径ｄ＝１ｍｍのシリンダーの場合と、直径Ｄ＝２００マイクロメートル（μｍ）の管内に配置される直径ｄ＝１００μｍのシリンダーの場合とで事実上同じである。反対に、測定を行うために用いられる流体の体積は１００分の１になる。
【０１７１】
この小型化は、本発明の粘度計が可動部品を全く有しないことによっても可能となる。
【０１７２】
この小型化は、本発明の粘度計を作製するコストを低減することにさらに寄与する。
【０１７３】
さらに、本発明の粘度計を、流体のチャネル及び細長要素（いずれも使い捨て専用である）を画定する部品によって作製することを想定することができる。
【０１７４】
これは、チャネル及びその内部に配置される細長要素を形成する部品に損傷を与える可能性が高い流体において測定を行うことになる場合に特に有利であり得る。この利点はまた、多くの場合に使い捨て装置を使用する必要がある医療分野において極めて有用であり得る。
【０１７５】
上記説明は本質的にニュートン流体に関する。それにもかかわらず、本発明の方法は非ニュートン流体にも適用可能である。
【０１７６】
非ニュートン流体の場合、粘度計の細長要素に加わる力は、流れの速度に伴って非線形に変化し、流体の粘度はせん断速度に依存する。
【０１７７】
非ニュートン流体に関して、本方法は、流れの速度を変更することなくチャネルの開口内で細長要素を動かしながら力を測定することによって、又は細長要素の位置を静止させたまま異なる流量で力を測定することによって、流れの速度に関して様々な値において連続的に適用される。
【０１７８】
関与する速度は空間平均され、時間内で先験的に瞬間的である。いずれにしても、その後、結果を使用しやすくなるのであれば、時間平均をとることが可能である。この場合、粘度計は、測定された力の変化を流量の関数として、又は開口内の位置の関数として示す曲線を提供するのに役立つ。
【０１７９】
次いで、得られる情報を、既知の数学モデルを用いて従来のやり方で処理することができる。
【０１８０】
中性多糖であるスクレログルカン２５０ｐｐｍ（百万分率）を含有する水溶液を用いて試験を行った。
【０１８１】
図１２に示されるようなＬｏｗ−ｓｈｅａｒ粘度計を用いて得られたレオロジー曲線は、流体の流動性（rheofluidizing nature）をはっきりと示す。
【０１８２】
測定点（●）は、カロー（Carreau）型流体に関しては以下の式を用いて当てはめることができる：
【０１８３】
【数１２】
0005474076.tif 000014
【０１８４】
式中、η、η０及びη∞はそれぞれ、動的粘度、ニュートンプラトーにおける粘度、及び無限であるせん断速度
【０１８５】
【数１３】
0005474076.tif 000015
【０１８６】
における粘度を表す。
【０１８７】
【数１４】
0005474076.tif 000016
【０１８８】
は、２つの流れ様式（ニュートンプラトー及びべき乗則レジーム（power law regime））間の遷移に対応するせん断速度である。Ｌｏｗ−ｓｈｅａｒ粘度計によって得られるこれらの係数の値を表２に示す。
【０１８９】
【表２】
0005474076.tif 000017
【０１９０】
続いて、式（６）によって示されるレオロジー方程式をＦｒｅｅｆｅｍ＋＋ソフトウェアに再び入力して、シリンダーの内部に対して加わる摩擦力の数値推定値を得た（式（５）を参照されたい）。
【０１９１】
さらに、１ｍｍに等しい直径ｄ及び７ｃｍの液浸部分ｌを有する細長要素を有する、図６に示される粘度計を用いて、同じ温度で同じ流体を用いて２つの試験を行った。
【０１９２】
図１３は、本発明の粘度計を用いて行った試験に対応する曲線Ｃ２（−＋−）及びＣ３（−▼−）と共に、図１２に示されるレオロジー曲線から得られた曲線Ｃ１（−●−）をプロットする。
【０１９３】
それらの３つの曲線間の比較によって、本発明の粘度計を用いて得られた結果がＦｒｅｅｆｅｍ＋＋を用いて数値計算された力と良好に相関することが示される。数値計算による値と実験値との最大の差は７％である。
【０１９４】
したがって、本発明の粘度計は、例えば１ｍＰａ・ｓ未満から１００ｍＰａ・ｓの範囲にある粘度でニュートンであるか又は非ニュートンである流体に対して測定を行うことを可能にする。このシステムによって、正確かつ迅速な（数秒を要する）測定を連続的に行うことが容易となる。
【０１９５】
さらに、本発明の方法は、測定チャネル又はセルのサイズに関係なく使用可能であるという利点を有する。すなわち、力は流れの速度に比例する（流れはチャネル又は容器の断面に反比例する）。したがって、粘度計を小型化することによって、所与の流量で、試験に供した粘度計と少なくとも同程度の力を得ることが可能である。この場合、小型化された装置において用いられる流体の量は、ミリメートル（ｍＬ）オーダーであり得る。
【０１９６】
小型化された粘度計の用途は、大量に得るには高価であるか又は大量に得るのが難しい流体、場合によっては非ニュートン流体の測定であり得る。
【０１９７】
本発明の方法及び粘度計は本質的に、液体の粘度の測定を意図する。
【０１９８】
本発明の方法は、別の利点、すなわち、測定を連続的にかつ非常に短い時間間隔で行うことができ、したがって粘度の変化を測定することができるという利点を呈する。この利点は、或る特定の用途において（特にプロセスエンジニアリングにおいて）必要不可欠である。反対に、毛細管粘度計、球粘度計、又はカップ粘度計を用いる場合、連続的な測定は非常に困難である。
【０１９９】
行った実験において、応答時間は約０．２秒である。制御／記録装置を有するポンプ及び天秤の接続（interfacing）によって応答時間が本質的に制限されることが分かっている。それにもかかわらず、測定値を取得する速度を増加させることが完全に可能であり、また、チャネルの開口内での細長要素の移動を制御する働きをする移動テーブルと接続することが完全に可能である。
【０２００】
概して、或る特定の流体の粘度は経時にわたって変化し得ることが知られている。本発明の文脈では、経時にわたるこの変化は、流体が流れる瞬間から測定することができる。
【図１】
0005474076.tif 000019
【図２】
0005474076.tif 000020
【図３】
0005474076.tif 000021
【図４】
0005474076.tif 000022
【図５】
0005474076.tif 000023
【図６】
0005474076.tif 000024
【図７】
0005474076.tif 000025
【図８】
0005474076.tif 000026
【図９】
0005474076.tif 000027
【図１０】
0005474076.tif 000028
【図１１】
0005474076.tif 000029
【図１２】
0005474076.tif 000030
【図１３】
0005474076.tif 000031