overset 小间隙-多浮体算例 计算结果明显有误
-
最近在做同轴双浮体装置,目前是内浮体不动(后续会动,所以做成了overset),外浮体初始位移为-1做自衰减
内浮体的overset较大,将外浮体overset包裹,四面消波处理
目前计算结果,外浮体在最初的几秒存在负速度,这是完全不合理的,之前尝试过单独外浮体/单独内浮体/外浮体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; } } // ************************************************************************* //
请大家帮忙看看是不是哪里设置有问题
-
计算结果
/*---------------------------------------------------------------------------*\ | ========= | | | \\ / 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
-
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