关于rhoCentralFoam计算高超声速问题时温度出现负值的原因

[Foamer24]

@zxyf1 rhoCentralFoam改了好多边界条件和参数，最后一直报相同的温度错误，不知道出错的地方在哪里，不知道这种热力学边界要怎么调比较好，是不是这个求解器本身鲁棒性不好。

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib64/libc.so.6"
#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5  ? at ??:?
#6  __libc_start_main in "/lib64/libc.so.6"
#7  ? at ??:?
Floating point exception (core dumped)

[zxyf1]

@Foamer24 就是求解器本身鲁棒性不好，你可以尝试调换flux limiter，应该在fvScheme里面。

[Foamer24]

@zxyf1 rhoSonicFoam and rhopSonicFoam have been superceded by rhoCentralFoam. 换成rhoPimpleFoam、pisoFoam或者sonicFoam可能会好一些

请问flux limiter 具体指的是%(#ff0000)[fluxScheme Kurganov;]？

[Foamer24]

@zxyf1 ε=(´ο｀*)))唉，将fluxScheme Kurganov; 改为fluxScheme Tamdor;算的时间会变长一点，但是也会出现上面的问题。

[adcpk]

请问插值格式用的什么

[Foamer24]

@adcpk 您好！插值格式用的是linear

interpolationSchemes
{
    default         linear;
}

[adcpk]

@Foamer24 应该不是这个default吧，后面主要物理量都覆盖了

[adcpk]

@Foamer24 抱歉，是divSchemes

[Foamer24]

@adcpk 在 关于rhoCentralFoam计算高超声速问题时温度出现负值的原因 中说：

divSchemes

您好！divSchemes用的是Euler

ddtSchemes
{
    default         Euler;
}

[youmengtian]

@Foamer24 你可以贴一个fvSchemes完整的设置看一下，很有可能是离散格式的问题。

[Foamer24]

@youmengtian 您好！我完整的格式是这样的

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluxScheme          Tadmor;// Kurganov;

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
         default cellMDLimited Gauss linear 1.0;
}

divSchemes
{
    default         none;
    div(tauMC)      Gauss linear;
    div(phi,nuTilda) Gauss linearUpwind limited;
}



laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}


snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}