CFD中文网

qjh888

你应该是安装 swak4foam
注意，swak4foam需要使用bison 2.7。
所以在安装以前，先确定一下bison的版本：

$ bison -V

如果bison的版本大于2.7，比如 3.0.4，运行以下：

$     cd $FOAM_RUN/swak4Foam
$     ./maintainanceScripts/compileRequirements.sh
$     ./Allwmake

安装完成之后运行一下：

$   funkySetFields

如果正确安装的话会出现类似于运行icoFoam的界面。

Janry

qjh888

这个是foam-extend 中的一个bug. 我已经repot了。
请点击这儿->Bug#308
有兴趣的朋友可以看看

qjh888

问题已经得到解决。
原因是Foam extended 3.0 使用的是mixingPlane RC1 的版本，如果使用Density based solver 的话，会出现一系列的错误。原因参考我在上一条回复中po出的Cfd-online的连接。
解决方法：
安装foam extend 3.1 和 foam extend 3.2，使用foam extend 3.2 中的mixingPlane文件（RC2版本）更新3.1中的C文件，然后：
$ cd $FOAM/foam-extend-3.1/src
$ wmake libso finiteVolume

Janry

qjh888

问题已经解决。
这个问题确实是由于网格过于稀疏导致的。
由于在OpenFOAM中网格是6面体结构网格，一个环形会离散成多边体，滑移网格每转一个角度，两个接触面之间的距离就会发生波动。当网格密度过低时，在特定的时间步下就会造成间距过大，GGI interpolation就会报错，找到没有被完全覆盖的面。

对于这个问题的解决方法是：细化网格。
Janry

qjh888

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

qjh888

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

qjh888

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

qjh888

此问题请管理员关闭。。已解决

qjh888

为了方便，我把我的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;
}

qjh888

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

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

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

qjh888

我想请教一下各位大神，在计算Laval喷管的时候，p,T,U的边界条件分别如何设置？
我使用rhoCentralFoam, sonicFoam, transonicMRFDyMFoam 均没有正确的得到激波。
p进口使用过： fixedValue, totalPressure;
p的出口使用： waveTansmissive;
T进口使用过：fixedValue, totalTemperature
U进口使用： zeroGradient,
对称面使用slip.
不过目前为止，均没有得到任何激波，而同样的网格，使用Ansys可以轻松得到激波。
请问，OF算laval喷管，是否有什么特定的要求？
谢谢

qjh888

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

qjh888

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

qjh888

各位朋友：

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

我的设置如下：

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

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

qjh888

@qjh888

Best posts made by qjh888

Latest posts made by qjh888