有算过水锤的吗

[李东岳]

水锤这面有参考的sci么，发过来瞅瞅？我做个标准算例

[shrine]

@李东岳 DOI 10.1007/s13369-013-0942-1

[李东岳]

几何有点复杂，没法画网格..

[shrine]

@李东岳 https://doi.org/10.1115/1.4000557

[胡同小巷]

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

[胡同小巷]

@李东岳 report-water-hammar.pdf

[李东岳]

@胡同小巷 很久很久以前了 当时没有备份 过几天我算一下这个 https://doi.org/10.1115/1.4000557

[胡同小巷]

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

[胡同小巷]

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

[李东岳]

@胡同小巷 用codedFixedValue 参考 https://www.cfd-china.com/topic/3499 里面有模板

[胡同小巷]

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

[胡同小巷]

@李东岳 老师您好，我用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

[李东岳]

@胡同小巷 类似这个 https://www.cfd-china.com/topic/5031/用rhopimplefoam计算t型管冷热流混合-计算几步就停掉/1 不涉密的话你清理下算例 发在论坛我给你调一下

[胡同小巷]

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

[胡同小巷]

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

[胡同小巷]

@李东岳 老师您好，我这个水击可以计算了，但是感觉结果不对阿，和我用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

[李东岳]

@胡同小巷 你那个算例我下载了之后没法运行。好像不是openfoam8的。不过你这个压力的图，看起来比较正常。