摘要：

      舰船自航评估和研究通常是由拖曳水池自航试验实现，实施周期长，耗资大。船舶三维多面体OPENFOAM

潜艇外流场分析报告

项目概述

该项目展示了如何模拟潜艇周围的外部流体流动的示例

舰船自航评估和研究通常是由拖曳水池自航试验实现，实施周期长，耗资大。采用流体动力数值仿真分析技术，能够系统分析和研究各种水下航行器的航行性能，指导水下航行器设计，提出优化设计方案和改进航行性能的措施。

本项目对潜艇的水动力航行性能进行了分析和研究。

模型简化

本算例使用几何来源simscale官网，几何模型如图1所示。

1. 潜艇几何简化模型

网格设置

本算例采用snappyHexMesh对几何模型进行网格划分（如图2），详细参数设置如下所示：

geometry

{solid_0

{type triSurfaceMesh;

file "solid_0.stl";

}

innerbox1

{type searchableBox;

min (-10 -15 -10);

max ( 80 15 15);

}

innerbox2

{type searchableBox;

min (-75 -15 -10);

max ( 90 15 15);

}

};

castellatedMeshControls

{maxLocalCells 2000000;

maxGlobalCells 4000000;

minRefinementCells 1;

maxLoadUnbalance 0.2;

nCellsBetweenLevels 4;

resolveFeatureAngle 30;

allowFreeStandingZoneFaces false;

features

(

{file "solid_0.extendedFeatureEdgeMesh";

level 4;

}

);

refinementSurfaces

{solid_0

{level (4 4);

patchInfo

{ type wall;

inGroups (solid_0Group);

}

}

}

refinementRegions

{solid_0

{mode distance;

levels ((1 4));

}

innerbox1

{mode inside;

levels ((1E15 2));

}

innerbox2

{mode inside;

levels ((1E15 1));

}

}

locationInMesh (-120 0 0);

}

snapControls

{nSmoothPatch 3;

tolerance 3.0;

nSolveIter 300;

nRelaxIter 5;

nFeatureSnapIter 10;

implicitFeatureSnap true;

explicitFeatureSnap false;

multiRegionFeatureSnap false;

}

addLayersControls

{relativeSizes true;

layers

{

"solid_0_.*"

{nSurfaceLayers 5;

}

}

expansionRatio 1.3;

finalLayerThickness 0.3;

minThickness 0.001;

nGrow 0;

featureAngle 180;

slipFeatureAngle 75;

nRelaxIter 5;

nSmoothSurfaceNormals 1;

nSmoothNormals 3;

nSmoothThickness 10;

maxFaceThicknessRatio 0.5;

maxThicknessToMedialRatio 0.3;

minMedianAxisAngle 90;

nBufferCellsNoExtrude 0;

nLayerIter 50;

nRelaxedIter 20;

}

meshQualityControls

{maxNonOrtho 65;

maxBoundarySkewness 20;

maxInternalSkewness 4;

maxConcave 80;

minVol 1e-13;

minTetQuality 1e-15;

minArea -1;

minTwist 0.01;

minDeterminant 0.001;

minFaceWeight 0.02;

minVolRatio 0.01;

minTriangleTwist -1;

nSmoothScale 4;

errorReduction 0.75;

relaxed

{maxNonOrtho 75;

}

}

writeFlags

(scalarLevels // write volScalarField with cellLevel for postprocessing

layerSets // write cellSets, faceSets of faces in layer

layerFields // write volScalarField for layer coverage

);

mergeTolerance 1E-6;

1. 潜艇网格划分

2. 网格信息参数及网格质量参数如下所示：

3. 网格信息参数

网格总数

---

数量 1544724 4346105 1402876 
网格类型 
类型 hexahedra prisms tet wedges polyhedra 数量 1280641 21250 25 100957 网格质量 
评价指标 最大纵横比 最小体积 最大非正交性 最大歪斜率 值 19.795 6.38e-005 64.11 5.33

边界条件参数

该算例湍流模型选用K-Omega SST模型。外流场壁面应用对称平面边界条件，进出口边界条件见表1。

表1 边界条件参数

边界条件 参数

---

速度进口 15m/s 压力出口 0Pa

求解器参数设置

fvSchemes设置

ddtSchemes {

default steadyState;

}

gradSchemes {

default Gauss linear;

grad(U) Gauss linear;

grad(p) Gauss linear;

}

divSchemes {

default none;

div(phi,U) bounded Gauss upwind;

div(phi,omega) bounded Gauss upwind;

div((nuEff*dev2(T(grad(U))))) Gauss linear;

div(phi,k) bounded Gauss upwind;

}

laplacianSchemes {

default Gauss linear corrected;

laplacian(nuEff,U) Gauss linear corrected;

laplacian((1|A(U)),p) Gauss linear corrected;

laplacian(nu,U) Gauss linear corrected;

}

interpolationSchemes {

default linear;

interpolate(HbyA) linear;

}

snGradSchemes {

default corrected;

}

fluxRequired {

default no;

p ;

}

wallDist

{method meshWave;

}

fvSolution设置

solvers {

p {

relTol 0.005;

cacheAgglomeration on;

nPreSweeps 0;

agglomerator faceAreaPair;

nCellsInCoarsestLevel 10;

solver GAMG;

nPostSweeps 0;

mergeLevels 1;

smoother GaussSeidel;

tolerance 1e-06;

}

k {

relTol 0.01;

preconditioner DILU;

tolerance 1e-06;

solver PBiCG;

}

U {

relTol 0.01;

preconditioner DILU;

tolerance 1e-05;

solver PBiCG;

}

omega {

relTol 0.01;

preconditioner DILU;

tolerance 1e-06;

solver PBiCG;

}

}

relaxationFactors {

fields {

p 0.3;

}

equations {

omega 0.7;

k 0.7;

U 0.7;

}

}

SIMPLE {

pRefValue 0.0;

residualControl {

p 1e-06;

k 1e-06;

U 1e-06;

omega 1e-06;

}

nNonOrthogonalCorrectors 0;

pRefCell 0;

}

controlDict设置

application simpleFoam;

startFrom startTime;

stopAt endTime;

purgeWrite 0;

writeFormat ascii;

writePrecision 12;

writeCompression on;

timeFormat general;

timePrecision 12;

runTimeModifiable yes;

deltaT 1.0;

startTime 0.0;

writeInterval 1000;

adjustTimeStep no;

endTime 2000.0;

writeControl timeStep;

functions {

}

计算结果

本算例计算总步长为1272步时云图如下所示：