Core Concepts
바셋 기억항을 포함하는 맥시-라일리 방정식을 효율적으로 해결하기 위한 수치 해법을 제안한다. 기존의 직접 적분 방식과 프라사스 등이 제안한 편미분 방정식 변환 방식을 개선하여 정확도와 계산 효율성을 향상시켰다.
Abstract
이 논문은 유한 크기의 구형 입자가 유체 내에서 운동하는 과정을 모델링하는 맥시-라일리 방정식(MRE)을 효율적으로 해결하는 수치 해법을 제안한다.
MRE에는 입자의 과거 궤적에 의한 영향을 모델링하는 바셋 기억항이 포함되어 있어, 이를 직접 다루기 어렵다. 최근 프라사스 등은 이 항을 분수 미분으로 표현하여 편미분 방정식 형태로 변환하는 방법을 제안했다.
이 논문에서는 프라사스 등의 방법을 기반으로 하여, 유한 차분 기법을 이용한 새로운 수치 해법을 제안한다. 2차 및 4차 정확도의 유한 차분 이산화를 수행하고, 암시적-명시적 룽게-쿠타 시간 적분법을 적용하여 효율성을 높였다.
제안된 유한 차분 기반 해법과 프라사스 등의 다항식 전개 방식, 그리고 기존의 직접 적분 방식을 다양한 유동장 문제에 적용하여 비교 분석하였다. 정확도와 계산 효율성 측면에서 각 방법의 장단점을 확인하였다.
Stats
입자의 밀도가 유체 밀도의 2/3배일 때(R=7/9), 정상 와류 유동장에서 입자 궤적의 최대 상대 오차는 약 10^-12이다.
진동 유동장에서 입자 밀도가 유체 밀도와 같을 때(R=1), 정상 상태에 도달하는 입자 위치의 최대 상대 오차는 약 10^-13이다.
정지 유동장에서 입자 밀도가 유체 밀도의 3/2배일 때(R=4/3), 입자 속도 감소 과정의 최대 상대 오차는 약 10^-12이다.
Quotes
"바셋 기억항을 포함하는 맥시-라일리 방정식을 효율적으로 해결하기 위한 수치 해법을 제안한다."
"기존의 직접 적분 방식과 프라사스 등이 제안한 편미분 방정식 변환 방식을 개선하여 정확도와 계산 효율성을 향상시켰다."