移植sprayFoam液相功能

[流体中二少年]

新手改求解器，由于老板要求引入液相功能，最近尝试把sprayFoam的液相功能移植到rhoReactingCentralFoam（rhoCentralFoam添加了reactingFoam的反应功能）上，编译求解器时，报的错误意思是方程数据类型不一致，我也大概知道了问题出在哪里，

solve
(
    fvm::ddt(rhoU)
  + fvc::div(phiUp)
  - rho()*g
  - parcels.SU(U)
);

这是动量方程的代码，我的理解是里面既有标量场，又有矢量场，所以数据类型不一致，希望老师们指点以下，给点思路，其实在能量方程那里遇到同样的问题

[Tens]

@流体中二少年 报错信息贴出来才好看

[流体中二少年]

@tens 不好意思，隔这么久回复你，求解器已经改好了，问题在于矩阵的形式不同，接下来准备验证求解器的可靠性