[논문 리뷰] Boundary conditions for a capillary fluid in contact with a wall

Henri Gouin, Witold Kosiński|arXiv (Cornell University)|2008. 02. 14.

Lattice Boltzmann Simulation Studies참고 문헌 7인용 수 19

한 줄 요약

이 논문은 제2계층 연속체 모형과 가상일 원리(가상일 원리)를 사용하여 경계면에서 경계 조건을 유도한다. 이는 경계면에서 유체 밀도, 그 법선 미분, 벽의 곡률이 함께 상호작용하여 경계면 응력을 결정하며, 유체가 벽에 축적되거나 감소하는 데 있어 분기점이 되는 임계 밀도 ρᵢ = γ₁/γ₂를 규명한다.

ABSTRACT

Contact of a fluid with a solid or an elastic wall is investigated. The wall exerts molecular forces on the fluid which is locally strongly nonhomogeneous. The problem is approached with a fluid energy of the second gradient form and a wall surface energy depending on the value of the fluid density at the contact. >From the virtual work principle are obtained limit conditions taking into account the fluid density, its normal derivative to the wall and the curvature of the surface

연구 동기 및 목표

• 제2계층 연속체 역학을 확장하여 고체 벽에서의 캐필라리 유체에 대해 일관된 경계 조건을 수립하기 위해.
• 유체 밀도에 의존하는 벽 표면 에너지를 변분 프레임워크에 통합하기 위해.
• 유체 경계면의 기계적 평형에서 경직성 벽과 탄성 벽의 반응을 모두 고려하기 위해.
• 유체 밀도, 그 법선 기울기, 벽 곡률 간의 관계를 유도하기 위한 조건을 도출하기 위해.
• 가상일 원리를 통해 동적 경우로 모델을 일반화하여 핵심 경계 조건을 유지하기 위해.

제안 방법

• 유체 에너지를 밀도 ρ와 그 기울기에 의존하는 제2계층 기능으로 기술하며, 벽 표면에 표면 에너지 B(ρ)를 포함한다.
• 전체 에너지(유체, 위치 에너지, 벽 표면 기여)에 대해 가상일 원리를 적용한다.
• 체적 적분의 변분을 통해 본체 평형 방정식을 도출하여 −div σ + ρ ∂Ω/∂x = 0을 얻는다.
• 경직성 벽의 경우 운동학적 제약 조건(n*ζ = 0)을 적용하고, 탄성 벽의 경우 법선 방향 이동을 허용한다.
• 표면 에너지의 변분을 통해 경계 조건을 유도하며, 밀도의 법선 미분과 응력 균형 조건을 도출한다.
• 경직성 벽의 경우 라그랑주 승수를 도입하고, 탄성 벽의 경우 응력 벡터 T를 도입하여 경계면 힘과 벽의 변형을 연결한다.

실험 결과

연구 질문

• RQ1 제2계층 유체 모형과 벽 표면 에너지가 함께 유체-벽 경계면에서 경계 조건을 어떻게 결정하는가?
• RQ2 유체 밀도와 그 법선 기울기가 경계면 응력과 적촉 거동을 결정하는 데 어떤 역할을 하는가?
• RQ3 벽 곡률이 유체-벽 경계면에서 응력 균형에 어떻게 영향을 미치는가?
• RQ4 경계 조건이 경직성 벽과 탄성 벽 간에 캐필라리 유체 모형에서 어떻게 다름을 보이는가?
• RQ5 임계 유체 밀도 ρᵢ = γ₁/γ₂가 벽에서의 유체 축적 또는 감소 맥락에서 물리적으로 어떤 의미를 갖는가?

주요 결과

• 경직성 벽에서 경계 조건 G − A = 0는 유체와 벽 에너지 기여 간의 일관성을 보장한다.
• 응력 균형 조건 2(A−H)/Rₘ n + gradₜg(A−H) + σ n = λ n은 곡률, 표면 에너지, 경계면 응력을 연결한다.
• 탄성 벽의 경우 법선 응력 성분은 Tₙ = P + 2B/Rₘ + (1/C)(B′ₚ)²이며, 탄성 응력 성분은 Tₜg = 0이다.
• 조건 C dρ/dn = γ₁ − γ₂ρ는 ρₛ < ρᵢ 또는 ρₛ > ρᵢ 여부에 따라 벽 근처에서 유체 밀도가 증가하거나 감소하는지를 결정한다.
• 임계 밀도 ρᵢ = γ₁/γ₂는 분기점이 되며, 이 이하일 경우 유체 감소가 발생하고, 이 이상일 경우 유체 축적이 발생한다.
• 경계 조건의 변화 없이 동적 경우로 일반화할 수 있으며, 이로 인해 해의 물리적 구조가 유지된다.