mc-levha

Monte Carlo dilim taşınımı

Monte Carlo

ÖLÇÜM SÜRÜYOR

TOHUM00003039

Dilim kalınlığı L = 3.0 cm

Toplam tesir kesiti Σ_t = 1.00 /cmSaçılma oranı c = Σ_s/Σ_t = 0.70Geliş kosinüsü μ₀ = 1.00

Geçmiş sayısı N = 20 000

⚅ AYNI TOHUM = AYNI İZLER

ÇALIŞIYOR

GEÇİŞ T

———

YANSIMA R

———

YUTULMA A

———

TARİHÇE

ƒ μ′(μ, ξ) — CANLI v1kodu düzenle

scatter(mu, xi, xi2) → μ′● CANLI v1

function scatter(mu, xi, xi2) {

}

Çarpışmada yeni yön kosinüsünü sen yazarsın. İzotropik: 2*xi−1; ileri-yönlü: mu ağırlıklı. Albedo ve nüfuzu değiştir.

Lewis & Miller · Lux & Koblinger · kaynakarka uç: TS · analog MC · track-length

Uçuş örneklemesi

s=-\dfrac{\ln\xi}{\Sigma_t}

Bir sonraki çarpışmaya kadarki mesafe, üstel serbest-yol dağılımından çekilir.

Soğurma olasılığı

P_{abs}=1-c,\quad c=\dfrac{\Sigma_s}{\Sigma_t}

Her çarpışmada nötron 1−c olasılıkla soğurulur, c olasılıkla saçılır.

Çarpışmasız iletim

T_0=e^{-\Sigma_t L/\mu_0}

Saf soğurucuda (c=0) analog yürüyüş tam olarak bu Beer–Lambert sonucunu verir.

Space çalıştır/duraklat · . adım · R sıfırla (tut) · 1/2/3 seviye · A onay

k → 1.000Kritikliğe yaklaşılıyor… / Approaching criticality…