Differences

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

--- 물리:소성_사건_plastic_event_의_전파_인자_propagator [2023/04/23 08:48] – minwoo
+++ 물리:소성_사건_plastic_event_의_전파_인자_propagator [2023/09/07 06:57] (current) – minwoo
@@ Line 1: / Line 1: @@
 ====== Oseen tensors ======
 ===== 레이놀즈 수 (Reynolds number) =====
-'비압축성'을 갖는 흐름(incompressible flow, $\nabla \cdot \boldsymbol{u}=0$)에 대한 나비에-스토크스 방정식(Navier-Stokes equation)은 다음과 같다.
+'비압축성'을 갖는 흐름(incompressible flow, $\nabla \cdot \mathbf{u}=0$)에 대한 나비에-스토크스 방정식(Navier-Stokes equation)은 다음과 같다.
 $$
-\rho \frac{\partial \boldsymbol{u}}{\partial t} + \rho \boldsymbol{u} \cdot \nabla \boldsymbol{u} = - \nabla p + \mu \nabla ^2 \boldsymbol{u} + f_{ext}
+\rho \frac{\partial \mathbf{u}}{\partial t} + \rho \mathbf{u} \cdot \nabla \mathbf{u} = - \nabla p + \mu \nabla ^2 \mathbf{u} + f_{ext}
 $$
-$p$는 압력이다. $\rho \boldsymbol{u} \cdot \nabla \boldsymbol{u}$의 항은 관성 항(inertial term)이며, $\mu \nabla ^2 \boldsymbol{u}$는 점성 항(viscous term)으로서
+$p$는 압력이다. $\rho \mathbf{u} \cdot \nabla \mathbf{u}$의 항은 관성 항(inertial term)이며, $\mu \nabla ^2 \mathbf{u}$는 점성 항(viscous term)으로서
 특성 속도의 크기가 $U_0$이며 그 속도에 대한 특성 길이의 크기를 $L$이라고 할 때, 각 항들은 다음과 같이 근사가 될 수 있다.
 $$
-\rho \boldsymbol{u} \cdot \nabla \boldsymbol{u} \approx \rho \frac{U_0 ^2}{L}, \ \  \mu \nabla ^2 \boldsymbol{u} \approx \frac{\mu U_0}{L^2}
+\rho \mathbf{u} \cdot \nabla \mathbf{u} \approx \rho \frac{U_0 ^2}{L}, \ \  \mu \nabla ^2 \mathbf{u} \approx \frac{\mu U_0}{L^2}
 $$
@@ Line 25: / Line 26: @@
 레이놀즈 수가 큰 값일 경우에는 점도가 낮은 경우이고, 작은 값인 경우에는 점도가 높은 경우이다.
-즉, $Re \ll 1 $으로서 '스토크스 흐름(Stokes flow)'을 가정한다면, $\rho \boldsymbol{u} \cdot \nabla \boldsymbol{u}$의 항은 근사적으로 무시할 수 있다.
+즉, $Re \ll 1 $으로서 '스토크스 흐름(Stokes flow)'을 가정한다면, $\rho \mathbf{u} \cdot \nabla \mathbf{u}$의 항은 근사적으로 무시할 수 있다.
-또한, 그 경우에 유속의 변화율인 $\rho \frac{\partial \boldsymbol{u}}{\partial t} \approx \rho \frac{U_0 ^2}{L}$의 항도 무시할 수 있으므로
+또한, 그 경우에 유속의 변화율인 $\rho \frac{\partial \mathbf{u}}{\partial t} \approx \rho \frac{U_0 ^2}{L}$의 항도 무시할 수 있으므로
 나비에-스토크스 방정식은 다음과 같이 표현된다.
 $$
--\nabla p(\boldsymbol{r}) + \mu \nabla ^2 \boldsymbol{u(r)} = -\boldsymbol{F}\delta(\boldsymbol(r)), \\
+-\nabla p(\mathbf{r}) + \mu \nabla ^2 \mathbf{u(r)} = -\mathbf{F}\delta(\mathbf(r)), \\
-\nabla \cdot \boldsymbol{u} = 0.
+\nabla \cdot \mathbf{u} = 0.
 $$
-여기서 $\boldsymbol{F}$는 점성이 있는 액체(viscous liquid)에 담겨져 있는 하나의 점 입자에 작용되는 힘이다.
+여기서 $\mathbf{F}$는 점성이 있는 액체(viscous liquid)에 담겨져 있는 하나의 점 입자에 작용되는 힘이다.
 $$\\$$
-$\boldsymbol{u(r)}$에서 $\boldsymbol{r}$이 무한하면 $0$의 값으로 사라져야 한다면 ($\lim_{r \to \infty} |\boldsymbol{u(r)}| = 0$), 해 $\boldsymbol{u(r)}$은 유일하게 결정되어야 한다.
+$\mathbf{u(r)}$에서 $\mathbf{r}$이 무한하면 $0$의 값으로 사라져야 한다면 ($\lim_{r \to \infty} |\mathbf{u(r)}| = 0$), 해 $\mathbf{u(r)}$은 유일하게 결정되어야 한다.
 $$\\$$
-===== $p(\boldsymbol{r}), \ \boldsymbol{u(r)}$ =====
+===== $p(\mathbf{r}), \ \mathbf{u(r)}$ =====
-위에서 $f_{ext}=\boldsymbol{F}\delta(\boldsymbol{r})$라고 설정한 것에 따라, 해의 형태를 다음과 같이 가정하자.
+위에서 $f_{ext}=\mathbf{F}\delta(\mathbf{r})$라고 설정한 것에 따라, 해의 형태를 다음과 같이 가정하자.
-$$ p(\boldsymbol{r}) = \frac{\boldsymbol{F} \cdot \boldsymbol{P(r)}}{8\pi \mu}, \ \ \boldsymbol{u(r)} = \frac{\mathbb{G}(\boldsymbol{r}) \cdot \boldsymbol{F}}{8\pi \mu} $$
+$$ p(\mathbf{r}) = \frac{\mathbf{F} \cdot \mathbf{P(r)}}{8\pi \mu}, \ \ \mathbf{u(r)} = \frac{\mathbb{G}(\mathbf{r}) \cdot \mathbf{F}}{8\pi \mu} $$
-이때 $\boldsymbol{P(r)}$은 어떤 벡터장이며, $\mathbb{G}(\boldsymbol{r})$은 힘 $\boldsymbol{F}$(point force)에 적용되는 텐서이다.
+이때 $\mathbf{P(r)}$은 어떤 벡터장이며, $\mathbb{G}(\mathbf{r})$은 힘 $\mathbf{F}$(point force)에 적용되는 텐서이다.
 $$\\$$
@@ Line 55: / Line 56: @@
 $$
-p(\boldsymbol{r}) = \frac{P_j F_j}{8\pi \mu}, \ \ u_i(\boldsymbol{r}) = \frac{\mathbb{G}_{ij} F_j}{8\pi \mu}
+p(\mathbf{r}) = \frac{P_j F_j}{8\pi \mu}, \ \ u_i(\mathbf{r}) = \frac{\mathbb{G}_{ij} F_j}{8\pi \mu}
 $$
@@ Line 64: / Line 65: @@
 $$
-\mathcal{F}\{f \} = \hat{f} (\boldsymbol{k}) = \int_{\mathbb{R}^3} d\boldsymbol{r} f(\boldsymbol{r}) e^{-i\boldsymbol{k}\cdot \boldsymbol{r}}, \\
+\mathcal{F}\{f \} = \hat{f} (\mathbf{k}) = \int_{\mathbb{R}^3} d\mathbf{r} f(\mathbf{r}) e^{-i\mathbf{k}\cdot \mathbf{r}}, \\
-\mathcal{F}^{-1}\{\hat{f} \} = f(\boldsymbol{r}) = \frac{1}{(2 \pi)^3} \int_{\mathbb{R}^3} d\boldsymbol{k} \hat{f}(\boldsymbol{k}) e^{i\boldsymbol{k}\cdot \boldsymbol{r}}
+\mathcal{F}^{-1}\{\hat{f} \} = f(\mathbf{r}) = \frac{1}{(2 \pi)^3} \int_{\mathbb{R}^3} d\mathbf{k} \hat{f}(\mathbf{k}) e^{i\mathbf{k}\cdot \mathbf{r}}
 $$
 위의 변환 식에 따라서, 도함수의 푸리에 변환에 대한 간단한 결과를 다음과 같이 얻을 수 있다.
 $$
-\frac{d}{d\boldsymbol{r}} f(\boldsymbol{r}) = i\boldsymbol{k} \frac{1}{(2 \pi)^3} \int_{\mathbb{R}^3} d\boldsymbol{k} \hat{f}(\boldsymbol{k}) e^{i\boldsymbol{k}\cdot \boldsymbol{r}} = i\boldsymbol{k} \mathcal{F}^{-1}\{\hat{f} \} \\
+\frac{d}{d\mathbf{r}} f(\mathbf{r}) = i\mathbf{k} \frac{1}{(2 \pi)^3} \int_{\mathbb{R}^3} d\mathbf{k} \hat{f}(\mathbf{k}) e^{i\mathbf{k}\cdot \mathbf{r}} = i\mathbf{k} \mathcal{F}^{-1}\{\hat{f} \} \\
-\to \mathcal{F}\{\frac{d}{d\boldsymbol{r}} f(\boldsymbol{r}) \} = i \boldsymbol{k} \mathcal{F}\{f\}
+\to \mathcal{F}\{\frac{d}{d\mathbf{r}} f(\mathbf{r}) \} = i \mathbf{k} \mathcal{F}\{f\}
 $$
 $$\\$$
-따라서, 나비에-스토크스 방정식인 $ -\nabla p(\boldsymbol{r}) + \mu \nabla ^2 \boldsymbol{u(r)} = -\boldsymbol{F}\delta(\boldsymbol{r}) $을 다음과 같이 고칠 수 있다.
+따라서, 나비에-스토크스 방정식인 $ -\nabla p(\mathbf{r}) + \mu \nabla ^2 \mathbf{u(r)} = -\mathbf{F}\delta(\mathbf{r}) $을 다음과 같이 고칠 수 있다.
 $$
--i \boldsymbol{k} \hat{p} - \mu k^2 \hat{u} = -\boldsymbol{F}.
+-i \mathbf{k} \hat{p} - \mu k^2 \hat{u} = -\mathbf{F}.
 $$
-위의 식에, 가정했던 $p(\boldsymbol{r})$과 $u(\boldsymbol{r})$의 해의 형태를 대입하고 성분 형태로 고치면 다음과 같다.
+위의 식에, 가정했던 $p(\mathbf{r})$과 $u(\mathbf{r})$의 해의 형태를 대입하고 성분 형태로 고치면 다음과 같다.
 $$
@@ Line 95: / Line 96: @@
 $$
-여기에, 비압축성(incompressibility) 조건인 $\nabla \cdot \boldsymbol{u} = 0$을 다음과 같이 푸리에 공간에 대해서 적용할 수 있다.
+여기에, 비압축성(incompressibility) 조건인 $\nabla \cdot \mathbf{u} = 0$을 다음과 같이 푸리에 공간에 대해서 적용할 수 있다.
 $$ k_i \mathbb{G}_{ij}=0. $$
-(위의 식은 $\boldsymbol{u(r)} = \frac{\mathbb{G}(\boldsymbol{r}) \cdot \boldsymbol{F}}{8\pi \mu}$의 식을 참고하고 발산(divergence)의 공식을 참고하면 확인할 수 있다.)
+(위의 식은 $\mathbf{u(r)} = \frac{\mathbb{G}(\mathbf{r}) \cdot \mathbf{F}}{8\pi \mu}$의 식을 참고하고 발산(divergence)의 공식을 참고하면 확인할 수 있다.)
 $$ \\ $$
@@ Line 108: / Line 109: @@
 \to \frac{\hat{P}_j}{8\pi \mu} = -\frac{ik_j}{k^2}.
 $$
+이를 앞서 얻은 $-ik_i \frac{\hat{P}_j}{8\pi \mu} - k^2 \frac{\hat{\mathbb{G}_{ij}}}{8\pi} = -\delta_{ij}$ 의 식에 대입하여 다음과 같이 $\hat{\mathbb{G}}_{ij}$의 식을 얻는다.
+$$
+-\frac{k_ik_j}{k^2} - k^2 \frac{\hat{\mathbb{G}}_{ij}}{8\pi} =-\delta_{ij}\\
+\to \frac{\hat{\mathbb{G}}_{ij}}{8\pi}=\frac{\delta_{ij}}{k^2}
+-\frac{k_ik_j}{k^4}
+$$
+$$ \\ $$
+이때, 위의 결과에 대해서 다음을 'Oseen tensor'라고 부른다.
+$$
+\hat{\mathbf{O}}(\mathbf{k})=\frac{\hat{\mathbb{G}}}{8\mu\pi}=\frac{1}{\mu k^2}\left(\mathbf{I}-\frac{\mathbf{k}\mathbf{k}}{k^2}\right)
+$$
+편의상 $\mathbf{k}$를 $\mathbf{q}$로 표기하면, 다음과 같이 쓰인다.
+$$
+\hat{\mathbf{O}}(\mathbf{q}) =\frac{1}{\mu q^2}\left(\mathbf{I}-\frac{\mathbf{q}\mathbf{q}}{q^2}\right)
+$$
@@ Line 123: / Line 148: @@
 $$
-그의 r에 대한 기울기(gradient)를 계산하는 것과 같다.
+그의 \mathbf{r}에 대한 기울기(gradient)를 계산하는 것과 같다.
 $$ \\ $$