自由来流下的网框表面流速异常
-
各位老师好，
【问题描述】
本人在使用pimpleFoam模拟固定式网框在自由来流（u=0.5m/s）中的响应。整个流场的计算结果符合预期，可是网框表面处的流速竟然高达1e2量级，明显与现实不符，结果云图如下：（侧视图）
【分析】
网框截面原本比本帖中展示的要更小，起初正方形边长仅分布四个网格，为了排除网格数量对结果的影响，特意扩大了截面，增加边长网格数量。但是现在依然存在这个问题。【初始条件设置】
以下摘抄0文件夹设置：
p：/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "../varsSettings" //Include global variables dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value $internalField; } frontAndBack { type zeroGradient; } topAndBottom { type zeroGradient; } frame { type fixedFluxPressure; } } // ************************************************************************* //U:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volVectorField; location "0"; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "../varsSettings" //Include global variables dimensions [0 1 -1 0 0 0 0]; internalField uniform $flowVelocity; boundaryField { inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; inletValue uniform (0 0 0); value $internalField; } frontAndBack { type slip; value uniform (0 0 0); } topAndBottom { type slip; value uniform (0 0 0); } frame { /* type noSlip; */ type movingWallVelocity; value uniform (0 0 0); } } // ************************************************************************* //pointDisplacement:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 5 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class pointVectorField; object pointDisplacement; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 0 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type fixedValue; value uniform (0 0 0); } outlet { type fixedValue; value uniform (0 0 0); } frontAndBack { type fixedValue; value uniform (0 0 0); } topAndBottom { type fixedValue; value uniform (0 0 0); } frame { type calculated; value uniform (0 0 0); } } // ************************************************************************* //【算例设置】
functionalTest_NoPZRefine1.zip
生成网格createMesh，
跑算例runCaseOnly，
全流程runCase希望各位老师能给出建议~@李东岳
-
@李东岳 静态网格也是有这个问题的。按理说不应该啊...我再跑跑其他类似的原生算例
p：
#include "../varsSettings" //Include global variables dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type zeroGradient; } outlet { type fixedValue; value $internalField; } frontAndBack { type empty;//zeroGradient; } topAndBottom { type zeroGradient; } frame { type zeroGradient;//fixedFluxPressure; } }U：
#include "../varsSettings" //Include global variables dimensions [0 1 -1 0 0 0 0]; internalField uniform $flowVelocity; boundaryField { inlet { type fixedValue; value $internalField; } outlet { type inletOutlet; inletValue uniform (0 0 0); value $internalField; } frontAndBack { type empty;//slip; //value uniform (0 0 0); } topAndBottom { type slip; value uniform (0 0 0); } frame { type noSlip;//movingWallVelocity; //value uniform (0 0 0); } }算例文件 OpenfOAM2012 网格数4w
functionalTest_NoPZRefine2_2D_staticMesh.zip -
~