[shimakasaei]

Hi,
I am using olafoam to simulate my model with the inlet of irregular wave.
I did set all of the files such as U, P_rgh, alpha.water according to my model.
when start to run a simulation (running command: olaFlow), I got an error and I do not know how to solve it!!


this is the error:

PIMPLE: iteration 1
MULES: Solving for alpha.water
alpha.water BC on patch inlet
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib64/libc.so.6"
#3 StokesIFun::waveLength(double, double) at ??:?
#4 Foam::waveAlphaFvPatchScalarField::updateCoeffs() at ??:?
#5 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() at ??:?
#6 ? at ??:?
#7 ? at olaFlow.C:?
#8 ? at ??:?
#9 __libc_start_main in "/lib64/libc.so.6"
#10 ? at ??:?
Floating point exception (core dumped)


can anyone please help me through this??

Best,
shima

[Phicau]

Hi Shima,

sure, I can help you with this, can you upload your case so that I can take a look?

Just a couple of quick checks that can guide you here and may solve your problem:
- Make sure that the wave period for all your components is different from 0.
- Have you performed setFields?

Best,

Pablo

[shimakasaei]

Hi pablo,
Thanks for the reply, this is my alpha.water file.

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type waveAlpha;
waveDictName waveDict;
value uniform 0.0;
}
outlet
{
type zeroGradient;
}
leftwall
{
type zeroGradient;
}
rightwall
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}
top
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
A3b1
{
type zeroGradient;
}
A3b2
{
type zeroGradient;
}
A3b3
{
type zeroGradient;
}


}


-------------------------
and for the waveDict, yes I am sure that I have waveHeights, wavePeriods, wavePhases and waveDirs and none of the wavePeriods are zero...!

for setFields, no I did not use that...I used blockMesh, surfaceFeatureExtract and snappyHexMesh..... should I use that also???


thanks for the help
shima

[Phicau]

Hi Shima,

yes, setFields is needed to fill your mesh with water. Make sure the box in setFieldsDict covers the whole mesh at the level that you want to set water at.

The BC is failing because there is no water at the boundary!

Best,

Pablo

[shimakasaei]

Originally Posted by Phicau

Hi Shima,

yes, setFields is needed to fill your mesh with water. Make sure the box in setFieldsDict covers the whole mesh at the level that you want to set water at.

The BC is failing because there is no water at the boundary!

Best,

Pablo

I did run setFeilds and it showed me warning:

Reading setFieldsDict

Setting field default values
--> FOAM Warning :
From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
in file setFields.C at line 121
Field alpha.water not found

Setting field region values
Adding cells with center within boxes 1((0 0 0) (10 1.5 0.3))
--> FOAM Warning :
From function bool setCellFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
in file setFields.C at line 121
Field alpha.water not found

End

--------------
and the same error when I tried to run olaFlow!

the whole flume dimension is (10 1.5 0.75) and the box in setFeildDict is (10 1.5 0.30) which contains 30 cm water in the flume.

I do not know what is the BC problem...!!! I will appreciate your help.


Best,
shima

[Phicau]

Hi Shima,

the error in olaFlow is the same because nothing has changed. The warning in setFields is self-explanatory. It indicates that the file alpha.water has not been found in the 0 folder, so it has not been filled with water. Make sure that you have a 0 folder and that alpha.water is inside before running setFields.

I suggest that you try to understand the steps needed to run a case and the correct order. Take a look at the baseWaveFlume tutorial included in olaFlow. It has a runCase script with all the commands needed to run the case.

Best,

Pablo

[shimakasaei]

Originally Posted by Phicau

Hi Shima,

the error in olaFlow is the same because nothing has changed. The warning in setFields is self-explanatory. It indicates that the file alpha.water has not been found in the 0 folder, so it has not been filled with water. Make sure that you have a 0 folder and that alpha.water is inside before running setFields.

I suggest that you try to understand the steps needed to run a case and the correct order. Take a look at the baseWaveFlume tutorial included in olaFlow. It has a runCase script with all the commands needed to run the case.

Best,

Pablo

thank you so much, I think I find the problem:
defaultFieldValues
(
volScalarFieldValue alpha.water 0
volScalarFieldValue porosityIndex 0******
volVectorFieldValue U (0 0 0)****
);

regions
(
boxToCell
{
box ( 0 0.000 0.000 ) ( 10 1.5 0.30 );
fieldValues ( volScalarFieldValue alpha.water 1 );****
}


);

when I added the parts that are shown with *****, it does work!

thank you so much
shima

[Mariano_]

Dear colleges,
I am recently starting with OpenFoam and OlaFlow. I have already managed to compile and run the examples of the olaFlow tutorials. I ask if anyone can provide me with the files to plot the example breakwater and basewaveflume for python 3. It would help me a lot for these first steps. Thank you very much!


Mariano.

[Phicau]

Hi Mariano,


I have added python3 version of the postprocessing scripts in the reference folder. Check it out and download the latest version of the repository: https://github.com/phicau/olaFlow


Best,


Pablo

[Mariano_]

Hi Pablo,
thank you so much for your time and work!



Best,


Mariano.