@东岳 试了一下,等待结果中~
LiuZhuang
帖子
-
-
也烦请各路大神看到这个帖子,成功加载过此类边界条件的给点指导建议,小弟在此谢过了~~
-
原谅我又一次没发代码,因为不知道这种代码怎么写,在OF里是否有这种边界条件,该如何设置?
有一个圆柱面,我想在圆柱表面添加一个指向圆心的固定速度,大小为72m/s,请问这种条件在OF里有没有这种边界条件,我该怎么去设置呢?圆心的坐标为(-615,2892.5 0)圆柱的轴方向为(0 1 0)
PS:我的坐标系是直角坐标系
@东岳 -
@麦迪文 等着看云图
-
@东岳 @麦迪文
感谢东岳先生和大佬,问题已经解决了,phi1=phi+phi也实现了,浓度场也变化了;
现在想问一下关于声明里有一段代码理解的依然不是很好这个是我仿照U的声明做的一个Uslip,目的是为了读取0文件夹下的Uslip。
这个是我新生命的phi1场,其目的是为了让我的Uslip从体心插值到面上,这个思路是否正确呢?
这句代码的含义是将速度从体心插值到面上并且初始化(没理解错的话)
ps:我的网格文件实在是有点大,而且处理成5000以下计算效果很不好,就不麻烦东岳先生帮我验证了。感谢东岳先生百忙之中帮我研究我的问题,这些问题很基础,但对于小白来说也很棘手,感激东岳先生,拜谢!
-
@100yearsalone 新写的phit1里的U,这个U是从0文件夹下的U调用的么?
-
@东岳 确实在方程中添加了phi1=phi+phi,也影响了浓度场的分布,现在我遇到的问题是phi 和 phi1文件里的值还是一摸一样的,正常phi1的值也应该和phi 有所区别,我如何才能让phi和phi1的值有区别呢?问题是否出在声明文件中呢? 这是我写的phi1声明。
这是我写的浓度方程文件,在括号里加入了phi1 = phi +phi; -
@东岳 或者说我想在浓度方程的对流项中将phi改为phi1,phi1的表达式为phi+phi(第二个phi想改动成自己写的速度场文件(即Unew)&mesh.Sf)
-
@东岳 这个新的phi1我是想放在浓度方程 中去的,而且第二个phi想达到一个Unew乘以网格插值出来一个新的面通量,如果压力方程里做phi = 2*phi1的计算的话,会不会影响连续性方程和动量方程呢?我并不想去影响这两个方程
-
@东岳 果然是我的phi1写的方式不准确,我在createFields.H中声明了一个phi1,如下图
我的想法是创立一个新的phi1让phi1=2倍的phi,不知道这种写法错误在哪,做出来的phi和phi1两个文件值相同。
如何创立一个新的phi1呢?除了在createFields里填写还需要在哪注意呢? 我那个phi+phi的写法应该错误很大吧,应该怎么去实现phi1=phi+phi的感觉呢? -
@金石为开 请问您的driftFluxFoam关于气液两相流的问题有没有算例的设置呢?能不能给我分享一下?
-
@东岳 层流施密特数1,湍流施密特数800,而且入口速度非常小(0.1m/s)这样的湍流扩散项系数应该不会大吧,我也在怀疑对流项占优和扩散项占优的问题,如果这种取值很小的话应该就是我写的文件没有让openfoam明白我想给Unew添加进去的想法吧,苦笑~
-
@东岳 我做的测吹方式要比刘仲秋老师(东岳先生您发给我的文章)的模型简单,我是用浓度方程去代替含气率方程(他们在方程形式上有很高的类似度)而且以前有前人用混合单项流模型去模拟的E-E模型,现在我已经做好了滑移速度在模型中的分布文件Unew,和导师探讨后认为Unew文件应写入到我之前写的浓度方程的对流项中去。然后我发现pEqn.H中关于phi 进行了一系列眼花缭乱的操作,我进行了简单的复制后发现wmake报错(我的想法是按照pEqn.H简单的复写一个phi1,给phi1变成U+Unew的形式),这种操作对于我来说十分吃力,我弄不清楚如果添加这个新的phi1应该对除了createFields和浓度方程之外还有哪些文件进行改动(比如说或者pEqn.H)说实话CFD理论理解十分浅显,根本看不懂Weller的操作,希望东岳大神帮我指点一下迷津
ps:我将混合单项流模型的论文里含气率守恒公式和我做的简单的pEqn.H粘贴在下方。
这个是湍流nu
pEqn.tar 这个是我做简单改动后的pEqn.H
标蓝处我做的简单添加,在文件的第48行
我用的OF版本是2.4.0的,基于buoyantBousssinesqSimpleFoam求解器进行的改动 -
@东岳 对了东岳老师,我听说最近你也在做测吹的研究,请问您采用的例子有什么可供参考的么?我也在做这方面研究,现在就是速度太大导致计算很慢。希望能与您多加交流
-
我进行了测试,发现原来的phi和新的phi1做出来的浓度分布毫无变化。。。。。。。
-
@东岳 这个是我在createFields.H里做的phi1
这个是我做的新的速度矢量
这个是我新写的浓度方程
-
@东岳 我也在做测吹方面的研究,但是E-E有点算不起(因为我的吹入速度很大)东岳先生的结晶器图片我看过,很漂亮。
-
@东岳 Unew目前是一个矢量vector(1000 0 0)和U比起来大了很多,我用了我自己画的网格进行测试(浓度入口很小)算例使用的是buoyantBOussinesqSimpleFoam的例子,仿照的边界条件进行了改动。
-
或者说我的pEqn.H文件写的不全,导致了phi1虽然满足编译,但是并没有离散且参与到浓度方程的计算中去?现在不了解问题的错误在哪里?非常的头疼!!!
-
东岳老师及各位大神们早上好!
本人的课题遇到了一个新的瓶颈,在此发表希望大伙给与指导与点拨。
我想在浓度方程对流项里引入一个速度,并且通过简单的学习将原速度U和新速度Unew引入到phi中,将phi定义为phi1(编译成功),但是我将phi1放入浓度方程的对流项(动量方程与压力泊松方程的对流项仍保持为phi)中时,发现浓度的扩散并没有Unew的效果。
在这里我想问的是,我这样做是不是也需要考虑将动量方程和压力泊松方程的phi也改成phi1?并且对浓度方程里新的phi1应该怎样进行调整(是否需要在动量方程和压力泊松方程里也同样调整phi1)?才能使我加入的phi1能够对浓度的扩散有一个流动的效果?
下面我将改动的求解器代码发到这里,希望大神们能帮我搞清楚这个问题ps:一直对压力方程理解的不是很好。希望大神们给与指导,这些文件我是基于buoyantBoussinesqSimpleFoam的基础上加入浓度方程进行修改的。
-
/OpenFOAM/OpenFOAM-2.4.0/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/0/T.water
这是我想修改的温度文件的位置。
初始蒜粒里面wall给的边界条件是zeroGradient.我想做一个壁面散热的效果(壁面以热通量W/m2的形式向外散热)很明显zeroGradient不再适合。
我想求教:这种热通量向外散热的边界条件该如何设置?是采用fixedGradient还是其他形式的壁面函数?
如果是fixedGradient,那么它的具体值的单位是什么?如果采用其他方式,那么有没有能使twoPhaseEulerFoam运行起来的成功的壁面函数?如果有,还望各路大神不吝赐教!!!
最后还是麻烦东岳老师帮我看一下,或者号召一下大家伙,看看我的问题。拜谢拜谢!
@东岳 -
目前研究了原始的bubbleColumn,发现出口设置属于inletOutlet混合模式。仔细研究了一下初始算例,发现他们有一段空气域才如此设置。
现在我的问题是顶部没有那段空气域,可能所有的出口条件就不适用于inletOutlet混合形式。
现在对出口条件十分疑惑,我的想法是气泡在outlet部分自由逸出,液体在outlet部分作为自由液面(固定壁面)流动。不知按照这种思路设置的outlet采用zeroGradiet【fixed Value(0 0 0)】还是有其他的outlet边界条件设置。希望能与fomer们交流,求教
@东岳 -
东岳先生您好,我现在在用twoPhaseEulerFoam解决气孔侧吹的问题;
请问壁面润滑力中Hosokawa模型在OF里怎么实现调用?
还有一个问题是气体入口速度很大(200m/s)的情况下,气体的k和epsilon初始值就不会与液体的k和epsilon初始值相同而是差很大的数量级吧?(液体初始在计算域内静止不动)
另外,有什么比较好的计算程序能这气相和液相的k、epsilon定的准确一点?
多次叨扰,感愧于心,拜谢! -
@东岳 或者我是否可以这样理解,当value或者inletValue的值设为0的时候,他们对应的边界条件为zeroGradient,如果他们的值不为0的时候,对应的边界条件为value设置的fixedValue的值?上个楼层有问题,我在贴一段代码
这个type取zeroGradient还是fixedValue是根据phi是否为0确定,0则为zeroGradient,非0则为fixedValue,而fixedValue的值则取决于value值的大小,这种理解是否有误??
多有麻烦,不胜叨扰,在此拜谢! -
@东岳 或者我是否可以这样理解,当value或者inletValue的值设为0的时候,他们对应的边界条件为zeroGradient,如果他们的值不为0的时候,对应的边界条件为value设置的fixedValue的值?
-
@东岳 我研究了一下他的初始算例,对其中的pressureinletOutletVelocity边界条件不是很清晰。我想做的工况是在outlet边界上气体自由逸出,液体则为壁面(zeroGradient),这是我的OF自带的原始边界条件,请问我该如何设置?
PS:inletOut边界条件是否也需要做改变?如何让所有的边界条件满足这种工况? -
@气固两相流 请问您看过的一些文献资料都是哪些?求分享,我学习学习
-
东岳老师您的邮箱多少?关于twoPhaseEulerFoam的问题我有几点想跟您请教
-
我的OF2.4.0也遇到了这个问题导致计算发散,而且在计算中我不需要考虑相间传热,只考虑液体的壁面散热,我该如何解决这个问题呢?@东岳 @CFD中文网
-
@东岳 最近根据东岳先生的指导,对这个求解器不那么陌生了,也简单做了一个自己的小例子。目前程序能跑起来。拜谢东岳先生!!!
这是我做出来的一点东西:这个是两相流的含气率分布。但是我看了一些论文在顶部设置的边界条件叫做液相壁面采用无滑移边界条件。气泡达到顶部以上浮速度自由逸出。不太理解这句话对应的OF边界条件是什么?
PS:目前顶部采用的边界条件为:U.airU.water
p文件:
目前还是依据算例进行自己改动,仍然没达到气泡以上浮速度自由逸出的效果,请问东岳先生对此有什么好的建议么?
叨扰叨扰,再次拜谢!!! -
@东岳 https://github.com/OpenFOAM/OpenFOAM-2.4.x/blob/master/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/0/p
这个文件的代表了什么意义呢?为什么内部的值高达10的5次方 -
@东岳 https://github.com/OpenFOAM/OpenFOAM-2.4.x/blob/master/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/0/k.air 这个是多相流的k.air文件;
https://github.com/OpenFOAM/OpenFOAM-2.4.x/blob/master/tutorials/multiphase/twoPhaseEulerFoam/RAS/bubbleColumn/0/nut.air 这个是欧拉两相流的nut.air文件
我想问的是nut.air文件里面的object为什么是k,他们两个代表什么意义?nut.air是否是湍流雷诺数的意思? -
@东岳 也就是Theta文件,上一条回复说的有点含糊
-
@东岳 收到收到,我去拜读一下。如果不需要设置的话我能否删除这个文件。因为并行计算的时候总是提示我这个文件有错误。但是单核计算的时候却能够运行
-
@东岳 确实是我考虑欠妥。现在遇到的问题是在twoPhaseEluerFoam中遇到了一个Theta文件,我查询了量纲是m2/s2,这个文件的物理意义代表了什么?在做气液两相流的时候应该怎么去设置这个文件?
-
@东岳 或者说我做气液流动的问题能否直接删除Theta文件?
-
@东岳 气液不需要设置Theta的话,那么还有哪些也不需要改动?我现在不动Theta能算,但是并行计算的话无法计算
-
@金石为开 您好,您毕业了么?还做这一块关于滑移和气液两相流的问题么?如果有的话请加我QQ:694467562,我有一些问题想咨询您
-
@东岳 这个模型我很需要,请问东岳先生是用OF哪个版本的哪个求解器实现的?能否引入侧吹(也就是气泡从侧方吹入,而不是底部吹入)?
-
@东岳 感谢东岳先生的鼓励,目前我的理解是将这个model放入/home/dcy/OpenFOAM/OpenFOAM-2.4.0/tutorials/heatTransfer/buoyantBoussinesqPimpleFoam/hotRoom/constant/RASProperties中的RASmodel,将原始的KEpsilon更改为mixtureKEpsilon,我不知道这个思路是否正确,还需要做哪些改进??
-
@东岳 看了东岳老师对twoPhaseEluerFoam的解读,现在我想用mixture做气体侧吹进入液体中的例子。但是哪些算例里能够使用mixtureKEpsilon模型?mixtureKEpsilon是否考虑了曳力虚拟质量力等源项力呢?
-
查阅了一些资料,了解了phi是一个通量,含义为Uf*Sf,那么我在温度方程里的对流项中的Ux变成Ux+C,请问我该怎么做?
-
然后这个phi的意义是不是速度矢量u?这个phi是不是一个对象,它的位置在哪里?我想在这个phi的z方向上加一个固定的速度,请问我该怎么做?@东岳 -
@东岳 那个我添加成了,可能是我考虑多了,加和不加区别不大
-
@东岳 alphatJayatillekeWallFunction;这是0文件alphat文件夹下的边界条件,我仿照这个设置了A边界条件,其意义是将nut/Prt+nu/Pr的热扩散系数更换为nut/Sct+nu/Sc的浓度扩散系数给出浓度扩散的边界条件。不知我的思路是否正确
-
这个是按照温度方程改变的浓度方程
-
@东岳 原先的温度变化比较合理,自从我加入了自己编写的浓度方程和新的nut/Sct+nu/Sc的边界条件后,温度场就呈现出了浓度扩散那种一圈一圈的变化
-
@东岳 ,东岳老师,我做这个求解器的浓度方程添加进去了,程序也能算了。浓度的效果还不错。但是我现在遇到的问题是温度的变化也呈现出一种类似于浓度扩散的变化,请问是为什么?我需要把哪些代码发上去给您看看???
-
我目前能回答你第一个问题,这个是考虑热浮力情况下的计算,TRrf是参考温度,一般取液相相变温度。Pr一般可以取一个很大的数,Prt一般取1.p-rgh一般不用修改,如果想了解的话可以看看方程里的代码,他的目的是为了使计算更准确(p—rgh先加后减)。alphat的意思是nut/Prt+nu/Pr。也就是动量方程里的有效层流黏度和湍流黏度。最后一个是一种p+p_rgh混合的感觉,不用了解p的边界取法。
我也是新手,希望我这点浅见您能够去其非取其是,也希望对您有点小小的帮助 -
@regina 动力黏度在constant/transportpropretities那个文件夹下
想在边界添加一个指向圆心的速度
想在边界添加一个指向圆心的速度
想在边界添加一个指向圆心的速度
求助:液池热毛细对流算例,已编译好求解器和算例,但就是和文献自编程序不一样
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
自己写的一个空化求解器的气泡数,气泡质量,气泡半径和气泡半径速度的守恒方程的一些疑问
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
使用QBMM模拟气泡流
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
升力、曳力、虚拟质量力等
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
关于在浓度方程对流项中添加U+sudu(sudu添加无效果)的问题
of2.4.0下twoPhaseEulerFoam气液两相温度设定问题
twoPhaseEulerFoam下的气液两相流边界条件设置问题
升力、曳力、虚拟质量力等
twoPhaseEulerFoam下的bubbleColumn算例的问题
twoPhaseEulerFoam下的bubbleColumn算例的问题
twoPhaseEulerFoam下的bubbleColumn算例的问题
twoPhaseEulerFOAM
TwoPhaseEulerFoam求解器算例设置
twoPhaseEulerFoam 热传递
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
请教气液两相流动模拟,我采用了vof,mixture和欧拉欧拉多相流模型,但结果很不相同。请指教。
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
有用twoPhaseEulerFoam求解器计算管道内固液流动的算例吗?
多相流与矩方法的数学模型
OpenFOAM模拟管道流
请教气液两相流动模拟,我采用了vof,mixture和欧拉欧拉多相流模型,但结果很不相同。请指教。
请教气液两相流动模拟,我采用了vof,mixture和欧拉欧拉多相流模型,但结果很不相同。请指教。
buoyantBoussinesqSimpleFoam中TEqn.H的问题
buoyantBoussinesqSimpleFoam中TEqn.H的问题
关于buoyangBoussinesqSimpleFoam添加新的边界条件
关于buoyangBoussinesqSimpleFoam添加新的边界条件
关于buoyangBoussinesqSimpleFoam添加新的边界条件
关于buoyangBoussinesqSimpleFoam添加新的边界条件
关于buoyangBoussinesqSimpleFoam添加新的边界条件
buoyantBoussinesqSimpleFoam的使用
关于simpleFoam求解器中设置流体属性的问题