ダクト内間欠流の計算2

昨日の続きです。

    inlet
    {
    type             groovyBC;
    variables       "p0=0.5;"; // =(0.5*U0*U0)/4
    valueExpression "p0*(pow(sin(pi*time()/50),7))";    
    value        uniform 0;
    }
    outlet
    {
    type            groovyBC;
    variables       "p0=0.5;"; // =(0.5*U0*U0)/4
    valueExpression "p0*(-pow(sin(pi*time()/50),7))";
    value            uniform 0;
    }

という感じで、圧力差がついた後に、しばらく圧力差がない時間帯を設けてみた。

その結果がこちら。圧力差は対称なはずなんだが、結果は非対称で不思議な結果です。