Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. pimplefoam在瞬态模拟过程中库朗数突然增大然后报错停止

pimplefoam在瞬态模拟过程中库朗数突然增大然后报错停止

已定时 已固定 已锁定 已移动 OpenFOAM
5 帖子 3 发布者 3.2k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • G 离线
    G 离线
    Guoguo
    写于 最后由 编辑
    #1

    最近在使用openFOAM对下肢血管进行血流动力学模拟,为了更准确所以在入口处使用了随时间变化的边界条件,但计算一段时间后,库朗数会突然增大,然后程序报错停止,有没有大神可以帮忙看看应该如何修改设置。
    版本:openfoam 10
    求解器:pimpleFoam
    U文件,入口速度方向是z轴的负方向

    boundaryField
    {
        WALL
        {
            type            noSlip;
        }
    
        INLET
        {
            type			uniformFixedValue;
    		uniformValue	
    		{
    			type		tableFile;
    			file		"v.txt";
    		}
    		value			(0 0 -0.032180);
        }
    
        OUTLET1
        {
            type            zeroGradient;   
        }
    
       OUTLET2
        {
            type            zeroGradient;   
        }
    

    controlDict文件

    application     pimpleFoam;
    
    startFrom       startTime;
    
    startTime       0;
    
    stopAt          endTime;
    
    endTime         0.7;
    
    deltaT          1e-5;
    
    writeControl    adjustableRunTime;
    
    writeInterval   0.01;
    
    purgeWrite      0;
    

    fvSolution文件

    solvers
    {
        p
        {
            solver           GAMG;
            tolerance        1e-7;
            relTol           0.1;
    	preconditioner  DIC;
            smoother        symGaussSeidel;
    	nCellsInCoarsestLevel    10;
    	mergeLevels      1;
        }
    
        pFinal
        {
            $p;
            relTol          0.1;
        }
    
        U
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-05;
            relTol          0.1;
        }
    
        UFinal
        {
            $U;
            relTol          0.1;
        }
    }
    

    fvSchemes文件

    ddtSchemes
    {
        default         Euler;
    }
    
    gradSchemes
    {
        default         Gauss linear;
    }
    
    divSchemes
    {
        default             none;
    
        div(phi,U)          Gauss linearUpwind grad(U);
    
        div((nuEff*dev2(T(grad(U))))) Gauss linear;
    }
    
    laplacianSchemes
    {
        default         Gauss linear corrected;
    }
    
    interpolationSchemes
    {
        default         linear;
    }
    
    snGradSchemes
    {
        default         corrected;
    }
    
    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 编辑
    #2

    把你的时间步长调小一点

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    G 1 条回复 最后回复
  • G 离线
    G 离线
    Guoguo
    在 中回复了 李东岳 最后由 编辑
    #3

    @李东岳 李老师你好,我将时间步长调小后依然会出现库朗数增大的情况,但是我发现每次都在迭代到0.23s左右出现这种情况,这是否与我入口速度的边界条件有关?
    Figure_1.png
    这是我入口速度的波形图

    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 编辑
    #4

    pimplefoam不至于如此,估计还是边界条件不太稳定。另外随时间变化的入口,最有可能发散应该在你那个0.08秒左右而不是0.23秒。不知道网格质量如何

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 学流体的小明学 离线
    学流体的小明学 离线
    学流体的小明 神
    写于 最后由 编辑
    #5

    fvSolution里面的pFinal,relTol应该调到0吧。
    矩阵最终求解的这些*Final一般都强制要求满足容差。
    应该是这个原因导致你某一步没收敛,或者误差逐渐积累,导致发散。

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]