CFD中文网

请问有人在of中成功使用动态增厚火焰模型吗

请问您成功了吗，我也有这方面的需要

以及对于入口气态氧化剂超声速射入燃烧室，入口液态燃料低速进入燃烧室进行非预混燃烧，最后速度为低马赫数的用sprayFoam进行燃烧仿真合适吗
请各位给小白点建议

@cccrrryyy 我想请问一下您对这个“如果采用低马赫数近似，密度随压力不变，那是不是就相当于认为声速无穷大？这样假设会导致什么呢？”有个结论了吗，对于发动机来说，燃料进入燃烧室通常是速度较快的，但燃料开始燃烧之后速度又是正常低马赫数，请问这种用reactingFoam合适吗

同问，我也出现了这个问题，不知道该不该去调整

@King_RVM 方便问您一些具体的网格导入和修改边界条件的细节吗，因为我现在导入网格改边界条件总是失败

求问大佬如何绘制楔形模型的，是在icem里画完了导入openfoam的嘛，再改的wedge的边界条件嘛，我是要算一个轴对称圆柱燃烧室

@李东岳 好的，非常感谢东岳老师

原模型是圆柱燃烧室，为了简化计算只画了1/12的部分
我个人想法是将两个面改为wedge，但是checkMesh的时候发现识别不了，会报错，求大神指教

报错如下：
Checking geometry...
Overall domain bounding box (0 -6.363451432e-19 -1.169089492e-12) (0.2372166655 0.01928534508 0.009642672539)
Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
Mesh has 3 solution (non-empty) directions (1 1 1)
Wedge WALL_RANSHAOSHI2 with angle 0 degrees
***Cannot find opposite wedge for wedge WALL_RANSHAOSHI2
Boundary openness (-1.075387266e-17 -1.69211467e-15 -1.2396473e-14) OK.
Max cell openness = 4.03570353e-16 OK.
Max aspect ratio = 7.562500293 OK.
Minimum face area = 4.890642405e-09. Maximum face area = 2.500000033e-07. Face area magnitudes OK.
Min volume = 6.010518271e-13. Max volume = 9.284234891e-11. Total volume = 1.511515734e-05. Cell volumes OK.
Mesh non-orthogonality Max: 63.85295712 average: 15.3592327
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.8283279386 OK.
Coupled point location match (average 0) OK.

Failed 1 mesh checks.

又又又遇到了新问题，一步都没算下去，请问这种问题一般出在哪里
Create time

Create mesh for time = 0

PIMPLE: No convergence criteria found

PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode

Using LTS
Reading thermophysical properties

Selecting thermodynamics package
{
type heRhoThermo;
mixture multiComponentMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model.

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Prt 0.85;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
Creating reaction model

Selecting combustion model EDC
Selecting chemistry solver
{
solver ode;
method standard;
}

standardChemistryModel: Number of species = 544 and reactions = 2446
Selecting ODE solver seulex
Creating field dpdt

Creating field kinetic energy K

No MRF models present

No fvModels present
No fvConstraints present

Starting time loop

Time scales min/max:
Flow = 9.728569846e-08, 0.0001
Temperature = 4.494232837e+307, 4.494232837e+307
Composition = 4.494232837e+307, 4.494232837e+307
Overall = 9.728569846e-08, 0.0001
Time = 1e-07

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0001092756986, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.006862564884, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.0001167898495, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::standardChemistryModel<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::perfectGasFoam::specie >, Foam::sensibleEnthalpy> > >::solve<Foam::Field<double> >(Foam::Field<double> const&) at ??:?
#4 Foam::combustionModels::EDC::correct() at ??:?
#5 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#6 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception (core dumped)

@李东岳 是的东岳老师，发散的原因还在寻找，目前找到一个原因是有可能是我把出口当成壁面了，给的出口压力边界条件zeroGradient 所以压力出不去

@vbcwl 请问你调用成功了嘛，我也在尝试firefoam的这里，我想改成正庚烷的

好巧我也是！！跑了好几天然后出来浮点数溢出的报错，我的问题好像是压力一直无法收敛，一直再增长，这是我的部分报错，可以一起交流一下，我小白一只也在搞reactingfoam

@bestucan 好的好的，十分感谢！

@李东岳 好的，谢谢东岳老师，我去看看

这个问题解决了，但是出现了新的问题

No MRF models present

No fvModels present
No fvConstraints present
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::nutkWallFunctionFvPatchScalarField::nut() const at ??:?
#4  Foam::nutWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#5  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#6  Foam::RASModels::kEpsilon<Foam::CompressibleMomentumTransportModel<Foam::dynamicTransportModel> >::correctNut() at ??:?
#7  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
#8  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  ? in "/opt/openfoam9/platforms/linux64GccDPInt32Opt/bin/reactingFoam"
Floating point exception (core dumped)

目前研一刚入门阶段，针对reactingfoam进行一些改动，用的是自己的构型，边界条件也相应做了一些调整，原本的参考案例是用CH4，我自己改成了C7H16，也把化学机理加了进来，现在遇到了这个问题不知道怎么办（报错如下所附），请求各位大神指点一下

--> FOAM FATAL ERROR: 
Sum of mass fractions is zero for species 
544
(
各种组分
）
    From function void Foam::basicSpecieMixture::correctMassFractions()
    in file mixtures/basicSpecieMixture/basicSpecieMixture.C at line 205.

FOAM exiting

@香柏树 请问你用ReactingFoam ras改les成功了吗，想和您请教一下

@myheart 好的好的，十分感谢，我再去研究一下

一枚小白，请问有用过reactingfoam的大佬吗，有很多问题都想请教一下