A Volume of Fluid Immersed Boundary Method for Industrial Polymer Mixing

본 논문은 부분적으로 채워진 산업용 압출기 내 자유 표면 고분자 혼합을 정확하게 시뮬레이션하기 위해 고점도 대비로 인한 수치적 불안정성을 극복하는 오픈FOAM에 구현된 새로운 블록 결합 체적-유체 침수 경계 (BC-VOF-IB) 솔버를 제시한다.

핵심

거대한 용기 안에서 두껍고 끈적한 꿀을 공기와 섞으려 한다고 상상해 보세요. 이것이 바로 산업용 고분자 혼합에서 일어나는 일입니다. 피렐리와 같은 기업들은 타이어, 의료 기기, 자동차 부품 등을 만들기 위해 녹은 플라스틱에 첨가제를 혼합해야 합니다. 목표는 모든 것이 완벽하게 혼합되어 최종 제품이 강하고 균일하도록 하는 것입니다.

그러나 이 과정을 컴퓨터로 시뮬레이션하는 것은 수학자와 엔지니어들에게 악몽과 같습니다. 그 이유와 이 논문이 어떻게 단순한 비유를 통해 이를 해결하는지 살펴보겠습니다.

문제: "두꺼운 꿀 대 얇은 공기"의 싸움

이러한 기계 안에는 두 가지 매우 다른 유체가 존재합니다.

1. 고분자 용융물: 매우 두껍고 끈적하며 느리게 움직입니다 (차가운 꿀처럼).
2. 공기: 매우 얇고 빠르게 움직입니다.

회전하는 나사가 있는 기계 안에서 이 두 가지가 어떻게 상호작용하는지 시뮬레이션하려 할 때, 표준 컴퓨터 프로그램들은 혼란에 빠집니다. 마치 같은 트랙에서 달리는 달팽이와 레이싱 카의 움직임을 동일한 규칙 세트로 계산하려 하는 것과 같습니다. 컴퓨터는 "달팽이"(두꺼운 플라스틱) 가 너무 빠르게 움직이지 않도록 하기 위해 아주, 아주 작은 단계를 취하려 시도합니다. 이로 인해 시뮬레이션이 극도로 느려져, 실제 몇 초간의 혼합을 완료하는 데 며칠이 걸리기도 합니다.

더욱이, 이러한 기계에는 고정된 용기 안에서 움직이는 복잡한 회전 부품 (나사) 이 있습니다. 전통적으로 이를 시뮬레이션하려면 회전하는 나사를 완벽하게 감싸는 디지털 메쉬 (작은 상자들의 격자) 를 구축해야 합니다. 나사가 회전함에 따라 이 격자는 끊임없이 재형성되어야 하는데, 이는 마라톤을 뛰는 사람이 입고 있는 스웨터를 뜨개질하듯 하려는 것과 같습니다. 이는 messy(지저분하고) 어렵고 오류가 발생하기 쉽습니다.

해결책: 새로운 "스마트 격자"와 "팀 접근법"

이 논문의 저자들은 OpenFOAM 이라는 소프트웨어를 사용하여 이러한 시뮬레이션을 실행하는 새로운 방법을 개발했습니다. 그들은 두 가지 강력한 기법을 결합했습니다.

1. 침수 경계법 (The "Ghost Wall" Trick)
회전하는 나사에 맞춰 격자를 재형성하는 대신, 그들은 격자를 고정되고 단단하게 유지했습니다 (단단한 얼음 덩어리처럼). 그리고 컴퓨터에게 "이 얼음 덩어리 안에 회전하는 나사가 있어"라고 알려주었습니다.

• 비유: 바닥에 고정된 타일 격자가 있는 수영장을 상상해 보세요. 수영객에게 맞춰 타일을 움직이는 대신, 물에게 "수영객을 통과하지 마"라고 말하기만 하면 됩니다. 컴퓨터는 수학적으로 나사 주위에 "유령 벽"을 만들어 유체가 나사를 우회하도록 강제하며, 격자를 다시 구축할 필요가 없습니다. 이로 인해 복잡하고 움직이는 형태를 처리하는 것이 훨씬 쉬워집니다.

2. 유체 부피법 (VOF, The "Tracking Paint" Trick)
두꺼운 플라스틱이 어디에서 끝나고 공기가 시작되는지 파악하기 위해, 그들은 셀을 채우는 "페인트"를 사용합니다.

• 비유: 컴퓨터 격자를 3 차원 체스판이라고 상상해 보세요. 어떤 칸은 100% 플라스틱이고, 어떤 칸은 100% 공기이며, 어떤 칸은 혼합된 상태입니다. 컴퓨터는 각 칸에 들어 있는 "플라스틱 페인트"의 양을 추적하여 액체의 표면을 파악합니다.

3. 블록 결합 방식 (The "Team Huddle")
이것이 가장 중요한 돌파구입니다. 표준 시뮬레이션에서는 컴퓨터가 유체의 X, Y, Z 방향 속도를 하나씩 차례로 계산합니다. 마치 세 사람이 차례로 이야기하는 것과 같습니다. 유체가 매우 두꺼울 때 (고분자와 같이), 이러한 "차례를 기다리는" 접근 방식은 시뮬레이션이 충돌하거나 극도로 느려지게 만듭니다. 두꺼운 유체는 모든 방향을 밀접하게 결합시키기 때문입니다.

저자들은 이를 블록 결합 (Block-Coupled) 방식으로 변경했습니다.

• 비유: 세 사람이 차례로 이야기하는 대신, 모두 한데 모여 동시에 문제를 해결합니다. 모든 방향의 움직임을 하나의 거대하고 상호 연결된 팀으로 취급함으로써, 컴퓨터는 두꺼운 플라스틱과 얇은 공기 사이의 거대한 차이를 걸림 없이 처리할 수 있습니다.

결과: 몇 시간에서 몇 분으로

연구팀은 두 가지 시나리오에서 새로운 방법을 테스트했습니다.

1. 도그본 (Dog-Bone) 모양의 채널: 두꺼운 플라스틱이 좁고 비틀어진 채널로 주입되는 테스트 사례입니다.

   • 기존 방식: 표준 컴퓨터 프로그램은 아주 작은 단계를 강요받았기 때문에 충돌하거나 몇 초를 시뮬레이션하는 데 7 시간이 걸렸습니다.
   • 새로운 방식: 그들의 새로운 "팀 허들" 방법은 같은 작업을 16 분 만에 완료했으며, 플라스틱이 극도로 두꺼워졌을 때도 충돌하지 않았습니다.

2. 실제 산업용 기계: 그들은 플라스틱 펠릿을 만드는 데 사용되는 실제 단일 나사 및 이중 나사 압출기를 시뮬레이션했습니다.

   • 그들은 플라스틱이 기기를 어떻게 채우는지, 압력이 어떻게 상승하는지, 그리고 공기가 어떻게 밀려나는지 성공적으로 보여주었습니다.
   • 그들은 그들의 "유령 벽" 방법이 격자를 재형성하는 기존의 어렵고 복잡한 방법만큼 효과적이지만, 훨씬 빠르고 설정이 쉽다는 것을 입증했습니다.

다음 단계는?

이 논문은 이것이 산업계에 큰 진전이라고 결론지었습니다. 이는 학술적 수학과 실제 공장 요구 사항 사이의 간극을 메워줍니다. 그러나 저자들은 현재 모델이 온도가 일정하다고 가정 (등온) 한다는 점을 지적합니다. 실제로는 플라스틱을 혼합하면 열이 발생하여 플라스틱의 두께가 변합니다. 온도 효과와 더 복잡한 "신축성" 있는 플라스틱 거동을 추가하는 것이 향후 연구의 다음 단계입니다.

요약하자면: 그들은 회전하는 기계 안에서 두꺼운 플라스틱이 공기와 혼합되는 모습을 컴퓨터로 관찰하는 더 빠르고 안정적인 방법을 개발했습니다. 나사가 회전할 때마다 디지털 세계를 다시 구축할 필요 없이, 몇 시간이 걸리던 과정을 몇 분으로 단축했습니다.

기술 요약

기술 요약: 산업용 고분자 혼합을 위한 부피 기반 유체 침수 경계 방법

문제 제기
고분자 혼합 공정의 수치 시뮬레이션은 비뉴턴 유체 역학, 다상 자유 표면 역학, 그리고 이동하는 복잡한 기하학적 구조물이 복잡하게 상호작용하기 때문에 상당한 도전에 직면해 있습니다. 단일 스크류 (SSE) 및 투 스크류 압출기 (TSE) 와 같은 산업용 믹서는 종종 고점도 고분자 용융물과 공기가 공존하는 부분 충전 조건 하에서 작동합니다. 이러한 시나리오를 시뮬레이션하려면 회전하는 고체 경계를 처리하면서도 상 간의 계면을 포착해야 합니다.

이러한 맥락에서 표준 전산 유체 역학 (CFD) 접근법은 두 가지 주요 장애물에 직면합니다:

1. 기하학적 복잡성: 복잡한 기하학적 구조를 가진 회전 스크류에 대한 몸체 적합 (conforming) 메쉬를 생성하는 것은 계산 비용이 많이 들며, 종종 메쉬 왜곡을 초래할 수 있는 동적 메쉬 기법을 필요로 합니다.
2. 수치적 불안정성: 고분자 용융물과 공기 사이의 극심한 점도 차이 (종종 $10^5$에서 $10^8$을 초과) 는 표준 분리형 부피 기반 유체 (VOF) 솔버에서 심각한 수치적 불안정성을 유발합니다. 구체적으로, 점성 확산 항에서 전치된 속도 기울기의 명시적 처리는 엄격한 시간 단계 제약 ($\Delta t \propto h^2$) 을 부과하여 고점도 유동의 시뮬레이션을 계산적으로 불가능하게 만들거나 불안정하게 만듭니다.

방법론
본 연구는 OpenFOAM 유한 체적 라이브러리 내에서 구현된 부피 기반 유체 (VOF) 계면 포착 접근법과 비정합 침수 경계 (IB) 방법을 통합한 새로운 수치 프레임워크를 제안합니다.

• 침수 경계 (IB) 접근법: 이 방법은 회전하는 스크류와 같은 고체 경계가 삼각화된 STL 표면으로 표현되는 고정된 배경 격자를 활용합니다. 가중 최소 제곱 (WLS) 보간 연산자를 사용하여 침수 경계와 교차하는 셀에 경계 조건 (속도에 대한 디리클레, 부피 분율에 대한 뉴만) 을 부과합니다. 이는 몸체 적합 메쉬 생성의 필요성을 제거하고 복잡하며 회전하는 기하학적 구조를 따라 이동하는 접촉선을 처리할 수 있게 합니다.
• VOF 공식화: 계면은 상 부피 분율 ($\alpha$) 에 대한 이송 방정식을 풀어서 추적됩니다. 이 방법은 수치적 확산을 상쇄하기 위한 압축 속도 항을 포함하여 계면의 날카로움과 유계를 유지하기 위해 MULES (Multidimensional Universal Limiter with Explicit Solution) 알고리즘을 사용합니다.
• 블록 결합 (BC) 방식: 높은 점도 차이로 인한 안정성 문제를 해결하기 위해 저자는 운동량 방정식에 대한 블록 결합 방식을 도입합니다. 속도 성분을 순차적으로 처리하는 표준 분리형 솔버와 달리, BC 방식은 전체 점성 확산 항 $\nabla \cdot (\mu(\nabla \mathbf{u} + \nabla^\top \mathbf{u}))$을 암시적으로 처리합니다. 이는 단일 선형 시스템 내에서 세 가지 속도 성분을 결합하는 것을 포함합니다. 접선 기울기 항의 이산화는 탄성 응용 분야에서 차용된 [11] 확장된 스텐실 (면 이웃뿐만 아니라 점 이웃을 포함하는) 에 대한 WLS 보간과 결합된 유한 영역 (FA) 방법을 활용합니다.
• 해법 알고리즘: 솔버는 압력 - 속도 결합을 처리하기 위해 수정된 PIMPLE 알고리즘 (SIMPLE 및 PISO 의 결합) 을 사용합니다. IB 제약 조건은 대수적 시스템에 통합되어 행렬 계수와 소스 항을 수정함으로써 침수 표면에서의 경계 조건을 만족시키면서도 질량 보존을 유지합니다.

주요 기여

1. BC-VOF-IB 솔버 개발: 본 논문은 비정합 IB 방법과 블록 결합 VOF 접근법을 결합한 UCInterIbFoam 솔버를 소개합니다. 이는 OpenFOAM-10 프레임워크 내에서 이상 유동에 대한 이러한 결합 방법의 첫 번째 구현입니다.
2. 블록 결합을 통한 안정성 향상: 이 연구는 점성 확산의 완전한 암시적 처리가 시간 단계 안정성 제약을 크게 완화한다는 것을 보여줍니다. 이로 인해 코란트 수 위반으로 인해 시뮬레이션이 중단되는ことなく 점도 비율이 $10^8$까지인 고점도 고분자 유동을 시뮬레이션할 수 있게 되었습니다.
3. 이동 접촉선 처리: 이 프레임워크는 부분 충전 혼합 장치의 필수 조건인 침수 경계와 교차하는 두 상의 계면이 이동하는 접촉선의 기하학적 및 알고리즘적 어려움을 성공적으로 해결합니다.

수치 결과
제안된 방법은 두 가지 범주의 문제에 대해 검증 및 테스트되었습니다:

• 벤치마크 사례 (도그본 사출 성형):

  • 고점도 파워-로 법칙 유체 주입 사례에 대해 분리형 (interFoam), 블록 결합 정합 (UCInterFoam), 그리고 블록 결합 비정합 (UCInterIbFoam) 솔버 간의 비교가 수행되었습니다.
  • 표준 분리형 솔버는 엄격한 시간 단계 제한 ($\Delta t \approx 10^{-6}$초) 으로 인해 합리적인 시간 내에 안정적인 결과를 생성하지 못했습니다.
  • 블록 결합 솔버는 최대 $10^{-3}$초의 시간 단계로 안정적인 해를 달성하여 계산 시간을 7 시간에서 정합의 경우 16 분, 비정합의 경우 58 분으로 단축했습니다.
  • 비정합 IB 해는 기하학적 제한 부근의 점도 피크에서 격자의 비정합 특성으로 인한 미세한 불일치가 있었으나, 정합 해와 정성적으로 일치하는 것을 보여주었습니다.

• 산업 응용 (SSE 및 TSE):

  • 단순화된 단일 스크류 압출기 (SSE) 와 실제적인 투 스크류 압출기 (TSE) 에 대해 완전 충전 및 부분 충전 (스타브 피드) 조건 모두에서 시뮬레이션이 수행되었습니다.
  • 솔버들은 자유 표면의 진화, 충전 비율, 그리고 압력 장을 성공적으로 포착했습니다.
  • 결과는 물리적 기대와 일치했습니다: 완전 충전 시나리오 ($Q_{in} > Q_0$) 는 입구에서 출구로 향하는 음의 압력 구배를 보인 반면, 부분 충전 시나리오 ($Q_{in} < Q_0$) 는 출구 근처에서 압력 증가를 보여주었습니다.
  • IB 방법은 동적 메쉬 재생성 없이 TSE 의 복잡한 메쉬링 기하학을 처리하는 데 효과적이었습니다.

의의 및 주장
본 논문은 제안된 BC-VOF-IB 프레임워크가 학술적 CFD 연구와 산업용 고분자 가공 요구 사항 사이의 격차를 해소하는 의미 있는 단계라고 주장합니다. 높은 점도 차이와 관련된 수치적 불안정성을 극복하고 회전 기하학적 구조물에 대한 복잡한 몸체 적합 메쉬 생성의 필요성을 제거함으로써, 이 방법은 산업적으로 관련성이 있는 부분 충전 혼합 장치에서 속도와 압력 장의 물리적으로 일관된 예측을 가능하게 합니다.

저자들은 현재 연구가 등온 공정 및 일반화된 뉴턴 유체 역학으로 제한되어 있음을 인정합니다. 그들은 향후 개발이 산업용 고분자 혼합의 물리학을 완전히 포착하기 위해 열적 효과 (에너지 방정식) 와 점탄성 구성 모델을 포함해야 한다고 명시합니다. 그러나 현재 프레임워크는 이러한 확장을 위한 엄격한 기초를 마련하고 복잡하며 이동하는 기하학적 구조물에서의 자유 표면 유동을 시뮬레이션하기 위한 강력한 도구를 제공합니다.