Vortex의 모습(유체역학 교과서에 꼭 나오는 사진)

비행기가 지나가면 날개로 인해 vortex가 만들어 진다



Vortex sheet의 운동방정식인 Birkhoff-Rott방정식은 singular integral의 형태여서 수치해석하기가 쉽지 않다. regular한 미분방정식들은 대단한 수치해석기법없이 정밀한 numerical solution을 얻지만 Birkhoff-Rott과 같은 singular integral들은 수치해석에 상당히 다양한 기법이 들어간다. 일단 point vortex method가 한 방법인데 이 방법을 적용하기 위해 많은 테크닉들이 들어간다. 일단 desingularizaition parameter인 Krasny delta와 Fourier filtering technique을 활용해서 roll-up 문제를 해결한 Krasny는 이 공로로 ICM에서 초청연설을 한 적이 있고 이 분야의 세계적인 대가로 이름을 날리게 되었다. 하지만 이것만으로는 여전히 Long time computation에 어려움이 많았는데 이유는 컴퓨터 계산은 횟수가 누적될 수록 에러도 누적되기 때문이고 vortex sheet자체가 가진 물리적인 불안정성(Kelvin-Helmholtz instability)때문인데 이 에러를 줄이고 빠르게 계산하기 위해 많은 기법이 개발되었다. Birkhoff-Rott 방정식을 Krasny delta를 사용해 regular하게 만든 후 컴퓨터 연산을 빠르게 하기 위해 변수변환을 통해 식을 변형하고 이것을 direct computation하지 않고 quadrature들을 활용해 빠르게 계산하는 방법을 사용한다. 물론 대부분 Runge-Kutta 4th order method를 활용하고.


즉, 기본적으로 다음과 같은 기법들을 사용한다.


0) point vortex method

1) Desingularization by Krasny delta

2) Runge-Kutta 4th order method

3) Change of variables

4) Quadrature

5) Fourier filtering (periodic boundary condition인 경우)


이것들은 일단 기본이고 좀 더 정확하고 오랜시간에 걸쳐 계산하기 위해서 김선철 교수님의 논문(이준엽,손성익 공저)에서는 point insertion과 redistribution 기법을 활용했는데 이를 위해 cubic spline method까지 사용했다.


참고로 이론적인 배경은 생각하지 않고 계산기법만 따진다면 수학과 학부 2학년을 마친 학생들은 대부분 이해할 수있는 내용들이다.


첨부한 파일은 Krasny의 당시 논문인데, 이 분야의 중요한 논문이다. 아무도 해결못한 roll-up을 해결했으니 당연하다.

A study of singularity formation in a vortex sheet by the point-


이 논문의 특징은 친절함인데, vortex sheet 시뮬레이션의 어려움을 상세히 설명하고 그 어려움들을 하나하나 어떻게 하면 극복할 수 있는지 자세하게 설명한다.



해양에서 속도가 다른 두 해류가 만나 경계에서 vortex sheet 생기고 roll-up이 일어나면서 섞이고 있다




  1. *저녁노을* 2014.02.13 08:40 신고

    잘 보고갑니다.
    ㅎㅎ
    즐거운 하루 되세요

  2. 주영 2014.04.28 18:27 신고

    좋은 블로그네요...
    지금 singular integral 배우고있는데 여기에 대한 기초적인 교재같은걸 올려주시면 정말 감사하겠습니다.
    이를 이용한 미분방정식 풀이에 대해 포스팅해주시면 더 좋구요...
    그럼 좋은 하루 되세요...

    • '주영'님 죄송하게도 제가 지금 singular integral을 해석학적으로 깊이 공부하지 않았기 때문에 별로 도움을 드릴 수 없습니다. 제가 아직 배움이 짧다보니 이런 일이 생기는군요.. 앞으로 더 분발하겠습니다.ㅠㅠ

  3. 주영 2014.04.28 18:36 신고

    그 미분방정식..Cauchy 's equation에 대해서 올려주시면 좋겠습니다.
    이해하기 좋게 설명을 잘 해주시는것같아 부탁드립니다.

  4. 주영 2014.04.28 18:46 신고

    그리고 singular integral을 한국어로는 뭐라고 할수 있을가요?

    • singular integral은 특이적분이라고 합니다

      singularity를 특이성이라고 부르는 것을 보면 singular하다는 건 좀 특이하다는 거라고 표현하는 것 같습니다

+ Recent posts