Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新
    1. Home
    2. lv1995
    L
    • Profile
    • Following
    • Followers
    • Topics
    • Posts
    • Groups

    lv1995

    @lv1995

    30
    Posts
    803
    Profile views
    0
    Followers
    0
    Following
    Joined Last Online

    lv1995 Follow

    Best posts made by lv1995

    This user hasn't posted anything yet.

    Latest posts made by lv1995

    • overPimpleDyMFoam计算中出现报错

      在使用overPimpleDyMFoam计算旋转机械 时,出现了如下报错:

      Create time
      
      Create mesh for time = 0
      
      Selecting dynamicFvMesh dynamicOversetFvMesh
      Selecting motion solver: rigidBodyMotion
      Selecting rigidBodySolver Newmark
      
      PIMPLE: Operating solver in PISO mode
      
      Reading field p
      
      Reading field U
      
      Reading/calculating face flux field phi
      
      Creating cellMask field to block out hole cells
      
      inverseDistance : detected 2 mesh regions
          zone:0 nCells:29272  voxels:(193 193 1) bb:(-300.001077033 -200.001077033 1.99892296657) (700.001077033 200.001077033 3.00107703343)
          zone:1 nCells:8145  voxels:(193 193 1) bb:(-50.0001414007 -49.9939730307 1.99985859933) (50.0001414007 49.9720506607 3.00014140067)
      --> FOAM Warning : 
          From function bool Foam::oversetPolyPatch::master() const
          in file oversetPolyPatch/oversetPolyPatch.C at line 149
          The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation.
      Overset analysis : nCells : 37417
          calculated   : 37186
          interpolated : 204 (interpolated from local:204  mixed local/remote:0  remote:0)
          hole         : 27
      
      Creating interpolatedCells field 
      
      Selecting incompressible transport model Newtonian
      Selecting turbulence model type RAS
      Selecting RAS turbulence model kOmegaSST
      Selecting patchDistMethod meshWave
      RAS
      {
          RASModel        kOmegaSST;
          turbulence      on;
          printCoeffs     on;
          alphaK1         0.85;
          alphaK2         1;
          alphaOmega1     0.5;
          alphaOmega2     0.856;
          gamma1          0.555555555556;
          gamma2          0.44;
          beta1           0.075;
          beta2           0.0828;
          betaStar        0.09;
          a1              0.31;
          b1              1;
          c1              10;
          F3              false;
          decayControl    false;
          kInf            0;
          omegaInf        0;
      }
      
      Reading/calculating face velocity Uf
      
      No MRF models present
      
      No finite volume options present
      Courant Number mean: 0.0195783147649 max: 21.9831544632
      
      Starting time loop
      
      Courant Number mean: 0.0195783147649 max: 21.9831544632
      Time = 0.005
      
      forces forces:
          Not including porosity effects
      
      
      
      --> FOAM FATAL ERROR: 
      Could not find rho:rho
      
          From function void Foam::functionObjects::forces::initialise()
          in file forces/forces.C at line 204.
      
      FOAM exiting
      

      使用的动网格文件如下:

      motionSolverLibs    ("librigidBodyMeshMotion.so");
      
      dynamicFvMesh       dynamicOversetFvMesh;
      
      motionSolver        rigidBodyMotion;
      
      // rigidBodyMotionCoeffs  - not needed any more
      //  {
          report          on;
      
          solver
          {
              type Newmark;
              gamma   0.1;    // Velocity integration coefficient
              beta    0.1;    // Position integration coefficient
          }
      
          accelerationRelaxation 1.0;
      
          //- prescribedRotation requires some sub-iterations to converge faster
          //  to desired value.
          nIter   3;
      
          bodies
          {
              rotor
              {
                  type            rigidBody;
                  parent          root;
                  centreOfMass    (-0.265 0 0);
                  mass            1;
                  rho             1;
                  inertia         (0.001 0 0 0.001 0 0.001);
                  transform       (1 0 0 0 1 0 0 0 1) (-0.265 0 0);
                  patches        (BLADE1);
                  innerDistance   1000;
                  outerDistance   2000;
                  joint
                  {
                      type            Rz;
                  }
              }
          }
      
          restraints
          {
              rotorRotation
              {
                  type                    prescribedRotation;
                  body                    rudder;
                  referenceOrientation    (1 0 0 0 1 0 0 0 1);
                  axis                    (0 0 1);    // Axis of rotation
                  omega                   sine;       // Function1 entry
                  omegaCoeffs
                  {
                      frequency   4;
                      amplitude   0.2;
                      scale       (0 1 0);
                      level       (0 1 0);
                  }
              }
             
          }
      

      但是如果把动网格文件改为:

      dynamicFvMesh       dynamicOversetFvMesh;
      
      dynamicOversetFvMeshCoeffs
      {
      //    layerRelax 0.3;
      }
      
      solver          multiSolidBodyMotionSolver;
      
      multiSolidBodyMotionSolverCoeffs
      {
          rotor
          {
              solidBodyMotionFunction rotatingMotion;
              rotatingMotionCoeffs
              {
                  origin      (0 0 0);
                  axis        (0 0 1);
                  omega       1.08;
              }
          }
      }
      

      就可以正常计算了。请问这是哪里出现错误了,跪求各位大神解答。

      posted in OpenFOAM
      L
      lv1995
    • RE: OpenFOAM重叠网格

      @wwzhao 你好,请问mergeMesh具体怎么用呢?需要添加那些信息或者有没有算例啊?

      posted in OpenFOAM
      L
      lv1995
    • openfoam中两套网格合并

      openfoam手册中写到mergeMeshes可以将两套网格合并,那么具体是怎么操作的呢?controlDict中需要添加什么信息?

      posted in OpenFOAM
      L
      lv1995
    • RE: bluecfd运算出错

      @dzw05
      k:

      dimensions      [0 2 -2 0 0 0 0];
      internalField   uniform 0.081;
      boundaryField
      { 
          inlet
          {
              type            zeroGradient;
          }
          outlet
          {
              type            zeroGradient;
          }
          frontAndBack
          {
              type            symmetry;
           }
          inter2
          {
              type            cyclicAMI;
           }
          inter1
          {
              type            cyclicAMI;
           }
      
          blades
          {
              type            kqRWallFunction;
              value           $internalField;
          }
      }
      

      nut

      dimensions      [0 2 -1 0 0 0 0];
      
      internalField   uniform 0;
      
      boundaryField
      {
              
          inlet
          {
              type            fixedValue;
              value           $internalField;
          }
      
          outlet
          {
              type             calculated;
              value           uniform 0;
      
          }
          frontAndBack
          {
              type            symmetry;
           }
          inter2
          {
              type            cyclicAMI;
           }
          inter1
          {
              type            cyclicAMI;
           }
           blades
          {
              type            nutkWallFunction;
              value           uniform 0;
          }
      }
      

      omega

      dimensions      [0 0 -1 0 0 0 0];
      
      internalField   uniform 0.095;
      
      boundaryField
      {
          
      
          outlet
          {
              type            zeroGradient;
          }
          blades
          {
              type            omegaWallFunction;
              value           $internalField;
          }
          inlet
          {
              type            zeroGradient;
          }
              frontAndBack
          {
              type            symmetry;
           }
          inter2
          {
              type            cyclicAMI;
           }
          inter1
          {
              type            cyclicAMI;
           }
      }
      

      p

      dimensions      [0 2 -2 0 0 0 0];
      
      internalField   uniform 0;
      
      boundaryField
      {
          inlet
          {
              type            zeroGradient;
          }
      
          outlet
          {
              type            fixedValue;
              value           uniform 0;
         }
              frontAndBack
          {
              type            symmetry;
           }
          inter2
          {
              type            cyclicAMI;
           }
          inter1
          {
              type            cyclicAMI;
           }
      
      
          blades
          {
              type            zeroGradient;
      
          }
      }
      

      u

      dimensions      [0 1 -1 0 0 0 0];
      
      internalField   uniform (14 0 0);
      
      boundaryField
      {
          
          inlet
          {
              type            fixedValue;
      
              value      uniform (14 0 0);
          }
      
          outlet
          {
              type            zeroGradient;
          }
              frontAndBack
          {
              type            symmetry;
           }
          inter2
          {
              type            cyclicAMI;
           }
          inter1
          {
              type            cyclicAMI;
           }
      
      
           blades
          {
              type            rotatingWallVelocity;
      
             origin           (0 0 0);
      
             axis             (0 0 1);
      
             omega           0;
          }
      }
      

      boundary

      inlet
          {
              type            patch;
              nFaces          400;
              startFace       6377766;
          }
          outlet
          {
              type            patch;
              nFaces          400;
              startFace       6378166;
          }
          frontAndBack
          {
              type            symmetry;
              inGroups        1(symmetry);
              nFaces          4800;
              startFace       6378566;
          }
          blades
          {
              type            wall;
              inGroups        1(wall);
              nFaces          80334;
              startFace       6383366;
          }
          inter1
          {
              type            cyclicAMI;
              inGroups        1(cyclicAMI);
              nFaces          13994;
              startFace       6463700;
              matchTolerance  0.0001;
              transform       noOrdering;
              neighbourPatch  inter2;
          }
          inter2
          {
              type            cyclicAMI;
              inGroups        1(cyclicAMI);
              nFaces          11572;
              startFace       6477694;
              matchTolerance  0.0001;
              transform       noOrdering;
              neighbourPatch  inter1;
          }
      

      fvSchemes

      ddtSchemes
      {
          default Euler;
      }
      
      gradSchemes
      {
          default         Gauss linear;
          grad(p)         Gauss linear;
          grad(U)         Gauss linear;
      }
      
      divSchemes
      {
          default         none;
          div(phi,U)      Gauss linearUpwind grad(U);
          div(phi,k)      Gauss limitedLinear 1;
          div(phi,omega)  Gauss limitedLinear 1;
          div((nuEff*dev2(T(grad(U))))) Gauss linear;
      }
      
      laplacianSchemes
      {
          default         Gauss linear limited corrected 0.5;
      }
      
      interpolationSchemes
      {
          default         linear;
      }
      
      snGradSchemes
      {
          default         corrected;
      }
      
      wallDist
      {
          method meshWave;
      }
      

      fvSolution

      solvers
      {
          "pcorr.*"
          {
              solver           GAMG;
              tolerance        0.02;
              relTol           0;
              smoother         GaussSeidel;
          }
      
          p
          {
              $pcorr;
              tolerance        1e-7;
              relTol           0.01;
          }
      
          pFinal
          {
              $p;
              tolerance        1e-7;
              relTol           0;
          }
      
          "(U|k|omega)"
          {
              solver          smoothSolver;
              smoother        symGaussSeidel;
              tolerance       1e-06;
              relTol          0.1;
          }
      
          "(U|k|omega)Final"
          {
              $U;
              tolerance       1e-06;
              relTol          0;
          }
      
          cellDisplacement
          {
              solver          GAMG;
              tolerance       1e-5;
              relTol          0;
              smoother        GaussSeidel;
          }
      }
      
      PIMPLE
      {
          correctPhi          yes;
          nOuterCorrectors    2;
          nCorrectors         2;
          nNonOrthogonalCorrectors 5;
          pRefCell  0;
          pRefPoint 0;
          pRefValue 0;
      }
      
      relaxationFactors
      {
          fields
          {
              p               0.2;
          }
          equations
          {
              "(U|k|omega)"   0.7;
              "(U|k|omega)Final" 0.7;
          }
      }
      
      cache
      {
          grad(U);
      }
      
      posted in OpenFOAM
      L
      lv1995
    • RE: bluecfd运算出错

      @东岳 请问还有可能是哪里的设置有问题呢?

      posted in OpenFOAM
      L
      lv1995
    • RE: bluecfd运算出错

      @东岳 还是一样的错误

      posted in OpenFOAM
      L
      lv1995
    • RE: bluecfd运算出错

      openfoam报错如下:

      #0  Foam::error::printStack(Foam::Ostream&) at ??:?
      #1  Foam::sigFpe::sigHandler(int) at ??:?
      #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
      #3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
      #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
      #5  ? at ??:?
      #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
      #7  ? at ??:?
      
      posted in OpenFOAM
      L
      lv1995
    • bluecfd运算出错

      使用pimpleDyMFoam求解器计算风机所受转矩,使用kωSST模型,运算第一步就出现如下错误:

      Create mesh for time = 0
      
      Selecting dynamicFvMesh dynamicMotionSolverFvMesh
      Selecting motion solver: solidBody
      Selecting solid-body motion function rotatingMotion
      Applying solid body motion to cellZone inter
      
      PIMPLE: no residual control data found. Calculations will employ 2 corrector loops
      
      Reading field p
      
      AMI: Creating addressing and weights between 14390 source faces and 12742 target faces
      AMI: Patch source sum(weights) min/max/average = 0.570448, 1.45816, 1.00126
      AMI: Patch target sum(weights) min/max/average = 0.000958814, 1.18605, 0.993037
      Reading field U
      
      Reading/calculating face flux field phi
      
      Selecting incompressible transport model Newtonian
      Selecting turbulence model type RAS
      Selecting RAS turbulence model kOmegaSST
      Selecting patchDistMethod meshWave
      RAS
      {
          RASModel        kOmegaSST;
          turbulence      on;
          printCoeffs     on;
          alphaK1         0.85;
          alphaK2         1;
          alphaOmega1     0.5;
          alphaOmega2     0.856;
          gamma1          0.555556;
          gamma2          0.44;
          beta1           0.075;
          beta2           0.0828;
          betaStar        0.09;
          a1              0.31;
          b1              1;
          c1              10;
          F3              false;
      }
      
      No MRF models present
      
      Reading/calculating face velocity Uf
      
      No finite volume options present
      
      Courant Number mean: 0.0006176 max: 0.225383
      
      Starting time loop
      
      forces Forces:
          Not including porosity effects
      forces forceCoeffs:
          Not including porosity effects
      forceCoeffs forceCoeffs:
          Not including porosity effects
      Courant Number mean: 0.0006176 max: 0.225383
      Time = 0.001
      
      AMI: Creating addressing and weights between 14390 source faces and 12742 target faces
      AMI: Patch source sum(weights) min/max/average = 0.570366, 1.45766, 1.00147
      AMI: Patch target sum(weights) min/max/average = 0, 1.18298, 0.993152
      GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 0.0220577, No Iterations 4
      GAMG:  Solving for pcorr, Initial residual = 0.230843, Final residual = 0.0159595, No Iterations 1
      GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
      GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
      GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
      GAMG:  Solving for pcorr, Initial residual = 0.0189167, Final residual = 0.0189167, No Iterations 0
      time step continuity errors : sum local = 2.90421e-009, global = -5.82674e-012, cumulative = -5.82674e-012
      PIMPLE: iteration 1
      smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.0228862, No Iterations 1
      smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.0203193, No Iterations 1
      smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.0226345, No Iterations 1
      Generating stack trace...
      
      
      Backtrace:
              ZN10StackTraceC1Ev [0x705c1465+0x25]
                       module: D:\blueCFD-Core-2017\ThirdParty-5.x\platforms\mingw_w64GccDPInt32\lib\libstack_trace.dll
              ZN4Foam5error10printStackERNS_7OstreamE [0x18d1c88+0x218]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
              ZN4Foam6sigFpe13sigFpeHandlerEi [0x18d2af3+0x33]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
              (No symbol) [0x4047cd]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              _C_specific_handler [0x7ffe9fcc7c58+0x98]
                       module: C:\WINDOWS\System32\msvcrt.dll
              0_chkstk [0x7ffe9ff8f7dd+0x11d]
                       module: C:\WINDOWS\SYSTEM32\ntdll.dll
              RtlWalkFrameChain [0x7ffe9fefd856+0x13f6]
                       module: C:\WINDOWS\SYSTEM32\ntdll.dll
              KiUserExceptionDispatcher [0x7ffe9ff8e70e+0x2e]
                       module: C:\WINDOWS\SYSTEM32\ntdll.dll
              ZN4Foam6divideERNS_5FieldIdEERKdRKNS_5UListIdEE [0x1877fa4+0x24]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\lib\libOpenFOAM.dll
              (No symbol) [0x42c1e7]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              (No symbol) [0x4361db]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              (No symbol) [0x4363dd]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              (No symbol) [0x4546ca]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              (No symbol) [0x4013f7]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              (No symbol) [0x40152b]
                       module: D:\blueCFD-Core-2017\OpenFOAM-5.x\platforms\mingw_w64GccDPInt32Opt\bin\pimpleDyMFoam.exe
              BaseThreadInitThunk [0x7ffe9d393dc4+0x14]
                       module: C:\WINDOWS\System32\KERNEL32.DLL
              RtlUserThreadStart [0x7ffe9ff63691+0x21]
                       module: C:\WINDOWS\SYSTEM32\ntdll.dll
      

      求教各位大神,这是哪里的设置出现问题?

      posted in OpenFOAM
      L
      lv1995
    • RE: 关于snappy与cfmesh

      @wayne 你好,我现在也开始做船,但是因为刚起步,然后船舶的网格比较复杂,我不太知道应该怎么划分,可以借鉴您的网格参考一下么

      posted in OpenFOAM
      L
      lv1995
    • RE: 如何输出某点随时间变化的坐标?

      @东岳 这个模型是根据实际的物体建的,所以模型重心就直接取了实际物体的重心坐标,没有数学公式。

      posted in OpenFOAM
      L
      lv1995