CFD中文网

如果局部网格尺寸为1e-5m时，由于vof中库朗数的限制，时间步长也需要设置的很小，为1e-6或者1e-7，此时计算时间过长，该如何解决呢？

@星星星星晴 好的 我搞了一个

三维的类似情况，我最后一个打错字了。

teoria.png
如图所示，这是非常经典的Le Méhauté (1976)波浪理论适用图，但我有点迷糊，主要有两方面疑问：

图中红色部分看起来是Stokes三阶波/四阶波与五阶波的交接处，那么我应该选择用三阶波/四阶波还是五阶波呢？同样地，图中浅蓝色和五阶波也有交集，那么在既是浅蓝色又是五阶波的区域应该用cnoidal还是五阶波呢？ 我还注意到图中五阶波下面还有一个流函数理论，但我看olaFlow的用户手册，用流函数的话需要输入一大堆参数（如下图所示），这些参数应该怎么选取呢？本菜鸡在此请教一下各位:135:
2024-07-08 21-55-00 的屏幕截图.png
2024-07-08 22-01-59 的屏幕截图.png

上传一下我的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; // ************************************************************************* //

可以把正确的UDF发来参考一下吗？我也做类似的复现研究，谢谢！

畸形波无论是在fluent还是在openfoam中模拟方法都未有明确作法，请大家指导一下模拟方法。

居然还发公众号了

基于OpenFOAM开发了一个简单的程序，做的是气液两相流模拟

看起来更符合 international journal of multiphase flow。不过这俩个都是好期刊，能中哪个都行。在我眼里没什么区别。

已经收到了

https://www.cfd-online.com/Forums/openfoam-meshing/133850-snappyhexmesh-problems.html

@李东岳 安装ansys 时就可以同时安装cfx么？按照我经验，目前的ansys学术版里面可以直接安装cfx，不需要任何额外的操作。

没有低版本的openfoam，你可以试下下面的方法。
理论上来说在applications/solvers/lagrangian/DPMFoam/DPMTurbulenceModels/DPMTurbulenceModels.C中加上其他湍流模型应该就可以，比如：

#define makeLESModel(Type) \ makeTemplatedTurbulenceModel \ (singlePhaseTransportModelPhaseIncompressibleTurbulenceModel, LES, Type) #include "kEpsilon.H" makeRASModel(kEpsilon); #include "realizableKE.H" makeRASModel(realizableKE); #include "Smagorinsky.H" makeLESModel(Smagorinsky);

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

我只是有空看下论坛，随缘回复一下我能帮助的问题。平时没什么时间。:135:

在pEqn.H中对U进行更新时有

const surfaceScalarField phig(-rhorAUf*ghf*fvc::snGrad(rho)*mesh.magSf()); U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf);

我想请问其中的p_rghEqn.flux()代表什么意思，以及为什么phig中要先乘rhorAUf再在U中除以rhorAUf；
以及之后的根据密度差判断收敛性中，我发现compressibleContinuityErrs.H代码如下，只计算了误差值，但没有判断是否符合要求，想请问各位大佬是在哪里判断的呢

dimensionedScalar totalMass = fvc::domainIntegrate(rho); scalar sumLocalContErr = (fvc::domainIntegrate(mag(rho - thermo.rho()))/totalMass).value(); scalar globalContErr = (fvc::domainIntegrate(rho - thermo.rho())/totalMass).value(); cumulativeContErr += globalContErr; Info<< "time step continuity errors : sum local = " << sumLocalContErr << ", global = " << globalContErr << ", cumulative = " << cumulativeContErr << endl; }

谢谢大家！

@AppleKiller 感谢，我再看看，以前也接触过这些资料，感觉用欧拉法求解器为基础修改的拉格朗日求解器不大好，没怎么深入去学习。我再看看这些思路，感谢

any official links?

我最近出现了这个情况。还有人有这个情况么，很奇怪