@李东岳 从另外一个思路来做,能不能指定内场中一区域的速度分布满足某一个函数关系式呢?
液固两相流
帖子
-
-
@yhdthu 前辈能够分享一下解决的经验吗?
-
二维计算域示意图如上图所示,请问大家在openFoam中能否指定计算域中的某一块区域的涡量为某一值吗? -
@fubianhanshu 需要调用particleTrack函数哟
-
有没有大佬可以提点意见呀!
-
背景:在使用欧拉-拉格朗日算法求解固液两相流动时,使用颗粒追踪的后处理函数,如图所示:
在结果中将会输出颗粒的属性:
现在我希望增加输出项,简单的如颗粒的体积,我修改了kinematicParcelIO.C中的writefields函数:添加了如下代码:
但是计算结果中还是没有颗粒体积的输出?
请问大家有碰到过类似的问题吗?
谢谢各位了! -
@bubbleEL 谢谢你
-
计算区域为旋转区域与非旋转区域,颗粒由非旋转区域入射,但是后期处理颗粒轨迹时候发现,颗粒无法通过动静交界面(cyclicAMI).
在颗粒求解文件中设置交界面为none
请问大家有碰到过这问题吗?有什么好的解决方法吗? -
@bubbleEL 我目前也碰到了颗粒无法进入旋转区域的问题,我的交界面设置的也时none,请问楼主是进行了怎样的修改,使得颗粒能够进入旋转区域呢?
-
请问楼主这个问题解决了吗?
-
@东岳
我的网格有多个交界面(用的是cyclicAMI),交界面对于拉格朗日的计算会有影响吗? -
目前我在采用稳态的欧拉拉格朗日求解器,求解水轮机的磨损问题,颗粒数为50000,网格数为1100w左右。颗粒文件的设置如下:```
code_textsolution { 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 手动划分网格吗?
-
-
反馈以下:
p的初始条件,进口设置为totalPressure
u的初始条件,进口可设置为pressureinletuniformvelocity -
已成功植入emoij
-
基于MRF模型计算液固两相流动时,在旋转区域,颗粒会受到科式力的作用,在openfoam中颗粒受力并没有该项。
请教各位大佬应该如何把这个力加上去呀? -
背景:在水轮机的流场计算中,通常是水头固定的,因此需要在水头已知的情况下计算流量。
请教以下大家在openfoam中能够怎样实现呢? -
@宝丁 请问这个边界条件中指定的压力值是总压吗?
-
压力梯度力的计算公式如下图1所示:
图1
虚拟质量力的计算公式如下图2所示:
图2
在openFOAM中的计算代码如下所示:
压力梯度力:
value.Su() = mass*td.rhoc()/p.rho()*DUcDt;
虚拟质量里:
value.Su() = Cvm*PressureGradient
OpenFOAM里面的计算公式与实际的计算公式有一定的差别,请教一下大家这种处理是基于什么考虑呢? -
-
@lllwonderliquid 请问你的问题目前解决了吗?
我目前也在使用MRF以及AMI边界条件进行计算,但是碰到了与你类似的问题 -
@sjlouie91
请教一下您,应该怎样定义AMI边界条件呢?可以直接在boundary中定义吗? -
@东岳 固体密度是2650kg/m3 液体是1000kg/m3,没有轨迹交叉现象,颗粒的体积分数小于1%
-
@Ant 由于我计算的工况,颗粒浓度很小,所以我没有考虑相分数,以及颗粒之间的碰撞
-
@东岳 目前我的认知是,在单相耦合的情况下,对于颗粒的受力计算可能存在一些问题,在openfoam中,将颗粒简化为一个点,如果颗粒大于网格尺寸,可能对于颗粒受力的计算存在较大的误差,虽然颗粒小于网格的话貌似这种简化也会带来误差!
@刘雄国 您所说的bug指的是哪一块啊? -
如果在不考虑体积分数的情况下,颗粒粒径大于网格单元,会带来哪些影响呢?
基于DPMFoam自定义求解器,取消了对于体积分数的求解,颗粒粒径为0.45mm,计算模型为弯管(边界层网格0.2mm),入口采用patchinjection,颗粒数为10000,从inlet入射,结果进行颗粒求解时,50%的颗粒直接由inlet escape?
-
如果不考虑相分数的情况下,颗粒的粒径可以大于网格尺寸吗?@东岳
求助帖—OpenFOAM中初始场的指定
关于函数初场设定
求助帖—OpenFOAM中初始场的指定
关于颗粒属性的输出
关于颗粒属性的输出
关于颗粒属性的输出
有没有人用颗粒求解器计算过带AMI面的问题的?
液固两相流动+滑移网格颗粒运动求解的问题
有没有人用颗粒求解器计算过带AMI面的问题的?
有没有人用颗粒求解器计算过带AMI面的问题的?
关于稳态欧拉-拉格朗日计算速度的探讨
关于稳态欧拉-拉格朗日计算速度的探讨
关于OF并行分区的问题
关于OF并行分区的问题
关于OF并行分区的问题
OpenFoam中如何设置总压进出口
多重坐标系下如何考虑颗粒在旋转区域的科式力?
多重坐标系下如何考虑颗粒在旋转区域的科式力?
OpenFoam中如何设置总压进出口
如何在给定总压入口的同时给定速度的方向
openFOAM中颗粒压力梯度力于虚拟质量力的计算问题
超算上采用openfoam计算时,mpi报错
对旋转机械进行计算时的错误提示
交界面处拓扑结构不同,如何定义边界条件
请问DPMFoam不能用于非常精细的网格的模型么?
请问DPMFoam不能用于非常精细的网格的模型么?
请问DPMFoam不能用于非常精细的网格的模型么?
不考虑体积分数的情况下,颗粒粒径能否大于网格单元
请问DPMFoam不能用于非常精细的网格的模型么?