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

hachikoi

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

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

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

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

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;
    }
}


// ************************************************************************* //

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

hachikoi

计算结果

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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

hachikoi

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

hachikoi

背景网格的插值区，计算区和hole

中心柱overset的插值区，计算区和hole

外环柱overset的插值区，计算区和hole

hachikoi

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

队长别开枪

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

hachikoi

@队长别开枪
内柱：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之间相互接触就会发散，不知问题到底出在哪

王长宏

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

zbxia

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

CFD中文网

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

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