麻烦哪位大佬遇到过类似的问题,帮忙解决一下,谢谢!
问题描述:
利用inter(DyM)Foam一个圆柱射流进入大空间:
通过blockMesh进行网格划分,利用codeFixedValue设置圆柱形液柱的入口面。
U边界:
type codedFixedValue;
value uniform (0 0 0);
name liquidColumnVelocity;
code #{
const scalar columnVelocity = 10;
const scalar radiusLiquidColumn = 4e-4;
const vectorField& Cf = patch().Cf();
const scalar yc = 0.01;
const scalar zc = 0.0;
forAll(Cf,faceI)
{
const scalar deltay = Cf[faceI].y()-yc;
const scalar deltaz = Cf[faceI].z()-zc;
const scalar radius = sqrt(deltay*deltay+deltaz*deltaz);
if(radius<=radiusLiquidColumn)
{
(*this)[faceI] = vector(columnVelocity,0,0);
}
}
#};
alpha.water边界:
type codedFixedValue;
value uniform 0;
name liquidColumnAlpha1;
code #{
const scalar radiusLiquidColumn= 4e-4;
const vectorField& Cf = patch().Cf();
const scalar yc = 0.01;
const scalar zc = 0.0;
forAll(Cf,faceI)
{
const scalar deltay = Cf[faceI].y()-yc;
const scalar deltaz = Cf[faceI].z()-zc;
const scalar radius = sqrt(deltay*deltay+deltaz*deltaz);
if(radius<=radiusLiquidColumn)
{
(*this)[faceI] = 1.0;
}
}
#};
报错描述
当不使用自适应网格时,算例正常计算;当采用自适应网格功能interDyMFoam,网格未加密前(未达到设置加密条件),正常计算;网格开始加密计算中断,报错如下:
Courant Number mean: 4.98466e-06 max: 0.00220318
Interface Courant Number mean: 0 max: 0
deltaT = 6.18034e-08
Time = 3.11186e-07
Selected 0 cells for refinement out of 384000.
Selected 0 split points out of a possible 0.
smoothSolver: Solving for alpha.water, Initial residual = 0.061072, Final residual = 2.16265e-15, No Iterations 2
Phase-1 volume fraction = 3.73802e-07 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 3.73802e-07 Min(alpha.water) = -8.38315e-28 Max(alpha.water) = 1
DICPCG: Solving for p_rgh, Initial residual = 0.109645, Final residual = 0.00458278, No Iterations 5
time step continuity errors : sum local = 8.25061e-11, global = 4.48406e-13, cumulative = 6.73588e-12
DICPCG: Solving for p_rgh, Initial residual = 0.0059725, Final residual = 0.000280082, No Iterations 19
time step continuity errors : sum local = 4.64407e-12, global = 1.03878e-12, cumulative = 7.77466e-12
DICPCG: Solving for p_rgh, Initial residual = 0.000280406, Final residual = 9.32371e-08, No Iterations 86
time step continuity errors : sum local = 1.54319e-15, global = 1.34038e-17, cumulative = 7.77467e-12
ExecutionTime = 38.03 s ClockTime = 46 s
Courant Number mean: 5.97925e-06 max: 0.0026174
Interface Courant Number mean: 8.02695e-08 max: 0.0026174
deltaT = 7.41181e-08
Time = 3.85304e-07
Selected 12 cells for refinement out of 384000.
Refined from 384000 to 384084 cells.
Selected 0 split points out of a possible 12.
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 8.07683e-06, No Iterations 124
time step continuity errors : sum local = 7.03583e-13, global = -1.22968e-15, cumulative = 7.77345e-12
smoothSolver: Solving for alpha.water, Initial residual = 0.00350774, Final residual = 8.77793e-10, No Iterations 1
Phase-1 volume fraction = 3.96963e-07 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 3.962e-07 Min(alpha.water) = -2.57151e-25 Max(alpha.water) = 1
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in “/lib64/libc.so.6”
#3 double Foam::sumProd(Foam::UList const&, Foam::UList const&) at ??:?
#4 Foam::PCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const at ??:?
#5 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) at ??:?
#6 Foam::fvMatrix::solve(Foam::dictionary const&) at ??:?
#7 Foam::fvMatrix::solve() at ??:?
#8 ? at ??:?
#9 __libc_start_main in “/lib64/libc.so.6”
#10 ? at ??:?
Floating point exception (core dumped)
请问论坛上的各位大佬怎么解决?万分感谢!