关于两相流vof求解器边界条件



  • 大家好,我用interMixingFoam求解器求解水渠流动,即采用相对运动法设置水流从左向右运动,空气的速度总是不对,相分数和压力场感觉没有什么问题
    rr.png
    ss.png
    p.png

    水槽两侧和底面都是壁面边界条件,进口、出口、大气边界压力条件分别为:

        inlet
        {
            type            fixedFluxPressure;
            value           uniform 0;		
        }
    
        outlet
        {
    	type            zeroGradient;			
        }
    
        atmosphere
        {
            type            totalPressure;
            p0              uniform 0;
            U               U;
            phi             phi;
            rho             rho;
            psi             none;
            gamma           1;
            value           uniform 0;	
        }
    

    进口、出口、大气边界速度条件分别为:

        inlet
        {
            type            fixedValue;
            value           uniform ($velocity 0 0);
        }
    
        outlet
        {
            type            inletOutlet;
            inletValue      uniform (0 0 0);
            value           $internalField;
        }
    
        atmosphere
        {
            type            pressureInletOutletVelocity;
            value           $internalField;		
        }
    

    fvSolution文件:

    solvers
    {
        "alpha.*"
        {
            nAlphaCorr      1;
            nAlphaSubCycles 2;
            cAlpha          1;
    
            solver          smoothSolver;
            smoother        GaussSeidel;
            tolerance       1e-06;
            relTol          0;
            nSweeps         1;
        }
    
        pcorr
        {
            solver          PCG;
            preconditioner  DIC;
            tolerance       1e-10;
            relTol          0;
        }
    
        p_rgh
        {
            solver          PCG;
            preconditioner  DIC;
            tolerance       1e-07;
            relTol          0.05;
        }
    
        p_rghFinal
        {
            solver          PCG;
            preconditioner  DIC;
            tolerance       1e-07;
            relTol          0;
        }
    
        U
        {
            solver          smoothSolver;
            smoother        GaussSeidel;
            tolerance       1e-06;
            relTol          0;	
    		nSweeps         1;	
        }
    }
    
    PIMPLE
    {
        momentumPredictor no;
        nCorrectors     3;
        nNonOrthogonalCorrectors 0;
    }
    
    relaxationFactors
    {
        fields
        {
            p       0.3;
        }
        equations
        {
            ".*"           0.7;
        }
    }
    

    不知道哪里设置的不合理,或者是按照我设的边界条件,可能某些连续性条件不满足、不合理。
    谢谢大家!



  • @chengan-wang 你看一下相对速度的方向对不对?并且看起来好像是非结构网格,还比较糙



  • @东岳 谢谢,速度方向没有问题,网格就是用blockMesh画的结构化网格,刚开始算,网格没有画太多,比较稀疏


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2