放射线动画 — 熱放射的可视化
什么是辐射视角系数
🎓
简单来说，你可以把它想象成“看见”的程度。想象一下，你站在一个房间里，面前有一扇窗户。你能“看见”多少天空，就代表有多少来自天空的光线能到达你的眼睛。视角系数 $F_{12}$ 就是描述从表面1发出的热辐射，有多少比例能“直接看到”并到达表面2。它是一个介于0和1之间的数，只和两个面的形状、大小以及它们怎么摆放有关。在实际工程中，比如计算汽车发动机舱里某个零件的辐射散热，第一步就是要算这个系数。
🙋
诶，真的吗？那它和材料本身有关系吗？比如一个是亮闪闪的金属，一个是黑色的橡胶。
🎓
这是个好问题！视角系数本身和材料颜色、粗糙度这些都没关系，它纯粹是个“几何”问题。就像刚才说的，只取决于“能不能看见”。材料本身的属性，比如你提到的发射率 $\varepsilon$，是另一个参数，它决定了表面“愿意”发射或吸收多少辐射能。在我们的模拟器里，你可以分别调整 $\varepsilon_1$ 和 $\varepsilon_2$ 的滑块，你会发现它们不影响 $F_{12}$ 的计算，但会直接影响最终计算出的辐射换热量。你可以马上试试看！
🙋
原来是这样！那模拟器里说的“互易定理” $A_1F_{12}= A_2F_{21}$ 又是什么？感觉好神奇。
🎓
这个定理是视角系数计算的核心，它保证了能量交换的对称性。简单理解就是：从面1“看”面2，和从面2“看”面1，虽然“看见”的比例（$F_{12}$ 和 $F_{21}$）可能不同，但乘以各自的面积后，这个交换量是相等的。比如，一个巨大的天花板和一个小小的桌面，天花板“看”桌面的比例 $F_{12}$ 很小，但桌面“看”天花板的 $F_{21}$ 几乎接近1。你可以在工具里选择“平行矩形板”，然后只改变边长a和b，观察 $F_{12}$ 和 $F_{21}$ 如何变化，但下面验证互易定理的数值始终是相等的，这就是几何的魔力！
物理模型与关键公式
辐射视角系数 $F_{12}$ 的基本定义，来源于对两个表面之间辐射能交换的几何积分。它计算从表面1的微元 $dA_1$ 发出并到达表面2的辐射能比例。
$$F_{12}= \frac{1}{A_1}\int_{A_1}\int_{A_2}\frac{\cos \theta_1 \cos \theta_2}{\pi r^2}dA_2 dA_1$$
其中，$A_1$ 和 $A_2$ 是表面面积，$r$ 是两个微元之间的距离，$\theta_1$ 和 $\theta_2$ 分别是连线与两个表面法线的夹角。这个双重积分对于复杂形状很难直接计算，因此本工具内置了常见几何（如平行板、同轴圆盘）的解析解公式。
基于视角系数，两个漫灰表面之间的净辐射换热量由以下公式给出：
$$Q_{12}= A_1 F_{12}\sigma \varepsilon_{\text{eff}}(T_1^4 - T_2^4)$$
其中，$\sigma$ 是斯蒂芬-玻尔兹曼常数，$T$ 是绝对温度，$\varepsilon_{\text{eff}}$ 是考虑两个表面发射率 $\varepsilon_1$ 和 $\varepsilon_2$ 后的有效发射率。这个公式是工程热设计中计算辐射散热的核心。
现实世界中的应用
电子设备散热：在手机或服务器芯片的散热设计中，需要计算芯片表面（面1）到散热鳍片或外壳（面2）的辐射传热。通过计算视角系数，工程师可以优化布局，判断是否需要增加辐射涂层（改变发射率 $\varepsilon$）来强化散热。
航空航天热控：卫星在太空中主要依靠辐射与外部环境换热。计算卫星太阳能板、仪器舱等不同表面之间的视角系数至关重要，用于精确控制卫星各部分的温度，防止仪器过热或过冷。
工业炉设计：在钢铁或玻璃熔炼炉中，炉壁、炉顶和被加热物料之间存在着复杂的辐射换热。计算各表面间的视角系数，是设计炉膛形状、布置燃烧器以及评估热效率的基础。
建筑节能与日照分析：在绿色建筑设计中，需要计算窗户与室内各墙面、地板之间的辐射视角系数，以分析冬季的热量损失和夏季的太阳得热。这同样适用于计算城市峡谷中建筑物表面的日照率。
常见误解与注意事项
开始使用此工具时，有几个需要注意的要点。首先是“角系数与温度无关”。滑块上没有温度参数对吧？这是因为F₁₂纯粹是决定“几何可见性”的系数。实际的传热能力需要乘以发射率和温差后才能得到。例如，即使F₁₂=0.5，如果表面发射率较低，实际传递的热能也会较小。其次是“计算结果仅包含‘直接’辐射”这一限制。工具计算的是表面1直接可见于表面2的比例。在实际封闭空间（例如电子设备机箱内部）中，还存在经壁面多次反射后最终到达的“间接”部分。因此，最稳妥的做法是将本工具的结果视为“辐射传热量的下限值”。最后要注意参数的“无量纲化”。例如在“同轴圆盘”计算中，半径与“间距c”的比值（R₁/c, R₂/c）才是本质影响因素。直径50mm的圆盘与间距10mm（比值=5）的情况，和直径5m与间距1m（比值=5）的情况，其F₁₂值是完全相同的。请记住：关键不在于具体尺寸，而在于“形状比例”。
// Radiation view factor animation
(function() {
const el = document.getElementById('radViewAnimCanvas');
const ctx = el.getContext('2d');
let animT = 0;
const rays = [];
const NRAYS = 40;
function getParams() {
const sel = document.getElementById('geom');
const config = sel ? sel.value : 'parallel_rect';
const F12El = document.getElementById('statF12');
const F12 = F12El ? parseFloat(F12El.textContent) || 0.5 : 0.5;
return { config, F12 };
}
function resize() {
const dpr = window.devicePixelRatio || 1;
const w = el.parentElement.clientWidth - 36;
const H = Math.round(w * 0.38);
if (Math.abs(el.width - w * dpr) > 2 || el.height !== H * dpr) {
el.width = w * dpr;
el.height = H * dpr;
el.style.height = H + 'px';
ctx.setTransform(dpr, 0, 0, dpr, 0, 0);
rays.length = 0;
for (let i = 0; i < NRAYS; i++) {
rays.push({ life: Math.random(), speed: 0.008 + Math.random() * 0.01, angle: 0, ox: 0, oy: 0, hits: false });
}
}
return { W: w, H };
}
function frame() {
const { W, H } = resize();
const { config, F12 } = getParams();
animT += 1/60;
ctx.clearRect(0, 0, W, H);
ctx.fillStyle = '#001F3F';
ctx.fillRect(0, 0, W, H);
const padL = 20, padR = 20;
const surf1X = padL + 10;
const surf2X = W - padR - 10;
const surf1Y1 = H * 0.2, surf1Y2 = H * 0.8;
const surf2Y1 = H * 0.2, surf2Y2 = H * 0.8;
const cx = W / 2;
// Surface 1 (emitter) — hot
ctx.fillStyle = '#d63031';
ctx.fillRect(surf1X - 4, surf1Y1, 8, surf1Y2 - surf1Y1);
ctx.fillStyle = 'rgba(255,100,50,0.15)';
const grad = ctx.createLinearGradient(surf1X, 0, cx, 0);
grad.addColorStop(0, 'rgba(255,80,0,0.2)');
grad.addColorStop(1, 'rgba(0,0,0,0)');
ctx.fillStyle = grad;
ctx.fillRect(surf1X, surf1Y1, cx - surf1X, surf1Y2 - surf1Y1);
// Surface 2 (receiver) — blue
ctx.fillStyle = '#007BFF';
ctx.fillRect(surf2X - 4, surf2Y1, 8, surf2Y2 - surf2Y1);
// Labels
ctx.fillStyle = '#e17055';
ctx.font = 'bold 11px Roboto Mono';
ctx.textAlign = 'center';
ctx.fillText('面1 (放熱)', surf1X, surf1Y1 - 8);
ctx.fillStyle = '#007BFF';
ctx.fillText('面2 (受熱)', surf2X, surf2Y1 - 8);
ctx.fillStyle = 'rgba(255,255,255,0.6)';
ctx.font = '10px Roboto Mono';
ctx.fillText('F₁₂ = ' + F12.toFixed(3), cx, H - 8);
// Animate rays
for (let r of rays) {
r.life += r.speed;
if (r.life >= 1) {
r.life = 0;
// Originate randomly from surf1
r.ox = surf1X + 4;
r.oy = surf1Y1 + Math.random() * (surf1Y2 - surf1Y1);
// Random direction toward right hemisphere
const baseAngle = -Math.PI/2 + Math.random() * Math.PI;
r.angle = baseAngle;
// Hits surf2 if within F12 probability (more rays hit if F12 is higher)
r.hits = Math.random() < F12;
}
const t = r.life;
// Ray travels from ox,oy at angle
const len = W * 0.85;
const dx = Math.cos(r.angle) * len * t;
const dy = Math.sin(r.angle) * len * t;
const rx = r.ox + dx;
const ry = r.oy + dy;
// Check if ray has gone past surf2
if (rx >= surf2X - 4 && r.hits) {
// hit — stop at surf2
const hitAlpha = (1 - r.life) * 0.8;
ctx.strokeStyle = `rgba(255,200,0,${hitAlpha})`;
} else if (rx >= surf2X && !r.hits) {
// miss — continue past
ctx.strokeStyle = `rgba(255,80,0,${(1 - r.life) * 0.3})`;
} else {
ctx.strokeStyle = `rgba(255,${Math.round(150 - t*100)},0,${(1-t)*0.6})`;
}
ctx.lineWidth = 1;
ctx.beginPath();
ctx.moveTo(r.ox, r.oy);
ctx.lineTo(rx, ry);
ctx.stroke();
}
requestAnimationFrame(frame);
}
frame();
})();