抱歉说漏了，rho、rho.water、rho.air这三个文件都要放进0文件夹。funkySetFields是能够使用函数表达式的，用来设置场内分布；边界上的值可以用funkySetBoundaryField来设置。在system下添加funkySetFieldsDict和funkySetBoundaryDict来使用。里面的内容，随便举个例子：

//funkySetFields
expressions
 (
    pressureWater
    {
     field p_rgh; //需要设置的物理量
     expression "100000 - 99800*pos().y";  //这行写你要的函数
     condition  "(pos().x <= 5) && (pos().x >= 0)"; //起作用的范围、可以注释掉，不影响使用
     keepPatches 1; //填1就好，要不然可能把边界条件改成zeroGradient
    }
 );

//funkySetBoundaryField
pressureWater
 {
    field p_rgh;
    expressions
    (
        {
            target value;
            patchName INLET; //选择边界面
            //variables "maxY=max(pts().y);";  //可以注释掉，不影响使用
            expression "100000 - 99800*pos().y"; //这行写你要的函数
        }
         //可以再expression下同时写多个边界的设置
        {
            target value;
            patchName OUTLET;
            //variables "maxY=max(pts().y);";
            expression "100000 - 99800*pos().y";
        }
    );
}

更多细节参考这里链接文本，以及swak4foam的tutorials
另外，你的这个算例里面压力也应当要设置，需要和密度相匹配。

CFD中文网

Jacobian