CFD中文网

本人用sprayFaom进行一些模拟，请问如何获得液滴的蒸发速率？
质量方程：
代码：

 solve
(
    fvm::ddt(rho)
  + fvc::div(phi)
  - parcels.Srho(rho)
);

论文中说Smass是蒸发速率，那我需要把代码中的parcels.Srho(rho)给输出出来，如何实现呢？

前辈，你好！请问sprayFoam中的液滴蒸发率怎么获得？

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

刚接触sprayfoam，请教三个个问题：

1. 初始场指定区域怎么布置液滴？setFields？
2. 液滴的喷射位置和喷射速度怎么设置为跟随inlet边界？
3. 我模拟的是预气化的液滴，比如C7H16既有液滴也有气体，我在算例里面该怎么区分他们？

谢谢前辈，希望解答或者提供一些解决问题的方向

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

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

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

@strive_dut 非常感谢，我试试

@strive_dut 大哥，你最后是怎么解决的，我也遇到和你相同的问题

if (p<5 )               
			{
                    	field[faceI] = vector(0,0,0);
                	}
		else
                	{
                    	field[faceI] = vector(U_0*t,0,0);
                	}
		}

@chien 请问用forall如何比较，能否提供下code，万分感谢

小白一枚，在下使用codedFixedValue自定义边界条件，根据边界处不同的压力，给定入口对应的速度
网上查到提取时间的代码是scalar t = this->db().time().value();
那压力要如何提取，请各位前辈不吝赐教。

时隔多年，没想到在这遇到孟前辈，老师让我复现你的RDE数值模拟的论文，入口边界用swak4Foam完美解决了，出口边界你论文中提到的无反射边界条件该怎么设置，waveTransmissive这个边界条件可以自动实现超音速和亚声速转换吗