@李东岳
老师,高尔夫球和离散气泡模型中的那种效果叫什么名字?在paraview之中要如何实现啊?
100yearsalone
帖子
-
-
@tidedrinker
pimpleDyMFoam是新版本快? -
在非压缩的求解器中,我查了4.1之下的创建通量的
createPhi
文件,发现有两个一个里面写的是fvc::flux(U)
,还有一个是linearInterpolate(hU) & mesh.Sf()
.说明里面写了前一个是相对relative.这两者有什么区别呢?我試了非压缩下的标量运输,fvc::flux(U)
和linearInterpolate(U) & mesh.Sf()
(和上面比少了一个h)两个的结果是一样的,请问下这两者是不是同一种东西两个写法呢? -
@李东岳
ngrd就是两个标量的乘积,第一个式子求了ng,第二个式子求ngrd,最后相除得rd
不过听您这么一说我可能又想错了
这图上面ng0是气泡数,mg0是气泡质量,rg0是气泡半径,rhog0是根据这三个计算的密度。
在运输之前只有一个格子有值,运行运输方程后变成三个格子有值了,可是从密度来看另外两个格子都变成了极大的数,
我就是找个方法比较稳定的传输这几个值,让他们运输以后不偏差密度值太多。 -
@李东岳
老师,您之前给我写了一个气泡数量密度的运输方程是图中的第一个式子
对于有些标量,例如a我想跟随这个N做移动,我设另外一个标量N乘以a等于Na做同一样的运输方程,如图中式子2
求完N和Na以后,再相除求出a,这种方法在OF里面应该是没问题的吧?
我试了一下,好像不太符合我的预期。
我设了ng0作为第一个式子的N,得出下图:
而rg0作为a,a为常值1,ngrd的结果:
我很不解为什么ngrd的运输方程计算后会出现这种结果? -
@李东岳
读取每个cell的值 -
@李东岳
我按照你给我那篇论文里的式子
不含有扩散项和源项
其中mg21就是气泡数量密度,上面的phis2是通量
气泡个数发生的测试条件是这样,我觉得我设置得没有什么问题,但是有一点不太明白
就是在结果里面为什么会出现负值? -
@youmengtian 我当时一直找不到答案,都没管了。。。。。谢谢你
-
@李东岳
你好,李老师
我目前就想试试气泡数量密度的传输方程,我试了几个进行了测试。
你的回答公式无法正式显示,所以我自己查了一下这样定义了。
其中ng2是气泡数量密度,phis2是通量,我感觉应该是没问题的,但是实际运行后结果是
ng2没有发生移动。
我试了一下另外两种定义rg0和rg1,phis0=rhog0U,phis1=rhoU
成功流动,结果一样
rhog和rho
rho和rhog不一样,结果是一样,代表我以后传输方程只用随便定义一个密度是吗?如果不是的话,为什么一开始的式子没发生流动呢?
-
@yhdthu 没处理
-
@mohui
我是在自己写的空化代码中求 Rayleigh–Plesset 方程,得出气相的数据和气体和液体alpha。把原方程求的rho当成混合流的密度来看,再求出液相的密度。 -
大家好,我一直是在基于rhoPimpleFOAM去做一个空化的求解器。
基本思路是;
自己写的空化求解代码 -> rhoPimpleFOAM的代码 -> 自己定义的各个量在网格的守恒方程从第一步的空化求解代码得到四个自定义的值,分别是气泡数,气泡质量,气泡半径和气泡半径速度。
气泡质量 单位kg/m3
气泡数 单位1/m3
气泡半径 实际上是半径乘数量 单位 1/m2
气泡半径速度 实际上是半径速度乘数量 单位 1/m2s
以上四个就是自定义量的守恒方程,用以计算在格子中生成后,跟随流动的变化。
其中自定义通量phit1和密度rhot的定义分别为
在creatfieds处定义了rhot,然后在0文件夹中见了一个rhot文件作如下设置:
在creatfields中定义了phit1,我是仿照通量phi去写的
我自己简单的试了一下好像可以,但是在真的运算中气泡数缺出现了负值如图
不知各位有没有类似定义过守恒方程经验的同学能够提供点意见? -
@李东岳
级测试,如您所说改为NO_READ后不再出现警告,谢谢 -
我自己定义了几个通量,在计算时出现了以下的警告。
我按照提示改成了mustread可是还是出现了相同信息。各位同学有遇过相同的情况吗?
-
@wwzhao
谢谢赵老师的回答,我尝试下写入不同的值看是否能正常写入文档 -
我最近在试着一个OF里面的一个简单的功能实现。
就是在任意一个求解器中加入一段这样的代码:
我想对网格上的一个值alphal进行写入1.0,然后看下会不会在0之类的时间文件夹之下的文档输出。
alphal的定义如下:
这是边界条件的文档:我进行计算后输出的却没有像速度文件之类的输出:
还是像原来一样:
我查了log文件都像这样的写入了,可是为什么没有在时间文件里写入呢?
添加了write的代码还是一样: -
@影川风
我在研究的是rhoPimpleFoam,虽然不一样,但是在rhopimpleFoam也出现了psi,我也想对其进行修改,
我之前也是找不到修改的地方所以我是直接thermo的代码给//掉,自己写了一个状态方程,
但在pEqn里面求解时还是有psi,觉得自己修改得不完全。 -
@李东岳
的确奇怪,我现在都不知道为什么纯粹进行值之间的交换会出现问题,又没有显示命名的重复,
不过我把发生错误的变量名字改了以后,错误就不在出现,虽然解决了,但还是对于错误的原因存疑。 -
我一直都在试,好像错误总出现在cell的值的交换当中。
-
你好,请问你这个问题解决了吗?
-
大家好,
我在调试自己的程序时,好像在读取速度矢量的值时出现了浮点溢出。我是第一次遇到这种浮点溢出,如何解决完全没有任何头绪。
这是我的求解器发生错误的区域,我使用了Info去调试。
在其他网格处是可以运行的,但是好像在某一个网格处在读取u【cell】时就立刻报错了。
我的边界条件也是沿用以前的计算,格子也应该是没问题,请问下各位大神以前有没有遇到过类似的错误呢? -
今天具体查了下书。
在最后的矩阵出有u+c和u-c。
但是我在东岳流体上看可压缩和不可压缩求解器的不同,以及看了下代码。
发现都找不到类似音速c的出现,请问有谁比较清楚音速c在可压缩求解器的应用的嘛? -
一个关于离散格式的问题,如果在可压缩流体的计算当中,流体的音速代表了信息传递的速度,所以在求解的时候根据流体的速度和音速的比来选择离散格式。我查了书有写关于可压缩的计算当中,气体的情况音速c=根号下gammaRT的公式来定义的,请问下这个公式的在哪个文件夹处呢?
-
@李东岳
可是在使用OF本来就有的求解器时不会出现不收敛的情况,而自己修改了求解器就会出现这种情况,会不会是自己不正确的修改令到计算难以收敛呢? -
@李东岳
就是二维只有xy,三次元就是三维旋转机械。
我个人说法有点难以理解。。。。。。
老师,是否清楚在滑移网格中计算出现PIMPLE: not converged within 10 iterations的原因? -
我在二次元的情况下速度的壁面设为fixedValue 0, 二次元没问题但是三次元就像上面一样发散。
printStack的信息代表了除以0发生发散,结果我试了下把壁面设为zeroGradient结果就没问题了。
不知道为什么,但是解决了,以上给大家参考。 -
你是照着openfoamwiki上的来装吗?
-
我使用rhoPimpleDyMFoam去试运行了一下滑移网格,原生的算例annularThermalMixer计算后没有任何问题。
但是我使用自己的模型或者其他算例中的propellar无任何修改的用这个Foam设置也和算例一样但都不行,但是propellar在其他的DyM算例当中运行没有任何问题。
因为计算还没有进到Foam.C里面,所以我也不知道该如何debug,请问下大家有没有类似的处理经验。
-
自己挖的坑自己填回去。。。。。。。。。。。。。
我今天发现是我自己定义的TELF函数出了问题,所以计算出来的结果有问题。
但是这种方法加入状态方程我实践过的确是可行的。
就是换了以后,连续性误差的地方会变的非常怪,查了continuityErr这个文档,发现里面计算连续性误差还是使用rho-thermo.rho所以会出现误差值十分大,
目前还在想着这么去改thermo下的东西,因为自己不能把黏性系数和普朗数设成固定值,
怎样把这些thermo文档下设置函数还在研究中。 -
@李东岳
这几天我一直在试着这个方程,一直都是焓he无法收敛,最后是发散。
所以我改了下方程,让他变成了能量的方程。
其中状态方程是自己定义的,如前面所说
把thermo.rho()换成了这个图里的状态方程,rhoG和rhoL的函数都是自己定义的,这样上述四个式子可以求出内部能量e,EEqn我是这样改了下 原来是thermo.he()的,直接使用理想气体状态方程时是没问题的,但一用这种方式就无法收敛以及发散。
我觉得我是气体状态方程的地方出现了问题,但具体有什么问题我不是很清楚。 -
@李东岳
老师,请问下e是在哪个地方可以获取? -
@李东岳
现在是在做空化 -
@李东岳 谢谢:laughing: ,根据你的指导的确是把Tl成功的输出了。
上面这个是添加在TlEqn下面的代码。
但是输出到Tl的时间步文件中的各个格子的值似乎没有被计算,而是直接被输出了。
我查了下时间步的phi和自己定义的phil输出的也不一样。
下面是phil
而且输出的文档也有些异样
在运算是控制台输出时最后一步TlEqn的计算特别慢,好像卡了一下,从这个文档来看,最后的iterations也到了1000啊,好像是有些问题。 -
大家好,:happy:
我在使用rhoPimpleFoam来计算空化,在其中植入自己的代码,目前第一步是在改了这个压缩求解器中的状态方程.原来的设置是perfectGas就是理想气体状态方程来计算的,我把它改成了一个这样的方程式:
在pEqn的头部
和尾部的地方把rho=thermo.rho()的地方换成了这个.
my_Eqn_of_State的代码,就是把通过气相和液相求出综合的密度,其中RHOBF和RHOL2F是由自定义数据来实现的,根据T和P的不同决定密度.
计算后发现是可以计算的,类那边也看过但觉得对自己改来说太复杂了就用这样的方法去改,不知道方法正确与否希望大家给点意见.
下图是整个主题的代码,加了两行TEqn和ElEqn
因为一开始计算的是综合的密度,所以得出的也是混合体的压力和速度什么的,由混合物的速度可以求出液体的速度.其中混合体的温度T是通过焓方程求出来的.我老师让我把液体的速度再代入焓方程求液体的温度.
所以我在pimple的循环之后加了一个ElEqn的来像计算出液体的温度.
因为现在只是想知道方法,所以直接让混合体速度等于液体速度,这一步根据计算结果来看没有问题.然后我把EEqn文档中的数值符号全部加了l在后面,volScalarField也全部在creatField里面声明了.但计算结果来看,好像有问题,而且也没有输出Tl的值.
请问一下大家我对于这个求解器的修改有什么不完美的地方吗?
还望大家赐教. -
我在使用pimpleDymFoam时出现了以下错误。但是我除了0文件里的边界条件名称意外,完全没有修改过其他的文档。原来文档cellZone处写的是innerCylinder,我打开网格生成后的cellZone文件发现里面有一项wallrotatetip,我就把innerCylinder替换成这个。然后输入pimpleDymFoam,却出现了第一个错误。
1.为什么代码会显示没有solidBodyMotionFvMesh的呢?明明在原example中也能运转啊。
我按照代码所示,把solidBodyMotionFvMesh换成了dynamicMotionSolverFvMesh。后出现了第二个错误。
2.提示所指的keyword solver到底指的是什么啊?如何修改才能使这个Foam正常运行呢?
谢谢大家
-
问题已解决,但是还是不懂为什么,用fluentMeshToFoam会出现warning而且时间特别久,fluent3DMeshToFoam不用很久得可以了而且没warning
-
跟大家分享下解决方法,把controldict里的write mode改成ascii再进行格式转换就可以了!!
-
请问下如何转换格式,谢谢啦
-
有人在转mesh后出现这种乱吗? -
@cfd-china 你好,今天下午我试了一下按网上openFoam滑移网格教程去做,内容是
‘本次教程主要说明一下openFoam中的滑移网格,本次的算例的网格采用的是simCafe.org网站上的网格文件,将网格文件拷到新建的AMI文件夹中,创建system文件夹和system文件夹下的controlDict文件,然后执行fluentMeshToFoam命令将fluent网格转换为openFoam格式网格,命令输入如下图所示,注意fluentMeshToFoam命令后要添加一个选项-writeZones,则转换完后openFoam网格单元组cellZone文件中会有各个旋转域的信息,省去后续的网格操作,此外该命令只能转换ascii码的fluent网格格式,需将2进制网格文件转换为ascii码格式。’
我按照这个步骤,下载了msh文件导入后,提示没有controlDict,我吧pimpleDymFoam的propeller的controlDict文件复制过来,但再运行了发现这个文件与原来文件夹的surface和force文件都有关系,我打开看内容好像和原来的propeller都有关系。请问如果自己到入ICEM模型时,controlDict要如何处理呢?
-
不同的zones指的是什么呢?我在这个模型是可以在FLUENT进行旋转运行的,画的时候已经分开成不同的zones了吧,管理员有没有这方面的资料可以参考一下呢
-
大家好,我是用OpenFOAM做流体旋转机械的研究的。想把自己做的ICEM模型导入Dym的求解器里面看下能否运行。下图是我格式转换时的出现的代码, 小弟第一次转换不知道是不是哪里出了错。
然后我就打开paraFoam
我运行原来的算例builtin是有好几个的。但是现在导入ICEM模型后就只有一个了,中间旋转的Zone怎么去设置呢?如何实现转动呢?忘不吝赐教。
酷炫的流体动图,堪比电影艺术
OpenFOAM5.0计算速度比4.1更慢?
相对通量和通量的区别?
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
可压缩求解器中的音速的定义和离散格式
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
自定义守恒方程时出问题
自定义守恒方程时出问题
关于计算后的数值写入的问题
关于计算后的数值写入的问题
怎么才能在代码里快速找到函数的定义?
读取速度矢量时出现浮点溢出
读取速度矢量时出现浮点溢出
Foam::sigSegv::sigHandler(int)错误
读取速度矢量时出现浮点溢出
可压缩求解器中的音速的定义和离散格式
可压缩求解器中的音速的定义和离散格式
使用rhoPimpleDyMFoam计算在开头就无法进行下去错误
使用rhoPimpleDyMFoam计算在开头就无法进行下去错误
使用rhoPimpleDyMFoam计算在开头就无法进行下去错误
如何在centOS/Redhat6.2上编译OpenFOAM4.1
使用rhoPimpleDyMFoam计算在开头就无法进行下去错误
关于在rhoPimpleFoam里面修改EEqn和状态方程的问题
关于在rhoPimpleFoam里面修改EEqn和状态方程的问题
关于在rhoPimpleFoam里面修改EEqn和状态方程的问题
关于在rhoPimpleFoam里面修改EEqn和状态方程的问题
关于在rhoPimpleFoam里面修改EEqn和状态方程的问题
关于在rhoPimpleFoam里面修改EEqn和状态方程的问题
关于dynamicMesh文件的几个问题
ICEM的旋转模型导入pimpleDymFOAM算例中运行
mesh转化之后出现乱码
mesh转化之后出现乱码
mesh转化之后出现乱码
ICEM的旋转模型导入pimpleDymFOAM算例中运行
ICEM的旋转模型导入pimpleDymFOAM算例中运行
ICEM的旋转模型导入pimpleDymFOAM算例中运行