静電気学会誌，44, 1（2020）26-31

J. Inst. Electrostat. Jpn.

論文

ガソリンエンジン用点火プラグにおける流動により延伸した火花放電の直径に関する考察
傅 建華 ，今村 宰＊，秋濱 一弘 ，山﨑 博司＊＊, 1＊＊＊＊（2019年9月12日受付；2019年12月16日受理）

A Study on the Diameter of Elongated Spark Discharge Chanel in Flow Fieldin the Ignition Plug for Gasoline Engine＊＊＊, 1Chien-hua FU , Osamu IMAMURA＊＊, Kazuhiro AKIHAMA＊＊

and Hiroshi YAMASAKI（Received September 12, 2019; Accepted December 16, 2019）
From the global CO2 reduction trend, improvement in thermal efficiency is required for the powertrain of automobiles. A
promising direction for gasoline engine is lean burn. In the lean burn, since the laminar burning velocity decreases, theturbulent flow is intensified to shorten the burning time. Stable ignition in the high velocity and turbulent flow is required in
combination with the advance ignition timing. At this time, since the discharge path at the spark plug gap is elongated by theflow and it is assumed that the elongation of the discharge path affects the behavior of the ignition to flame propagation, it isimportant to understand the phenomena and to model the effect of such discharge path elongation on the ignition behavior.
From this background, the discharge characteristics at spark plug gap in the flow were investigated using a small wind tunnel.
The deformation of the discharge path in the flow was observed using a high-speed camera, and the current and voltage at thespark plug were measured. The relationship between them was investigated and compared with previous works. Especially inthis report, the relations between the electric current and diameter of discharge path were discussed analytically.

1．はじめにている．流動による放電路の変形は放電路の抵抗値の増

パリ協定による世界的な CO2 の低減の方向性から，内加をもたらし，電流，電流プロファイルが変化する他，放燃機関においても熱効率の向上が求められている．ガソ電路の変形により，点火位置や初期火炎の拳動に影響をリンエンジンにおいては，リーンバーンが熱効率の向上及ぼすことが考えられる．ガソリンエンジンにおける筒策として有望視されている．リーンバーンでは層流燃焼内の燃焼をモデル化する上で，このような放電路の延伸速度が低下するため乱流燃焼速度を高めて燃焼期間を短と再放電を含めた点火のモデリングが必要であり，その縮する必要があるが，空気過剰率の増加に伴って要求さモデル化を目指して，定容容器を主として用いた実験

れる流動が強くなり，あわせて点火時期が進角することや，小型風洞を用いた実験 ，渦流式の燃焼室を用いたから，安定して点火することが困難になってくる．特に

実験

流動がある場合には，放電した放電路が流動によって引においては，Kim の式 と呼ばれる電極間の放電路長をき伸ばされ，場合によっては放電が維持できなくなり，加味した以下のような電圧，電流の関係式が用いられて放電路が短絡したり，再放電するような現象も観察されいることが多い．

2）3）などが近年実施されている．他方で点火モデル

4）5）- 0.32 0.51

Vgc = 40.46lspkis
キーワード：火花放電，着火，プラズマ＊
日本大学大学院生産工学研究科機械工学専攻（〒275-8575 千葉県習志野市泉町 1-2-1）
Department of Mechanical Engineering, Graduate Schoolof Industrial Technology, Nihon University, Chiba, 1-2-1,
Japan＊＊
日本大学生産工学部環境安全工学科（〒275-8575 千葉県習志野市泉町 1-2-1）
Department of Sustainable Engineering, College ofIndustrial Technology, Nihon University, Chiba, 1-2-1,
Japan1
imamura.osamu@nihon-u.ac.jpDOI：https://doi.org/10.34342/iesj.2020.44.1.26

1）

p

（1）

Vgc，lspk，is，p は，各々，ギャップ間の電圧，放電路長さ，電流，圧力を表す．Kim の式は，そもそもは自動車用の点火プラグ周りの流動について調べるために，点火プラグ間の電圧と電流を測定することで点火プラグ周りの流動を見積もるために作られた実験式であるが，今村らの実験の結果 もこの式に沿うものとなっている．他方で2）（1）式については，電圧の増加とともに電流がおよそ 1/3乗で減じていくことに関しての物理的な解釈がない．また（1）式から放電路にどの程度のジュール熱が加わっているかについての示唆は得られるが，放電路の太さの情ガソリンエンジン用点火プラグにおける流動により延伸した火花放電の直径に関する考察（傅 建華ら）

27

報がなくそのエネルギー密度が不明である．エンジン筒されているものと同様なものである．本実験では，10 個内では放電により雰囲気の予混合気が点火し火炎伝播すのコイルに同時に 4［ ms ］の充電を行っている．点火プるが，球状の火炎が外側の未燃部に向かって伝播するにラグはプラグのギャップがちょうど直径 25［ mm ］の流あたり，未燃部を火炎面からの伝熱で予熱する必要があ路の中心になるように設置されている．オリフィス間のるため，火炎面の曲率が大きい時は火炎伝播に関しては距離は 360［ mm ］であり，点火プラグは上流側のオリフ不利な状況となる．放電路を円柱状とし，このプラズマィスから 230［ mm ］の位置に設置されている．また点火の周りから火炎の伝播が始まるとすると，プラズマが細プラグの中心から 18［ mm ］下流の内壁に圧力孔を設け，く曲率が大きい場合に比べて，プラズマが太く曲率が小そこでの圧力をテストセクションの圧力として計測した．

さい場合の方が初期火炎の曲率が小さいと想像され，火電磁弁が開き流動が生じてからテストセクション部の圧炎伝播しやすいと考えらえる．このような観点から放電力が立ち上がるのに 1 秒程度かかるため，電磁弁の開の路の径を点火モデルのパラメータとして加えることは意信号から 1.7 秒後にイグニションコイルに充電し，放電，義があるものと思われる．

燃焼の様子を観察した．電磁弁の開の信号から 2.0 秒後

以上の観点から，本報では小型風洞を用いて，放電路に電磁弁を閉止し，実験を終了している．放電の様子はの延伸の様子について，特に放電路の太さについて着目

高速度カメラ（ Phantom M310）にて，35000 fps でモノクし検討を行ったものである．その概要を以下に報告する．

ロ撮影を行った．解像度は 320 × 240 で視野は 12.6［ mm ］× 9.5［ mm ］程度である．電流，電圧の計測については，

2．実験装置および方法

テクトロニクス社の TCPA300 および P6015A をプローにブとして用い，オシロスコープを用いて計測した．本報

示されるものと同様である．ガスボンベからバッファー

では予混合気としてはプロパン空気の予混合気であり，タンクに混合気が一度貯められた後，電磁弁を開くこと

当量比 φ は 0 から 1 の間で変化させているが，予混合気でバッファータンクから定断面積の流路内に混合ガスがでは放電路からの発光の分布が強く，本報で紹介する結流れ込む．その流路内で火花放電を行い，その様子を観果の多くは空気を作動流体として用いた場合である．テ

察できる実験装置となっている．流路は直径 25［ mm ］

ストセクションの圧力は最大で 2.5［ MPa ］程度である．

実験装置の全体図を図 1 に示す．実験装置は文献

2）の円形断面積であり，点火プラグ前後にオリフィスを挿入することで流量の調整を行っている．火花放電が行わ

3．実験結果および考察れるテストセクションでの平均の流速はオリフィスの面

3.1 電圧および電流計測の結果積が一定であれば圧力によらず，理想的にはオリフィス

図 2 には典型的な電圧および電流波形を示す．雰囲気径を定めることで圧力によらず一定流速下における試験流体の条件は，当量比 φ = 0.8 のプロパン空気予混合気が実施できる．流速の分布および絶対値については，で平均流速が Um = 18.5［ m ⁄ s ］である．放電時の圧力はPIV を用いた計測により検討を行っている ．本報では

図 2 に示されるようにおよそ 0.8［ MPa ］である．時間

流速として，オリフィス径を基に算出された流速 Um，軸の 0 はプラグギャップ間の放電が始まった時刻を 0を用いて記述する．なお本報で示す実験結果は主に Um［ ms ］として取っており，放電時間はおよそ 1［ ms ］程= 18.5［ m ⁄ s ］の場合である．

度である．時刻 0［ ms ］において，電流，電圧とも大き6）

点火プラグは NGK 社製の PSPE プラグ，点火コイルは10 個のコイルを 2 直 5 並列に接続しており，文献

7）に示

く立ち上がっているが，この部分には本報では着目をしていないため時間分解能が十分ではない．その後，電流は緩やかに減少する一方で電圧は急上昇，急降下を繰り返している．これは，放電路の短絡や再放電に起因しており，その時刻で電流値にも変化が見られる．電圧の上昇は放電路が伸長することによって放電路の抵抗値が上昇することに起因している．その後，放電路が短絡もしくは再放電による放電路長さの減少によって，放電路の抵抗値が低下することから，電圧は急激に減少する．電流が 0［ mA ］となったところで放電は終了し，その時刻

図 1 実験装置概略図
Fig.1 Experimental apparatus.

以降，高速度カメラでも放電路は観察されない．それ以降はプラグ間の電荷はコイルに流れるなどして電圧が低28

静電気学会誌 第44巻 第 1 号（2020）

図 2 典型的な電流，電圧，圧力のプロファイル（当量比 0.8,平均流速 18.5［ m / s ］）
Fig.2 Typical profile of electric current, voltage and pressure in）.
the test section（φ = 0.8，Um = 18.5［ m ⁄ s ］

下している．図 2 には点火プラグから 18［ mm ］下流の地点に設けた圧力孔で測定した圧力も示している．圧力孔は下流オリフィスから上流に 112［ mm ］の位置にあり，点火プラグ近傍の圧力を代表していると考えらえる．
3［ ms ］付近から徐々に圧力が上昇しており，これは点火によるものである．

3.2 放電路径に関する検討
次に得られた放電路の写真から，放電路の径（太さ）を算出することを試みた．図 3 に検討の結果を示す．実験 条 件 は 当 量 比 φ = 0.66， プ ラ グ 付 近 の 圧 力 が 1.57
［ MPa ］の場合である．図 3 の（a）（，b）は撮影した時間が異なっており，各々上段の左側に示したものが高速度カメラにより撮影された画像である．図 3 においては，流動は左手から右手に向かって流れており，それに伴って放電路形状が右手側に向かって大きく変形している様子がわかる．この放電路形状を評価するため，画像を 8bitのビットマップとして保存後，輝度を基に二値化し，そ

図 3 放電路の長さと面積に及ぼす閾値の影響
Fig.3 Effect of threshold on the discharge path lengths and areas.の閉曲線で囲まれた部分の面積と周囲の長さの 1/2 を計測した．二値化の閾値を変えて，閾値がその面積と周囲電路の面積については閾値に対して大きく変化してい長さの 1/2 に及ぼす影響について検討したものが図 3 のる．また図 （3 a）（，b）からわかるように同一の実験であグラフである．図中の上部の右側 2 つの画像は二値化しっても，時刻によって閉曲線を与える閾値の範囲が異なた後の画像であり，二値化の際の閾値は，下記のグラフることがわかる．の点線部分に対応している．図 （3 a）および（b）において以上の見積もりから，本報においては，撮影された画

放電路の形状は大きく異なっており，面積は閾値を下げ像から，放電路の周囲長の 1/2 を放電路長の長さとし，ることにより大きく上昇している．周囲長さの 1/2 につ放電路の面積を放電路長の長さで除したものを放電路径いてみてみると，図 （3 a）
において閾値が 200 以下のとこ（太さ）として定義することとした．図 3 に示す通り，ろでは周囲長さの 1/2 が急落しているが，これは上部と放電路長については，二値化の閾値の影響を受け，また下部の湾曲した部分は結合してしまうことによる．図 3

放電開始時と放電の終了時では閉曲線を生じる閾値が異（b）における閾値が 140 以下のところの急落も同様な理なることから，各実験において一つの閾値を定め，その由である．図 （3 b）においては閾値が 210 以上では単一実験における放電開始から放電終了までの放電路径につの閉曲線が得られず，計測が行えなかった．閉曲線が得いて比較することを試みた．図 4 にその例を示す．放電られている閾値の範囲では，広い範囲で周囲の長さの時間は 1［ ms ］程度であり，35000 fps で測定しているの1/2 は大きな変化をしていないことがわかる．他方で放で，30 枚強の撮影画像が得られており，それらを放電ガソリンエンジン用点火プラグにおける流動により延伸した火花放電の直径に関する考察（傅 建華ら）

29開始からの時刻とともに，上段に並べた．これらを閾値

電路が結合しているために除外，また時刻が 600.0［ μs ］を輝度 70 として二値化した画像が下段である．のものは，放電路が切断されているために除外した．図

二値化後の画像から，適切な閉曲線が得られているも4 においては時刻の前に◎を付したものが，適切な閉曲のを抽出した．例えば時刻が 171.4［ μs ］のものは，放線が得られたと判断したものである．図 4 では適切な閉

図 4 撮影された放電路画像（上段）と閾値 70 で二値化された放電路画像（下段）（空気，平均流速 18.5［ m / s ］，圧力 0.43 MPa（ abs ））Fig.4 Observed original images（ Upper rows ）and binarized
images at 70 a.u. threshold（ Lower rows ）of discharge path
, p = 0.43 MPa（ abs ））.
（ Air, Um = 18.5［ m ⁄ s ］

図 5 各圧力における電流と放電路径との関係
Fig.5 Relationships between electric current and diameter ofdischarge path for various pressure.

静電気学会誌 第44巻 第 1 号（2020）30

曲線が得られていると判断されたものは，6 枚である．
各実験において，閉曲線が得られる画像が多くなるように閾値を調整した．なお空気でなく，当量比が 0.5～1.0

またポアソンの式は，2dne2 (V ) =dx
ε（4）の場合についての検討も行ったが，放電の発光の分布がであり，V は電圧，ε は誘電率である．速度 v は，ここで強く，二値化によって撮影された画像で適切な閉曲線がは簡単に電子の衝突を考えず，電子の質量を m として得られていると判断できるが画像が取得できなかったことから，本報では作動ガスが空気の場合についてのみ述
べる．

12mv = eV2（5）これから，v と ne を消去してこのような手法を用いて，空気の圧力が 0.24［ MPa ］（ abs ）～0.82［ MPa ］（ abs ）の条件における放電路の径を，2d2 (V) = dx

( εj )( 2em ) V
12（6）

-12

その時刻における電流値の関数として示したものが図 5j が一定として，x = 0 で dV / dx = 0，V = 0 の境界条件をである．各グラフにおいては，採用した閾値をあわせて与えて，積分すると，示している．図 4 の画像は，図 （5 b）に対応している．
図 5 では電流を横軸とした．放電路長 lspk，電極間電圧j=-( 49ε )( 2em ) V x
1232（7）-2Vgc，電流 i は時刻により，また雰囲気条件により異なるの関係が得られる（チャイルド・ラングミュアの式）．

が，これらの間には実験式である Kim の式

式（7）から，5）が成り立つとする．Kim の式からプラズマが放電経路に沿って一様であるとして，単位長さあたりの電気抵抗である線抵
抗 R を算出すると，R = Vgc /（lspk i ）∝ i–1.33となり，線抵抗は電流のみの関数となる．線抵抗 R が雰囲気の組成や流速を変化させても，電流のみ関数として整理できる
ことは，著者らの過去の実験

2）

でも実験事実として確かめられている．線抵抗 R に寄与する要素としては，プラズマの断面積と電気伝導率が考えられるが，線抵抗 Rが電流のみ関数として整理できることという実験結果から，ここではプラズマの断面積が電流 i だけの関数であると仮定し，解析を行った．図 5 からは，電流値に対して放電路径をプロットすると電流値が大きいほど，放電

3j ∝ V2（8）である．電流 i は電流密度に面積 S を乗じればよいから，（9）i = jS

面積 S は放電路径 d を用いて，S = ( 4 )d と書けるとすると，
π2

式（2）および（8）を用いると，式（9）から，
3i ∝ V2 ⊗ i

2a（10）（2 1–2a ）
3（11）

すなわち，
V∝i

路径が大きい傾向が見て取れ，電流値と放電路径の間に

a = 0.5 で 0，a = 0.75 で式（11）における i の指数は，（ –1）は相関があることが示唆される．図 5 の各グラフは両対V と i との関係は式/ 3 となる．仮に a = 0.75 であれば，（1）数グラフで示しており，図から電流 i と放電路径 d とのによく似ていることとなる．

間の関係を以上から，式（7）に示されるように電流密度は電圧のd ∝ ia（2）増加に伴って増加するが，式（1）の Kim 式では電流は電圧の増加に伴って減少している．この一見相反する効果と 書 く と，a は 両 対 数 グ ラ フ の 傾 き で あ り，0.5 か ら

は，式（2）
において実験結果から a > 0 であったこと，す0.75 の範囲にあることがわかる．

なわち放電路の直径は電流の増加に伴って増加する効果

3.3 考察を考慮することで，定性的にではあるが説明可能である．

3.2 にて示されたように，放電路径の径と電流値の間上記のモデルが定量的に適応できるかどうかについてには相関があることが示唆されたので，これをもとに（1）考察する．電流密度 j を式（7）を用いて算出してみる．

式の Kim 式について考察を行う．先ずは電流密度 j につ

図 2 を参考に V = 2［ kV ］，x = 1［ mm ］とし，電流密度いて考える ．ここでは x 座標の 1 次元のみを考える．

が j = 2.1 ⊗ 10 ［ A ⁄ m ］であり，図 5 から d = 0.5［ mm ］

電流密度 j は以下で与えられる．

を代入すると，i = 41［ mA ］となる．オーダーとしては8）

j = - nev

（3）n は電子の密度，e は電気素量，v は電子の速度である．

52

図 5 と合致しているものの，電流値としては図 5 よりも小さい値となっている．上記では電子は無衝突としているが，実際には電子の速度は分子などとの衝突によりこガソリンエンジン用点火プラグにおける流動により延伸した火花放電の直径に関する考察（傅 建華ら）

31の見積もりよりも小さくなっているであろうから，このとし，投影面積を周囲長さで除したものを放電路径（太部分については今後修正が必要である．また本報では高さ）として定義した．放電路径 d はその時刻の電流 i と速カメラの撮影画像から閾値を用いて放電路径を算出しの間に d ∝ i の関係が認められ，a は 0.5 から 0.75 の間たが，3.2 で述べた通り，輝度の分布の差が大きく，単の値であった．チャイルド・ラングミュアの式を用いる

純に一様な放電路径となっているかについては，疑問のと a = 0.75 の場合に電圧が電流の（ –1）⁄ 3 乗の比例する余地がある．そのため，放電路径の見積もりについては，ことになり，点火プラグ間の放電における実験式である

今後異なる手法を検討し，高精度化を図っていくことがKim の式を定性的に説明することができた．なお電流の必要と思われる．

絶対値については，誤差が大きく，放電路径の見積もり

a

最後に放電路長について考察する．式（1）の Kim 式での高精度化が必要である．これに加えて Kim の式の放は，放電路が延伸すると，それに伴って抵抗値が上昇す電路長が電圧に比例するという関係についても，改良のることから，放電路長 lspk は電圧 Vgc に比例するモデル式余地があることが示唆された．となっている．上記で式（11）を導出した過程が正しいけ
れば，V と x の関係も同様に考えると，式（7）からは，4V ∝ x3（12）

謝辞
本研究の一部は，総合科学技術・イノベーション会議の SIP（戦略的イノベーション創造プログラム）「革新となり，必ずしも放電路長の長さは電圧に比例しないこ的燃焼技術」（管理法人 JST ）の支援により実施されたとが想定される．

ものであり，ここに記して謝意を表す．

Kim の式（式（1））で算出に用いている放電路長さ は5）点火プラグのギャップから放電路の最遠点を計測し，ギ参考文献

ャップから最遠点までの距離の 2 倍にギャップ長を加え

1） 白石泰介，寺地た，いわゆる放電路長を矩形にて近似したものである．
過去の報告

1-4）

もこの手法を踏襲しているが，図 3 に示しているとおりで，実際の放電路は湾曲しており，この淳，森吉泰生：スパーク放電チャンネ

ル形成に関する点火環境および放電波形特性の影響解析．自動車技術会論文集 46［2］（2015）283-288
2） 今村宰，文鉉太，岩田和也，秋濱一弘：小型風洞を方法では放電路長の見積もりに誤差があることが予測さ用いた高速気流中における火花放電の放電特性に関するれる．また湾曲の効果の他にも，図 3 は一方向から観察

研究．自動車技術会論文集，50［3］（2019）737-742

しているため，観察の奥行方向への変化については観測

3） 佐山勝悟，木下雅夫，政所良行，増田糧，冬頭孝之：できておらず，これも放電路長の見積もりに誤算が生じ

高速流条件での放電経路の短縮化現象のモデル化（第 1る一因であると想像される．本報で展開した議論は，

報）－点火現象の可視化による定量解析－．第 28 回内Kim の式（式（1））における電流の指数を合理的に説明

燃機関シンポジウム講演番号 94，94-20178077（2017）

できる可能性があるものであるが，放電路長の指数につ4） 堀司：火花点火機関における放電経路伸長と再放電のいては Kim の式（式（1））と合致していない．これから放モデリング．自動車技術会論文集，48［3］（2017）641-647

電路長については，湾曲部，3 次元変形を含めて精度よ

5） Jaehong Kim and Richard W. Anderson, Spark Anemometry of

く計測し，吟味し直す必要があることを示している．

Bulk Gas Velocity at the Plug Gas of a Firing Engine, SAE
952459（1995）

4．まとめ

6） 今村宰，文鉉太，岩田和也，秋濱一弘：小型風洞を一様流速を生じさせることができる実験装置を用い用いた高流気流中における火花放電経路の流速追従性にて，流動中における火花放電の電流，電圧プロファイル関する研究．2018 年度日本機械学会年次大会予稿集，の取得，および高速度カメラを用いて放電路の変形挙動

J0710104（2018）について実験的に把握を行い，特に放電路径の変化に着

7） 鄭 棟元，佐々木耕作，菅田健志，横森 剛，飯田訓正：

目して検討を行った．その結果，明らかになったことは火花放電のパターンがリーンバーン SI エンジンにおけ以下のとおりである．

る燃焼のサイクル変動に及ぼす影響．第 27 回内燃機関

流動によって延伸した放電路の写真を一定の閾値を決め二値化し，閉曲線が得られる場合について周囲長さ，および投影面積を算出した．周囲長さの 1/2 を放電路長

シンポジウム講演予稿集，講演番号 6（2016）8） 高村秀一：プラズマ理工学入門，pp145-149，森北出版社（1997）