群组信息 私有

administrators

管理员

  • RE: 算例初始条件设置求助!

    其中气相匀速0.2m/s

    interFoam不能只设置气相匀速0.2m/s,因为两相公用一个速度。

    发布在 OpenFOAM
  • RE: 算例初始条件设置求助!

    你用的什么求解器

    发布在 OpenFOAM
  • RE: paraview显示面场?

    um. I dont understand. 我这是个2为网格,每个网格单元应该有4个内部面面,但是两个empty边界面上却具有通量?

    发布在 OpenFOAM
  • RE: paraview显示面场?
    void Foam::QMOM::test()
    {
    	const labelUList& owner = mesh_.owner();
    	const labelUList& neighbour = mesh_.neighbour();
    
    	forAll(M_[0], cellI)
    	{
    		forAll(mesh_.cells()[cellI], faceI)   
    		{
    			if (cellI == 1897)
    			{
    				Info << "face" << faceI << nl;
    				label faceIndex = mesh_.cells()[cellI][faceI];
    				//Info << "mFlux[" << faceIndex << "] = " 
    				//		<< mFlux_[faceIndex].value() << nl;
    				forAll(mFlux_, kth)
    				{
    					Info << "mFlux[" << kth << "][" 
    						<< faceIndex << "] = " 
    						<< mFlux_[kth][faceIndex] << nl;
    				}
    			}
    		}
    	}
    }
    

    一个显示某网格(如第1897个网格点)周围的面(如6个面)的通量的代码。

    如下:

    PIMPLE: iteration 1
    MOMsolve() 
    momentsAdvection() function
    updateAbs function
    Wheeler function
    updateAbs function end
    KBFVM() FUNCTION 
    face0
    mFlux[0][3783] = -1.2927335654110290391e-12
    mFlux[1][3783] = -6.2675976060777894205e-12
    mFlux[2][3783] = -3.1150947399844006554e-11
    mFlux[3][3783] = -1.5863902768703573911e-10
    mFlux[4][3783] = -8.2750268353901142759e-10
    mFlux[5][3783] = -4.4202305898578507575e-09
    mFlux[6][3783] = -2.417473810022032185e-08
    mFlux[7][3783] = -1.3535140721940442922e-07
    face1
    mFlux[0][3784] = 2.0540870520051721362e-13
    mFlux[1][3784] = 9.4661825056687303207e-13
    mFlux[2][3784] = 4.4532871356728447274e-12
    mFlux[3][3784] = 2.1369999282078349297e-11
    mFlux[4][3784] = 1.0452769480664259678e-10
    mFlux[5][3784] = 5.2076058382076055416e-10
    mFlux[6][3784] = 2.6403668330037581177e-09
    mFlux[7][3784] = 1.3610985929982985127e-08
    face2
    mFlux[0][11360] = -0.38831761037105733125
    mFlux[1][11360] = 3.2471357929156595401e-12
    mFlux[2][11360] = 0.055793628070796436047
    mFlux[3][11360] = 8.30049086908404522e-11
    mFlux[4][11360] = 1.4292276668538757853e-05
    mFlux[5][11360] = 7.1461383336000636559e-05
    mFlux[6][11360] = 0.00036534632227100779416
    mFlux[7][11360] = 0.0019098593165646248534
    face3
    mFlux[0][16315] = 5.2242945423600630761
    mFlux[1][16315] = 1.101005821401192256e-10
    mFlux[2][16315] = 0.31179646759308882187
    mFlux[3][16315] = 2.812383025709799438e-09
    mFlux[4][16315] = 1.4292279304637762273e-05
    mFlux[5][16315] = 7.1433737066938304697e-05
    mFlux[6][16315] = 0.00036507857573279831271
    mFlux[7][16315] = 0.0019078661663797318558
    face4
    mFlux[0][3604] = 6.9764003936297054674e-12
    mFlux[1][3604] = 3.4017174538427634926e-11
    mFlux[2][3604] = 1.6992537663665707908e-10
    mFlux[3][3604] = 8.693051567457588884e-10
    mFlux[4][3604] = 4.5534209594944540047e-09
    mFlux[5][3604] = 2.4416268967863339806e-08
    mFlux[6][3604] = 1.340108550200319568e-07
    mFlux[7][3604] = 7.5279189328388737134e-07
    face5
    mFlux[0][3781] = -9.33250737219599897e-12
    mFlux[1][3781] = -4.5911180782797612225e-11
    mFlux[2][3781] = -2.3141179782345137439e-10
    mFlux[3][3781] = -1.1946980508076885345e-09
    mFlux[4][3781] = -6.3157644055822882014e-09
    mFlux[5][3781] = -3.4182357847142154226e-08
    mFlux[6][3781] = -1.8937287564175264068e-07
    mFlux[7][3781] = -1.0737815840220180914e-06
    kth is 0
    gMin is 4.8282703298717213403e-05
    kth is 1
    gMin is 4.8215082927350118335e-05
    kth is 2
    gMin is 4.8203426534192527083e-05
    kth is 3
    gMin is 4.8197530295908814742e-05
    deltaT is 5.000000000000000409e-06
    SetdeltaT  END
    KBFVM() FUNCTION END 
    momentsAdvection() function END
    updateAbs function
    Wheeler function
    1897's b[i] < 0
    m = 8.9981758148512276923e-06
    m = 4.1459899559787042548e-05
    m = 0.00019500283102905703906
    m = 0.00093553135820788324833
    m = 0.0045747020117296623154
    m = 0.022783875013194655745
    m = 0.11547491828732039942
    m = 0.59499954912512731386
    neigh cells are 4(1807 1896 1898 1987)
    
    
    --> FOAM FATAL ERROR: 
    
    
        From function b < 0
        in file momentsInversion/wheeler/wheeler.C at line 91.
    
    FOAM exiting
    
    
    

    我这一个月,一直在和1897网格单元作斗争

    发布在 OpenFOAM
  • RE: owner neighbour 网格寻址

    注意在变量赋值的时候要使用faceIndex而不是faceI,因为p[faceI]只是一直在1 2 3 4 5…中循环,p[faceIndex]是真正的面址

    发布在 OpenFOAM
  • RE: 有限体积法的强弱形式

    不是很懂有限元。
    据我了解有限差分和有限体积你可以选择任意一种形式的方程去离散。主要看你写的代码来自的方程的形式。

    发布在 基本理论
  • RE: 最近国内有没有组织OpenFOam的培训?

    @cq
    嗯,可以。

    市面上还有一些其他家的,不推荐。曾经有人参加过别的家的,然后又参加我的,他说特别水。连讲师做什么的都不知道。

    发布在 OpenFOAM
  • RE: outlet设为速度出口,压力为固定值结果很奇怪

    你用的什么进口条件?你的设置符合常规设置么?see Ferziger and Peric page 315

    发布在 OpenFOAM
  • RE: 算例初始条件设置求助!

    匀速0.2m/s

    0/U.air:

    internalField uniform (0.2 0 0);
    
    boundaryFields
    
    发布在 OpenFOAM
  • RE: 粘度模型,求解随温度变化的粘度?

    这都是非牛顿流体的本构方程,实际上非牛顿流体的书上都有的,我建议你参考Non-Newtonian Fluids: An Introduction。有问题再回帖 :big_mouth:

    发布在 OpenFOAM

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待