翼型摆动cfd模拟Word文档下载推荐.docx

1、表2翼型贴体网格信息参数4物性参数分析所涉及流场介质为空气，其相关物性参数如表3所示。表3物性参数5边界条件自由流的流速（Air）设定为102m / s，分析对象为低马赫数，空气选用不可压缩介质，密度为1kg/m3，湍流模型选用kOmegaSST模型，初始边界条件设置如下表4所示：表4初始边界条件设置6、openfoam求解器设置本项目为求解摆动翼型流场，湍流模型选用kOmegaSST，需分别设置对应fvSchemes离散方法，fvSolution方程求解方法，指定场函数setFieldsDict ，动网格运动参数及求解控制参数。1.1 离散方法fvSchemes设置ddtSchemes /时

2、间离散格式，该项目瞬态计算采用欧拉离散Eulerdefault Euler;gradSchemes/梯度离散，采用高斯方法，有界线性插值，为二阶离散 cellLimited Gauss linear 1;divSchemes/散度离散，对流项U采用带限制器高斯线性插值，为二阶离散，k、epsilon采用高斯迎风格式，为一阶离散。default Gauss linear;div（rhoPhi,U） Gauss linearUpwindV grad（U）;div（phi,alpha） Gauss vanLeer;div（phirb,alpha） Gauss linear;div（phi,k） Ga

3、uss upwind;div（phi,epsilon） Gauss upwind;div（rho*nuEff）*dev2（T（grad（U） Gauss linear;laplacianSchemes/拉普拉斯项离散，扩散项及压力方程离散均采用高斯理论线性插值，并带有正交修正default Gauss linear corrected;interpolationSchemes/插值格式，默认线性插值default linear;snGradSchemes/梯度法向分量，默认带有非正交修正default limited 1;oversetInterpolation/嵌套网格插值格式 method

4、 inverseDistance;wallDist method meshWave;1.2 方程求解方法fvSolution设置solvers cellDisplacement solver PCG; preconditioner DIC; tolerance 1e-06; relTol 0; maxIter 100; p PBiCGStab; DILU; 1e-6; 0.01; minIter 2; pFinal （U|k|epsilon|omega） 1e-8;0.0;（U|k|epsilon|omega）Final $U;PIMPLE momentumPredictor true; co

5、rrectPhi false; oversetAdjustPhi nOuterCorrectors 1; nCorrectors 3; nNonOrthogonalCorrectors 2; ddtCorr checkMeshCourantNo yes;relaxationFactors fields equations.*1.3 域设置defaultFieldValues（ volScalarFieldValue zoneID 123）;regions cellToCell set c0; fieldValues （ volScalarFieldValue zoneID 0 ）; set c

6、1; volScalarFieldValue zoneID 11.4 动网格参数设置dynamicFvMesh dynamicOversetFvMesh;dynamicOversetFvMeshCoeffssolver multiSolidBodyMotionSolver;multiSolidBodyMotionSolverCoeffs movingZone solidBodyMotionFunction tabulated6DoFMotion; CofG （0.33 -0.5 0）; timeDataFileName $FOAM_CASE/constant/6DoF.dat;1.5 求解控制

7、参数设置application pimpleDyMFoam;startFrom startTime;startTimestopAt endTime;endTime 10;deltaTwriteControl runTime;writeInterval 0.2;purgeWritewriteFormat ascii;writePrecisionwriteCompression off;timeFormat general;timePrecision 6;runTimeModifiable true;adjustTimeStepmaxCo 4;maxDeltaTfunctions7、结果分析 图4. 残差收敛曲线10s内各时间段速度场如下图所示：图5. 10s内翼型摆动速度云图10s内各时间段压力场如下图所示：图6. 10s内翼型摆动压力云图