Differences

This shows you the differences between two versions of the page.

--- 전산물리학:변분법 [2025/03/14 11:25] – admin
+++ 전산물리학:변분법 [2025/03/14 13:11] (current) – [예제] admin
@@ Line 1: / Line 1: @@
 ======전자기학에서======
+경계 $\Gamma$로 둘러싸인 공간 $\Omega$에서 푸아송 방정식
+$$\nabla^2 u = -f$$
+이 정의되어 있다. 디리클레 경계조건 $\left. u \right|_\Gamma = u_0$ 혹은 노이만 경계조건 $\left. \partial u / \partial n \right|_\Gamma = \nabla u \cdot \hat{n} = 0$이 만족된다고 가정한다. 이러한 $u$를 찾는 것은 다음 적분을 최소화하는 해 $u = \arg\min_{w} J\left[w\right]$를 찾는 것과 동치이다.
+$$J\left[w\right] \equiv \frac{1}{2} \int_\Omega \left| \nabla w \right|^2 - \int_\Omega f w ~dA.$$
+우리의 추정해인 $w$가 그러한 답 $u$에서 벗어난 정도가 $v$라고 하자. 경계조건을 만족해야 하므로 디리클레라면 $\left. v \right|_\Gamma = 0$일 것이다.
+적분을 해보면
+\begin{eqnarray*}
+J\left[ u+v \right] &=& \frac{1}{2} \int_\Omega \left| \nabla u \right|^2 dA - \int_\Omega fu ~dA\\
+&&+\int_\Omega \nabla u \cdot \nabla v ~dA - \int_\Omega fv ~dA\\
+&& + \frac{1}{2} \int_\Omega \left| \nabla v \right|^2 dA
+\end{eqnarray*}
+여기에서 세 번째와 네 번째 항을 따로 적어보자. 발산정리를 활용하면 다음처럼 고쳐 쓸 수 있다:
+\begin{eqnarray*}
+\delta J &=& \int_\Omega \nabla u \cdot \nabla v ~dA - \int_\Omega fv ~dA\\
+&=& \int_\Omega \nabla \cdot \left( v\nabla u \right) dA - \int_\Omega v \nabla^2 u ~dA - \int_\Omega fv ~dA\\
+&=& \int_\Gamma v \nabla u \cdot \hat{n} dl + \int_\Omega v \left(\nabla^2 u - f\right) ~dA.
+\end{eqnarray*}
+두 항 모두 0임을 알 수 있다. 따라서 다음 결론을 얻는다:
+\begin{eqnarray*}
+J\left[ u+v \right] &=& \frac{1}{2} \int_\Omega \left| \nabla u \right|^2 dA - \int_\Omega fu ~dA + \frac{1}{2} \int_\Omega \left| \nabla v \right|^2 dA\\
+&\ge& \frac{1}{2} \int_\Omega \left| \nabla u \right|^2 dA - \int_\Omega fu ~dA\\
+&=& J\left[ u \right].
+\end{eqnarray*}
+=====예제=====
+푸아송 방정식 $\nabla^2 V = -\rho_0 = const.$가 2차원 공간 $\Omega \equiv [-1,1] \times [-1,1]$에서 정의되어 있다. 경계조건은 $V(\pm 1, y) = V(x, \pm 1)=0$이다.
+시험 함수 $\Phi$를 아래의 기저로 나타내자:
+$$u_{mn} = \left(1-x^2 \right) \left(1-y^2 \right) \left( x^{2m} y^{2n} + x^{2n} y^{2m} \right).$$
+여기에서 $m,n = 0,1,2,\ldots$이다.
+차 근사로서 $m=n=0$인 하나의 기저만 사용해보자. 그러면 시험 해는 $\Phi = a_1 u_1$이며, 변분법을 활용하기 위해 이를 적분에 집어 넣으면
+\begin{eqnarray*}
+J &=& \frac{1}{2} \int_\Omega \left| \nabla \Phi \right|^2 dA - \int_\Omega \rho_0 \Phi ~dA\\
+&=& \frac{a_1^2}{2} \int_{-1}^1 \int_{-1}^1 \left[ 4x^2 \left(1-y^2\right)^2 + \left(1-x^2\right) 4y^2 \right] dx dy - \rho_0 a_1 \int_{-1}^1 \int_{-1}^1 \left(1-x^2\right) \left(1-y^2\right) dx dy\\
+&=& \frac{128}{45} a_1^2 - \frac{16}{9} \rho_0 a_1.
+\end{eqnarray*}
+$a_1$에 대한 이 이차함수의 최소를 구해보면 $a_1 = \frac{5}{16} \rho_0$임을 얻는다.
+$m=n=1$인 기저까지 활용해서
+$\Phi = a_1 \left(1-x^2 \right) \left(1-y^2 \right) +  a_2 \left(1-x^2 \right) \left(1-y^2 \right) \left( x^2 y^2 + x^2 y^2 \right)$로도 연습하고 해의 모양이 얼마나 바뀌었는지 확인해볼 수 있다.
+답은 $a_1 = \frac{1295}{4432} \rho_0$와 $a_2 = \frac{525}{8864} \rho_0$이다.
 ======양자역학에서======
@@ Line 92: / Line 131: @@
 ======함께 보기======
-[[수학:범함수]]
+  * [[수학:범함수]]
+  * [[전산물리학:모멘트법]]
+  * [[전산물리학:유한요소법]]
 ======참고문헌======
   * http://www.fisica.uniud.it/~giannozz/Corsi/MQ/LectureNotes/mq.pdf
+  * https://ocw.mit.edu/courses/16-920j-numerical-methods-for-partial-differential-equations-sma-5212-spring-2003/b08aec3e6f8cd87891d84d05b6d7a5eb_lecs13_14_notes.pdf
+  * Karl E. Lonngren, Sava V. Savov, and Randy J. Jost, //Fundamentals of Electromagnetics with MATLAB//, (SciTech Publishing, 2007).
+  * Matthew N.O. Sadiku, //Computational Electromagnetics with MATLAB//, (CRC Press, Boca Raton, FL, 2019).