光线折射公式推导：Snell‘s Law

光线折射公式推导：Snell’s Law

设 入射光 $\mathbf{i}$ 和 折射光 $\mathbf{t}$ 均为单位向量，由Snell’s law可知 ${\eta }_i\sin {\theta }_i={\eta }_t\sin {\theta }_t$ 。

其中 ${\eta }_i,{\eta }_t$ 为反射系数（typically air = 1.0, glass = 1.3–1.7, diamond = 2.4）。

因为 $\sin {\theta }_i=\frac{|{\mathbf{i}}_{\parallel }|}{|\mathbf{i}|},\sin {\theta }_t=\frac{|{\mathbf{t}}_{\parallel }|}{|\mathbf{t}|}$ ，故：
$$\begin{array}{c}{\eta }_i\frac{|{\mathbf{i}}_{\parallel }|}{|\mathbf{i}|}={\eta }_t\frac{|{\mathbf{t}}_{\parallel }|}{|\mathbf{t}|}\\ {\eta }_i|{\mathbf{i}}_{\parallel }|={\eta }_t|{\mathbf{t}}_{\parallel }|\end{array}$$
又因 $\mathbf{i}$ 和 $\mathbf{t}$ 平行且方向相同：
$$\begin{array}{cccc}{\eta }_i\cdot {\mathbf{i}}_{\parallel }& =& {\eta }_t\cdot {\mathbf{t}}_{\parallel }& \\ {\mathbf{t}}_{\parallel }& =& \frac{{\eta }_i}{{\eta }_t}\cdot {\mathbf{i}}_{\parallel }& \\ & =& \frac{{\eta }_i}{{\eta }_t}(\mathbf{i}+\cos {\theta }_i\cdot \mathbf{n})& \end{array}$$
其中 ${\mathbf{i}}_{\parallel }=\frac{\mathbf{i}\cdot \mathbf{n}}{|\mathbf{i}||\mathbf{n}|}\mathbf{n}=\cos {\theta }_i\cdot \mathbf{n}$.

又因 $\mathbf{t}={\mathbf{t}}_{\perp }+{\mathbf{t}}_{\parallel }$ 且 $|\mathbf{t}{|}^2=|{\mathbf{t}}_{\perp }{|}^2+|{\mathbf{t}}_{\parallel }{|}^2$ ：
$$\begin{array}{ccc}|{\mathbf{t}}_{\perp }{|}^2& =& |\mathbf{t}{|}^2-|{\mathbf{t}}_{\parallel }{|}^2\\ |{\mathbf{t}}_{\perp }|& =& \sqrt{1-|{\mathbf{t}}_{\parallel }{|}^2}\\ {\mathbf{t}}_{\perp }& =& -\sqrt{1-|{\mathbf{t}}_{\parallel }{|}^2}\cdot \mathbf{n}\end{array}$$
从而：
$$\mathbf{t}={\mathbf{t}}_{\parallel }+{\mathbf{t}}_{\perp }$$

// r_in 为射入物体的光线
vec3 refract(const vec3& r_in, const vec3& outward_normal, double etai_over_etat) {
    double cos_thetai = dot(-r_in, outward_normal);
    vec3 r_out_parallel =  etai_over_etat * (r_in + cos_thetai * outward_normal);
    vec3 r_out_perp = -sqrt(1.0 - r_out_parallel.length_squared()) * outward_normal;
    return r_out_perp + r_out_parallel;
}

也可进一步化简：
$$\begin{array}{ccc}\mathbf{t}& =& {\mathbf{t}}_{\parallel }+{\mathbf{t}}_{\perp }\\ & =& \frac{{\eta }_i}{{\eta }_t}\cdot (\mathbf{i}+\cos {\theta }_i\cdot \mathbf{n})-\sqrt{1-|{\mathbf{i}}_{\parallel }{|}^2}\cdot \mathbf{n}\\ & =& \frac{{\eta }_i}{{\eta }_t}\cdot (\mathbf{i}+\cos {\theta }_i\cdot \mathbf{n})-\sqrt{1-{\sin }^2{\theta }_t}\cdot \mathbf{n}\\ & =& \frac{{\eta }_i}{{\eta }_t}\cdot \mathbf{i}+\left(\frac{{\eta }_i}{{\eta }_t}\cos {\theta }_i-\sqrt{1-{\sin }^2{\theta }_t}\right)\cdot \mathbf{n}\end{array}$$
由Snell’s law公式可得 $\sin {\theta }_t=\frac{{\eta }_i}{{\eta }_t}\cdot \sin {\theta }_i$，
$$\begin{array}{ccc}\mathbf{t}& =& \frac{{\eta }_i}{{\eta }_t}\cdot \mathbf{i}+\left(\frac{{\eta }_i}{{\eta }_t}\cos {\theta }_i-\sqrt{1-(\frac{{\eta }_i}{{\eta }_t}{)}^2{\sin }^2{\theta }_i}\right)\cdot \mathbf{n}\\ & =& \frac{{\eta }_i}{{\eta }_t}\cdot \mathbf{i}+\left(\frac{{\eta }_i}{{\eta }_t}\cos {\theta }_i-\sqrt{1-(\frac{{\eta }_i}{{\eta }_t}{)}^2(1-{\cos }^2{\theta }_i)}\right)\cdot \mathbf{n}\end{array}$$
其中 $\cos {\theta }_i=\frac{\mathbf{i}\cdot \mathbf{n}}{|\mathbf{i}||\mathbf{n}|}=\mathbf{i}\cdot \mathbf{n}$ 。

全反射现象 Total internal reflection

对于Snell’s law 的公式 $\sin {\theta }_t=\frac{{\eta }_i}{{\eta }_t}\cdot \sin {\theta }_i$ ，当空气中的反射率大于物体内的反射率的时候，有可能会出现 $\frac{{\eta }_i}{{\eta }_t}\cdot \sin {\theta }_i>1$ 的情况。

但是 $\sin {\theta }_t$ 的值一定是小于1的，故等式不可能成立。这个时候光线应该被反射，而不是折射。