LaTex用的人也挺多的
OF九月九
帖子
-
-
感觉自带的跟好用,键盘输入,我用mathtype编号
-
好吧,既然没有人回答,我就把我的笨办法贴出来
foamToVTK-->paraview-->打开算例文件和粒子Cloud文件-->File/Save Data-->得到每个时间步每个粒子信息.csv-->用python提取有用信息 -
如题,在一个lagrangian求解器中,设置一个面的粒子碰撞为escape
patchInteractionModel standardWallInteraction; standardWallInteractionCoeffs { type escape; }
那么怎么得到粒子escape时的速度?用ParaView能实现吗?或者用cloudFunctions?求知情人士告知,不甚感谢!
-
我觉得这本书不错Multiphase Flows with Droplets and Particles (Second Edition)-Clayton T. Crowe,网上有电子版,里面介绍离子的受力运动和粒子传热相变比较详细。
-
有没有用LES算一下,结果怎么样?
-
把我们办公室的另两位拉下水了,终于不是我一个苦逼的在学of了。
-
@chpjz0391 谢谢!应该怎么改进?增加kinematicCloudPositions吗?还是改变nParticle?
文章 http://xueshu.baidu.com/s?wd=02Sea+surface+drag+and+the+role+of+spray&rsv_bp=0&tn=SE_baiduxueshu_c1gjeupa&rsv_spt=3&ie=utf-8&f=8&rsv_sug2=1&sc_f_para=sc_tasktype%3D{firstSimpleSearch}&rsv_n=2 -
首先,谢谢您的回复!的确,算例是不考虑颗粒之间的碰撞。按您说的修改成
这样算正确吗?计算结果质量分数太小了,总觉得和论文中的对不上。论文中是这样描述的:难道是我的颗粒数50000设的太少了吗?而且如果我想要改变颗粒的质量分数,直接改变颗粒的密度就可以吗?谢谢!
-
在一个lagrangian算例中,kinematicCloudProperties如下
solution { active true; coupled true; //false; transient yes; cellValueSourceCorrection on; interpolationSchemes { rho cell; U cellPoint; mu cell; } integrationSchemes { U Euler; } sourceTerms { schemes { U semiImplicit 0.5; } } } constantProperties { parcelTypeId 1; rhoMin 1e-15; //parcel的最小密度 minParticleMass 1e-15; //parcel的最小质量 rho0 1250; //单个颗粒的密度 youngsModulus 6e8; poissonsRatio 0.35; constantVolume true; } subModels { particleForces { sphereDrag; gravity; } injectionModels { model1 { type manualInjection; massTotal 0; parcelBasisType fixed; nParticle 1; //每个parcel的颗粒数,设成1就是DEM SOI 0; positionsFile "kinematicCloudPositions"; U0 ( 0 0 0 ); sizeDistribution { type fixedValue; fixedValueDistribution { value 0.00125; } } } } dispersionModel none; //stochasticDispersionLES; patchInteractionModel standardWallInteraction; heatTransferModel none; surfaceFilmModel none; collisionModel none; stochasticCollisionModel none; radiation off; standardWallInteractionCoeffs { //type stick; type rebound; e 1; mu 0; } } cloudFunctions { }
kinematicCloudPositions如下:
50000 ( ...... )
计算域为35m×35m×5m的长方体。Euler相密度为1.205kg/m3,则粒子的质量负载为
请问这样计算正确吗?谢谢?
@散漫守望2016 -
@cfd-china 谢谢!我找到了,好像是主文件中的parcel.evolve()。
-
在lagrangian求解器中只能看到连续相的方程:
tmp<fvVectorMatrix> UEqn ( fvm::ddt(rho, U) + fvm::div(phi, U) + MRF.DDt(rho, U) + turbulence->divDevRhoReff(U) == rho.dimensionedInternalField()*g + parcels.SU(U) + fvOptions(rho, U) );
请问离散相方程怎么找到?谢谢!
-
@OF九月九 好吧,有个sphereToCell
-
@wwzhao 你好,如果我想设置一个球形区域,应该在setFields中怎么设置?
-
@wwzhao 明白了,我思维定势了,谢谢!
-
在两相流中,如果设置成不同高度的粒子浓度不同,我在setFields是这样设置的
defaultFieldValues ( volScalarFieldValue alpha.air 1 volScalarFieldValue alpha.particles 0 ); regions ( boxToCell { box (0 0 0) (35 35 1); fieldValues ( volScalarFieldValue alpha.air 0.8 volScalarFieldValue alpha.particles 0.2 ); box (0 0 1) (35 35 3); fieldValues ( volScalarFieldValue alpha.air 0.9 volScalarFieldValue alpha.particles 0.1 ); } );
即设置了三层不同浓度的粒子,但是setFileds后
Reading setFieldsDict Setting field default values Setting internal values of volScalarField alpha.air Setting internal values of volScalarField alpha.particles Setting field region values Adding cells with center within boxes 1((0 0 1) (35 35 3)) Setting internal values of volScalarField alpha.air Setting internal values of volScalarField alpha.particles End
只默认设置了第二层的浓度
这是怎么回事?是不是这样设置不行?应该怎么办?求教? -
@cfd-china :happy: THANK YOU!
-
如题所述,如果我运行
mpirun -np 4 ***Foam -parallel
在处理器设置为1或4的时候都没有出现错误。是否算例中processor的数目与虚拟机处理器的数量设置要一致?与主机的处理器数量又有什么关系?主机的处理器数量是否会限制虚拟机处理器数量的设置?求教!
-
@队长别开枪 果然如此,明白了,非常感谢!!!:cheeky:
-
感谢两位的回答,但是controlDict不是算例里面的控制文件吗?为什么用Tecplot读取它。而且我用了一个算例读取出现错误
又是怎么回事呢? -
@队长别开枪 新手一枚,你能说的具体一点吗,我用的是windows版的tecplot360 2015版,怎么读取ubuntu下openfoam的计算结果?非常感谢!!!
-
http://blog.sina.com.cn/s/blog_4a0a8b5d01014rmp.html 按这篇博客上的方法编译不出来,是不是版本问题?求指教。谢谢!
-
我知道是这个,我以为.A()和.H()包含什么重要东西呢。谢了!:happy:
-
大家好,在icoFoam解析中
volScalarField rAU(1.0/UEqn.A()); volVectorField HbyA("HbyA", U); HbyA = rAU*UEqn.H();//公式(15)
UEqn.A()、UEqn.H()是什么意思?
-
@cfd-china 你在羡慕别人的同时,别人也在羡慕你,凭LZ的条件,为什么不继续你的科研呢?博士毕业去重点高校当个教授,工资虽然没你同学高,但待遇好,日子过得肯定比他好,本来走的就不是同一条路,不能完全用工资去比。
-
如上式,添加方程的最后两项在OpenFOAM中如何实现?
Word里面你们怎么打公式
Word里面你们怎么打公式
Lagrangian怎么求得面上的粒子速度?
Lagrangian怎么求得面上的粒子速度?
DEM和DPM
有没有较好的讲解DES和LES的中文资料?
如何在tecplot中后处理OpenFoam计算结果
OpenFOAM教程 初级入门建议(2020年更新版)
lagrangian两相流中粒子质量负载的计算。
lagrangian两相流中粒子质量负载的计算。
lagrangian两相流中粒子质量负载的计算。
请教大家lagrangian离散相方程问题
请教大家lagrangian离散相方程问题
如果要在垂直方向上粒子浓度分层,setFields应该怎么设置?
如果要在垂直方向上粒子浓度分层,setFields应该怎么设置?
如果要在垂直方向上粒子浓度分层,setFields应该怎么设置?
如果要在垂直方向上粒子浓度分层,setFields应该怎么设置?
在虚拟机上用OpenFOAM并行运算应如何设置处理器数量?
在虚拟机上用OpenFOAM并行运算应如何设置处理器数量?
有谁知道怎么在openfoam3.0.1上编译foamToTecplot360?求教。
有谁知道怎么在openfoam3.0.1上编译foamToTecplot360?求教。
有谁知道怎么在openfoam3.0.1上编译foamToTecplot360?求教。
有谁知道怎么在openfoam3.0.1上编译foamToTecplot360?求教。
icoFoam求解器解析中问题。
icoFoam求解器解析中问题。
凭什么我辛苦二十年,现在却比别人差那么多(看化学和计算机专业的区别)
LES的N-S方程中加入科氏力项和重力项在OpenFOAM中如何实现?