파티클 해석 : 흐름을 따르지 않는 입자를 해석하는 방법

1. 서론

현대 전산유체역학(CFD)에서 유체는 일반적으로 공간에 고정된 요소를 기준으로 물리량을 계산하는 오일러(Eulerian) 접근법을 통해 연속체(Continuum)로 정의됩니다. 그러나 실제 산업 현장의 유동장 내에는 고체 분진과 같은 이산체(Discrete Phase)가 빈번하게 혼재되어 있습니다. 이러한 이산체들은 유체와 독립적인 질량과 관성을 지니고 있어, 유체의 흐름에 온전히 순응하지 않고 궤적을 이탈하는 비평형 거동을 보입니다.

이처럼 개별적인 파티클의 거동을 정밀하게 모사하기 위해, 본 자료에서 다루는 '파티클 해석'은, 유동장은 오일러 방식으로 처리하되 파티클은 뉴턴의 제2법칙에 기반하여 개별적으로 추적하는 오일러-라그랑주(Euler-Lagrange) 접근법을 채택합니다. 이는 연속체 가정을 보완하여 유체와 파티클 간의 운동량 교환 및 에너지 전달을 수치적으로 구현하는 핵심 방법론입니다.

1.1. 연속체와 이산체의 공존
: 오일러(Eulerian)와 라그랑주(Lagrangian) 좌표계

오일러-라그랑주 접근법은 파티클 해석을 수행하는 방법론입니다. 파티클의 거동을 정확히 해석하기 위해서는 공간과 개체를 바라보는 두 가지 물리적 관점의 차이를 명확히 이해하는 것이 좋습니다.

• 1.1.1. 오일러(Eulerian) 관점

공간 중심의 기술 (연속체) : 

강가에서 서서 특정 지점을 통과하는 물의 속도와 압력을 측정하는 방식에 비유할 수 있습니다. CFD에서는 고정된 계산 요소를 기준으로 나비에-스토크스(Navier-Stokes) 방정식을 풀어 전체적인 유동장을 도출합니다. 이 유동장은 파티클이 이동할 수 있는 공간적 배경을 제공하며, 파티클에 작용하는 항력(Drag Force) 등의 외력을 결정하는 기반이 됩니다.

• 1.1.2. 라그랑주(Lagrangian) 관점

개체 중심의 기술 (이산체) : 

강물에 띄운 나뭇잎 하나를 끝까지 쫓아가며 그 위치와 속도를 관찰하는 방식입니다. 파티클 해석에서 솔버는 파티클을 물리적 부피가 아닌 질량점(Point Mass)으로 인식합니다. 오일러 유동장에서 제공받은 정보를 바탕으로 파티클에 작용하는 외력을 뉴턴의 운동 제2법칙(F= ma)에 대입하여 개별 파티클의 궤적을 시간 적분으로 추적합니다.

1.2. 가정 및 유효 한계

오일러-라그랑주 모델링은 매우 유용하지만, 모든 조건에서 완벽하게 작동하는 것은 아닙니다. 해석의 신뢰성을 위해 솔버가 전제로 하는 다음의 공학적 가정을 인지하는 것이 좋습니다.

• 1.2.1. 구형 파티클 가정

연산의 효율성을 위해 모든 파티클은 기본적으로 완벽한 구형으로 간주됩니다. 불규칙한 형상의 파티클을 해석할 경우, 형상 계수(Shape Factor)를 적용하여 항력 계수(Drag Coefficient)를 보정을 통해 물리적 오차를 줄일 수 있습니다.

• 1.2.2. 저 농도 분산 유동 한계

본 모델은 파티클이 차지하는 체적 점유율이 통상 0.1%~5% 미만의 저농도 분산 유동에서 물리적 타당성을 갖습니다. 여기서 '희박(Dilute)'이라는 용어는 압력이 낮은 '희박 기체(Rarefied Gas)'와는 무관하며, 파티클 간의 직접 충돌을 무시할 수 있을 만큼 농도가 낮다는 것을 의미합니다. 만약 파티클 농도가 높아 충돌이 지배적인 환경이라면 DEM(Discrete Element Method)과 같은 접촉 역학 기반의 수치 기법을 병행하는 것을 권장합니다. 

2. 파티클 해석의 수치적 연동 체계와 물리적 매커니즘

제 1장에서 두 좌표계의 개념을 정립했다면, 제 2장에서는 실제 midas NFX 솔버 내부에서 이들이 어떻게 데이터를 주고 받으며 물리적인 위치를 확정하는지 그 연동 메커니즘을 상세히 분석합니다.

2.1. 오일러리안과 라그랑지안의 상호작용 3단계

파티클 해석은 공간 정보를 기반으로 개체의 이동을 계산하는 일련의 수치해석적 방식을 따릅니다. 솔버는 다음의 3단계를 거쳐 하나의 파티클 궤적을 완성합니다.

[1단계] Eulerian 좌표계 : 유동장의 물리정보 (환경 정보)

솔버는 요소 내에서 유체의 속도, 압력, 밀도 등을 먼저 계산합니다. 이는 파티클이 비행하게 될 유동장을 만다는 과저이며, 파티클을 밀어내는 힘인 항력의 기초 데이터가 됩니다.

[2단계] 물리적 보간(Interpolation) : 데이터의 가교

유동 정보는 요소의 절점에 존재하지만, 파티클은 요소 내부의 임의 좌표(x, y, z)에 위치합니다. 솔버는 요소의 절점의 정보를 읽어와 파티클의 현재 위치에서의 유동 속도를 실시간으로 추산합니다. 이를 통해 파티클이 이동하는 유동의 속도를 인식합니다.  

 😥 미리 보기는 여기까지!
내용을 이어서 보고 싶다면,
아래 정보를 입력해 주세요.