[논문 리뷰] A trust-region SQP method for the numerical approximation of viscoplastic fluid flow

Timm Treskatis, Miguel Moyers-González외 1명|arXiv (Cornell University)|2015. 04. 30.

Rheology and Fluid Dynamics Studies참고 문헌 34인용 수 3

한 줄 요약

이 논문은 점성플라스틱 유체 유동을 시뮬레이션하기 위해 응력 기반 이중 변수로 문제를 재정의함으로써 더 높은 정규성과 뛰어난 수렴성을 가능하게 하는 신뢰영역 순차적 정수계획법(TR-SSQP) 방법을 제안한다. 이 방법은 인위적 정규화나 메쉬 의존성 없이 허쉬엘-불크리 유동에서 최대 140배 빠른 CPU 시간을 달성하며, 널리 사용되는 증강 라그랑주(ALG2) 방법보다 뛰어난 성능을 보인다.

ABSTRACT

We present a new approach to the problem of stationary viscoplastic duct flow as modelled by the Herschel-Bulkley model, with Bingham fluids included as a special case. While the mathematical formulation of this problem is conventionally based on a variational inequality, or equivalently, on a nonsmooth minimisation problem for the flow velocity, we suggest an alternative approach. Considering the Lagrangian dual in terms of the stress, rather than the velocity, turns out to be advantageous in numerous ways. The objective functional possesses higher regularity, which ensures applicability of second order methods. Our numerical experiments with a trust-region SQP algorithm also demonstrate clearly superior performance compared to the widely used augmented Lagrangian method, although no artificial regularisation is introduced into the problem. Hence, besides providing a new theoretical angle to a classical problem, our results also pave the way for an entirely new class of numerical approaches to simulating flows of viscoplastic fluids.

연구 동기 및 목표

• 점성플라스틱 유체 유동을 시뮬레이션하는 데 발생하는 수치적 과제, 특히 유효 응력 거동으로 인한 비연속성과 악조건의 조건을 해결하기 위해.
• 기존 증강 라그랑주 방법의 한계를 극복하기 위해, 이는 느린 수렴성과 알고리즘 매개변수에 대한 민감성으로 인해 악영향을 받기 때문이다.
• 점성플라스틱 유동 문제의 이중 공식화에서 더 높은 정규성을 활용하는 안정적이고 이阶 수치 방법을 개발하기 위해.
• 인위적 점성도나 정규화 없이도 플라스틱 거동을 정확히 포착하는 새로운 수치 프레임워크를 제공하기 위해.
• 빙하 및 허쉬엘-불크리 흐름 시나리오에서 ALG2에 비해 상당한 계산 성능 향상을 입증하기 위해.

제안 방법

• 기존의 속도 기반 변분부등식 대신 응력 기반 이중 공식화를 사용하여 점성플라스틱 흐름 문제를 재정의함으로써 더 부드러운 목적 함수 기능을 도출한다.
• 유한요소법을 사용하여 이중 문제를 이산화하여 유한차원 최적화 문제를 도출한다.
• 결과로 얻어진 부드럽고 제약 조건이 있는 최적화 문제를 해결하기 위해 신뢰영역 순차적 정수계획법(TR-SSQP) 알고리즘을 적용한다.
• 신뢰영역 프레임워크를 사용하여, 헤시안 행렬 근사가 열악한 경우에도 전역 수렴성과 강건성을 확보한다.
• 정확한 선 탐색과 신뢰영역 갱신을 통합하여 내부 반복 없이 초선형 수렴성을 유지한다.
• 인위적 정규화나 점성을 피하여 점성플라스틱 거동의 물리적 정확성을 유지한다.

실험 결과

연구 질문

• RQ1 응력 변수를 기반으로 한 이중 공식화가 기존의 속도 기반 공식화에 비해 점성플라스틱 흐름 문제에서 더 부드러운 최적화 문제를 이끌 수 있는가?
• RQ2 신뢰영역 SQP 방법이 점성플라스틱 흐름 문제에서 증강 라그랑주(ALG2) 방법보다 더 빠른 수렴성과 우수한 성능을 보일 수 있는가?
• RQ3 Bingham 및 허쉬엘-불크리 유체에서 메쉬 정밀도 향상과 유효 응력 변화에 따라 제안된 TR-SSQP 방법의 성능은 어떻게 변화하는가?
• RQ4 인위적 정규화나 점성을 도입하지 않더라도 이 방법이 강건성과 효율성을 유지할 수 있는가?
• RQ5 유체의 유변학적 거동(예: 점성 감소 지수 α)이 TR-SSQP 알고리즘의 수렴 행동에 미치는 영향은 무엇인가?

주요 결과

• 허쉬엘-불크리 흐름에서 α = 1.5일 경우, TR-SSQP 방법은 ALG2 대비 최대 140배의 속도 향상을 달성하였다. 특히 더 미세한 메쉬에서 두드러진 성능 향상을 보였다.
• α = 1.5 및 τ₀ = 0.1 조건에서, TRS는 4308개 노드 메쉬에서 17초의 CPU 시간을 소비하였고, ALG2는 116.62초를 소비하여 114배의 속도 향상을 기록하였다.
• 메쉬 정밀도 향상과 낮은 유효 응력 조건에서 TRS 방법은 반복 횟수가 감소하는 경향을 보이며, 더 도전적인 조건에서도 효율성이 향상됨을 나타냈다.
• 비대칭 기하구조에서도 우수한 성능 유지를 보였으며, 키위 파이프 기하구조에서 최대 6.7배의 속도 향상을 기록하였다.
• TRS와 ALG2 해의 상대적 차이는 항상 1e-3 이하로 유지되어, 더 빠른 수렴에도 불구하고 해의 정확성이 확인되었다.
• TRS 방법은 메쉬 독립성과 강건한 수렴성을 유지하였으며, ALG2와 달리 매개변수 조정이나 내부 반복이 필요로 하지 않았다.