[논문 리뷰] Micro-Macro Modelling of an Array of Spheres Interacting Through Lubrication Forces

Aline Lefebvre-Lepot, Bertrand Maury외 1명|arXiv (Cornell University)|2008. 02. 19.

Adhesion, Friction, and Surface Interactions인용 수 10

한 줄 요약

이 논문은 라우브리케이션 힘을 통해 상호작용하는 강체 구의 밀도 농도에서 거시적 구성 모델을 제안하며, 국소 유체 분율에 반비례하는 확장 점성도를 가지는 뉴턴형 유변학적 방정식을 유도한다. 주요 기여는 입자 접촉과 군집 내 무한한 국소 점성도를 允허하는 상황에서도 이산 해가 연속 편미분방정식에 약한 수렴함을 엄밀히 증명한 것이다.

ABSTRACT

We consider here a discrete system of spheres interacting through a lubrication force. This force is dissipative, and singular near contact: it behaves like the reciprocal of interparticle distance. We propose a macroscopic constitutive equation which is built as the natural continuous counterpart of this microscopic lubrication model. This model, which is of the newtonian type, relies on an elongational viscosity, which is proportional to the reciprocal of the local fluid fraction. We then establish the convergence in a weak sense of solutions to the discrete problem towards the solution to the partial differential equation which we identified as the macroscopic constitutive equation.

연구 동기 및 목표

• 강체 구가 라우브리케이션 힘을 통해 상호작용하는 이산 시스템의 효과적 유변학적 거동을 기술하는 연속 거시적 모델을 유도하기 위해.
• 접촉 근처에서 라우브리케이션 힘이 특이해지는 문제를 다루며, 입자 간격이 0에 수렴할수록 점성도가 발산하는 상황을 고려하기 위해.
• 효과적 점성도가 군집에서 무한대가 되는 경우에도, 이산 해가 연속 편미분방정식의 해로 약한 위상에서 수렴함을 확립하기 위해.
• 기존의 균질화 프레임워크를 일반화하여, 무한한 점성도(완전 밀도 영역)를 가진 비유계 확산 계수를 允허함으로써, 고전적 균질화 이론의 표준 가정을 위반하는 상황을 다루기 위해.
• 시간에 의존하는 입자 구성으로의 모델 확장을 통해 속도장과 고체 분율에 대한 운반 방정식을 결합하기 위해.

제안 방법

• 이웃하는 구 사이의 상호작용을 입자 간격의 역수에 비례하는 라우브리케이션 힘으로 모델링하며, 단순화를 위해 비례 계수(κ = 1)를 사용한다.
• 직선상에 위치한 N+1개의 구로 이루어진 이산 시스템을 정의하며, 위치 qᵢ(t), 속도 uᵢ(t), 고체 분율 ρᵢ(t)를 고려한다. 여기서 ρᵢ = 1은 접촉을 의미한다.
• 상대 속도와 입자 간격을 이용해 이산 기울기를 정의함으로써, 조각별 상수 속도장 u^ε 및 그의 약한 의미에서의 도함수를 구성한다.
• 외부 체력 f가 존재하는 형태의 연속 근사인 힘 평형을 나타내는 거시적 편미분방정식 −∂ₓ[(1−ρ)⁻¹ ∂ₓu] = ρf 를 도출한다.
• 입자 운동에 따른 고체 분율의 진화를 기술하기 위해 운반 방정식 ∂ₜρ + ∂ₓ(ρu) = 0 을 도입한다.
• 속도 기울기와 시간 도함수에 대한 균일한 유계성 조건을 이용해, 이산 속도 및 밀도 장의 약한 수렴을 증명한다.

실험 결과

연구 질문

• RQ1 라우브리케이션 힘이 특이한 상황에서 상호작용하는 이산 구 입자 시스템으로부터 연속 거시적 모델을 어떻게 도출할 수 있는가?
• RQ2 입자 간격이 0에 수렴하고 라우브리케이션 힘이 지배적인 상황에서, 매크로스코픽 륨성 거동을 기술하는 효과적 구성관계는 무엇인가?
• RQ3 고체 분율이 1이 되어 국소 점성도가 무한대가 되는 상황에서도 이산 해가 연속 편미분방정식으로 수렴할 수 있는가?
• RQ4 이 모델은 입자 군집화와 접촉을 어떻게 처리하며, 결과적으로 유도된 특이한 편미분방정식의 수학적 구조는 어떠한가?
• RQ5 관성 항을 포함하거나 유체 점성도가 0으로 수렴할 경우, 시스템의 거동은 어떻게 되는가?

주요 결과

• 거시적 모델은 특이한 타원형 편미분방정식 −∂ₓ[(1−ρ)⁻¹ ∂ₓu] = ρf 로 표현되며, 이는 효과적 확장 점성도가 국소 유체 분율 (1−ρ)에 반비례함을 나타낸다.
• 입자 접촉(ρ = 1)을 허용하여 국소 점성도가 무한대가 되며, 연속적 수식은 매크로스코픽 군집을 자연스럽게 수반한다.
• 속도 및 기울기의 유계성 조건 하에, 이산 해의 수열이 연속 편미분방정식의 해로 약한 수렴함이 엄밀히 증명된다.
• 상대 속도 표현식에 (1−ρᵢ) 인자가 존재함으로써, 접촉 근처에서도 속도 기울기가 ε에 대해 균일하게 유계로 유지되어 폭발을 방지한다.
• 시간에 의존하는 시스템으로의 확장은 ∂ₜρ + ∂ₓ(ρu) = 0 과 ∂ₜ(ρu) + ∂ₓ(ρu²) − ∂ₓ[(1−ρ)⁻¹ ∂ₓu] = ρf 의 쌍방정식을 유도하며, 특이한 점성도를 가진 1차원 압력 없는 나비에-스토크스 시스템과 유사하다.
• 유체 점성도가 0으로 수렴하는 극한(ε→0)에서, 이 모델은 압력처럼 행동하는 장과 과거 역사를 추적하는 장 γ를 포함하는 접착성 입자 시스템을 시사하며, 저점성 영역에서의 접촉 피하기 이론과 일관된다.