19.27 격자 생성과 수치 해법

전산 유체역학에서 격자(grid 또는 mesh)는 연속적인 해석 도메인을 이산 요소의 집합으로 분할한 것이며, 지배 방정식의 이산화와 해의 정확도를 결정짓는 핵심 요소이다. 격자 생성(grid generation)은 기하학적 형상을 표현하는 작업인 동시에, 해석의 수렴성과 계산 효율성을 좌우하는 수치 해법의 사전 단계에 해당한다. 본 절에서는 격자의 분류, 생성 기법, 품질 지표, 그리고 이산 방정식의 수치 해법을 체계적으로 정리한다.

1. 격자의 분류

격자는 위상 구조에 따라 구조화 격자(structured grid)와 비구조 격자(unstructured grid), 그리고 이들의 장점을 결합한 혼합 격자(hybrid grid)로 분류된다.

구조화 격자는 (i, j, k)와 같은 논리적 인덱스를 이용하여 각 절점을 고유하게 식별할 수 있는 격자이며, 인접 정보가 암시적으로 주어져 메모리 효율이 높고 해법의 구현이 단순하다. 그러나 복잡한 형상에 대한 대응이 제한적이라는 단점이 있다. 비구조 격자는 삼각형, 사면체, 다면체 등 임의 형상의 요소로 구성되며, 인접 정보를 명시적으로 저장해야 한다. 복잡한 기하학적 형상에 자유롭게 대응할 수 있는 유연성이 가장 큰 장점이다.

멀티 블록 구조화 격자(multi-block structured grid)는 해석 도메인을 여러 블록으로 분할하고 각 블록 내부에는 구조화 격자를 구성하는 방식이다. 이를 통해 구조화 격자의 효율성과 복잡한 형상 대응 능력을 동시에 확보할 수 있다. 오버셋 격자(overset 또는 Chimera grid)는 독립적으로 생성된 격자 블록들이 서로 중첩되어 해석 도메인을 구성하며, 상대 운동이 있는 다중 물체 해석에 적합하다.

2. 구조화 격자 생성 기법

구조화 격자 생성의 대표적 기법으로는 대수적 방법(algebraic method), 편미분 방정식 기반 방법, 변분적 방법이 있다. 대수적 방법은 경계에서의 점 분포를 내부로 직접 보간하는 방식이며, 대표적으로 초한 보간(transfinite interpolation, TFI)이 사용된다.

편미분 방정식 기반 기법에서는 격자 좌표를 라플라스 또는 푸아송 방정식의 해로 얻는다. Thompson, Thames, Mastin이 제안한 TTM 방식은 타원형 방정식을 이용하여 매끄럽고 직교성이 우수한 격자를 생성한다. 쌍곡형 방법은 경계에서 출발하여 외부 방향으로 격자를 전진적으로 생성하며, 외부 유동 해석에 유용하다.

3. 비구조 격자 생성 기법

비구조 격자 생성의 핵심 알고리즘은 두 가지로 구분된다. 전진 경계법(Advancing Front Method, AFM)은 경계면에서 출발하여 내부 방향으로 요소를 순차적으로 생성하는 방식이며, 경계 인접 영역의 품질이 우수하다. Delaunay 삼각화(Delaunay triangulation)는 주어진 점 집합에 대해 최소각을 최대화하는 삼각화를 생성하는 기법으로, Bowyer-Watson 알고리즘이 대표적이다.

2차원 삼각형 격자는 3차원 사면체 격자로 확장되며, 최근에는 6면체와 다면체 격자 생성 기법도 발달하였다. 다면체 격자는 셀당 이웃 수가 많아 구배 계산의 정확도가 높고 수렴성이 우수한 장점이 있다.

4. 경계층 격자

점성 유동 해석에서는 벽면 근처의 경계층을 정확히 해석하기 위해 조밀한 격자가 요구된다. 이를 위해 벽면에 평행한 방향으로는 성긴 격자를, 벽면에 수직인 방향으로는 조밀한 격자를 생성하는 경계층 격자(prism layer 또는 boundary layer mesh)가 사용된다. 첫 격자점의 무차원 거리 y^+는 난류 모델의 요구 사항에 따라 결정되며, 저레이놀즈 해법의 경우 y^+ < 1을, 벽 함수 해법의 경우 30 < y^+ < 300 범위가 권장된다.

5. 격자 품질 지표

격자의 수치적 품질을 정량적으로 평가하기 위한 지표로는 다음이 널리 사용된다.

이들 품질 지표가 악화되면 이산화 오차가 커지고 반복 해법의 수렴성이 저하된다. 격자 생성 소프트웨어는 품질 지표의 분포를 시각화하여 부적절한 영역을 식별하는 기능을 제공한다.

6. 적응 격자 세분화

해의 구배가 큰 영역이나 충격파, 와류, 경계층과 같은 국소 특징이 존재하는 영역에서는 격자를 집중적으로 세분화해야 해석 정확도가 확보된다. 적응 격자 세분화(Adaptive Mesh Refinement, AMR)는 해 기반의 오차 지표를 사용하여 자동으로 격자를 세분화 또는 성기화하는 기법이다. 구조화 격자에 대해서는 Berger-Oliger 방식의 블록 구조 AMR이, 비구조 격자에 대해서는 셀 분할 기반의 h-세분화(h-refinement)가 주로 사용된다.

7. 대수 방정식 체계와 직접 해법

이산화 결과 얻어지는 대수 방정식 체계는 일반적으로 \mathbf{A}\mathbf{x} = \mathbf{b}의 형태를 가지며, 행렬 \mathbf{A}는 희소 행렬(sparse matrix)이다. 격자 크기 N이 수백만에서 수억에 이르는 현대의 CFD 문제에서는 행렬의 크기가 막대하므로 적절한 해법의 선택이 필수적이다.

직접 해법(direct solver)은 가우스 소거법이나 LU 분해를 통해 엄밀한 해를 구하는 방식이며, 소규모 문제나 잘 정제된 희소 행렬에 대해 효과적이다. 그러나 대규모 3차원 CFD 문제에서는 메모리 요구량과 계산 비용이 과도하여 일반적으로 반복 해법이 선호된다.

8. 반복 해법과 전처리

반복 해법(iterative solver)은 초기 추정해에서 출발하여 수렴 기준을 만족할 때까지 해를 점진적으로 갱신한다. 고전적 반복법으로는 Jacobi, Gauss-Seidel, Successive Over-Relaxation(SOR) 등이 있으며, Krylov 부분공간 방법으로는 Conjugate Gradient(CG), Generalized Minimal Residual(GMRES), BiCGSTAB 등이 널리 사용된다.

반복 해법의 수렴 속도를 가속하기 위한 전처리(preconditioning)는 필수적이며, 대표적 전처리 기법으로는 불완전 LU 분해(Incomplete LU, ILU), Jacobi 전처리, Algebraic Multigrid(AMG) 등이 있다. 이 중 AMG는 대규모 CFD 문제에서 격자 크기에 거의 무관한 수렴성을 제공하므로 현대 상용 해석 코드에서 표준 전처리로 채택되고 있다.

9. 멀티그리드 기법

멀티그리드(multigrid) 기법은 서로 다른 공간 해상도의 격자 위에서 교대로 해를 갱신함으로써, 저주파 오차와 고주파 오차를 효율적으로 감쇠시키는 해법이다. 기하학적 멀티그리드(geometric multigrid)는 실제 격자 계층을 이용하며, 대수적 멀티그리드(algebraic multigrid)는 행렬 구조만을 이용하여 가상의 격자 계층을 구성한다. V-사이클, W-사이클, F-사이클 등 다양한 사이클 전략이 존재하며, CFD의 수렴성 향상에 결정적 기여를 한다.

10. 비선형 해법과 의사 시간 적분

나비에-스토크스 방정식은 비선형이므로, 각 반복 단계에서 비선형 항을 선형화하여 해를 구한 뒤 갱신하는 방식이 사용된다. 뉴턴 방법과 그 변형인 Newton-Krylov 기법은 2차 수렴 속도를 제공하나 초기 추정해에 민감하다. 정상 상태 해석에서는 가상의 시간 미분 항을 추가하여 비정상 해석과 유사한 방식으로 수렴시키는 의사 시간 적분(pseudo-time stepping) 기법이 널리 사용된다.

11. 수렴 판정과 잔차 모니터링

해법의 수렴 여부는 방정식의 잔차(residual) 감소량과 관심 물리량(예: 양력, 항력, 질량 유량)의 안정화 여부를 기준으로 판정한다. 잔차는 대개 L^2 노름 또는 최댓값 노름으로 정의되며, 초기 잔차 대비 10^{-4}에서 10^{-6} 이하로 감소할 때 수렴으로 판단하는 것이 일반적이다. 비정상 해석에서는 각 시간 단계 내부의 반복 수렴과 전체 시간 적분의 안정성을 함께 점검해야 한다.

12. 로봇 공학에서의 응용

격자 생성과 수치 해법은 로봇 공학의 유체역학 해석 전반에 직접적으로 적용된다. 무인 항공기의 공력 해석에서는 동체와 날개의 복잡한 형상을 표현하는 혼합 격자가 사용되며, 경계층 해석을 위한 조밀 격자가 필수적이다. 회전 로터의 해석에는 오버셋 격자와 슬라이딩 메시(sliding mesh) 기법이 활용된다. 수중 로봇의 프로펠러 해석에서는 캐비테이션 모델링과 결합된 비정상 반복 해법이 요구되며, 소프트 로봇과 유체-구조 상호작용 해석에서는 격자 변형과 재격자화 기법이 핵심적으로 작용한다. 본 절에서 정리한 원리는 이러한 응용의 이론적 토대를 제공한다.

13. 출처

• Thompson, J. F., Soni, B. K., and Weatherill, N. P., eds., Handbook of Grid Generation, CRC Press, 1999.
• Frey, P. J., and George, P. L., Mesh Generation: Application to Finite Elements, 2nd ed., Wiley-ISTE, 2008.
• Ferziger, J. H., Perić, M., and Street, R. L., Computational Methods for Fluid Dynamics, 4th ed., Springer, 2020.
• Saad, Y., Iterative Methods for Sparse Linear Systems, 2nd ed., SIAM, 2003.
• Trottenberg, U., Oosterlee, C. W., and Schüller, A., Multigrid, Academic Press, 2001.
• Versteeg, H. K., and Malalasekera, W., An Introduction to Computational Fluid Dynamics: The Finite Volume Method, 2nd ed., Pearson, 2007.

14. 버전

v1.0