CFD中文网

@hustzzw 你这是*号都被格式化吃掉了吧。。看起来奇奇怪怪的

可以加些输出语句，看看y_h2这些量算的对不对。（小心除零）

好的，谢谢

我在使用solidDisplacementFoam求解器进行热应力计算时，需要对位移向量场重新设置边界条件：其中的x,y方向为zeroGradient边界，z方向为fixedValue边界。请问OpenFOAM里有类似的边界条件或者知道该怎么设置的吗？？

咦，这个工作我应该是2015年做的，怎么2017年才发帖..

用STARCCM做船舶阻力分析，考虑升沉和纵倾，之前流域底部设置为速度入口，结果会收敛。考虑浅水效应即流域底部的影响时，我把底部的边界条件修改为壁面 结果不收敛 且提示如图警告-the mesh may have non-positive volume cells，请问前辈是壁面条件不能修改为壁面而是其他，还是修改为壁面后需要再修改其他的地方，谢谢

@李东岳 我发了帖子之后才发现你之前发过通知了，不是催更，抱歉抱歉~:chitang:

背景是模拟固体火箭的射流问题。
考虑到固体火箭产生的射流中存在$Al_2O_3$颗粒，需要采用多相流的方式去模拟流动。目前我采用Fluent中的欧拉-拉格朗日方法初步开展了流动仿真，对比有无颗粒，流场还是有变化的。但是接下来采用DO方法引入颗粒辐射后，计算得到壁面处的辐射热流值很小，似乎都不影响总热流结果。
有人有这方面的经验不？

@东岳 原来是连续方程的引用啊，谢谢李老师啦

@李东岳 谢谢李老师

央企中国电子信息集团所属企业长沙军民先进技术研究有限公司招聘
数量：6人
工作地点：长沙市
岗位名称：流体计算软件研发工程师
岗位要求：
1、流体力学相关专业；
2、熟悉GPU和CUDA编程；
3、熟悉C++编程；
以上三项条件满足二项即可，满足三项者优先录取，流体力学强相关专业者也可面谈。

基本月薪待遇8K-15K，其他待遇面谈，特别优秀的硕士博士待遇还可以商谈。
联系人：郭先生，18908480199
bf037a90-2b9b-464c-83da-3b501a4e6431-image.png

@冰冰玩 好的，谢谢！我之后尝试一下。

@yan 我想问的是壁面边界上的粘性应力，您说的这个应该不是吧？

按照twoPhaseEulerFoam中的EEqns.H写出能量方程如下（算例文件中能量形式是sensibleInternalEnergy）大神们先看看这么写对吗？

fvScalarMatrix E1Eqn ( fvm::ddt(alpha1, rho1, he1) + fvm::div(alphaRhoPhi1, he1) - fvm::Sp(contErr1, he1) + fvc::ddt(alpha1, rho1, K1) + fvc::div(alphaRhoPhi1, K1) - contErr1*K1 + ( he1.name() == thermo1.phasePropertyName("e") ? fvc::div(fvc::absolute(alphaPhi1, alpha1, U1), p) + p*fvc::ddt(alpha1) : -alpha1*dpdt ) - fvm::laplacian ( fvc::interpolate(alpha1) *fvc::interpolate(thermo1.alphaEff(phase1.turbulence().mut())), he1 ) ); E1Eqn.relax(); E1Eqn -= ( Kh*(thermo2.T() - thermo1.T()) + Kh*he1/Cpv1 - fvm::Sp(Kh/Cpv1, he1) + alpha1*rho1*(U1&g) + fvOptions(alpha1, rho1, he1) );

a8753d25-8ea0-4cf2-8fae-10a71aa5cdfd-image.png
那如果不考虑传热的话，只把导热系数设为0和温度设为一样，这样做对吗？
我为了不考虑传热，是直接把twoPhaseEulerFoam.C中的EEqns.H给注释掉了，我这种做法对吗？

@netbuger :huahua: :huahua:

新人发帖，还不懂如何编辑帖子，请各位老师见谅。
捕获.JPG
该公式是安德森的空气动力学基础一书中对等熵条件下准一维流动推导出的速度面积关系式。根据该关系式，当管道流通面积A不变时，速度u是不变的，则根据质量守恒关系可以知道密度和压力都是不变的。但是我在观察真实条件下的抽气管道流动中，在管道流通面积不变的情况下，上游管道静压大于下游管道静压（下游安装有抽气机组）。请问这是否说明等熵条件下准一维流动速度面积关系式不适用于真实条件下的均匀管道流动？
流体力学小白，还希望有老师能指点一下。

$\frac{\mathrm{d} A}{A}=\left ( M^{2}-1 \right )\frac{\mathrm{d} u}{u}$

分别将小球和圆柱建立两个body，再采用结构化网格划分球体，再设置交界面节点对应生产圆柱体网格

2021-09-17 15-46-22屏幕截图.png

WechatIMG47.png

请问一下, 如果得到这个bubble的projection area和horizontal diameter呢?
我现在尝试两种方法, 第一种对气泡旁的alpha数据在python中进行3D griddata处理, 得到alpha < 0.5的数值, 然后将所有坐标投影到一个平面(比如 z=0), 然后使用CV2进行图像分析,得到半径, 圆心和投影面积
但是因为需要计算量巨大, 我们在griddata中的space number取30-50, 这样造成误差比较大.

第二种办法是使用openfoam的后处理sample获取 z=z_bubble的constant plane然后进行2D griddata, 得到的数据比较好, 但是问题是当遇到图2中的气泡时, 这个方法就不好用了, 请问openfoam有没有一种后处理方法可以将气泡表面捕获然后给出projection的坐标呢?

ICEM用block画的多个区域，怎么各自设置成单独的part?