关于稳态欧拉-拉格朗日计算速度的探讨
-
目前我在采用稳态的欧拉拉格朗日求解器,求解水轮机的磨损问题,颗粒数为50000,网格数为1100w左右。颗粒文件的设置如下:```
code_textsolution { active true; calcFrequency 15; maxTrackTime 6; maxCo 0.6; coupled true; cellValueSourceCorrection off; //suggested to be off interpolationSchemes { Urel cellPoint; //U for water rho cell; //rho of water mu cell; //mu of water DUcDt cell; //DUcDt of water } integrationSchemes { U Euler; } sourceTerms { resetOnStartup false; schemes { U semiImplicit 1; } } } constantProperties { rho0 2650; //the density for the particle } subModels { particleForces { sphereDrag; //drag forces gravity; // gravity force pressureGradient //pressure gradient force { U Urel; } virtualMass //virtual mass force { Cvm 0.5; } MRF { omega (0 0 -104.72); centreOfRotation (0 0 0); Ctemp MRFC; } } injectionModels { model1 { type patchInjection; parcelBasisType fixed; patchName INLET_SV; U0 (0 0 0); nParticle 420; duration 1; parcelsPerSecond 50000; flowRateProfile constant 1; massFlowRate 0; //positionsFile "kinematicCloudPositions"; sizeDistribution { type fixedValue; fixedValueDistribution { value 0.0002;//m } } } } dispersionModel stochasticDispersionRAS; patchInteractionModel UserLocalInteraction; //using the tabakoff rebound model UserLocalInteractionCoeffs { patches ( "(.*)" { type rebound; reboundModel 2; } OUTLET_SC { type none; } INLET_SV { type escape; } OUTLET_SV { type none; } INLET_GV { type none; } OUTLET_GV { type none; } INLET_RUN { type none; } OUTLET_RUN { type none; } INLET_DT { type none; } INLET_SC { type none; } OUTLET_DT { type escape; } ); } heatTransferModel none; surfaceFilmModel none; collisionModel none; stochasticCollisionModel none; radiation off; }
流场迭代15步之后,会进行颗粒的迭代,颗粒计算平均耗时5min。
但每次计算至210步时,颗粒就算不动了(前面的颗粒计算均是正常的)。
我如果把颗粒数降至10000,计算至750步,颗粒就算不动了。请问一下大家又碰到这类问题吗?有什么比较好的解决办法吗?