電気めっきシミュレーション
理論と物理
電気めっきシミュレーションとは
電気めっきのシミュレーションって何を予測するんですか？
ざっくり言うと、めっき膜厚の均一性を予測する技術だ。ワーク（被めっき体）の表面全体に同じ厚さのめっきを付けたいのに、実際は場所によって厚い・薄いが出てしまう。これを事前に計算で見極めるのがめっきシミュレーションだよ。
なんで場所によって厚さが変わるんですか？
電流は「電気的に近い道」を通りたがるからだ。例えば自動車のバンパーにクロムめっきをする場合、アノード（陽極）に近い凸部には電流が集中してめっきが厚くなる一方、凹部や穴の奥は電流密度が低くなり膜が薄くなる。この「スロウイングパワー（均一電着性）」の問題が、めっきシミュレーションが生まれた最大の動機だよ。
具体的にどんな場面で使われてるんですか？
代表的な用途を挙げると——
- PCBのビアホール内めっき：直径0.1〜0.3 mmの穴の奥まで均一にCuを付ける。スルーホールのアスペクト比が高くなるほど困難
- 半導体のダマシン配線：幅10 nm以下のトレンチをCuで完全充填（スーパーフィリング）
- 自動車の装飾めっき：複雑な3D形状のバンパーやグリルに均一なクロムめっき
- 航空機部品の硬質クロムめっき：耐摩耗性が要求される部品での膜厚管理
- めっき治具・ラック設計：補助陰極や遮蔽板の配置を最適化して均一性を向上
結構身近な製品に使われてるんですね！ じゃあ、理論的にはどうやって計算するんですか？
めっきシミュレーションは精度に応じて3段階のモデルがある。一次電流分布、二次電流分布、三次電流分布だ。段階が上がるほど物理を精密に取り込む代わりに計算コストが増える。順番に見ていこう。
一次電流分布（ラプラス方程式）
まず一次電流分布ってどういうものですか？
一次電流分布は電極反応の抵抗を完全に無視して、電解液中の電位分布だけで電流を求めるモデルだ。電解液の導電率 $\kappa$ が一様で化学反応の影響がないと仮定すると、電位 $\phi$ はラプラス方程式に従う：
電流密度は電位の勾配から得られる。これがオームの法則の微分形だ：
え、これだけ？ かなりシンプルですね。
そう、方程式自体はシンプルだ。でも一次分布は最も不均一な分布を与える「最悪ケース」なんだ。電極面を等電位面として扱うため、電流はひたすら幾何学的に最短距離を通ろうとする。凸部のエッジに電流が極端に集中する。実際のめっきでは電極反応に抵抗があるため、これよりは均一化する方向に働く。だから一次分布は「均一性の下限値」として使える。
一次電流分布の境界条件
- カソード（陰極）面：$\phi = \phi_{\text{cathode}}$（等電位）
- アノード（陽極）面：$\phi = \phi_{\text{anode}}$（等電位）
- 絶縁壁面：$\frac{\partial \phi}{\partial n} = 0$（電流の法線成分がゼロ）
電極/電解液界面での過電圧 $\eta$ を無視していることが一次分布の本質的な仮定。
二次電流分布（バトラー・ボルマー方程式）
二次になると何が変わるんですか？
二次電流分布では、電極表面での電荷移動反応（電極反応速度論）を考慮する。つまり、金属イオンが電子を受け取って金属原子になる反応には一定の「抵抗」がある。この関係を記述するのがバトラー・ボルマー方程式（Butler-Volmer equation）だ：
記号がたくさん出てきました…。ひとつずつ教えてください。
| 記号 | 意味 | 典型値・補足 |
|---|---|---|
| $i$ | 電極面の電流密度 | [A/m²] |
| $i_0$ | 交換電流密度 | 反応系に依存（Cu: 約1〜10 A/m²、Ni: 約0.1〜1 A/m²） |
| $\alpha_a, \alpha_c$ | アノード側・カソード側の移行係数 | 通常 $\alpha_a + \alpha_c = 1$（単電子反応の場合） |
| $F$ | ファラデー定数 | 96,485 C/mol |
| $\eta$ | 過電圧（活性化過電圧） | $\eta = \phi_m - \phi_s - E^{\text{eq}}$（金属電位 - 溶液電位 - 平衡電位） |
| $R$ | 気体定数 | 8.314 J/(mol·K) |
| $T$ | 絶対温度 | [K] |
過電圧 $\eta$ がキーパラメータなんですね。これが大きいほど反応が速い？
その通り。過電圧 $\eta$ は平衡状態からの駆動力で、$|\eta|$ が大きいほど電流が増える。でも重要なのは、この式が非線形境界条件として機能する点だ。電解液内のラプラス方程式は一次分布と同じだが、電極面で $\phi$ を等電位と置く代わりにバトラー・ボルマー式で $i$ と $\eta$（つまり $\phi$）を結びつける。この非線形性のおかげで、凸部に集中しようとする電流が「反応抵抗の壁」にぶつかって、結果として一次分布より均一化する。
なるほど！ 反応が遅い（$i_0$ が小さい）めっき液のほうが、逆に均一に付くってことですか？
するどいね。その通りで、$i_0$ が小さいほど電極反応が律速になり、電流分布は均一化する方向に向かう。シアン系の銅めっきが均一性に優れるのは、この反応速度論的な効果が大きいからだ。
ターフェル近似（高過電圧領域）
$|\eta| \gg RT/F$（約25 mV at 25°C）のとき、バトラー・ボルマー式のカソード側は指数関数1つに簡略化できる：
$$ i \approx -i_0 \exp\left(-\frac{\alpha_c F \eta}{RT}\right) $$対数を取ると線形関係（ターフェルプロット）が得られ、実験データから $i_0$ と $\alpha_c$ を求める定番の方法だ：
$$ \eta = a + b \log|i| $$ここで $a = -\frac{RT}{\alpha_c F}\ln i_0$、$b = -\frac{2.303RT}{\alpha_c F}$（ターフェル勾配）。
三次電流分布（ネルンスト・プランク方程式）
三次電流分布ではさらに何を考えるんですか？
三次電流分布では物質輸送（マスストランスポート）を考慮する。めっきが進むと電極近傍の金属イオン濃度が低下し（濃度分極）、供給が追いつかなくなる。イオンの移動を支配するのがネルンスト・プランク方程式だ：
右辺の3つの項はそれぞれ：
- $-D_i \nabla c_i$：拡散（濃度勾配によるフィックの法則）
- $-\frac{z_i F D_i c_i}{RT} \nabla \phi$：泳動（電場による荷電イオンの移動）
- $c_i \mathbf{v}$：対流（液の流れによる輸送）
ここで $D_i$ は拡散係数、$c_i$ は種 $i$ の濃度、$z_i$ は電荷数、$\mathbf{v}$ は流速ベクトルだ。
対流も入ってくるということは、流体の計算も必要になるんですか？
そうだ。三次分布を厳密にやるなら、ナビエ・ストークス方程式で流れ場を解き、それをネルンスト・プランク方程式の対流項に入れる。かなり重たい計算になるよ。ただ実務では、電極近傍の拡散境界層をネルンスト拡散層モデルで近似して、層の厚さ $\delta_N$（典型値10〜500 $\mu$m）だけをパラメータとして与える簡略法がよく使われる。
例えばPCBのビアホールだと、穴の奥は液が動かないから濃度が下がりやすい、ってことですよね？
まさにその通り。ビアホールの奥はほぼ停滞状態で対流がない。拡散だけでイオンが補給されるため、アスペクト比（穴の深さ/直径）が大きくなるほど濃度分極が厳しくなり、膜厚が薄くなる。この効果を正確に予測するのが三次電流分布の出番だ。
限界電流密度
電極表面のイオン濃度がゼロになると、それ以上電流は流せない。この上限が限界電流密度$i_L$：
$$ i_L = \frac{nFD c_{\infty}}{\delta_N} $$$c_{\infty}$ はバルク濃度、$\delta_N$ はネルンスト拡散層厚さ。$i > i_L$ になると析出品質が急激に劣化し、焦げ・樹枝状析出が生じる。
電気的中性条件
電解液のバルク領域では局所的な電荷は中性と仮定：
$$ \sum_i z_i c_i = 0 $$この条件と各種のネルンスト・プランク方程式、電極面のバトラー・ボルマー式を連立させるのが三次電流分布の完全な定式化。
Wagner数と均一電着性
一次・二次・三次って選ぶ基準はあるんですか？ 全部三次でやれば安全？
計算コストが全然違うから、必要十分な精度のモデルを選ぶ判断が重要だ。その判断に使えるのがWagner数（ワグナー数）Wa だ：
$\kappa$ は電解液の導電率、$L$ はワークの代表寸法、$\partial\eta/\partial i$ は分極曲線の傾き（電極反応の微分抵抗）。物理的には「電極反応抵抗 / 電解液抵抗」の比だ。
- Wa ≪ 1：電解液抵抗が支配的 → 一次分布に近い（不均一）
- Wa ≫ 1：電極反応抵抗が支配的 → 均一な電流分布
- Wa ≈ 1：両方が同程度の影響 → 二次分布の計算が必須
つまりWa数を見れば、一次で済むか二次が必要か判断できるんですね。現場ではどれくらいの値ですか？
| めっき浴 | 典型的なWa | 均一性の傾向 |
|---|---|---|
| 硫酸銅めっき（酸性） | 0.1〜1 | 低〜中 |
| シアン化銅めっき | 5〜50 | 高（優れた均一性） |
| ワットニッケルめっき | 0.5〜5 | 中 |
| 硬質クロムめっき | 0.01〜0.1 | 非常に低 |
| 錫めっき | 1〜10 | 中〜高 |
硬質クロムめっきはWaが0.01って…すごく不均一になりそうですね。
その通り。だからクロムめっきでは補助陰極（ロッバー）や遮蔽板を使って電流分布を制御する。これらの治具配置の最適化こそ、シミュレーションの最も得意とするところだ。
膜厚成長モデル
電流分布が分かったら、膜厚はどう計算するんですか？
ファラデーの電解法則で、局所電流密度から膜の成長速度を求める：
- $\delta$：膜厚 [m]
- $M$：析出金属のモル質量 [kg/mol]（Cu: 0.0636, Ni: 0.0587, Cr: 0.0520）
- $i_n$：電極面の法線電流密度 [A/m²]
- $n$：電子移動数（Cu²⁺→Cu: $n=2$, Cr⁶⁺→Cr: $n=6$）
- $\rho$：析出金属の密度 [kg/m³]
- $\eta_{\text{CE}}$：電流効率（副反応分を差し引く。クロムめっきは約10〜25%と極めて低い）
例えば硫酸銅めっきで $i = 300$ A/m²、電流効率100%の場合、成長速度は約3.7 $\mu$m/min になる。
スマホの中のめっきシミュレーション
あなたのスマートフォンのプロセッサには、幅10 nm以下の銅配線が数十億本張り巡らされています。このすべてが電気めっき（ダマシンプロセス）で形成されています。特に難しいのが、アスペクト比の高い（細くて深い）トレンチの底から銅を成長させる「スーパーフィリング」です。PEG（ポリエチレングリコール）などの抑制剤、SPS（ビス(3-スルホプロピル)ジスルフィド）などの加速剤、JGB（ヤヌスグリーンB）などの均一化剤——3種類の有機添加剤の吸着・拡散・消費の競争を物質移動モデルで計算し、ボトムアップ充填を実現するプロセス条件を探索します。半導体1チップの製造に必要なめっきシミュレーションの計算量は、自動車バンパー1個分の数千倍とも言われます。
数値解法と実装
FEMとBEMの使い分け
めっきシミュレーションではどんな数値解法が使われるんですか？ 有限要素法ですか？
主に2つの方法がある。
- 有限要素法（FEM）：電解液領域全体を要素で分割する。二次・三次電流分布や非線形材料（導電率の濃度依存性など）に強い。COMSOLやAnsysが使用
- 境界要素法（BEM）：電極と絶縁壁の表面だけをメッシュ化する。一次分布（ラプラス方程式）なら内部の体積メッシュが不要でメッシュ生成が楽。3D複雑形状のめっき治具設計に有利。ElSyCA PlatingMasterが代表的
BEMのほうが楽そうに聞こえますけど、デメリットはないんですか？
BEMは密行列（full matrix）になるので、要素数の2乗でメモリが増える。数万要素を超えると計算が重くなる。あと、電解液の導電率が場所によって変わる場合（濃度依存など）は適用が困難だ。実務では一次・二次分布まではBEM、三次分布が必要ならFEM、という使い分けが多い。
移動境界の取り扱い
めっきが進むと膜が厚くなって形状が変わりますよね？ それはどう計算するんですか？
いい質問だ。めっきの堆積によって電極形状が時々刻々変わる——いわゆる移動境界問題だ。主なアプローチは3つ：
- ALE法（Arbitrary Lagrangian-Eulerian）：メッシュを膜成長に合わせて変形させる。COMSOLの Deformed Geometry モジュールが使用。メッシュの歪みが大きくなるとリメッシングが必要
- レベルセット法：界面を暗黙的に追跡する。トポロジー変化（穴の充填完了など）に強い
- 準定常近似：膜厚の変化が遅い場合、電流分布を「その瞬間のフリーズ形状」で解いて膜厚を更新し、形状を進める。多くの実務ケースで十分に精度が出る
準定常近似が使えるかどうかはどう判断するんですか？
膜厚の変化量がワークの代表寸法に比べて十分小さければOK。例えば膜厚30 $\mu$m、ワークの凹凸が数mm以上なら、形状変化は0.1%以下で準定常が成り立つ。一方、半導体のトレンチ充填（穴幅100 nm、膜厚50 nm）では形状が劇的に変わるから、ALE法やレベルセット法が必須だ。
非線形ソルバーと収束戦略
バトラー・ボルマー式は非線形ですよね。収束が難しいイメージがあるんですけど…
その通り。バトラー・ボルマー式の指数関数は過電圧に対して急激に変化するから、ニュートン法の初期値が悪いとあっさり発散する。実務的なテクニックをまとめると：
- 初期値には一次分布の解を使う：ラプラス方程式の解をウォームスタートにする
- ダンピング（緩和）付きニュートン法：更新量に0.3〜0.7の緩和係数をかけて発散を防ぐ
- 電流の段階的増加：最終的な印加電流を一気にかけず、0.1倍→0.5倍→1.0倍と段階的に上げる（擬似時間ステッピング）
- ターフェル近似から切り替える：まずターフェル近似（対数関数＝穏やかな非線形性）で初期解を得て、バトラー・ボルマー式に切り替える
なるほど、いきなり本番条件で解こうとしちゃダメなんですね。
そう。電気化学の非線形方程式は「助走」が大事だ。COMSOLだと Auxiliary Sweep でパラメータを段階的に変える機能がある。Ansys Fluent なら Pseudo Transient を有効にすると安定しやすい。
弱形式の定式化
電解液中のラプラス方程式の弱形式は、試験関数 $w$ を用いて：
$$ \int_\Omega \kappa \nabla \phi \cdot \nabla w \, d\Omega = \oint_{\Gamma} i_n \cdot w \, d\Gamma $$右辺の境界積分がバトラー・ボルマー式で与えられるため、全体として非線形な連立方程式系：
$$ \mathbf{K} \boldsymbol{\phi} = \mathbf{f}(\boldsymbol{\phi}) $$ニュートン法で $\mathbf{J} \Delta\boldsymbol{\phi} = -\mathbf{r}$（$\mathbf{J}$はヤコビアン、$\mathbf{r}$は残差）を反復的に解く。ヤコビアンにはバトラー・ボルマー式の微分 $\partial i / \partial \eta$ が含まれる。
実践ガイド
解析ワークフロー
実際にめっきシミュレーションをやるとき、どんな手順で進めるんですか？
標準的なワークフローはこうだ：
- CADモデルの準備：めっき槽、アノード、カソード（ワーク）、治具（補助電極・遮蔽板）の3Dジオメトリ。実際に電解液が満たされる領域を抽出する
- 電解液ドメインの定義：導電率 $\kappa$、反応種の拡散係数 $D_i$、バルク濃度 $c_\infty$ を設定
- 電極反応パラメータの設定：交換電流密度 $i_0$、移行係数 $\alpha$、平衡電位 $E^{\text{eq}}$、電子移動数 $n$
- メッシング：電極近傍の境界層を細かく、バルク領域は粗く
- ソルバー設定と実行：上述の非線形収束戦略を適用
- 後処理：電流密度分布、電位コンター、膜厚分布を可視化し、均一性を評価
メッシングの勘所
めっきシミュレーション特有のメッシュの注意点ってありますか？
めっき解析のメッシュで気をつけるべきポイントは構造解析とかなり違う。重要な順に：
- 電極エッジの細分化：カソードのエッジ（角、穴のフチ）には電流が集中する。エッジから数mmの範囲で要素サイズ0.1〜0.5 mmが目安
- 拡散境界層の解像：三次分布では電極面から数十$\mu$mの拡散層内に3〜5層の要素が必要。バウンダリレイヤーメッシュ（プリズム要素）が有効
- アノード-カソード間のギャップ：両電極の距離方向に最低10要素。ここの解像度が電解液の電位降下に直結する
- バルク領域は粗くてOK：電極から離れた領域は電位勾配が小さいので、要素サイズを10倍にしても結果が変わらないことが多い
構造解析と違って、応力集中部じゃなくて電極のエッジに集中するんですね。
正確に言うと物理は同じで、「ラプラス方程式の解が特異的になる場所」にメッシュを集中させるんだ。90度の角なら理論上 $r^{-1/3}$ の特異性がある（$r$ はエッジからの距離）。構造解析の応力特異点と本質的に同じ数学だよ。
境界条件の設定
境界条件はどう設定するのが一般的ですか？
| 境界 | 一次分布 | 二次分布 | 三次分布 |
|---|---|---|---|
| カソード面 | $\phi = \phi_c$（等電位） | Butler-Volmer式 | Butler-Volmer式 + 濃度依存 |
| アノード面 | $\phi = \phi_a$（等電位） | Butler-Volmer式 or 全電流指定 | 同左 + 溶解反応速度 |
| 絶縁壁 | $\partial\phi/\partial n = 0$ | 同左 | 同左 + $\partial c/\partial n = 0$ |
| 液面（開放面） | $\partial\phi/\partial n = 0$ | 同左 | $c = c_\infty$（バルク濃度） |
実務上のポイント：アノードは溶解型（Cu, Ni）と不溶性（Pb-Sn, Ti-Pt, DSA）に大別され、不溶性アノードでは酸素発生反応の速度論も含める必要がある。
検証と妥当性確認
シミュレーション結果の正しさをどうやって確認するんですか？
めっきシミュレーションのV&Vで定番の手法を3段階で紹介するよ。
- 解析解との比較（Verification）：同心円筒電極（ワグナーの解析解あり）やハル・セル（Hull Cell）の理論解で一次分布を検証。全電流の積分値が入力電流と一致するかも必ず確認
- ハル・セル試験との比較（Validation）：267 mLハル・セルは標準化された実験装置で、アノードからカソードまでの距離が直線的に変わる。実測の膜厚分布とシミュレーション結果を比較する
- 実機のトライめっきとの比較：最終的には実際のめっき浴で試作品のめっきを行い、蛍光X線（XRF）やクロスセクションで膜厚を実測して比較する
ハル・セルって実験だけじゃなくシミュレーションの検証にも使えるんですね！
そう。ハル・セルは「一次電流分布の解析解が知られている簡単な形状」だから、シミュレーションコードの基本的な正しさを確認するのに最適なんだ。ASTM B571にも規定されている標準的な試験法だよ。
境界条件の考え方——「水槽のたとえ」
めっき槽の境界条件を設定するとき、「水槽に電極を沈めた状況」をイメージしよう。電極面は「電流を生む場所」（バトラー・ボルマー式）、槽壁面は「電流を通さない壁」（絶縁条件）、液面は「何事もなく通り過ぎる場所」（自然境界）。補助電極を追加するのは「水槽の中に堰を作って水流を制御する」のと同じだ。形状が電流の流れを支配するから、治具設計は流路設計に似ている。
ソフトウェア比較
対応ツールと機能比較
めっきシミュレーションに使えるソフトって、何がありますか？
| ソフトウェア | 解法 | 電流分布 | 移動境界 | 添加剤モデル | 強み |
|---|---|---|---|---|---|
| COMSOL Multiphysics （Electrodeposition Module） | FEM | 一次〜三次 | ALE / レベルセット | カスタムPDE | マルチフィジクスの柔軟性。電気化学分野で最も実績豊富 |
| Ansys Fluent （Electrochemistry Model） | FVM | 二次・三次 | 動的メッシュ | UDF | 流れとの連成が得意。大規模並列計算 |
| STAR-CCM+ （Electrochemistry） | FVM | 二次・三次 | モーフィングメッシュ | Java API | ポリヘドラルメッシュの自動生成 |
| ElSyCA PlatingMaster | BEM | 一次・二次 | 表面メッシュ更新 | 限定的 | 3D複雑形状の治具設計に特化。高速 |
| ELSYCA CuPlate | BEM+FEM | 一次〜三次 | 対応 | 対応 | PCBめっき専用。ビアホール充填に強い |
結局どれを使えばいいんでしょうか？
用途次第だ。
- 研究・論文用：COMSOLが圧倒的。カスタム方程式を自由に書けるし、文献でも最も引用される
- 自動車・装飾めっきの治具設計：PlatingMasterがBEM特有の高速さで定評あり
- 流れの影響が大きい場合（噴流めっき・バレルめっき）：Ansys FluentやSTAR-CCM+のCFD連成が強い
- 半導体/PCBの穴埋め：CuPlateやCOMSOLのレベルセット法
もし所属企業が既にAnsysのライセンスを持っているなら、まずFluentの電気化学モデルで始めるのが初期投資ゼロで現実的だ。
先端技術
添加剤モデリング
実際のめっき液には添加剤が入ってますよね？ あれもシミュレーションできるんですか？
できる。特に半導体のCuめっきでは添加剤のモデリングが不可欠だ。3種類の添加剤が競争的に電極表面に吸着して析出速度を制御する：
- 抑制剤（Suppressor）：PEGなど。電極面に吸着して析出を遅くする。凸部に多く吸着
- 加速剤（Accelerator）：SPSなど。表面に吸着して析出を速める。凹部の底に蓄積
- 均一化剤（Leveler）：JGBなど。対流依存的に吸着し、凸部の加速剤の効果を打ち消す
これらの表面被覆率 $\theta$ を追跡する表面拡散・吸着モデルと、バルク中の輸送方程式を連立させる。CUSOMOLではCOMSOLの Surface Reactions インターフェースで実装できる。
パルスめっきシミュレーション
パルスめっきって聞いたことがあるんですけど、シミュレーションは難しくなりますか？
パルスめっきはON/OFF（またはフォワード/リバース）を周期的に繰り返す手法で、DC（直流）めっきより均一性と結晶性が向上する。シミュレーション上は：
- ON期間：通常の電流分布計算（ミリ秒〜秒のオーダー）
- OFF期間：電流ゼロだが拡散層内の濃度が回復（拡散方程式を解く）
- リバースパルス：逆電流でカソード表面の凸部を優先的に溶解
時間スケールがミリ秒と秒で大きく異なるので、適応的時間ステッピングが重要だ。ON/OFF切り替え直後は $\Delta t$ を極めて小さく（$\mu$s）、定常に近づいたら大きくする。
デジタルツインとプロセス最適化
最近はデジタルツインっていう言葉もよく聞きますけど、めっきでも使われてるんですか？
先端的な工場では、めっきラインのデジタルツインが実用化されつつある。具体的には：
- リアルタイム濃度モニタリング：液分析センサーのデータをシミュレーションモデルに入力し、膜厚分布をリアルタイム予測
- 治具形状の自動最適化：遮蔽板の形状・位置をパラメータとして、膜厚の標準偏差を最小化する最適化ループ
- サロゲートモデル：高精度FEMを大量に走らせた結果から機械学習モデルを構築し、推論時間をミリ秒単位に短縮
- プロセスウィンドウの可視化：電流密度・温度・液流速の組み合わせマップで、品質OK/NGの領域を明示
特に自動車業界では、新車種の複雑な外装部品のめっき治具を試作なしで設計するワークフローが広がっている。
トラブルシューティング
よくあるエラーと対策
めっきシミュレーションでよくハマるポイントを教えてください！
| 症状 | 原因 | 対策 |
|---|---|---|
| ニュートン法が発散する | バトラー・ボルマー式の初期値が悪い | 一次分布の解を初期値に使用。ダンピング係数0.3〜0.5で開始 |
| エッジ付近で電流密度が異常に大きい | 一次電流分布の特異性（幾何学的効果） | エッジ細分化（要素サイズ0.1mm以下）。丸みを付ける（実形状に忠実に） |
| 膜厚の合計がファラデー則と合わない | 電流効率の未設定、またはメッシュの電流積分誤差 | 電流効率 $\eta_{\text{CE}}$ を正しく設定。境界積分の精度をメッシュ収束で確認 |
| 三次分布で負の濃度が出る | 拡散境界層の要素が粗すぎる | 境界層メッシュを追加（最低3層）。上流化スキームの適用 |
| ALE法でメッシュが潰れる | 膜厚成長に対してリメッシュ頻度が不足 | リメッシュ間隔を短くする。変形量が要素サイズの30%を超えたらリメッシュ |
| 計算時間が異常に長い | パルスめっきで時間ステップが小さすぎる | ON/OFF切り替え後のみ微小ステップ。定常部は大ステップに自動切り替え |
エッジの特異性って、モデルの角を丸めるだけで解消するんですか？
実は実物のエッジも完全なシャープエッジではないから、微小な丸み（フィレット R=0.1〜0.5 mm）を付けるのが物理的にも正しい。シミュレーション上は特異性が緩和されてメッシュ収束性が大幅に改善する。構造解析で応力集中部にフィレットを付けるのと同じ発想だよ。
まずはハル・セルの一次分布から自分で解いてみるといい。COMSOLのチュートリアルに例題がある。手を動かすのが一番の近道だからね。
初心者が陥りやすい落とし穴
「全部三次分布で計算すれば正確だろう」——この考えは危険だ。三次分布には拡散係数、バルク濃度、拡散層厚さなど追加のパラメータが必要で、これらが不正確なら二次分布より信頼性が下がることもある。まず一次・二次でWagner数を確認し、物質輸送が支配的な場合にのみ三次に進むのが正しい順序。「パラメータが増える = 不確実性が増える」ことを常に意識しよう。
関連トピック
なった
詳しく
報告