Abstract
本稿では、マイクロ流体デバイス内の不均一性を測定するために、混合マイクロシュリーレン技術の使用を紹介します。マイクロシュリーレンシステムはスリット板を除去し、ナイフエッジで変調器を置き換えることにより、対物レンズの後側焦点面への容易なアクセスを提供するホフマンモジュレーションコントラスト顕微鏡、から構成されています。マイクロシュリーレン技術の動作原理は、屈折率1-3のばらつきに起因する光偏向を検出することに依存しています。偏向光が逃げるか、それぞれ、明るいまたは暗い帯を生成するためにナイフエッジによって妨害されますか。混合物の屈折率は、その組成物を用いて直線的に変化する場合、像面における光強度の局所的な変化は、光軸に対して垂直な濃度勾配に比例します。マイクロシュリーレン画像は、三次元不均質性によって生成される乱れた光の2次元投影を提供します。
定量的な分析を達成するために、我々は、T-マイクロチャネルで二つの流体を混合する較正手順を説明します。我々は、対応するマイクロシュリーレン画像と密接に相関するT-マイクロチャンネル内の濃度勾配を得るために、数値シミュレーションを行っています。比較すると、マイクロシュリーレン画像のグレースケールの読み出し及びマイクロ流体デバイスにおいて提示濃度勾配との間の関係が確立されます。この関係を用いて、準マイクロシュリーレン画像から混合の不均一性を分析し、マイクロ流体発振器4での測定でマイクロシュリーレン技術の能力を実証することができます。光学的に透明な流体は、マイクロシュリーレン法は、混合プロセスの三次元の特徴を保持する瞬間フルフィールド情報を提供するための魅力的な診断ツールです。
Introduction
流体混合は、多くの工業プロセス及び生物系に見られる重要な問題です。マイクロフルイディクスの出現により、マイクロスケールでの混合による質量輸送機構の間で拡散支配でその課題に多くの注意をもたらしています。定量的な検証が必要な有効マイクロミキサーを設計しているので、いくつかの測定方法は、5-7を開発しました。それにもかかわらず、一般的に効率的なマイクロミキサ5に見られる三次元構造は、一般的な測定技術は、配信に失敗濃度フィールドのより正確な表現を必要とします。による角度8または反応速度6の視聴制限のために、上記の方法は、正確に、混合物の均質性を考慮していない誤った結果を生成することができます。
光学的に透明な微細構造中での混合、光学的に透明な流体、マイクロシュリーレン技術3,9-14用9-13、15または位相勾配16を可視化するために使用されています。簡単な光学レイアウトと高感度の両方からマイクロスケールシュリーレン技術の利点と、光障害を引き起こすが、混合を評価するのに使用するのに適している特定のフロー機能の非侵襲的な調査だけでなく、を可能にします。本稿では、顕微鏡の対物レンズの後焦点面でのナイフエッジを挿入することにより、マイクロスケールシュリーレンシステムを構築し、定量的な分析を実現するためのキャリブレーション手順を説明し、マイクロ流体発振器4での検証測定を報告しています。混合流体の屈折率は、組成物に対して直線的に変化し、ターゲットのマイクロ流体デバイスの厚さは、上と同じになるように測定を実施するために、作動流体が適切に選択されますeは、キャリブレーションに使用されます。種の濃度に加えて、マイクロシュリーレン法は、温度や塩分などの線形屈折率に相関している他のスカラ量の傾きを測定するために拡張することができます。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
マイクロ流体デバイスの作製1。
- T-マイクロチャネルの輪郭を描画する図形レイアウトソフトウェア（ 例えば、AutoCADの）を使用します。 T-マイクロチャンネルでは、2つの供給チャネルは、90ミクロン、幅2,500μmの長さ、および合流チャネルは180ミクロン、幅3000ミクロンの長さです。 1100ミクロンの直径を有する個々の円に、各チャネルの端を接続します。
- マーク「クリア」と暴露し、対象地域のために「暗い」、それぞれ。ネガ型フォトレジストの場合（ 例えば、SU-8）は、T-マイクロチャネルの形状は、「クリア」であり、その周辺で「暗いです」。
- クロム·オン·ガラスマスク上にT-マイクロチャネルのパターンを転写するために442ナノメートルの波長2μmの最小特徴サイズを有するレーザパターン発生器を使用します。
- フォトマスクは、基板（ 例えば、片面研磨されたシリコンウェーハ）と永久エポキシフォトレジストを使用して（ 例えば、、SU-8）は、標準的なリソグラフィプロセスを介して金型を作ります。フォトレジスト層は55.2ミクロンの厚さです。一般的に、フォトレジストの厚さは、対物レンズ17-19の相関の深さよりも薄くする必要があります。
- T-マイクロチャネル20を製造するための金型及びポリジメチルシロキサン（PDMS）などの透明な材料を使用してください。
- 流体接続のために、パンチの外径2mmのステンレス鋼管を使用する貫通孔PDMS上の円形パターンに整列されます。
- PDMSと30秒60 Wの酸素プラズマによるスライドガラスの表面を処理します。スライドガラスにPDMSを取り付けます。 2つの材料の酸化表面は強固な結合を作成します。 120℃で5分間ホットプレート上に接着PDMS構造を配置します。
- 流体接続するためのパンチ穴にテフロンチューブを挿入します。
2.実験セットアップ
- Hoffmaからマイクロスケールシュリーレンシステムを構築凝縮器の前側焦点面にスリット板を除去し、5倍の対物レンズ3の後側焦点面でのナイフエッジと変調器を置き換えることにより、n個の変調コントラスト顕微鏡。対物レンズ17-19の開口数に依存する相関関係の深さは、マイクロ流体デバイスの全体の深さをカバーするのに十分であるべきです。ナイフエッジの表面は、その反射率を低減するために、陽極酸化アルミニウムにより黒化されます。
- Cマウントアダプタを介して、顕微鏡の三眼鏡筒にハイスピードカメラをマウントします。カメラは、ビームスプリッタを介して顕微鏡の光路に面しています。イーサネットケーブルを介して、デスクトップコンピュータにカメラを接続します。そのグレースケールの読み出しは入力輝度に比例するように、カメラを1にガンマ補正を設定します。
- 光源をオンにします。過剰な熱を避けるために、LED（発光ダイオード）照明を使用。
- 画像処理ソフトウェア（ 例えば、<使用/ EM>はMATLAB関数の関数imread）が取得した画像の階調値を取得します。ナイフエッジを削除し、画像の平均階調読み出しが最大値よりも約10％小さくなるように照明、絞りと露光時間を調整します。これは、0％のカットオフのバックグラウンド強度を示し、我々は、8ビットの画像のために230の値を使用します。
- 完全に入射光を遮断するナイフエッジを挿入します。画像の平均グレースケールの読み取りを記録します。これは、100％のカットオフのためのバックグラウンド強度を示し、値は約15 8ビットの画像のためにあります。
- ナイフエッジの位置を調整するように、取得された画像の平均階調読み出し、0％と100％のカットオフ値の中間にあります。ここで遮断の程度は50％に設定されています。
- 構成成分として、互いに完全に混和性であることが知られている屈折率21を有する2つの透明な液体を準備します。屈折INDの依存性を評価するために、混合物の濃度のEXは、文献21を確認するか、グラッドストーン-デール式22を使用しています。曲線が全範囲にわたって非線形であれば、他の流体コンポーネントを選択します。その後、溶液の屈折率は濃度に比例して変化するの下に指定された組成物を選択します。例えば、作動流体として0.05と水の質量分率でエタノール水溶液を希釈してください。
- 試料ステージ上にT-マイクロチャネルを配置します。ナイフエッジにコンフルエントチャンネルパラレル（ 図1）となるようT-マイクロチャンネルを配置します。
- シリンジAは、基準流体（水）として機能する作動流体で満たされ、そしてシリンジBは、他の作動流体（水、エタノールを希釈）が充填されている二つの同一の注射器を準備します。注射器の大きさは、所望の流量Qと、シリンジポンプの仕様に依存する：Q =πD2 V / 4、dは SYRの内径インゲとVは、プランジャの速度です。脈動は、通常、V 23を増大させるために小さな注射器を選択することによって防止することができます。
- ビーカー内のT-マイクロチャネルからの出口流体を収集します。確認テフロン出口管をビーカーの壁に固定され、その端部が、液滴切離しに起因する振動を避けるために、ビーカー中の液体レベルより下です。
3.キャリブレーション
- 混合2流体および参照画像の画像を取得します。
- シリンジポンプの流量を設定された所定のレイノルズ数Reにおいて、Q。Qは、Q =μ（W + D）のRe /μ、ρは作動流体の粘度および密度である4ρ、およびWから計算されます及びDは、それぞれ、T-マイクロチャネルの合流流路の幅および深さです。
- シリンジAとSYと他のポンプと負荷つのポンプringe Bのは、テフロンチューブを介してAとシリンジBを注射器T-マイクロチャネルの2つの入口を接続します。同一の体積流量でT-マイクロチャネルに作動流体を供給するようにシリンジポンプを起動します。
- 定常流が確立するまで待ちます。定常流状態が定常シュリーレンパターンの出現によって定義されます。
- 30fpsのフレームレートで流体混合の20フレームを記録するカメラ制御ソフトウェアを使用してください。
- のみ一定の速度でT-マイクロチャネルの合流チャネルに一つの入口を介して基準流体（水）をシリンジポンプBにロードされているポンプを停止します。
- 定常流状態に到達し、何のシュリーレンパターンが観察されなくなるまで待ちます。
- 光学的な不均一性は、T-マイクロチャネル内に存在しない場合には、参照画像を撮るためにカメラ制御ソフトウェアを使用してください。 30fpsのフレームレートで20フレームを記録します。
- 異なるレイノルズ数で3.1.7に3.1.1を繰り返します。= 1 再、5、10、20及び50複雑な流れ構造は、T-マイクロチャネル24の合流領域に出現しないようにします。
- iおよびjは画素インデックスである参照画像I 0（i、j）は 25で取得された画像I（i、j）を分割する画像処理ソフトウェアを使用します。
- T-マイクロチャンネル内で指定された流体の混合をシミュレートするCFD（計算流体力学）パッケージを採用しています。
- T-マイクロチャネルのジオメトリの三次元モデルを構築します。構造格子への流れドメインを離散化。精度を高めるために、合流及びT-マイクロチャネルの中央領域に細かいメッシュを使用します。
- 流体の物理的性質を割り当て、フロードメインに境界条件を確立します。解決プロセスの間、最後に得られた濃度から濃度依存性の拡散係数を決定します局所濃度を更新するために反復26。
- グリッド研究27を行うことにより、演算結果の感度を調べます。
- =（X I、Y jの ）W：各ノードの（x i、y i）とのxy -planeに、台形則によって、チャネルの深さを横切る濃度フィールドの平均値を取るようにCFDポスト処理ツールを採用{Σkの [（X I、Y jを 、Z K）W + W（X I、Y jを 、Z K +1）]∙（Z kを 1 - ZのK）/ 2} Dが / D 28、チャネル深さ。 derivaを計算するために、中央差分スキームを使用してクロスストリーム方向に対する濃度の的：（W /∂∂Y）I、J = [（X I、Y jを+1）W - W（X I、Yの J -1）]（/ Y jを 1 - Y jを -1）。
- 正と負の両方の勾配のために、 私は / I 0（3.2）で得られたグレースケール値の比率や、Y = 0（中心線、流れ方向として指定された場所で（3.4）で得られたW /∂∂yの質量分率の勾配を抽出方向）または様々な指定されたX（クロスストリーム方向）。
- 結果をプロットし、I / I 0との間の関係を決定し、W /∂∂Y、I / I 0 /∂Y + C 2、W C 1∂（C 1、C 2は定数）、線形回帰と。
4.定量
- 繰り返して、ターゲットマイクロ流体デバイス内で混合するために3.1〜3.2を繰り返します。ターゲットマイクロ流体デバイスの深さは同一またはT-マイクロチャンネルのそれに近いものであるべきです。非定常現象が予想される場合、代わりにステップ3.1.4でビデオクリップ（一連の画像）を取得。露光時間が2.4、2.5、3.1.4と3.1.7で使用される値と同じである必要がありながら、フレームレートは、明らかに過渡流れのダイナミクスを解決するのに十分高くなければなりません。
- ターゲットマイクロ流体デバイス内の質量分率の勾配にグレースケール値の比を変換するために、ステップ3.6で得られた関係を使用してください。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
質量分率の正および負の勾配の両方について異なるレイノルズ数下のグレースケール比率I / I 0は、T-マイクロ流路の途中に現れる対称バンドと（ 図2）が示されています。低レイノルズ数では、シュリーレンバンドの尾が原因ミキシングインターフェースを介して分散液に拡大し、ぼやけます。レイノルズ数が増加するにつれて、拡散長は、より狭い帯域につながる短縮します。異なる下流の場所で、強度変化ΔIの変動は/ Iは、クロスストリーム方向0を定量（ 図3）が示されています。較正プロセスの結果は、（ 図4Aに示されています および4B）。 I / I 0〜W /∂∂yの関係が線形とレイノルズNUMとは無関係ですBER。回帰分析から、I / I 0 = -110∂∂Y +∂W /∂Y> 0、I / I 0 = -160∂用1.03 / W /∂Y + 0.83 W∂W /∂Y <0のため、∂∂Y / Wは μmで-1です。相対的な不確実性が±3.8％であり、それぞれ、 図4に3.2％、±します。データ点は平準どこ検出限界に達します。これは、正および負の勾配の傾きのずれが3珍しいことではないことに留意されたいです。これらの式を使用して、マイクロ流体発振器4の時間質量分率の勾配の変化（ 図5）に見られます。混合インタフェースは、キャビティ領域内に偏向され、不安定cを流すれますommences。このビデオ図は明らかに、マイクロ流体発振器内の流れの振動の性質を明らかにし、マイクロ流体デバイスの時間分解全視野の濃度勾配を捕捉するためにマイクロシュリーレン法の能力を実証します。
光学装置の図1の回路図。ナイフエッジの向きは、屈折率の正の勾配で暗いバンドが生成されます。光は、屈折率を増加させる方向に偏向します。対物レンズは、画像を反転させているため、ブロッキング- Y領域は、歪んだ光を遮蔽し、暗いバンドが生成されます。
unde T-マイクロチャンネル内で混合するためのグレースケールの読み出し図2.レシオ異なるフロー構成をrは正および負の勾配は、それぞれ、暗い部分と明るいバンドをもたらします。レイノルズ数が増加すると、帯域がより濃縮なります。
図3の正および負の両方のグラデーションのクロスストリーム方向の強度変化の変化。 再 = 1及びRe = 5。
質量分率の勾配とグレースケール比と図4の関係。正と負の両方の勾配については、グレースケールの比率が質量分率の勾配で直線的に変化します。
915 / 52915fig5.jpg "/>
図5（ビデオ図）。 再 = 250フロー振動を介して混合特性でマイクロ流体発振器における質量分率勾配の進化に成功マイクロシュリーレン法によって捕捉されます。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
マイクロ流体デバイス内の流体の混合のために、マイクロスケールのシュリーレン法は、光強度の変化を定量化することを介して濃度勾配の大きさを測定することができます。この技術の原理は、光伝播の交代を検出することに依存しているため、作動流体とマイクロ流体デバイスは、入射光に対して透明でなければなりません。作動流体の予備的評価が不可欠となるように加えて、プロトコルは、溶液およびその組成の屈折率との間の線形関係を必要とします。本明細書において実証エタノール水溶液に加えて、マイクロシュリーレン法が正常塩分勾配29とsolutocapillary対流30を測定するために適用されます。較正手順で使用される正確な測定値、絞り範囲、照明レベル、暴露時間、対物レンズとマイクロチャネルの深さについて定量手順で使用されるものと同一でなければなら。また、対物レンズの相関関係の深さは、マイクロ流体デバイスの全体の深さをカバーするのに十分に大きくなければなりません。
T-マイクロチャンネル内での混合のキャリブレーション処理は、マイクロスケールシュリーレン技術の正確な定量において最も重要なステップです。提案手法の実装を成功させるために、ユーザは、適切にチューブ接続を揃える余分な熱を減少させるためにLED光源を使用し、フロー発振23を回避するために、流体送達のための小さな注射器や空気圧を利用する、低レイノルズ数でのキャリブレーション手順を実施する必要があります24、および高次光学効果31を除去するために、焦点にマイクロ流体デバイスを配置します。最低測定勾配（明るいパターン、Y <0は、∂/∂W）W /∂∂yの最高測定可能な勾配（明暗パターン、一方、カメラのダイナミックレンジにリンクされています> 0）が到達したときにナイフエッジが完全に遮断偏向光を。濃度勾配の広い範囲を検出するために、高ISO値であれば黒つぶれや白とびが発生しないように有利です。マイクロシュリーレンシステムが識別することができない、それ以下の検出限界は、カメラが解決することができる最小の強度変化に依存します。最小限の強度変化は、ノイズの程度や色調、階調のレベルによって制限されます。したがって、大きなピクセル深度で高感度カメラは、低信号用途に望まれます。
マイクロシュリーレン法の重要性は、二つの折り目です。一方で、それは簡単な光学構成でリアルタイムに非定常フルフィールド測定を可能にします。いかなる異物が流れ場を乱すために導入されないように、他方では、非侵襲的です。マイクロシュリーレン法はmicrof三次元不均質性の二次元投影を生成するので既存の方法によって包まれたままluidicデバイス、複雑な混合現象が明瞭に見ることができます。この技術の将来のアプリケーションは、電気化学的処理中に濃度勾配を定量化やマイクロ流環境中の微生物の走化性を研究するために栄養勾配を決定することを含みます。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
著者らは、開示することは何もありません。
Acknowledgments
この作品は、認可番号101から2221-E-002から064-MY3の下で台湾の科学技術省によってサポートされていました。
Materials
| Name | Company | Catalog Number | Comments |
| Permanent Epoxy Negative Photoresist | MicroChem | SU-8 2150 | |
| single side polished silicon wafer | Light Technology | S4W1PP5SABUP1 | p-type, diameter: 100±0.5 mm, thickness: 525±25 mm, orientation: (1 0 0) |
| syringe pump | kdScientific | kds210 | |
| syringe, i.e. 10 ml | Terumo | SS-10L2138 | |
| Hoffman modulation contrast microscope | Leica Microsystems | DM IL LED | |
| 5X objective lens | Leica Microsystems | N PLAN | NA = 0.12 |
| knife-edge | custom made part | ||
| camera, i.e. high speed | Integrated Design Tools | NX7-S1 | |
| C-mount adapter HC 0.63x | Leica Microsystems | 541537 | |
| camera operating software | Integrated Design Tools | MotionPro X Studio 2.02.01 | |
| polydimethylsiloxane (PDMS) | Dow Corning | Sylgard-184 silicone elastomer kit | |
| Teflon tubing, i.e. O.D. x I.D. 1/16 in. x 0.031 in. | Supelco | 58700-U | |
| micro glass slide | Matsunami Glass | S2215 | |
| hot plate | Yeong-Shin | HP-303DN | |
| distilled water, i.e. HPLC grade | Alps Chemicals | ||
| ethyl alcohol, i.e. reagent grade | Nihon Shiyaku Reagent | EA448652 | |
| image processing software | Mathworks | MATLAB R2009a | |
| computational fluid dynamics package | ESI group | CFD-ACE+ 2008 |
References
- Merzkirch, W. Flow Visualization. , Academic Press. New York. (1974).
- Settles, G. S. Schlieren and Shadowgraph Techniques. , 2nd edn, Springer. New York. (2001).
- Sun, C. -l, Hsiao, T. -h Quantitative analysis of microfluidic mixing using microscale schlieren technique. Microfluidics and Nanofluidics. 15 (2), 253-265 (2013).
- Sun, C. -l, Sun, C. -Y. Effective mixing in a microfluidic oscillator using an impinging jet on a concave surface. Microsystem Technologies. , (2011).
- Strook, A. D.
- Wheat, P. M., Posner, J. D. Quantifying mixing using equilibrium reactions. Physics of Fluids. 21 (3), 037101 (2009).
- Liu, R. H. Passive mixing in a three-dimensional serpentine microchannel. Journal of Microelectromechanical Systems. 9 (2), 190-197 (2000).
- Munson, M. S., Yager, P. Simple quantitative optical method for monitoring the extent of mixing applied to a novel microfluidic mixer. Analytica Chimica Acta. 507 (1), 63-71 (2004).
- Bradfield, W. S., Sheppard, J. J. Microschlieren-a technique for the study of details in compressible flow. Aero/Space Engineering. 5 (18), 37-56 (1959).
- Scroggs, S. D., Settles, G. S. An experimental study of supersonic microjets. Experiments in Fluids. 21, 401-409 (1996).
- Phalnikar, K. A., Alvi, F. S., Shih, C. 31st AIAA Fluid Dynamics Conference & Exhibit. , Anaheim, California. (2001).
- Phalnikar, K. A., Kumar, R., Alvi, F. S. Experiments on free and impinging supersonic microjets. Experiments in Fluids. 44 (5), 819-830 (2008).
- Micro-schlieren characterization of a high momentum micro-fluidic actuator. Naughton, J. W., Bishop, D. S., Walrath, D. E., Lindberg, W. R. 22rd AIAA Aerodynamic Measurement Technology and Ground Testing Conference, St. Louis, MS, USA, , The American Institute of Aeronautics and Astronautics (AIAA). (2002).
- Chen, H. -M., Wu, H. -W., Lee, J. -N., Ting, C. -C. 23th National Conference of the Chinese Society of Mechanical Engineers. , YongKang, Tainan, Taiwan. (2006).
- Davidson, G. P., Emmony, D. C. A schlieren probe method for the measurement of the refractive index profile of a shock wave in a fluid. Journal of Physics E: Scientific Instruments. 13, 92-97 (1980).
- Xie, H. Schlieren confocal microscopy for phase-relief imaging. Optics Letters. 39, 1238-1241 (2014).
- Inoué, S., Spring, K. R. Video Microscopy: The Fundamentals. , 2nd edn, Plenum Press. New York. (1997).
- Meinhart, C. D., Wereley, S. T., Gray, M. H. B. Volume illumination for two-dimensional particle image velocimetry. Measurement Science and Technology. 11 (6), 809-814 (2000).
- Olsen, M. G., Adrian, R. J. Out-of-focus effects on particle image visibility and correlation in microscopic particle image velocimetry. Experiments in Fluids. 29 (1), S166-S174 (2000).
- Friend, J., Yeo, L. Fabrication of microfluidic devices using polydimethylsiloxane. Biomicrofluidics. 4 (2), 026502 (2010).
- Wohlfarth, C. Landolt-Börnstein: Numerical Data and Functional Relationships in Science and Technology. Optical Constants: Refractive Indices of Pure Liquids and Binary Liquid Mixtures (Supplement to III/38). Lechner, M. D. 47, Springer. Berlin. (2008).
- Heller, W. Remarks on refractive index mixture roles. Journal of Physical Chemistry. 69 (4), 1123-1129 (1965).
- Korczyk, P. M., Cybulski, O., Makulska, S., Garstecki, P. Effects of unsteadiness of the rates of flow on the dynamics of formation of droplets in microfluidic systems. Lab on a Chip. 11 (1), 173-175 (2011).
- Dreher, S., Engler, M., Kockmann, N., Woias, P. Theoretical and experimental investigations of convective micromixers and microreactors for chemical reactions. Micro and Macro Mixing: Analysis, Simulation and Numerical Calculation. Bockhorn, H., Mewes, D., Peukert, W., Warnecke, H. .-J. , Springer-Verlag. Berlin, Germany. (2010).
- Huang, C. Y., Gregory, J. W., Sullivan, J. P. A Modified schlieren technique for micro flow visualization. Measurement Science & Technology. 18 (5), N32-N34 (2007).
- Tyn, M. T., Calus, W. F. Temperature and concentration dependence of mutual diffusion coefficients of some binary liquid systems. Journal of Chemicaland Engineering Data. 20 (3), 310-316 (1975).
- Celik, I. B. Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. Journal of Fluids Engineering. 130 (7), 078001 (2008).
- Tasić, A. Ž, Djordjević, B. D., Grozdanić, D. K., Radojković, N. Use of mixing rules in predicting refractive indices and specific refractivities for some binary liquid mixtures. Journal of Chemical and Engineering Data. 37 (3), 310-313 (1992).
- A comparison of coarse-resolution numerical simulation with experimental measurements of wind turbine aerodynamic performance. Hsiao, P. -J., Chen, S. -T., Hsiao, T. -H., Sun, C. -l 37th National Conference on Theoretical and Applied Mechanics & The 1st International Conference on Mechanics, , ScienceDirect. Hsinchu, Taiwan. (2013).
- Sun, C. -l, Huang, C. -Y. Microscale schlieren visualization of near-bubble mass transport during boiling of 2-propanol/water mixtures in a square capillary. Experiments in Fluids. 55 (7), 1778 (2014).
- Panigrahi, P. K., Muralidhar, K. Schlieren and Shadowgraph Methods in Heat and Mass Transfer. , Springer. New York. (2012).
