CFD中文网

李博的桌面？

@李东岳 李老师，您好，我按照您的建议进行了相关的测试。结果发现当更换为GAMG时计算比PCG慢非常多，如果我在计算一开始取残差与相对残差值较大时，完全不收敛的情况基本上到500步左右就改善了。但是我又发现一个特别奇怪的现象：
首先通过pimpleFoam (nOuterCorrectors=1) 计算到一定步时计算结果收敛很好，残差基本下降至1e-5左右；但是当我在运行时更改nOuterCorrectors=2时，初始残差立即上升至0.4左右，计算异常缓慢，并且继续运行算例残差也下降不了。请问这个可能是什么原因导致的呢？

@bestucan 厉害，第一次知道这种用法！:146:

@是嘿嘿牛啊 密度是被动地，是根据流体种类、温度和压力变化的，你可以设置上下两层温差或不同的流体种类

计算的是气液混输泵的气液两相流计算，air-water(25℃)的混合介质，air为离散相给定体积分数和固定的气泡直径，water为连续项，湍流模型为SST k-w，欧拉双流体模型。计算出来之后总感觉残差是不是波动有点大了。
有没有这方面的大佬教教我，是不是设置有问题还是其他的问题。下面是计算得出的残差
47a7ac58-69fe-4db7-a71d-e8ba6fd07ebc-image.png
800b8c63-4b22-47a5-9131-cc03980b1291-image.png
05ca6207-15dd-45ab-bb5e-9628346b8a9e-image.png
e252f689-2fd6-4d82-b46a-77c7453a624a-image.png

想解决IV型激波，入射激波如何通过边界条件设置产生

@东岳 老师，没理解您的意思。后面的$\alpha_q\lambda_q\nabla·\vec{v_q}I$中单位二阶张量是表示体形变率仅出现在法相应力中，红框内等于$-\frac{2}{3}\lambda_q$时才能与法相应力对应啊

导出时候选择了如下这些选项
image.png

udf有问题，你试试这个

#include "udf.h" DEFINE_PROFILE(inlet_velocity, t, i) { real x[ND_ND],y,vel; face_t f; begin_f_loop(f,t) { F_CENTROID(x, f, t); y = x[1]; if(y>0.0393) vel = 29.5; else vel = 29.5 * pow((y / 0.0393), 1/7); F_PROFILE(f, t, i) = vel; } end_f_loop(f,t) }

moment method对于nucleation和condensation的模拟简单吗？

即为简单。但是溶解/蒸发数值上要复杂些。2017年我在HZDR的时候同事用fixed pivot method模拟增长，我用moment method，相比来看moment method要更简单，并且不需要处理相空间的drift：https://www.sciencedirect.com/science/article/pii/S1004954118301174/pdfft?md5=84216a8842ce3fd6d1a99d4f327c1ee1&pid=1-s2.0-S1004954118301174-main.pdf

@HITSC30 icem不是很熟练，我去试一下

上传一下我的snappyhexMeshDict文件，请各位老师帮忙看看```
code_text

========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 10 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { format ascii; class dictionary; object snappyHexMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #includeEtc "caseDicts/mesh/generation/snappyHexMeshDict.cfg" castellatedMesh true; snap true; addLayers true; geometry { cylinder { type triSurfaceMesh; file "cylinder.stl"; } refinementBox { type searchableBox; min (0.0 0.0 0.08); max (1.2 0.2 0.22); } }; castellatedMeshControls { maxLocalCells 200000; maxGlobalCells 5000000; minRefinementCells 10; maxLoadUnbalance 0.10; nCellsBetweenLevels 6; features ( { file "cylinder.eMesh"; level 3; } ); refinementSurfaces { cylinder { level (3 3); patchInfo { type wall; } } } refinementRegions { refinementBox { mode inside; level 1; } } insidePoint (0.4 0.1 0.15); } snapControls { nSmoothPatch 3; tolerance 2.0; nSolveIter 30; nRelaxIter 5; nFeatureSnapIter 10; explicitFeatureSnap true; implicitFeatureSnap false; multiRegionFeatureSnap false; } addLayersControls { layers { "(cylinder_stl_face<stlunit=MM>|cylinder_stl_top<stlunit=MM>)" { nSurfaceLayers 10; } } relativeSizes true; expansionRatio 1.05; finalLayerThickness 0.8; minThickness 0.5; } meshQualityControls {} writeFlags ( // scalarLevels // layerSets // layerFields ); mergeTolerance 1e-6; // ************************************************************************* //

@xpqiu
.................
25T硬盘确实花了不少钱。:sunglasses:

您好，请问这个问题您是否有好的解决办法了吗？
强制设为0可能不太对，反映出来的可能不是压实效果。。

@gemini 谢谢您，服务fluent版本跟调试用的是一致的，我去试试修改成您说的这样

退Fortran用c++保平安，虽然没用过Fortran，但看知乎上专门论述Fortran的缺陷，发现用c++一下就好了。但是c++的学习周期长，东西多，也是问题。

@东岳 谢谢！

@Wayne

我已联系你:sunglasses:

@ibelief 在 颗粒表面化学反应的文献 中说：

基于简单碰撞理论煤粉燃烧动力学模型的研究

非常感谢

各位老师好！
报错信息：
Screenshot from 2021-11-20 14-31-59.png

solidParticleCloud.C的代码：

solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); Cloud<solidParticle>::move(*this, td, mesh_.time().deltaTValue()); if(mesh_.time().value()> td.cloud().tInjStart_ && // mesh_.time().value()< td.cloud().tInjEnd_) {this->inject(td);} void Foam::solidParticleCloud::inject(solidParticle::trackingData &td) //. { //. label cellI = -1; //. label tetFaceI = -1; label tetPtI = -1; mesh_.findCellFacePt(td.cloud().posB_, cellI, tetFaceI, tetPtI); if(cellI > 0){ solidParticle* ptr1 = new solidParticle(mesh_, td.cloud().posB_, cellI, tetFaceI, tetPtI,td.cloud().dB_, td.cloud().UB_); Cloud<solidParticle>::addParticle(ptr1); } }