A Jacobian-free Newton-Krylov method for cell-centred finite volume solid mechanics

기술 요약: 셀 중심 유한 체적 고체 역학을 위한 자코비안 프리 뉴턴-크릴로프 방법

문제 정의
고체 역학을 위한 전통적인 유한 체적(FV) 방법은 주로 분리형 암시적 접근 방식(Picard 반복법)이나 완전 명시적 기법에 의 eyes 의존합니다. 분리형 방법은 구현의 단순성과 메모리 효율성을 제공하지만, 기하학적 비선형성, 큰 회전 또는 복잡한 결합이 포함된 문제에서 수렴 속도가 느려지는 경우가 많습니다. 반대로, 유한 요소 커뮤니티에서 흔히 사용되는 뉴턴-라프슨 유형의 방법은 이차 수렴(quadratic convergence)을 제공하지만, 전체 자코비안 행렬을 명시적으로 조립, 저장 및 역행렬 계산하는 과정이 필요합니다. 이러한 요구 사항은 계산 비용과 메모리 집약적이며, 기존의 분리형 FV 프레임워크를 블록 결합 뉴턴 방법으로 적응시키는 것은 매우 까다롭습니다. 결과적으로, 대대적인 코드 리팩토링 없이도 기존의 분리형 FV 프레임워크에 통합될 수 있는 견고하고 효율적인 블록 결합 솔버가 필요합니다.

방법론
저자들은 셀 중심 유한 체적 고체 역학에 맞춤화된 자코비안 프리 뉴턴-크릴로프(JFNK) 접근 방식을 제안합니다. 핵심 방법론은 유한 차분법을 사용하여 벡터($v$)에 대한 자코비안 행렬($J$)의 작용을 근사함으로써 자코비안 행렬을 명시적으로 형성하는 것을 피합니다:
$$Jv \approx \frac{R(u + \epsilon v) - R(u)}{\epsilon}$$
여기서 $R$은 잔차 함수, $u$는 해 벡터, $\epsilon$은 섭동 스칼라입니다. 이 근사법은 뉴턴 반복 내에서 선형 시스템을 해결하기 위해 일반화 최소 잔차(GMRES) 크릴로프 부공간 솔버와 결합됩니다.

효율성과 견고성을 보장하기 위해, 본 방법은 프리컨디셔닝(preconditioning) 전략을 채택합니다. 정확한 자코비안을 사용하는 대신, 저자들은 관성 항과 확산 유사 항(분리형 솔버에서 사용되는 Rhie-Chow 안정화에 기반함)으로부터 유도된 컴팩트 스텐실 근사 자코비안을 활용합니다. 이 프리컨디셔너($P$)는 개념적으로 기존의 분리형 접근 방식에서 사용되는 행렬과 유사하며, 이를 통해 새로운 방법이 OpenFOAM의 solids4foam 툴박스와 같은 프레임워크 내의 기존 행렬 조립 및 저장 루틴을 활용할 수 있게 합니다.

구현에는 비선형 솔버(SNES) 및 선형 대수를 위한 PETSc 라이브러리가 사용되었으며, solids4 foam FV 솔버와 인터페이스되었습니다. 주요 알고리즘 특징은 다음과 같습니다:

• 프리컨디셔너: LU 직접 솔버(MUMPS), 불완전 LU(ILU), 그리고 대수적 다중 격자(HYPRE BoomerAMG)의 비교.
• 전역화(Globalisation): 해를 이차 수렴 영역으로 유도하기 위한 라인 서치(line search) 및 과도 지속(transient continuation, 예측자-교정자)의 사용.
• 안정화: 제로 에너지 모드(checkerboarding)를 억제하기 위해 잔차에 Rhie-Chow 유형의 항이 추가되었으며, 수렴 속도를 위해 튜닝된 전역 스케일링 인자($\alpha$)가 적용되었습니다.

주요 기여

1. 유한 체적 고체 역학에 대한 첫 적용: 본 연구는 유한 체적 고체 역학을 위해 특별히 제안된 최초의 JFNK 방법으로, 유체 역학 응용 분야와 고체 역학 사이의 간극을 메웁니다.
2. 원활한 통합: 컴팩트 스텐실 근사 자코비안을 프리컨디셔너로 사용함으로써, 최소한의 코드 수정만으로 기존의 분리형 FV 프레임워크 내에서 뉴턴 유형의 수렴을 채택할 수 있게 합니다.
3. 포괄적인 벤치마킹: 차원(2-D/3-D), 재료 거동(선형 탄성, 초탄성, 탄소성), 물리 현상(정적, 동적, 유체-구조 상호작용)이 다양한 일련의 케이스를 통해 본 방법을 엄격하게 테스트했습니다.

결과
본 연구는 계산 시간, 메모리 사용량 및 견고성 측면에서 JFN키(JFNK) 방법을 분리형 접근 방식과 벤치마킹합니다:

• 효efficiency: JFNK 방법은 모든 선형 및 비선형 탄성 케이스에서 분리형 접근 방식을 능가했습니다. 많은 사례에서 1~2 자릿수의 속도 향상이 관찰되었으며, 2-D 선형 탄성 Cook's membrane 케이스에서는 최대 341배의 속도 향상이 나타났습니다.
• 메모리: 자유도가 $10^4$ 미만인 문제의 경우, 메모리 사용량은 분리형 방법과 유사합니다. 더 큰 문제의 경우, JFNK의 메모리 사용량은 증가하여 구형 공동(spherical cavity) 케이스에서 분리형 방법보다 최대 3.6배에 달했는데, 이는 주로 GMRES 솔버와 다중 격자 프리컨디셔너의 저장 요구 사항 때문입니다.
• 견고성 및 한계:
  • JFNK 방법은 분리형 방법이 수렴에 실패한 초탄성(큰 변형/회전) 케이스에서 우수한 견고성을 입증했습니다.
  • 그러나 JFNK 방법은 탄소성(소성 변형이 포함된 Cook's membrane) 케이스에서 발산했습니다. 저자들은 이를 소성 변형이 큰 경우에 대해 소형 변형 탄성에 기반한 컴팩트 스텐십 프리컨디셔너가 자코비안의 불충분한 근사치가 되었기 때문이라고 분석했습니다.
• 프리컨디셔너 선택:
  • LU 직접 솔버는 소-중규모 크기의 문제와 2-D 케이스에서 가장 빨랐습니다.
  • **다중 격자(AMG)**는 더 큰 3-D 문제에서 더 효과적이었으며, LU보다 더 나은 확장성과 낮은 메모리 증가율을 보였습니다.
  • **ILU(k)**는 일반적으로 가장 느리고 견고성이 떨어졌습니다.
• 안정화: Rhie-Chow 안정화 인자($\alpha$)를 높이는 것은 일반적으로 수렴 속도와 견고성을 개선했지만, 이는 전-점근적(pre-asymptotic) 정확도의 희생을 수반했습니다.

의의 및 주장
본 논문은 자코비안 프리 뉴턴-크릴로프 방법이, 특히 최소한의 수정만으로 도입이 가능한 기존의 분리형 프레임워크를 위한 유한 체적 고체 역학의 유망한 발전임을 주장합니다. 연구는 JFNK가 전체 자코비안 조립의 과도한 메모리 비용 없이도 탄성 문제에서 상당한 속도 향상과 개선된 견고성을 달 수 있음을 보여줍니다.

저자들은 현재의 프리컨디셔닝 전략 하에서 해당 방법이 아직 탄소성에 대해 견고하지 않다는 점을 명시하며 현재의 한계에 대해 겸허한 태도를 유지하고 있습니다. 그들은 향상된 소성 프리컨디셔너 개발과 비직교 보정(non-orthogonal corrections)의 포함을 향후 과제로 식별했습니다. 본 구현은 커뮤니티의 탐색과 확장을 장려하기 위해 solids4foam 툴박스에 공개되었습니다.