CFD中文网

@chengan-wang 你好，请问你说少了一个定义是啥？我好像遇到了和你一样的问题，先谢谢啦。

sed -i "s/oldString/newString/g" `grep oldString -rl /path`

谢谢，已更新

这个是当我把水的物性从常数改为多项式后遇到的问题
Courant Number mean: 0.134868655 max: 0.237479859
deltaT = 3.4820232e-06
Time = 0.00378807

PIMPLE: iteration 1
MULES: Solving for alpha.gas
MULES: Solving for alpha.liquid
gas fraction, min, max = 0.00155200072 1e-15 0.0672465558
liquid fraction, min, max = 0.998447999 0.932753444 1
Phase-sum volume fraction, min, max = 1 0.999994574 1.00026614
MULES: Solving for alpha.gas
MULES: Solving for alpha.liquid
gas fraction, min, max = 0.00155951826 1e-15 0.0674835707
liquid fraction, min, max = 0.998440482 0.932516429 1
Phase-sum volume fraction, min, max = 1 0.999994794 1.00025528
MULES: Solving for alpha.gas
MULES: Solving for alpha.liquid
gas fraction, min, max = 0.00156703579 1e-15 0.0677203557
liquid fraction, min, max = 0.998432964 0.932279644 1
Phase-sum volume fraction, min, max = 1 0.999995014 1.00025493
Constructing momentum equations
Tf.gasAndLiquid: min = 350.131001, mean = 456.035361, max = 549.778465
iDmdt.gasAndLiquid: min = -0.959612341, mean = 11.8287784, max = 447.419943, integral = 0.000182163188
wDmdt.gasAndLiquid: min = 0, mean = 138.094307, max = 5727.35791, integral = 0.00212665234
smoothSolver: Solving for e.gas, Initial residual = 0.104953394, Final residual = 1.40390231e-09, No Iterations 2
smoothSolver: Solving for e.liquid, Initial residual = 0.282448508, Final residual = 6.30884149e-18, No Iterations 2
[stack trace]

#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4 void Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#6 Foam::heatTransferModels::RanzMarshall::K(double) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#7 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::BlendedInterfacialModelFoam::heatTransferModel::evaluate<double, Foam::fvPatchField, Foam::volMesh, double>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > (Foam::heatTransferModel::*)(double) const, Foam::word const&, Foam::dimensionSet const&, bool, double) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#8 Foam::BlendedInterfacialModelFoam::heatTransferModel::K(double) const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#9 Foam::ThermalPhaseChangePhaseSystem<Foam::PhaseTransferPhaseSystem<Foam::TwoResistanceHeatTransferPhaseSystem<Foam::MomentumTransferPhaseSystemFoam::multiphaseSystem > > >::correctInterfaceThermo() in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactingMultiphaseSystem.so
#10 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingMultiphaseEulerFoam
#11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#12 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingMultiphaseEulerFoam

shurufa diu shi
Could you please check the other terms? e.g., lift force

针对一个高度为3mm，宽度为1mm的矩形微通道，用fluent mesh 画网格后，在入口流速为0.3 0.4 m/s时，计算都可以收敛。但是当入口流速为0.5m/s时，计算收敛不了，这可能是什么原因导致的，尝试过增加网格数量也无效

看了一下你的模型，之所以你在相贯线处网格画不好，这主要问题在于你的几何模型——相贯线与底面相切，出现了三角区域。从而导致网格在该三角区域网格质量奇差无比。如下图所示

该问题导致出现了如下问题：

解决的办法是将底面拉伸一段距离，避免出现相贯线与底面相切即可，如下图所示
0_1536142189833_4ccfa31a-f8b8-4435-87b6-c322423cd7e1-image.png
如此，便可愉快的进行分块和关联
0_1536142304777_9fb3a5ee-3695-47f4-971f-61fe9c119d23-image.png
接着对其进行O型块划分
0_1536142342178_11fb8b28-9d22-441a-978f-e0f2d4ed1c7c-image.png
最终网格效果图如下
0_1536142436479_986e33bd-8587-4ee0-a010-5f914675e718-image.png

画周期网格不行吗，取出1/10画好就行了啊

@wwxdli 谢谢，我试一试

可以，例子google一下都有。

语法参考http://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf

@Wayne 对，电磁场里面会用到

您好，请问这里的Re_thetac如果不修改，即为none时，Re_thetac的默认值是多少？

倒是没有在研究这个，最近在用有限体积法求解翼型绕流流函数的问题，用到一点点。谢谢你，我找来看一下:happy:

各位大佬，在编译udf出现问题，为什么会出现解码失败呢，flunet不是使用UTF-8 解码，为什么会出现gbk' code，从而导致Error at host: open_udf_library: couldn't get udf_data
???????????????

UnicodeDecodeError: 'gbk' codec can't decode byte 0x8f in position 155: illegal multibyte sequence:

注：路径检查没有中文和非字符

@小葵 :threaten:

@李东岳 老师你好，我想请问下，设置了重力加速度g后，dynamicDict下不设置初速度，物体不能下落，这是哪里没设置好吗？

本人小白，有个疑问：为什么在做流体仿真的时候设置完边界条件后，还要进行初始化后才能开始计算。而做固体有限元分析的时候设置完边界条件后就可以直接开始计算了？照理说微分方程有边界条件后就可以确定唯一的解了，想不明白初始化的目的是什么，各位老师可不可以深入浅出的给我讲解一下？谢谢

@liqiwen
Hello, 欢迎加入！

流体理论中层流中的边界层的确定使用的是普朗特的理论，而ansys flunt中加入边界层只需要把流体域边界处的网格加密，无需其他设定就可以计算了，请问普朗特的理论在cfd中是怎么引入的；

普朗特什么理论？虽然边界层分为层流边界层和湍流边界层。但是在计算中主要调查湍流边界层。CFD软件把网格加密或者加密到合适的高度，是因为边界层内速度和网格距离存在一定的关系。这个关系成为壁面法则。要实施壁面法则，壁面处网格需要满足距离上的一些要求，因此需要加密网格。湍流神秘y+

湍流和层流的边界层理论在flunt中的设定是否一致？

湍流可以使用壁面函数。层流边界层只是速度变化比较大。不能使用壁面函数。因此需要细化网格。

计算火焰对铜管的加热，火焰仿真已经解决，现在需要考虑流固耦合对流换热和铜的热传导。
方案1.网格目前只考虑流体网格，暂时不考虑固体网格。已知fluent中有类似的薄壳边界可以实现，openfoam中pyrolysisZones貌似能实现。
方案1存在问题：1）设置指导；2）pyrolysisZones能考虑三维传热吗？还是只能一维径向？
方案2.划分流体区域和固体区域，耦合计算。
方案2存在问题：固体部分不熟悉，网格耦合和区域划分设置求指导还有就是tut参照。:140: