[Stephy]

Hello all,
I am trying to run the wigley hull tutorial on OpenFoam 2.0.0.
I have built the mesh using snappyHexMesh and I ran the LTSInterFoam computations. No error, eveything went well.
But then, when I run paraFoam, I can see the mesh, but when playing and displaying any vector like U or p, they would mysteriously disappear (I cannot select them instead of "Solid color"), and when clicking "apply" again, paraView closes automatically and I receive the following error message :
"
created temporary 'wigleyHull.OpenFOAM'
--> FOAM Warning :
From function polyMesh::readUpdateState polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 204
Number of patches has changed. This may have unexpected consequences. Proceed with care.


--> FOAM FATAL IO ERROR:
size 12000 is not equal to the given value of 182759

file: /home/alex/OpenFOAM/stephy-2.0.0/run/tutorials/multiphase/LTSInterFoam/wigleyHull/100/nut from line 18 to line 12043.

From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
in file lnInclude/Field.C at line 236.

FOAM exiting
"

I checked the files mentionned, but I couldn't find where this value of 182759 comes from, I couldn't find it in any file.

Could you please help me if you have any idea of what could solve the problem ?

Thank you very much,

Stephy

[Curro5150]

Try renaming the "0" folder to something like "0.org" and then renaming the "2" folder to "0". Otherwise openFoam (or maybe paraFoam) gets somehow confused with the original mesh before using snappyHexMesh. That solved the problem for me.

Regards,

Francisco

[jameel]

Originally Posted by Curro5150

Try renaming the "0" folder to something like "0.org" and then renaming the "2" folder to "0". Otherwise openFoam (or maybe paraFoam) gets somehow confused with the original mesh before using snappyHexMesh. That solved the problem for me.

Regards,

Francisco

Dear Francisco,
I have the same problem above, I tried to use your solution and I changed the zero file name to 0.org. but what do you mean by( the"2" folder ). please write to me
atheel,

[jorkolino]

Originally Posted by jameel

Dear Francisco,
I have the same problem above, I tried to use your solution and I changed the zero file name to 0.org. but what do you mean by( the"2" folder ). please write to me
atheel,

Hi,

maybe I can help. When constructing the mesh, snappy does it in several steps (1st castellating the mesh, then 2. snapping the mesh, then 3. adding layers). Each stage is written into a numbered folder, starting from 1. (0 goes for the blockMesh I suppose). So your final mesh should be the latest number (it can be number 3 or 1 instead of 2, depending on which steps you employed in the meshing process). So follow the instructions above with your latest folder number instead of number 2. Regards

[sandy13]

Originally Posted by jorkolino

Hi,

maybe I can help. When constructing the mesh, snappy does it in several steps (1st castellating the mesh, then 2. snapping the mesh, then 3. adding layers). Each stage is written into a numbered folder, starting from 1. (0 goes for the blockMesh I suppose). So your final mesh should be the latest number (it can be number 3 or 1 instead of 2, depending on which steps you employed in the meshing process). So follow the instructions above with your latest folder number instead of number 2. Regards

Dear jorkolino,
Thank you so much for this explanation. It was very helpful,

[jovenitta]

Hi all. I am a newbie and trying to run wigley hull in OpenFoam 2.3. I get the same error message. I did exactly what has been told but the end of the message doesn't change. When I try to apply the volume fields in paraview (4.1). It suddenly exits and I get the error message below:
--> FOAM FATAL IO ERROR:
size 12000 is not equal to the given value of 182759

file: /home/itu/OpenFOAM/itu-2.3.0/run/tutorials/multiphase/LTSInterFoam/wigleyHull/2000/nut from line 18 to line 12047.

From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
in file /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.C at line 292.

FOAM exiting

Segmentation fault (core dumped)

This is quite an old thread but is there still anybody can help me with this issue? Thanks..

[jovenitta]

Somehow I couldn't solve wigleyHull turorial using OF 2.3. I installed OF 2.2.2 instead and now it works.

[Sachin m]

Hi Jovenitta,

Try the following steps and it will work.

1) first use the blockMesh command
2) snappyHexMesh -overwrite
3) renumberMesh -overwrite
4) foamToVTK
5) LTSInterFoam

[kubashmuba]

Hi,

I tried what Sachin m suggested, and it works once I have copied alpha1.org into alpha1. I have OpenFOAM 2.1.1 installed and took the wigleyHull tutorial from 2.0.0 I believe, since it doesn't seem to feature in versions after that. Running paraFoam and creating a contour plot for alpha1 at 0.5 to visualize the free surface gives nothing though, and the info in the contour plot for [min, max] of alpha1 says [0, 0]. So there seems to be a problem somewhere. I am new at openFOAM so I'm not even sure where to look for an error... What additional information could I provide to get help?

Thank you

[Sachin m]

I am using 2.1.0

once you have completed the calculations, open parafoam. tick the tab on the left bottom corner which says alpha1. then on the active variable control tool bar , change the tab showing 'p' to alpha1.
Then on the common filters tab click on contour. You ll get the free surface.

[kubashmuba]

Ah I see. I'm using OpenFOAM 2.1.1 because I'm on Ubuntu 14.04, and that was the closest version to 2.1.0, which can only be installed on Ubuntu 10 or 11.

I'm using a virtual box so I guess I can start from scratch again so I can install 2.1.0. First I'll try to compile 2.1.0 so I don't have to reinstal Ubuntu, maybe that will work, doesn't anybody know if it's possible to compile 2.1.0 for Ubuntu 14.04?

But I'm wondering why a minor version bump (2.1.0 -> 2.1.1) would break the wigley hull tutorial that badly, and why is it that I get [0,0] for alpha1. The computation didn't seem to go wrong or crash, just the result seems incorrect. Was there a change in format or to LTSInterFoam that could explain the difference in outcome between 2.1.0 and 2.1.1?

[Sachin m]

Hi all,
there is a slight mistake in the steps that i mentioned. i forgot to mention the setField command.
So here is the steps :

1) first use the blockMesh command
2) snappyHexMesh -overwrite
3) renumberMesh -overwrite
4) foamToVTK
5) setFields
6) LTSInterFoam

[Sachin m]

Hi kubashmuba,

the reason u dont get a free surface is because you did not type the command setFields. Which converts the domain into air and water.In your case you have only one value of alpha; that is zero.
Try the steps that i have mentioned just now.

[kubashmuba]

Hi Sachi m,

I downloaded and compiled 2.1.0 instead, it took me a little while to recompile everything.

- I then copied alpha1.org into alpha1 in the 0 folder
- I also decompressed wigley.stl and wigley-scaled-oriented.stl from $FOAM_TUTORIALS/resources/geometry into wigleyHull/constant/triSurface

I then ran blockMesh with no visible problem
Then I ran snappyHexMesh -overwrite, no apparent problem
Running renumberMesh -overwrite however gave me this error:

-->FOAM FATAL IO ERROR:
keyword hull_wall is undefined in dictionary "/home/.../OpenFOAM/.../run/tutorials/multiphase/LTSInterFoam/wigleyHull/0/nut::boundaryField"

file /home/.../OpenFOAM/.../run/tutorials/multiphase/LTSInterFoam/wigleyHull/0/nut::boundaryField from line 26 to line 29

From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 461.

FOAM exiting

So then I tried rerunning snappyHexMesh and I started getting jovenitta's errors:

size 12000 is not equal to the given value of 182759

So I ran ./Allclean and ran ./Allrun, which doesn't have some of the steps Sachi m listed, but I think this proves that you can get in a bad state and running ./Allclean can help you recover from that.

It did work in the end, so thanks a lot for your help Sachi m.

I see some artefacts when large cells can be seen in the alpha data, so I'm going to try again with your additional renumberMesh and foamToVTK commands. By the way, what do they do?

[vava10]

Hey
I am new to OpenFoam and I am using OpenFoam v-2006
I am working with wigley hull case from https://github.com/OpenFOAM/OpenFOAM...oam/wigleyHull

I followed the following steps

1. blockMesh
2. snappyHexMesh -overwrite
3. transformPoints -scale '(2 1 1)'
4. cp 0/alpha1.org 0/alpha1
5. setFields

when I run setFields I get the following error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2006 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _b45f8f6f58-20200629 OPENFOAM=2006
Arch : "LSB;label=32;scalar=64"
Exec : setFields
Date : Sep 24 2020
Time : 17:54:26
Host : LAPTOP-A0NOGRUJ
PID : 3337
I/O : uncollated
Case : /home/sam/OpenFOAM/OpenFOAM-v2006/tutorials/wigleyHull
nProcs : 1
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

--> FOAM Warning :
From static Foam::IOstreamOption::compressionType Foam::IOstreamOption::compressionEnum(const Foam::word&, Foam::IOstreamOption::compressionType)
in file db/IOstreams/IOstreams/IOstreamOption.C at line 115
Unknown compression specifier 'compressed', using compression off
Create mesh for time = 0

Reading setFieldsDict

Setting field default values
Setting internal values of volScalarField alpha1


--> FOAM FATAL IO ERROR:
Cannot find patchField entry for hull

file: /home/sam/OpenFOAM/OpenFOAM-v2006/tutorials/wigleyHull/0/alpha1.boundaryField at line 26 to 49.

From void Foam::GeometricField<Type, PatchField, GeoMesh>::Boundary::readField(const Foam:imensionedField<TypeR, GeoMesh>&, const Foam::dictionary&) [with Type = double; PatchField = Foam::fvPatchField; GeoMesh = Foam::volMesh]
in file /home/pawan/OpenFOAM/OpenFOAM-v2006/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 172.

FOAM exiting

I don't know why I am getting this error is coming
I will greatly appreciate if anyone could help me with this.

Regards
vava10

[Antimony]

Hi,

Is "hull" present in both alpha1 and constant/polyMesh/boundary? If it is present in one but not the other you will need to ensure that the same name is present in both places.

Based on the link you posted, in the 0 folder, all of them have "hull_wall". But what you will get after snappyHexMesh is "hull" (look at snappyHexMeshDict). So to be consistent either:

a. "hull_wall" should be in place of "hull" in snappyHexMeshDict

Or

b. "hull" should be in place of "hull_wall" in 0/*

Hope this helps.

Cheers,
Antimony

[vava10]

Hey Antimony,

Thank you so much.
It helped. That was the issue. I changed everything to hull
that step was cleared .

but when I ran interform the following happened:


sam@LAPTOP-A0NOGRUJ:~/OpenFOAM/OpenFOAM-v2006/tutorials/wigleyHull$ interFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2006 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _b45f8f6f58-20200629 OPENFOAM=2006
Arch : "LSB;label=32;scalar=64"
Exec : interFoam
Date : Sep 26 2020
Time : 10:38:14
Host : LAPTOP-A0NOGRUJ
PID : 881
I/O : uncollated
Case : /home/sam/OpenFOAM/OpenFOAM-v2006/tutorials/wigleyHull
nProcs : 1
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

--> FOAM Warning :
From static Foam::IOstreamOption::compressionType Foam::IOstreamOption::compressionEnum(const Foam::word&, Foam::IOstreamOption::compressionType)
in file db/IOstreams/IOstreams/IOstreamOption.C at line 115
Unknown compression specifier 'compressed', using compression off
Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Using LTS
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

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.555556;
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
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 1.59522e-06, No Iterations 5
time step continuity errors : sum local = 8.07781e-09, global = 7.43151e-10, cumulative = 7.43151e-10

Starting time loop

Flow time scale min/max = 0.000258129, 0.0900016
Smoothed flow time scale min/max = 0.000258129, 0.00551245
Time = 1

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 4.05791e-06, Final residual = 1.90272e-09, No Iterations 1
Phase-1 volume fraction = 0.925716 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.925716 Min(alpha.water) = -1.14545e-07 Max(alpha.water) = 1


--> FOAM FATAL ERROR:

request for volScalarField alpha1 from objectRegistry region0 failed
available objects of type volScalarField are

17
(
alpha.water_0
interfaceProperties:K
nut
yWall
alpha.water
rho
k
p_rgh
nu
gh
nu1
rDeltaT
p
rho_0
nu2
alpha.air
omega
)


From const Type& Foam:bjectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>]
in file /home/pawan/OpenFOAM/OpenFOAM-v2006/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 463.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::exitOrAbort(int, bool) at ??:?
#2 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam:bjectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&, bool) const at ??:?
#3 Foam:utletPhaseMeanVelocityFvPatchVectorField::u pdateCoeffs() at ??:?
#4 Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<doub le>, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#5 ? at ??:?
#6 ? at ??:?
#7 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#8 ? at ??:?
Aborted (core dumped)]



In the case file which I send you before I did change the name of alpha1 to alpha.water.

I will attach all the files where I made change here.

from this portion of error "request for volScalarField alpha1 from objectRegistry region0 failed
available objects of type volScalarField are" I do get its looking for volScalarField alpha1. I did change the setFieldDict

But I believe I changed everything correctly.I changed all the names by following a youtube tutorial for VOF method. May be I missed something. But I couldn't find it.

I really appreciate your help





I couldn't attach the alpha.water file from 0 folder so I will copy the content below.

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2006 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];


internalField nonuniform List<scalar>
182759
(
1
1
1
1
1
1
1
1
1
1
1
1
....
...
..
1
1
1
1
1
1
1
1
1
)
;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
sides
{
type symmetry;
}
hull
{
type zeroGradient;
}
}


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

[Antimony]

Hi,

The short answer is that alpha1 does not seem to be looked for in the code.

Look at the error message that you got:

--> FOAM FATAL ERROR:

request for volScalarField alpha1 from objectRegistry region0 failed
available objects of type volScalarField are

17
(
alpha.water_0
interfaceProperties:K
nut
yWall
alpha.water
rho
k
p_rgh
nu
gh
nu1
rDeltaT
p
rho_0
nu2
alpha.air
omega
)

It lists what are the available ones - do you have alpha1 or alpha.water? You can make a copy of alpha1 as alpha.water and see if that solves the problem.

If it doesn't then you need to mention what error message you get.

Cheers,
Antimony

[vava10]

Dear Antimony,


I tried what you said. but it is still not working and I get the same error. it is not reading alpha1 even after renaming the alpha file in the 0 folder.

this is how I named alpha in all the files:

• in the transport properties

phases (water air)

water
{
...
}

air
{
..
}

• in set fields

volScalarFieldValue alpha.water 0

• in the 0 folder

alpha.water

( I tried renaming this to alpha1 asper your suggestion, but still got the same error)

• in fvSolution

alpha.water.*


Do you find any mistake in it?

and why are they asking for alpha1 when I changed every alpha files to alpha.water?? which Part of the entire case set up is responsible for asking for volScalarFieldValue alpha1??. If I know if which part of the case is responsible for volScalarFieldValue "alpha1" maybe it become fixable


I am thankful for your help.
I really appreciate it
Thank you in advance for your time.

Kind regards
vava10

[Yann]

Hello vava10,

You're trying to run a tutorial made for OpenFOAM-2.2.x on OpenFOAM-v2006. There is a 7 years time gap between these versions, and a lot of changes in the code.

Is there a reason for you to use this specific tutorial instead of the tutorial from the OpenFOAM-v2006 version ?

Regards,

Yann