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. 数组变量边界条件的设置

数组变量边界条件的设置

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

    单独计算和通过数组循环计算的结果分别如下(取相同方向):
    0_1520941821583_QQ截图20180313194929.png
    0_1520941836934_QQ截图20180313194947.png
    循环计算时在四个顶点处有数据错误。

    在计算各方向的方程时为了方便使用了数组进行循环,所设的变量如下:

        PtrList<volScalarField> PHI(nRay_);
        forAll(PHI, nRay_)
        {
    	PHI.set
    	(
    		nRay_,
    		new volScalarField
    		(
              		    IOobject
    			    (
    				"PHI"+i,
    				runTime.timeName(),
    				mesh,
    				IOobject::MUST_READ,
    				IOobject::AUTO_WRITE
    			    ),
    			    mesh
    		)
    	);
        } 
    

    计算方程如下:

        while (simple.loop())
        {      
              Info<< "Time = " << runTime.timeName() << nl << endl;
              i = 0; 
              forAll(PHI,i)
              {
                 while (simple.correctNonOrthogonal())
                 {  
                    fvScalarMatrix TEqn
                    (
                         fvm::ddt(PHI[i])
                       + fvc::div(phiA[i],PHI[i])
                       ==
                       - fvm::Sp(v*sigmaT,PHI[i])
                       + fvOptions(PHI[i])
                    );
                    TEqn.relax();
                    fvOptions.constrain(TEqn);
                    TEqn.solve() ;
                    fvOptions.correct(PHI[i]);
                 }
              }
    

    变量的边界条件设置如下:

    dimensions      [0 -1 0 0 0 0 0];
    
    internalField   uniform 0;
    
    boundaryField
    {
        movingWall1
        {
            type            empty;
        }
        movingWall2
        {
            type            zeroGradient;
        }
        fixedWall1
        {
            type            zeroGradient;
        }
        fixedWall2
        {
            type            empty;
        }
        fixedWall3
        {
            type            zeroGradient;
        }
        fixedWall4
        {
            type            zeroGradient;
        }    
    }
    

    单独调出每个变量的边界条件均与所设相同。求教这是什么问题?

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

    并没有看出什么错误 :quwan:

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yang山青
    在 中回复了 李东岳 最后由 编辑
    #3

    @李东岳 谢谢。我将原代码放入了另外一个求解器中就没有了这种问题,可能是路径的问题

    1 条回复 最后回复

  • 登录

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