[treem22]

I'm attempting to create a simple surface wave and watch it develop from inlet to outlet, but am running across an issue when the wave model is attempting to update. The following is my interDyMFoam log file:

Courant Number mean: 0 max: 0
 70 Interface Courant Number mean: 0 max: 0
 71 Time = 0.001
 72
 73 PIMPLE: iteration 1
 74 Selected 34033 cells for refinement out of 625000.
 75 Refined from 625000 to 863231 cells.
 76 Selected 0 split points out of a possible 34033.
 77 Execution time for mesh.update() = 8.93 s
 78 Selecting waveModel StokesII
 79 Wave model: patch inlet
 80     Type : StokesII
 81     Velocity field name : U
 82     Phase fraction field name : alpha.water
 83     Transformation from local to global system : (1 -0 -0 0 1 0 -0 -0 1)
 84     Number of paddles: 1
 85     Reference water depth : 0.7576
 86     Active absorption: 1
 87     Wave height : 0.1
 88     Wave angle  : 0
 89     Ramp time : 3
 90     Wave period : 3
 91     Wave length : 7.71495
 92     Wave phase : 4.71239
 93     Wave type: Intermediate
 94 Updating StokesII wave model for patch inlet
 95 Selecting waveModel shallowWaterAbsorption
 96 Wave model: patch outlet
 97     Type : shallowWaterAbsorption
 98     Velocity field name : U
 99     Phase fraction field name : alpha.water
100     Transformation from local to global system : (-1 -0 -0 0 -1 0 -0 -0 1)
101     Number of paddles: 1
102     Reference water depth : 0.756
103     Active absorption: 1
104 Updating shallowWaterAbsorption wave model for patch outlet
105 DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 9.37693e-07, No Iterations 223
106 time step continuity errors : sum local = 1.73663e-14, global = -2.5422e-16, cumulative = -2.5422e-16
107 MULES: Solving for alpha.water
108 Phase-1 volume fraction = 0.75613  Min(alpha.water) = 0  Max(alpha.water) = 1
109 Updating StokesII wave model for patch inlet
110 MULES: Solving for alpha.water
111 Phase-1 volume fraction = 0.75613  Min(alpha.water) = -3.14696e-06  Max(alpha.water) = 1
112 Updating StokesII wave model for patch inlet
113 MULES: Solving for alpha.water
114 Phase-1 volume fraction = 0.75613  Min(alpha.water) = -6.29392e-06  Max(alpha.water) = 1
115 Updating StokesII wave model for patch inlet
116 DICPCG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.0993804, No Iterations 1
117 time step continuity errors : sum local = 0.0368612, global = 1.74905e-08, cumulative = 1.74905e-08
118 GAMG:  Solving for p_rgh, Initial residual = 0.0484055, Final residual = 9.00315e-08, No Iterations 7
119 time step continuity errors : sum local = 6.85691e-08, global = -9.55399e-09, cumulative = 7.93653e-09
120 ExecutionTime = 39.19 s  ClockTime = 39 s

Courant Number mean: 0.00132011 max: 0.21225
123 Interface Courant Number mean: 0 max: 0
124 Time = 0.002
125
126 PIMPLE: iteration 1
127 Selected 149175 cells for refinement out of 863231.
128 Refined from 863231 to 1907456 cells.
129 Selected 0 split points out of a possible 149175.
130 Execution time for mesh.update() = 24.82 s
131 Updating StokesII wave model for patch inlet
132 [0] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/lib/libOpenFOAM.so"
133 [0] #1  Foam::sigSegv::sigHandler(int) in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/lib/libOpenFOAM.so"
134 [0] #2  ? in "/lib64/libc.so.6"
135 [0] #3  Foam::waveModel::waterLevel() const in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/lib/libwaveModels.so"
136 [0] #4  Foam::waveModel::correct(double) in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/lib/libwaveModels.so"
137 [0] #5  Foam::waveVelocityFvPatchVectorField::updateCoeffs() in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/lib/libwaveModels.so"
138 [0] #6  Foam::fvPatchField<Foam::Vector<double> >::evaluate(Foam::UPstream::commsTypes) in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/bin/interDyMFoam"
139 [0] #7  Foam::correctUphiBCs(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&) in "/opt/Software/openfoam/161    2+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/lib/libfiniteVolume.so"
140 [0] #8  ? in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/bin/interDyMFoam"
141 [0] #9  ? in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/bin/interDyMFoam"
142 [0] #10  __libc_start_main in "/lib64/libc.so.6"
143 [0] #11  ? in "/opt/Software/openfoam/1612+/OpenFOAM-v1612+/platforms/linux64IccDPInt32Opt/bin/interDyMFoam"

Here is my fvSolution file:

solvers
 20 {
 21     "alpha.water.*"
 22     {
 23         nAlphaCorr      1;
 24         nAlphaSubCycles 3;
 25         cAlpha          1;
 26     }
 27
 28     pcorr
 29     {
 30         solver          PCG;
 31         preconditioner  DIC;
 32         tolerance       1e-6;
 33         relTol          0;
 34     }
 35
 36     pcorrFinal
 37     {
 38         solver          PCG;
 39         preconditioner  DIC;
 40         tolerance       1e-6;
 41         relTol          0;
 42     }
 43
 44     p_rgh
 45     {
 46         solver          PCG;
 47         preconditioner  DIC;
 48         tolerance       1e-6;
 49         relTol          0.1;
 50     }
 51
 52     p_rghFinal
 53     {
 54         solver          GAMG;
 55         smoother        DIC;
 56         tolerance       1e-7;
 57         relTol          0;
 58     }
 59
 60     U
 61     {
 62         solver          PBiCG;
 63         preconditioner  DILU;
 64         tolerance       1e-6;
 65         relTol          0.1;
 66     }
 67
 68     UFinal
 69     {
 70         solver          PBiCG;
 71         preconditioner  DILU;
 72         tolerance       1e-6;
 73         relTol          0;
 74     }
 75 }
 76
 77 PIMPLE
 78 {
 79     momentumPredictor no;
 80     nCorrectors     2;
 81     nNonOrthogonalCorrectors 0;
 82 }

Here is my fvSchemes file:

ddtSchemes
 20 {
 21     default         Euler;
 22 }
 23
 24 gradSchemes
 25 {
 26     default         Gauss linear;
 27 }
 28
 29 divSchemes
 30 {
 31     div(rhoPhi,U)   Gauss linearUpwind grad(U);
 32     div(phi,alpha)  Gauss vanLeer;
 33     div(phirb,alpha) Gauss linear;
 34     div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
 35 }
 36
 37 laplacianSchemes
 38 {
 39     default         Gauss linear orthogonal;
 40 }
 41
 42 interpolationSchemes
 43 {
 44     default         linear;
 45 }
 46
 47 snGradSchemes
 48 {
 49     default         orthogonal;
 50 }

The domain is a structured rectangular prism that is 1m x 1m x 5m. The two 1m x 1m faces on the ends are the inlet and outlet. The initial water/air interface is at y=0.75m from the ground. An image is attached.


Any help will be most appreciated. Thanks!