∞ Joumalo ft h eCombustionS o c i 巴t yo fJ a p a n . l47No.139( 2 5 )5 8 6 3 Vo 日本燃焼学会誌第 47巻 139号 ( 2 5年) 5 8 6 3 ∞ -原著論文/ORIG町 ALPAPER. 火花点火希薄燃焼の着火限界のシミュレーション NumericalC a l c u l a t i o no ft h eI g n i t i o nL im i tonLean-Burn8 1Engines 太田晴也 1* ・橋本公太郎 2 ・新井充 2 ・田村昌三 2 * l 2 I T A， Haruya ， HASHIMOTO， Kohtaro， ARAI， Mitsuru2， andTAMURAMasamitsu3 Ol l ∞ 東京ガス縁式会社フロンティア研究所〒2 3 0 - 45横浜市鶴昆区末広町1 7 7 TokyoGasCo L t d . .1 7 7 .S u e h i r o c h o .T s附・附n i k u .Yok οhama.Kallagawa.2300 0 4 5 .Japan 叫 2 東京大学大学院新領域創成科学研究科環境学専攻〒1/3-0077文京区本郷7ふ I TokyoU l l i v e r s i t y .7 ふ1 .H O l l g o .B l l l z k y o・k u .Tokyo1 1 3心0 3 3 .Japan ∞ ∞ ∞ 3日受理/ R e c e i v e d6J u l y .2004;A c c e p t e d1 3J a n u a r y2 5 2 4年 7月 6日受付 ;2 5年 1月 1 g n i t i o np r o b a b i l i t yo ft h elean-bumS Ie n g i n ei nt h ev i c i n i t yo ft h el e a n b u ml i m i twas A b s t r a c t:A modelp r e d i c t i n gt h巴 i c o n s t r u c t e du s i n gt h ep e r c o l a t i o n model i nc o n t i n u o u ss p a c e .C a l c u l a t e dp e r c o l a t i o np r o b a b i l i t i e sr e p r o d u c e dt h e e x p e r i m e n t a lr e s u l t swhent h ef o l l o w i n gc o n d i t i o n swerea p p l i e d .c o n s i d e r i n gt h ea c t u a li g n i t i o nonS Ie n g i n e .F i r s t .t h eh e a t g e n e r a t i o no ft h es p a r kp l u gs h o u l dbel a r g e r( 1 5t o20t i m e s )t h a nt h a to ft h ee x o t h e r r n i c陀 a c t i o no ft h ef u e lp a r t i c l e.S e c o n d . sr e s u l t t h es c a l eo ft h ep e r c o l a t i o ns h o u l dbemoret h a nt h et w i c eo ft h er a n g eo ft h eh e a tg e n e r a t i o noft h es p a r kp l u g.Thi s u g g e s t e dt h a tt h es i z eo ft h ef l a m ec o r ewouldh a v ea b o u tt h et w i c es i z ea st h ee f f e c t i v er a n g eo fh e a t i n gbyt h es p a r kp l u g KeyWords:S le n g i n e，1 9ni t i o nl i r n i t .N u m e r i c a lC a l c u l a t i o n， Combustion 値計算による，火炎核成長のモデル化の研究もなされてい 1. 緒言る[ 5 ] . これらの数値モデルでは，初期条件によって着火す希薄燃焼とは，燃料混合気が化学量論比よりも希薄な条るか否かが決定する. しかしながら，実際の火花点火エン件において燃焼することである.希薄燃焼の長所は， 1サジンの希薄限界近傍では，ある希薄限界で着火確率が lかイクルあたりの燃料の噴射量が少ないために低燃費であるら 0に変化するのではなく，ある濃度の幅をもって lからこと，燃焼温度が低いことから NO xの生成量も抑制される Oに変化する.従って，着火限界近傍の着火確率を求めることが挙げられる.このため，自動車の内燃機関に希薄燃ためには，決定論的モデルではなく，確率モデルを用いる焼を用いるリーンパーンエンジンの開発が進められてい必要がある.また，一般的なエンジンのシミュレーシヨンる.しかし，希薄燃焼は燃焼限界に近いために制御が難しでは，化学反応計算や流体力学計算により，化学反応や熱し不完全な燃焼により未燃炭化水素 (THC)やすす (soot) 力学過程による発熱や吸熱を表現する. しかし燃料分子のの排出が増加するため，燃焼限界についての理解を深める酸化反応には，多くの化学反応が関与する.燃料分子に含ことは，リーンパーンエンジンの開発に貢献するものと期まれる炭素数が増えると，反応式の数は指数関数的に増加待される. するため，エンジンで用いられる燃料の酸化反応を網羅す火花点火エンジンの希薄限界に関する研究は，数多くなることは困難である. そこで，本研究では，より簡便に希されてきている [1-3]. 火花点火エンジンの着火機構につい薄燃焼の着火限界をシミュレートする確率モデルを考ても研究がなされており，点、火プラグによる火花が火炎核えた. を形成し，この火炎核が乱流場で燃焼反応熱と未燃混合気希薄燃焼における着火限界付近での着火確率を再現するや点火プラグへの熱損失のバランスで成長するかどうかためには，燃料粒子が酸化反応を起こして発熱し，その熱で，着火が決まると考えられている [ 4 ] . そして，詳細な数が他の燃料粒子に伝わり新たに酸化反応を起こす，反応の連鎖を表現することが重要である *Correspondingauthor .E m a i l :h a r u y a @ t o k y o g a s . c o . j p rもののつながり方」を表現するモデルとしては，パーコレーションモデルがよく知られている.パーコレーシヨンモデルは，「つながり」 ( 5 8 ) 太田晴也ほか:火花点火希薄燃焼の着火限界のシミュレーション 59 という一般的な概念を取り扱うため，素粒子や磁石，金属，高分子，地震発生のメカニズムや渦状星雲の生成まで，広 6， 7 ] . 太田らは格子状ノ号ーコレい分野に応用が可能である [ ーションモデルを用いて燃焼発生機構の考察を行っている [ 8 ] . 本研究では，燃料粒子の位置は任意として，希薄燃焼における着火の連鎖をパーコレーション理論に則して調べた.ここで，パーコレーションモデルを用いたのは，火炎核が成長するか否かを確率論的に表すためであり，混合気の不均一性を議論する訳ではないことを記す. F i g . lC o n n e c t i o no fi g n i t i o nr e a c t i o n s 2 . 計算方法および実験装置・実験方法 2 . 1 . 計算方法火花点火エンジンではスパークプラグにおける電位差か逆に，加熱しても次のステップで酸化反応を起こす燃料粒ら生じるスパーク一点のみが反応の元となり，圧縮着火エ子がない場合には，失火と判定した.以上の手続きを，燃ンジンのように複数の燃料粒子が同時に反応を起こすので料濃度を変化させて行うことで，燃料濃度と着火確率とのはない.放電による熱エネルギーは周りの燃料粒子を加熱関係を求めた.ある燃料濃度における着火確率は，燃料粒し，一定温度を超えた燃料粒子は反応を起こし発熱する. 子の初期配置をランダムに設定した 2 0 0パターンに対するこの発熱が他の燃料粒子を加熱する，という繰り返しによ計算結果から算出した. り発熱反応が伝播し，着火の元となる火炎核を形成するも通常のパーコレーションモデルとの関係で考えると，比例定数 k は連続空間で、のパーコレーションモデルにおけのと考えられる. 作成したモデルでは，連続的な空間に燃料粒子を設定した個数ランダムに分布させるる，つながりの半径 rに対応する粒子の初期温度は 300K とした.次に連続空間の中心に配置したスパークプラグにおここでは k=7 ∞とすることで，つながりの範囲が半径 lの球である連続空間ノ号ーコレーションモデルとほぼ同じものとなる.パーコレーいて，発熱を起こす.発熱は，その時点において，すべてシヨンは，燃料粒子の発熱により，他の燃料粒子が酸化反の燃料粒子が温度上昇を起こすことで表現する.この研究応を起こすとし寸反応の連鎖に相当する.つまり，火炎核では，加熱の影響の伝播則を式 ( 2 )の形に仮定して考えた. の形成は「連鎖が起きる半径 r内に次の燃料粒子が存在する」という連鎖によって，スパークプラグから距離 L の地 L 1 T i= kIf i / 点まで到達することによって表現される(図 1 ) .これはパー ( 2 ) コレーシヨンモデルでのクラスターの形成と対応する. L 1 T iは，一回のステップにおける i番目の粒子の温度上昇， j=Oとした場合の通常のパーコレーションモデルと異なる点は，つながり f i jは，スパークプラグ(j=0 )と i番目のの範囲の外にある燃料粒子に対しでも，粒子聞の距離に依粒子との距離. kは比例定数である. kが大きいほど粒子存する影響(温度上昇)を与えることである.燃料粒子の熱の温度上昇が大きく. kが小さいほど温度上昇が小さいこの伝達では，ある距離以上離れたとしても，熱が全く伝わとを示す.つまり k は発熱の大きさを表すパラメータでらなくなることはないためである.従って，発熱した粒子ある. 温度上昇により叩0 0K を超えた燃料粒子は，次のステからの距離が r よりわずかに遠い場所にあった粒子でも，再度わずかな温度上昇を受けることにより，酸化反応を起ップで発熱反応を起こすものとする.この発熱によって他こすことになる.もう一つの相違点は，火花着火エンジンの燃料粒子も式 ( 2 ) と同様の加熱を受けることになるが，でのスパークプラグにおける発熱を，燃料粒子 1つの反応スパークプラグの放電による加の方が，燃料粒子の発熱反による発熱よりも大きいと設定することである.これは，すなわちスパークプラグの発スパークプラグの発熱によるつながりの半径 R を，燃料粒熱における比例定数 kを k，.燃料の発熱における比例定数応による加熱よりも大きい子の発熱によるつながりの半径 rよりも大きくすること， kを k 2とおくと. k，>k 2とすることで発熱量の大小が表現すなわちスパークプラグにおける比例定数 k，を燃料の発 ∞OKを超えた燃料粒子が複数ある場合は，前のされる. 1 熱における比例定数 k 2よりも大きくすることで表現さステップで発熱した粒子から最も近い燃料粒子が発熱反応れる . シミュレーションを行うにあたって決定すべきパラメーを起こす. この繰り返しによって燃料粒子の酸化反応は伝播し，スタは 2つある. 一つは火炎核の大きさの設定で，スパークノ号ークプラグから一定距離 L 以上離れた燃料粒子が反応をプラグからどれだけ離れた距離にある粒子が反応を起こし起こしたら火炎核が形成したと見なし，着火と判定した. たら着火と見なすか，という距離 L によって表される.も ( 5 9 ) ∞ 第47巻 1 3 9号 ( 2 5年) 日本燃焼学会誌 60 Ta b l e1 E n g i n eo p e r a t i n gc o n d i t i o n う一つは，スパークプラグによる発熱が，粒子一つの発熱に比べてどれだけ大きな規模であるかという比である.こ F u e l iso.octane EngineSpeed 1200rpm される.今回は k 2=700とすることで粒子の発熱によるつ F u e !I n j e c t i o nPressure 12MPa ながりの半径 rを l とし，スパークプラグの発熱によるつ Thr o t t l eOpening 100% I g n i t i o n百 mmg 18 BTDC SwirlValveOpening 0，100 れは，スパークプラグでの発熱でのつながりの半径 R と粒子の発熱におけるつながりの半径 rとの比 R l rによって表ながりの半径 R によって表すことにした.実際の計算においては，発熱における比例定数 k，と k 2の比によって決定したこの 2つのパラメータ Lと Rを変化させた場合の影 0 響から，実際の着火現象を左右する要素が何かを調べる. 2 . 2 . 実験方法火花点火エンジンを用いて均一燃焼での燃料希薄条件において運転を行い，燃料の濃度と着火確率の関係を測定し I n t a k eAir 一一ーでi=-- た.測定に用いたエンジンはトヨタ D・4エンジンをリカルド社が単気筒エンジンに改造したものである.装置図を図 2，運転条件を表 lに示す.エンジンは均一燃焼条件で運転し，空気過剰j 率 λを変えることで燃料濃度を変化させた. 希薄着火限界付近の燃料濃度において運転を行い，測定された失火率から着火確率を調べた.ある空気過剰率毎に 100サイクル運転を行い，圧力上昇が 100kPa以下であったサイクルを失火とみなし，失火したサイクル数から失火率を算出した. 3 . 結果および考察 3 . 1 . 実験による燃料濃度と着火確率の関係エンジンでの測定によって得られた失火率から着火確率を算出した結果を図 3に示す.横軸は空気過剰率入の逆数 F i g . 2 Exp 巴r i m e n t a le n g i n 巴組dm e a s u r e m e n ta p p訂 . a t u s e s で，これは燃料の濃度に比例する.空気過剰率 λは量論混 • 。合比の空気量を l とした場合の空気量であり，燃料の体積 swirll00 swirlO w i r l100がスワール弁の開度を 100 は小さく無視できる. s (全開)とした場合， s w i r l0がスワール開度を o(全閉)とし e e ・ 8. た場合である.スワール弁を聞くと燃焼室内の空気の流れができるため，燃料噴霧の混合が促進される.従って，〉、 s w i r l100は燃料粒子の分布が比較的均一な場合， swirlOは。。分布が不均一な部分が多い場合を表すと考えられる. 結果から，着火確率はある着火限界濃度において突然 0 206 0 から lになるのではなく，ある濃度範囲で徐々に上昇して園。 E0.4 + ' • 旬いることがわかる.これは，希薄着火限界に近い濃度での着火現象が，濃度により決定論的に決まっているのではな。・ • 0 0.8 + ' c く，確率現象であることを示している.着火確率曲線が 0 から立ち上がる部分を比較すると，スワール弁を開いた場ー0.2 bD 。。合では燃料濃度が 0 . 4付近で立ち上がるのに対し，閉じた場合には 0 . 3付近で立ち上がっている.また前者は着火確率のカープが比較的なだらかだが，後者の方が角度が急で O~ 識位 0.2 0. 4 0 . 6 0.8 1 /λ ある.結果として，同じ濃度における着火確率を比較した場合には，スワールがない方が着火確率が高い.これはス 3E x p e r i r r 即日a li g n i t i o na b i l i t yf o rd e f e r e n tf u e lc o n c e n t r a t i o n F i g. ワールがない場合に，部分的に燃料濃度が高い部分に着火するケースが出てくるためだと考えられる. ( 6 0 ) 太田晴也ほか:火花点火希薄燃焼の着火限界のシミュレーション 6 1 総和を想定している空間の体積で割った値であり，図4 (および図 1 1)の横軸の濃度 p と同義である. ーコレーシヨンモデルにおいては，無限に大きなクラノf スターが存在し始める臨界濃度，つまり浸透関値は次元数のみによって決まり，二次元の場合には面積分率 0. 4 5，三 . 1 6のときであることが知ら次元の場合では体積分率が約 0 P ( p ) れている [ 2 ] . 図 4 と図 5の比較において，二次元と三次元の違いがよく見られる.図 5では，着火確率が Oから立ち上がる部分，つまり大きなクラスターが発生し始める体積分率は 0 . 1と0 . 2の間であり，上記の値と一致している. 以上より， R = 1の場合には，通常のパーコレーションモ。。デルと同じ特徴を持つことが確認された.しかし図 5の結 pc 果は，実験結果と比較した場合に，着火確率の上昇が緩や p かであり，着火確率が lへなかなか到達しないためこのず F i g. 4P e r c o l a t i o np r o b a b i l i t y( P )f o rc o n c e n t r a t i o ni n2 d i m e n t i o n a l g r i d p e r c o l a t i o nm o d e l れを解消する必要がある. 次にパラメータ R を 3 ，5 ，1 0， 1 5，2 0 とし，それぞれの Rの値に対して L を 8 ，1 0， 1 5，2 0，3 0，4 0( R=2 0のみ)と変化させ，体積分率と着火確率の関係を調べた.こ〉、 + ' ••• コ 0 . 6 4コ」伺。 c . . 。 0.4 E 二 z h z a グからの発熱の規模を，燃料粒子一つの発熱よりどの程度大きくするか，ということを意味する.この部分が，火花点火エンジンにおける着火の特徴を表現している.まず R =3の場合について，・・ . •• 、 • 、』 + ' •• • こで L>Rである . Rを大きくすることは，スパークプラ • • •• . . 。 0 . 2 • 0~ 0 . 5 L を 8， 1 0， 1 5，2 0 と変化させた結果を図 6，R =5の場合を図 7，R = 1 0の場合を図 8，R= 1 5の場合を図 9に， R=20の場合を図 1 0に示す. R を lから 3へ増やしたことにより，臨界濃度が約 0 . 1 程度に低くなった.また，着火確率が lへ近づくまでの立ち上がりが格段に早くなった. R =1の場合は，着火確率と臨界濃度が対応していた. しかし， R =3 となったこと 1 . 5 で，着火確率は「半径 lの任意の粒子が大きなクラスター z に含まれる確率」ではなく，「半径 3の球の一部が大きな V o l u m ef r a c t i o n クラスターに含まれる確率」へと変化した.この体積 2 7 倍への変化により，大きなクラスターが存在していれば， F i g . 5I g n i t i o np r o b a b i l i t yf o rv o l u m ef r a c t i o n スパークプラグの発熱がそこに届きやすくなったものと考えられる . Rをさらに大きくしてしベと，着火確率は，大きなクラスターが存在する確率そのものに近づいて ~j く. 3 . 2 . シミュレーションによる燃料濃度と着火確率の関係この系の両端を結ぶクラスターの出現確率を図 1 1 に示す 3 . 2 . 1 . Rが着火確率に及ぼす影響 [ 2 ] . R が大きくなると，系の端まで達する大きなクラスタ R=I の場合は，スパークプラグは特別な点ではなく，他ーがある場合に，そのクラスターにスパークプラグからのの燃料粒子と同等である.この場合は，モデルは単純なパ発熱が届く確率が高くなる.つまり， R を lから大きくすーコレーションモデルとなる.着火の可否は，スパークプ 1に近づいて行くことることで，着火確率は図 4から図 1 ラグが，距離 L より遠くまで続く燃料粒子のクラスターにが予想される.図 3の実験結果を見ると，明らかに図 4よ含まれているかどうかで決まる. R =1では，スパークプ I に近い. これはスパークプラグの発熱が大きいりも図 I ラグと燃料粒子の区別はないため，任意の粒子が巨大なクという，火花点火エンジンの特徴が表れているものと推測ラスターに含まれている確率と同値である. これは，パーされるコレーションモデルでいう浸透確率，すなわち選んだ任意の点が無限の大きさのクラスターに含まれる確率と同義で 3 . 2 . 2 . Lの変化が着火確率に及ぼす影響ある.二次元正方格子サイトにおける浸透確率 P と濃度 p 着火と判定する距離 L を変化させた影響を調べる.パーとの関係 [ 2 ]を図 4に， R=1 ，L=8とした場合の計算結果コレーシヨンモデルにおいては，系の大きさを変えると，を図 5に示す.図 5横軸の体積分率は，燃料粒子の体積の無限の大きさのクラスターが生成する確率の，濃度に対す ( 6 1 ) 6 2 ∞ 第4 7巻 1 3 9号 ( 2 5年) 日本燃焼学会誌ふ守9.\.~ 量得‘畠h q_~う aco 一H 一戸﹄回目 h H一一一﹄伺DO﹂ azo-H 一戸﹄凶日 AO ﹂﹀一三一一門芯ら実' 0 . 8 Xo 令官h 0 . 6 〈ム車白 4 0. × ミ〉。 0 . 2 0 . 1 中 0 . 6 o v 。 4 0. + ' 0 0 v o 0 . 8 可。 + 。 0 . 2 • 0 . 2 0 . 3 4 0. 0 . 5 。。 0 . 6 0 . 0 5 Volume fraction × ﹄凶︼ F 十~十十日口 E 平十口千士取。。 0 . 2 十口 L FH判 0 . 2 正証﹃勺 ~ QUFOA P aco 一判一巳切︻ OLazo 一一三一一一 D 叩﹄ h H一一一心句SO ﹄﹄ Z J nunununυ p~ 0. 4 0 . 1 5 F i g . 8I g n i t i o np r o b a b i l i t yf o rv o 1 u m ef r a c t i o n ..，~ 0 . 6 0 . 1 Volume fraction F i g . 6I g n i t i o np r o b a b i l i t yf o rv o 1 u m ef r a c t i o n 0 . 8 8 十 O 口苦 0 . 0 50 . 10 . 1 50 . 20 . 2 50 . 30 . 3 5 0. 4 . 0 40.060.08 0 . 1 0.120.14 V o l u m ef r a c t i o n Volume fraction F i g . 9I g n i t i o np r o b a b i l i t yf o rv o 1 u m ef r a c t i o n F i g . 7I g n i t i o np r o b a b i l i t yf o rv o l u m ef r a c t i o n る依存性が変化する.パーコレーションモデルにおける，系のサイズよりも燃料の濃度によって左右され，系はある 1 において，系の両端を結ぶクラスターの出現確率の図 1 程度の大きさがあればそれ以上の拡大は不要でらあることが系の大きさが有限である場合(細線)には，着火確率はなだ示唆された. らかに変化していくが，系の大きさが無限になると(太線) 臨界濃度 p cを超えた濃度で、は必ず大きさが無限大のクラ 3 . 2 . 3 . 計算と実験との比較スターが出現する.従って，パラメータ L を変えることは，各条件の計算結果を見ると，「濃度 0から立ち上がり始大きなクラスターの出現確率，つまり着火確率の濃度によめるまでの濃度変化」と，「立ち上がり始めてから着火確る影響が変化するものと予想される.具体的には，濃度変率が lに達するまで、の濃度変化」の比率が R によって変化化に対する着火確率の立ち上がりが急になっていくことがすることがわかる.実験結果では両者の比率が 1:1から 1 考えられる. 5または R= :2の間であるのに対し，計算結果では R=1 2 0の場合にこの範囲内の結果が得られたが， R が 10より，5 ，1 0， 1 5，2 0( 図 6・図 1 0 )それぞれにつしかし R=3 いて Lを変化させた場合において，着火確率が lに近くな小さい場合は， 1 :3以上となった.従って，スパークプラる部分での傾きに多少変化が見られるものの，立ち上がりグの放電による加熱が，燃料の発熱反応による加熱よりもの角度に変化は見られなかった.発熱反応の伝播の可否は，大きいことを表現するためには， R がある程度以上の大き ( 6 2 ) 太田晴也ほか:火花点火希薄燃焼の着火限界のシミュレーション一一 - E﹂ B一 L = 圃・ -nunu- -34・・ム - ~.~ 63 . 4 . 結言 h H一一一心司 SOLaco 一一主 ZM- • 0.8 火花点火エンジンの希薄着火限界における着火現象を，連続空間でのパーコレーシヨンモデルを用いることで，簡 ~ 便にシミュレートした.火花点火エンジンによる着火を表 0.6 現するため，以下の条件を適用したところ，実際のエンジ • 0. 4 ンでの試験結果をよくシミュレートできることがわかった. 1 . スパークプラグの発熱によるつながりの範囲を，燃料 A 0. 2 の発熱によるつながりの範囲よりも大きく設定する ( 1 520倍) 2 . 到達距離 L はスパークプラグの発熱によるつながりのム 0 0 範囲 R の倍以上とする ~ ~ト」 0 . 0 2 0 . 0 4 0 . 0 6 0 . 0 8 0 . 1 前者はスパークプラグの放電による加熱が燃料の反応に 0 . 1 2 よる加熱よりも大きいことを表現し，後者は反応がここま Volume fraction で進めば着火を起こす，という火炎核のサイズがスパーク F i g . l OI g n i t i o np r o b a b i l i t yf o rv o l u m ef r a c t i o n プラグによる加熱の有効範囲の倍程度であることを示して一一 85 戸官。﹄叫53uE。OHO-目白一山ω￡﹄。b s-AU し当る. 謝辞本研究を行うにあたり，パーコレーションモデルについて助言をいただいた東京大学地震研究所小屋口剛博教授，実験に関してご助力頂いた JOMOテクニカルリサーチセンターの方々に謝意を表す. 1 . Hamamoto，Y . . Wakisaka，T .，and O h i g a s h i，S .，P r o c . S i x t e e n t hI n t .C o n g r .FISITA2・5 3( 19 7 6 ) . ハ H u UEES&︿ -A司﹄宮且U References pc 2 . Wakisaka，T .，a n dHamamoto，Y .，B u l l e t i noft h eJSME2 2 : l 282( 19 7 9 ) . p F i g . l lA p p e a r a n c ep r o b a b i l i t yo ft h el a r g ec l u s t e rw h i c hc o n n e c t sb o t h 陀 a (2D) e n d so f t h ea 3 . P e t e r s，B .D .，1 .Mec h .E .C9 2 :6 3( 19 7 9 ) . 4 .A r i c i，0 .，Tabaczynski，R .J .，andA r p a c i，V .S .，Combust. .3 0 :3 1( 19 8 3 ) . S c i .Tech .，andHeywood，J .B .，P r o c .23rdS y m p .( i n t . ) 5 .P i s c h i n g e r，S 0 3 3( 19 9 0 ) C o n b u s t .1 さであることが必要で、ある.計算結果では， Rの変化が， .， I n t r o d u c t i o n t i n 6 . Odagaki， T oP e r c o l a t i o n P h y s i c s ( J a p a n e s e )( 19 9 2 ) . 着火確率が立ち上がる濃度(闘値)に大きな影響を与えてい ∞ るが，これは実際のエンジンにおいてスパークプラグのエ 7 . Odagaki，T .，S c i e n c eo f C o n n e c t i o n( i nJ a p a n e s e )( 2 0 ) . ネルギーが着火確率に大きな影響を与えることと類似性を .，andUmemura，A .，A p p l i c a t i o nofP e r c o l a t i o n 8 . Takamori，S ，見せている.また，今回の計算結果では， R と L が 1:2 ，J .C o m b u s t .S o c .Japan4 5 : 9 5 T h e o r yt oSprayCombustion すなわちスパークプラグによる加熱の有効範囲が，火炎核 1 0 2( 2 0 0 3 ) . のサイズの半分以上の場合に実験結果に近い着火確率曲線が得られた.これは火炎核のサイズは，スパークプラグによる加熱の有効範囲の倍程度であることを示唆している. ここでは，スワールが強い場合と弱い場合の影響と R，L との関係を明確に示すには至らなかった. ( 6 3 )