점근적 경계층 구조는 강력하게 불균일한 격자를 요구합니다: 속도 및 스칼라 구배가 가파른 고체 벽면 법선 방향으로는 매우 작은 셀을 사용하고, 구배가 약한 외부 유동 영역으로 갈수록 점진적으로 큰 셀을 배치하십시오. 필요한 벽면 근처 해상도는 솔버와 벽면 처리 방식에 따라 다릅니다. DNS는 점성 및 난류 스케일을 직접 해석하고, LES는 벽면 근처 구조물에 대한 적절한 해상도나 벽면 모델을 필요로 하며, RANS는 일반적으로 (y^+)와 같은 지정된 벽면 좌표를 목표로 합니다. 파일럿 플랜트에서는 측정된 유량, 압력 강하, 온도 및 벽면 조건이 해당 격자 전략을 검증하는 데 필요한 증거를 제공합니다.
높은 레이놀즈 수에서 유동은 균일하게 해석하기 어렵지 않습니다. 경계층은 벽면 전단, 압력 손실, 박리, 열 또는 물질 전달을 제어하는 얇은 점성 영역을 생성하므로, 격자 설계는 국부적인 경계층 스케일에 따라 해상도를 할당하고 파일럿 플랜트 측정값을 통해 이를 검증해야 합니다.
경계층이 격자 설계를 결정하는 이유
외부 유동에는 균일한 고해상도가 필요하지 않음
높은 레이놀즈 수에서 유동은 일반적으로 대략 비점성인 외부 영역과 벽면에 인접한 얇은 점성 경계층으로 개념적으로 분리될 수 있습니다.
벽면 법선 방향의 속도 구배는 외부 영역의 대부분에서 작습니다. 그곳에 균일하게 미세한 격자를 사용하는 것은 관심 있는 물리량의 개선 없이 계산 자원만 낭비하는 결과를 초래합니다.
벽면 영역에 중요한 구배가 존재함
점착 조건(no-slip condition)에 따라 벽면에서의 속도는 0이 되어야 하며, 이후 경계층을 가로질러 외부 유동 속도를 향해 급격히 상승합니다.
이는 벽면 전단 응력을 직접 결정하는 큰 (\partial u/\partial y) 값을 생성합니다:
[ \tau_w = \mu \left.\frac{\partial u}{\partial y}\right|_w ]
벽면 법선 방향으로 너무 거친 격자는 이 구배를 수치적으로 평탄화합니다. 이로 인한 오차는 잘못된 마찰 계수, 압력 강하, 항력, 체류 시간 분포 또는 열전달 계수로 나타날 수 있습니다.
경계층 두께는 레이놀즈 수에 따라 감소함
특성 스케일링은 다음과 같습니다.
[ \frac{\delta}{L} \sim \frac{1}{\sqrt{Re}} ]
(적절한 조건 하의 층류 경계층 기준). 따라서 레이놀즈 수가 증가하면 일반적으로 장비 길이에 비해 경계층이 더 얇아집니다.
실질적인 결과는 명확합니다. 적절한 벽면 함수나 벽면 모델 전략을 사용하지 않는 한, 더 높은 레이놀즈 수 시뮬레이션일수록 더 공격적인 벽면 법선 방향 세분화가 필요합니다.
구조를 격자로 변환하기
벽면 법선 방향부터 세분화
최우선 순위는 벽면 법선 방향 간격입니다. 경계층을 가로질러 여러 개의 셀이 배치되어야 하며, 속도 구배가 가장 큰 벽면 근처에는 특히 주의 깊게 배치해야 합니다.
일반적인 벽면 근처 메시는 벽면에서 멀어질수록 두께가 점진적으로 커지는 인플레이션(inflation) 또는 프리즘(prism) 층을 사용합니다. 아주 작은 벽면 셀과 큰 외부 셀 사이의 급격한 변화는 정확도를 떨어뜨리고 수렴 문제를 일으킬 수 있습니다.
벽면에서 멀어질수록 격자 확장
격자는 벽면에서 외부 유동을 향해 부드럽게 확장되어야 합니다. 이는 속도 프로파일이 외부 유동 값에 가까워짐에 따라 해상도 필요성이 감소함을 반영합니다.
확장 비율은 이산화 품질을 유지할 수 있을 만큼 적당해야 합니다. 과도한 확장은 보간 오차, 불량한 구배 추정 및 솔버 민감도를 유발할 수 있습니다.
유동 방향 발달 해석
경계층은 하류로 갈수록 성장하므로, 필요한 벽면 법선 방향 범위와 셀 분포는 장비를 따라 변화합니다.
유동 방향 격자는 다음을 해석해야 합니다:
- 입구 영역 발달.
- 압력 구배의 변화.
- 박리 및 재부착.
- 혼합, 굴곡부, 수축부 및 팽창부.
- 열 또는 물질 전달의 국부적 변화.
입구 경계층만을 기준으로 설계된 메시는 굴곡부, 밸브, 임펠러, 배플 또는 갑작스러운 면적 변화 근처에서 부적절해질 수 있습니다.
3차원 효과가 중요한 곳에 스팬 방향 해상도 추가
직선형의 이상적인 유동에서는 스팬 방향 해상도가 벽면 법선 방향보다 덜 중요할 수 있습니다. 그러나 파일럿 플랜트에서는 3차원 효과가 흔합니다.
굴곡부, 매니폴드, 분배판, 튜브 번들, 회전 장비 및 비대칭 피팅은 2차 유동과 교차 유동 구배를 생성할 수 있습니다. 이는 단순히 2차원 경계층 가정에 의존하기보다 횡단 방향으로 충분한 해상도를 요구합니다.
CFD 방법과 격자 맞추기
RANS 시뮬레이션
RANS 모델은 종종 벽면 해석(wall-resolved) 메시나 벽면 함수를 사용합니다.
벽면 해석 RANS 메시는 첫 번째 셀을 벽면 근처 속도 프로파일을 해석하기에 충분히 가깝게 배치합니다. 벽면 함수 메시는 첫 번째 셀을 적절한 로그 층(logarithmic-layer) 범위에 배치하고 벽면 근처 영역을 모델링합니다.
올바른 목표는 단순히 "가능한 한 미세하게"가 아닙니다. 선택한 난류 모델과 벽면 처리에 필요한 (y^+) 범위입니다. 해당 처리 방식과 일치하지 않는 메시는 신뢰할 수 없는 벽면 전단 및 열전달 예측을 생성할 수 있습니다.
LES 시뮬레이션
LES는 큰 난류 구조를 해석하고 작은 구조는 모델링합니다. 고체 벽면 근처에서는 가장 작은 난류 스케일을 해석하는 데 비용이 매우 많이 듭니다.
결과적으로 LES는 다음 중 하나가 필요할 수 있습니다:
- 벽면 해석 LES를 위한 매우 미세한 벽면 근처 해상도.
- 높은 레이놀즈 수의 장비 규모 시뮬레이션을 위한 검증된 벽면 모델.
격자는 또한 해석된 난류 구조를 표현하기 위해 유동 방향 및 스팬 방향으로 충분히 미세해야 합니다. 벽면 법선 방향만 세분화하는 것으로는 충분하지 않습니다.
DNS 시뮬레이션
DNS는 난류 모델 없이 모든 관련 공간 및 시간 스케일을 해석합니다. 격자 요구 사항이 훨씬 엄격하며, 특히 벽면 근처와 높은 레이놀즈 수에서 그렇습니다.
DNS는 일반적으로 기초 연구나 작은 대표 형상에 적합하며, 산업 운영 조건의 전체 규모 파일럿 플랜트 장비에는 적합하지 않습니다. 그럼에도 불구하고 축소 차수 모델이나 벽면 처리를 검증하기 위한 참조 데이터를 제공할 수는 있습니다.
파일럿 플랜트 측정값을 사용하여 메시 정의 및 검증
측정 가능한 운영 조건에서 시작
격자 설계는 명목상의 형상뿐만 아니라 실제 운영 범위를 기반으로 해야 합니다. 중요한 입력값은 다음과 같습니다:
- 유량 및 속도.
- 유체 밀도, 점도 및 열적 특성.
- 입구 난류 또는 교란 수준.
- 벽면 온도 또는 열유속.
- 표면 거칠기.
- 운영 압력 및 온도.
- 과도 상태의 시동 또는 순환 조건.
이러한 요소들이 국부적 레이놀즈 수, 경계층 두께 및 세분화가 가장 중요한 영역을 결정합니다.
벽면 해상도에 의존하는 물리량 검증
파일럿 플랜트 데이터는 경계층에 민감한 시뮬레이션 출력값과 비교해야 하며, 여기에는 다음이 포함됩니다:
- 배관 및 장비를 따른 압력 강하.
- 마찰 계수.
- 벽면 전단 응력.
- 속도 프로파일.
- 박리 및 재부착 위치.
- 열전달 계수.
- 농도 또는 물질 전달 프로파일.
- 혼합 및 체류 시간 거동.
전체 유량의 일치만으로는 벽면 근처 격자가 적절함을 입증할 수 없습니다.
격자 수렴 연구 사용
격자 연구는 거친(coarse), 중간(medium), 미세한(fine) 메시와 같이 체계적으로 관련된 최소 세 가지 메시를 비교해야 합니다. 비교는 잔차 감소뿐만 아니라 공학적 출력값에 초점을 맞춰야 합니다.
시뮬레이션이 수치적으로는 수렴하더라도 격자에 의존적일 수 있습니다. 압력 강하, 벽면 전단 또는 열전달 예측값이 세분화에 따라 크게 변한다면, 메시는 아직 솔루션 독립성을 입증하지 못한 것입니다.
수치적 및 실험적 불확실성 결합
파일럿 플랜트 측정값에는 센서, 교정, 유동 교란 및 운영 변동성으로 인한 불확실성이 포함되어 있습니다. 수치 결과에는 이산화, 모델링 및 입력 불확실성이 포함됩니다.
적절한 질문은 시뮬레이션과 실험이 점 대 점으로 동일한지가 아닙니다. 차이가 결합된 불확실성과 일치하는지, 그리고 동일한 경향이 운영 조건 전반에 걸쳐 지속되는지입니다.
경계, 박리 및 과도 상태 고려
부착 유동과 박리 유동 영역 분리
경계층 가정은 점진적인 기하학적 변화가 있는 주로 부착된 유동에서 가장 신뢰할 수 있습니다. 강한 역압력 구배는 박리를 유발할 수 있으며, 여기서 경계층은 벽면 부착력을 잃고 유동은 매우 불안정해집니다.
박리 영역은 하나 이상의 방향에서 국부적인 세분화가 필요합니다. 후류(wake), 전단층, 재순환 영역 및 재부착 지점은 모두 압력 손실과 혼합을 제어할 수 있습니다.
기하학적 교란 주변 세분화
파일럿 플랜트에는 "평평한 벽면을 따른 얇은 층"이라는 단순한 그림을 무효화하는 특징들이 포함되어 있습니다. 다음 주변을 세분화하십시오:
- 입구 및 출구.
- 엘보우 및 티(tee).
- 밸브 및 피팅.
- 급격한 팽창 및 수축부.
- 배플 및 내부 장치.
- 튜브 입구.
- 임펠러 및 회전 부품.
- 다공성 분배기.
이러한 영역은 장비 규모의 결과를 지배할 수 있는 강한 구배, 2차 유동 또는 재순환을 생성합니다.
관련 시간 스케일 해석
시동, 정지, 맥동 또는 주기적 운영의 경우 경계층은 시간에 따라 진화합니다. 시간 단계(time step)는 관련 과도 상태를 해석해야 하며 공간 메시와 호환되어야 합니다.
스트루할 수(Strouhal number) 및 레이놀즈 수와 같은 무차원량은 공정 변화가 유동 발달에 비해 느린지, 아니면 뚜렷한 비정상 경계층 반응을 생성할 만큼 빠른지 결정하는 데 도움이 됩니다.
트레이드오프 이해
미세한 격자는 비용을 급격히 증가시킴
벽면 근처 세분화는 셀 수, 메모리 요구량 및 계산 시간을 증가시킵니다. 3차원 인플레이션 층은 긴 배관, 열교환기 및 복잡한 공정 장비에서 이러한 비용을 특히 크게 만들 수 있습니다.
따라서 목표는 모든 곳의 최대 해상도가 아닌 표적 해상도(targeted resolution)입니다.
벽면 함수는 비용을 줄이지만 가정을 추가함
벽면 함수는 해석되지 않은 벽면 근처 영역을 모델링함으로써 높은 레이놀즈 수 시뮬레이션을 실용적으로 만들 수 있습니다. 유동이 강한 박리, 급격한 압력 구배 변화, 거칠기 효과, 강한 곡률 또는 상당한 열 및 물질 전달을 가질 때는 신뢰성이 떨어집니다.
사용 시 메시 위치 및 물리적 영역과 일치해야 합니다.
경계층 단순화는 보편적으로 유효하지 않음
축소된 경계층 방정식은 적절한 높은 레이놀즈 수의 주로 하류로 흐르는 유동에서 포물선형으로 취급될 수 있습니다. 이는 정보가 주로 하류로 전파되기 때문에 효율적인 행진 절차(marching procedures)를 가능하게 합니다.
그러나 상류 영향, 재순환, 강한 3차원성, 박리 또는 타원형 압력 결합이 중요한 곳에서는 전체 나비에-스토크스 처리가 필요합니다. 파일럿 플랜트를 단순히 전체 레이놀즈 수가 높다는 이유만으로 단순화해서는 안 됩니다.
형식적 차수가 공학적 정확도를 보장하지 않음
고차 이산화는 절단 오차를 줄일 수 있지만, 부적절한 벽면 배치, 불량한 형상 표현, 부적절한 난류 모델링 또는 박리에 대한 불충분한 해상도를 보상할 수는 없습니다.
마찬가지로 수렴 연구를 수행하지 않고 셀을 추가하는 것만으로는 결과가 신뢰할 수 있음을 확립할 수 없습니다.
프로젝트에 적용하는 방법
실용적인 워크플로우는 운영 범위와 단위 공정에 중요한 벽면 민감 출력값을 중심으로 격자를 설계하는 것입니다.
- 주요 초점이 압력 강하나 벽면 전단인 경우: 벽면 법선 방향 세분화를 우선시하고 일관된 (y^+) 전략을 선택하며, 측정된 압력 손실 및 마찰 계수와 비교하여 검증하십시오.
- 주요 초점이 열 또는 물질 전달인 경우: 두께가 다를 수 있는 유체 역학적 경계층과 열적 또는 농도 경계층을 모두 해석하고, 예측된 전달 계수를 파일럿 플랜트 데이터와 비교하십시오.
- 주요 초점이 혼합이나 체류 시간인 경우: 직선 벽면에만 해상도를 집중하기보다 후류, 재순환 영역, 입구, 내부 장치 및 3차원 교차 유동 영역을 세분화하십시오.
- 주요 초점이 박리나 비정상 운영인 경우: 국부적인 다방향 세분화, 충분히 작은 시간 단계, 그리고 과도 상태의 압력, 속도 또는 재부착 거동을 포착하는 측정값을 사용하십시오.
- 주요 초점이 계산 효율성인 경우: 계층화된 국부 세분화 메시를 사용하고, 모든 곳에 가장 비싼 방법을 적용하기보다 필요한 물리학에 따라 RANS, LES 또는 DNS를 선택하십시오.
- 주요 초점이 모델 검증인 경우: 단일 파일럿 플랜트 데이터 포인트가 아닌, 메시 수렴 연구를 수행하고 여러 운영 조건을 비교하십시오.
점근적 경계층 구조에 기반한 격자는 물리학과 공학적 결정이 가장 민감한 곳에 계산 노력을 집중시킵니다.
요약 표:
| 격자 설계 우선순위 | 권장 접근 방식 | 검증할 주요 출력값 |
|---|---|---|
| 벽면 근처 해상도 | 프리즘 또는 인플레이션 층을 사용하고 목표 y+와 일치하는 벽면 처리 선택 |
벽면 전단, 마찰 계수, 압력 강하 |
| 외부 유동 메시 | 속도 및 스칼라 구배가 약한 곳에서 셀 크기를 점진적으로 증가 | 전체 유동 발달 및 효율 |
| 기하학적 교란 | 굴곡부, 밸브, 입구, 출구, 내부 장치, 후류 및 재순환 영역 세분화 | 박리, 혼합, 체류 시간 |
| CFD 방법 | 메시 해상도를 RANS, LES 또는 DNS 요구 사항에 맞춤 | 난류 및 벽면 모델 신뢰성 |
| 열 및 물질 전달 | 유체 역학적, 열적 및 농도 경계층 해석 | 전달 계수 및 프로파일 |
| 검증 | 거친, 중간, 미세 메시를 파일럿 플랜트 측정값과 비교 | 격자 독립성 및 불확실성 일치 |
공정 장비를 위한 신뢰할 수 있는 CFD 연구를 계획 중이십니까? 요구 사항을 논의하려면 LABPARK에 문의하십시오. LABPARK는 화학 공학, 바이오 공정 및 생명 공학, 환경 및 수처리 응용 분야를 위한 교육 및 직업 훈련용 단위 조작 파일럿 플랜트를 제공합니다. 당사의 파일럿 플랜트는 대학, 연구소 및 기업이 측정 가능한 조건 하에서 압력 강하, 열 및 물질 전달, 혼합, 체류 시간, 박리 및 과도 운영 거동을 평가하도록 돕습니다. 관련 장비 선택, 운영 조건 정의, 경계층 기반 격자 전략 지원, 그리고 시뮬레이션 결과와 실제 실험 검증을 연결하는 것을 도와드릴 수 있습니다.