유한체적법 전산유체 (Finite Volume Method in CFD)
쉽게 풀면
큰 방 안의 온도 변화를 계산한다고 상상해봅시다. 방 전체를 한 번에 계산하기는 어려우니, 방을 여러 개의 작은 상자로 나누고 각 상자의 벽을 통해 얼마나 열이 들어오고 나가는지를 따로 계산합니다. 유한체적법은 바로 이런 방식으로 유체가 흐르는 공간을 작은 체적들로 나누고, 각 체적의 경계면을 통과하는 질량, 운동량, 에너지의 흐름을 계산합니다. 이 방법은 각 상자에 들어오고 나간 양이 정확히 맞아떨어지도록 설계되어 있어서, 질량이나 에너지가 계산 중에 사라지거나 생겨나지 않는다는 장점이 있습니다.
왜 중요한가
유한체적법은 질량, 운동량, 에너지 보존을 자연스럽게 만족시키는 구조 덕분에 상용 및 오픈소스 CFD 소프트웨어 대부분이 채택하고 있는 표준적인 이산화 방법입니다. 복잡한 형상과 비구조 격자에도 적용이 가능해 실제 공학 문제에 폭넓게 쓰입니다.
논문에서는 이렇게 쓰입니다
수치해석 절차에서 유한체적법이 기본 이산화 기법으로 채택된 예문입니다.
복잡한 형상 해석에 유한체적법의 유연성이 활용된 예문입니다.
조금 더 깊게 보면
유한체적법은 지배방정식의 적분형을 각 제어체적에 적용하고, 발산 정리를 이용해 체적 내부의 변화율을 경계면을 통과하는 플럭스의 합으로 표현합니다. 대류항과 확산항을 이산화하는 방식(상류차분, 중심차분 등)과 시간 적분 기법의 선택에 따라 정확도와 안정성이 달라지며, 압력-속도 연성을 위한 SIMPLE 계열 알고리즘이 함께 사용되는 경우가 많습니다.
주의할 점
격자 품질이 낮거나 왜곡이 심한 경우 수치 확산이나 오차가 커질 수 있으므로, 유한체적법 자체의 보존 특성만으로 해석 정확도가 자동으로 보장되는 것은 아닙니다.