Development and application of a multiphase Lagrangian structure function model in anisotropic turbulence

기술 요약: 비등방성 난류에서의 다상 라그랑지안 구조 함수 모델의 개발 및 적용

문제 제기
관성 입자가 난류 흐름에 반응하는 에너지적 응답은 다상 역학을 이해하고 입자 분산을 위한 차수 축소 모델을 개발하는 데 매우 중요하다. 많은 연구가 균질 등방성 난류(HIT)에서의 입자 에너지론을 조사해 왔으나, 실제 흐름은 특히 고체 경계 근처에서 공간적 불균질성과 비등방성을 나타내는 경우가 많다. 난류 경계층(TBL)과 같은 벽면 결합 난류 흐름에서는 평균 전단(mean shear)과 경계 조건에 의해 대규모 비등방성이 생성된다. 핵심적인 미해결 과제는 이러한 대규모 비등병성이 입자상으로 어떻게 전달되는지, 그리고 이것이 형식적으로는 등방성인 소규모 영역까지 지속되는지 여부이다. 또한, 이러한 환경에서 입자의 관성, 중력 침강, 그리고 흐름의 비등방성 사이의 복잡한 결합 관계를 정확하게 모델링하기 위해서는 추가적인 조사가 필요하다.

연구 방법론
저자들은 순 침강(net deposition) 조건 하의 난류 경계층 내에서 스케일별 입자 에너지론을 조사하기 위해 이론적 및 수치적 접근 방식을 결합하여 사용하였다.

1. 이론적 프레임워크:

   • 본 연구는 임의의 스케일에서 입자상의 2차 구조 함수(SoSF)에 대한 정확하지만 닫히지 않은(unclosed) 보존 법칙을 유도하기 위해 운동론적 접근 방식을 활용한다.
   • 입자 역학은 입자-유체 밀도비가 큰 극한에서의 작고 구형이며 상호작용하지 않는 입자에 대한 Maxey-Riley 방정식을 사용하여 모델링된다.
   • 저자들은 종방향 및 수직 2차 구조 함수의 차이에 기초하여 입자상의 비등방성 지표인 $\eta^{(p)}$를 정의한다.
   • 에너지 플럭스와 비등방성을 지배하는 항들의 균형을 규명하기 위해, 상대 속도 모멘트와 반경 방향 분포 함수의 멱법칙(power-law) 거동을 가정하여 유도된 연속 방정식에 대한 점근적 분석을 수행한다.

2. 수치 시뮬레이션:

   • 일방향 결합(one-way coupled) 입자 함유 벽면 결합 난류의 직접 수치 시뮬레이션(DNS)이 중간 마찰 레이놀즈 수($Re_\tau = 315$)에서 수행되었다.
   • 시뮬레이션은 TBL의 로그층과 버퍼층을 다룬다.
   • 입자의 관성을 평가하기 위해 스토크스 수($St^+$)를, 중력 표류의 영향을 평가하기 위해 침강 파라미터($Sv^+$)를 변화시킨다.
   • 오일러적(유체) 관점과 라그랑지안(입자) 관점이 모두 활용되며, 특히 국부적인 수직 역학을 격리하기 위해 수직 분리가 zero인 "슬래브 단위(slab-wise)" 통계에 초점을 맞춘다.

주요 기여

• 보존 법칙 유도: 본 논문은 수송, 1차 모멘트와의 결합, 그리고 에너지 플럭스를 위한 항들을 명시적으로 식별함으로써, 임의의 스케일에 대한 입자상의 에너지론에 관한 정확한 보존 방정식을 유도하였다. 이는 비등방성 다상 흐름에서 스케일별 에너지 전달을 분석하기 위한 연속체 프레임워크를 제공한다.
• 비등방성 전달의 정량화: 본 연구는 대규모 유체의 비등방성이 입자상으로 어떻게 "기여(donated)"되는지, 그리고 입자 관성과 침강이 이러한 전달을 어떻게 수정하는지를 체계적으로 정량화한다.
• 점근적 분석: 저자들은 DNS 데이터를 해석하기 위해 점근적 논거를 제공하며, 이를 통해 관찰된 비등방성 지표의 비단조적(non-monotonic) 거동을 대규모 흐름의 비등방성과 소규모 간 스케일 에너지 전달 사이의 경쟁 관계와 연결시킨다.

결과

• 관성과 비등방성: 작은 스토크스 수에서 입자상은 유체를 모사하여 소규모에서 거의 등방성을 보인다. 그러나 관성이 증가함에 따라($St^+ = 10$), 입자상은 소규모에서도 상당한 비등방성을 유지하며, 이는 입자 관성이 대규모 흐름 비등방성의 기억을 보존함을 나타낸다.
• 침강 효과: 침강 파라미터($Sv^+$)의 영향은 유동 영역에 따라 크게 달라진다. 유동 자체가 본질적으로 비등방적인 버퍼층($z^+ = 30$)에서는 $Sv^+$가 증가함에 따라 수직 성분은 증가하고 종방향 성분은 감소하는 복잡한 구조 함수 변화가 나타난다. 소규모 흐름이 등방적인 로그층($z^+ = 90$)에서는 높은 침강률($Sv^+$$)이 입자상을 소규모에서 등방적 상태로 유도하는 경향이 있다.
• 스케일 비국소성(Scale Non-Locality): 입자 구조 함수는 유체의 기대되는 $r^2$ 스케일링에 비해 소산 범위에서 더 완만한 $r^{2/3}$ 스케일링에 접근하는 경shallow한 스케일링을 보인다. 이는 관성과 침강에 의해 주로 구동되는 관성 입자의 스케일 비국소성을 반영한다.
• 비등방성 지표 거동: 입자 비등방성 지표 $\eta^{(p)}$는 분리 스케일 $r$이 감소함에 따라 비단조적 거동을 보인다. $r/z \to 0$일 때 다시 증가하기 전에 국부 최솟값에 도달한다. 이는 대규모 흐름에 의해 부과된 비등방성과 입자상 내부의 간 스케일 에너지 전달에 의해 생성된 비등방성 사이의 경쟁을 시사한다.

의의 및 주장
본 논문은 비등방성 난류 내 관성 입자의 스케일별 에너지 플럭스에 대한 중요한 통찰을 제공한다고 주장한다. 입자 관성이 형식적으로 등방적인 소규모 영역까지 대규모 비등방성을 보존할 수 있음을 입증함으로써, 저자들은 입자상의 국부 등방성을 가정하는 모델들의 한계를 지적한다. 유도된 보존 법칙과 확인된 비등방성 전달 메커니즘은 관성 침강 입자에 대한 연속체 모델 개발, 특히 이러한 비국소적 비등방성 효과가 이방향 결합(two-way coupled) 역학에 어떻게 영향을 미칠지에 대해 직접적인 시사점을 준다. 본 연구는 정확한 모델링을 위해 입자 관성과 중력 표류를 동시에 고려하는 것이 필수적임을 강조하는데, 이는 두 효과가 단순히 가산적인 것이 아니라 입자 역학 내에서 암묵적으로 결합되어 있기 때문이다.