在OpenFOAM的sprayFoam案例求解中，怎么设置入射液滴的温度？

mingyang

是在OpenFOAM中的aachenBomb案文件夹下constant/cloudProperties文件里设置入射液滴的温度吗？下面截图的constantProperties内T0是什么意思？路过的大佬指点一二。
QQ截图20231117215707.png

wangfei9088

@mingyang 这个T0是入射燃料的温度，也用来初始化parcel（液滴）的温度，具体见ThermoParcel.H和ThermoParcelI.H。

 template<class ParcelType>
 inline Foam::ThermoParcel<ParcelType>::constantProperties::constantProperties
 (
     const dictionary& parentDict
 )
 :
     ParcelType::constantProperties(parentDict),
     T0_(this->dict_, "T0"),           //T0在这里
     TMin_(this->dict_, "TMin", 200.0),
     TMax_(this->dict_, "TMax", 5000.0),
     Cp0_(this->dict_, "Cp0"),         //Cp0在这里
     epsilon0_(this->dict_, "epsilon0"),
     f0_(this->dict_, "f0")
 {}

顺便说一句，每个时间步的parcel温度是通过calcHeatTransfer函数更新的，但还是赋值给这个T0，也就是说它也是更新的。虽然喷雾燃烧涉及传热相变，但是这个过程是与液滴表面温度相关的，液滴的温度和环境温度相对来说影响小一些。所以，个人觉得第一个时间步的T0（初始）值的影响应该不大。

另外，给定的密度rho0在KinematicCloud.C和KinematicParcel.H里有用到。比热容Cp0也要给定，在ThermoParcel.H和ThermoCloud.C里。constantVolume体现在ReactingParcel.H里。

mingyang

谢谢大佬的讲解，我大概明白了，如果在设置的时候T0/Cp0/rho0三者的时候，我如果没有正确给定在对应T0温度下的Cp0或者rho0，计算过程是不是会有发散的情况，因为我计算时出现过下面截图的FOAM Warning，在给定物性的therm.dat文件中温度范围是200-6000。但是计算过程好像出现低于200k的情况。最后，再次感谢大佬的讲解。
QQ截图20231117140409.png

wangfei9088

@mingyang 初始T0的值可能会有warning，但应该不会导致发散，因为时间步推进一次就正常了。图里有很多网格里的T都低于200K，还特别低。没有温升，没有燃烧，应该有能量方程方面相关的问题。

mingyang

好的，谢谢大佬的讲解，我再看看相应求解器下的能量方程。感谢

CFD中文网

在OpenFOAM的sprayFoam案例求解中，怎么设置入射液滴的温度？