[{"data":1,"prerenderedAt":131},["ShallowReactive",2],{"$f1fug44q4cq41n":3},{"slug":4,"publishedAt":5,"category":6,"title":7,"summary":8,"image":9,"imageAlt":14,"readingMinutes":15,"serviceId":16,"keywords":17,"imageCaption":22,"sections":23,"sources":121},"mesh-independence-study","2026-10-12","methods","Исследование сеточной независимости: пример расчёта GCI","Пошаговое исследование сеточной независимости: три сетки, коэффициент измельчения, наблюдаемый порядок и индекс сеточной сходимости, с итоговой фразой для отчёта.",{"src":10,"thumbnail":11,"width":12,"height":13},"\u002Fassets\u002Fcfd-cases\u002Fshell\u002Fimages\u002F06.webp","\u002Fbrand\u002Fstudies\u002Fshell.svg",1448,1086,"Иллюстрация DEEB-FLOW: график потерь давления от массового расхода",4,"review",[18,19,20,21],"сеточная независимость","индекс сеточной сходимости","GCI CFD","сеточная сходимость","Концептуальная иллюстрация из библиотеки кейсов DEEB-FLOW. Она не является результатом из этой статьи.",[24,30,36,48,98,105,115],{"id":25,"title":26,"paragraphs":27},"two-meshes","Почему две сетки ничего не доказывают",[28,29],"«Мы посчитали на двух сетках, и поля выглядели одинаково». Эта фраза встречается в очень многих отчётах и означает не то, на что надеялся автор. За двумя похожими картинками полей может скрываться 10% разницы в том единственном числе, которое нужно заказчику.","Исследование сеточной независимости отвечает на вопрос, насколько изменился бы результат на бесконечно мелкой сетке, и отвечает числом. В статье такое исследование разобрано от трёх сеток до фразы для отчёта. Используется методика, опубликованная Челиком с соавторами в 2008 году и основанная на индексе сеточной сходимости Роуча.",{"id":31,"title":32,"paragraphs":33},"what-you-find-out","Что вы хотите выяснить?",[34,35],"Любое решение методом конечных объёмов содержит погрешность дискретизации, которая уменьшается с уменьшением ячеек. Если сетка достаточно мелкая, погрешность убывает упорядоченно: при схеме второго порядка уменьшение размера ячейки вдвое снижает её примерно вчетверо.","Имея решения на трёх сетках, можно измерить, как быстро меняется результат, экстраполировать его к нулевому размеру ячейки и сказать, насколько результат на мелкой сетке далёк от этого предела. Индекс сеточной сходимости GCI и есть это расстояние с коэффициентом запаса. Он ограничивает только численную погрешность и ничего не говорит о том, верна ли модель турбулентности.",{"id":37,"title":38,"paragraphs":39},"four-rules","Четыре правила до начала",[40,47],{"list":41,"ordered":46},[42,43,44,45],"**Выберите величину, которая важна.** Потери давления, тепловая мощность, максимальная температура. Проводите исследование на том, что войдёт в отчёт.","**Измельчайте сетку везде и в одинаковое число раз.** Меняйте глобальный размер ячейки. Если измельчать только там, где вы предполагали градиенты, вы проверяете собственное предположение.","**Сохраняйте одинаковый подход к пристеночной области.** При пристеночных функциях первая ячейка должна оставаться в допустимом диапазоне на всех трёх сетках. Если пограничный слой разрешается, измельчайте его вместе с остальным.","**Доводите каждый расчёт до сходимости.** Итерационная погрешность должна быть намного меньше различий, которые вы пытаетесь измерить.",true,"Используйте коэффициент измельчения не менее 1,3 по размеру ячейки. В трёх измерениях это означает рост числа ячеек примерно в 2,2 раза на каждом шаге, а коэффициент 1,5 даёт рост примерно в 3,4 раза.",{"id":49,"title":50,"paragraphs":51},"worked-example","Разбор на примере",[52,53,72,73,76,77,81,82,85,86,89,90,94,95],"Трубопроводный фитинг, стационарное течение, искомая величина: потери давления. Три сетки с постоянным коэффициентом измельчения 1,5. Числа подобраны для иллюстрации методики.",{"table":54},{"head":55,"rows":59},[56,57,58],"Сетка","Число ячеек","Потери давления (Па)",[60,64,68],[61,62,63],"3, грубая","0.40 × 10⁶","1297",[65,66,67],"2, средняя","1.35 × 10⁶","1262",[69,70,71],"1, мелкая","4.56 × 10⁶","1245","**Шаг 1, коэффициент измельчения.** По числу ячеек в трёх измерениях он равен 1,5 для каждой пары сеток.",{"formula":74},[75],"r = (N₁ \u002F N₂)^(1\u002F3) = (4.56 \u002F 1.35)^(1\u002F3) = 1.5","**Шаг 2, разности.** Обе положительны, вторая меньше первой, значит, результат сходится монотонно.",{"formula":78},[79,80],"ε₃₂ = 1297 − 1262 = 35 Pa","ε₂₁ = 1262 − 1245 = 17 Pa","**Шаг 3, наблюдаемый порядок.** Формально схема имеет второй порядок, а наблюдается 1,78. Это нормальный результат: реальные сетки с пограничными слоями и ограничителями редко дают ровно 2.",{"formula":83},[84],"p = ln(ε₃₂ \u002F ε₂₁) \u002F ln(r) = ln(35 \u002F 17) \u002F ln(1.5) = 1.78","**Шаг 4, экстраполяция к нулевому размеру ячейки.** Это оценка того, что дала бы бесконечно мелкая сетка.",{"formula":87},[88],"φ_ext = φ₁ − ε₂₁ \u002F (r^p − 1) = 1245 − 17 \u002F 1.06 = 1229 Pa","**Шаг 5, индекс сеточной сходимости.** Относительное изменение между мелкой и средней сеткой составляет 1,37%. С коэффициентом запаса 1,25, рекомендованным для исследования на трёх сетках, GCI мелкой сетки равен 1,6%, средней 3,3%.",{"formula":91},[92,93],"GCI_fine = 1.25 × 1.37% \u002F (r^p − 1) = 1.6%","GCI_medium = 1.25 × 2.77% \u002F (r^p − 1) = 3.3%","**Шаг 6, проверка асимптотической области.** Значение, близкое к 1, говорит о том, что все три сетки находятся там, где погрешность убывает упорядоченно, и экстраполяции можно доверять.",{"formula":96},[97],"GCI_medium \u002F (r^p × GCI_fine) = 3.27 \u002F (2.06 × 1.61) = 0.99",{"id":99,"title":100,"paragraphs":101},"report","Что написать в отчёте",[102,104],{"quote":103},"Потери давления на мелкой сетке составляют 1245 Па при численной неопределённости 1,6% (индекс сеточной сходимости, три сетки, коэффициент измельчения 1,5, наблюдаемый порядок 1,78). Экстраполированное значение равно 1229 Па.","И одно решение. Если для задачи допустимы 3,3%, для остальных вариантов достаточно средней сетки, а это меньше трети затрат. В этой экономии и состоит практический смысл проводить исследование в начале проекта.",{"id":106,"title":107,"paragraphs":108},"when-numbers-misbehave","Когда числа ведут себя плохо",[109],{"list":110},[111,112,113,114],"**Разности меняют знак.** Это осциллирующая сходимость, и экстраполяция недействительна. Сначала проверьте итерационную сходимость, затем добавьте более мелкую сетку.","**Наблюдаемый порядок далёк от формального.** Ниже примерно 0,5 или намного выше порядка схемы. Относитесь к GCI с осторожностью и скажите об этом. Частые причины: грубая сетка пропускает особенность течения, или пристеночные функции на одной из сеток вышли из допустимого диапазона.","**Разности очень малы.** Если все три результата совпадают до четырёх знаков, проверьте, действительно ли сетки различаются там, где это важно.","**Результат всё ещё меняется, а на четвёртую сетку нет ресурсов.** Приведите результат на мелкой сетке с его GCI, укажите, что асимптотическая проверка не пройдена, и считайте число предварительным.",{"id":116,"title":117,"paragraphs":118},"proportionate-effort","Какой объём работы соразмерен задаче?",[119,120],"Полное исследование на трёх сетках с указанием GCI оправдано, когда важно абсолютное значение: гарантия, сравнение с испытаниями, публикация. Для быстрого сравнения двух вариантов часто достаточно одного измельчения на интересующей величине, чтобы показать, что порядок вариантов не меняется. Называйте это проверкой чувствительности: оценку неопределённости она не даёт.","Чего я не принял бы, так это отсутствия проверки вообще. У результата на одной сетке погрешность неизвестна, а неизвестная хуже большой.",[122,125,128],{"href":123,"label":124},"https:\u002F\u002Fdoi.org\u002F10.1115\u002F1.2960953","Celik et al. (2008), Journal of Fluids Engineering: procedure for estimation and reporting of uncertainty due to discretization in CFD applications",{"href":126,"label":127},"https:\u002F\u002Fdoi.org\u002F10.1115\u002F1.2910291","Roache (1994), Journal of Fluids Engineering: a method for uniform reporting of grid refinement studies",{"href":129,"label":130},"https:\u002F\u002Fwww.grc.nasa.gov\u002Fwww\u002Fwind\u002Fvalid\u002Ftutorial\u002Fspatconv.html","NASA \u002F NPARC: examining spatial (grid) convergence",1791759706180]