광학 두께

편집 역사 토론
계산물리 물리 마지막 수정: 2026-08-02 05:33:48

1. 개요[편집]

매질의 두께를 미터로 재면 아무 의미가 없다. 광자가 몇 번이나 부딪히는지로 재야 한다.

광학 두께(optical depth / optical thickness, τ\tau)는 복사가 매질을 통과할 때 겪는 소산의 누적량을 무차원으로 잰 양이다. 소산계수 χν=κν+σν\chi_\nu = \kappa_\nu + \sigma_\nu (흡수 + 산란)에 대해 경로 ss 를 따라

τν=∫0sχν ρ ds′\tau_\nu = \int_0^s \chi_\nu\,\rho\,ds'

로 정의하고, 방출·산란이 없다면 비어-람베르트 법칙 Iν=Iν,0e−τνI_\nu = I_{\nu,0}e^{-\tau_\nu} 로 세기가 줄어든다. 이 문서의 요지는 정의 자체가 아니라 τ\tau 의 크기에 따라 지배 물리가 완전히 다른 방정식으로 바뀐다는 것, 그리고 그 두 극한에서 서로 다른 방식으로 평균한 불투명도를 써야 한다는 것이다.

복사 전달 방정식이 방정식의 형태와 이산화를, 흑체복사가 방출항 Bν(T)B_\nu(T) 의 크기를 다룬다면, 이 문서는 τ\tau 라는 하나의 무차원수가 해석 전략을 어떻게 갈라놓는지에 한정한다.

2. 얇음과 두꺼움, 그리고 광구[편집]

τ≪1\tau \ll 1 이면 광학적으로 얇다(optically thin). 매질을 지나는 광자의 대다수가 한 번도 상호작용하지 않고 빠져나가고, 따라서 우리는 매질 내부를 통째로 들여다본다. 방출된 복사가 재흡수 없이 탈출하므로 광학적으로 얇은 냉각이 성립하고, 총 방출량이 매질 전체 부피에 대한 단순 적분으로 계산된다. 성운의 방출선, 희박한 화염, 초기 단계의 플라스마 냉각이 여기 해당한다.

τ≫1\tau \gg 1 이면 광학적으로 두껍다(optically thick). 광자가 나오기 전에 수없이 흡수·재방출되어 국소 열평형에 가까워지고, 우리가 보는 것은 표면뿐이다. 이때 복사장은 국소적으로 등방에 가까워지고 Iν→Bν(T)I_\nu \to B_\nu(T) 로 수렴한다. 표면이 어디냐를 정하는 관례가 τ≈1\tau \approx 1 인 면이며, 항성에서 이것을 광구(photosphere)라 부른다.1 태양이 뚜렷한 원반으로 보이는 이유가 이것이다 — 물리적 경계면이 있어서가 아니라 τ=1\tau=1 을 넘는 순간 더 이상 안 보이기 때문이다.

산란이 섞이면 척도가 하나 더 필요하다. 단일 산란 알베도

ω=σsσa+σs\omega = \frac{\sigma_s}{\sigma_a + \sigma_s}

는 상호작용 한 번당 광자가 살아남을 확률이다. ω→0\omega \to 0 은 순수 흡수(빛이 사라진다), ω→1\omega \to 1 은 보존 산란(방향만 바뀐다)이다. 우유·구름·대리석·피부는 전부 ω\omega 가 0.99를 넘는 고알베도 매질이고, 그래서 광자가 수백 번 튕기며 특유의 부드러운 빛이 나온다.

3. 두꺼운 극한 — 확산 근사와 로세란트 평균[편집]

τ≫1\tau \gg 1 에서 복사 전달은 확산 방정식으로 환원된다. 국소적으로 Iν≈Bν+I_\nu \approx B_\nu + 작은 이방성 보정을 넣고 1차 모멘트를 취하면 복사 열유속이 온도 기울기에 비례하는 꼴로 나온다.

F=−16σT33 κR ρ ∇TF = -\frac{16\sigma T^3}{3\,\kappa_R\,\rho}\,\nabla T

괄호 앞이 곧 복사 열전도도 kradk_{\mathrm{rad}} 다. 즉 아주 불투명한 매질 안에서 복사는 열전도와 수학적으로 구별되지 않으며, 그냥 전도 솔버에 kk 를 더해 주면 된다. 유리 용융로·단열재·항성 내부 계산이 이 근사 위에 서 있다.

여기서 결정적인 것이 κR\kappa_R 의 정의다.

1κR=∫0∞1κν∂Bν∂T dν∫0∞∂Bν∂T dν\frac{1}{\kappa_R} = \frac{\displaystyle\int_0^\infty \frac{1}{\kappa_\nu}\frac{\partial B_\nu}{\partial T}\,d\nu}{\displaystyle\int_0^\infty \frac{\partial B_\nu}{\partial T}\,d\nu}

κ\kappa 가 아니라 1/κ1/\kappa 를 평균한다. 이유는 물리적으로 명백하다. 두꺼운 매질에서 에너지를 실어 나르는 것은 가장 잘 통과하는 주파수이고, 저항이 병렬로 붙은 회로처럼 가장 투명한 창이 흐름을 지배한다. 반면 얇은 극한에서는 총 방출량이 중요하므로 플랑크 평균

κP=∫0∞κνBν dν∫0∞Bν dν\kappa_P = \frac{\int_0^\infty \kappa_\nu B_\nu\,d\nu}{\int_0^\infty B_\nu\,d\nu}

즉 κ\kappa 자체의 산술평균(플랑크 가중)을 쓴다. 두 평균은 조화평균과 산술평균의 관계라 항상 κR≤κP\kappa_R \le \kappa_P 이고, 흡수선이 촘촘한 실제 기체에서는 한 자릿수 이상 벌어질 수 있다. 두꺼운 문제에 플랑크 평균을 넣으면 매질이 실제보다 훨씬 불투명해져 복사 열유속을 수십 % 이상 과소평가한다. 이 실수가 아직도 논문 리뷰에서 잡힌다.2

4. 수치적 함의 — 격자와 반복[편집]

격자 설계 규칙부터. 복사 문제의 셀 크기 기준은 길이가 아니라 셀당 광학 두께 Δτ\Delta\tau 다. Δτ≫1\Delta\tau \gg 1 인 셀은 셀 내부의 지수 감쇠를 전혀 분해하지 못해 1차 상류 도식이 통째로 틀리고, Δτ≪1\Delta\tau \ll 1 인 셀만 잔뜩 만들면 계산을 낭비한다. 경험칙은 관심 영역에서 Δτ∼1\Delta\tau \sim 1 이하, 경계층 근처에서는 0.1까지 줄이는 것이다.

반복 수렴은 더 고약하다. 산란항이 있는 복사 전달 방정식을 푸는 고전적 방법인 소스 반복(source iteration)은 물리적으로 “산란 nn 번 겪은 광자까지 누적”하는 것이라, 반복당 오차 감소율이 대략 ω\omega 다. 푸리에 해석으로 얻는 스펙트럴 반경은 매질이 두꺼워질수록 위쪽 극한

ρSI⟶ω(τ≫1)\rho_{\mathrm{SI}} \longrightarrow \omega \qquad (\tau \gg 1)

에 붙는다. 가장 느린 것은 파장이 계 크기만 한 저주파 오차 성분이다. ω=0.999\omega = 0.999 인 구름 셀에서는 오차가 반복당 0.1%씩만 줄어 수천 번을 돌려야 한다. 표준 처방이 확산 합성 가속(DSA, diffusion synthetic acceleration)이다 — 수송 스윕 한 번마다 저차 확산 문제를 풀어 저주파 오차를 통째로 걷어낸다. 잘 이산화된 DSA는 ω\omega 와 무관하게 반복당 감소율을 0.2 근처로 눌러 준다. 다만 저차 방정식과 고차 방정식의 이산화가 어긋나면 오히려 발산하는 것으로 악명 높다.

방법 선택도 τ\tau 로 갈린다. τ≫1\tau \gg 1 이면 P1/확산 근사가 싸고 정확하다. 중간 영역(τ∼1\tau \sim 1)이 가장 어려워서 이산종좌표법(SN_N)이나 유한체적 복사법을 쓰는데, 각을 유한개 방향으로 자른 대가로 레이 효과(ray effect)라는 줄무늬 아티팩트가 생긴다. τ≪1\tau \ll 1 이거나 기하가 복잡하면 몬테카를로 방법 기반 광선 추적이 사실상 유일한 실용해다. 광자를 하나씩 날리고 −ln⁡ξ-\ln \xi 로 다음 상호작용까지의 τ\tau 를 뽑는 방식이라 기하 복잡도에 거의 무관하지만, 통계 잡음이 N−1/2N^{-1/2} 로만 줄어든다.

5. 다른 동네의 같은 개념[편집]

τ\tau 는 사실 평균자유행로로 잰 거리다. ℓ=1/(χρ)\ell = 1/(\chi\rho) 로 두면 τ=s/ℓ\tau = s/\ell 이니, “광자가 평균 몇 번 부딪히는가”가 곧 광학 두께다. 이 관점이 여러 분야를 하나로 묶는다.

  • 중성자 수송: 같은 양을 광학거리(optical distance) 또는 mean free path 수로 부르고, Σt\Sigma_t 는 거시 전단면적이다. DSA도 원래 이쪽에서 나왔다.
  • 의료 영상: CT의 선감쇠계수 μ\mu 가 곧 χρ\chi\rho 이고, 투영값 ln⁡(I0/I)\ln(I_0/I) 가 정확히 τ\tau 다. 라돈 변환이 재구성하는 것이 이 적분이다.
  • 렌더링: 참여 매질(안개·연기·피부)의 볼륨 렌더링 적분에서 투과율 T=e−τT = e^{-\tau} 가 곧 가중치다. 레이 트레이싱과 물리 기반 렌더링의 불균질 매질 샘플링(델타 트래킹)은 최대 소산계수로 상한을 잡고 가상 충돌을 섞어 τ\tau 적분을 편향 없이 처리하는 기법이다.
  • 대기 CFD·기후 모델: 대역별 τ\tau 를 미리 계산해 두는 상관 kk-분포 방식이 표준이며, 여기서도 두꺼운 대역과 얇은 대역에 다른 근사를 쓴다.

구조적으로 가장 닮은 무차원수는 크누센수 Kn=ℓ/L\mathrm{Kn} = \ell/L 이다. 실제로 τ≈1/Kn\tau \approx 1/\mathrm{Kn} 이며, ”Kn≪1\mathrm{Kn}\ll1 이면 연속체(확산) / Kn≫1\mathrm{Kn}\gg1 이면 자유분자”라는 판정이 ”τ≫1\tau\gg1 이면 확산 근사 / τ≪1\tau\ll1 이면 자유 탈출”과 정확히 같은 논리다. 두 분야가 서로 모르는 채 같은 극한 이론을 각자 만들어 쓴 셈이다.34

6. 관련 문서[편집]

7. Footnotes[편집]

  1. 엄밀하게는 τ=2/3\tau = 2/3 이 표준이다. 에딩턴 근사에서 관측되는 유효온도가 τ=2/3\tau=2/3 지점의 온도와 같아진다는 에딩턴-바비에 관계에서 나온 값인데, 실무에서는 그냥 “1 근처”로 말한다. 2/3냐 1이냐로 싸우는 건 항성대기 하는 사람들끼리만 하자. ↩

  2. 외우는 요령: 로세란트는 “빠져나가는 놈”을 센다 — 그래서 가장 투명한 창이 이기고, 조화평균이다. 플랑크는 “방출하는 놈”을 센다 — 그래서 가장 불투명한 선이 이기고, 산술평균이다. 헷갈리면 흡수선 하나가 무한히 불투명한 극단을 상상하면 된다. 로세란트는 꿈쩍도 안 하고 플랑크는 무한대로 간다. ↩

  3. 그래서 학회에서 “우리는 광학적으로 두꺼운 극한을 확산으로 근사했습니다”와 “우리는 연속체 극한에서 나비에-스토크스를 씁니다”가 같은 말이라는 걸 서로 30분쯤 지나서야 알아차리는 광경이 벌어진다. ↩

  4. 참고로 τ\tau 는 시간이 아니다. 기호가 겹쳐서 응력·이완시간·전단응력과 한 논문 안에서 충돌하는 일이 잦은데, 복사 전달 논문에서 τ\tau 는 거의 항상 광학 두께다. 아닌 경우에는 저자가 반드시 각주를 단다. 안 달았으면 그 논문을 의심하는 게 맞다. ↩