CFD中文网

@xpqiu 老是您好，同求一份，感谢。834749050@qq.com

@尚善若水 对，用的PeleC。很抱歉，这个机理也是从别人那里来的，我没办法直接分享给你，不过我的经验是无论Li的机理还是Burke的机理，都能算出这个效果。

@钟余 不好意思，最近比较忙，没想到你已经做了这么多的工作。是这样的，计算域长点，到后面会稳定两个胞格，但是我用了取巧的方法，你稍微给大点宽度，现在是60mm，你给个61.8mm试试，然后给4个点火源。这是我最近用5阶weno算的结果： ，和openfoam算的差不多。

@钟余 你是想映射一段ZND爆轰作为初始流场？这个可以用sdtoolbox算一段一维爆轰数据，然后set到流场里面，这个不难。其次，这个两步反应模型，我没具体了解过，不过你可以参考下这个https://github.com/risinyoung/rhoCentral1StepFoam。

@钟余 原因比较复杂，大概有：网格尺寸(一般诱导区长度或者半反应区长度需要有10个以上的网格来解析)，化学机理(适合高温高压)，还有就是detonationFoam的rhoY并不没有进行通量计算，当然我不知这个影响大不大(这个你可以参考deepflame的通量格式)，还有就是和点火个数有关(这个你得自己试试换下点火个数)。我没跑过20kpa的算例，你可以试试经典Oran得6670pa算例。

@李东岳 在 移植sprayFoam液相功能 中说：

rhoReactingCentralFoam

感谢李老师的回复，之前用过rhoReactingCentralFoam算爆轰胞格结构，效果还可以，但是这个求解器不支持二维网格的AMR并且KNP格式耗散大了点，后面我用一个HLLC-type的黎曼求解器和用上AMR，对惰性/反应性多组分激波管，旋转爆轰，斜爆轰进行验证，效果也非常好，链接文本

@xuluxi18 我现在不搞两相爆轰了，但是据我所知现在用of算爆轰燃烧已经非常成熟了，github上面有很多开源的求解器，例如Deepflame，blastFoam，detonationFoam，这些求解器有显式RK时间推进，耗散更小的黎曼求解器以及支持二维网格的AMR等等，比这个rhoReactingCentralFoam功能强大很多，你可以试试，祝好！

本人用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这个边界条件可以自动实现超音速和亚声速转换吗