foam-extend求解稳态问题,如何才能使最终解残差下降
-
各位朋友,最近我在计算一个稳态的Laval喷管的稳态问题。可是最后的马赫盘的位置一直在非常微小幅度的振荡,导致最后的残差一直居高不下。我尝试过如下的几种方法:
1。 waveTransmissive改为稳态适用版本;
2。 减少co数和松弛因子。
3。 加密网格。我用的是kOmegaSST, 求解器比较小众,是transonicMRFDyMFoam. 请问各位在稳定稳态求解方面有什么好的经验?
谢谢各位! -
为了方便,我把我的fvScheme贴在这儿:
ddtSchemes { default none;//EulerLocal physDeltaT CoDeltaT; ddt(rho) EulerLocal physDeltaT CoDeltaT; ddt(rhoU) EulerLocal physDeltaT CoDeltaT; ddt(rhoE) EulerLocal physDeltaT CoDeltaT; ddt(rho,omega) EulerLocal physDeltaT CoDeltaT; ddt(rho,k) EulerLocal physDeltaT CoDeltaT; p steadyState; // Needed for waveTransmissive; phiMin steadyState; // Needed for waveTransmissive; phiMax steadyState; // Needed for waveTransmissive; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(rhoSlope) Gauss linear; grad(pSlope) Gauss linear; grad(USlope) Gauss linear; grad(TSlope) Gauss linear; snGradCorr(U) Gauss linear; snGradCorr(h) Gauss linear; snGradCorr(omega) Gauss upwind phi; snGradCorr(k) Gauss upwind phi; grad(omega) Gauss upwind phi; grad(k) Gauss upwind phi; } divSchemes { default none; div(phi,U) Gauss linearUpwind grad(U); div(phi,k) Gauss QUICK;//linearUpwind; div(tau) Gauss linear; div((tau&U)) Gauss linear; div(phi,omega) Gauss QUICK;//linearUpwind; } laplacianSchemes { default Gauss linear limited 0.33; } interpolationSchemes { default none; interpolate(rho) linear corrected; } snGradSchemes { default orthogonal; } fluxRequired { default no; }
2017年7月6日 00:21
1/3
2017年7月6日 08:33