CFD中文网

Update: compressibleInterFoam是“使用基于VOF（流体体积）相分数的界面捕获方法求解两种可压缩的非等温不混溶流体，并进行可选的网格运动和网格拓扑更改（包括自适应重新网格划分）。” 并不适合上述问题。

@东岳 谢谢李老师。最近在研究一个求解器，compressibleInterFoam. 我把您推荐的这个也研究一下，看看是否可以求解我想要的问题。

RT
本人想计算两种气体的混合工质（如He + CH4）气动问题，要求可以可压、湍流、捕捉组份变化。请教有没有不错的求解器可以使用。
谢谢！

为了方便，我把我的fvScheme贴在这儿：

ddtSchemes
{
    default          none;//EulerLocal physDeltaT CoDeltaT;
    ddt(rho)       EulerLocal physDeltaT CoDeltaT;
    ddt(rhoU)      EulerLocal physDeltaT CoDeltaT;
    ddt(rhoE)      EulerLocal physDeltaT CoDeltaT;
    ddt(rho,omega) EulerLocal physDeltaT CoDeltaT;
    ddt(rho,k)     EulerLocal physDeltaT CoDeltaT;
    p              steadyState;   // Needed for waveTransmissive;
    phiMin           steadyState;   // Needed for waveTransmissive;
    phiMax          steadyState;    // Needed for waveTransmissive;
    
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
    grad(rhoSlope)  Gauss linear;
    grad(pSlope)    Gauss linear;
    grad(USlope)    Gauss linear;
    grad(TSlope)    Gauss linear;
    snGradCorr(U)   Gauss linear;
    snGradCorr(h)   Gauss linear;
    snGradCorr(omega) Gauss upwind phi;
    snGradCorr(k)   Gauss upwind phi;
    grad(omega)     Gauss upwind phi;
    grad(k)         Gauss upwind phi;
}

divSchemes
{
    default                         none;
    div(phi,U)                      Gauss linearUpwind grad(U);
    div(phi,k)                      Gauss QUICK;//linearUpwind;
    div(tau)                        Gauss linear;
    div((tau&U))                    Gauss linear;
    div(phi,omega)                  Gauss QUICK;//linearUpwind;

}

laplacianSchemes
{
    default         Gauss linear limited  0.33;
}

interpolationSchemes
{
    default               none;
    interpolate(rho) linear corrected;
}

snGradSchemes
{
    default         orthogonal;
}

fluxRequired
{
    default         no;
}

各位朋友，最近我在计算一个稳态的Laval喷管的稳态问题。可是最后的马赫盘的位置一直在非常微小幅度的振荡，导致最后的残差一直居高不下。我尝试过如下的几种方法：

1。 waveTransmissive改为稳态适用版本；
2。 减少co数和松弛因子。
3。 加密网格。

我用的是kOmegaSST, 求解器比较小众，是transonicMRFDyMFoam. 请问各位在稳定稳态求解方面有什么好的经验？
谢谢各位！

各位，我想问一下你们使用janafThermo 和surthlandTaransport的经验。我之前的一个用constant Thermo和ideal gas equation 下收敛的算例，现在换成janaf之后直接发散，完全没有办法去让它收敛。请问这种情况要怎么办呢？
谢谢

题主可以搜索一下 nasa gas polynomial.
不知道这个是不是你要的。

各位朋友：

我想问一下，为什么在出口处使用inletOutlet 边界条件并不能限制回流呢？
如图所示，有一个cell出现了巨大的回流，而且并没有被 inletOutlet 设置成为 （0 0 0）。

我的设置如下：

patch_outlet
{
    type     inletOutlet;
    inletValue  uniform(0 0 0); 
}

各位是否有人遇到过这种情况？
谢谢！

@李东岳

对的，您比较一下就知道了。Extend什么也没有。。。

@李东岳
谢谢东岳大神。

我看了一下OpenFOAM的代码，它的wallShearStress 可以求解可压流动。
只是再确认一下您说的最后一句话：“壁面剪切力，既然是壁面，只求解壁面，因此你说的是正确的，类似的处理在边界条件和壁面函数上很常见”
也就是说虽然求解所有的patch,但是只有wall的边界条件是有意义的。是否是这样呢？

目前知道第一句，就是
RASModel->devReff() 是//- Return the effective stress tensor including the laminar stress
那么问题来了，这个wallShearStress 后处理器是否只能适用于incompressible flow?

有人知道怎么使用不？

此问题可以关闭了。
答案在这儿：
https://www.cfd-online.com/Forums/openfoam-programming-development/121114-foamtocgns-cgnstofoam-openfoam-2-1-1-a.html

各位好！
最近正在使用wallShearStress这个后处理器。
其中有一些代码并不知道是什么作用，还望各位大神不吝赐教。

RASModel->devReff()//这一句确实不知道是用来干什么的


//下面这几句分别是什么作用呢？
singlePhaseTransportModel laminarTransport(U, phi);

autoPtr<incompressible::RASModel> RASModel
(
    incompressible::RASModel::New(U, phi, laminarTransport)
);


//下面的是对所有的patch, 也就是说不仅仅是wall，进行求解？
forAll(wallShearStress.boundaryField(), patchi)
    {
        wallShearStress.boundaryField()[patchi] =
        (
            -mesh.Sf().boundaryField()[patchi]
            /mesh.magSf().boundaryField()[patchi]
        ) & Reff.boundaryField()[patchi];
    }

非常感谢！

@random_ran

谢谢！
我研究一下，然后反馈

@random_ran
谢谢，不过我需要大量的切面，所以只能用编程来实现…

@赵一铭

:crying:
看来这个sampleDict是不能完成这种工作了…唉～

@赵一铭 在 sampleDict 是否可以切取圆柱面采集数据？ 中说：

定义你的圆环面在这里

嗯，谢谢，但是有一个问题，这个integrationPlan.stl如何去定义呢？是要自己选取一个圆环面么？

谢谢！

各位朋友，

请问一下谁用过这个转换器？或者，有什么方法可以将OpenFOAM格式的网格转换成Cgns标准格式的网格？
谢谢！

各位好！

我现在在尝试着使用OpenFOAM切取一个圆环面上的数据。我看sampleDict中有一项功能是：

    distance
    {
        // Isosurface from signed/unsigned distance to surface
        type            distanceSurface;
        signed          true;

        // Definition of surface
        surfaceType     triSurfaceMesh;
        surfaceName     integrationPlane.stl;
        // Distance to surface
        distance        0.0;

        interpolate     false;
        regularise      false;       // Optional: do not simplify
        // mergeTol        1e-10;    // Optional: fraction of mesh bounding box
                                     // to merge points (default=1e-6)
}

不知道如何使用，是否有朋友使用过这个功能？

谢谢！