@李东岳 李老师好像对这个比较懂哈哈
VamosZhi
帖子
-
-
各位老板好:
在openfoam里面基于interPhaseChangeDyMFoam计算船用螺旋桨的空化问题。想要得到类似实验结果和其他人计算结果的非常清晰的空化尾流:
进口速度,旋转速度,根据cavitation number计算得到的vapor pressure是正确的。但是我的结果就是随着计算的进行,空化contour变得不明显,尤其到了计算后期:
我的网格:
我用的是SchnerrSauer模型,该部分我设的参数如下:
我了解到这些空话模型里面很多参数也是经验参数,不知道各位大佬有什么意见?就是想获得很清晰的空化轨迹。感谢!!!项目卡在这里4个多月了 -
请问在snappyhex中如何调整参数准确解析结构的边缘,比如下图圆柱端的直角边,以及翼型的尾缘。
结构抓取的就非常的难看。我看其他网站也有类似的介绍。
但是好像改了参数效果还是不好。各位有什么建议?目测应该不是继续加密的问题。 -
U#4.92Re100_v0.0001.zip 各位好,我现在想验证一篇论文,Re100, reduced order velocity 4.92,两个振动自由度上刚度相同,无机械阻尼。几篇论文显示最大相对振幅只有0.6。但是经过转换后,放在openfoam里面计算,我的相对振幅总是奔着3到4去了。我试了overset网格,多种不同参数组合(但保证无量纲参数一样),改时间步长,都不行。各位知道是什么原因吗?在这个简单的验证上拖了10天了。。如果帮助,不甚感激。case得压缩包以及论文附在后面。
-
您好,您现在还有做SIXDOF得问题吗?有的话我们可以交流一下,谢谢啦
-
@wwzhao
老师,我之前一直有个小问题没有解决。我在SIXDOF里面定义的时候,要求了输出motion state。angularMomentum (0 0 0); g (0 0 0);//-9.81 rho rhoInf; rhoInf 1;//2375.9;//840;//77.2;//1.924; report on; reportToFile on; solver { type symplectic;//GrankNicolson;Newmark; }
但是只有总有LOG文件里面有,没有输出单独的文件。请问有什么方法,单独把MOTION state输出吗?像我时间步长小,周期大,文件太大了,log文件没有办法处理。我用的V1712
-
@wwzhao
嗯嗯,我想用它先算一个初始场,再放到overpimple里面算。没事,我自己再查查哪里有问题哈 -
@wwzhao
Re100-fr1.0-simple.zip
老师您好,您能看一下我这个oversimplefoam哪里设置的不对吗?
感觉就是一直计算不了,识别不出oversetpatch, 我条件以及参数的设置就是按照tutorial里面那个改的。谢谢!
-
@wwzhao 老师,我问一下,TOPOSET 和 setfields 在OVERSET 定时的作用是什么?我把TOPOSET 和 setfields删掉的时候,好像计算时并不能识别出oversetpatch的重叠网格
-
@netbuger 哥,目前我还没有到后处理的那一步,到了跟你继续交流!
-
@wwzhao 好了老师,问题男孩又来了。我想先用Oversimplefoam算一个基本场,用的tutorial里面那个基础。但是算的时候发现,只在oversetpatch的区域里面进行计算,如图。
我开始以为的U,P这些边界条件没有overset的设置
加进去后发现反而一点都不算了。。如图。。
-
@wwzhao 老师!你太厉害啦!
是这个问题,都是这种小问题注意不到,哎。应该是:
type calculated;
value uniform (0 0 0);
之前是FIXED,我在加拿大这边,这两天晚上找原因找的头发都快掉了哈哈哈 -
各位好,我在做涡激振动,想用overset来处理6dof的网格变形问题。但是计算发现,圆柱的patch在计算时不移动,但周围的其他网格移动,导致第一层网格被挤压,计算崩掉。
我觉得可能是6dof 的coefficient定义里面,只能选择patch的原因?因为像其他solid motion solver,可以直接选择定义的cellzone,进行整体移动,各位有遇到类似问题吗?或者有更好的tutorials,我比较着急,用手机打了一下,代码后面上传。
![0_1574344909239_D51BD371-9335-40E8-998F-ACD1888DF2AF.png](正在上传 100%)
![0_1574344919737_2477F3C6-A10A-478B-A67D-92521F192D0B.png](正在上传 100%) -
@kennyadapter 您好,您用的是2D overset网格吗?
是的话,SNAPPY之后不是要extrudemesh一下,这个有什么要注意的问题吗?
因为我merge之后,两套网格不在一个平面上。。 -
好的老师,谢谢!如果有什么相关的资料方便提前发给我的,如果能发给我,不胜感激。邮箱:z66cheng@uwaterloo.ca, 看您有时间的时候。
-
东岳老师官网那个PimpleDyMFoam解析链接失效了,有谁有的话能发我一份吗?还有就是想知道各种solid motion solver的控制方程是什么,尤其是sixDoFRigidBodyMotion的控制方程。谢谢!!
-
您好!PimpleDyMFoam解析能发我一份吗?东岳老师官网那个链接失效了,谢谢。2316721479@qq.com
-
@VamosZhi 忘说了,是OpenFOAM,用的pimpleDyMFoam
-
有人知道检测这种弹性系统振幅用的工具命令是什么吗?我觉得这个很常用,应该有一个简单的检测命令,只是暂时没有找到。有人知道不?万分感谢。之间只算过固定圆柱的绕流情况。 -
@wwzhao 听说这方面您很懂
-
请问如何实现open foam里面两套网格。两套网格分别采用不同的求解器求解,其中一套网格求解时需要用到另一套网格事先插值进来的参数。谢谢!
-
您好!请问您还有当时无反射边界条件那段的代码吗?我想借看一下,理解后改一改。谢谢
-
上次的问题解决后,面临这个新问题,最近也找了三角函数的定义,但都是关于边界条件的定义,没有明确的类似案例。
总的公式仍然是:只是公式8中的源项换了,换成了:
我已经写好程序去找到这个网格,label cellI; label cellIob; DynamicList<label> p1I; scalar value(1000); vector p1loc(1,0.5,0.5); forAll(mesh.cells(),cellI) { if(mag(U.mesh().C()[cellI]-p1loc)<=value) { cellIob=cellI; value=mag(U.mesh().C()[cellI]-p1loc); } } p1I.append(cellIob); Info<<"location"<<U.mesh().C()[cellIob];
但是如何在定义三角函数源项以及在网格中添加该源项仍然难。如有了解,希望提示一下。谢谢!
-
各位老哥,是我蠢了,我昨天把时间步长降到0.00001,居然算过去了。。。谢谢各位大佬帮助,我现在新问题是,如何在场内某一点上加上随时间三角函数变化的波动项。我自己查查,各位了解也可以提示我用哪种函数
-
@wwzhao 哦哦好的,我再想想看。不过这种方法的流场部分本来是按照不可压缩流算得,这个应该是没有问题的,我的工作是把这部分方程加到计算流场的大循环中。
-
@wwzhao 您好,声场之前的流场计算是不可压缩的。还有,可能是我不太了解,我这段程序发散:
fvScalarMatrix PpEqn ( (rho/c21)*fvm::ddt(Pp) +fvc::div(F) ); Info<<"point4.1"<<endl; solve(PpEqn==-(rho/c21)*fvc::ddt(p));
有哪些可能呢?
-
@李东岳 嗯嗯是的老师,求解公式8
-
经过试算发现,应该是求解Pp的那段程序导致的发散,但是为什么我还不清楚,有什么别的写法使它稳定吗?
-
@李东岳 嗯嗯,那我现在能改吗?公式带‘的项都是波动项,对应于我代码中加了p的变量,如u'对应于Up,p'对应于Pp,也是我要求解的东西,公式中的g相当于一个中间变量,没有什么实际意义,对应于代码中F。
我通过计算出的那几步调用paraview发现,波动项在几步计算后便达到了e23次量级,应该是这个原因出现“浮点数除外”
我把带有F的那几项换成FVM隐式后,求解器就编译不成功了。 -
“但是当把求解F的方程和求解Pp的方程其中一个去掉的时候,计算便能顺利进行。”这句话当我没说。。。。。去掉其中一个后虽然会计算的更久,但还是出现“浮点数例外,核心已转储”。
-
我是年初接触OpenFoam, 之前一直做的是实验,又是硕士联培只有一年,东西从头做,所以时间很紧,还希望了解的人提点一二。做的是风力机气动噪声,现在要在icoFOAM 中把声波扰动方程APE加进去。
这段代码跟在PISO后面写,如下:Info<<"point1"<<endl; volScalarField c21(1.4*(p+Pp+p0)); volScalarField c20(1.4*(p+Pp+p0)); Info<<"point2"<<endl; while (piso.correct()) { Info<<"point3"<<endl; fvVectorMatrix FEqn ( fvm::ddt(F) +fvc::div(F*(U+Up)) +rho*fvc::div(U*Up) -rho*fvc::laplacian(nu,Up) -(rho)*fvc::div(nu*dev2(T(fvc::grad(Up)))) ); Info<<"point3.1"<<endl; solve(FEqn == -(rho)*fvc::grad(Pp)); Info<<"point4"<<endl; fvScalarMatrix PpEqn ( (rho/c21)*fvm::ddt(Pp) +fvc::div(F) ); Info<<"point4.1"<<endl; solve(PpEqn==-(rho/c21)*fvc::ddt(p)); Info<<"point5"<<endl; rhos=rhos+(rho)*(Pp+p-(Pp.oldTime()+p.oldTime()))/(1.5*c21-0.5*c20); Up=(F-rhos*U)/(rho+rhos); volScalarField c22(1.4*(p+Pp+p0)); c20=c21; c21=c22; Info<<"point6"<<endl; }
但是运行一定循环次数后便出现“浮点数例外,核心已转储”
但是当把求解F的方程和求解Pp的方程其中一个去掉的时候,计算便能顺利进行。
不知各位大佬有没有出现过这种情况。
谢谢!
interPhaseChangeDyMFoam 计算水下螺旋桨propeller:叶尖涡空化 不明显
interPhaseChangeDyMFoam 计算水下螺旋桨propeller:叶尖涡空化 不明显
snappyhexmesh结构边缘抓取问题
Re100层流VIV涡激振动振幅不对
OF4.1这行代码释义,动网格问题
overset来处理6dof的网格变形问题
overset来处理6dof的网格变形问题
overset来处理6dof的网格变形问题
overset来处理6dof的网格变形问题
OpenFOAM重叠网格
overset来处理6dof的网格变形问题
overset来处理6dof的网格变形问题
overset来处理6dof的网格变形问题
overset来处理6dof的网格变形问题
OpenFOAM重叠网格
PimpleDyMFoam解析, SIXDOF motion的govern equation
PimpleDyMFoam解析, SIXDOF motion的govern equation
pimpleDyMFoam只有adaptive mesh refine,不存在网格变形,那Ub=0吗?
涡激振动振幅监测
涡激振动振幅监测
Openfoam重叠网格插值
Openfoam重叠网格插值
non-reflecting boundary conditions的定义和初始化问题
在某指定网格加入源项,源项为随时间变化的三角函数
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”
自定义Solver出现“浮点数例外,核心已转储”