2011년 7월 10일 일요일

원통 좌표계(圓筒座標系, circular cylindrical coordinate system)

[경고] 아래 글을 읽지 않고 "원통 좌표계"를 보면 바보로 느껴질 수 있습니다.
1. 텐서 미적분학
2. 직교 좌표계 텐서 미적분학


[그림 1] 원통 좌표계의 표현(출처: wikipedia.org)

직교 좌표계(orthogonal coordinate system)에서 데카르트 좌표계(Cartesian coordinate system) 다음으로 많이 쓰이는 좌표계가 [그림 1]의 원통 좌표계(圓筒座標系, circular cylindrical coordinate system)이다. 원통 좌표계 $(\rho, \phi, z)$가 2차원 $(\rho, \phi)$로만 표현되면 극좌표계(polar coordinate system)라고도 한다. [그림 2]는 원통 좌표계 좌표값의 변화 모습을 보여준다.
[그림 2] 원통 좌표계 좌표값의 변화 모습(출처: wikipedia.org)

[그림 1]과 [그림 2]의 좌표계 구성으로 인해 데카르트 좌표계 $X$에서 원통 좌표계 $U$로 가는 좌표변환(coordinate transform)은 아래와 같다.

                  (1)

식 (1)은 [그림 3]과 삼각 함수(trigonometric function)를 이용해 쉽게 유도할 수 있다.

[그림 3] 데카르트 좌표계와 원통 좌표계의 상호 변환(출처: wikipedia.org)

식 (2)의 계량 텐서(metric tensor)와 관련된 척도 인자(尺度因子, scale factor) $h_i$는 식 (3)처럼 계산된다.

                         (2)

                         (3)

식 (1)을 식 (3)에 대입해 계산하면 원통 좌표계의 척도인자를 얻을 수 있다.

                         (4)

식 (4)를 직교 좌표계에 대한 벡터 연산자 공식에 대입해보자.

                       (5)

        (6)

                       (7)

                       (8)

그러면, 원통 좌표계에 대한 벡터 연산자를 모두 정의할 수 있다.

                       (9)

                       (10)

                       (11)

                       (12)

다음으로 원통 좌표계를 구성하는 단위 벡터(unit vector)를 데카르트 좌표계 관점으로 계산해보자. 벡터 $\bar r$이 $\rho, \phi, z$방향으로 바뀔 때 생기는 벡터 방향이 $\rho, \phi, z$방향 단위 벡터가 된다. 이를 수식으로 표현하면 다음과 같다.

                       (13)

예를 들면 $\phi$방향 단위 벡터는 $\phi = 0$에서는 $y$방향으로, $\phi = 90^\circ$는 $-x$방향으로 형성된다.
식 (13)이 표현하는 것처럼 단위 벡터 $\rho, \phi$는 방위각(方位角, azimuth) $\phi$에 따라 변한다. 단위 벡터 $\rho, \phi$를 미분할 때는 방위각의 변화율을 다음처럼 반드시 고려해야 한다.

                       (14)

식 (13)을 이용하면 데카르트 좌표계를 원통 좌표계로 바꾸어주는 공식을 행렬(matrix) 형태로 만들 수 있다.

                       (15)

식 (15)에서 재미있는 것은 행렬 ${\bf T}_{cr}$과 ${\bf T}_{rc}$가 서로 전치 행렬(transpose) 관계라는 것이다.

식 (9)에서 (12)에 제시한 원통 좌표계 벡터 연산자는 무식한 방법으로 구할 수도 있다. 데카르트 좌표계(Cartesian coordinate system)에서 유도한 결과와 완전 미분(exact differential)을 결합하면 식 (9)부터 (12)까지 있는 원통 좌표계 벡터 연산자를 동일하게 유도할 수 있다.
예를 들어 식 (16)의 데카르트 좌표계 구배 연산자(gradient operator)는 원통 좌표계에서 식 (17)처럼 유도된다.

                         (16)

                        (17)

식 (17) 유도에서 식 (15)에 증명한 벡터 관계식을 사용하였다.
식 (18)의 데카르트 좌표계 발산 연산자(divergence operator)는 원통 좌표계에서 식 (19)처럼 유도된다.

                        (18)

   

                  (19)

무식하게 계산한 식 (17), (19)는 신기하게도 텐서 미적분학(tensor calculus)을 기반으로 계산한 식 (9), (10)과 일치한다. 사실 텐서 미적분학도 위에서 연쇄법칙(chain rule)을 이용해 유도한 방식을 수학적으로 우아하게 만든 것일 뿐이다.

[다음 읽을거리]
1. 구 좌표계

댓글 16개 :

  1. :-) 도움이 되었다니 저도 좋네요!

    답글삭제
  2. 잘보고가네요!! 내용도 도움이되었고! 위쪽에 조금은 느리게 살자라는 글을 보면서 또 힘을 얻고 가네요ㅎ 문제푸는 방법만 외우는것 대신에 내용을 이해하려고 하다보니 진도가 더디고 그래서 마음도 흔들리고 있었는데ㅎ 속도보단 이렇게 꾸준히 해야되는 거겠죠?

    답글삭제
  3. 감사합니다. 윗글은 저의 각오이기도 합니다. 공부 계속 할거면 모든 부분들에 의심이 없도록 철저하게 이해해야겠지요.

    답글삭제
  4. 식 (19)에서 중간에 (로 분의 에프 로) 오타 같네요

    답글삭제
  5. 댓글 감사합니다. 제가 식 (19)를 다시 확인해봤는데 틀린 부분은 없습니다. Fρ/ρ가 나오는 부분은 위쪽에 있는 ∂/∂φ(·)를 적용할 때 생기는 부분입니다.

    답글삭제
  6. 아이고 도움되는 유익한 블로그네요 ^^감사합니다

    답글삭제
  7. 질문좀 드릴게요..

    cylindrical a_pi vector를 rectangular 형태로 변환시키면 -sin(pi)a_x + cos(pi)a_y 가 맞는지 해서요...

    a_pi vector를 rectangular vector인 a_x와 외적을 해야하는데 바로 안되서 말입니다.

    답글삭제
    답글
    1. 예, 맞습니다. 유도는 식 (13)으로 하면 됩니다.
      벡터 r이 φ방향으로 바뀔 때 생기는 벡터 방향이 φ방향 단위벡터입니다. 예를 들면 φ = 0에서는 y방향으로, φ = 90도는 -x방향으로 φ방향 단위벡터가 형성됩니다.

      삭제
    2. 감사합니다. 대전에서 전파 관련 일 하시나봅니다? 저도 대전에서 전파공학전공중인 학생이거든요

      학생이신지 혹시?? 여기 좋은 자료들이 너무나도 많네요

      삭제
  8. 질문이 있습니다
    17에 a로우/ax가 cos(pi)로 되는 이유가 궁금합니다! api/ax도 마찬가지로요..

    답글삭제
    답글
    1. 식 (1)을 직접 미분하시면 나옵니다.

      삭제
  9. 저 시간이 많이 지났지만 질문드려도 될까요... 식(18)과 식(19)사이에 d(Fx)/dx 랑 d(Fy)/dy 구하는 부분에서 d로/dx 는 cos(pi)가 맞는데 dpi/dx 는 왜 sin(pi)/로 인가요?? 텐서미적분은 잘 모르고 cartesian 좌표계에서 cylindrical 좌표계로 divergent를 손으로 유도해보고 있는데 항상 저기서 막힙니다... 다른 책들 찾아봐도 그냥 그렇다고만 나오고 이유를 못찼겠어요

    답글삭제
  10. 위에 17년 9월 29일에 질문글 썼던 사람인데요 궁금한걸 보충설명 해야될 것 같아서요.
    식(13)에 로와 파이방향 벡터를 X와 Y방향 벡터로 표현 하잖아요? 로에서 cos(pi)가 d로/dx, 파이에서 -sin(pi)가 d파이/dx (d를 흔히말하는 '라운드 디'대신 쓰겠습니다)일텐데 식(19)에서 d파이/d가 -sin(pi)/로 이렇게 나와서 이 부분이 궁금합니다..

    답글삭제
    답글
    1. $\phi = \tan^{-1} (y/x)$를 $x$ 혹은 $y$에 대해 미분한 결과가 식 (19)에 쓰였습니다. 간단하게 생각하세요, 익명님. ^^

      삭제

욕설이나 스팸글은 삭제될 수 있습니다. [전파거북이]는 선플운동의 아름다운 인터넷을 지지합니다.