Асимптотическая структура пограничного слоя определяет сильно неравномерную сетку: используйте очень малые ячейки в направлении нормали к твёрдым стенкам, где градиенты скорости и скалярных величин велики, и постепенно увеличивайте размер ячеек во внешнем потоке, где градиенты слабее. Требуемое разрешение вблизи стенки зависит от решателя и способа обработки пристеночной области: DNS непосредственно разрешает вязкие и турбулентные масштабы, LES требует достаточного разрешения пристеночных структур или использования пристеночной модели, а RANS обычно ориентируется на заданную безразмерную координату стенки, например (y^+). На пилотной установке измеренные расходы, перепады давления, температуры и условия на стенке обеспечивают необходимые данные для валидации такой стратегии построения сетки.
При высоких числах Рейнольдса течение не является одинаково сложным для разрешения во всей области. Пограничный слой образует тонкую вязкую область, которая определяет касательное напряжение на стенке, потери давления, отрыв, а также тепло- или массообмен; поэтому при проектировании сетки следует распределять разрешение в соответствии с локальными масштабами пограничного слоя и проверять его по измерениям на пилотной установке.
Почему пограничный слой определяет проектирование сетки
Внешнему потоку не требуется равномерно мелкое разрешение
При высоких числах Рейнольдса течение часто можно концептуально разделить на внешнюю область, которая является приблизительно невязкой, и тонкий вязкий пограничный слой, прилегающий к стенке.
Градиенты скорости в направлении нормали к стенке во значительной части внешней области невелики. Равномерно мелкая сетка в этой области приводит к неоправданным затратам вычислительных ресурсов, не обязательно повышая точность расчёта интересующих параметров.
В пристеночной области сосредоточены наиболее важные градиенты
Условие прилипания требует, чтобы скорость достигала нуля на стенке. Затем она быстро возрастает до скорости внешнего потока в пределах пограничного слоя.
Это создаёт большие значения (\partial u/\partial y), которые непосредственно определяют касательное напряжение на стенке:
[ \tau_w = \mu \left.\frac{\partial u}{\partial y}\right|_w ]
Слишком грубая сетка в направлении нормали к стенке численно сглаживает этот градиент. Возникающие ошибки могут проявляться в виде неправильных коэффициентов трения, перепадов давления, сопротивления, распределений времени пребывания или коэффициентов теплоотдачи.
Толщина пограничного слоя уменьшается с ростом числа Рейнольдса
Характерный масштаб имеет вид
[ \frac{\delta}{L} \sim \frac{1}{\sqrt{Re}} ]
для ламинарного пограничного слоя при соответствующих условиях. Таким образом, увеличение числа Рейнольдса обычно делает пограничный слой более тонким относительно длины оборудования.
Практическое следствие очевидно: для моделирования течений с более высоким числом Рейнольдса обычно требуется более интенсивное уточнение сетки в направлении нормали к стенке, если только не используется соответствующая стратегия пристеночных функций или пристеночной модели.
Преобразование структуры пограничного слоя в сетку
В первую очередь уточняйте сетку по нормали к стенке
Первостепенное значение имеет шаг сетки в направлении нормали к стенке. Через пограничный слой должно проходить несколько ячеек, причём особенно тщательно их следует размещать вблизи стенки, где градиент скорости максимален.
В типичной пристеночной сетке используются слои с постепенным расширением или призматические слои, толщина которых плавно увеличивается по мере удаления от стенки. Резкие переходы от очень малых пристеночных ячеек к крупным внешним ячейкам могут снизить точность и вызвать проблемы со сходимостью.
Растягивайте сетку по мере удаления от стенки
Сетка должна плавно укрупняться от стенки к внешнему потоку. Это отражает уменьшение потребности в разрешении по мере приближения профиля скорости к значению скорости внешнего потока.
Коэффициент расширения должен оставаться достаточно умеренным, чтобы сохранять качество аппроксимации. Чрезмерное растяжение может привести к ошибкам интерполяции, неточным оценкам градиентов и повышенной чувствительности решателя.
Разрешайте развитие пограничного слоя вдоль потока
Пограничные слои увеличиваются вниз по потоку, поэтому требуемая протяжённость сетки по нормали к стенке и распределение ячеек изменяются вдоль оборудования.
Сетка в продольном направлении также должна разрешать:
- Развитие течения во входной области.
- Изменения градиента давления.
- Отрыв и повторное присоединение потока.
- Смешение, повороты, сужения и расширения.
- Локальные изменения тепло- или массообмена.
Сетка, спроектированная только на основе входного пограничного слоя, может оказаться недостаточной вблизи поворота, клапана, рабочего колеса, перегородки или резкого изменения площади сечения.
Добавляйте разрешение в поперечном направлении при наличии трёхмерных эффектов
В прямолинейных идеализированных течениях разрешение в поперечном направлении может быть менее требовательным, чем разрешение по нормали к стенке. Однако на пилотных установках трёхмерные эффекты встречаются часто.
Повороты, коллекторы, распределительные решётки, трубные пучки, вращающееся оборудование и несимметричные соединения могут создавать вторичные течения и поперечные градиенты. Для них требуется достаточное разрешение в поперечных направлениях, а не только опора на двумерное приближение пограничного слоя.
Согласуйте сетку с методом CFD
Моделирование RANS
Модели RANS обычно используют либо сетки с разрешением пристеночной области, либо пристеночные функции.
Сетка RANS с разрешением пристеночной области располагает первую ячейку достаточно близко к стенке, чтобы разрешить профиль скорости вблизи стенки. Сетка с пристеночными функциями располагает первую ячейку в соответствующем диапазоне логарифмического слоя и моделирует область, расположенную ближе к стенке.
Правильная цель заключается не просто в том, чтобы сделать сетку «как можно более мелкой». Необходимо использовать диапазон (y^+), соответствующий выбранной модели турбулентности и способу обработки пристеночной области. Сетка, не согласованная с этим способом, может привести к ненадёжным прогнозам касательного напряжения на стенке и теплоотдачи.
Моделирование LES
LES разрешает крупные турбулентные структуры и моделирует более мелкие. Вблизи твёрдых стенок разрешение наиболее мелких турбулентных масштабов становится особенно затратным.
Поэтому LES может потребовать:
- Очень высокого разрешения вблизи стенки для LES с разрешением пристеночной области.
- Валидированной пристеночной модели для крупномасштабных расчётов оборудования с высоким числом Рейнольдса.
Сетка также должна быть достаточно мелкой в продольном и поперечном направлениях, чтобы представлять разрешаемые турбулентные структуры. Одного уточнения только в направлении нормали к стенке недостаточно.
Моделирование DNS
DNS разрешает все значимые пространственные и временные масштабы без использования модели турбулентности. Требования к сетке значительно выше, особенно вблизи стенок и при больших числах Рейнольдса.
DNS обычно подходит для фундаментальных исследований или небольших типовых геометрий, но не для полноразмерного оборудования пилотных установок при промышленных режимах работы. Тем не менее DNS может предоставлять эталонные данные для валидации моделей пониженного порядка или пристеночных подходов.
Используйте измерения на пилотной установке для определения и валидации сетки
Начинайте с измеряемых рабочих условий
Проектирование сетки должно основываться на фактическом диапазоне рабочих режимов, а не только на номинальной геометрии. Важные исходные данные включают:
- Расход и скорость потока.
- Плотность, вязкость и теплофизические свойства жидкости.
- Уровни турбулентности или возмущений на входе.
- Температуру стенки или тепловой поток.
- Шероховатость поверхности.
- Рабочие давление и температуру.
- Переходные условия запуска или циклической работы.
Эти параметры определяют локальные числа Рейнольдса, толщины пограничных слоёв и области, в которых уточнение сетки наиболее важно.
Валидируйте величины, зависящие от разрешения у стенки
Данные пилотной установки следует сопоставлять с результатами моделирования, чувствительными к пограничному слою, включая:
- Перепад давления вдоль труб и оборудования.
- Коэффициент трения.
- Касательное напряжение на стенке.
- Профили скорости.
- Положения отрыва и повторного присоединения потока.
- Коэффициенты теплоотдачи.
- Профили концентрации или массообмена.
- Смешение и характеристики времени пребывания.
Одного совпадения общего расхода недостаточно, чтобы доказать адекватность пристеночной сетки.
Проводите исследования сходимости по сетке
Исследование сетки должно сравнивать как минимум три систематически связанных варианта, например грубую, среднюю и мелкую сетки. Сравнение следует сосредоточить на инженерных выходных параметрах, а не только на уменьшении невязок.
Расчёт может численно сойтись, сохраняя зависимость от сетки. Если прогнозы перепада давления, касательного напряжения на стенке или теплоотдачи существенно изменяются при уточнении, сетка ещё не демонстрирует независимость решения от сетки.
Учитывайте численную и экспериментальную неопределённость совместно
Измерения на пилотной установке содержат неопределённость, связанную с датчиками, калибровкой, возмущениями потока и изменчивостью рабочих условий. Численные результаты содержат неопределённость дискретизации, моделирования и исходных данных.
Правильный вопрос заключается не в том, совпадают ли расчёт и эксперимент в каждой точке. Следует определить, согласуется ли разница с совокупной неопределённостью и сохраняются ли одинаковые тенденции при различных рабочих режимах.
Учитывайте границы, отрыв и переходные процессы
Разделяйте области присоединённого и отрывного течения
Приближения пограничного слоя наиболее надёжны для преимущественно присоединённого течения при плавных изменениях геометрии. Сильные неблагоприятные градиенты давления могут вызвать отрыв, при котором пограничный слой теряет контакт со стенкой, а течение становится сильно нестационарным.
Зоны отрыва требуют локального уточнения более чем в одном направлении. След, слой сдвига, область рециркуляции и точка повторного присоединения могут определять потери давления и смешение.
Уточняйте сетку вблизи геометрических возмущений
Пилотные установки содержат элементы, делающие простую картину «тонкого слоя вдоль плоской стенки» неприменимой. Уточняйте сетку вблизи:
- Входов и выходов.
- Колен и тройников.
- Клапанов и соединительных элементов.
- Резких расширений и сужений.
- Перегородок и внутренних устройств.
- Входных участков труб.
- Рабочих колёс и вращающихся частей.
- Пористых распределителей.
Эти области создают сильные градиенты, вторичные течения или рециркуляцию, которые могут определять результат на уровне всего оборудования.
Разрешайте значимые временные масштабы
При запуске, остановке, пульсациях или периодической работе пограничный слой изменяется во времени. Шаг по времени должен разрешать соответствующий переходный процесс и быть согласованным с пространственной сеткой.
Безразмерные величины, такие как число Струхаля и число Рейнольдса, помогают определить, являются ли изменения процесса медленными по сравнению с развитием потока или достаточно быстрыми для формирования отдельной нестационарной реакции пограничного слоя.
Понимание компромиссов
Более мелкие сетки быстро увеличивают стоимость расчёта
Уточнение сетки вблизи стенки увеличивает количество ячеек, требования к памяти и вычислительное время. Трёхмерные пристеночные слои могут особенно сильно повышать стоимость расчёта в длинных трубах, теплообменниках и сложном технологическом оборудовании.
Поэтому целью должно быть целевое разрешение, а не максимальное разрешение во всей области.
Пристеночные функции снижают стоимость, но добавляют допущения
Пристеночные функции могут сделать расчёты при высоких числах Рейнольдса практически осуществимыми за счёт моделирования неразрешённой области вблизи стенки. Однако они менее надёжны при сильном отрыве, быстрых изменениях градиента давления, влиянии шероховатости, сильной кривизне, а также значительном тепло- и массообмене.
Их применение должно соответствовать положению первой ячейки сетки и физическому режиму.
Упрощения пограничного слоя применимы не всегда
В подходящих течениях с высоким числом Рейнольдса, преимущественно развивающихся вниз по потоку, приведённые уравнения пограничного слоя можно рассматривать как параболические. Это позволяет использовать эффективные пошаговые процедуры, поскольку информация распространяется главным образом вниз по потоку.
Однако там, где важны влияние вышележащих областей, рециркуляция, сильная трёхмерность, отрыв или эллиптическая связь по давлению, требуется полное описание уравнениями Навье—Стокса. Пилотную установку не следует упрощать только потому, что её общее число Рейнольдса велико.
Формальный порядок аппроксимации не гарантирует инженерную точность
Аппроксимация высокого порядка может уменьшить ошибку усечения, но не компенсирует неправильное расположение ячеек у стенки, неточное представление геометрии, неподходящее моделирование турбулентности или недостаточное разрешение отрыва.
Аналогично, добавление ячеек без проведения исследования сходимости не доказывает достоверность результата.
Как применить это к вашему проекту
Практический рабочий процесс заключается в проектировании сетки с учётом диапазона рабочих режимов и чувствительных к стенке выходных параметров, важных для конкретной технологической операции.
- Если основной интерес представляет перепад давления или касательное напряжение на стенке: уделите приоритетное внимание уточнению в направлении нормали к стенке, выберите согласованную стратегию по (y^+) и проведите валидацию по измеренным потерям давления и коэффициентам трения.
- Если основной интерес представляет тепло- или массообмен: разрешайте как гидродинамический, так и тепловой или концентрационный пограничный слой, поскольку их толщины могут различаться, и сравнивайте прогнозируемые коэффициенты переноса с данными пилотной установки.
- Если основной интерес представляет смешение или время пребывания: уточняйте сетку в следах, зонах рециркуляции, на входах, во внутренних устройствах и в трёхмерных областях поперечного течения, а не концентрируйте разрешение только у прямых стенок.
- Если основной интерес представляет отрыв или нестационарная работа: используйте локальное многомерное уточнение, достаточно малые шаги по времени и измерения, отражающие переходные изменения давления, скорости или поведения точки повторного присоединения.
- Если основной интерес представляет вычислительная эффективность: используйте слоистую локально уточнённую сетку и выбирайте RANS, LES или DNS в соответствии с требуемой физикой, а не применяйте самый затратный метод повсеместно.
- Если основной интерес представляет валидация модели: проводите исследования сходимости по сетке и сравнивайте несколько рабочих режимов, а не только одну точку данных пилотной установки.
Сетка, построенная с учётом асимптотической структуры пограничного слоя, концентрирует вычислительные ресурсы там, где физика процесса и инженерные решения наиболее чувствительны.
Сводная таблица:
| Приоритет при проектировании сетки | Рекомендуемый подход | Ключевые выходные параметры для валидации |
|---|---|---|
| Разрешение вблизи стенки | Используйте призматические или растущие слои и выбирайте пристеночную модель, согласованную с целевым значением y+ |
Касательное напряжение на стенке, коэффициент трения, перепад давления |
| Сетка во внешнем потоке | Постепенно увеличивайте размер ячеек там, где градиенты скорости и скалярных величин слабее | Общее развитие потока и эффективность |
| Геометрические возмущения | Уточняйте сетку в поворотах, клапанах, входах, выходах, внутренних устройствах, следах и зонах рециркуляции | Отрыв, смешение, время пребывания |
| Метод CFD | Согласуйте разрешение сетки с требованиями RANS, LES или DNS | Надёжность моделей турбулентности и пристеночных моделей |
| Тепло- и массообмен | Разрешайте гидродинамические, тепловые и концентрационные пограничные слои | Коэффициенты переноса и профили |
| Валидация | Сравнивайте грубую, среднюю и мелкую сетки с измерениями на пилотной установке | Независимость от сетки и согласованность неопределённостей |
Планируете надёжное CFD-исследование для вашего технологического оборудования? Свяжитесь с LABPARK, чтобы обсудить ваши требования. LABPARK предоставляет учебные и профессиональные пилотные установки для изучения единичных операций в химической инженерии, биопроцессах и биотехнологиях, а также в сфере экологической и водоочистной обработки. Наши пилотные установки помогают университетам, исследовательским институтам и предприятиям оценивать перепад давления, тепло- и массообмен, смешение, время пребывания, разделение и переходное поведение при измеряемых условиях эксплуатации. Мы поможем подобрать подходящее оборудование, определить рабочие условия, разработать стратегии построения сетки с учётом пограничного слоя и связать результаты моделирования с практической экспериментальной валидацией.