SRF求解器中如何用fvOptions给定恒Urel？

东方白杨

各位大佬，最近想在旋转状态下算一个流动，这个流动也是通过fvOptions给定压力梯度来达到指定Urel的。在静止状态下直接使用fvOptions文件指定就行了，但是在旋转状态下因为U到Urel不同，fvOptions如何指定恒Urel呢？而不是U？

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

momentumSource
{
    type            meanVelocityForce;

    selectionMode   all;

    fields          (Urel);
    Ubar            (1.3 0 0);
}

东方白杨

@ 李东岳老师

李东岳

你那个代码，应该可以直接施加到Urel上，不过我不太确定，从代码的原理来看没问题，但我没实操过，你试试

东方白杨

不太行，直接改成Urel运行出来速度都3000多了，很离谱。。。。感觉只能改方程加手动压力梯度吧

李东岳

那就是不能这么加Urel源项，代码起作用了，但是算法不对

东方白杨

@李东岳排除了一下问题，在静止状态下rpm=0用SRF求解器，写Fvoptions（Urel）是可以正常算的。问题主要出现在周期性边界条件这里，当我进出口不用周期性边条Urel就是正常的，一旦用了周期性边条，旋转状态下Urel就不对劲了。

李东岳

当我进出口不用周期性边条Urel就是正常的，一旦用了周期性边条，旋转状态下Urel就不对劲了。

这个看起来是openfoam的bug了。

东方白杨

@李东岳老师我又想了一下，感觉旋转状态下用cyclic周期性边条有问题，因为如果是一个直通道绕轴旋转的话，进出口都不在一个旋转半径上，存在离心力加速，如果此时用周期性边条的话，算出来显示通道内中心点速度Urel越来越大了。

CFD中文网

SRF求解器中如何用fvOptions给定恒Urel？