JP2010182106A - 撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム - Google Patents

撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム Download PDF

Info

Publication number
JP2010182106A
JP2010182106A JP2009025148A JP2009025148A JP2010182106A JP 2010182106 A JP2010182106 A JP 2010182106A JP 2009025148 A JP2009025148 A JP 2009025148A JP 2009025148 A JP2009025148 A JP 2009025148A JP 2010182106 A JP2010182106 A JP 2010182106A
Authority
JP
Japan
Prior art keywords
solid
liquid
model
gas
molecule
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2009025148A
Other languages
English (en)
Inventor
Hiroshi Teranishi
浩 寺西
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Motor Corp
Original Assignee
Toyota Motor Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Toyota Motor Corp filed Critical Toyota Motor Corp
Priority to JP2009025148A priority Critical patent/JP2010182106A/ja
Publication of JP2010182106A publication Critical patent/JP2010182106A/ja
Pending legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

【課題】撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法、及びコンピュータ上でシミュレーションを行うシミュレーションプログラムを提供する。
【解決手段】固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、分子動力学法及び／又は分子軌道法を用いて、固気モデル及び／又は固液界面モデルの初期状態における、固体基質への、気体分子又は液体分子の吸着エネルギーを計算するステップと、初期状態の計算結果を解析するステップと、計算結果の解析を基に、固気モデル及び／又は固液界面モデルの撥液性の良否を判定するステップと、判定に基づいて、固気モデル及び／又は固液界面モデルを最適な撥液性被膜を有するモデルとして選択するステップと、を有することを特徴とする、撥液性被膜の探索方法。
【選択図】図２

Description

本発明は、分子動力学法及び／又は分子軌道法に基づいた、撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法、及びコンピュータ上でシミュレーションを行うシミュレーションプログラムに関する。
エンジンの内部には、燃料やエンジンオイル、又はそれらの不完全燃焼生成物が徐々に堆積する。この堆積物のことをエンジンデポジットと呼び、過度のデポジットは自動車の排気や運転性に悪影響を及ぼす。このようなエンジンデポジットを抑制するためには、その生成メカニズムを明らかにする必要がある。しかし、エンジンデポジットの性質や組成物の混合割合等を把握するためには、エンジンデザイン、運転条件及び堆積部位等の複雑な要因を複合的に解析しなければならないため、エンジン全域のデポジットの生成メカニズムを詳細に解析する手法は未だ確立されていない。
エンジンデポジットを除去、抑制するための技術はこれまでにも開発されている。特許文献１には、燃焼変動を回避して、ピストン頂面等の燃焼室内部に生成されたデポジットを適切に洗浄することを目的とした、少なくとも筒内噴射用インジェクタを備えた、アルコール等の代替燃料を同時に使用可能なシステムを搭載した自動車などの車両における、内燃機関の燃料噴射制御装置の技術が開示されている。
特許文献２には、燃焼室内に直接燃料を噴射する直噴ガソリンエンジンをスプレーガイド燃焼モードで運転するように構成した際に、運転時間の増大に伴い、噴霧成層混合気の形状が変化し、燃焼安定度が低下した場合に、燃焼安定度の回復を図り、安定した燃焼を継続させることを課題とした、直噴形エンジンの制御装置の技術が開示されている。
特許文献３には、エンジンデポジットの予防および抑制のため、主に燃料組成物中に使用することを目的とする、ポリアルキルフェノキシアルカノールのアミノカルバメートに関する技術が開示されている。
特開２００８−２７４８０１号公報 特開２００８−２６７２５２号公報 特開平１０−１３９８１１号公報
特許文献１に開示された技術は、あくまでもアルコール含有燃料を用いた結果発生するエンジンデポジットの洗浄のみをその課題としているため、その他の燃料や、エンジンオイル等を用いた際に発生するエンジンデポジットを洗浄する手段として、そのまま応用できるとは考え難い。また当該技術は、エンジンデポジットを洗浄することを目的としているため、エンジンデポジットの予防・抑制などには適していない。
また、特許文献２に開示された技術は、あくまでも燃焼室内に直接燃料を噴射する直噴ガソリンエンジンにおいて発生するエンジンデポジットの抑制のみをその課題としているため、その他のエンジンを用いた際にそのまま応用できる技術ではない。また当該文献には、ガソリンエンジンのエンジンデザインに基づいた物理的考察についての言及はあるものの、エンジンデポジットに含有されている化学成分やその含有比、及び当該化学成分が各々どのような悪影響をエンジンにもたらすかについての化学的考察については全く述べられていない。
また、特許文献３は、実施例において、単気筒エンジンを用い、非加鉛のレギュラーオクタン価のガソリンを基材燃料として、当該文献中に開示されたポリアルキルフェノキシアルカノールのアミノカルバメートを、当該基材燃料に添加する前後のエンジンデポジットの量の比較実験結果について開示している。しかし、当該文献には、単気筒エンジン以外のエンジンを用いた場合や、非加鉛のレギュラーオクタン価のガソリン以外の基材燃料（例えば、ディーゼル燃料等）を用いた場合において、エンジンデポジット量が低減するとの実験結果は何ら開示されておらず、当該実施例において採用されなかったエンジンや基材燃料を用いる場合には、最適な燃料添加剤を得るために、当該文献に開示された技術的範囲を超えて、更なる試行錯誤が必要であると考えられる。
このように、エンジンデポジットによって引き起こされる諸問題を課題とする従来技術は、比較的狭い技術的範囲において適用されるに過ぎないものであった。当該従来技術はいずれも、あらゆる生成要因におけるエンジンデポジットの生成メカニズムについて統一的に考察がなされているものではないため、当該技術的範囲を超えて当該従来技術を適用しようとすると、材料の最適化等において試行錯誤が必要であると予想される。
本発明は、上記実状を鑑みて成し遂げられたものであり、分子動力学法及び／又は分子軌道法に基づいた、撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法、及びコンピュータ上でシミュレーションを行うシミュレーションプログラムを提供することを目的とする。
本発明の撥液性被膜の探索方法は、分子動力学法及び／又は分子軌道法に基づいた、撥液性被膜の探索方法であって、固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、前記判定に基づいて、前記固気モデル及び／又は前記固液界面モデルを最適な撥液性被膜を有するモデルとして選択するステップと、を有することを特徴とする。
このような構成の撥液性被膜の探索方法は、前記分子動力学法及び／又は前記分子軌道法を用いることによって、多角的に前記固気モデル及び／又は前記固液界面モデルについて計算・解析し、当該モデルの良否を判定し、さらに、当該判定結果を基に最適な撥液性被膜を選択することができる。
本発明の撥液性被膜の探索方法は、前記撥液性判定の結果に基づいて、前記固気モデル、及び／又は、前記固液界面モデルの改善仕様を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの改善仕様における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記改善仕様の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記改善仕様の前記計算結果の解析を基に、前記改善仕様の撥液性の良否を判定するステップと、必要に応じて繰り返される、前記改善仕様作成ステップ、前記分子動力学法及び／又は前記分子軌道法を用いた計算ステップ、前記計算結果解析ステップ並びに前記撥液性判定ステップからなるステップ群と、を有することが好ましい。
このような構成の撥液性被膜の探索方法は、前記撥液性判定の結果を反映させて、新たに改善仕様を作成し、さらに当該改善仕様についても前記分子動力学法及び前記分子軌道法を用いて計算・解析を行うことによって、モデルの最適化を図り、その結果最適な撥液性被膜を選択することができる。また、本発明の撥液性被膜の探索方法は、必要に応じて繰り返される前記ステップ群を有することによって、前記改善仕様の更なる改善を図ることができる。
本発明の被膜の撥液性を予測するためのシミュレーション方法は、分子動力学法及び／又は分子軌道法に基づいて、被膜の撥液性を予測するためのシミュレーション方法であって、固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、を有することを特徴とする。
このような構成の被膜の撥液性を予測するためのシミュレーション方法は、上述した本発明に係る撥液性被膜の探索方法を具体的に実施することができる。
本発明の被膜の撥液性を予測するためのシミュレーション方法は、前記撥液性判定の結果に基づいて、前記固気モデル、及び／又は、前記固液界面モデルの改善仕様を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの改善仕様における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記改善仕様の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記改善仕様の前記計算結果の解析を基に、前記改善仕様の撥液性の良否を判定するステップと、必要に応じて繰り返される、前記改善仕様作成ステップ、前記分子動力学法及び／又は前記分子軌道法を用いた計算ステップ、前記計算結果解析ステップ並びに前記撥液性判定ステップからなるステップ群と、を有することが好ましい。
このような構成の被膜の撥液性を予測するためのシミュレーション方法は、上述した本発明に係る撥液性被膜の探索方法の好ましい形態を具体的に実施することができる。
本発明の被膜の撥液性を予測するためのシミュレーションプログラムは、分子動力学法及び／又は分子軌道法に基づいて、被膜の撥液性を予測するためのシミュレーションプログラムであって、コンピュータ上で実行されることにより、固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、を実行することを特徴とする。
このような構成の被膜の撥液性を予測するためのシミュレーションプログラムは、上述した本発明に係るシミュレーション方法を具体的に実施することができる。
本発明の被膜の撥液性を予測するためのシミュレーションプログラムは、コンピュータ上において、さらに、前記撥液性判定の結果に基づいて、前記固気モデル、及び／又は、前記固液界面モデルの改善仕様を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの改善仕様における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記改善仕様の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記改善仕様の前記計算結果の解析を基に、前記改善仕様の撥液性の良否を判定するステップと、必要に応じて繰り返される、前記改善仕様作成ステップ、前記分子動力学法及び／又は前記分子軌道法を用いた計算ステップ、前記計算結果解析ステップ並びに前記撥液性判定ステップからなるステップ群と、を実行することが好ましい。
このような構成の被膜の撥液性を予測するためのシミュレーションプログラムは、上述した本発明に係るシミュレーション方法の好ましい形態を具体的に実施することができる。
本発明によれば、前記分子動力学法及び前記分子軌道法の両方を用いることによって、多角的に前記固気モデル及び／又は前記固液界面モデルについて計算・解析し、当該モデルの良否を判定し、さらに、当該判定結果を基に最適な撥液性被膜を選択することができる。
本発明の探索方法に用いる固気モデル又は固液界面モデルの典型例の概略断面図である。 本発明の撥液性被膜の探索方法の典型例を示すフローチャートである。 実施例、比較例１及び比較例２の解析用モデルを用いた吸着エネルギーの計算結果をまとめたものである。
本発明の撥液性被膜の探索方法は、分子動力学法及び／又は分子軌道法に基づいた、撥液性被膜の探索方法であって、固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、前記判定に基づいて、前記固気モデル及び／又は前記固液界面モデルを最適な撥液性被膜を有するモデルとして選択するステップと、を有することを特徴とする。
上記探索方法を具体的に実施するための、本発明の被膜の撥液性を予測するためのシミュレーション方法は、分子動力学法及び／又は分子軌道法に基づいて、被膜の撥液性を予測するためのシミュレーション方法であって、固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、を有することを特徴とする。
上記シミュレーション方法を具体的に実施するための、本発明の被膜の撥液性を予測するためのシミュレーションプログラムは、分子動力学法及び／又は分子軌道法に基づいて、被膜の撥液性を予測するためのシミュレーションプログラムであって、コンピュータ上で実行されることにより、固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、を実行することを特徴とする。
本発明においては、「撥液性」とは、液体をはじく性質のことをいい、その具体例としては、撥水性及び撥油性が挙げられる。
実験済みの材料以外の材料を利用したいが、当該技術分野においては、実験上のノウハウが存在しないという場合には、最適な材料を確実に選択することが不可能であった。例えば、基板上に被膜を作製し、当該被膜の撥水性・撥油性等を調べる場合には、当該被膜に最適な材料を選び出すことは難しかった。特に、新規材料で被膜を作製しようとする場合は、実験計画法等に基づいた評価計画の立案、データ出力及び評価が必要となり、莫大な試験工数をかけなければならなかった。
また、従来技術においては、上記実験をたとえ経たとしても、新規材料の撥水性・撥油性に対して○又は×というような大ざっぱな定性的評価レベルでの選択には役立つが、定量的な評価レベルでの選択に役立つほどの情報のフィードバックができなかった。このことにより、実証実験によって、最適な選択であったと評価を下すことができるような新規材料を得るまでに、試行錯誤しながら実験を実施しなければならず、上記同様の莫大な試験工数をかけなければならなかった。
本発明は、このような試行錯誤により生じる無駄を省き、多角的に固気モデル及び／又は固液界面モデル等の実験モデルを解析し、当該モデルの良否を判定し、さらに、当該判定結果を基に最適な撥液性被膜を選択することができる、分子動力学法及び／又は分子軌道法に基づいた、撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法、及びコンピュータ上でシミュレーションを行うシミュレーションプログラムに関する。
本発明の前提としては、分子動力学法及び／又は分子軌道法により計算・解析することが可能な実験モデル（後述するような固気モデル、固液界面モデル等）を作成することができるソフトウェアを準備することが必要である。
このようなソフトウェアとしては、具体的には、分子動力学法においては、Ｍａｔｅｒｉａｌｓ Ｅｘｐｌｏｒｅｒ（商品名。富士通株式会社製）、Ｍａｔｅｒｉａｌｓ Ｓｔｕｄｉｏ Ｄｉｓｃｏｖｅｒ、Ｍａｔｅｒｉａｌｓ Ｓｔｕｄｉｏ Ｆｏｒｃｉｔｅ（いずれも商品名。アクセルリス株式会社製）等が挙げられる。
このようなソフトウェアは、作成したモデルを分子動力学法及び／又は分子軌道法による計算、及び当該計算結果の解析においても用いることができる。すなわち、当該ソフトウェアは、モデルの作成・計算・解析のためのプラットフォームとすることができる。
以下、本発明の探索方法に用いる固気モデル又は固液界面モデルの典型例について詳細に述べる。図１は、本発明の探索方法に用いる固気モデル又は固液界面モデルの典型例の概略断面図である。
本発明においては、特に、エンジン内部に付着したエンジンデポジットを低減させるための最適なエンジン内機構開発をその主要な課題としているため、固液界面モデル及び固気界面モデルの作成においては、エンジンの内壁面の微細構造や、エンジン作動時の構造変化等を考慮する必要がある。
固気モデル又は固液界面モデル１００は、固体基質１、被膜２、液体分子又は気体分子３からなる。被膜２は、固体基質１の一部又は全部に接触している複数の有機系分子２ａから構成されている。
解析用モデルの作成方法としては、具体的には、市販品（ソフトウェアメーカー等により販売されているもの等）の原子／分子モデル作成ソフトによりモデル作成する方法が例示できる。このような市販品のソフトにおいては、検討する材料の原子の種類・位置・原子間の結合情報等を、該ソフトが指定するフォーマットに従って入力することで、望みの解析用モデルが作成可能である。
このような固気モデル又は固液界面モデルに用いられる固体基質の例としては、マルテンサイト系ステンレス、Ｃｒ２Ｏ３等の無機固体基質等が挙げられる。
また、このような固気モデル又は固液界面モデルに用いられる表面被膜の例としては、フルオロアルキルシラン膜等の有機化合物被膜等が挙げられる。
また、このような固液界面モデルに用いられる液体の例としては、エンジンデポジット中に含まれる物質、具体的には、トリメチルベンゼン等の有機化合物等が挙げられる。
分子動力学法とは、多粒子系、すなわち原子や分子の集団の構造や動的挙動を調べるために、粒子間の相互作用が既知であるとして、計算機によって古典的なニュートンの運動方程式を数値的に解く方法のことである（化学事典（第一版、丸善）より）。分子動力学法においては、一般に巨視的な原子、分子集団のモデルとして数百から数万個の粒子系を選び、互いに相互作用する粒子の運動を１〜数十ｆｓ間隔で追跡して各粒子の運動方程式を解き、位置と運動量を決定する。運動方程式を解くためのアルゴリズムとしては、Ｖｅｒｌｅｔ法、Ｌｅａｐ−ｆｌｏｐ法、ＳＨＡＫＥ法、ＲＡＴＴＬＥ法、ＲＯＬＬ法、ＲＥＳＰＡ法、Ｔｒｅｅ法、カットオフ法、Ｅｗａｌｄ法、ＰＭＥ（Ｐｅｒｔｃｌｅ Ｍｅｓｈ Ｅｗａｌｄ）法、ＦＭＭ法等がある。
分子軌道法とは、分子内の電子配置やエネルギー準位を研究する理論的方法のことである（化学事典（第一版、丸善）より）。分子軌道法は、その具体的手法としては、分子全体に広がった軌道（分子軌道）に電子を詰めていき、その軌道に対する波動方程式を解くことによって、分子内の電子配置やエネルギー準位を決定する。
本発明においては、固気モデル及び固液界面モデルを分子動力学法及び／又は分子軌道法により解析して、固体基質の表面に接近しようとする気体又は液体の吸着エネルギーを計算する。
吸着エネルギーは、以下の式により求めることが可能である。
Figure 2010182106
 （ただし、Ｅａｄｓｏｒｐｔｉｏｎは吸着エネルギー、ＥＡｌｌＭｏｄｅｌは基板、被膜及び吸着分子からなる解析用モデルの全エネルギー、ＥＣｏａｔ＋ＢａｓｅＭａｔｅｒｉａｌは基板及び被膜のエネルギー、ＥＭｏｌｅｃｕｌｅは吸着分子のエネルギーである。）
上述した分子動力学法を用いた計算、及び、分子軌道法を用いた計算は、どちらを先に行ってもよく、また、同じ計算を複数回繰り返し行ってもよい。ただし、分子動力学法を用いた計算、又は、分子軌道法を用いた計算のいずれか一方は、少なくとも１回は行うこととする。
ここで、「改善仕様」とは、分子軌道法や分子動力学法などの計算において新たに得られた情報に基づいて改善が期待される仕様のことである。改善仕様は、具体的には、モデルの初期状態からの材料変更や、パラメータ変更等を適宜行うことによって作成する。
以下、本発明の撥液性被膜の探索方法の典型例について詳細に述べる。図２は、本発明の探索方法の典型例を示すフローチャートである。なお、以下は、被膜解析用の固液界面モデルの選択について説明するが、固気モデルの選択も、下記の説明同様に行うことができる。
まず、被膜解析用の固液界面モデルの現状仕様を作成する。ここでいう「現状仕様」とは、既知の情報に基づいて、技術者（撥液性被膜の設計者等）が性能向上を期待できる仕様のことである。なお、現状仕様としては、固液界面モデルの初期状態を用いてもよい。
固液界面モデルとしては、図１に示したような、固体基質、当該固体基質を被覆する表面被膜、及び液体からなる典型例を例示することができる。
次に、現状仕様の吸着エネルギーを、分子動力学法及び分子軌道法を用いて計算する。なお、このときの計算の手順としては、まず、分子動力学法及び分子軌道法を用いて、基板、被膜及び吸着分子からなる解析用モデルの全エネルギーであるＥＡｌｌＭｏｄｅｌ、基板及び被膜のエネルギーであるＥＣｏａｔ＋ＢａｓｅＭａｔｅｒｉａｌ、吸着分子のエネルギーであるＥＭｏｌｅｃｕｌｅをそれぞれ計算し、その後にこれらの計算結果に基づいて、式（１）に示したような吸着エネルギーＥａｄｓｏｒｐｔｉｏｎについて計算するという方法を例示することができる。
このとき、吸着エネルギーＥａｄｓｏｒｐｔｉｏｎが０未満である場合（図２において「Ｙｅｓ」の場合；すなわち、固体基質への液体分子の反発が起こる場合）には、当該モデルが良い撥液性を示すことを判定したのち、当該モデルを最適な撥液性被膜を有するモデルとして選択し、必要であれば、当該モデルを固液界面モデルとして用いた固体基質、表面被膜及び液体を用いた実験の準備を始めることができる。
なお、当該吸着エネルギーが０以上である場合（図２において「Ｎｏ」の場合；すなわち、固体基質への液体分子の吸着が起こる場合）には、解析結果を反映させた上で、新たに被膜解析用の固液界面モデルの改善仕様を作成する。
新たに作成された改善仕様の吸着エネルギーＥａｄｓｏｒｐｔｉｏｎを、再び分子動力学法及び分子軌道法を用いて計算する。このとき、当該吸着エネルギーが０未満である場合（図２において「Ｙｅｓ」の場合）には、当該モデルが良い撥液性を示すことを判定したのち、当該モデルを最適な撥液性被膜を有するモデルとして選択し、必要であれば、当該モデルを固液界面モデルとして用いた固体基質、表面被膜及び液体を用いた実験の準備を始めることができる。なお、この判定の際に、既に計算した現状仕様の吸着エネルギーと、改善仕様の吸着エネルギーとを比較することにより、撥液性の良否の判定に反映させてもよい。なお、吸着エネルギーが０以上である場合（図２において「Ｎｏ」の場合）には、解析結果を反映させた上で、再び改善仕様を作成する。
本発明によれば、分子動力学法及び分子軌道法の両方を用いることによって、多角的に固気モデル及び／又は固液界面モデルを用いて計算・解析し、当該モデルの良否を判定し、さらに、当該判定結果を基に最適な撥液性被膜を選択することができる。
本発明のシミュレーション方法を用いて、解析用モデルを作成し、計算検討を行った。
１．解析用モデルの作成
［実施例］
本実施例の解析用モデルは、固体基質として、マルテンサイト系ステンレスの一種である、ＳＵＳ４４０Ｃ、表面被膜としてフルオロアルキルシラン膜、液体としてトリメチルベンゼンを含んだエンジンデポジットを有する解析用モデルである。
エンジンデポジットは、当該エンジンデポジット中に含まれている原子及び分子を分析して、可能な限り近似したクラスタモデルを作成した。実際の分析結果より、エンジンデポジット中にはトリメチルベンゼンを確認したため、トリメチルベンゼンを含んだクラスタモデルを作成した。
また、フルオロアルキルシラン膜については、当該膜中に含まれているフルオロアルキルシラン分子を分析して、可能な限り近似したクラスタモデルを作成した。
さらに、固体基質に関しては、目的とする結晶モデルを厳密に作成することは困難であったため、元素組成が近いＣｒ２Ｏ３の結晶をモデルとして代用した。
［比較例１］
比較例１の解析用モデルは、固体基質として、マルテンサイト系ステンレスの一種である、ＳＵＳ４４０Ｃ、液体としてトリメチルベンゼンを含んだエンジンデポジットを有する解析用モデルである。すなわち、実施例の解析用モデルのような表面被膜は、比較例１の解析用モデルには含まれない。
エンジンデポジット及び固体基質に関しては、実施例と同様にモデルを作成した。
［比較例２］
比較例２の解析用モデルは、固体基質として、マルテンサイト系ステンレスの一種である、ＳＵＳ４４０Ｃ、表面被膜として金薄膜、液体としてトリメチルベンゼンを含んだエンジンデポジットを有する解析用モデルである。
エンジンデポジット及び固体基質に関しては、実施例と同様にモデルを作成した。
２．解析用モデルを用いた吸着エネルギーの計算
分子軌道法を利用した分子シミュレーションプログラムを用いて、実施例、比較例１及び比較例２の解析用モデルについて、上記式（１）の吸着エネルギーを計算した。なお、分子軌道法としては密度汎関数法（ＤＦＴ）を用いた。
図３は、実施例、比較例１及び比較例２の解析用モデルを用いた吸着エネルギーの計算結果をまとめたものである。
図３に示すように、計算結果から、比較例１の解析用モデルについての吸着エネルギー値は正の値を算出した。したがって、エンジンデポジット中に含まれるトリメチルベンゼンが、固体基質であるＳＵＳ４４０Ｃに対して吸着しやすい成分であることがシミュレーションできた。
また、計算結果から、比較例２の解析用モデルについての吸着エネルギー値は正の値を算出した。したがって、固体基質（ＳＵＳ４４０Ｃ）を金薄膜によって被覆しても、エンジンデポジット中に含まれるトリメチルベンゼンは固体基質に対して依然吸着することがシミュレーションできた。
しかし、計算結果から、実施例の解析用モデルについての吸着エネルギー値は負の値を算出した。したがって、固体基質（ＳＵＳ４４０Ｃ）をフルオロアルキルシラン膜によって被覆することで、フルオロアルキルシラン膜がエンジンデポジット中に含まれるトリメチルベンゼンに対して反発し、その結果、トリメチルベンゼンが固体基質に対して吸着し難くなることがシミュレーションできた。
以上より、実施例、比較例１及び比較例２の解析用モデルを用いた吸着エネルギーの計算結果を比較することにより、実施例に用いたフルオロアルキルシラン膜は、被膜が存在しない場合（比較例１）、又は金薄膜を用いた場合（比較例２）と比較して、エンジンデポジットをより吸着し難くする効果があるというシミュレーション結果が得られた。
１ 固体基質
２ 被膜
２ａ 有機系分子
３ 液体分子又は気体分子
１００ 固気モデル又は固液界面モデル

Claims (6)

  1. 分子動力学法及び／又は分子軌道法に基づいた、撥液性被膜の探索方法であって、
    固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、
    前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、
    前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、
    前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、
    前記判定に基づいて、前記固気モデル及び／又は前記固液界面モデルを最適な撥液性被膜を有するモデルとして選択するステップと、を有することを特徴とする、撥液性被膜の探索方法。
  2. 前記撥液性判定の結果に基づいて、前記固気モデル、及び／又は、前記固液界面モデルの改善仕様を作成するステップと、
    前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの改善仕様における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、
    前記改善仕様の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、
    前記改善仕様の前記計算結果の解析を基に、前記改善仕様の撥液性の良否を判定するステップと、
    必要に応じて繰り返される、前記改善仕様作成ステップ、前記分子動力学法及び／又は前記分子軌道法を用いた計算ステップ、前記計算結果解析ステップ並びに前記撥液性判定ステップからなるステップ群と、を有する、請求項１に記載の撥液性被膜の探索方法。
  3. 分子動力学法及び／又は分子軌道法に基づいて、被膜の撥液性を予測するためのシミュレーション方法であって、
    固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、
    前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、
    前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、
    前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、を有することを特徴とする、被膜の撥液性を予測するためのシミュレーション方法。
  4. 前記撥液性判定の結果に基づいて、前記固気モデル、及び／又は、前記固液界面モデルの改善仕様を作成するステップと、
    前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの改善仕様における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、
    前記改善仕様の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、
    前記改善仕様の前記計算結果の解析を基に、前記改善仕様の撥液性の良否を判定するステップと、
    必要に応じて繰り返される、前記改善仕様作成ステップ、前記分子動力学法及び／又は前記分子軌道法を用いた計算ステップ、前記計算結果解析ステップ並びに前記撥液性判定ステップからなるステップ群と、を有する、請求項３に記載の被膜の撥液性を予測するためのシミュレーション方法。
  5. 分子動力学法及び／又は分子軌道法に基づいて、被膜の撥液性を予測するためのシミュレーションプログラムであって、コンピュータ上で実行されることにより、
    固体基質と気体分子を有する固気モデル、及び／又は、固体基質と液体分子を有する固液界面モデルの初期状態を作成するステップと、
    前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの初期状態における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、
    前記初期状態の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、
    前記計算結果の解析を基に、前記固気モデル及び／又は前記固液界面モデルの撥液性の良否を判定するステップと、を実行することを特徴とする、被膜の撥液性を予測するためのシミュレーションプログラム。
  6. コンピュータ上において、さらに、
    前記撥液性判定の結果に基づいて、前記固気モデル、及び／又は、前記固液界面モデルの改善仕様を作成するステップと、
    前記分子動力学法及び／又は前記分子軌道法を用いて、前記固気モデル及び／又は前記固液界面モデルの改善仕様における、前記固体基質への、前記気体分子及び前記液体分子の少なくともいずれか一方の吸着エネルギーを計算するステップと、
    前記改善仕様の、前記分子動力学法及び／又は前記分子軌道法の計算結果を解析するステップと、
    前記改善仕様の前記計算結果の解析を基に、前記改善仕様の撥液性の良否を判定するステップと、
    必要に応じて繰り返される、前記改善仕様作成ステップ、前記分子動力学法及び／又は前記分子軌道法を用いた計算ステップ、前記計算結果解析ステップ並びに前記撥液性判定ステップからなるステップ群と、を実行する、請求項５に記載の被膜の撥液性を予測するためのシミュレーションプログラム。
JP2009025148A 2009-02-05 2009-02-05 撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム Pending JP2010182106A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009025148A JP2010182106A (ja) 2009-02-05 2009-02-05 撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009025148A JP2010182106A (ja) 2009-02-05 2009-02-05 撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム

Publications (1)

Publication Number Publication Date
JP2010182106A true JP2010182106A (ja) 2010-08-19

Family

ID=42763659

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009025148A Pending JP2010182106A (ja) 2009-02-05 2009-02-05 撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム

Country Status (1)

Country Link
JP (1) JP2010182106A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109448798A (zh) * 2018-11-30 2019-03-08 湘潭大学 一种分子动力学模拟二元掺杂离子液体润滑薄膜结构特性的方法
JP2021504831A (ja) * 2017-11-30 2021-02-15 １キュービー インフォメーション テクノロジーズ インコーポレイテッド１Ｑｂ Ｉｎｆｏｒｍａｔｉｏｎ Ｔｅｃｈｎｏｌｏｇｉｅｓ Ｉｎｃ． 量子古典コンピューティングハードウェア用いた量子コンピューティング対応の第一原理分子シミュレーションのための方法とシステム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JPN6013036365; 中島章: '中島ナノウエッティングプロジェクト' インターネット , 2005 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2021504831A (ja) * 2017-11-30 2021-02-15 １キュービー インフォメーション テクノロジーズ インコーポレイテッド１Ｑｂ Ｉｎｆｏｒｍａｔｉｏｎ Ｔｅｃｈｎｏｌｏｇｉｅｓ Ｉｎｃ． 量子古典コンピューティングハードウェア用いた量子コンピューティング対応の第一原理分子シミュレーションのための方法とシステム
JP7312173B2 (ja) 2017-11-30 2023-07-20 グッド ケミストリー インコーポレイテッド 量子古典コンピューティングハードウェア用いた量子コンピューティング対応の第一原理分子シミュレーションのための方法とシステム
CN109448798A (zh) * 2018-11-30 2019-03-08 湘潭大学 一种分子动力学模拟二元掺杂离子液体润滑薄膜结构特性的方法
CN109448798B (zh) * 2018-11-30 2021-09-03 湘潭大学 一种分子动力学模拟二元掺杂离子液体润滑薄膜结构特性的方法

Similar Documents

Publication Publication Date Title
Saliba et al. Comparison of gasoline direct-injection (GDI) and port fuel injection (PFI) vehicle emissions: emission certification standards, cold-start, secondary organic aerosol formation potential, and potential climate impacts
Raza et al. A review of particulate number (PN) emissions from gasoline direct injection (GDI) engines and their control techniques
Zhao et al. Intermediate volatility organic compound emissions from on-road gasoline vehicles and small off-road gasoline engines
Zhao et al. Intermediate volatility organic compound emissions from on-road diesel vehicles: chemical composition, emission factors, and estimated secondary organic aerosol production
Mehl et al. Detailed chemistry promotes understanding of octane numbers and gasoline sensitivity
Naik et al. Applying detailed kinetics to realistic engine simulation: The surrogate blend optimizer and mechanism reduction strategies
Yao et al. A hybrid mechanism for n-dodecane combustion with optimized low-temperature chemistry
Zhang et al. A computational investigation of fuel chemical and physical properties effects on gasoline compression ignition in a heavy-duty diesel engine
Battistoni et al. Combustion CFD modeling of a spark ignited optical access engine fueled with gasoline and ethanol
JP2010182106A (ja) 撥液性被膜の探索方法、被膜の撥液性を予測するためのシミュレーション方法及びシミュレーションプログラム
Zheng et al. Trends of full-volatility organic emissions in China from 2005 to 2019 and their organic aerosol formation potentials
Xiao et al. Kinetic modeling study of the ignition process of homogeneous charge compression ignition engine fueled with three-component diesel surrogate
Vella et al. Prediction of the distillation curve and vapor pressure of alcohol–gasoline blends using pseudocomponents and an equation of state
Zhao et al. Comprehensive assessment for the Impacts of S/IVOC emissions from mobile sources on SOA formation in China
Fleitmann et al. Molecular Design of Fuels for Maximum Spark-Ignition Engine Efficiency by Combining Predictive Thermodynamics and Machine Learning
Osipova et al. Effect of addition of methyl hexanoate and ethyl pentanoate on the structure of premixed n-heptane/toluene/O2/Ar flame
Maio et al. ECFM-LES modeling with AMR for the CCV prediction and analysis in lean-burn engines
Kamp et al. Perspectives on current and future requirements of advanced analytical and characterization methods in the automotive emissions control industry
St. John et al. Development of a Data-Derived Sooting Index Including Oxygen-Containing Fuel Components
Hakeem et al. Characterization and speciation of fuel oil dilution in gasoline direct injection (di) engines
Jatoth et al. Impact of thermophoresis factor on soot particle trajectories near the in-cylinder wall in a diesel engine
Ismail et al. Application of Adaptive Local Mesh Refinement (ALMR) Approach for the Modeling of Reacting Biodiesel Fuel Spray Using OpenFOAM
Mariani et al. Neural network-based prediction of liquid-phase diffusion coefficient to model fuel-oil dilution on engine cylinder walls
Gao et al. Steady-state Calibration of a Diesel Engine in CFD using a GPU-based Chemistry Solver
Puduppakkam et al. Predicting emissions using CFD simulations of an e30 gasoline surrogate in an HCCI engine with detailed chemical kinetics

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20111129

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20120703

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130723

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20131126