Mixing of miscible liquids: Dimensionless scaling for intermediate-to-large density differences in a stirred tank

본 연구는 혼합된 액체가 포함된 교반 탱크의 수치 시뮬레이션을 활용하여, 혼합 시간이 리처드슨 수와 양의 상관관계를 보이지만, 파워, 프루드 및 리처드슨 수를 기반으로 유도된 지수 스케일링 법칙은 중간에서 큰 밀도 차이 범위의 모든 데이터를 단일 마스터 곡선으로 수렴시킨다는 것을 입증한다.

핵심

거대한 두 개의 서로 다른 액체를 섞는 상황을 상상해 보세요. 바닥에는 무겁고 끈적한 시럽이 있고, 위에는 가볍고 묽은 주스가 있습니다. 여기에 거대한 회전식 paddles(임펠러) 를 넣어 두 액체를 섞으려 합니다.

실제 세계에서는 의약품부터 하수 처리에 이르기까지 다양한 공장에서 이러한 작업이 흔히 이루어집니다. 하지만 함정이 하나 있습니다. 액체들의 밀도, 즉 무게가 다르기 때문에 무거운 액체는 바닥에 가라앉으려 하고 가벼운 액체는 위로 떠오르려 합니다. 이로 인해 액체를 섞으려 하는 회전식 paddles 와 액체를 분리시키려 하는 중력 사이에 '싸움'이 발생합니다.

이 논문은 마치 탐정 이야기와 같습니다. 과학자들이 강력한 컴퓨터 시뮬레이션을 사용하여 두 액체가 완벽하게 섞이는 데 정확히 얼마나 시간이 걸리는지, 그리고 매번 물리적인 탱크를 제작하고 값비싼 실험을 수행하지 않고도 그 시간을 어떻게 예측할 수 있는지 규명했습니다.

설정: 디지털 실험실

연구진은 표준 산업용 혼합 탱크의 가상 버전을 구축했습니다.

• 탱크: 벽과 액체가 게으른 강처럼 원형으로만 돌지 않도록 막아주는 네 개의 수직 판 (배플) 이 있는 큰 원통입니다.
• paddle: 중앙에 있는 회전하는 날개입니다.
• 액체: 무거운 액체와 가벼운 액체를 50 대 50 으로 섞은 것을 시뮬레이션했습니다. 실제 화학 물질을 사용한 것이 아니라, 점도 (끈기) 가 동일한 '무거운' 유체와 '가벼운' 유체로만 취급했습니다.
• 방법: 표준 수학 방정식을 대신하여 **격자 볼츠만 방법 (Lattice Boltzmann Method)**이라는 영리한 트릭을 사용했습니다. 이를 액체를 연속적인 덩어리가 아니라, 수조 개의 보이지 않는 당구공들이 튀어 오르고 충돌하는 것으로 시뮬레이션하는 것이라고 생각하세요. 이를 통해 난류 (혼란스러운 소용돌이) 가 어떻게 행동하는지 정확히 파악할 수 있었습니다.

핵심 질문: 얼마나 빨리 섞을 수 있을까?

주요 목표는 혼합 시간을 예측하는 '마법 공식'을 찾는 것이었습니다.

• 변수: 두 가지 주요 요소를 변경했습니다.
  1. paddle 의 회전 속도 (레이놀즈 수): 회전 속도가 빠를수록 일반적으로 난류가 증가하고 혼합이 빨라집니다.
  2. 무게 차이 (리처드슨 수): 액체들의 무게가 거의 같으면 쉽게 섞입니다. 하지만 하나가 훨씬 무거우면 중력이 혼합을 방해하여 깨기 어려운 층을 형성합니다.

발견: '중력 대 회전'의 싸움

연구진은 몇 가지 흥미로운 패턴을 발견했습니다.

1. 중력이 영향을 미치지 않을 때 (동일한 무게):
   두 액체의 무게가 정확히 같다면, 혼합 시간은 놀라울 정도로 일정합니다. paddle 을 얼마나 빠르게 회전시키든 (특정 범위 내에서), '무차원 혼합 시간' (즉, 'paddle 이 몇 바퀴 회전해야 하는가'라는 fancy 한 표현) 은 약 20 바퀴로 일정하게 유지됩니다. 이는 마치 자연의 법칙과 같습니다. 물이 충분히 교반되면, paddle 회전 수로 볼 때 더 빠르게 회전한다고 해서 혼합이 더 빨라지지 않는다는 것입니다.

2. 중력이 반격할 때 (서로 다른 무게):
   액체들의 무게가 다르면 무거운 액체는 바닥에 머물고 싶어 합니다. 무게 차이가 클수록 혼합하기 어려워집니다.

   • 경향: 무게 차이가 클수록 혼합에 더 많은 시간이 걸립니다.
   • 놀라운 반전: '무게 차이'를 일정하게 유지한 채 paddle 만 더 빠르게 회전시킨다면, 혼합 시간이 항상 줄어드는 것은 아닙니다. 때로는 더 빠르게 회전할수록 특정 혼합 지점에 도달하는 데 오히려 더 오래 걸립니다.
   • 이유는 무엇일까요? 무거운 액체를 두꺼운 담요라고 상상해 보세요. paddle 을 너무 빠르게 회전시키면 많은 에너지가 발생하지만, 무거운 액체는 가벼운 액체가 침투할 수 없는 안정적인 '덮개'를 형성합니다. 에너지는 상층부를 소용돌이치게 하는 데 낭비되고 하층부는 고립된 채로 남습니다. 마치 무거운 야채들이 바닥에 단단한 덩어리로 가라앉은 수프 냄비를 저으려 할 때, 숟가락을 더 빠르게 저어도 야채 덩어리를 부수는 대신 위쪽 국물만 튀기는 것과 같습니다.

해결책: 새로운 '마스터 곡선'

이 팀의 가장 큰 업적은 이러한 모든 요소를 결합한 단일하고 간단한 공식을 만드는 것이었습니다. 그들은 혼합 시간을 세 가지 특정 수치 (동력, 프루드 수, 리처드슨 수) 의 관점에서 바라보면, 모든 복잡한 데이터 포인트가 하나의 매끄러운 지수 곡선 위에 모인다는 것을 발견했습니다.

이렇게 생각하세요. 이전에는 엔지니어들이 새로운 액체가 어떻게 섞일지 예측하기 위해 추측하거나 수백 번의 실험을 수행해야 했습니다. 이제 그들은 '레시피'를 갖게 되었습니다. 무게 차이와 회전 속도를 알려주기만 하면 이 공식이 높은 정확도로 혼합 시간을 예측해 줍니다.

결론

이 논문은 이러한 특정 산업용 탱크에 대해 다음과 같이 결론 내립니다.

• 난류가 핵심입니다: 액체가 완전히 교반되면 혼합 행동은 예측 가능합니다.
• 중력이 지배적입니다: 액체들의 밀도가 다르다면 중력이 혼합을 저항하는 '층화 (stratification)'를 생성합니다.
• 빠른 것이 항상 좋은 것은 아닙니다: 밀도 차이가 큰 시스템에서는 단순히 모터 속도를 높인다고 해서 혼합이 빨라지는 것이 보장되지 않습니다. 때로는 오히려 더 안정적인 분리를 초래할 뿐입니다.

저자들은 이 새로운 공식을 제공하여 엔지니어들이 값비싼 프로토타입을 먼저 제작하지 않고도 더 나은 혼합 공정을 설계할 수 있도록 돕습니다. 그들은 향후 이 공식을 다른 탱크 모양과 paddle 유형에 대해 테스트할 계획이지만, 현재로서는 시뮬레이션한 표준 탱크에는 완벽하게 작동합니다.

기술 요약

기술 요약: 교반 탱크 내 중간~대규모 밀도 차이를 위한 무차원 스케일링

문제 제기
혼화성 액체의 혼합은 제약 생산에서 폐수 처리에 이르기까지 산업 공정에서 중요한 단위 조작입니다. 단일 상 균질 시스템에서의 난류 혼합은 잘 규명되어 있지만, 유체 간 밀도 차이가 존재하면 부력에 의한 성층화가 발생합니다. 이 현상은 중력력이 대류 수송을 상쇄하는 수직 층상 영역을 형성하여 수직 혼합을 억제하고 균일성을 달성하는 데 필요한 시간을 증가시킬 수 있습니다. 단일 상 추적자 혼합 및 특정 성층화 사례에 대한 경험적 상관관계는 존재하지만, 다양한 유동 체제와 기하학적 구조에 걸쳐 중력, 관성, 점성 및 부력 효과를 일관되게 고려하는 일반적으로 유효한 프레임워크는 여전히 부재합니다. 본 연구는 중간~대규모 밀도 차이를 가진 두 가지 혼화성 액체가 포함된 격벽이 있는 교반 탱크에서의 혼합 시간을 예측하는 과제를 다룹니다.

방법론
저자들은 격벽이 있는 교반 탱크 내 난류 혼합을 시뮬레이션하기 위해 격자 볼츠만 방법 (LBM) 을 사용한 전산 유체 역학 (CFD) 을 활용했습니다. 시스템은 동일한 동점도를 가지지만 밀도 ($\rho_l$ 및 $\rho_h$) 가 다른 두 가지 혼화성 액체를 50/50 부피 비율로 구성되었으며, 여기서 가벼운 유체는 처음에 무거운 유체 위에 성층화되어 있었습니다.

• 수치적 접근법: 시뮬레이션은 D3Q19 격자 체계를 사용하여 비압축성 나비에 - 스토크스 방정식을 풀었습니다. 높은 레이놀즈 수 ($Re > 10^5$) 를 처리하기 위해 서브 그리드 난류를 모델링하기 위해 스마고린스키 - 릴리 대형 와류 시뮬레이션 (S-L LES) 모델을 구현했습니다.
• 스칼라 수송: 가벼운 액체의 농도를 추적하기 위해 질량 플럭스 기반 스칼라 수송 모델을 채택했습니다. 이 모델은 농도 의존성 중력 체력을 통해 스칼라 장을 운동량 방정식과 결합하여 국부적 밀도 변화로 인한 부력 효과를 시뮬레이션할 수 있게 했습니다.
• 시뮬레이션 공간: 임펠러 회전 속도 ($N$) 와 가벼운 상의 밀도 ($\rho_l$) 를 변화시켜 체계적인 매개변수 연구를 수행했습니다. 이로 인해 레이놀즈 수 ($Re$) 가 50,000 에서 125,000 까지, 리처드슨 수 ($Ri$) 가 0 에서 14 까지 광범위한 무차원 수 범위가 도출되었습니다.
• 검증: 수치적 정확도는 격자 수렴 지수 (GCI) 분석을 통해 검증되었으며, 이는 생산 격자에서 약 2.39% 의 수치적 불확실성을 추정했습니다. 또한, 시뮬레이션 결과는 문헌의 실험적 추적자 데이터와 비교하여 검증되었으며, 프로브 응답 및 혼합 시간 정의에서 좋은 일치를 보였습니다.

주요 결과
본 연구는 90% 균일성에 도달하는 데 필요한 시간으로 정의된 혼합 시간 ($t_m$) 을 정량화하고, 이를 무차원 형태 ($t^*_m = N \cdot t_m$) 로 표현했습니다.

1. 균질 기준선 ($Ri = 0$): 밀도 차이가 없는 경우, 무차원 혼합 시간은 조사된 전체 레이놀즈 수 범위에서 약 20 의 일정 값으로 수렴했습니다. 이는 유동이 완전히 난류 상태이며 관성 혼합이 지배적임을 확인한 것으로, 격벽이 있는 탱크에 대한 기존 문헌과 일치합니다.
2. 성층화의 영향 ($Ri > 0$): 리처드슨 수가 증가함에 따라 밀도 차이가 커짐을 나타내며, 무차원 혼합 시간은 체계적으로 증가했습니다. 데이터는 모든 레이놀즈 수에 걸쳐 명확한 지수 추세를 따랐습니다.
3. 제안된 스케일링 법칙: 저자들은 모든 시뮬레이션 데이터를 단일 지수 스케일링 관계로 수렴시키는 마스터 곡선을 도출했습니다:
   $$t^*_m = 20 \cdot e^{0.58 \cdot N_p \cdot Fr \cdot Ri}$$
   여기서 $N_p$는 동력수, $Fr$은 프루드 수,$Ri$는 리처드슨 수입니다. 계수 20 은 관성 기준선을 나타내며, 지수 항은 부력의 안정화 효과를 포착합니다.
4. 비단조적 거동: 일정 $Ri$에서 차원 있는 혼합 시간 ($t_m$) 을 분석할 때, 연구는 높은 리처드슨 수 ($Ri \geq 8$) 에서$Re$에 대한 비단조적 의존성을 관찰했습니다.$Re$(및 따라서 임펠러 속도) 를 증가시키는 것은 초기에 난류를 강화하여 혼합 시간을 단축시키지만, 추가적인 증가는 시뮬레이션 설정에서 매개변수의 결합으로 인해 더 강한 밀도 성층화를 초래하여 수직 수송을 억제합니다. 이러한 상호 작용은 중간 레이놀즈 수에서 최소 혼합 시간을 초래합니다.
5. 유동 구조: 속도장 분석은 높은 $Ri$에서 가벼운 상이 임펠러의 동력 입력으로부터 "봉쇄"될 수 있으며, 수직 구조가 주로 더 밀도가 높은 액체에서 나타날 수 있음을 보여주었습니다. 파워 스펙트럼 밀도 (PSD) 분석은 부력이 유동의 공간적 조직에 영향을 미치지만, 임펠러 구동 역학이 여전히 에너지 캐스케이드를 지배함을 나타냈습니다.

의의 및 주장
본 논문은 완전히 난류 영역 내에서 초기에 밀도 성층화된 시스템의 난류 혼합 시간에 대한 간결하고 예측 가능한 상관관계를 제공한다고 주장합니다. 동력수, 프루드 수, 리처드슨 수의 결합된 영향에 기반한 도출된 스케일링 법칙은 엔지니어 및 응용 전문가를 위한 간단한 분석적 표현을 제공합니다.

저자들은 완전히 난류 영역의 성층화 시스템에 대한 비교 가능한 상관관계가 이전에 보고된 바가 없다고 지적합니다. 이 연구의 의의는 값비싼 실제 테스트에 대한 의존도를 줄여 혼합 공정 설계를 간소화할 수 있는 능력에 있습니다. 제안된 모델은 부력 효과가 중요한 시스템에서 혼합 시간을 추정할 수 있게 하며, 강화된 난류와 강화된 성층화 사이의 복잡한 경쟁을 포착합니다. 저자는 겸손하게 향후 연구가 제안된 스케일링의 견고성을 평가하기 위해 각도 교반기 등을 포함한 다양한 탱크 기하학과 임펠러 유형에서 이 상관관계를 테스트하는 것을 포함할 것이라고 제안합니다.