浅水方程出入口条件设置

[xiaolin]

求问怎么设置出口边界h值，每个单元不同，（因为h0不同）

[李东岳]

出口边界的h值是可以计算的，这样限定不会对方程体系构成问题么？

[xiaolin]

@东岳 原例子就是给定了出口的h和入口的hu，模型是出入口h0分别统一，当我设置h0不一致时程序就发散了

[李东岳]

我弄混了，我说的是hTotal，你试试 https://www.cfd-china.com/topic/3499/openfoam小代码 里面的codedFixedValue，功能类似

[xiaolin]

@东岳
李老师，怎么用到h0中的数据呢，报了这个错误

    outlet
    {
        type            codedFixedValue;
        value           uniform 0.01; //default value
        redirectType    linearTBC; //name of new BC type
        code
        #{
        const vectorField& Cf = patch().Cf(); // get face center coordinate;
        scalarField& h = *this;
        forAll(Cf,faceI)
        {
        	h[faceI] = 0.01-h0[faceI]; 
        }
        #};

[xiaolin]

@东岳 顶部#include已加