CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    关于稳态欧拉-拉格朗日计算速度的探讨

    OpenFOAM
    2
    4
    783
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 液
      液固两相流 最后由 编辑

      目前我在采用稳态的欧拉拉格朗日求解器,求解水轮机的磨损问题,颗粒数为50000,网格数为1100w左右。颗粒文件的设置如下:```
      code_text

      solution
      {
          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步,颗粒就算不动了。

      请问一下大家又碰到这类问题吗?有什么比较好的解决办法吗?

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        发代表性网格截面看一下

        线上CFD课程开始报名:http://www.dyfluid.com/class.html

        CFD高性能服务器 http://dyfluid.com/servers.html

        液 1 条回复 最后回复 回复 引用
        • 液
          液固两相流 @李东岳 最后由 编辑

          @东岳
          我的网格有多个交界面(用的是cyclicAMI),交界面对于拉格朗日的计算会有影响吗?

          1 条回复 最后回复 回复 引用
          • 李东岳
            李东岳 管理员 最后由 编辑

            很有可能,cyclicAMI有可能把粒子算丢了

            线上CFD课程开始报名:http://www.dyfluid.com/class.html

            CFD高性能服务器 http://dyfluid.com/servers.html

            1 条回复 最后回复 回复 引用
            • First post
              Last post