CFD中文网

@东岳
老师这是我目前算出来的一些结果：
链接: https://pan.baidu.com/s/163fY1uXF_sCCcN9nt1MiCQ
提取码: w3gp
网格量很少的，我现在只是想通过调试把路给走通，方便以后的使用

@东岳 不算是风场吧，是超声速管流，因为我模型前面还有很长一段模型，但前面的部分我不关注，并且为了节省计算资源，我就先算好前面这一步的出口，然后赋值给我这部分的出口，以后不同的工况都直接用这个边界层进口就可以了，也不需要每次都算前面一部分，这是我的本科毕设，先算直管道，再算一定曲率的和变截面的管道，后面再算一些比较成熟的工程模型，从二维到三维。我现在用的用的求解器是rhoCentralFoam,湍流模型是komegaSST,但我发现算出来的的结果和Fluent的差距有点大，一直在调试，现在能算出边界层，但感觉不是很理想。昨天看了老师的rhoCentralFoam解析，发现rhoCentralFoam是很适合于求解超声速可压缩流体的，但还是没有调试好，我这边有一些算好的结果，要是要是有空的话，我可以发给老师老师帮我参考一下~

@疾风GAVIN
我的问题解决了，问题出在文件类型上，文件类型应该是C源代码，而不是纯文本文档，现在我已经实现给入口添加边界层的功能了，谢谢各位老师的解答

谢谢老师前面的讲解，但我按照顺序给0文件夹中的各个参数的进口赋值了nonuniform List<scalar>和nonuniform List<vector>的边界条件,但我在运行和时候出现了以下错误：

--> FOAM FATAL IO ERROR: 
Expected a ')' while reading binaryBlock, found on line 21 the doubleScalar 4.7

file: /home/jifeng/OpenFOAM/jifeng-6/GraduateDesign/Straight_20_0_200_0.98U1BL/0/p at line 21.

    From function Foam::Istream& Foam::Istream::readEnd(const char*)
    in file db/IOstreams/IOstreams/Istream.C at line 109.

FOAM exiting

我参考了tutorials中用到同样设置的0文件夹，我没有找到设置的不正确的地方，下面是我的p文件：

FoamFile
{
	version             2.0;
	format              binary;
	arch                "LSB;label=32;scalar=64";
	class               volScalarField;
	location            "0";
	object              p;
}


dimensions	[1 -1 -2 0 0 0 0];
internalField				uniform 100000.0;
boundaryField
{
	inlet
	{
		type			fixedValue;
		value			nonuniform List<scalar> 
59
(
105559.3
106142.5
102728.5
104058.5
105988.4
104673.2
104496.1
104798.2
104722
104710.9
104719.2
104720.1
104730.8
104733.2
104734
104727.2
104725.6
104750.9
104773.8
104843.1
104878.7
104986.6
105056.3
105129.8
105207.8
105255.9
105315.9
105348.5
105349.4
105378.3
105346.6
105345.5
105313.3
105254.7
105207.5
105124.8
105053.1
105003.3
104907.8
104863.5
104812.3
104775
104754.5
104753.6
104754.6
104755.9
104752.2
104744.7
104736.2
104719.5
104742.3
104709.9
104497.9
105188.1
105279
103137.2
104778
106991.6
103429.3
)
;
	}
	outlet
	{
        	type            	waveTransmissive;
        	field           	p;
        	psi             	thermo:psi;
        	gamma           	1.4;
        	fieldInf        	100000.0;
        	lInf            	1;
        	value           	uniform 100000;
	}
	wall
	{
		type			zeroGradient;
	}
	frontAndBackPlanes
	{
		type			empty;
	}
}

我并没有发现和tutorials里面不一样的地方，希望大家能帮我看一下原因，谢谢

@东岳 谢谢老师！！！

@cccrrryyy 老师，我这里为什么 会显示 symbol lookup error呀

jifeng@jifeng-GL552VW:~/OpenFOAM/jifeng-6/yinJun/yinJun_2$ writeCellCentres 
writeCellCentres: symbol lookup error: writeCellCentres: undefined symbol: _ZNK4Foam11regIOobject11writeObjectENS_8IOstream12streamFormatENS1_13versionNumberENS1_15compressionTypeE

@cccrrryyy 老师您好，我觉得您的第二个方法能有帮助，而且从思路上我觉得是能实现的
但是老师我对于具体的操作方法还不太清楚，我想问一下老师说的Paraview里面的Cellcenters我会，但具体的writeCellCentres怎么实现，怎么输出我还是不太清楚，我想输出一个界面有U、T、k、omega这些参数然后再将它插值到进口面，这个过程我还打算修改其X坐标为0（为了与进口匹配），同时关于0文件的里面inlet的设置该用什么边界条件具体什么格式我还不知道，有点手足无措的感觉，UserGuide上面讲解的也很粗造，不知道具体该怎么设置，想向老师请教一下老师具体是怎么实现的？谢谢老师~

嗯嗯我觉得老师说的这个方法很有启发，我准备先尝试一下第一种方法，谢谢老师

我计算的对象是超声速可压缩管流，因为实际情况中入口是有已经发展好的边界层的，所以如果入口条件直接给定速度或者压力，算出来的结果和真实的有边界层的入口是有很大差别的，所以希望能在计算开始的时候就在入口添加一定管道百分比厚度的边界层，我有两个思路，但不知道怎么在OpenFOAM里面实现：
1.是直接给一定厚度白分比厚度的边界层，我在一片韩国人的论文里看到过这样的设置

摘自Numerical study on the start and unstart phenomena in a scramjet inlet-isolator model
但是我不知道他是怎么设置的，向各位老师和同学求助
2.截取已经计算好的一个有边界层的Profile直接赋给进口，这个是我使用Fluent里面的用到的操作，相比自己设置入口边界层而言很方便，不需要自己给定壁面温度和压力，直接根据计算的Profile里面的 结果来就可以了，但是这个操作我也不知道在OpenFOAM里面应该怎么实现，我知道好像有一个setFieldsDict什么的文件，但不知道具体怎么设置，要是各位老师和同学知道怎么用的，希望教我一下，非常感谢！

@东岳 谢谢东岳老师

使用OpenFOAM时间不长，小白一枚，请问大家OpenFOAM使用Ctrl+C暂停后，该怎么继续计算？暂停对计算后再计算对计算结果是否会有影响？
非常感谢

哦哦，那应该就对网格没什么影响，可能是Pointwise这个软件里面才有的设置，谢谢李老师

如图所示

在使用Pointwise网格绘制软件绘制完OpenFOAM网格，并设置好Boundary Conditions后，下面有一个Set Volume Conditions，里面可以选中体网格并设置图中所示的条件，之前使用的时候不清楚，一直是随便选的最下面的那个Volume Condition，因为自己对OpenFOAM的Volume Conditions并不清楚，所以所以想在论坛上请教一下大家：OpenFOAM网格里面的Volume Conditions分别代表什么意思？适用于什么情况？是否会影响计算的收敛？对计算结果有什么影响？
非常感谢！！！

@cccrrryyy 嗯嗯是的，symmetryPlane要求对称面的normal必须一样，有些时候就算差的很小不影响计算，但是它也会报错的，我后面修改成symmetry后，发现并行计算也可以进行了，问题应该就是在symmetry和symmetryPlane的区别上。谢谢您～

刚刚我发现原因了，对称面不能命名为symmetryPlane,不然在运行的时候会报错：

inconsistent patch and patchField types for 
    patch type symmetry and patchField type symmetryPlane

我把对称面改称sym后就可以开始计算了！

我把它改成3之后还是一样的报错，后面我是了一下不并行进行计算，它显示的报错是
：

--> FOAM FATAL IO ERROR: 
inconsistent patch and patchField types for 
    patch type symmetry and patchField type symmetryPlane

file: /home/jifeng/OpenFOAM/jifeng-6/run/Isolater-Huang/Isolater-Huang_1/0/alphat.boundaryField.symmetryPlane from line 59 to line 59.

    From function static Foam::tmp<Foam::fvPatchField<Type> > Foam::fvPatchField<Type>::New(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&, const Foam::dictionary&) [with Type = double]
    in file /home/jifeng/OpenFOAM/OpenFOAM-6/src/finiteVolume/lnInclude/fvPatchFieldNew.C at line 160.

但是我的alphat以及其他的初始条件以及polyMesh文件夹里面的boundary里面的symmetryPlane边界面都是用的是symmetry，不知道为什么它会说alphat里面会

inconsistent patch and patchField types for 
    patch type symmetry and patchField type symmetryPlane

因为模型的网格量比较大，所以我选择并行计算，但在过程中遇到了一些问题，之前也没有遇到过，在论坛上向大家请教一下为什么这样？应该怎么进行修正？谢谢！

decomposeParDict

numberOfSubdomains  3;

/*
    Main methods are:
    1) Geometric: "simple"; "hierarchical", with ordered sorting, e.g. xyz, yxz
    2) Scotch: "scotch", when running in serial; "ptscotch", running in parallel
*/

method              scotch;

simpleCoeffs
{
    n               (4 2 1); // total must match numberOfSubdomains
    delta           0.001;
}

hierarchicalCoeffs
{
    n               (4 2 1); // total must match numberOfSubdomains
    delta           0.001;
    order           xyz;
}

运行的代码

mpirun -np 2 rhoSimpleFoam parallel

报错

--> FOAM FATAL ERROR: 
Wrong number of arguments, expected 0 found 1

--> FOAM FATAL ERROR: 
Wrong number of arguments, expected 0 found 1


FOAM exiting


FOAM exiting

-------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

  Process name: [[65460,1],1]
  Exit code:    1
--------------------------------------------------------------------------

不是的，后面我用komegaSST算了一下这个模型，吧压力求解器改成了PCG之后就收敛了，而且结果比较满意，希望对遇到同样问题的朋友也有参考价值

@东岳
恩恩，好的，我按照老师的思路再来来设置一下，谢谢李老师

@东岳
谢谢李老师，我检查了一下速度方向，是正确的，但mu是0是因为我参照的tutorials里面的给定方法，自己也不太清楚具体应该怎么修改，这个mu和之前算incom里面的nu是一个意思吗？
我刚刚在自己电脑上重新设置了一下，修改该了一些边界条件，用的是rhoCentralFoam,现在没有发散并且正在计算了，时间步长我修改了，可以计算也没有发散
P

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 101325;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 101325;
    }

    outlet
    {
        type            zeroGradient;
        phi             phi;
        rho             rho;
        psi             thermo:psi;
        U               U;
        gamma           1.4;
        inletOutlet     false;
        correctSupercritical true;
        lInf            0;
        fieldInf        101325;
        value           uniform 101325;
    }

    sym
    {
        type            symmetryPlane;
    }

    wall
    {
        type            zeroGradient;
    }
}

U

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (208.3 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (208.3 0 0);
    }

    outlet
    {
        type           zeroGradient;
        inletValue     uniform (0 0 0);
        value          $internalField;
    }

    sym
    {
        type            symmetryPlane;
    }

    wall
    {
        type            noSlip;
    }

T

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 300;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 300;
    }

    outlet
    {
        type            zeroGradient;
    }

    sym
    {
        type            symmetryPlane;
    }

    wall
    {
        type            zeroGradient;
    }

thermophysicalProperties

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture // air at room temperature (293 K)
{
    specie
    {
        molWeight   28.9;
    }
    thermodynamics
    {
        Cp          1005;
        Hf          0;
    }
    transport
    {
        mu          0;
        Pr          0.71;
    }
}

turbulenceProperties

simulationType laminar;

我下面想进行有粘可压计算，但我不知道0文件里面增加的alphat、k、nut、omega的初始条件和边界值具体应该给多少，请问李老师这个有什么可以参考的吗？