水の粘度 (mpa s) y R φd p τ 図 3. 円管内流速分布の解析モデル 応力を τ とすると 円筒の全側面に作用するせん 断応力による力は となる 定常状態 では この力と圧力による力 () 式が釣り合うので これから τ を求めると () (3) 壁から円筒側面までの距離を y とす

Similar documents
伝熱学課題

第 3 章二相流の圧力損失

Microsoft PowerPoint - Š’Š¬“H−w†i…„…C…m…‰…Y’fl†j.ppt

2 図微小要素の流体の流入出 方向の断面の流体の流入出の収支断面 Ⅰ から微小要素に流入出する流体の流量 Q 断面 Ⅰ は 以下のように定式化できる Q 断面 Ⅰ 流量 密度 流速 断面 Ⅰ の面積 微小要素の断面 Ⅰ から だけ移動した断面 Ⅱ を流入出する流体の流量 Q 断面 Ⅱ は以下のように

Microsoft PowerPoint - 第7章(自然対流熱伝達 )_H27.ppt [互換モード]

Microsoft PowerPoint - 12_2019裖置工�榇諌

<4D F736F F F696E74202D C CC89C88A B8CDD8AB B83685D>

ギリシャ文字の読み方を教えてください

第1章 単 位

untitled

7 章問題解答 7-1 予習 1. 長方形断面であるため, 断面積 A と潤辺 S は, 水深 h, 水路幅 B を用い以下で表される A = Bh, S = B + 2h 径深 R の算定式に代入すると以下のようになる A Bh h R = = = S B + 2 h 1+ 2( h B) 分母の

Microsoft PowerPoint - zairiki_3

20年度一次基礎略解

伝熱学課題

Taro-水理計算.$td

横浜市環境科学研究所

Microsoft PowerPoint - 発表II-3原稿r02.ppt [互換モード]

1

Microsoft PowerPoint - 流体力学の基礎02(OpenFOAM 勉強会 for geginner).pptx


技術資料 JARI Research Journal OpenFOAM を用いた沿道大気質モデルの開発 Development of a Roadside Air Quality Model with OpenFOAM 木村真 *1 Shin KIMURA 伊藤晃佳 *2 Akiy

Microsoft Word 卒業論文2.doc

フィードバック ~ 様々な電子回路の性質 ~ 実験 (1) 目的実験 (1) では 非反転増幅器の増幅率や位相差が 回路を構成する抵抗値や入力信号の周波数によってどのように変わるのかを調べる 実験方法 図 1 のような自由振動回路を組み オペアンプの + 入力端子を接地したときの出力電圧 が 0 と

問題 2-1 ボルト締結体の設計 (1-1) 摩擦係数の推定図 1-1 に示すボルト締結体にて, 六角穴付きボルト (M12) の締付けトルクとボルト軸力を測定した ボルトを含め材質はすべて SUS304 かそれをベースとしたオーステナイト系ステンレス鋼である 測定時, ナットと下締結体は固着させた

<4D F736F F D208D5C91A297CD8A7793FC96E591E631318FCD2E646F63>

構造力学Ⅰ第12回

Microsoft PowerPoint EM2_15.ppt

パソコンシミュレータの現状

<4D F736F F D20824F B CC92E8979D814696CA90CF95AA82C691CC90CF95AA2E646F63>

物性物理学 I( 平山 ) 補足資料 No.6 ( 量子ポイントコンタクト ) 右図のように 2つ物質が非常に小さな接点を介して接触している状況を考えましょう 物質中の電子の平均自由行程に比べて 接点のサイズが非常に小さな場合 この接点を量子ポイントコンタクトと呼ぶことがあります この系で左右の2つ

T ダイの流動解析 HASL 社 FlatCAD を使用した池貝製 T ダイの流動解析事例 各種の樹脂粘度を考慮した T ダイの流路設計 Rich Green on Land Deep Blue in Sky and Sea 株式会社池貝開発室横田新一郎

粒子画像流速測定法を用いた室内流速測定法に関する研究

Microsoft PowerPoint - ‚æ3‘Í [„Ý−·…‡†[…h]

<4D F736F F D2097CD8A7793FC96E582BD82ED82DD8A E6318FCD2E646F63>

Microsoft PowerPoint - 熱力学Ⅱ2FreeEnergy2012HP.ppt [互換モード]

Microsoft PowerPoint - 夏の学校(CFD).pptx

<4D F736F F D208D5C91A297CD8A7793FC96E591E631308FCD2E646F63>

Microsoft PowerPoint - siryo7

19年度一次基礎科目計算問題略解

座標軸以外の直線のまわりの回転体の体積 ( バウムクーヘン分割公式 ) の問題の解答 立体の体積の求め方 図 1 の立体の体積 V を求める方法を考えてみる 図 1 図 1 のように 軸の から までの長さを 等分する そして とおく とすると となる 図 1 のように のときの 軸に垂直な平面 に

PowerPoint プレゼンテーション

<4D F736F F D208D7E959A82A882E682D18F498BC78BC882B B BE98C60816A2E646F63>

第 2 章 構造解析 8


Microsoft Word - NumericalComputation.docx

微分方程式による現象記述と解きかた

Microsoft PowerPoint - 混相流学会2011_髙木.ppt

Brookfield回転粘度計

(162) 岡 山本 渡辺まず Jeffery(2) と Starkey(3) は粘性流体中に一個の剛体粒子があるとき 流動による散逸エネルギーを極少にするようにふるまうと考え 軸集中の説明をこころみている 山本と大木は同様の原理を仮定し 定常状態にある懸濁液内での粒子の分布を統計的に論じた 一方

実験題吊  「加速度センサーを作ってみよう《

PowerPoint プレゼンテーション

Microsoft Word - thesis.doc

第1章 単 位

Microsoft PowerPoint - シミュレーション工学-2010-第1回.ppt

Microsoft Word - 1B2011.doc

Microsoft PowerPoint - 第3回MSBS研究会.pptx

Microsoft Word - Chap17

円筒歯車の最適歯面修整の設計例 適正な歯面修整で負荷容量の増大を目指すー 目次 1. はじめに 2 2. ` 解析例 ( はすば歯車の例 ) 歯車諸元 (C 面取り ) 歯車諸元 (R 面取り ) 最適歯面修整 歯先修整 + 歯先

破壊の予測

Microsoft PowerPoint - 第5章(対流熱伝達)講義用_H27.ppt [互換モード]

研究成果報告書

A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical

例 e 指数関数的に減衰する信号を h( a < + a a すると, それらのラプラス変換は, H ( ) { e } e インパルス応答が h( a < ( ただし a >, U( ) { } となるシステムにステップ信号 ( y( のラプラス変換 Y () は, Y ( ) H ( ) X (

Kumamoto University Center for Multimedia and Information Technologies Lab. 熊本大学アプリケーション実験 ~ 実環境における無線 LAN 受信電波強度を用いた位置推定手法の検討 ~ InKIAI 宮崎県美郷

Microsoft Word - 中村工大連携教材(最終 ).doc

PowerPoint プレゼンテーション

第6章 実験モード解析

Brookfield回転粘度計

Microsoft Word - 資料9 水理計算参考

Microsoft PowerPoint - H21生物計算化学2.ppt

コンクリート工学年次論文集 Vol.32

Microsoft Word - 第5章.doc

Microsoft PowerPoint - 基礎電気理論 07回目 11月30日

Microsoft PowerPoint - elast.ppt [互換モード]

品 名 ホロコーンパターンノズル ホロコーンノズル 単孔式 KSC 多孔式 KSC ー, ホロコーンアトマイジングノズル QC ノズル 単孔式 多孔式 型 KSN 式 KSWC ー QC ー T, KSWC ー QC KSFC ー, KSWC ー QC ー EE C. C.5 C.7 ホロコーンパタ

<4D F736F F F696E74202D E94D58B9393AE82F AC82B782E982BD82DF82CC8AEE E707074>

03マイクロ波による光速の測定

伝熱学課題

<4D F736F F D208D5C91A297CD8A7793FC96E591E6398FCD2E646F63>

報道関係者各位 平成 26 年 5 月 29 日 国立大学法人筑波大学 サッカーワールドカップブラジル大会公式球 ブラズーカ の秘密を科学的に解明 ~ ボールのパネル構成が空力特性や飛翔軌道を左右する ~ 研究成果のポイント 1. 現代サッカーボールのパネルの枚数 形状 向きと空力特性や飛翔軌道との

ニュートン重力理論.pptx

水理学(水理学Ⅰ、水理学Ⅱの復習)

<4D F736F F D20824F B834E835882CC92E8979D814690FC90CF95AA82C696CA90CF95AA2E646F63>

RMS(Root Mean Square value 実効値 ) 実効値は AC の電圧と電流両方の値を規定する 最も一般的で便利な値です AC 波形の実効値はその波形から得られる パワーのレベルを示すものであり AC 信号の最も重要な属性となります 実効値の計算は AC の電流波形と それによって

Microsoft Word - 実験テキスト2005.doc

ÿþŸb8bn0irt

DVIOUT

<4D F736F F F696E74202D AB97CD8A E631318FCD5F AB8D5C90AC8EAE816A2E B8CDD8AB B83685D>

<4D F736F F D208D5C91A297CD8A7793FC96E591E6328FCD2E646F63>

Microsoft PowerPoint - 第9回電磁気学

材料強度試験 ( 曲げ試験 ) [1] 概要 実験 実習 Ⅰ の引張り試験に引続き, 曲げ試験による機械特性評価法を実施する. 材料力学で学ぶ梁 の曲げおよびたわみの基礎式の理解, 材料への理解を深めることが目的である. [2] 材料の変形抵抗変形抵抗は, 外力が付与された時の変形に対する各材料固有

軸受内部すきまと予圧 δeff =δo (δf +δt ) (8.1) δeff: 運転すきま mm δo: 軸受内部すきま mm δf : しめしろによる内部すきまの減少量 mm δt: 内輪と外輪の温度差による内部すきまの減少量 mm (1) しめしろによる内部すきまの減少量しめしろを与えて軸受

Microsoft Word - cavitation.doc

国土技術政策総合研究所 研究資料

2014年度 名古屋大・理系数学

線積分.indd

数値流体解析 (CFD) によるスプレー性能の最適化ブリテン No.J955A 数値流体解析 (CFD) による スプレー性能の最適化

平成 23 年度 JAXA 航空プログラム公募型研究報告会資料集 (23 年度採用分 ) 21 計測ひずみによる CFRP 翼構造の荷重 応力同定と損傷モニタリング 東北大学福永久雄 ひずみ応答の計測データ 静的分布荷重同定動的分布荷重同定 ひずみゲージ応力 ひずみ分布の予測 or PZT センサ損

Chap2.key

<4D F736F F F696E74202D2091E6328FCD E9F8CB392E88FED944D936093B1298D758B F E291E892C789C1292E B8CDD8

Transcription:

圧損型粘度計による日用品流体の粘度測定 望月聡 小椋勝仁 山田重良 桜木俊一 概要 : 私達の日常生活で利用される 醤油やソース ドレッシングやマヨネーズなどの食品流体 さらに シャンプーやリンス 液体洗剤などの日用品流体は 充填機を利用して容器に自動充填され製品として出荷されている しかし この充填プロセスにおいて 日用品流体の流体物性 ( 粘度または粘性係数 ) に起因する様々な技術課題が存在している 一般的に 日用品流体は非ニュートン流体の性質を保有するものが多く 起動時にポンプに非常に大きな起動トルクを要求する場合や 充填時の発泡現象による容器からの吹きこぼれ問題の発生 ゴマや玉ねぎなどの具材入り食品流体の粘度測定法の確立など未解決な問題が多数存在する これまで これらの日用品流体は 回転トルク式の粘度計を利用して粘度の測定が実施されてきたが 回転円筒と固定円筒の間の間隙が mm 以下と非常に狭く 具材入りのタレなどの測定では具材が間隙に挟まってしまい正確な測定値が得られないという問題も発生していた 筆者らは 具材などの混入粒子の影響を受けることのない新たな粘度測定法として 円管の異なった 点間で発生する圧力損失を利用した粘度計を製作し その測定値と回転トルク式粘度計の測定値を比較検討し圧力損失型粘度計の精度を検証した またこの研究により得られた結果から 充填機から容器に噴出される食品流体の噴流構造の制御に有効な知見を得ることができた キーワード : 非ニュートン流体 粘度 ハーゲン ポアズイユの法則 レイノルズ数 乱流ジェット. 緒言ドレッシングやゴマダレ マヨネーズなどの食品流体 シャンプーやリンス 液体洗剤などの日用品流体は充填機を利用して容器に自動充填され製品として出荷されている これらの日用品流体は非ニュートン流体の性質を持つものが多く 流動特性も未解明なものが多い 筆者らは これまでの回転トルク式粘度計 ( 図 参照 ) では対応できなかった特殊な性状の流体にも対応できるように充填機流体回路の一部に圧損型粘度計を組み込み リアルタイムでの粘度特性の測定を試みた. 圧損型粘度計の理論的背景. 管路内速度分布と圧力損失図 は充填機に装着された圧損型粘度計の模式図である 距離 l ( 本実験では l = m とした ) だけ離れた 点間の圧力 P と P を測定することにより 流体の粘度計測が可能になる 図 3 は 管路内流速分布を求めるための解析モデルを示す () 同図より 半径 長さ の微小円筒流体要素の 左右底面に作用する圧力による力の合計は p ( p ) () 一方 微小円筒流体要素の側面に作用するせん断 図. 回転トルク式粘度計 図. 圧損型粘度計の概略図

水の粘度 (mpa s) y R φd p τ 図 3. 円管内流速分布の解析モデル 応力を τ とすると 円筒の全側面に作用するせん 断応力による力は となる 定常状態 では この力と圧力による力 () 式が釣り合うので これから τ を求めると () (3) 壁から円筒側面までの距離を y とすると = R - y であるから d = -dy が成り立ち ニュートン流体 に対するせん断応力の式は μ を粘性係数として du du () dy d () 式を (3) 式に代入して変形すると du d (5) 両辺を積分すると u p + d C (6) 積分定数 C は 管壁 = R において u = が成立 することから定められ 以下の管路内速度分布の 式が得られる u ( R ) (7) (7) 式より 速度分布は回転放物面になっているこ とが分かる 上式中の / は 下流方向への圧 力勾配であるから / < であることに注意が 必要である. ハーゲン ポアズイユの法則を利用した流体 の粘度測定 流速分布の式 (7) を用いて流量 Q を計算する 半 径が と +d に囲まれた環状部分の面積 πd を通過する流体の流量は d u である したが V u max x d って 管を流れる全流量 Q は R R Q d u ( R R 8 ) d (8) 長さが l の区間における圧力の降下量を Δp とする と - /=Δp / l であるから 上式は R p d p Q 8 l 8 l (9) また 管路の断面積を A 断面平均流速を V とす れば 連続の式より Q AV ( d / ) V こ の式を (9) 式に代入すると次の関係式を得る p d /(3lV ) () すなわち 圧力差 Δp と断面平均流速 V を知るこ とにより 流体の粘度を測定することができる ただし () 式は () 式の層流せん断応力の仮定か ら層流状態の流れに対してのみ成立することに注 意を要する 3. 実験結果と考察 3. 圧損型粘度計の検定 本研究で使用する圧損型粘度計の測定精度を検 定するために 粘度の知られている流体として水 を用いて粘度測定を実施した 図 に測定結果を 示す 水はニュートン流体であるため流れのせん 断速度に依存せず一定の粘度を持つ 図中では は一定値.33 mpa s の値を示して いる 一方 圧損型粘度計は Re 数 3 の値から 算出された臨界流量.5 cc/s 以上の流量で大き 5 3 層流 乱流 6 流量 (cc/s) (kpa センサ ) ( 測定値.33 mpa s) 臨界流量 (.5 cc/s) 理科年表掲載粘度 :. mpa s ( ) φ d = 6 mm 図. 水の粘度測定による粘度計の検定

リンス粘度 (mpa s) 醤油粘度 (mpa s) く逸脱した値となっている また cc/s 以下の 微小流量領域では 圧力センサーの発生電圧が小 さく誤差の増大が見て取れる なお 本実験では 圧損測定部の管内径は d = 6 mm である 3. 圧損型粘度計による醤油粘度の測定 充填機によるボトル充填用流体の例として 醤 油の粘度測定を実施した 図 5 に結果を示す 醤 油は水や空気などと同様のニュートン流体である ため 流れのせん断速度に依存せず一定の粘度を 持つ では一定値の.66 mpa s を示 している 一方 圧損型粘度計では 種類の圧力 センサー ( 定格圧力 kpa と 3kPa) による測定 結果を示している この場合も 臨界流量を超え た乱流領域では測定値の大きな逸脱が見られる 充填機を用いた醤油の充填では 充填流量を増 やしていくとボトル内で発泡現象が発生し 吹き こぼれによる不良品の発生が起こるが 今回の研 6 5 3 層流 3 5 流量 (cc/s) 乱流 図 5. 醤油粘度測定結果 ( 測定値.66 mpa s) 臨界流量 (.93 cc/s) (kpa) (3kPa) φ d = 6 mm 究で 臨界流量を境界として発泡現象の有無の領 域に分かれることが明らかになった 図 6 は醤油 の噴出流量と発泡現象の関係を撮影したものであ るが 臨界流量を超えると急激な発泡現象が現れ る これは ノズルから噴出される液体ジェット の構造が 臨界流量を超える流量では乱流ジェッ トに遷移し 周囲の空気を大量にジェット中に吸 引することにより発泡現象を引き起こすものと推 察される 圧損型粘度計ではグラフの形状から臨 界流量を明確に知ることができるので便利である 3.3 非ニュートン流体 ( リンス ) の粘度測定 日用品流体では 流れのせん断速度に依存して 粘度が変化する非ニュートン流体の性質を持つも のが多い では 回転円筒の回転速度 から流れのせん断速度 du/dy が容易に計算でき るが 圧損型粘度計では管壁面での流れの速度勾 配 ( 微分係数 ) を知る必要がある 前述の (7) 式を で微分すると du () d また - /=Δp / l と = R ( 壁面 ) の関係を () 式に代入して 次の壁面せん断速度の関係を得る du dy y du d R p d p R l l () 図 7 は 測定流体にリンスを使用したときのせん 断速度と粘度の関係の測定値を示す 同図より との測定結果は非常に良く (a) cc/s (b) 5 cc/s (c) 5 cc/s ( 層流 ) ( 臨界領域 ) ( 乱流 ) 図 6. ノズルからの噴出流量と発泡現象の発生 5 5 35 3 5 5 5 5 5 5 3 せん断速度 du/dy (/s) (3 kpa) (CPA 型 ) φ d = mm 図 7. 非ニュートン流体 ( リンス ) の粘度特性

玉ねぎジュース粘度 (mpa s) 液体洗剤粘度 (mpa s) 一致している このように 非ニュートン流体は ポンプ起動時の低流速時 ( 低せん断速度時 ) に粘 度が非常に高く ( 水の数千倍 ) 流速の増加とともに 急速に粘度の低下が起こる また このリンスの ような高粘度流体はポンプの最大吐出流量におい てもレイノルズ数が非常に小さく層流状態に保た れている 3. 具材入り食品流体の粘度測定 図 8 は玉ねぎの微粒子が混入している玉ねぎジ ュースの粘度測定結果である この図より 計測 流体はせん断速度に依存して粘度が変化する非ニ ュートン流体であることが分かる また 回転粘 度計 (CPA 型 ) の場合 具材と回転面の接触干渉 と思われる不規則な変動が見られる 一方 圧損 粘度計 ( 管径 mm) と隙間ギャップを大きく設定 できる Spindle 型は滑らかなデータの 接続ができており 具材の接触干渉の影響が微小 であることが窺える 5 5 35 3 5 5 5 6 8 せん断速度 du/dy (/s) 図 8. 玉ねぎジュースの粘度特性 3.5 その他の特異な粘度特性を持つ流体 (3 kpa) (CPA 型 ) (Spindle 型 ) φ d = mm 図 9 は液体洗剤の粘度測定結果であるが 特異 な性質を示した による第 回目の計 測では の測定値が示すようにせん断 速度に依存しない一定の粘度を持つニュートン流 体の性質を示した 次に同じ計測サンプルを 週 間後に測定した第 回目の測定では非ニュートン 流体の性質を示す結果となった さらに再度確認 のために 回目の測定から 週間後に 3 回目の 測定を行ったが同じく非ニュートン流体の性質を 6 8 6 図 9. 液体洗剤の粘度特性 示した この原因として 管路内の流動を繰り返 すことにより分子構造などの変化を誘発したとも 考えられるが 定かではない. 結言 ハーゲン ポアズイユの法則に基づく圧損型粘 度計を製作し 充填機流体回路に組み込み 各種 日用品流体の粘度測定を試みた その結果 回転 式粘度計とほぼ同等の精度で粘度測定が可能であ ることを確認した また 充填ボトル内での発泡 現象と噴流ジェット構造の相関も明らかになった 参考文献 資料 () 菊山功嗣 佐野勝志 流体システム工学 共 立出版 7, pp. 76~78. 著者紹介 3 5 せん断速度 du/dy (/s) 望月聡静岡理工科大学理工学部機械工学科 年 小椋勝仁靜甲株式会社 IT 推進室 山田重良靜甲株式会社専務付 桜木俊一静岡理工科大学理工学部機械工学科教授博士 ( 工学 ) E-mail: sakuagi@me.sist.ac.jp (3 kpa) 回目 (3 kpa) 回目 (3 kpa) 3 回目 (Spindle 型 ) φ d = mm

本論文は 日本技術士会中部本部 7 年 度研究業績発表会で 優秀論文 発表賞を受賞しました