Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. overset 小间隙-多浮体算例 计算结果明显有误

overset 小间隙-多浮体算例 计算结果明显有误

已定时 已固定 已锁定 已移动 OpenFOAM
10 帖子 5 发布者 4.7k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • H 离线
    H 离线
    hachikoi
    写于 最后由 编辑
    #1

    最近在做同轴双浮体装置,目前是内浮体不动(后续会动,所以做成了overset),外浮体初始位移为-1做自衰减

    内浮体的overset较大,将外浮体overset包裹,四面消波处理

    7b20a5fc-282a-4b85-a71c-7137b1a657e7-V@Y4DVKX`BUHWPMRE725GU5.png

    目前计算结果,外浮体在最初的几秒存在负速度,这是完全不合理的,之前尝试过单独外浮体/单独内浮体/外浮体overset+内浮体背景扣除都是没有问题的。怀疑是overset插值引起的错误结果,请大家帮忙看看是什么问题引起的呢?:chouchou:

    fvSchemes:

    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        object      fvSchemes;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    ddtSchemes
    {
        default                             Euler;
    }
    
    gradSchemes
    {
        default                             Gauss linear;
    }
    
    divSchemes
    {                    
        div(rhoPhi,U)                       Gauss limitedLinearV 1;
        div(U)                              Gauss linear;
        div(phi,alpha)                      Gauss vanLeer;
        div(phirb,alpha)                    Gauss linear;
        div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
        div(phi,k)                          Gauss upwind;
        div(phi,epsilon)                    Gauss upwind;
        div(phi,omega)                      Gauss upwind;
    }
    
    laplacianSchemes
    {
        default                             Gauss linear corrected;
    }
    
    interpolationSchemes
    {
        default                             linear;
    }
    
    snGradSchemes
    {
        default                             corrected;
    }
    
    oversetInterpolation
    {
        method                              inverseDistance;
    }
    
    fluxRequired
    {
        default                             no;
        p_rgh;
        pcorr;
        alpha.water;
    }
    
    wallDist
    {
        method                              meshWave;
        //method                              Poisson;
    }
    
    oversetInterpolationSuppressed
    {
        grad(p_rgh);
        surfaceIntegrate(phiHbyA);
    }
    
    /*
    oversetInterpolationRequired
    {
        k;
        epsilon;
        omega;
        nut;
        alpha.water;
    }
    */
    
    
    // ************************************************************************* //
    

    fvSolution:

    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        object      fvSolution;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    solvers
    {
        "cellDisplacement.*"
        {
            solver               PCG;
            preconditioner       DIC;
            tolerance            1e-06;
            relTol               0;
            maxIter              100;
        }
    
        "alpha.water.*"
        {
            nAlphaCorr           2;
            nAlphaSubCycles      1;
            cAlpha               1;
            icAlpha              0;
            MULESCorr            yes;
            nLimiterIter         5;
            alphaApplyPrevCorr   no;
            solver               smoothSolver;
            smoother             symGaussSeidel;
            tolerance            1e-08;
            relTol               0;
        }
    
        "pcorr.*"
        {
            solver               PCG;
            preconditioner       DIC;
            tolerance            1e-09;
            relTol               0;
        }
    
        p_rgh
        {
            solver               PBiCGStab;
            preconditioner       DILU;
            tolerance            1e-09;
            relTol               0.01;
        }
    
        p_rghFinal
        {
            $p_rgh;
            relTol               0;
        }
    
        "(U|k|omega|epsilon).*"
        {
            solver               smoothSolver;
            smoother             symGaussSeidel;
            tolerance            1e-08;
            relTol               0;
        }
    
        yPsi
        {
            solver               PBiCGStab;
            preconditioner       DILU;
            tolerance            1e-6;
            relTol               0.0;
        }
    
        yPsiFinal
        {
            $yPsi;
            relTol               0;
        }
    }
    
    PIMPLE
    {
        momentumPredictor        no;
        nOuterCorrectors         5;
        nCorrectors              2;
        nNonOrthogonalCorrectors 1;
        ddtCorr                  yes;
        correctPhi               no;
        moveMeshOuterCorrectors  no;
        turbOnFinalIterOnly      no;
        oversetAdjustPhi         yes;
    }
    
    relaxationFactors
    {
        fields
        {
        }
        equations
        {
            ".*"                 1;
        }
    }
    
    cache
    {
       grad(U); 
    }
    
    
    // ************************************************************************* //
    
    

    dynamicMeshDict

    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       dictionary;
        object      dynamicMeshDict;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dynamicFvMesh                    dynamicOversetFvMesh;
    
    solvers
    {
        dynamicPlatform
        {
            motionSolverLibs         (libfvMotionSolvers);
            motionSolver             multiSolidBodyMotionSolver;
            multiSolidBodyMotionSolverCoeffs
            {
                zonePlatform
                {
                    solidBodyMotionFunction  linearMotion;
                    linearMotionCoeffs
                    {
                        velocity             (0 0 0);
                    }
                }
            }
        }
    
        dynamicWECMiddle
        {
            motionSolverLibs         (librigidBodyMeshMotion);
            motionSolver             rigidBodyMotion;
            report                   on;
            cellSet                  cellWECMiddle;
            cOfGdisplacement         CofG;
            bodyIdCofG               1;
            solver
            {
                type                 Newmark;
            }
            accelerationRelaxation   0.8;
            accelerationDamping      0.9;
            nIter                    3;
            bodies
            {
                dynamicWECMiddle
                {
                    type             rigidBody;
                    parent           root;
                    mass             2.0126e+06;
                    centreOfMass     (0 0 0); // COG in the local coordinate system, based on transform
                    inertia          (1.3312e+08 0 0 1.3312e+08 0 2.1384e+08); // (Ixx Ixy Ixz Iyy Iyz Izz)
                    transform        (1 0 0 0 1 0 0 0 1) (0 0 -1);
                    joint
                    {
                        type         Pz; 
                    }
                    patches
                    (
                        WECMiddle
                    );
                    innerDistance    100;
                    outerDistance    101;
                }
            }
            restraints
            {
            }
    
            /*
            motionSolverLibs   (sixDoFRigidBodyMotion);
            motionSolver       sixDoFRigidBodyMotion;
            cellSet            cellWECMiddle;
            patches            (WECMiddle);
            innerDistance      100;
            outerDistance      101;
            centreOfMass       (0 0 -1);
            cOfGdisplacement   CofG;
            mass               2.0126e+06;
            momentOfInertia    (1.3312e+08 1.3312e+08 2.1384e+08);
            report             on;
            accelerationRelaxation 0.8;
            accelerationDamping    0.9;
            solver
            {
                type           Newmark;
            }
            constraints
            {
                zAxis
                {
                    sixDoFRigidBodyMotionConstraint line;
                    direction  (0 0 1);
                }
            }
            */
        }
    
    }
    
    
    // ************************************************************************* //
    

    0:

    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volScalarField;
        object      alpha.water;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 0 0 0 0 0 0];
    
    internalField   uniform 0;
    
    boundaryField
    {
        #includeEtc "caseDicts/setConstraintTypes"
    
        inlet
        {
            type            zeroGradient;
        }
        
        outlet
        {
            type            zeroGradient;
        }
    
        "(rightSide|leftSide)"
        {
            type            zeroGradient;
        }    
    
        lowerSide
        {
            type            zeroGradient;
        }  
    
        atmosphere
        {
            type            inletOutlet;
            inletValue      $internalField;
            value           $internalField;
        }
    
        oversetWECMiddle
        {
            type            overset;
        }
    
        WECMiddle
        {
            type            zeroGradient;
        }
    
        oversetPlatform
        {
            type            overset;
        }
    
        platform
        {
            type            zeroGradient;
        }
    }
    
    
    // ************************************************************************* //
    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volScalarField;
        object      p_rgh;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [1 -1 -2 0 0 0 0];
    
    internalField   uniform 0;
    
    boundaryField
    {
        #includeEtc "caseDicts/setConstraintTypes"
    
        inlet
        {
            type            fixedFluxPressure;
            value           $internalField;
        }
        
        outlet
        {
            type            fixedFluxPressure;
            value           $internalField;
        }
    
        "(rightSide|leftSide)"
        {
            type            fixedFluxPressure;
            value           $internalField;
        }  
    
        lowerSide
        {
            type            fixedFluxPressure;
            value           $internalField;
        }  
    
        atmosphere
        {
            type            totalPressure;
            p0              $internalField;
        }
    
        oversetWECMiddle
        {
            type            overset;
        }
    
        WECMiddle
        {
            type            fixedFluxPressure;
        }
    
        oversetPlatform
        {
            type            overset;
        }
    
        platform
        {
            type            fixedFluxPressure;
        }
    }
    
    
    // ************************************************************************* //
    
    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       pointVectorField;
        object      pointDisplacement;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 1 0 0 0 0 0];
    
    internalField   uniform (0 0 0);
    
    boundaryField
    {
        #includeEtc "caseDicts/setConstraintTypes"
    
        inlet
        {
            type            fixedValue;
            value           $internalField;
        }
        
        outlet
        {
            type            fixedValue;
            value           $internalField;
        }
    
        "(rightSide|leftSide)"
        {
            type            fixedValue;
            value           $internalField;
        }  
      
        lowerSide
        {
            type            fixedValue;
            value           $internalField;
        }  
    
        atmosphere
        {
            type            fixedValue;
            value           $internalField;
        }
    
        oversetWECMiddle
        {
            patchType       overset;
            type            zeroGradient;
        }
    
        WECMiddle
        {
            type            calculated;
            value           $internalField;
        }
    
        oversetPlatform
        {
            patchType       overset;
            type            zeroGradient;
        }
    
        platform
        {
            type            fixedValue;
            value           $internalField;
        }
    }
    
    
    // ************************************************************************* //
    
    
    /*--------------------------------*- C++ -*----------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  v2306                                 |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    FoamFile
    {
        version     2.0;
        format      ascii;
        class       volVectorField;
        object      U;
    }
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    
    dimensions      [0 1 -1 0 0 0 0];
    
    internalField   uniform (0 0 0);
    
    boundaryField
    {
        #includeEtc "caseDicts/setConstraintTypes"
    
        inlet
        {
            type            waveVelocity;
            value           $internalField;
        }
        
        outlet
        {
            type            waveVelocity;
            value           $internalField;
        }
    
        "(rightSide|leftSide)"
        {
            type            waveVelocity;
            value           $internalField;
        }  
        
        lowerSide
        {
            type            fixedValue;
            value           $internalField;
        }  
    
        atmosphere
        {
            type            pressureInletOutletVelocity;
            value           $internalField;
        }
    
        oversetWECMiddle
        {
            type            overset;
        }
    
        WECMiddle
        {
            type            movingWallVelocity;
            value           $internalField;
        }
    
        oversetPlatform
        {
            type            overset;
        }
    
        platform
        {
            type            fixedValue;
            value           $internalField;
        }
    }
    
    
    // ************************************************************************* //
    
    

    请大家帮忙看看是不是哪里设置有问题

    1 条回复 最后回复
  • H 离线
    H 离线
    hachikoi
    写于 最后由 编辑
    #2

    计算结果

    /*---------------------------------------------------------------------------*\
    | =========                 |                                                 |
    | \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
    |  \\    /   O peration     | Version:  2306                                  |
    |   \\  /    A nd           | Website:  www.openfoam.com                      |
    |    \\/     M anipulation  |                                                 |
    \*---------------------------------------------------------------------------*/
    Build  : _fbf00d6bf2-20230626 OPENFOAM=2306 version=v2306
    Arch   : "LSB;label=32;scalar=64"
    Exec   : overInterDyMFoam -parallel
    Date   : Nov 24 2023
    Time   : 21:42:20
    Host   : DESKTOP-KVQEVMV
    PID    : 617
    I/O    : uncollated
    Case   : /mnt/d/Tutorials/UMaine/decay/decayRASkOmegaSST/decayWECMiddleSpar/decayHeave1V0/background
    nProcs : 8
    Hosts  :
    (
        (DESKTOP-KVQEVMV 8)
    )
    Pstream initialized with:
        floatTransfer      : false
        nProcsSimpleSum    : 0
        nonBlockingExchange: 0 (tuning: 0)
        exchange algorithm : -1
        commsType          : nonBlocking
        polling iterations : 0
    trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
    fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
    allowSystemOperations : Allowing user-supplied system call operations
    
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Create time
    
    Create mesh for time = 0
    
    Selecting dynamicFvMesh dynamicOversetFvMesh
    Selecting motion solver: multiSolidBodyMotionSolver
    Applying motion to entire mesh
    Selecting solid-body motion function linearMotion
    Applying solid body motion linearMotion to 1449278 points of cellZone zonePlatform
    Selecting motion solver: rigidBodyMotion
    Applying motion to cellSet: cellWECMiddle
    Selecting rigidBodySolver Newmark
    
    PIMPLE: no residual control data found. Calculations will employ 5 corrector loops
    
    Reading field p_rgh
    
    Reading field U
    
    Reading/calculating face flux field phi
    
    Creating cellMask field to block out hole cells
    
    --> FOAM Warning : 
        From bool Foam::oversetPolyPatch::master() const
        in file oversetPolyPatch/oversetPolyPatch.C at line 145
        The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation.
    --> FOAM Warning : 
        From bool Foam::oversetPolyPatch::master() const
        in file oversetPolyPatch/oversetPolyPatch.C at line 145
        The master overset patch is not the first patch. Generally the first patch should be an overset patch to guarantee consistent operation.
    Creating interpolatedCells field 
    
    Reading transportProperties
    
    Selecting incompressible transport model Newtonian
    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 g
    
    Reading hRef
    Calculating field g.h
    
    No MRF models present
    
    No finite volume options present
    Reading/calculating face velocity Uf
    
    Courant Number mean: 0 max: 0
    
    Starting time loop
    
    Courant Number mean: 0 max: 0
    Interface Courant Number mean: 0 max: 0
    deltaT = 0.012
    Time = 0.012
    
    PIMPLE: iteration 1
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999950388828)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 0.00826852870788)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503453812) (11.2436789924 17.5000299924 10.2504506036)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    Selecting waveModel shallowWaterAbsorption
    Wave model: patch inlet
        Type : shallowWaterAbsorption
        Velocity field name : U
        Phase fraction field name : alpha.water
        Transformation from local to global system : (1 -7.77501308186e-32 7.17951403026e-23 7.77501308186e-32 1 0 -0 -0 1)
        Number of paddles: 1
        Reference water depth : 199.999999317
        Active absorption: 1
    Updating shallowWaterAbsorption wave model for patch inlet
    Selecting waveModel shallowWaterAbsorption
    Wave model: patch outlet
        Type : shallowWaterAbsorption
        Velocity field name : U
        Phase fraction field name : alpha.water
        Transformation from local to global system : (-1 -3.16385813721e-32 7.17951403181e-23 3.16385813721e-32 -1 0 -0 -0 1)
        Number of paddles: 1
        Reference water depth : 199.999999317
        Active absorption: 1
    Updating shallowWaterAbsorption wave model for patch outlet
    Selecting waveModel shallowWaterAbsorption
    Wave model: patch rightSide
        Type : shallowWaterAbsorption
        Velocity field name : U
        Phase fraction field name : alpha.water
        Transformation from local to global system : (-1.61558713389e-32 1 1.05282388693e-18 -1 -1.61558713389e-32 -0 -0 -0 1)
        Number of paddles: 1
        Reference water depth : 199.999999317
        Active absorption: 1
    Updating shallowWaterAbsorption wave model for patch rightSide
    Selecting waveModel shallowWaterAbsorption
    Wave model: patch leftSide
        Type : shallowWaterAbsorption
        Velocity field name : U
        Phase fraction field name : alpha.water
        Transformation from local to global system : (2.15411617852e-32 -1 -5.88975123258e-18 1 2.15411617852e-32 0 -0 -0 1)
        Number of paddles: 1
        Reference water depth : 199.999999317
        Active absorption: 1
    Updating shallowWaterAbsorption wave model for patch leftSide
    Execution time for mesh.update() = 31.64 s
    smoothSolver:  Solving for alpha.water, Initial residual = 2.61985281096e-07, Final residual = 9.3161532224e-09, No Iterations 23
    Phase-1 volume fraction = 0.66655233285  Min(alpha.water) = 0  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.66655233262  Min(alpha.water) = -4.06575814682e-23  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1, Final residual = 0.00934508785126, No Iterations 132
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.93805455326e-05, Final residual = 9.81032812434e-08, No Iterations 157
    time step continuity errors : sum local = 1.31199443232e-09, global = 5.31630376358e-10, cumulative = 5.31630376358e-10
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.05167563797e-05, Final residual = 1.04029750059e-07, No Iterations 144
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.3336587163e-06, Final residual = 6.33395425015e-10, No Iterations 233
    time step continuity errors : sum local = 6.1751513569e-12, global = 5.28990371798e-12, cumulative = 5.36920280076e-10
    smoothSolver:  Solving for omega, Initial residual = 0.000102569828088, Final residual = 2.86343248998e-09, No Iterations 3
    smoothSolver:  Solving for k, Initial residual = 0.999999998453, Final residual = 9.33528038491e-09, No Iterations 66
    PIMPLE: iteration 2
    smoothSolver:  Solving for alpha.water, Initial residual = 1.27402311839e-07, Final residual = 8.81590646774e-09, No Iterations 2
    Phase-1 volume fraction = 0.666552379444  Min(alpha.water) = -3.83779065969e-28  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552381689  Min(alpha.water) = -1.621967531e-06  Max(alpha.water) = 1.00000000013
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.61393401599e-05, Final residual = 2.95529990406e-07, No Iterations 137
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.32733975555e-05, Final residual = 8.98829957804e-07, No Iterations 58
    time step continuity errors : sum local = 2.97727638562e-09, global = 1.12545136858e-09, cumulative = 1.66237164866e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.25486710901e-05, Final residual = 7.9528642603e-07, No Iterations 111
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.1598896674e-05, Final residual = 2.69740404024e-10, No Iterations 215
    time step continuity errors : sum local = 8.87609898546e-13, global = 2.19947839059e-13, cumulative = 1.6625915965e-09
    smoothSolver:  Solving for omega, Initial residual = 1.05706805697e-06, Final residual = 3.06165398729e-09, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 2.17813199337e-05, Final residual = 8.83104952045e-09, No Iterations 11
    PIMPLE: iteration 3
    smoothSolver:  Solving for alpha.water, Initial residual = 1.08796870261e-07, Final residual = 6.88901360938e-09, No Iterations 3
    Phase-1 volume fraction = 0.666552350255  Min(alpha.water) = -9.78968376404e-18  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552348822  Min(alpha.water) = -4.03306452835e-06  Max(alpha.water) = 1.00000000001
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000260128355201, Final residual = 2.38866841716e-06, No Iterations 80
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.97760690631e-05, Final residual = 3.74084676736e-07, No Iterations 98
    time step continuity errors : sum local = 6.45059786074e-10, global = 6.27376936037e-10, cumulative = 2.28996853253e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.18893127841e-05, Final residual = 4.93880984564e-07, No Iterations 140
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.00033869064e-05, Final residual = 9.88435829622e-10, No Iterations 215
    time step continuity errors : sum local = 1.79260009404e-12, global = 8.66436780833e-14, cumulative = 2.29005517621e-09
    smoothSolver:  Solving for omega, Initial residual = 2.53541232093e-07, Final residual = 3.5196014854e-09, No Iterations 2
    smoothSolver:  Solving for k, Initial residual = 2.94650853902e-05, Final residual = 7.53900830036e-09, No Iterations 11
    PIMPLE: iteration 4
    smoothSolver:  Solving for alpha.water, Initial residual = 7.59424341536e-08, Final residual = 8.47687199918e-09, No Iterations 1
    Phase-1 volume fraction = 0.666552350369  Min(alpha.water) = -1.7236282171e-09  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552350436  Min(alpha.water) = -0.000416460503671  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000167246351824, Final residual = 1.40182075997e-06, No Iterations 107
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.16507599585e-05, Final residual = 4.37468479875e-07, No Iterations 64
    time step continuity errors : sum local = 6.00063571489e-10, global = -2.63397570915e-10, cumulative = 2.0266576053e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.29697246248e-05, Final residual = 8.11491883613e-07, No Iterations 151
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.44922341191e-05, Final residual = 6.85329760483e-10, No Iterations 217
    time step continuity errors : sum local = 1.22184650642e-12, global = -9.18162980737e-13, cumulative = 2.02573944232e-09
    smoothSolver:  Solving for omega, Initial residual = 1.52703236419e-07, Final residual = 7.73512960224e-09, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 8.02738033721e-06, Final residual = 9.37089259035e-09, No Iterations 23
    PIMPLE: iteration 5
    smoothSolver:  Solving for alpha.water, Initial residual = 5.62034902794e-08, Final residual = 8.1720263608e-09, No Iterations 1
    Phase-1 volume fraction = 0.666552346355  Min(alpha.water) = -3.66106998989e-10  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552344225  Min(alpha.water) = -0.000289714267251  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000179222117898, Final residual = 1.69745284627e-06, No Iterations 123
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.19874219811e-05, Final residual = 1.89807848936e-07, No Iterations 148
    time step continuity errors : sum local = 3.71990181733e-10, global = 2.12074694904e-11, cumulative = 2.04694691181e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.41850750136e-05, Final residual = 3.01352648671e-07, No Iterations 126
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.45426477763e-05, Final residual = 9.45301356671e-10, No Iterations 178
    time step continuity errors : sum local = 1.78223171988e-12, global = -2.78459422187e-13, cumulative = 2.04666845238e-09
    smoothSolver:  Solving for omega, Initial residual = 1.60789297939e-07, Final residual = 8.69888544542e-09, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 4.85804248814e-06, Final residual = 9.50862177947e-09, No Iterations 7
    ExecutionTime = 557.59 s  ClockTime = 560 s
    
    Courant Number mean: 1.63749449277e-05 max: 49.9913602954
    Interface Courant Number mean: 6.95403940679e-07 max: 3.55036474521
    deltaT = 0.000240041477749
    Time = 0.01224
    
    PIMPLE: iteration 1
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999948351772)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 0.00870400248066)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503433442) (11.2436789924 17.5000299924 10.2504526406)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    Updating shallowWaterAbsorption wave model for patch inlet
    Updating shallowWaterAbsorption wave model for patch outlet
    Updating shallowWaterAbsorption wave model for patch rightSide
    Updating shallowWaterAbsorption wave model for patch leftSide
    Execution time for mesh.update() = 40.12 s
    smoothSolver:  Solving for alpha.water, Initial residual = 8.36900182683e-09, Final residual = 8.36900182683e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552343446  Min(alpha.water) = -0.000289714267251  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552342656  Min(alpha.water) = -0.000289584326027  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.0228028713103, Final residual = 0.000158123063473, No Iterations 103
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000476023249144, Final residual = 4.19682603135e-06, No Iterations 132
    time step continuity errors : sum local = 8.93614961201e-11, global = 7.31703459824e-11, cumulative = 2.11983879837e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000192656249381, Final residual = 1.36166337649e-06, No Iterations 163
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 6.70954837961e-05, Final residual = 9.678157846e-10, No Iterations 265
    time step continuity errors : sum local = 1.58991434156e-14, global = -7.99361856081e-15, cumulative = 2.11983080475e-09
    smoothSolver:  Solving for omega, Initial residual = 1.65654710594e-06, Final residual = 9.11101653564e-10, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 9.84944919853e-05, Final residual = 7.48999418113e-09, No Iterations 2
    PIMPLE: iteration 2
    smoothSolver:  Solving for alpha.water, Initial residual = 6.06658567684e-09, Final residual = 6.06658567684e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552342656  Min(alpha.water) = -0.000289584326027  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552341998  Min(alpha.water) = -0.000289373569159  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 6.09071653377e-05, Final residual = 5.86498406611e-07, No Iterations 135
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.67956672479e-05, Final residual = 2.01508297651e-07, No Iterations 149
    time step continuity errors : sum local = 3.3067714853e-12, global = 3.30648753556e-12, cumulative = 2.12313729228e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.18188041175e-05, Final residual = 7.73661385616e-08, No Iterations 156
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.55315884803e-06, Final residual = 6.90107537102e-10, No Iterations 213
    time step continuity errors : sum local = 1.07514813361e-14, global = 1.01727697345e-14, cumulative = 2.12314746505e-09
    smoothSolver:  Solving for omega, Initial residual = 1.82307346832e-09, Final residual = 1.82307346832e-09, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 6.07724465011e-08, Final residual = 3.64698977772e-09, No Iterations 1
    PIMPLE: iteration 3
    smoothSolver:  Solving for alpha.water, Initial residual = 4.95615350895e-09, Final residual = 4.95615350895e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552341998  Min(alpha.water) = -0.000289373569159  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552341433  Min(alpha.water) = -0.00028921797597  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.34942515551e-06, Final residual = 2.89216757541e-08, No Iterations 153
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.71613055767e-06, Final residual = 1.69845429434e-08, No Iterations 139
    time step continuity errors : sum local = 2.85144328781e-13, global = 2.82486392054e-13, cumulative = 2.12342995145e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 7.57358817177e-07, Final residual = 4.68743757682e-09, No Iterations 148
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.30816263038e-07, Final residual = 7.98377737781e-10, No Iterations 155
    time step continuity errors : sum local = 1.24636818568e-14, global = 1.09792420598e-14, cumulative = 2.12344093069e-09
    smoothSolver:  Solving for omega, Initial residual = 1.92226875551e-09, Final residual = 1.92226875551e-09, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 1.34277181319e-08, Final residual = 2.5244507753e-09, No Iterations 1
    PIMPLE: iteration 4
    smoothSolver:  Solving for alpha.water, Initial residual = 4.14180619703e-09, Final residual = 4.14180619703e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552341433  Min(alpha.water) = -0.00028921797597  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.66655234095  Min(alpha.water) = -0.000289070037242  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.69585590378e-07, Final residual = 8.74528894708e-09, No Iterations 140
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.81105228656e-07, Final residual = 2.72909179203e-09, No Iterations 126
    time step continuity errors : sum local = 4.56821656802e-14, global = 1.54082829784e-14, cumulative = 2.12345633897e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.25226124652e-08, Final residual = 6.67713245408e-10, No Iterations 124
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.2978743237e-08, Final residual = 9.74587014884e-10, No Iterations 109
    time step continuity errors : sum local = 1.66374816607e-14, global = 1.64971601989e-14, cumulative = 2.12347283613e-09
    smoothSolver:  Solving for omega, Initial residual = 1.95628819863e-09, Final residual = 1.95628819863e-09, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 5.62125672229e-09, Final residual = 5.62125672229e-09, No Iterations 0
    PIMPLE: iteration 5
    smoothSolver:  Solving for alpha.water, Initial residual = 3.60189757848e-09, Final residual = 3.60189757848e-09, No Iterations 0
    Phase-1 volume fraction = 0.66655234095  Min(alpha.water) = -0.000289070037242  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552340533  Min(alpha.water) = -0.000288909525708  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.67847330072e-07, Final residual = 1.98395076727e-09, No Iterations 149
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.06975221441e-07, Final residual = 1.02026137309e-09, No Iterations 136
    time step continuity errors : sum local = 1.74929356065e-14, global = 1.46083619252e-14, cumulative = 2.12348744449e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.31574420273e-08, Final residual = 9.65472946334e-10, No Iterations 119
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.49455168275e-08, Final residual = 8.73609993614e-10, No Iterations 102
    time step continuity errors : sum local = 1.40323018836e-14, global = 2.66664498002e-15, cumulative = 2.12349011114e-09
    smoothSolver:  Solving for omega, Initial residual = 1.95616404291e-09, Final residual = 1.95616404291e-09, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 5.67015345545e-09, Final residual = 5.67015345545e-09, No Iterations 0
    ExecutionTime = 1116.64 s  ClockTime = 1124 s
    
    Courant Number mean: 1.23753094205e-06 max: 2.67476268045
    Interface Courant Number mean: 1.22084714043e-08 max: 0.040419192207
    deltaT = 8.97430936596e-05
    Time = 0.0123298
    
    PIMPLE: iteration 1
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999947523286)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 0.00975950626086)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503425157) (11.2436789924 17.5000299924 10.2504534691)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    Updating shallowWaterAbsorption wave model for patch inlet
    Updating shallowWaterAbsorption wave model for patch outlet
    Updating shallowWaterAbsorption wave model for patch rightSide
    Updating shallowWaterAbsorption wave model for patch leftSide
    Execution time for mesh.update() = 39.64 s
    smoothSolver:  Solving for alpha.water, Initial residual = 1.89506393044e-09, Final residual = 1.89506393044e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552340293  Min(alpha.water) = -0.000288909525708  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552339977  Min(alpha.water) = -0.000288846889594  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.00695412581841, Final residual = 6.58530505703e-05, No Iterations 148
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000208251042114, Final residual = 2.05743333992e-06, No Iterations 125
    time step continuity errors : sum local = 5.895261321e-11, global = -3.50991476391e-11, cumulative = 2.0883909635e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.00010380483611, Final residual = 9.76125551947e-07, No Iterations 112
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.2687177802e-05, Final residual = 9.50115732418e-10, No Iterations 213
    time step continuity errors : sum local = 2.05769533909e-14, global = 1.23631828244e-14, cumulative = 2.08840332668e-09
    smoothSolver:  Solving for omega, Initial residual = 6.1051758878e-07, Final residual = 4.82103625934e-10, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 3.01714859593e-05, Final residual = 9.37540222535e-09, No Iterations 1
    PIMPLE: iteration 2
    smoothSolver:  Solving for alpha.water, Initial residual = 1.63714691561e-09, Final residual = 1.63714691561e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552339977  Min(alpha.water) = -0.000288846889594  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552339706  Min(alpha.water) = -0.000288806485034  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.23761401106e-05, Final residual = 2.08661810278e-07, No Iterations 134
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.2229647521e-06, Final residual = 8.84903438452e-08, No Iterations 146
    time step continuity errors : sum local = 1.93794820853e-12, global = 1.93793640713e-12, cumulative = 2.09034126309e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.31258175956e-06, Final residual = 3.09958221974e-08, No Iterations 129
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.16173305739e-06, Final residual = 7.61858990072e-10, No Iterations 182
    time step continuity errors : sum local = 1.79845063561e-14, global = -7.97949298118e-15, cumulative = 2.0903332836e-09
    smoothSolver:  Solving for omega, Initial residual = 6.86324608392e-10, Final residual = 6.86324608392e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 2.59199319034e-08, Final residual = 3.72155986805e-09, No Iterations 1
    PIMPLE: iteration 3
    smoothSolver:  Solving for alpha.water, Initial residual = 1.464298727e-09, Final residual = 1.464298727e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552339706  Min(alpha.water) = -0.000288806485034  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552339468  Min(alpha.water) = -0.000288763304213  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.11368632855e-06, Final residual = 9.12572869945e-09, No Iterations 145
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.51282511339e-07, Final residual = 4.39016946545e-09, No Iterations 148
    time step continuity errors : sum local = 9.74366513151e-14, global = -9.73777380055e-14, cumulative = 2.09023590586e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.01376950811e-07, Final residual = 9.97729577993e-10, No Iterations 132
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.03614929149e-08, Final residual = 6.77515000163e-10, No Iterations 146
    time step continuity errors : sum local = 1.65113572261e-14, global = -1.60186025025e-14, cumulative = 2.09021988725e-09
    smoothSolver:  Solving for omega, Initial residual = 6.93072732889e-10, Final residual = 6.93072732889e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 4.6403143789e-09, Final residual = 4.6403143789e-09, No Iterations 0
    PIMPLE: iteration 4
    smoothSolver:  Solving for alpha.water, Initial residual = 1.34365518119e-09, Final residual = 1.34365518119e-09, No Iterations 0
    Phase-1 volume fraction = 0.666552339468  Min(alpha.water) = -0.000288763304213  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.66655233926  Min(alpha.water) = -0.000288719352487  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 7.66014318067e-08, Final residual = 8.71802548845e-10, No Iterations 148
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.2919023473e-08, Final residual = 9.70813176141e-10, No Iterations 97
    time step continuity errors : sum local = 2.14135270487e-14, global = -1.99862439702e-14, cumulative = 2.09019990101e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.08856570606e-08, Final residual = 9.76022773547e-10, No Iterations 58
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.07307604617e-09, Final residual = 9.29257136903e-10, No Iterations 21
    time step continuity errors : sum local = 1.67526610845e-14, global = 8.81445888723e-15, cumulative = 2.09020871547e-09
    smoothSolver:  Solving for omega, Initial residual = 6.93390357342e-10, Final residual = 6.93390357342e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 4.67499439824e-09, Final residual = 4.67499439824e-09, No Iterations 0
    PIMPLE: iteration 5
    smoothSolver:  Solving for alpha.water, Initial residual = 1.2520957575e-09, Final residual = 1.2520957575e-09, No Iterations 0
    Phase-1 volume fraction = 0.66655233926  Min(alpha.water) = -0.000288719352487  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552339079  Min(alpha.water) = -0.000288674615659  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.50640815229e-08, Final residual = 9.99499573606e-10, No Iterations 31
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.01438500186e-09, Final residual = 9.86510742554e-10, No Iterations 10
    time step continuity errors : sum local = 1.8597559213e-14, global = 9.98907756168e-15, cumulative = 2.09021870455e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.37064254853e-09, Final residual = 9.87017098349e-10, No Iterations 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.9024984935e-10, Final residual = 9.9024984935e-10, No Iterations 0
    time step continuity errors : sum local = 1.88156236059e-14, global = 9.85210599887e-15, cumulative = 2.09022855665e-09
    smoothSolver:  Solving for omega, Initial residual = 6.93401046968e-10, Final residual = 6.93401046968e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 4.67636350062e-09, Final residual = 4.67636350062e-09, No Iterations 0
    ExecutionTime = 1551.51 s  ClockTime = 1567 s
    
    Courant Number mean: 7.44659378168e-07 max: 0.33617200472
    Interface Courant Number mean: 5.92610839418e-09 max: 0.0133063005972
    deltaT = 0.000107691712392
    Time = 0.0124375
    
    PIMPLE: iteration 1
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999946925817)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 0.00133640649968)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503419182) (11.2436789924 17.5000299924 10.2504540666)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    
    
    1 条回复 最后回复
  • H 离线
    H 离线
    hachikoi
    写于 最后由 编辑
    #3
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999946925817)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 0.00133640649968)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503419182) (11.2436789924 17.5000299924 10.2504540666)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    Updating shallowWaterAbsorption wave model for patch inlet
    Updating shallowWaterAbsorption wave model for patch outlet
    Updating shallowWaterAbsorption wave model for patch rightSide
    Updating shallowWaterAbsorption wave model for patch leftSide
    Execution time for mesh.update() = 41.86 s
    smoothSolver:  Solving for alpha.water, Initial residual = 8.27297434572e-10, Final residual = 8.27297434572e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338974  Min(alpha.water) = -0.000288674615659  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338827  Min(alpha.water) = -0.000288619975557  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000659631763981, Final residual = 6.27878556034e-06, No Iterations 129
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.00981519206159, Final residual = 6.50074384625e-05, No Iterations 95
    time step continuity errors : sum local = 5.76258113924e-11, global = 3.76596506058e-11, cumulative = 2.12788820726e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000194046206327, Final residual = 1.08198900443e-06, No Iterations 137
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.43945734601e-05, Final residual = 8.53757949514e-10, No Iterations 262
    time step continuity errors : sum local = 3.60665856703e-15, global = -2.65475508792e-15, cumulative = 2.1278855525e-09
    smoothSolver:  Solving for omega, Initial residual = 7.30361462036e-07, Final residual = 3.97781225689e-10, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 3.49789830215e-05, Final residual = 3.4609685557e-09, No Iterations 2
    PIMPLE: iteration 2
    smoothSolver:  Solving for alpha.water, Initial residual = 8.16330553355e-10, Final residual = 8.16330553355e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338827  Min(alpha.water) = -0.000288619975557  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338696  Min(alpha.water) = -0.00028861608681  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.64320126049e-05, Final residual = 3.53085923555e-07, No Iterations 139
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.6562573215e-05, Final residual = 1.63348260399e-07, No Iterations 126
    time step continuity errors : sum local = 6.65506742196e-13, global = 6.65113179978e-13, cumulative = 2.12855066568e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 7.05305597495e-06, Final residual = 4.86049614563e-08, No Iterations 143
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.0062961958e-06, Final residual = 3.82850298196e-10, No Iterations 202
    time step continuity errors : sum local = 1.20177987263e-15, global = 2.18943311055e-16, cumulative = 2.12855088463e-09
    smoothSolver:  Solving for omega, Initial residual = 5.22965286878e-10, Final residual = 5.22965286878e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 8.764653366e-09, Final residual = 8.764653366e-09, No Iterations 0
    PIMPLE: iteration 3
    smoothSolver:  Solving for alpha.water, Initial residual = 8.34145135352e-10, Final residual = 8.34145135352e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338696  Min(alpha.water) = -0.00028861608681  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338578  Min(alpha.water) = -0.000288612326428  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.26546928799e-06, Final residual = 2.26376259356e-08, No Iterations 133
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 7.73727028489e-07, Final residual = 7.58193596374e-09, No Iterations 142
    time step continuity errors : sum local = 3.32815134752e-14, global = 3.32775001386e-14, cumulative = 2.12858416213e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.24877961418e-07, Final residual = 2.65040195196e-09, No Iterations 151
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.36702609294e-07, Final residual = 9.73556700389e-10, No Iterations 142
    time step continuity errors : sum local = 4.50942836308e-15, global = 3.24662375376e-15, cumulative = 2.12858740875e-09
    smoothSolver:  Solving for omega, Initial residual = 5.24948390428e-10, Final residual = 5.24948390428e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 8.92548068881e-09, Final residual = 8.92548068881e-09, No Iterations 0
    PIMPLE: iteration 4
    smoothSolver:  Solving for alpha.water, Initial residual = 8.39608653245e-10, Final residual = 8.39608653245e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338578  Min(alpha.water) = -0.000288612326428  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338472  Min(alpha.water) = -0.000288608472331  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.67244877187e-07, Final residual = 2.53807729896e-09, No Iterations 156
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.94717131277e-08, Final residual = 8.39554535244e-10, No Iterations 127
    time step continuity errors : sum local = 3.85428841404e-15, global = 2.11240953212e-15, cumulative = 2.12858952116e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.8066547318e-08, Final residual = 6.21581241679e-10, No Iterations 95
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 6.6868460246e-09, Final residual = 9.95083309543e-10, No Iterations 14
    time step continuity errors : sum local = 3.31155029848e-15, global = 2.63481371504e-16, cumulative = 2.12858978464e-09
    smoothSolver:  Solving for omega, Initial residual = 5.25158396946e-10, Final residual = 5.25158396946e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 8.94269523662e-09, Final residual = 8.94269523662e-09, No Iterations 0
    PIMPLE: iteration 5
    smoothSolver:  Solving for alpha.water, Initial residual = 8.49256382234e-10, Final residual = 8.49256382234e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338472  Min(alpha.water) = -0.000288608472331  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338375  Min(alpha.water) = -0.000288604542788  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.63410214157e-07, Final residual = 1.58076469785e-09, No Iterations 115
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.91057467543e-09, Final residual = 9.88772998935e-10, No Iterations 27
    time step continuity errors : sum local = 4.32969109332e-15, global = 3.26463782778e-15, cumulative = 2.12859304928e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 6.98750307525e-09, Final residual = 9.88518391293e-10, No Iterations 14
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.44234940454e-09, Final residual = 9.8025747263e-10, No Iterations 4
    time step continuity errors : sum local = 3.90688904253e-15, global = 3.55096694681e-15, cumulative = 2.12859660025e-09
    smoothSolver:  Solving for omega, Initial residual = 5.25212396357e-10, Final residual = 5.25212396357e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 8.94951618755e-09, Final residual = 8.94951618755e-09, No Iterations 0
    ExecutionTime = 2032.89 s  ClockTime = 2056 s
    
    Courant Number mean: 6.42988902653e-07 max: 0.749137867548
    Interface Courant Number mean: 7.49720100137e-09 max: 0.0185446989986
    deltaT = 0.000122067132065
    Time = 0.0125595
    
    PIMPLE: iteration 1
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999947421749)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 -0.00946196445189)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503424141) (11.2436789924 17.5000299924 10.2504535706)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    Updating shallowWaterAbsorption wave model for patch inlet
    Updating shallowWaterAbsorption wave model for patch outlet
    Updating shallowWaterAbsorption wave model for patch rightSide
    Updating shallowWaterAbsorption wave model for patch leftSide
    Execution time for mesh.update() = 43.03 s
    smoothSolver:  Solving for alpha.water, Initial residual = 5.62693014026e-10, Final residual = 5.62693014026e-10, No Iterations 0
    Phase-1 volume fraction = 0.66655233832  Min(alpha.water) = -0.000288604542788  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338241  Min(alpha.water) = -0.000288600000297  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.00610895566778, Final residual = 1.76062946381e-05, No Iterations 187
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000464626734755, Final residual = 4.31868001907e-06, No Iterations 135
    time step continuity errors : sum local = 8.24477020959e-11, global = 6.57112660894e-11, cumulative = 2.19430786634e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000453800900833, Final residual = 3.94151614699e-06, No Iterations 147
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000147155197186, Final residual = 5.92856342839e-10, No Iterations 256
    time step continuity errors : sum local = 7.18291390647e-15, global = -4.28128690922e-15, cumulative = 2.19430358505e-09
    smoothSolver:  Solving for omega, Initial residual = 8.2278501676e-07, Final residual = 3.51637527097e-10, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 3.84521991436e-05, Final residual = 9.69284044515e-09, No Iterations 1
    PIMPLE: iteration 2
    smoothSolver:  Solving for alpha.water, Initial residual = 5.18632871574e-10, Final residual = 5.18632871574e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338241  Min(alpha.water) = -0.000288600000297  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338169  Min(alpha.water) = -0.000288587367539  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.26591595287e-05, Final residual = 7.90344473439e-07, No Iterations 112
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.51930809664e-05, Final residual = 3.22380409956e-07, No Iterations 136
    time step continuity errors : sum local = 3.9496037615e-12, global = 3.83624615083e-12, cumulative = 2.1981398312e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.71819636959e-05, Final residual = 1.70602992744e-07, No Iterations 169
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 7.29985597309e-06, Final residual = 3.6900051247e-10, No Iterations 226
    time step continuity errors : sum local = 4.54190924639e-15, global = -2.53540221376e-15, cumulative = 2.1981372958e-09
    smoothSolver:  Solving for omega, Initial residual = 5.7544716098e-10, Final residual = 5.7544716098e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 2.66232952193e-08, Final residual = 3.81501056948e-09, No Iterations 1
    PIMPLE: iteration 3
    smoothSolver:  Solving for alpha.water, Initial residual = 5.38627943033e-10, Final residual = 5.38627943033e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338169  Min(alpha.water) = -0.000288587367539  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338103  Min(alpha.water) = -0.000288576915369  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.90201537023e-06, Final residual = 3.77572757143e-08, No Iterations 139
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.67020370565e-06, Final residual = 1.49001158022e-08, No Iterations 133
    time step continuity errors : sum local = 1.86485260925e-13, global = 1.62184431696e-13, cumulative = 2.19829948023e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 8.50239811927e-07, Final residual = 3.72560315988e-09, No Iterations 135
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.74883768517e-07, Final residual = 5.88727468347e-10, No Iterations 143
    time step continuity errors : sum local = 6.87566012445e-15, global = 2.61279932851e-15, cumulative = 2.19830209303e-09
    smoothSolver:  Solving for omega, Initial residual = 5.82349200687e-10, Final residual = 5.82349200687e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 5.00557298894e-09, Final residual = 5.00557298894e-09, No Iterations 0
    PIMPLE: iteration 4
    smoothSolver:  Solving for alpha.water, Initial residual = 5.67042728486e-10, Final residual = 5.67042728486e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338103  Min(alpha.water) = -0.000288576915369  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552338043  Min(alpha.water) = -0.000288566342977  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.27362568217e-07, Final residual = 2.24680912268e-09, No Iterations 148
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.01785649227e-08, Final residual = 9.24088431816e-10, No Iterations 131
    time step continuity errors : sum local = 1.10699528045e-14, global = 6.27749363094e-15, cumulative = 2.19830837052e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 4.51230608273e-08, Final residual = 9.6696452405e-10, No Iterations 117
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.92996272797e-08, Final residual = 7.57883021183e-10, No Iterations 102
    time step continuity errors : sum local = 8.10115202691e-15, global = 8.02059882017e-15, cumulative = 2.19831639112e-09
    smoothSolver:  Solving for omega, Initial residual = 5.82669063188e-10, Final residual = 5.82669063188e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 5.04746641437e-09, Final residual = 5.04746641437e-09, No Iterations 0
    PIMPLE: iteration 5
    smoothSolver:  Solving for alpha.water, Initial residual = 5.98420767334e-10, Final residual = 5.98420767334e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552338043  Min(alpha.water) = -0.000288566342977  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552337988  Min(alpha.water) = -0.000288555539222  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.69188395856e-08, Final residual = 9.84231165019e-10, No Iterations 51
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.26446402637e-09, Final residual = 9.89036434404e-10, No Iterations 4
    time step continuity errors : sum local = 1.00675233117e-14, global = 4.92300885785e-15, cumulative = 2.19832131413e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.20396438674e-09, Final residual = 9.80998673529e-10, No Iterations 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.89421021316e-10, Final residual = 9.89421021316e-10, No Iterations 0
    time step continuity errors : sum local = 1.01490028677e-14, global = 4.93624069931e-15, cumulative = 2.19832625037e-09
    smoothSolver:  Solving for omega, Initial residual = 5.82678554279e-10, Final residual = 5.82678554279e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 5.04987373213e-09, Final residual = 5.04987373213e-09, No Iterations 0
    ExecutionTime = 2526.09 s  ClockTime = 2556 s
    
    
        From virtual bool Foam::regIOobject::readIfModified()
        in file db/regIOobject/regIOobjectRead.C at line 271
            Re-reading object fvSchemes from file "/mnt/d/Tutorials/UMaine/decay/decayRASkOmegaSST/decayWECMiddleSpar/decayHeave1V0/background/system/fvSchemes"
    
        From virtual bool Foam::regIOobject::readIfModified()
        in file db/regIOobject/regIOobjectRead.C at line 271
            Re-reading object fvSolution from file "/mnt/d/Tutorials/UMaine/decay/decayRASkOmegaSST/decayWECMiddleSpar/decayHeave1V0/background/system/fvSolution"
    Courant Number mean: 4.99784351636e-07 max: 0.569730783454
    Interface Courant Number mean: 5.92138108728e-09 max: 0.0209198035939
    deltaT = 0.00014349253782
    Time = 0.012703
    
    PIMPLE: iteration 1
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    forces forces:
        rho: rho
        Not including porosity effects
    
    Rigid-body motion of the dynamicWECMiddle
        Centre of rotation: (0 0 -0.999948522777)
        Orientation: (1 0 0 0 1 0 0 0 1)
        Linear velocity: (0 0 -0.00588418016272)
        Angular velocity: (0 0 0)
    inverseDistance : detected 3 mesh regions
        zone:0 nCells:3302728  voxels:(92 92 92) bb:(-1250.00076811 -1250.00076811 -800.000768115) (-749.999231885 -749.999231885 -499.999231885)
        zone:1 nCells:1329243  voxels:(92 92 92) bb:(-17.5000567663 -17.5000567663 -25.0007917663) (17.5000557663 17.5000547663 -1.56244523367)
        zone:2 nCells:1783736  voxels:(92 92 92) bb:(-4.41261149239 5.32858310761 -12.2503435152) (11.2436789924 17.5000299924 10.2504524696)
    Overset analysis : nCells : 6415707
        calculated   : 5795510
        interpolated : 159835 (from local:606  mixed local/remote:100  remote:159129)
        hole         : 455019
    
    Updating shallowWaterAbsorption wave model for patch inlet
    Updating shallowWaterAbsorption wave model for patch outlet
    Updating shallowWaterAbsorption wave model for patch rightSide
    Updating shallowWaterAbsorption wave model for patch leftSide
    Execution time for mesh.update() = 41.31 s
    smoothSolver:  Solving for alpha.water, Initial residual = 4.31650457418e-10, Final residual = 4.31650457418e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552337956  Min(alpha.water) = -0.000288555539222  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552337912  Min(alpha.water) = -0.000288542572206  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.00236165217602, Final residual = 2.3271010307e-05, No Iterations 107
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 0.000412146953105, Final residual = 4.00471828788e-06, No Iterations 119
    time step continuity errors : sum local = 5.53488316472e-11, global = -3.70246236732e-11, cumulative = 2.1613016267e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.71606087261e-05, Final residual = 5.21097796183e-07, No Iterations 138
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.35125965116e-05, Final residual = 9.88359006305e-10, No Iterations 215
    time step continuity errors : sum local = 1.51332173524e-14, global = 4.52701018796e-15, cumulative = 2.16130615371e-09
    smoothSolver:  Solving for omega, Initial residual = 9.65248956233e-07, Final residual = 3.75230694262e-10, No Iterations 1
    smoothSolver:  Solving for k, Initial residual = 4.42745824449e-05, Final residual = 4.43358495994e-09, No Iterations 2
    PIMPLE: iteration 2
    smoothSolver:  Solving for alpha.water, Initial residual = 5.14219964389e-10, Final residual = 5.14219964389e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552337912  Min(alpha.water) = -0.000288542572206  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552337866  Min(alpha.water) = -0.000288456349238  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.58024318209e-05, Final residual = 1.50927926767e-07, No Iterations 135
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.73850293211e-06, Final residual = 5.74727092462e-09, No Iterations 165
    time step continuity errors : sum local = 6.69721692409e-14, global = -3.03450383446e-14, cumulative = 2.16127580867e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.49876301523e-06, Final residual = 1.39670720646e-08, No Iterations 123
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 6.76715177075e-07, Final residual = 7.98008738186e-10, No Iterations 148
    time step continuity errors : sum local = 1.22344536329e-14, global = -5.78068108562e-15, cumulative = 2.16127002799e-09
    smoothSolver:  Solving for omega, Initial residual = 6.32837304089e-10, Final residual = 6.32837304089e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 1.53271768735e-08, Final residual = 2.46952790962e-09, No Iterations 1
    PIMPLE: iteration 3
    smoothSolver:  Solving for alpha.water, Initial residual = 5.76591398138e-10, Final residual = 5.76591398138e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552337866  Min(alpha.water) = -0.000288456349238  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552337822  Min(alpha.water) = -0.000288366845858  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.77537562021e-07, Final residual = 5.48225168085e-09, No Iterations 102
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.56677431916e-07, Final residual = 1.4141887376e-09, No Iterations 140
    time step continuity errors : sum local = 2.23087432416e-14, global = -1.32089190452e-15, cumulative = 2.1612687071e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 6.59080063594e-08, Final residual = 9.9295509487e-10, No Iterations 132
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 3.25757162506e-08, Final residual = 9.37883485476e-10, No Iterations 103
    time step continuity errors : sum local = 1.46786743791e-14, global = -1.4653311533e-14, cumulative = 2.16125405379e-09
    smoothSolver:  Solving for omega, Initial residual = 6.49111444261e-10, Final residual = 6.49111444261e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 4.0118490653e-09, Final residual = 4.0118490653e-09, No Iterations 0
    PIMPLE: iteration 4
    smoothSolver:  Solving for alpha.water, Initial residual = 6.1973337999e-10, Final residual = 6.1973337999e-10, No Iterations 0
    Phase-1 volume fraction = 0.666552337822  Min(alpha.water) = -0.000288366845858  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.66655233778  Min(alpha.water) = -0.000288274846298  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 5.74398048665e-08, Final residual = 9.82136138533e-10, No Iterations 120
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 7.02405804242e-09, Final residual = 9.98479118902e-10, No Iterations 9
    time step continuity errors : sum local = 1.27292028318e-14, global = -7.9193313664e-16, cumulative = 2.16125326185e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 2.30850519151e-09, Final residual = 9.51027282588e-10, No Iterations 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.75182912812e-10, Final residual = 9.75182912812e-10, No Iterations 0
    time step continuity errors : sum local = 1.26570808677e-14, global = -4.64721120315e-16, cumulative = 2.16125279713e-09
    smoothSolver:  Solving for omega, Initial residual = 6.49171974338e-10, Final residual = 6.49171974338e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 4.0055169094e-09, Final residual = 4.0055169094e-09, No Iterations 0
    PIMPLE: iteration 5
    smoothSolver:  Solving for alpha.water, Initial residual = 6.63999797793e-10, Final residual = 6.63999797793e-10, No Iterations 0
    Phase-1 volume fraction = 0.66655233778  Min(alpha.water) = -0.000288274846298  Max(alpha.water) = 1
    MULES: Correcting alpha.water
    MULES: Correcting alpha.water
    Phase-1 volume fraction = 0.666552337738  Min(alpha.water) = -0.00028818033445  Max(alpha.water) = 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.61804964183e-08, Final residual = 9.64384311072e-10, No Iterations 8
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.73777456394e-09, Final residual = 9.66855572087e-10, No Iterations 1
    time step continuity errors : sum local = 1.11810878112e-14, global = 1.93801378722e-15, cumulative = 2.16125473514e-09
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 1.14956370447e-09, Final residual = 9.54095995662e-10, No Iterations 1
    DILUPBiCGStab:  Solving for p_rgh, Initial residual = 9.55825543332e-10, Final residual = 9.55825543332e-10, No Iterations 0
    time step continuity errors : sum local = 1.1122693211e-14, global = 1.9981078044e-15, cumulative = 2.16125673325e-09
    smoothSolver:  Solving for omega, Initial residual = 6.49156613018e-10, Final residual = 6.49156613018e-10, No Iterations 0
    smoothSolver:  Solving for k, Initial residual = 4.00623653355e-09, Final residual = 4.00623653355e-09, No Iterations 0
    ExecutionTime = 2914.54 s  ClockTime = 2959 s
    
    1 条回复 最后回复
  • H 离线
    H 离线
    hachikoi
    写于 最后由 编辑
    #4

    背景网格的插值区,计算区和hole
    image.png
    中心柱overset的插值区,计算区和hole
    image.png
    外环柱overset的插值区,计算区和hole

    image.png

    1 条回复 最后回复
  • H 离线
    H 离线
    hachikoi
    写于 最后由 编辑
    #5

    大佬有没有什么建议呀 @队长别开枪

    队长别开枪队 1 条回复 最后回复
  • 队长别开枪队 离线
    队长别开枪队 离线
    队长别开枪 超神
    在 中回复了 hachikoi 最后由 编辑
    #6

    @hachikoi 还没做过多浮体的例子,最近做了一个overset+AMR+PLIC的算例,结果还行。你的算例几何尺寸可以分享一下么,圣诞节放假的时候我可以用我的求解器试着算算。

    H 1 条回复 最后回复
  • H 离线
    H 离线
    hachikoi
    在 中回复了 队长别开枪 最后由 编辑
    #7

    @队长别开枪
    内柱:Cylinder z1(0 0 -20) z2 (0 0 15) 半径 6.25
    外环:z1(0 0 -7.25) z2(0 0 5.25) 内半径 7.5 外半径 12.5
    外环原先是关于水面对称,下降一米让其上浮

    decayHeave1V0.rar

    麻烦大佬给看看

    据另一位做六自由度双浮箱的同学说,他也遇到了类似的问题 只要overset之间相互接触就会发散,不知问题到底出在哪

    Z S 2 条回复 最后回复
  • 王 离线
    王 离线
    王长宏
    写于 最后由 编辑
    #8

    @hachikoi 您好,能帮忙问一下做六自由度双浮箱的同学,他做的双浮箱涉及两个浮体之间的耦合吗,想交流一下。

    1 条回复 最后回复
  • Z 离线
    Z 离线
    zbxia
    在 中回复了 hachikoi 最后由 编辑
    #9

    @hachikoi 你好,请问你说的外浮体overset➕内浮体背景网格扣除,内浮体通过什么工具在背景网格里面扣除呢?最近遇到类似的问题,卡住了,我的案例是浮体通过overset实现,背景网格内紧邻浮体有一个固定式结构,我用Toposet和sanppyHexMesh工具在背景网格里扣除固定式结构时,第一个时间步计算压力就发散,尝试修改网格尺寸,浮体与背景网格固定式结构距离,压力求解方式等均未解决,谢谢回答!

    1 条回复 最后回复
  • S 离线
    S 离线
    Sichan
    在 中回复了 hachikoi 最后由 编辑
    #10

    @hachikoi 您好,请问您的问题现在解决了吗?我现在也在做和您相似的双浮体模拟,但是在模拟时出现了第一部就会发散的问题,重叠网格挖的洞数量很大。

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]