CFD中文网

@李东岳 在 关于div的小问题 中说：

你要求的那一项方程形式是怎么写的

文献中给的方程是这样子的，下面这张图是我现在采用的方式，不知道对不对 （其中DIVUP是速度的散度）

$ \nabla\cdot(Ue)=U\cdot\nabla e + e \nabla \cdot U$ ，$\nabla \cdot U \neq 0$

关于这一项$ \nabla\cdot(Ue)$在openfoam中是直接写成div(phi,e)呢？还是div(phi,e)+ediv(U)呢？

目前测试出来差别比较大

@李东岳 从另外一个思路来做，能不能指定内场中一区域的速度分布满足某一个函数关系式呢？

@yhdthu 前辈能够分享一下解决的经验吗？

二维计算域示意图如上图所示，请问大家在openFoam中能否指定计算域中的某一块区域的涡量为某一值吗？

@fubianhanshu 需要调用particleTrack函数哟

有没有大佬可以提点意见呀！

背景：在使用欧拉-拉格朗日算法求解固液两相流动时，使用颗粒追踪的后处理函数，如图所示：

在结果中将会输出颗粒的属性：

现在我希望增加输出项，简单的如颗粒的体积，我修改了kinematicParcelIO.C中的writefields函数：添加了如下代码：

但是计算结果中还是没有颗粒体积的输出？
请问大家有碰到过类似的问题吗？
谢谢各位了！

@bubbleEL 谢谢你

计算区域为旋转区域与非旋转区域，颗粒由非旋转区域入射，但是后期处理颗粒轨迹时候发现，颗粒无法通过动静交界面（cyclicAMI）.
在颗粒求解文件中设置交界面为none
请问大家有碰到过这问题吗？有什么好的解决方法吗？

@bubbleEL 我目前也碰到了颗粒无法进入旋转区域的问题，我的交界面设置的也时none，请问楼主是进行了怎样的修改，使得颗粒能够进入旋转区域呢？

请问楼主这个问题解决了吗?

@东岳
我的网格有多个交界面（用的是cyclicAMI），交界面对于拉格朗日的计算会有影响吗？

目前我在采用稳态的欧拉拉格朗日求解器，求解水轮机的磨损问题，颗粒数为50000，网格数为1100w左右。颗粒文件的设置如下：```
code_text

solution
{
    active          true;
    
    calcFrequency   15;

    maxTrackTime    6;

    maxCo           0.6;

    coupled         true;

    cellValueSourceCorrection off; //suggested to be off

    interpolationSchemes
    {
    
       Urel      	cellPoint;		//U for water
       rho     		cell;	 		//rho of water
       mu      		cell;			 //mu of water
       DUcDt        cell;	 		//DUcDt of water 

    }

    integrationSchemes
    {
        U               Euler;
    }

    sourceTerms
    {
     resetOnStartup false;
        schemes
        {
            U semiImplicit 1;
        }
    }
}

constantProperties
{
    rho0 2650;	//the density for the particle
}

subModels
{
    particleForces
    {
        sphereDrag; //drag forces

        gravity;  // gravity force

		pressureGradient //pressure gradient force
		{
		  U   Urel;
		}
		
		virtualMass    //virtual mass force
		{
		  Cvm   0.5;
		}
		
		MRF
		{
		 omega (0 0 -104.72);
		 centreOfRotation (0 0 0);
		 Ctemp	MRFC;
		}
	
     }

    injectionModels
    {

        model1
        {
		 type					patchInjection; 		
         parcelBasisType   		fixed;					
         patchName        		INLET_SV;					
         U0                		(0 0 0);				
		 nParticle          	420;
		 duration            	1;					
		 parcelsPerSecond		50000;
         flowRateProfile 		constant 1;
		 massFlowRate     		 0;       			
		 //positionsFile   		"kinematicCloudPositions";		
            sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   0.0002;//m
                }
            }
        }
    }

    dispersionModel stochasticDispersionRAS;

    patchInteractionModel UserLocalInteraction;
    //using the tabakoff rebound model 

    UserLocalInteractionCoeffs
    {
        patches
        (
            "(.*)"
            {
                type rebound;
				reboundModel 2;
            }

            OUTLET_SC
            {
                type none;
            }
            INLET_SV
            {
                type escape;
            }
            OUTLET_SV
            {
                type none;
            }
            INLET_GV
            {
                type none;
            }
            OUTLET_GV
            {
                type none;
            }
            INLET_RUN
            {
                type none;
            }
            OUTLET_RUN
            {
                type none;
            }
            INLET_DT
            {
                type none;
            }

            INLET_SC
            {
                type none;
               
            }
            OUTLET_DT
            {
                type escape;
               
            }
        );
    }

   
    heatTransferModel none;

    surfaceFilmModel none;

    collisionModel none;

    stochasticCollisionModel none;

    radiation off;
}

流场迭代15步之后，会进行颗粒的迭代，颗粒计算平均耗时5min。
但每次计算至210步时，颗粒就算不动了（前面的颗粒计算均是正常的）。
我如果把颗粒数降至10000，计算至750步，颗粒就算不动了。

请问一下大家又碰到这类问题吗？有什么比较好的解决办法吗？

@xpqiu 您能够介绍一下上述方法，单看您发的文件不太懂

@xpqiu 手动划分网格吗？

@东岳 在 关于OF并行分区的问题 中说：

https://www.openfoam.com/releases/openfoam-v1712/parallel.php

基金会版本的用不了

反馈以下：
p的初始条件，进口设置为totalPressure
u的初始条件，进口可设置为pressureinletuniformvelocity

已成功植入emoij

基于MRF模型计算液固两相流动时，在旋转区域，颗粒会受到科式力的作用，在openfoam中颗粒受力并没有该项。
请教各位大佬应该如何把这个力加上去呀？