[논문 리뷰] Reply to "Comment on Sochi's variational method for generalised Newtonian flow" by Pritchard and Corson

Taha Sochi|arXiv (Cornell University)|2015. 05. 29.

Rheology and Fluid Dynamics Studies참고 문헌 8인용 수 2

한 줄 요약

이 논문은 티하 소치의 일반화된 뉴턴 유체의 원형 관 내 유동에 대한 변분법이 프리처드와 코슨의 비판에 대해 변론한다. 변분 형식이 총 응력 최소화에 기반하여 수학적으로 타당하며, 제곱법 유체뿐 아니라 모든 비점성플라스틱 일반화된 뉴턴 유체에 대해 유효하다고 주장한다. 이를 위해 오일러-라그랑주 및 딜리클레 원리를 사용하여 핵심 방정식을 재유도하고, 운동 방정식과의 동치성을 입증하며, 본질적 한계나 단순성에 대한 주장은 기각한다.

ABSTRACT

In this article we challenge the claim that the previously proposed variational method to obtain flow solutions for generalized Newtonian fluids in circular tubes and plane slits is exact only for power law fluids. We also defend the theoretical foundation and formalism of the method which is based on minimizing the total stress through the application of the Euler-Lagrange principle.

연구 동기 및 목표

• 프리처드와 코슨이 변분법이 제곱법 유체에 대해서만 정확하다고 주장한 데에 대응하기 위해.
• 1차원 유동계에서 응력 최소화에 기반한 변분 원리의 수학적 및 물리적 타당성을 확립하기 위해.
• 변분 형식(식 (8))이 제곱법 유체에 국한되지 않고 비점성플라스틱 일반화된 뉴턴 유체 전반에 대해 적용 가능하다는 것을 보여주기 위해.
• 변분 형식(오일러-라그랑주 대비 딜리클레) 간의 차이를 명확히 하고, 한 형식의 실패가 원리 전체를 무효화하지는 않는다는 것을 논증하기 위해.
• 식 (1)에서 부분 도함수를 사용하는 것이 본질적인 문제인지, 아니면 일관된 적용 맥락 하에서는 표기 문제일 뿐인지 밝히기 위해.

제안 방법

• 응력 균형 dτ/dr = G 에서 출발하여 오일러-라그랑주 체계를 사용해 핵심 방정식 d/dr(μ dγ/dr) = 0 을 재유도한다.
• 함수 f = μγ 에 오일러-라그랑주 방정식을 적용하여 d/dr(γ dμ/dr) = 0 을 도출하고, 이가 원래 형태와 수학적으로 동치임을 보여준다.
• 두 형태( d/dr(γ dμ/dr) = 0 과 d/dr(μ dγ/dr) = 0 )가 서로 동치임을 입증하며, 그 합이 일정하므로 각각이 시스템의 행동을 독립적으로 포괄함을 보여준다.
• 이전 연구 [2]에서 사용한 딜리클레 기반 변분 접근법을 적용하여 다양한 륨탈로지에 대해 정확한 수치해를 도출할 수 있음을 보여준다.
• 어떤 유도 과정에서 나타나는 명백한 단순성(예: 0=0)은 원리 자체의 결함이 아니라 중복된 연산의 결과임을 주장한다.
• 변분형식과 운동형식 간의 동치성을 입증하여, 둘 다 동일한 물리적 시스템을 표현하고 있음을 보여준다.

실험 결과

연구 질문

• RQ1 변분 형식 d/dr(μ dγ/dr) = 0 이 비점성플라스틱 일반화된 뉴턴 유체 전반에 대해 유효한가, 아니면 제곱법 유체에만 국한되는가?
• RQ2 응력 최소화에 기반한 변분 원리는 딜리클레 원리 또는 기타 방법을 통해 수학적으로 정당화될 수 있는가?
• RQ3 변분 방정식에서 부분 도함수를 사용하는 것이 원리의 타당성을 해치는가, 아니면 표기 문제일 뿐인가?
• RQ4 복잡한 륨탈로지에 대해서도 변분법이 정확한 해를 도출할 수 있는가, 비록 수치적 방법이 필요하더라도?
• RQ5 변분형식과 운동형식 간의 명백한 동치성은 수학적 우연인가, 아니면 더 깊은 물리적 일관성의 반영인가?

주요 결과

• 변분 형식 d/dr(μ dγ/dr) = 0 은 운동 방정식 dτ/dr = G 와 수학적으로 동치이며, 양자 모두 전체 시스템 역학을 표현한다.
• 두 형태( d/dr(γ dμ/dr) = 0 과 d/dr(μ dγ/dr) = 0 ) 간의 동치성은 그 합이 일정하기 때문에 발생하며, 이는 각각이 시스템 행동을 독립적으로 포괄함을 보장한다.
• 이 방법은 제곱법 유체에 국한되지 않으며, 논문은 원형 관과 평판 슬릿에서 9종의 다양한 륨탈로지 모델에 대해 그 타당성을 입증한다.
• 식 (1)에서 부분 도함수를 사용하는 것은 일관된 적용 맥락 하에서는 본질적 결함이 아니며, 방정식이 정확히 적용된다면 문제가 되지 않는다.
• 0 = 0 으로 이어지는 유도 과정은 타당하지만 생산적이지 않은 수학적 연산이며, 변분 원리의 비일관성의 징후는 아니다.
• 논문 [2]에서 사용한 딜리클레 기반 변분 접근법은 오일러-라그랑주 형식과 무관하게 일반화된 뉴턴 유체에 대해 정확한 수치해를 도출할 수 있다.