twoLiquidMixingFoam中添加能量方程温度出现发散



  • 请教各位CFDers
    我打算在twoLiquidMixingFoam中添加能量方程
    \begin{equation}
    \frac{\partial T}{\partial t }+\nabla \cdot (\mathbf{U}T)= \nabla \cdot (a \nabla T)
    \end{equation}
    我在原始求解器中加入相应的代码如下:

    {
        fvScalarMatrix TEqn
        (
            fvm::ddt(T)
            + fvm::div(phi, T)
            - fvm::laplacian(DT, T)
        );
    
        TEqn.solve();
    }
    

    本质上就是利用twoLiquidMixingFoam中的流场数据计算温度场,温度场不影响流场。计算结果发现温度场出现发散,不知道原因是什么?
    希望大家能提供一点思路建议。谢谢



  • 这是一个非常简单的方程,我觉得可能是你的边界条件不对,不应该出现发散



  • @东岳 谢谢,我也觉得这个方程很简单,这段代码放到icoFoam里面没有出现过问题。但是放在两种液体混合求解器中就不行了。我感觉边界条件没有问题,流体从一侧流入,另一侧流出。

    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volScalarField;
        object      T;
    }
    dimensions      [0 0 0 1 0 0 0];
    
    internalField   uniform 293;
    
        inlet
        {
            type            groovyBC;
            valueExpression "(pos().y>1) ? (295):(293)"; 			
        }
    
        outlet
        {
            type            zeroGradient;	
        }
    
        freeSurface
        {
            type            zeroGradient;
        }
    
        bottom
        {
            type            fixedValue;
            value           293;			
        }
    
        front
        {
            type            symmetryPlane;	
        }
    
        back
        {
            type            symmetryPlane;		
        }
    
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volVectorField;
        location    "0";
        object      U;
    }
    
    dimensions      [0 1 -1 0 0 0 0];
    
    internalField   uniform (1 0 0); 
    
    boundaryField
    {
        inlet
        {
            type            fixedValue;
            value           uniform (1 0 0);
        }
    
        outlet
        {
            type            inletOutlet;
            inletValue      uniform (0 0 0);
            value           $internalField;		
        }
    
        freeSurface
        {
            type            pressureInletOutletParSlipVelocity;
            value           uniform (1 0 0);
        }
    
        bottom
        {
            type            fixedValue;
            value           uniform (0 0 0);			
        }
    
        front
        {
            type            symmetryPlane;		
        }
    
        back
        {
            type            symmetryPlane;	
        }
    
    }
    

Log in to reply
 

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