• 研究課題をさがす
  • 研究者をさがす
  • KAKENの使い方
  1. 前のページに戻る

分子動力学法による超臨界状態における流体の動的液体構造

研究課題

研究課題/領域番号 06214223
研究種目

重点領域研究

配分区分補助金
研究機関法政大学

研究代表者

片岡 洋右  法政大学, 工学部, 教授 (30025407)

研究期間 (年度) 1994
研究課題ステータス 完了 (1994年度)
配分額 *注記
1,900千円 (直接経費: 1,900千円)
1994年度: 1,900千円 (直接経費: 1,900千円)
キーワード超臨界流体 / 分子動力学 / 相互拡散係数 / 自己拡散係数 / レナード・ジョーンズ流体 / 濃度の揺らぎ / 圧縮率
研究概要

2次元レナード・ジョーンズ(LJ)混合液体(Ar-Kr)における超臨界流体の分子動力学シミュレーションを行った。計算方法は分子動力学法で最も単純な体積と全エネルギー一定の分子動力学法である。単位はArのLJ単位である。濃度はモル分率x=0.50とした。
状態を確かめるためにまず、圧力Pを計算し、その密度ρ・温度T変化を確かめた。以下で調べる温度T0.67 密度ρ0.3が臨界点に近いことを確認できた。
輸送係数のうち特に自己拡散係数D1,D2と相互拡散係数Dを重点的に計算した。
相互拡散係数は熱力学因子Qと速度の時間相関関数の積分Lの積で書かれる。(D=Q*L)熱力学因子Qは密度ρが0.3付近で小さい。Qは化学ポテンシャルの微分で書かれ、その計算は動径分布関数の積分Gab(Kirkwood-Buffのパラメータ)より行われた。積分Gabの値が妥当なものであることはPVT関係から求めた圧縮率κ_TとKirkwood-Buffのパラメータより計算したものを比較して次のように確かめた。温度T0.67 密度ρ0.3が臨界点に近いのでこの付近で圧縮率は大きな値となる。この傾向は再現されており半定量的な一致が確かめられた。
時間相関関数の積分は重心の変異の2乗平均の長時間部分での時間tに対する傾きから計算した。
得られた相互拡散係数と自己拡散係数をT=0.67で密度ρに対してプロットして、相互拡散係数の異常な密度依存性が得られた。
濃度の揺らぎの大きい例の分子配置をみると粒子が密な部分と粗な部分がまず見える。さらに少しKrが集まりやすい傾向が見える。

報告書

(1件)
  • 1994 実績報告書
  • 研究成果

    (4件)

すべて その他

すべて 文献書誌 (4件)

  • [文献書誌] K.Yasuoka,M.Matsumoto and Y.Kataoka: "Evaporation and Condensation at a ligid Surface,I,Argon" J.Chem.Phys.101. 7904-7911 (1994)

    • 関連する報告書
      1994 実績報告書
  • [文献書誌] M.Matsumoto,K.Yasuoka and Y.Kataoka: "Evaporation and condensation at a liqid surface,II.Methanol" J.Chem.Phys.101. 7912-7917 (1994)

    • 関連する報告書
      1994 実績報告書
  • [文献書誌] Y.Kataoka and M.Fujita: "An Equation of state for the self diffusion coefficient in lennard-Jones fluid derived by molecular dynamics simulation" Bull.Chem.Soc.Jpn.68. 152-159 (1995)

    • 関連する報告書
      1994 実績報告書
  • [文献書誌] 片岡洋右: "分子動力学法とモンテカルロ法" 講談社, 179 (1994)

    • 関連する報告書
      1994 実績報告書

URL: 

公開日: 1994-04-01   更新日: 2018-06-07  

サービス概要 検索マニュアル よくある質問 お知らせ 利用規程 科研費による研究の帰属

Powered by NII kakenhi