fullseye

conv_gci — CONVERGE verification op

使い方

Grid convergence index of the finest solution: Celik et al. 2008 procedure, eqs. (1)-(7).

Takes the same inputs as :func:conv_observed_order (three solutions, finest first, and the refinement) and returns, with p from eq. (3):

phi_ext21 = (r21^p phi1 - phi2) / (r21^p - 1)                  (4)   (and phi_ext32)
e_a21     = |(phi1 - phi2) / phi1|                            (5)
e_ext21   = |(phi_ext12 - phi1) / phi_ext12|                   (6)
GCI_fine21 = 1.25 e_a21 / (r21^p - 1)                         (7)

(safety_factor = 1.25 as in eq. (7)). Fractions and percentages are both returned, and report reads “phi1 +- GCI”. Note what the GCI is not: it bounds the discretization error of the fine solution only — “this does not account for modeling errors” (Celik).

Refusal as in :func:conv_observed_order, plus phi1 = 0 or phi_ext = 0 (relative errors undefined). Defaults: Celik Table 1, first column (GCI 2.2 %).

詳しい使い方ガイド

参考(サンプルデータ・文献)

実行できる例(この op を実際に呼ぶ検証済みサンプル)

型が繋がる次の op(table を入力に取れる)

—

同カテゴリ(verification)

conv_observed_order


Provenance: converge.py — CONVERGE operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.