[논문 리뷰] Numerical simulation of salt migration -- Large deformation in viscoelastic solid bodies

I-Shih Liu, Rolci Cipolatti외 2명|arXiv (Cornell University)|2011. 08. 10.

Geological formations and processes참고 문헌 14인용 수 3

한 줄 요약

이 논문은 유동성 있는 점탄성 고체에서의 대변형을 정확하게 수치 모델링할 수 있도록, 업데이트된 기준 구성형태를 기반으로 한 순차적 선형 근사 방법을 제안한다. 이 방법은 현재 구성형태를 기준으로 하여 대변형을 작은 증분 단계의 집합으로 간주하며, 심지어 극한의 변형 조건에서도 최소한의 메쉬 재구성으로 장기간의 변형과 이동을 효과적으로 시뮬레이션할 수 있다.

ABSTRACT

We consider instability of a two layered solid body of a denser material on top of a lighter one. This problem is widely known to geoscientist in sediment-salt migration as salt diapirism. In the literature, this problem has often been treated as Raleigh-Taylor instability in viscous fluids instead of solid bodies. In this paper, we propose a successive linear approximation method for large deformation in viscoelastic solids as a model for salt migration.

연구 동기 및 목표

• 소금 이동을 고려한 점탄성 고체에서의 대변형을 시뮬레이션하는 데 도전하는 데 목적이 있다.
• 기존의 유체 기반 모델(예: 레일리-테일러 불안정성)이 소금 다이피르리즘에서 고체 특성을 반영하지 못하는 한계를 극복하는 데 목적이 있다.
• 실제 점탄성 재료 거동을 반영할 수 있는 안정적이고 정확한 다이피르 형성 시뮬레이션을 가능하게 하는 수치적 방법을 개발하는 데 목적이 있다.
• 극한의 변형 조건에서도 최소한의 메쉬 재구성으로 장기적인 다이피르 진화를 시뮬레이션할 수 있는지를 입증하는 데 목적이 있다.

제안 방법

• 각 시간 단계에서 현재 변형된 구성형태를 기준으로 선형화하는 상대적술적 기술을 사용한다.
• 각 단계에서 구성 함수는 현재 상태에서 평가되고, 소규모 증분 변형을 위해 선형화되며, 이로 인해 문제는 선형 경계값 문제의 열거로 변환된다.
• 유한요소법을 사용하여 각 단계에서 선형화된 방정식을 풀며, 경계 조건은 직접적으로 현재 구성형태에 적용된다.
• 모든 시간 단계에서 기준 상태를 업데이트함으로써 初기 구성형태 기반의 증분 데이터의 복잡성을 피한다.
• 이 방법은 2D 및 3D 영역 모두에 구현되었으며, 소금에는 무어니-리블린 유형의 점탄성 모델을, 침전물에는 탄성 매개변수를 사용하였다.
• 메쉬는 매 단계마다 업데이트되지만, 대변형 조건에서도 증분 접근의 안정성 덕분에 메쉬 재구성은 거의 필요로 하지 않는다.

실험 결과

연구 질문

• RQ1 순차적 선형 근사 방법은 소금 다이피르리즘과 관련된 점탄성 고체에서의 대변형을 효과적으로 시뮬레이션할 수 있는가?
• RQ2 제안된 방법은 기존의 유체 기반 레일리-테일러 불안정성 모델에 비해 다이피르 형성의 특성을 얼마나 잘 반영하는가?
• RQ3 극한의 변형 조건에서도 최소한의 메쉬 재구성으로 장기적인 다이피르 진화를 어느 정도까지 시뮬레이션할 수 있는가?
• RQ4 초기 불안정성과 기초 암석의 기울기가 점탄성 고체 모델에서 소금 구조의 기초 형성과 형태에 어떤 영향을 미치는가?

주요 결과

• 이 방법은 1억 년에 걸쳐 소금 다이피르의 형성을 성공적으로 시뮬레이션하였으며, 다이피르의 최대 높이에 도달하는 데 약 2,000만 년이 소요되었다.
• 시뮬레이션 결과에 따르면 밀도 역전에 기인한 부력이 다이피르리즘을 주도하며, 약 2,000만 년 후 다이피르는 성숙하고 안정된 상태에 이르렀다.
• n=300일 때(3,000만 년)의 변형된 메쉬는 실험적 실리콘 푸티 모델과 유사한 큰 변형 패tern을 보였다.
• 기초 암석이 1,000만 년 동안 점차 1도 기울여지자, 모델은 오른쪽에서 왼쪽으로 이르는 소금 베개와 다이피르의 연속적인 형성을 생성하였으며, 이는 독일 북부 지역의 지질학적 관측 결과와 일치하였다.
• 수치 결과는 독일 북부의 페르미아니안 소금 복합체에서 관측된 다중 다이피르 형성 구조를 정성적으로 재현하였으며, 모델의 지질학적 타당성을 검증하였다.
• 지질학적으로 校정되지 않은 편리한 재료 매개변수를 사용함에도 불구하고, 이 방법은 수치적 안정성을 유지하며 빈번한 메쉬 재구성 없이 물리적으로 타당한 변형 패턴을 생성하였다.