我目前正在开展一个涉及管道两相流的项目，我需要在MatLab中对其进行建模。即使我已经开始并得到一些相关的结果，我感觉有些事情是不对的，我认为这可能是我解决平流方程的方式。
我想解决的等式是
d(Hl x rhoL)/dt + d(Hl x rhoL x Vl)/dx = 0 ......其中：Hl =液体分数;
rhoL =液体密度;
V1 =液体速度;
液体不可压缩我没有考虑密度，因为它是恒定的。
我试图在时间t + dt（速度不恒定）的液体速度求解这个等式。
到目前为止我所做的是考虑Hl = U和Hl x Vl = V并且使用时间上的向前差异和空间的向后差异我最终得到以下公式：
Vl(i,1)=(((Hlold(i,1)-Hl(i,1))*(deltaX/deltaT))+(Vl(i-1,1)*Hl(i-1,1)))/Hl(i,1));
你们认为这是对的吗？有人可以帮助我理解如何通过节点将具有变化速度的PDE分散吗？
非常感谢！