Thin airfoil theory는 말 그대로 날개가 얇다고 가정을 한 상태에서 날개 주변의 유동에 대해 해석적인(analytic) 해를 구하는 것이다. 상당히 극단적인 가정이 많이 들어가기는 하지만 의외로 잘 맞아서(실험 데이터와 오차범위를 비교해놓은 것을 보면 정확하다는 말만 나온다.) 이전에 공학과 흑묘백묘론이라는 글을 쓰는 모티브가 되기도 했다. 물론 모티브는 내 엉망인 글 솜씨에 의해서 망해버렸지만.

근사를 하는 과정은 다음과 같다.

1. 점성항을 제거한다.
Navier-Stokes 공식을 Euler 공식으로 바꾸는 것이다. 이렇게 공식을 바꾸어주면 그나마 풀 수 있는 문제로 바뀌게 된다.[각주:1]

2. 회전하지 않는다고 가정한다.
속도 벡터장의 회전도(curl)를 항등적으로 0으로 취급한다는 뜻이다. 1번과 함께 이 조건이 만족되면 속도를 스칼라 함수의 물매(gradient)로 쓸 수 있다고 해서 potential flow라고 부른다. 컴퓨터 없이 유체역학을 공부하게 되면 이런 종류의 흐름만 배울 것이다.

3. 날개를 camber만 남긴다.
camber는 날개의 윗 단면과 날개의 아랫 단면의 정중앙을 지나는 곡선이다. 임의의 수직선을 생각했을 때 camber가 날개를 정확히 반으로 가른다고 생각해도 좋고, 윗 단면과 아랫 단면의 평균이라고 생각해도 좋다. 중요한 것은 날개에 이 녀석만 남겨서 두께를 0으로 만들어 버린다는 것이다.

4. x축에 vortex를 적절히 배치한다.[각주:2]
말 그대로 적절히 배치한다. 이 배치를 잘 조절해서 camber만 남긴 날개와 유체의 흐름이 평행하도록 만들어주는 것이다. 비행기가 날 때 공기가 날개를 뚫고 흐르지는 않는 것을 모사한다.

정확한 것은 책을 찾아보세요. 블로그는 언제까지나 블로그이다.

문제는 결국 이 방정식으로 줄어들게 된다. gamma는 x축 위의 vortex 분포, V는 무한지점에서 불어오는 속도(그래서 아래첨자가 무한대이다), alpha는 받음각(angle of attack), x는 좌표, z는 camber의 공식. 0에서부터 c까지 적분한다는 소리는 날개의 앞쪽 끝에서 뒤쪽 끝까지 적분한다는 의미이다.

\frac1{2\pi}\int_{0}^{c}\frac{\gamma(\xi)d\xi}{x-\xi}=V_\infty\left(\alpha-\frac{dz}{dx}\right)

Anderson의 Fundamentals of Aerodynamics 4th Ed.을 쓰고 있는데, 여기서는 날개를 완전히 평평한 판으로 근사했을 때(dz/dx=0)의 해를 다음과 같이 구해놓았다.

\gamma(\theta)=2\alpha V_\infty\frac{1+\cos(\theta)}{\sin(\theta)}\\x=\frac{c}2(1-\cos(\theta))

재미있는 것은 이 식이 경험법칙인 Kutta condition을 만족한다는 사실이다. Kutta 조건은 어떤 실험을 해 보아도 날개 끝 지점에서 유체의 움직임은 끝의 생김새와 평행하다는 것이다. 수식으로 바꾸면 gamma가 날개의 끝 지점에서는 0이라는 말이다.


http://galileo.phys.virginia.edu/classes/311/notes/aero/node4.html

아래 그림이 Kutta 조건을 만족할 때의 유동 모습이다. 일단 책에 나온 해의 그래프는 다음과 같이 그려진다.

c를 1로 잡았다. y축은 어차피 분포를 보여주기 위한 목적이 전부이니 무시.

자, 어떻게 다음 적분방정식(integral equation)의 해가 저렇게 깔끔한 성질을 보여줄까?

\frac1{2\pi}\int_{0}^{c}\frac{\gamma(\xi)d\xi}{x-\xi}=V_\infty\alpha

결론부터 말하자면, 당연히 깔끔할 수 밖에 없다. 일반적으로 적분방정식의 해는 여러개가 나오는데 그 중에서 우리가 원하는 조건에 맞는 녀석만 남도록 경계조건(boundary condition)을 적용했기 때문이다.[각주:3] 방정식 자체는 alpha와 V의 부호를 동시에 바꾸어도 동일하다는 점에 주목하길 바란다. 이 말은 위에서 구한 분포가 반대 방향에서 불어오고 있을 때에도 유효한 답이 된다는 말과 똑같은 소리이다.

\frac1{2\pi}\int_{0}^{c}\frac{\gamma(\xi)d\xi}{x-\xi}=(-V_\infty)(-\alpha)
의 해도\gamma(\theta)=2\alpha V_\infty\frac{1-\cos(\theta)}{\sin(\theta)}\\x=\frac{c}2(1-\cos(\theta))
가 된다.


다르게 표현한다면 다음 해도 사실은 유효하다는 것이다.부호 반대.

\gamma(\theta)=2\alpha V_\infty\frac{1-\cos(\theta)}{\sin(\theta)}\\x=\frac{c}2(1-\cos(\theta))

하늘색이 새로운 해의 분포

그래서 수업시간에 교수님께서 질문하신 "왜 답이 이렇게 깔끔한지 알기나 하냐?"에 대한 답은 "깔끔하도록 제한해주었으니까요"가 되겠다.



2012.11.14
좌우 대칭으로 바꾸어주면서 gamma의 부호도 바꾸어주어야 하는데 그 부분이 빠졌다.(시계방향 회전을 좌우대칭으로 바꾸면 반시계방향 회전이 된다.) 따라서 실제 위의 방정식을 만족하는 해는 다음 식이 된다.


\gamma(\theta)=2\alpha V_\infty\frac{\cos(\theta)-1}{\sin(\theta)}\\x=\frac{c}2(1-\cos(\theta))


양력은 음수가 되는 것을 알 수 있다. 전혀 쓰일 이유가 없는 해인 것이다.


  1. N-S 방정식을 완전히 풀어낸 해는 모두 합쳐도 손으로 꼽을 수 있을 정도밖에는 안된다. 그 중 하나가 원통 내부를 흘러가는 유체의 유동방정식인데, 이 녀석도 유체가 안정적으로(laminar) 흐를 때에만 유용한 녀석이라 좀 문제가 있다. 수도꼭지에서 물을 틀면 물이 매끈한 원기둥처럼 나올 때가 있고 이곳 저곳 울퉁불퉁한 흐름이 생기는 경우도 있는데 후자가 근처에서 가장 쉽게 관찰할 수 있는 불안정한(turbulent) 경우이다. 그나마 제대로 푼 식이라도 제한적으로만 의미가 있다는 뜻이다. CFD(Computational Fluid Dynamics)가 발전한 이유이기도 하고. [본문으로]
  2. 전자기학에 비유한다면, 적절한 전하밀도를 분포시키는 것이다. [본문으로]
  3. 양자역학의 산란 파트에서 Born approximation을 할 때에도 적분방정식을 푼다. 이때 나오는 독립적인 해가 두개인가 되는데, 그 중 하나는 경계조건으로 날려버린다. [본문으로]

댓글을 달아 주세요

  1.  댓글주소  수정/삭제  댓글쓰기

    비밀댓글입니다

    2010.05.05 16:40
    • Favicon of https://dexterstory.tistory.com BlogIcon 덱스터 2010.05.05 23:19 신고  댓글주소  수정/삭제

      원통좌표계를 사용하고 theta 방향의 속도 성분만 r에 대한 함수로 존재한다고 가정하면 됩니다.(theta에 대한 속도분포의 변화가 전혀 없다고 가정할 경우) 비압축성이니 속도 분포가 r에 비례하는 꼴 하나와 r에 반비례하는 꼴 하나가 나올 것이고(연속성에서 나옵니다) 해는 이 두 분포의 조합으로 나올 것 같네요. B.C은 no-slip으로 주면 됩니다. 일반적인 유체역학 책에서는 이런 방법으로 풀 겁니다.

1 ··· 9 10 11 12 13 14 15 16 17 ··· 23 

글 보관함

카운터

Total : 659,458 / Today : 45 / Yesterday : 123
get rsstistory!