[논문 리뷰] A Numerical Formulation for Simulating Free-Surface Hydrodynamics

Thomas T. O’Shea, Kyle A. Brucker외 2명|arXiv (Cornell University)|2014. 10. 07.

Fluid Dynamics Simulations and Interactions참고 문헌 7인용 수 23

한 줄 요약

이 논문은 선박 주변의 자유면 유동, 특히 테일스턴 흐름을 시뮬레이션하기 위해 카르테시안 격자 유한체적 방법과 침수체 기법, 체적분율(VOF) 기법을 결합한 방법을 제시한다. 이 방법은 파도의 붕괴, 스프레이, 공기 혼입을 정확하게 포착하며, 다양한 흐름 속도에서의 항력 및 파도 절단 측정치와 수치 예측 간 강한 일치를 보인다.

ABSTRACT

Cartesian-grid methods in combination with immersed-body and volume-of-fluid methods are ideally suited for simulating breaking waves around ships. A surface panelization of the ship hull is used as input to impose body-boundary conditions on a three-dimensional cartesian grid. The volume-of-fluid portion of the numerical algorithm is used to capture the free-surface interface, including the breaking of waves. The numerical scheme is implemented on a parallel computer. Recent improvements to the numerical scheme are discussed, including implementation of a new multigrid procedure and conversion to MPI communication. Numerical predictions are compared to laboratory measurements of a towed transom-stern model.

연구 동기 및 목표

• 파도 붕괴 및 공기 혼입을 포함한 복잡한 자유면 유동을 시뮬레이션하기 위한 견고하고 즉시 사용 가능한 CFD 프레임워크 개발
• 신체에 맞는 격자를 사용하지 않기 위해 침수체 경계 조건을 적용한 카르테시안 격자 사용
• 새로운 다중격자 솔버와 MPI 기반 통신을 통해 대규모 시뮬레이션에서 수치적 효율성과 확장성 향상
• 다양한 흐름 속도에서 테일스턴 모델의 세부 실험 측정치와 수치 예측의 검증
• 수십억 개 격자 포인트를 포함한 대규모 해양 유체역학 문제의 고해상도 시뮬레이션 가능화

제안 방법

• 압축성 없는 스토크스 방정식을 해결하기 위해 분할단계 시간 적분 기법을 사용하는 카르테시안 격자 스태그너드 격자
• 체적분율(VOF) 방법을 통해 셀 단위의 체적 분율 함수를 사용하여 자유면 계면을 추적하며, 날카운 액체-공기 계면을 표현하기 위해 부드러운 스텝 함수를 사용
• 표면 패널화를 입력으로 사용하여 체적력 방법을 통해 선체에 대한 슬립 없음 및 침투 금지 조건을 침수체 경계 조건으로 구현
• 분할단계 방법에서 압력 푸아송 방정식의 수렴 속도를 향상시키기 위해 새로운 다중격자 솔버 구현
• 고성능 컴퓨팅 시스템에서의 병렬 확장성 향상을 위해 MPI 기반 통신 루틴 통합
• 난류 및 점성 응력에 대해 미세구조 모델을 사용하며, SGS 응력 텐서는 동적 절차를 통해 모델링

실험 결과

연구 질문

• RQ1 VOF 및 침수체 기법을 적용한 카르테시안 격자 방법이 선박 모델 주변의 붕괴 파도 및 테일스턴 흐름을 정확하게 시뮬레이션할 수 있는가?
• RQ2 다양한 흐름 속도에서의 항력 및 파도 절단 프로파일 수치 예측이 실험 측정치와 얼마나 잘 일치하는가?
• RQ3 새로운 다중격자 솔버가 대규모 자유면 유동 CFD 시뮬레이션에서 수렴 속도를 얼마나 향상시키는가?
• RQ4 낙하 파도의 자유면 높이 및 수입자 속도 수치 예측이 실험 데이터와 얼마나 일치하는가?
• RQ5 고주파수 표면 스펙트럼 특징을 해석하는 데에는 어떤 제약이 있으며, 이는 어림표현 또는 모델 단순화와 어떤 관련이 있는가?

주요 결과

• 테일스턴 모델의 항력 예측값은 DTMB 실험 측정치와 네 가지 흐름 속도에서 뛰어난 일치를 보였다.
• NFA 코드로 예측한 파도 절단 프로파일은 네 가지 속도 및 횡방향 위치에서 실험 데이터와 매우 밀접하게 일치하였으며, 특히 후미부에서 일관된 정확도를 확보하였다.
• 낙하 파도 시뮬레이션의 자유면 높이 예측값은 형상과 크기 모두에서 실험 측정치와 잘 일치하였다.
• 이 모델은 특히 중간 속도(7–8노트)에서 파도 붕괴, 스프레이 형성, 공기 혼입과 같은 복잡한 유동 특징을 성공적으로 포착하였다.
• R/V Athena 뒤의 로스터테일 영역에서 평균 자유면 높이 예측값은 양호한 일치를 보였지만, 고주파수에서는 차이가 나타났으며, 이는 어림표현 또는 수치 모델의 부착물 누락 가능성과 관련이 있을 수 있다.
• 새로운 다중격자 솔버와 MPI 통신 루틴의 구현은 솔버 수렴 속도와 병렬 확장성 향상에 크게 기여하여, 수십억 개 셀에 가까운 격자 해상도의 시뮬레이션을 가능하게 하였다.