CFD中文网

@Alex 太厉害了，大神，我想学习icem，您还有别的类似的icem实例教程吗，不胜感激

列表

还有监测的一个点的密度变化
微信图片_20191214213757.png

整体划分，然后删除多余的block，看看好用不

这个啊，先排错，
看看是“只有20可以”，还是“数字小了才可以”。
比如试试19，21，试试200。

错误定的越准，描述才越准，也越容易搜到前人的经验

foamer新手，求指教如何往风场中加颗粒？

对的。第二张图中的总流量就是单个喷嘴的流量（20kg/s 除以 10）

论坛类似的问题好多。不太好回复。主要是这样：类似这种问题，只要能用函数来表示，在OpenFOAM里面植入通通不是问题。重要的是：

如何用函数表示 用函数表示之后，如何写进入

我觉得你们应该知道如何用函数表示，但如何在OpenFOAM中编程进去，这个几句话说不明白...

@Leon-Zhao 感谢大兄弟

OpenFOAM中如何给初始条件添加合适扰动呢？除了延长入口外（太耗费计算资源），有什么别的方法实现充分发展湍流初始场吗？

@ibelief 朋友您好，能否再分享一下？是同位网格版本的吗？

@东岳 Yu老师的工作重心最近两年主要放在把TEMOM应用到更多的场景，就是推广它。深挖的工作，我们这边团队小，没博士，老师越来越忙……老师是有意向的，只是条件不允许。对于其他矩方法的研究目的更多是为了对比验证，增强说服力，然后可以更好的推广。当然这些都是我猜的，毕竟老师不跟我汇报工作:chigua2:

是做巷道通风仿真么，带条缝的附壁风筒？

缩小模型只能保证一个无量纲数相同，一般的流动考虑雷诺数就行，缩放会带来一些非真实的物理效应，但是计算量不够就不要讲究了

@FluidGao 哦，我下来看看，谢谢

@d-benjamin 这个是正则表达式，是用来匹配字符串的。"."代表任意字符，"*"代表出现多次（包括0次），合起来".*"代表匹配任意字符。

"epsilon.*"可以匹配以epsilon开头的任意字符（包括本身），如epsilon，epsilonFinal等。其中epsilonFinal是用于PIMPLE算法最后一个迭代循环的控制参数。

以上规则同样适用于"U.*"和"k.*"。

是否为容差 matchTolerance 太大导致的衔接不上~

@小八 ICER，文章已经找不到了。化工方向的:+1:

简单的点线面操作撤销是没问题的，检查拓扑或者是操作包含的数量太多ICEM就没办法了，解决这个问题的最好办法就是你在做包含多步操作的问题时候先另存为一下是最保险的，问了一下学长他也是用这个方法的

@b占余文乐 当然以上这些都是我个人见解，根据使用感受来的，有没什么官方的消息

各位老师同学好。本人最近在用InterPhaseChangeDyMFoam计算旋转机械（介质是水）的空化问题，边界条件是入口速度fixedValue ，压力totalpressure；出口速度inletOutlet，压力fixedFluxPressure，湍流采用标准K-e模型，额定的转速是303rad/s，入口段直径0.125m，速度给定5m/s，给定的k和epsilon初值分别是0.0284和0.09，时间步长设定1e-5，网格40W左右（非结构），用AMI连接域。在计算几步后就回出现崩溃，但是通过大幅度缩小时间步长和减少转速（大概到50rad）可以正常收敛，但是这两者都不太能接受。所以想寻求各位大神帮忙看看有没有更好的方法能稳定计算，或者如果是我的设置有问题也欢迎指正~~~~附上log和pyfoam的图

Create time Create mesh for time = 0 Selecting dynamicFvMesh dynamicMotionSolverFvMesh Selecting motion solver: solidBody Selecting solid-body motion function rotatingMotion Applying solid body motion to cellZone impeller PIMPLE: No convergence criteria found PIMPLE: No corrector convergence criteria found Calclations will do 3 corrections Reading field p_rgh Reading field U Reading/calculating face flux field phi AMI: Creating addressing and weights between 589 source faces and 754 target faces AMI: Patch source sum(weights) min/max/average = 0.996846, 1, 0.999444 AMI: Patch target sum(weights) min/max/average = 0.99046, 1, 0.999443 AMI: Creating addressing and weights between 2752 source faces and 2836 target faces AMI: Patch source sum(weights) min/max/average = 0.999976, 1.00012, 1.00003 AMI: Patch target sum(weights) min/max/average = 0.99994, 1.00012, 1.00003 Creating phaseChangeTwoPhaseMixture Selecting phaseChange model SchnerrSauer Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model kEpsilon RAS { RASModel kEpsilon; turbulence on; printCoeffs on; Cmu 0.09; C1 1.44; C2 1.92; C3 0; sigmak 1; sigmaEps 1.3; } Reading g Reading hRef Calculating field g.h No finite volume options present Reading/calculating face velocity Uf Courant Number mean: 3.46663e-05 max: 0.00375 Starting time loop Courant Number mean: 3.46663e-05 max: 0.00375 Time = 1e-05 PIMPLE: Iteration 1 AMI: Creating addressing and weights between 589 source faces and 754 target faces AMI: Patch source sum(weights) min/max/average = 0.99691, 1, 0.999455 AMI: Patch target sum(weights) min/max/average = 0.990625, 1, 0.999455 AMI: Creating addressing and weights between 2752 source faces and 2836 target faces AMI: Patch source sum(weights) min/max/average = 0.999971, 1.00012, 1.00004 AMI: Patch target sum(weights) min/max/average = 0.999955, 1.00011, 1.00004 Execution time for mesh.update() = 0.28 s GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0046122, No Iterations 7 time step continuity errors : sum local = 7.2582e-05, global = -6.94718e-05, cumulative = -6.94718e-05 smoothSolver: Solving for alpha.water, Initial residual = 0.999973, Final residual = 1.1325, No Iterations 10 Phase-1 volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 1 Min(alpha.water) = 1 Max(alpha.water) = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00343267, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.00357127, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 4.53115e-06, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 9.27694e-10, No Iterations 48 PIMPLE: Iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 1, Final residual = 4.3514e-09, No Iterations 3 Phase-1 volume fraction = 0.979984 Min(alpha.water) = 0.821031 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.979982 Min(alpha.water) = 0.820777 Max(alpha.water) = 1 smoothSolver: Solving for Ux, Initial residual = 0.152882, Final residual = 0.000456751, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.158655, Final residual = 0.000456774, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.173638, Final residual = 0.000499279, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 0.00387877, Final residual = 9.67937e-10, No Iterations 8 PIMPLE: Iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 0.000217683, Final residual = 5.3264e-11, No Iterations 2 Phase-1 volume fraction = 0.995947 Min(alpha.water) = 0.351358 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.995942 Min(alpha.water) = 0.345137 Max(alpha.water) = 1 smoothSolver: Solving for Ux, Initial residual = 0.0528343, Final residual = 6.27726e-07, No Iterations 3 smoothSolver: Solving for Uy, Initial residual = 0.0567381, Final residual = 7.23559e-07, No Iterations 3 smoothSolver: Solving for Uz, Initial residual = 0.120368, Final residual = 8.2518e-08, No Iterations 4 GAMG: Solving for p_rgh, Initial residual = 3.49149e-06, Final residual = 5.80041e-10, No Iterations 5 smoothSolver: Solving for epsilon, Initial residual = 0.0986499, Final residual = 5.16839e-07, No Iterations 4 smoothSolver: Solving for k, Initial residual = 1, Final residual = 1.29045e-07, No Iterations 6 ExecutionTime = 8.23 s ClockTime = 9 s Courant Number mean: 0.0197059 max: 5.41516 Time = 2e-05 PIMPLE: Iteration 1 AMI: Creating addressing and weights between 589 source faces and 754 target faces AMI: Patch source sum(weights) min/max/average = 0.996998, 1, 0.999472 AMI: Patch target sum(weights) min/max/average = 0.990859, 1, 0.999472 AMI: Creating addressing and weights between 2752 source faces and 2836 target faces AMI: Patch source sum(weights) min/max/average = 0.999977, 1.00014, 1.00004 AMI: Patch target sum(weights) min/max/average = 0.999968, 1.00012, 1.00004 Execution time for mesh.update() = 0.26 s GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00504073, No Iterations 3 time step continuity errors : sum local = 0.000188015, global = -1.15051e-06, cumulative = -7.06224e-05 smoothSolver: Solving for alpha.water, Initial residual = 9.67477e-06, Final residual = 6.18883e-11, No Iterations 2 Phase-1 volume fraction = 0.995772 Min(alpha.water) = 0.211932 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.995769 Min(alpha.water) = 0.174408 Max(alpha.water) = 1 smoothSolver: Solving for Ux, Initial residual = 0.325971, Final residual = 0.0183609, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.309661, Final residual = 0.0165613, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.351942, Final residual = 0.0299024, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 1.43572e-07, Final residual = 7.23286e-10, No Iterations 5 PIMPLE: Iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 0.0155352, Final residual = 2.2274e-09, No Iterations 3 Phase-1 volume fraction = 0.952704 Min(alpha.water) = 0.277133 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.952698 Min(alpha.water) = 0.228159 Max(alpha.water) = 1.00002 smoothSolver: Solving for Ux, Initial residual = 0.104737, Final residual = 0.00347544, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.0999211, Final residual = 0.00465959, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.18458, Final residual = 0.0128277, No Iterations 2 GAMG: Solving for p_rgh, Initial residual = 0.0234743, Final residual = 7.77366e-10, No Iterations 5 PIMPLE: Iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 0.00178812, Final residual = 5.32355e-09, No Iterations 2 Phase-1 volume fraction = 0.985629 Min(alpha.water) = 0.0661184 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.985616 Min(alpha.water) = 0.0661918 Max(alpha.water) = 1.00001 smoothSolver: Solving for Ux, Initial residual = 0.0821993, Final residual = 9.95813e-07, No Iterations 383 smoothSolver: Solving for Uy, Initial residual = 0.092795, Final residual = 9.84964e-07, No Iterations 401 smoothSolver: Solving for Uz, Initial residual = 0.174958, Final residual = 9.94153e-07, No Iterations 455 GAMG: Solving for p_rgh, Initial residual = 0.00206972, Final residual = 3.9387e-10, No Iterations 6 smoothSolver: Solving for epsilon, Initial residual = 0.393518, Final residual = 8.58378e-07, No Iterations 27 smoothSolver: Solving for k, Initial residual = 0.951336, Final residual = 9.76279e-07, No Iterations 374 ExecutionTime = 24.93 s ClockTime = 25 s Courant Number mean: 0.0203512 max: 31.9041 Time = 3e-05 PIMPLE: Iteration 1 AMI: Creating addressing and weights between 589 source faces and 754 target faces AMI: Patch source sum(weights) min/max/average = 0.997111, 1, 0.999493 AMI: Patch target sum(weights) min/max/average = 0.991164, 1, 0.999493 AMI: Creating addressing and weights between 2752 source faces and 2836 target faces AMI: Patch source sum(weights) min/max/average = 0.999987, 1.00015, 1.00004 AMI: Patch target sum(weights) min/max/average = 0.999955, 1.00014, 1.00004 Execution time for mesh.update() = 0.22 s GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00576691, No Iterations 3 time step continuity errors : sum local = 0.000687359, global = -2.57908e-06, cumulative = -7.32014e-05 smoothSolver: Solving for alpha.water, Initial residual = 0.000155352, Final residual = 4.58115e-09, No Iterations 3 Phase-1 volume fraction = 0.985431 Min(alpha.water) = 0.00232962 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.985415 Min(alpha.water) = -0.0498938 Max(alpha.water) = 1.00005 smoothSolver: Solving for Ux, Initial residual = 0.556211, Final residual = 0.0456756, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.483701, Final residual = 0.0456852, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.54432, Final residual = 0.0380588, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 0.000953864, Final residual = 3.46933e-10, No Iterations 7 PIMPLE: Iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 7.33141e-05, Final residual = 3.27551e-09, No Iterations 3 Phase-1 volume fraction = 0.985373 Min(alpha.water) = 0.00519769 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.985357 Min(alpha.water) = -0.0226666 Max(alpha.water) = 1.00002 smoothSolver: Solving for Ux, Initial residual = 0.251887, Final residual = 0.0202833, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.190243, Final residual = 0.0150066, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.216529, Final residual = 0.0209459, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 2.22683e-06, Final residual = 4.46301e-10, No Iterations 1 PIMPLE: Iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 5.92958e-05, Final residual = 4.78988e-09, No Iterations 3 Phase-1 volume fraction = 0.985322 Min(alpha.water) = 0.00340138 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.985304 Min(alpha.water) = -0.0173934 Max(alpha.water) = 1.00005 smoothSolver: Solving for Ux, Initial residual = 0.300509, Final residual = 0.0025018, No Iterations 1000 smoothSolver: Solving for Uy, Initial residual = 0.398457, Final residual = 0.00194781, No Iterations 1000 smoothSolver: Solving for Uz, Initial residual = 0.353341, Final residual = 0.00647722, No Iterations 1000 GAMG: Solving for p_rgh, Initial residual = 2.82287e-06, Final residual = 4.84489e-10, No Iterations 3 smoothSolver: Solving for epsilon, Initial residual = 0.244345, Final residual = 8.33659e-07, No Iterations 29 smoothSolver: Solving for k, Initial residual = 0.933459, Final residual = 0.00774907, No Iterations 1000 ExecutionTime = 56.97 s ClockTime = 58 s Courant Number mean: 0.039008 max: 94.2267 Time = 4e-05 PIMPLE: Iteration 1 AMI: Creating addressing and weights between 589 source faces and 754 target faces AMI: Patch source sum(weights) min/max/average = 0.997235, 1, 0.999519 AMI: Patch target sum(weights) min/max/average = 0.99154, 1, 0.999518 AMI: Creating addressing and weights between 2752 source faces and 2836 target faces AMI: Patch source sum(weights) min/max/average = 0.999987, 1.00015, 1.00005 AMI: Patch target sum(weights) min/max/average = 0.999955, 1.00015, 1.00005 Execution time for mesh.update() = 0.28 s GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.00975609, No Iterations 5 time step continuity errors : sum local = 0.00591848, global = 3.92436e-06, cumulative = -6.92771e-05 smoothSolver: Solving for alpha.water, Initial residual = 0.000372595, Final residual = 7.54046e-09, No Iterations 8 Phase-1 volume fraction = 0.984937 Min(alpha.water) = 0.00110965 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.984913 Min(alpha.water) = -2.29641 Max(alpha.water) = 2.59561 smoothSolver: Solving for Ux, Initial residual = 0.541575, Final residual = 0.0121456, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.62719, Final residual = 0.0557338, No Iterations 1 smoothSolver: Solving for Uz, Initial residual = 0.462147, Final residual = 0.0348564, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 3.06236e-05, Final residual = 2.15551e-10, No Iterations 2 PIMPLE: Iteration 2 smoothSolver: Solving for alpha.water, Initial residual = 0.00637181, Final residual = 5.1554e-08, No Iterations 10 Phase-1 volume fraction = 0.984205 Min(alpha.water) = 1.05723e-20 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.982971 Min(alpha.water) = -1336.22 Max(alpha.water) = 97.6268 smoothSolver: Solving for Ux, Initial residual = 0.351934, Final residual = 0.0239104, No Iterations 1 smoothSolver: Solving for Uy, Initial residual = 0.318012, Final residual = 0.0115367, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.337556, Final residual = 0.03052, No Iterations 1 GAMG: Solving for p_rgh, Initial residual = 4.60749e-08, Final residual = 3.66729e-10, No Iterations 1 PIMPLE: Iteration 3 smoothSolver: Solving for alpha.water, Initial residual = 0.998656, Final residual = 9.85984e-09, No Iterations 7 Phase-1 volume fraction = 0.983918 Min(alpha.water) = -0.000242279 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Liquid phase volume fraction = 0.984043 Min(alpha.water) = -16.1095 Max(alpha.water) = 105.668 [1] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [1] #1 Foam::sigFpe::sigHandler(int) at ??:? [1] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" [1] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? [1] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:? [1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:? [1] #6 ? at ??:? [1] #7 ? at ??:? [1] #8 ? at ??:? [1] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" [1] #10 ? at ??:?

1_1537536681529_residual.PNG 0_1537536681528_continuity.PNG