CFD中文网

@李东岳 老师您好，我这个水击可以计算了，但是感觉结果不对阿，和我用fluent计算的结果差很多的。
第一个是阀门处的压力随时间变化
第二个是残差图
希望老师可以抽点时间帮看一下
file:///home/huliren/%E5%9B%BE%E7%89%87/2021-08-23%2020-19-33%20%E7%9A%84%E5%B1%8F%E5%B9%95%E6%88%AA%E5%9B%BE.png
file:///home/huliren/%E5%9B%BE%E7%89%87/2021-08-23%2020-19-25%20%E7%9A%84%E5%B1%8F%E5%B9%95%E6%88%AA%E5%9B%BE.png

@李东岳 水击，运动粘度0.00003967，进口压力500000，出口速度经历0.4s由0.12降至0，求解0--0.4s内的出口处压力随时间变化曲线
[链接文本](8076956a-b043-4046-8cf6-9b11d520d44e-shuiji.zip 链接地址)

@胡同小巷 老师您好，我把这个文件发到下面了，在您邮箱也发了一份。望老师帮忙调整一下。

@李东岳 老师您好，我用rhoPimpleFoam计算水击但是运行不了，不知道什么错误
p文件

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

internalField   uniform 5e+05;

boundaryField
{
    fixedWalls
    {
        type            zeroGradient;
    }

    inlet
    {
        type            fixedValue;
        value           uniform 5e+05;
    }

    outlet
    {
        type            zeroGradient;
    }

    frontAndBack
    {
        type            empty;
    }

    
}

U文件

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

internalField   uniform (0.12 0 0);

boundaryField
{
    fixedWalls
    {
        type            noSlip;
    }

    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            uniformFixedValue;
        uniformValue    scale;
        scale
        {
        type            linearRamp;
        start           0;
        duration        0.4;
        }
        value           2;
        
    }

    frontAndBack
    {
        type            empty;
    }

}

T文件

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 20;

boundaryField
{
    fixedWalls
    {
        type            fixedValue;
        value           $internalField;
    }

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            fixedValue;
        value           $internalField;
    }

    frontAndBack
    {
        type            empty;
    }
}

thermophysicalProperties

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          eConst;
    equationOfState rPolynomial;
    specie          specie;
    energy          sensibleInternalEnergy;
}

mixture
{
    specie
    {
        molWeight       100;
    }
    equationOfState
    {
        C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16);
    }
    thermodynamics
    {
        Cv              4185;
        Hf              0;
    }
    transport
    {
        mu              3.967e-05;
        Pr              1;
    }
}

controlDict

application     rhoPimpleFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         0.4;

deltaT          0.00001;

writeControl    timeStep;

writeInterval   100;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

下面是报错：

Create time

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calculations will do 2 corrections


PIMPLE: Operating solver in transient mode with 2 outer correctors


Reading thermophysical properties

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          eConst;
    equationOfState rPolynomial;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U



--> FOAM FATAL IO ERROR: 
Expected a '(' while reading VectorSpace<Form, Cmpt, Ncmpts>, found on line 43 the label 2

file: /home/huliren/OpenFOAM/huliren-8/run/SHUJI/shuiji/0/U/boundaryField/outlet/value at line 43.

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

FOAM exiting

@李东岳 老师codedFixedValue不是设置非均匀边界的吗？我这个问题是均匀边界，但是是随时间变化的。
codedFixedValue可以吗？

@李东岳
又来请教问题了。
我这是水击问题，开始管道速度为0.1209，我不想让出口瞬间变为0，想让出口设置为慢慢的变为0，比如经历2s，以一定的加速度慢慢的变为0。那出口速度这里应该怎么写？
谢谢老师

@西湖冷月 我理解了。我是做的水锤问题。
我现实中的问题是，从关闭阀门到阀门完全关闭耗时4s。开始流速2m/s，管子内径25mm，管子长度150mm，一共画了1.5w个网格。
还有下面两个问题：
1、由于我画的是个长方形流动区域，我看bilibili上它设置的出口瞬间关闭(0.000001s)
可是我想让出口均匀的关闭，以加速度为-0.5m/s2关闭。不知道能不能这样设置。
2、关于计算参数设置：计算10000步，每步的时间把原来的0.000001s设置成0.0004s，即(10000*0.0004=4s)，这样是不是就代表算出来的是实际物理时间的4s的过程？但是不知道这样是否对，是否可行？
望指点

@李东岳 好的，感谢老师，期待老师

@李东岳 report-water-hammar.pdf

@shrine 是不是环境配置的问题，可以留个QQ或微信探讨吗？有时候不是很方便

@李东岳 东岳老师可以把您之前算得水击的openfoam文件发一下吗？研究探究。

@李东岳 从2020版开始fluent可以汉化了

不太理解这个时间步数量和时间步长、最大迭代步数的关系。
我开始这样理解的：每步最多可以迭代20次，然后每次的时间是0.000001s，我设置的是10000步，这么算的话10000*0.000001还不到1s就算完了。
好像不太对呀。
希望可以指点一下

@李东岳 好的，谢谢老师，我研究一下

@李东岳 东岳老师，您之前用单相流做水锤是用的哪个求解器？可以指导一下吗？

有没有做过水锤的老师，自己在流体力学书上找到了连续性方程和运动方程，在openfoam中也找到了可压缩求解器，但是看不懂求解器啥意思

@shrine 您好，水锤计算的怎么样了？我导师也让我在openfoam下做水锤