잠근경계법

Immersed boundary method

계산 유체 역학에서, 잠근 경계 방법은 원래 1972년 찰스 페스킨이 유체 구조(섬유) 상호작용을 시뮬레이션하기 위해 개발한 접근방식을 가리켰다.[1] 구조물 변형과 유체 흐름의 결합을 처리하면 수치 시뮬레이션에 많은 난제가 발생한다(탄성 경계는 유체의 흐름을 변화시키고 유체는 탄성 경계를 동시에 이동한다). 잠근 경계법에서 유체는 오일러 좌표에 표시되고 구조는 라그랑 좌표에 표시된다. 압축할 수 없는 Navier에 의해 제어되는 뉴턴 유체용–스톡스 방정식, 유체 방정식은

그리고 압축할 수 없는 유체의 경우, 우리는 그 상태를 가지고 있다.

담근 구조물은 일반적으로 로 표기되는 1차원 섬유 모음으로 표현된다 섬유는 파라메트릭 곡선 ) 이고 t 이다. 섬유 물리학은 섬유 힘 F(, ) 를 통해 표현된다 이 용어에 스프링 힘, 굽힘 저항 또는 다른 유형의 행동이 내장될 수 있다. 그런 다음 구조물이 유체에 가하는 힘을 모멘텀 방정식의 소스 항으로 보간한다.

여기서 (는) Dirac Δ 함수다. 힘을 다차원까지 확장하여 탄성 표면이나 3차원 고형물을 모델링할 수 있다. 질량이 없는 구조를 가정할 때 탄성섬유는 국부 유체 속도에 따라 이동하며 델타 함수를 통해 보간될 수 있다.

여기서 전체 유체 영역을 나타낸다. 이러한 방정식의 탈바꿈은 유체에 오일러식 격자와 섬유에 별도의 라그랑지식 격자를 가정하여 수행할 수 있다. 평활기 기능에 의한 델타 분포의 근사치를 통해 우리는 두 그리드 사이를 보간할 수 있을 것이다. 기존의 모든 유체 용액 용액은 섬유 방정식을 위한 용액 용액에 결합하여 밀폐된 경계 방정식을 해결할 수 있다. 이 기본 접근법의 변형들은 유체 흐름과 상호작용하는 탄성 구조를 포함하는 다양한 기계적 시스템을 시뮬레이션하기 위해 적용되었다.

페스킨에 의한 이 방법의 본래의 개발 이후, 신체의 표면에 부합하지 않는 격자의 복잡한 몰입체 위의 흐름을 시뮬레이션하는 다양한 접근법이 개발되었다. 여기에는 몰입 인터페이스 방법, 데카르트 그리드 방법, 고스트 유체 방법, 컷 셀 방법 등의 방법이 포함된다. Mittal과 Iaccarino는[2] 이 모든 (및 기타 관련) 방법을 "몰입된 경계 방법"이라고 하며 이러한 방법의 다양한 범주화를 제공한다. 구현의 관점에서, 그들은 몰입된 경계 방법을 지속적 강제 방법과 이산적 강제 방법으로 분류한다. 전자의 경우 소멸하기 전에 연속적인 Navier-Stokes 방정식에 힘 용어가 추가되는 반면, 후자의 경우 소멸된 방정식에 힘을 가한다(명백히 또는 암묵적으로). 이 분류법 하에서 페스킨의 원래 방법은 연속적인 강제 방법인 반면, 카르테시아 그리드, 컷 셀, 고스트-유체 방법은 이산 강제 방법인 것이다.

참고 항목

소프트웨어: 숫자 코드

메모들

  1. ^ Peskin, Charles S (1972-10-01). "Flow patterns around heart valves: A numerical method". Journal of Computational Physics. 10 (2): 252–271. doi:10.1016/0021-9991(72)90065-4. ISSN 0021-9991.
  2. ^ 미탈 & 아이아카리노 2005.

참조