Hydrodynamic Behavior of Non-spherical Particles in Confined Vertical Flows: A Resolved CFD-DEM Study

본 연구는 해상된 CFD-DEM 시뮬레이션을 활용하여 비구형 다금속 결핵이 형상 유도 와류 비대칭성으로 인해 체적 등가 구체에 비해 현저히 향상된 항력과 감소된 종단 속도를 경험함을 입증하는 동시에, 입자 크기와 구속 조건이 수직 수리 수송 과정 중의 뚜렷한 항력 변화 거동을 어떻게 지배하는지를 밝혀낸다.

핵심

당신은 심해저의 기묘하고 울퉁불퉁한 돌덩이들(다금속 결절이라고 불리는)을 거대한 수직 빨대를 이용해 배 위로 끌어올리려 한다고 상상해 보십시오. 이것이 바로 심해 채굴이 작동하는 방식입니다. 엔지니어들의 핵심적인 질문은 이것입니다: 이 기묘한 모양의 돌들이 물이 빨대 속으로 치솟아 오를 때 어떻게 움직이는가?

이 시스템을 설계하기 위해 사용되는 대부분의 컴퓨터 모델은 모든 돌을 마치 완벽하고 매끄러운 구슬처럼 취급합니다. 하지만 실제로는 이 돌들이 울퉁불퉁하고 불규칙하며, 구슬과는 전혀 다르게 생겼습니다. 이 논문은 다음과 같은 질문을 던집니다: 울퉁불퉁한 돌을 매끄러운 구슬처럼 취급하는 것이 실제로 효과가 있을까요, 아니면 틀린 답을 주는 것일까요?

이를 알아내기 위해 연구진은 지름길을 쓰지 않는 초정밀 컴퓨터 시뮬레이션(마치 고성능 비디오 게임 물리 엔진과 같은)을 구축했습니다. 대신 물이 돌을 밀어내는 힘을 단순히 추측하는 대신, 돌의 모든 돌출부와 틈새에 가해지는 물의 압력을 일일이 계산했습니다.

연구 결과는 다음과 같이 쉽게 설명할 수 있습니다.

1. "울퉁불퉁한 돌" vs "매끄러운 구슬"

연구진은 정지된 물속에 이 울퉁불퉁한 돌들을 떨어뜨려 얼마나 빨리 가라앉는지 관찰했는데, 울퉁불퉁한 돌들은 무게와 크기가 정확히 같은 매끄러운 구슬보다 약 28% 더 느리게 가라앉았습니다.

• 비유: 물속에서 수영하는 것을 상상해 보세요. 만약 당신이 매끄럽고 유선형인 돌고래라면 쉽게 미끄러져 나갈 것입니다. 하지만 당신이 울퉁불퉁하고 거친 나무토막이라면, 내려가는 길에 물을 더 많이 걸리게 될 것입니다.
• 이유: 울퉁불퉁한 돌들은 더 큰 "전면 면적"(물을 더 많이 받음)을 가지고 있으며, 뒤쪽에 무질서하고 비대칭적인 항적(거품 궤적과 같은 혼란스러운 흔적)을 남깁니다. 이 추가적인 항력이 속도를 크게 늦춥니다.
• 주의점: 비록 더 느리게 가라앉지만, 돌을 밀어 올리는 총 부력(위로 떠오르는 힘)은 구슬과 동일합니다. 단지 그 힘과 균형을 맞추기 위해 더 느리게 움직여야 할 뿐입니다.

2. 파이프 속의 "교통 체증"

다음으로, 연구진은 돌을 운반하기 위해 파이프 속으로 물이 솟구쳐 오르는 상황을 시뮬레이션했습니다. 그들은 "작은" 돌과 "큰" 돌 두 가지 크기를 살펴보았습니다.

• 매끄러운 구슬: 물의 속도가 빨라질수록 매끄러운 구슬들은 예측 가능한 방식으로 행동했습니다. 낮은 속도에서는 흔들리며 가라앉았지만, 높은 속도에서는 고속도로에 합류하는 자동차들처럼 직선으로 빠르게 올라갔습니다.
• 울퉁불퉁한 돌: 이들은 훨씬 더 혼란스러웠습니다.
  • 낮은 속도에서: 작은 울퉁불퉁한 돌들은 파이프 위로 올라가지 못했습니다! 그들은 바닥 근처에서 맴돌며 제자리에서 회전하고 흔들렸으며, 중력을 극복하지 못했습니다. 반면 매끄러운 구슬들은 위로 이동할 수 있었습니다.
  • 높은 속도에서: 물이 충분히 빨라져서 돌을 운반할 수 있는 속도가 되었을 때도, 울퉁불퉁한 돌들은 꼭대기에 도달하는 데 더 오랜 시간이 걸렸고 훨씬 더 불규칙하게 회전하며 이동했습니다. 이들은 마치 에스컬레이터를 타고 올라가면서 제자리에서 빙글빙글 도는 사람들의 무리 같았고, 매끄러운 구슬들은 그냥 곧게 달려 올라가는 모습이었습니다.

3. "팽이" 효과

가장 큰 차이점은 돌이 회전하는 방식이었습니다.

• 매끄로운 구슬: 이들은 주로 위로 올라가기만 했습니다. 별로 회전하지 않았습니다.
• 울퉁불퉁한 돌: 돌의 표면이 울퉁불퉁하기 때문에, 부딪히는 물이 돌을 격렬하게 회전시켰습니다. 이 회전(회전 운동)은 위아래로 움직이는 운동과 밀접하게 연결되어 있었습니다.
• 비유: 매끄러운 구슬을 총에서 발사된 탄환이라고 생각하면 직선으로 나갑니다. 반면 울퉁불퉁한 돌은 바람 터널 속에 던져진 부메랑이나 팽이와 같습니다. 모양 때문에 계속 뒤틀리고 방향을 바꾸며 회전합니다. 이 회전은 물과의 추가적인 "마찰"을 만들어내어 운반을 더 어렵게 만듭니다.

4. "힘의 변동" (덜컹거리는 승차감)

연구진은 물이 돌에 가하는 "밀어내는 힘"(항력)을 측정했습니다.

• 작은 돌: 매끄럽든 울퉁불퉁하든 밀어내는 힘은 비교적 일정했습니다.
• 큰 돌: 여기서 놀라운 일이 벌어졌습니다.
  • 큰 매끄러운 구슬: 물이 스쳐 지나갈 때 힘이 약간씩 변하며 예측 가능한 패턴의 "요철"을 만들었습니다.
  • 큰 울퉁불퉁한 돌: 힘의 변화가 매우 예측 불가능했습니다. 돌이 회전하고 물에 대한 형태가 계속 변했기 때문에, 힘이 갑자기 급증했습니다. 이는 마치 매끄러운 도로(매끄러운 구슬)를 달리는 자동차와, 자동차가 기울어진 정도에 따라 매초마다 요철이 바뀌는 도로(울퉁불퉁한 돌)를 달리는 자동차의 차이와 같습니다.

결론

연구는 매끄러운 구슬 모델을 사용하여 이 돌들이 어떻게 행동할지에 대한 대략적인 아이디어를 얻을 수는 있지만, 세부 사항을 놓치고 있다고 결론짓습니다.

• 만약 매끄러운 구슬 모델을 사용한다면, 돌이 실제보다 더 빠르고 쉽게 파이프를 올라갈 것이라고 생각할 수 있습니다.
• 울퉁불퉁한 돌은 움직이기 위해 더 빠른 물의 속도가 필요하며, 일단 움직이기 시작하면 회전하고 흔들리기 때문에 더 불안정하고 제어하기 어렵습니다.

요약하자면: 자연은 무질서합니다. 실제로 작동하는 기계를 설계하고 싶다면, 단순히 들쭉날쭉한 돌을 완벽한 구형인 것처럼 가정해서는 안 됩니다. "울퉁불퉁함"은 단순한 모델이 무시하는 많은 양의 추가적인 항력과 혼돈을 불러옵니다.

기술 요약

기술 요약: 구속된 수직 유동 내 비구형 입자의 유체역학적 거동

문제 정의
심해 평원의 다금속 결절(PMN) 회수는 수직 라이저를 통해 거칠고 불규칙하며 비구형인 입자를 운송하는 유압 시스템에 의존한다. 이러한 시스템의 효율은 강한 중력 작용과 기하학적 구속 하에서의 입자 부상, 현탁 안정성 및 운송 거동을 예측하는 데 달려 있다. 현재의 설계 관행은 구형 입자에 맞춰 보정된 경험적 항력 상관관계나 구형도(sphericity)와 같은 스칼라 기술자로 매개변수화된 미해상 CFD-DEM(Computational Fluid Dynamics-Discrete Element Method) 모델에 의존하는 경우가 많다. 이러한 접근 방식은 입자-파이프 직경 비율($d/D \ge 0.2$)이 높은 구속 유동에서, 방향 의존적 항력, 비대칭적인 벽면 반사 후류(wake), 또는 비구형 입자 특유의 회전-병진 운동 결합을 포착하지 못하기 때문에 한계가 있다. 또한, 실제 PMN의 형태를 부피 등가 구체(volume-equivalent spheres)로 대체하는 것이 앙상블 수준의 운송 지표를 보존하는지, 아니면 부상 임계값 및 압력 강하 예측에 편향을 초-래하는지는 여전히 불분명하다.

연구 방법론
본 연구는 심해 채굴 라이저를 대표하는 수직 파이프 내 비구형 PMN의 유체역학을 조사하기 위해 완전히 해상된(fully resolved) CFD-DEM 프레임워크를 채택하였다.

• 수치 프레임워크: 본 접근 방식은 불압축성 나비에-스토크스 방정식(OpenFOAM을 통해 해결)을 이산 입자에 대한 뉴턴 운동 방정식(LIGGGHTS를 통해 해결)과 결합한다.
• 유체-고체 결합: 경험적 항력 법칙 없이 계면에서 유체-고체 상호작용을 직접 해결하기 위해 침투 경계(Immersed Boundary, IB)법을 사용한다. 무슬립(no-slip) 조건은 직접적인 강제력을 통해 적용되며, 유체 응력 텐서를 입자 표면에 대해 적분하여 유체역학적 힘을 계산한다.
• 입자 표현: X선 컴퓨터 단층 촬영(CT) 스캔으로부터 재구성된 복잡한 PMN 기하 구조는 다구체 근사법(중첩되는 구형 하위 입자들의 강체 집합체)을 사용하여 모델링된다. 이를 통해 계산의 용이성을 유지하면서도 불규칙한 표면을 고충실도로 표현할 수 있다.
• 검증: 본 프레임워크는 단일 구체의 침강 실험 데이터(Ten Cate et al., 2002)와 다구체 구체의 드래프팅-키싱-텀블링(DKT) 벤치마크(Glow-inski et al., 2001; Sharma and Patankar, 2005)를 통해 검증되었다. 또한, 구체의 다구체 근사 모델을 매끄러운 구체와 비교 검증하여 기하학적 표현의 정확성을 확인하였다.
• 시뮬레이션 설정: 시뮬레이션은 유효 직경이 각각 20mm와 54mm(부피 등가 직경은 16.4mm 및 44mm)인 PMN과 그에 대응하는 부피 등가 구체를 대상으로, 직경 $D=200$mm인 파이프 내에서 수행되었다. 연구 범위는 입자 레이놀즈 수($Re_p$) 98~2904, 구속비($d/D$) 0.082 및 0.22를 포함한다. 운송 거동은 유속 $1.0u_t$에서 $3.0u_t$($u_t$는 종말 침강 속도) 범위에서 분석되었다.

주요 결과

1. 침강 및 항력 증강:

   • 비구형 PMN은 질량과 부력이 동일함에도 불구하고 부피 등가 구체보다 27~29% 느리게 침강한다.
   • 이러한 종말 속도의 감소는 구체 대비 1.8~2.0배 높은 항력 계수의 증대로 이어진다.
   • 항력 증가는 두 가지 요인에 기인한다: (1) 더 큰 순시 전면적(외접 구체 비율로 인해 부피 등가 면적보다 최대 1.5배 큼)과 (2) 형태 유도 효과(비대칭 후류 구조 및 회전-병진 결합 포함)이다.
   • 결정적으로, 종말 항력은 두 형상 모두에서 잠긴 무게(submerged weight)와 동일하다. PMN의 낮은 속도는 더 낮은 슬립 속도에서 필요한 유체역학적 저항을 생성하기 때문이다.

2. 체류 시간 및 운송 영역:

   • 운송 거동은 속도에 따라 간헐적인 침강 중심 운동에서 안정적인 대류 중심 부상으로 전환되는 양상을 보인다.
   • 저속($u_f = 1.0u_t$)에서 작은 PMN은 입구 근처에서 경계 부상 상태를 유지하며, 순수 상향 운송 없이 진동 운동을 보이는 반면, 구체는 수직 분산을 보인다.
   • 속도가 증가함에 따라 모든 입자의 체류 시간 분포가 좁아진다. 그러나 PMN은 높은 항력과 강한 회전-병진 결합으로 인해 구체에 비해 일관되게 더 긴 평균 체류 시간과 높은 분산을 나타낸다.
   • 고속($u_f = 3.0u_t$)에서는 모든 입자가 공간적으로 균일한 분포를 달성하지만, PMN은 여전히 강화된 회전 운동(각속도 최대 80 rad/s)과 나선형 궤적을 유지한다.

3. 항력 힘 통계 및 변동:

   • 평균 정규화 항력($\hat{f}$)은 모든 조건에서 1에 가깝게 유지되어, 항력과 잠긴 무게 사이의 정상 상태 균형을 확인시켜 준다.
   • 작은 입자 ($d/D = 0.082$): 구체와 PMN 모두 좁은 항력 분포를 보인다. 작은 PMN의 경우, 빠른 재배향(reorientation)이 방향 의존적 투영 면적의 변동을 효과적으로 평균화하여 구체와 정량적으로 유사한 거동을 보인다.
   • 큰 입자 ($d/D = 0.22$):
     • 큰 구체는 후류의 불안정성과 구속 효과로 인해 속도가 증가함에 따라 항력 분포가 단조롭게 넓어지는 양상을 보인다.
     • 큰 PMN은 훨씬 더 넓고 왜곡된 분포와 확장된 고항력 꼬리(high-drag tails)를 나타낸다. 이는 후류 이력 효과와 방향 의존적 힘의 중첩에 기인한다. 작은 입자와 달리, 큰 PMN의 회전 시간 척기는 후류 진화 시간 척기와 유사하여 투영 면적 변화의 평균화를 방지한다.

의의 및 결론
본 연구는 구속된 수직 유동 내에서 입자 형상이 항력, 후류 역학 및 회전-병진 결합에 미치는 영향을 규명한다. 본 연구의 주요 기여는 부피 등가 구형 모델을 비구형 PMN 운송에 적용할 때의 정량적 한계를 설정한 것이다.

• 질적 유사성: 기초적인 운송 물리(침강에서 대류 중심 영역으로의 진행)는 PMN과 구체 간에 질적으로 유사하다. 이는 보정된 항력 법칙을 사용하는 구형 모델을 통해 1차적인 운송 거동을 포착할 수 있음을 시사한다.
• 정량적 차이: 정량적인 차이는 상당하다. PMN은 구체보다 4090% 더 긴 체류 시간을 보이며 항력 계수는 1.82.0배 더 높다. 이러한 차이는 심해 채굴 라이저의 최소 부상 속도 및 압력 강하를 예측하는 데 매우 중요하다.
• 모델링 시사점: 구형 근사 모델이 앙상블 거동의 축소 모델링에는 적합할 수 있으나, 높은 구속비를 가진 큰 비구형 입자와 관련된 특정 간헐성 및 힘의 변동을 포착하는 데는 실패한다. 본 연구에서 제시된 해상된 CFD-DEM 접근 방식은 이러한 모델들을 위한 형상 특화 보정법을 개발하는 데 필수적인 기준을 제공한다.

저자들은 이러한 결과가 희박한 앙상블(dilute ensembles)에 적용되며, 보정 입력값으로서의 역할을 한다는 점을 명시하였다. 침강 저해(hindered settling) 또는 플러그 유동(plug flow)이 포함된 고농도 슬러리로의 확장은 본 연구의 범위를 벗어난다.