关于MPM以及APIC中矩阵C_p如何推出4/dx^2的问题,以及为什么使用张量积的问题

对于公式:

C_p^{n+1} = \frac 4 {{\Delta x}^2} \sum_i w_{ip} \mathbf{v}_i^{n+1}{\left(x_i - x_p^n\right)}^T

第一个问题:4/dx^2

GAMES201似乎并没有讲这个 \frac 4 {{\Delta x}^2} 如何推出来的,只说了MPM中这一项是从APIC里用Quadratic得来的,但是在上一讲APIC时也没有推导这一项

在APIC论文原文中,C矩阵是两个矩阵计算得出的:

C_p^n=B_p^n(D_p^n)^{-1} \\ D_p^n=\sum_i w_{ip}^n{\left(x_i - x_p^n\right)} {\left(x_i - x_p^n\right)}^T \\ B_p^{n+1}=\sum_i w_{ip}^n \tilde{\mathbf{v}}_i^{n+1} {\left(x_i - x_p^n\right)}^T

论文在后面又说在Quadratic差值形式下很简单可以化为 D_p^n = \frac 1 4 \Delta x^2 I , 从而得出GAMES201中计算 C_p^{n+1} 的形式。但是论文原文也没有说如何推出这个 \frac 1 4 \Delta x^2 I

我试了两整天也没办法把这个 \frac 1 4 \Delta x^2 I 推出来

问题2:张量积

APIC论文原文中提到使用 v(x) = v_p + C_p(x - x_p) 来近似估计粒子p邻域的速度场,并且GAMES201中也使用了一张直观的图片来解释 C_p 矩阵。
照这样看 C_p 矩阵应该类似于速度关于位置的导数关系, 在应用中应该使用解线性方程 v(x) = v_p + C_p(x - x_p) 的方式来求解 C_p

但是实际上,应用中使用了张量积来近似这个导数关系,并且通过除以 \Delta x^2 来解决了量纲对不上的问题。这非常的奇怪

以上两个问题,希望有高人解答