有算过水锤的吗

[shrine]

@李东岳 任务中是要算出压力波的传递。但是有这个减压井应该会小很多。看来我还是得试一下。
thermophysicalProperties.water文件中

equationOfState
    {
        C (0.001278 -2.1055e-06 3.9689e-09 4.3772e-13 -2.0225e-16);
    }

这个是啥？

还有就是这个求解器的tutorial只有层流，是不是只能算层流？

[胡同小巷]

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

[胡同小巷]

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

[李东岳]

@胡同小巷 http://dyfluid.com/rhoPimpleFoam.html 我用这个求解器算的 你可以看一下

[shrine]

@胡同小巷 程序还没跑起来，东岳老师那个是没有自由液面的，我的算例有自由液面，准备用compressibleInterFoam

[胡同小巷]

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

[李东岳]

水锤这面有参考的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 不涉密的话你清理下算例 发在论坛我给你调一下

[胡同小巷]

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