Inertial effects on the interphase drag force and rheology of dilute suspensions of buoyant droplets at low Reynolds number

본 연구는 상호 관계 정리(reciprocal theorem)를 사용하여 저레이놀즈수에서 부력이 있는 부유 액적의 희박 현탁액 내 관성 효과가 계면 항력 및 연속상의 유효 응력에 상대 속도와 속도 분산에 대한 이차 의존성을 도입함을 입증한다.

핵심

모두가 같은 방향으로 걸어가려 노력하는 북적이는 방을 상상해 보세요. 하지만 어떤 사람들은 떠다니는 풍선(부유 액적)이고, 다른 사람들은 공기(연속체 유체)입니다. 보통 우리가 이 풍선들이 공기 속에서 어떻게 움직이는지 연구할 때는, 공기가 완벽하게 정지해 있고 풍선이 아주 느리게 움직인다고 가정합니다. 마치 걸쭉한 시럽 속을 움직이는 달팽리처럼 말이죠. 이 느린 세상에서는 규칙이 단순합니다. 풍선이 더 빨리 움직일수록 공기가 더 강하게 밀어냅니다.

하지만 현실 세계에서는 항상 이렇게 느리거나 단순하지는 않습니다. 때때로 공기는 약간의 "추진력(관성)"을 가지고 있으며, 풍선들은 단순히 직선으로 이동하는 것이 아니라 이리저리 흔들릴 수도 있습니다. 니콜라스 핀치(Nicolas Fintzi)와 장 루 피에르송(Jean-Lou Pierson)이 작성한 이 논문은 다음과 같은 구체적인 질문을 던집니다. 공기에 약간의 속도가 있고, 풍선들이 자신만의 작은 에너지를 가지고 이리저리 튀어 오를 때, 이 떠다니는 풍선들에 가해지는 힘에는 어떤 일이 일어날까?

다음은 일상적인 비유를 사용한 그들의 발견에 대한 설명입니다.

1. "상반 정리(Reciprocal Theorem)"는 마법의 거울

이를 해결하기 위해 저자들은 공기 입자 하나하나와 풍선 하나하나를 모두 시뮬레이션하는 방식을 사용하지 않았습니다. 그것은 해변의 파도가 어떻게 움직이는지 이해하기 위해 모래알 하나하나를 세는 것과 같습니다. 대신, 그들은 상반 정리라는 수학적 도구를 사용했습니다.

이것은 마법의 거울과 같습니다. 바람이 약간 부는 환경에서 움직이는 풍선의 복잡하고 무질서한 현실을 직접 들여다보는 대신, 그들은 규칙이 더 단순한(예: 완벽하게 정지된 방) "거울 이미지" 문제를 바라보았습니다. 이 실제 문제와 단순한 거울 이미지를 비교함으로써, 그들은 모든 힘든 과정을 거치지 않고도 복잡한 힘을 계산할 수 있었습니다. 이는 공기가 풍선을 어떻게 밀고 당기는지에 대한 숨겨진 세부 사항을 볼 수 있게 해주는 지름길입니다.

2. "떨림"이 중요하다 (속도 분산)

많은 기존 모델에서는 모든 풍선이 정확히 같은 속도로 움직인다고 가정했습니다. 하지만 현실에서는 어떤 풍선은 더 빨리 떠내려가고, 어떤 것은 더 느리게 가며, 어떤 것은 위아래로 흔들릴 수 있습니다. 이러한 "떨림" 또는 **속도 분산(velocity variance)**은 사람들이 걷는 군중과 같습니다. 만약 모든 사람이 정확히 같은 속도로 걷는다면 질서 정연하겠지만, 누군가는 질주하고 누군가는 산책한다면 그 군중은 다른 종류의 압력을 만들어냅니다.

저자들은 이 "떨림"이 추가적인 힘을 만들어낸다는 것을 발견했습니다.

• 항력 (Drag Force): 공기는 단순히 풍선의 평균 속도에 따라 밀어내는 것이 아닙니다. 또한 풍선들이 평균 속도를 중심으로 얼마나 떨리고 있는지에 따라서도 밀어냅니다.
• 응력 (Stress, "쥐어짜는 힘"): 전체 풍선 집단을 관찰할 때, 그들의 떨림은 주변 공기에 추가적인 "쥐어짜는 힘"이나 압력을 만들어냅니다. 이는 마치 사람들이 초조하게 발을 구르는 군중과 같습니다. 비록 그들이 달리고 있지는 않더라도, 그들의 안절부절못하는 움직임은 방 안에 압박감을 만들어냅니다.

3. "속도의 제곱" 효과

가장 중요한 발견 중 하나는 이러한 힘이 더 빠르게 움직일 때 어떻게 행동하는가 하는 점입니다.

• 매우 느린, 시럽 같은 세상에서는 힘이 속도에 정비례합니다 (속도가 두 배가 되면 밀어내는 힘도 두 배가 됩니다).
• 이 새로운, 약간 더 빠른 세상에서는 힘이 속도의 제곱에 의존하기 시작합니다.

쇼핑 카트를 미는 장면을 상상해 보세요. 살살 밀면 쉽습니다. 하지만 두 배로 세게 밀면, 단순히 두 배로 힘든 것이 아니라 공기 저항과 바닥과의 상호작용 때문에 훨씬 더 힘들게 느껴집니다. 저자들은 이러한 떠다니는 액적들의 경우, "밀어내는 힘"이 속도 자체보다 훨씬 더 빠르게 증가하며, 또한 액적들이 얼마나 떨리고 있는지에 크게 좌우된다는 것을 보여주었습니다.

4. 이것이 왜 유체의 "레시피"를 바꾸는가

이 논문은 만약 여러분이 공기와 떠다니는 액적의 혼합물(예: 기포탑 또는 부유 탱크)이 어떻게 행동하는지 설명하고 싶다면, 예전의 단순한 레시피를 사용할 수 없다고 결론짓습니다.

• 옛날 레시 recipe: "액적이 얼마나 있는지를 바탕으로 점성(끈적임)을 더한다."
• 새로운 레시피: "점성을 더하되, 액적이 공기에 비해 얼마나 빠르게 움직이는지에 의존하는 항을 더하고, 또 다른 항으로서 그들이 얼마나 떨리고 있는지에 의존하는 항을 더한다."

이는 이 혼합물이 단순한 끈적한 액체(예: 꿀)라기보다는, 움직임의 속도와 혼란 정도에 따라 행동이 변하는 스마트 재료처럼 행동한다는 것을 의미합니다.

요약

요약하자면, 핀치와 피에르송은 영리한 수학적 거울을 사용하여, 떠다니는 액적이 약간의 속도를 가진 유체 속을 이동할 때 다음과 같은 현상이 일어남을 보여주었습니다.

1. 관성이 중요하다: 유체의 "추진력(oomph)"이 항력을 변화시킵니다.
2. 떨림이 중요하다: 액적들 사이의 무작위적인 속도 차이가 추가적인 힘과 압력을 만들어냅니다.
3. 비선형적 행동: 힘은 단순히 속도에 따라 증가하는 것이 아니라, 속도의 제곱 및 떨림의 제곱에 따라 증가합니다.

이는 엔지니어들에게 이 혼합물(예: 산업용 분리 탱크)이 어떻게 흐르는지 예측하려면, 단순히 평균 속도뿐만 아니라 액적들의 "안절부절못하는 움직임(fidgeting)"까지 고려해야 한다는 점을 이해하도록 도와줍니다.

기술 요약

기술 요약: 희박 부유 액적 현탁액의 관성 효과가 계면 항력 및 유변학에 미치는 영향

문제 정의
본 연구는 저(低)하지만 유한한 레이놀즈 수($Re = aU\rho_f/\mu_f$) 영역과 희박한 체적 분율($\phi \ll 1$)을 가진 부유 액적으로 구성된 분산 이상 유동의 모델링 문제를 다룬다. 스토크스 유동(점성 지배) 및 포텐셜 유동(비점성) 한계에 대해서는 광범위한 문헌이 존재하지만, 부유 액적에 대해 저(低)하지만 유한한 $Re$ 영역에서의 고차 유체역학적 힘 모멘트에 대한 폐쇄 모델(closure models)은 부족한 실정이다. 구체적으로, 본 연구는 유체 관성이 계면 항력과 유체역학적 힘의 1차 및 2차 모멘트를 어떻게 수정하는지 정량화하고자 한다. 이러한 모멘트들은 두 유체 모델(two-fluid models)에서 평균화된 운동량 방정식을 닫기 위해 필수적이며, 연속상의 유효 응력과 상간의 운동량 교환에 기여한다. 저자들은 상대 속도가 관성을 유도하면서도 표면 장력이 구형을 유지하기에 충분한($Ca \ll 1$) 이동하는 구형 액적의 특정 구성을 중점적으로 다룬다.

방법론
저자들은 폐쇄 관계를 도출하기 위해 상호 역 정리(reciprocal theorem)와 결합된 엄격한 통계적 평균화 절차를 채택한다.

1. 평균 방정식: 연구는 분산상과 연속상을 위한 앙상블 평균 질량 및 운동량 방정식으로부터 시작된다. 폐쇄 문제는 앙상블 평균된 항들(예: 계면 힘 $F$ 및 유효 응력 $\Sigma^{\text{eff}}$)을 거시적 변수(수밀도, 체적 분율, 평균 속도)로 표현하는 문제를 포함한다.
2. 조건부 평균화: 폐쇄 문제를 해결하기 위해, 저자들은 조건부 평균화된 장(field)을 활용한다. 특정 질량 중심 속도 $w$를 가진 "테스트 액적"이 배경 유동 속에 잠겨 있는 상황을 고려한다. 이 액적 주변의 교란 장에 대한 지배 방정식은 나비에-스토크스 방정식을 조건부 평균화함으로써 유도된다.
3. 일반 상호 역 정리: 본 연구의 핵심적인 방법론적 기여는 액적에 적용 가능한 일반적인 상호 역 정리를 유도하는 것이다. 이 정리는 임의의 유동 내 액적에 작용하는 힘 모멘트를 알려진 보조 해(스토크스 유동 해 및 점원(point source)이나 점력(point force)과 같은 특이점 해)와 연결한다. 기존 연구들이 오직 힘이나 스트레스렛(stresslet)에만 집중했던 것과 달리, 이 공식화는 임의 차수의 모멘트를 체계적으로 계산한다.
4. 섭동 분석: 저자들은 이 정리를 균일한 정상 상태 유동 내의 액적에 적용한다. 이들은 항력, 1차 모멘트(스트레스렛), 그리고 힘의 2차 모멘트에 대한 $O(Re)$관성 보정치를 계산한다. 이는 유한한$Re$에서 교란 장의 비균일성을 처리하기 위해 내부(스토크스) 해와 외부(오신(Oseen)) 해를 매칭하는 과정을 포함한다.
5. 앙상블 평균: 마지막으로, 유도된 힘 모멘트들을 액적 속도 분포에 대해 앙상블 평균한다. 이 단계는 분산상의 속도 분산($\langle \delta_p u'_\alpha u'_\alpha \rangle$)에 의존하는 항들을 거시적 방정식에 명시적으로 도입한다.

주요 결과
본 연구는 유체역학적 힘 및 모멘트에 대한 관성 보정에 관한 명시적인 해석적 표현식을 제공한다:

• 항력(Drag Force): 항력에 대한 $O(Re)$보정치는$O(\rho_f \phi U^2)$에 비례한다. 이 결과는 이동하는 액적에 대한 항력을 다룬 고전적인 Taylor 및 Acrivos (1964)의 해를 회복하지만, 이는 상호 역 정리와 앙상블 평균의 맥락 내에서 도출되었다.
• 1차 모멘트(Stresslet): 힘의 1차 모멘트에 대한 관성 보정은 $O(\rho_f \phi U^2)$의 척도를 갖는다. 주목할 점은, 액적이 구형을 유지하는 한, 액적 밀도 비($\zeta = \rho_p/\rho_f$)가 1차 및 2차 모멘트의 선도 차수(leading-order) 관성 보정에는 나타나지 않는다는 것이다. 이는 다른 맥락에서 밀도가 변형을 통해 관여하는 항력의 경우와 대조된다.
• 2차 모멘트: 힘의 2차 모멘트에 대한 $O(Re)$보정은$O(a \rho_f \phi U^2)$의 척도를 갖는다. 저자들은 트레이스(trace)와 무영(traceless) 부분을 모두 포함한 완전한 텐서 표현식을 유도한다.
• 속도 분산 기여: 액적 속도 분포에 대해 힘 모멘트를 앙상블 평균하는 과정이 분산상의 속도 분산에 비례하는 추가적인 기여를 도입한다는 것이 핵심적인 발견이다. 구체적으로, 평균화된 항력과 연속상의 유효 응력은 상대 속도의 제곱에 비례하고 속도 분산에 선형적으로 의존한다.
• 유효 응력: 결과적으로 도출된 연속상의 유효 응력 텐서는 $\rho_f \phi U_r U_r$ 및 $\rho_f \phi \langle u' u' \rangle$에 비례하는 항들을 포함한다. 이러한 항들은 현탁액이 비뉴턴(non-Newtonian) 거동을 보임을 나타내며, 유효 응력이 상대 속도의 제곱과 체적 분율 및 상대 속도의 공간적 구배에 의존함을 보여준다.

의의 및 주장
저자들은 본 연구를 부유 액적과 난류가 공존하는 부유법(flotation)이나 액-액 추출과 같은 산업 공정에서 보다 정확한 두 유체 모델링을 위한 필수적인 단계로 규정한다.

• 프레임워크의 신규성: 본 논문은 유체역학적 1차 및 2차 힘 모멘트의 류론적(rheological) 기여를 명시적으로 포함하여, 평균화된 두 유체 흐름 방정식의 틀 안에서 저(低)하지만 유한한 $Re$를 갖는 부유 액적을 조사한 첫 번째 연구라고 주장한다.
• 방법론적 엄밀성: 일반적인 상호 역 관계를 유도함으로써, 저자들은 고체 입자나 포텐셜 유동에 사용되었던 기존 방법론을 확장하여 임의 차수의 힘 모멘트를 계산할 수 있는 체계적인 도구를 제공한다.
• 폐쇄에 대한 함의: 본 연구는 고립된 입자에 대해 도출된 표준 항력 법칙이 평균화된 프레임워크에 적용될 때 반드시 수정되어야 함을 강조한다. 입자 속도 분포에 대한 평균화 과정은 종종 무시되거나 경험적으로 모델링되는 고차 모멘트(속도 분산 항)를 자연스럽게 생성한다. 저자들은 이러한 항들이 평균화 절차의 물리적으로 근거가 있는 결과임을 보여준다.
• 유변학적 거동: 본 연구는 저(低) $Re$에서 부유 액적의 희박 현탁액이 속도 구배와 분산에 의존하는 2차 구배 유체(또는 특정 한계에서의 Reiner-Rivlin 유체)처럼 행동함을 입증한다. 그러나 저자들은 1차 모멘트의 관성 기여는 중요하지만, 2차 모멘트의 관성 기여는 척도 분리 가정($a/L \ll 1$)으로 인해 실제 모델에서는 무시될 수 있다고 주의를 준다.
• 한계점: 저자들은 속도 분산 항 자체에 대한 명시적인 모델을 제안하지 않았기에 방정식계가 완전히 닫히지는 않았다고 겸허히 언급하였으나, 이를 위해 기존의 운동론(kinetic theory) 및 난류 모델링 접근 방식을 참조하였다. 또한, 본 결과는 구형 액적과 깨끗한 계면을 가정하며, 마랑고니 응력(Marangoni stresses) 및 상당한 액적 변형은 제외한다.