[gautami]

Dear All,

I am using OpenFOAM version 2.2.0. on SuSE linux 12.2.
I have installed swak4Foam (Version for OpenFOAM 2.x).I am trying to simulate wave tank using groovyBC with solver interFOAM. Kindly can anyone let me know where I could find an example case for wave tank simulation.
Thank you.

Regards Gautami.

[imani]

hi
i want to modelling heat transfer in a hot room by groovyBC
in one of my room's wall there is windows
in this wall i have 4 regions with different temperature
how can i do this?
can i have different "fractionExpression" and "valueExpression"??????????
thanks

[nimasam]

you can define your different Region with its position in value expression and give them different Value or gradient

Also look at following syntax in groovyBC command:

( true/false value ? true : false )

[imani]

thanks alot
my mean is like bellow....how can i do this?

valueExpression "150"
valueExpression "45"
valueExpression "23"
valueExpression "36"
.
.
.
fractionExpression 1 "position1..... ? 1:0"
fractionExpression 2 "position1..... ? 2:0"
fractionExpression 1 "position1..... ? 3:0"
fractionExpression 1 "position1..... ? 4:0"

[gautami]

Dear All,

I am using OpenFOAM version 2.2.0. on SuSE linux 12.2. I have installed swak4Foam (Version for OpenFOAM 2.x).

I am trying to simulate wave tank using groovyBC with solver interFOAM. I have the simulation case of wave tank which I did earlier in version OpenFOAM 1-5. Now I have used the same case file and made changes for pd to p_rgh and gamma to alpha1.

I run the case and then I get the following error.


#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib64/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 void Foam::MULES::explicitSolve<Foam::geometricOneField , Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7
in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/interFoam"
#8 __libc_start_main in "/lib64/libc.so.6"
#9
at /home/abuild/rpmbuild/BUILD/glibc-2.15/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

Kindly can somebody help me with this issue.Thank you

Regards,
Gautami.

[imani]

valueExpression " pos().y < 0.18 ? 3000 : 0)","((pos().y > 2.2 && pos().y < 2.9) && ((pos().x > 4.498 && pos().x < 7.502) || (pos().x > 0.498 && pos().x < 3.502)) ? 1000 : 0)" ;






it is not true
how can i correct it?

[nimasam]

well you should write something like this:

valueExpression "pos().y < 0.18 ? 3000 : ((pos().y > 2.2 && pos().y < 2.9) && ((pos().x > 4.498 && pos().x < 7.502) || (pos().x > 0.498 && pos().x < 3.502)) ? 1000 : 0)",

[gschaider]

Originally Posted by gautami

Dear All,

I am using OpenFOAM version 2.2.0. on SuSE linux 12.2. I have installed swak4Foam (Version for OpenFOAM 2.x).

I am trying to simulate wave tank using groovyBC with solver interFOAM. I have the simulation case of wave tank which I did earlier in version OpenFOAM 1-5. Now I have used the same case file and made changes for pd to p_rgh and gamma to alpha1.

I run the case and then I get the following error.


#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib64/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 void Foam::MULES::limit<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double, int, bool) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 void Foam::MULES::explicitSolve<Foam::geometricOneField , Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7
in "/opt/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/interFoam"
#8 __libc_start_main in "/lib64/libc.so.6"
#9
at /home/abuild/rpmbuild/BUILD/glibc-2.15/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

Kindly can somebody help me with this issue.Thank you

Regards,
Gautami.

That doesn't look like a groovyBC-problem per se. Most likely what is happening is that you're specifying a boundary condition ("made changes") that is unphysical (for instance: violates continuity) and because of this the solution blows up.

[gautami]

Dear Gschaider,

How are you doing? I am happy to see a reply from you.

I am really not able to make out what change is causing the run to blow. I am really sorry for the trouble, kindly could you please have a look at my case files and guide.I am extremely sorry for the trouble,but as always your suggestions are invaluable.Thank you.

Please find attached folders:
0 constant and system of my run.

Regards,
Gautami.0.tar.gz

constant.tar.gz

system.tar.gz

[gkarlsen]

Originally Posted by gautami

Dear All,

I am using OpenFOAM version 2.2.0. on SuSE linux 12.2.
I have installed swak4Foam (Version for OpenFOAM 2.x).I am trying to simulate wave tank using groovyBC with solver interFOAM. Kindly can anyone let me know where I could find an example case for wave tank simulation.
Thank you.

Regards Gautami.

Maybe this post could be of use to you (runs with OF 2.2)

http://www.cfd-online.com/Forums/ope...ywavetank.html

[gschaider]

Originally Posted by gautami

Dear Gschaider,

How are you doing? I am happy to see a reply from you.

I am really not able to make out what change is causing the run to blow. I am really sorry for the trouble, kindly could you please have a look at my case files and guide.I am extremely sorry for the trouble,but as always your suggestions are invaluable.Thank you.

Please find attached folders:
0 constant and system of my run.

Regards,
Gautami.Attachment 31738

Attachment 31739

Attachment 31740

I won't have any time to look at other peoples cases in the next few weeks.

Sorry

[gautami]

Dear Gschaider,

I totally understand.

Thank you.

Regards,
Gautami.

[imani]

thank you very much nima
excuse me i have another questions
now i define my different regions with different temperatures in valueExpression but
1.what i difference of these two form ( my means is fractionExpression)


type groovyBC;

variables "Kwall=0.06;Kglazing=0.06;Qwall=-8.5;Qglazing=-40.8;";// taken from table 3

valueExpression"pos().y < 0.15 ? 37 : (pos().y > 0.198 && pos().y < 0.2) && ((pos().x > 4.498 && pos().x < 7.502) || (pos().x > 0.498 && pos().x < 3.502)) ? 18 : (pos().y > 0.2 && pos().y < 2.2) && ((pos().x > 0.498 && pos().x < 0.5) || (pos().x > 3.5 && pos().x < 3.502) || (pos().x > 4.498 && pos().x < 4.5) || (pos().x > 7.5 && pos().x < 7.502)) ? 18 : (pos().y > 2.2 && pos().y < 2.202) && ((pos().x > 4.498 && pos().x < 7.502) || (pos().x > 0.498 && pos().x < 3.502)) ? 18 0";

gradientExpression "(pos().y > 0.2 && pos().y < 2.2) && ((pos().x > 0.5 && pos().x < 3.5) || (pos().x > 4.5 && pos().x < 7.5)) ? Qglazing/Kglazing : Qwall/Kwall";// in external wall which containts windows, BC is constant negative heat flux

fractionExpression "1";




or






type groovyBC;

variables "Kwall=0.06;Kglazing=0.06;Qwall=-8.5;Qglazing=-40.8;";// taken from table 3

valueExpression"pos().y < 0.15 ? 37 : (pos().y > 0.198 && pos().y < 0.2) && ((pos().x > 4.498 && pos().x < 7.502) || (pos().x > 0.498 && pos().x < 3.502)) ? 18 : (pos().y > 0.2 && pos().y < 2.2) && ((pos().x > 0.498 && pos().x < 0.5) || (pos().x > 3.5 && pos().x < 3.502) || (pos().x > 4.498 && pos().x < 4.5) || (pos().x > 7.5 && pos().x < 7.502)) ? 18 : (pos().y > 2.2 && pos().y < 2.202) && ((pos().x > 4.498 && pos().x < 7.502) || (pos().x > 0.498 && pos().x < 3.502)) ? 18 0";

gradientExpression "(pos().y > 0.2 && pos().y < 2.2) && ((pos().x > 0.5 && pos().x < 3.5) || (pos().x > 4.5 && pos().x < 7.5)) ? Qglazing/Kglazing : Qwall/Kwall";// in external wall which containts windows, BC is constant negative heat flux

fractionExpression "0";



2. what is this error...


Time = 1

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00995599, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.0104898, No Iterations 1
DILUPBiCG: Solving for T, Initial residual = 1, Final residual = 0.0591834, No Iterations 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 7.1898, No Iterations 1001
time step continuity errors : sum local = 84.0209, global = 0.0235848, cumulative = 0.0235848
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 at kEpsilon.C:0
#6 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/buoyantBoussinesqSimpleFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/home/dell/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/buoyantBoussinesqSimpleFoam"
Floating point exception (core dumped)



i think my problem is that my mesh is very big....because my room is very large ,when i fine my meshs it can not make blockMesh successfully
what should i do
thanks alot

[gautami]

Dear Gkarlsen,

Thank you very much. I appreciate you providing the link for the wave tank case files.I was able to download the case files and run the simulation. It did not blow . Right now I am working on the post processing. Shall keep you posted. Thank you.

Regards,
Gautami.

[gautami]

Dear All,

I was able to run the wavetank simulation with groovyBC at the inlet(2nd order Stokes waves).Amplitude=0.1 and wavelength=5.The simulation does not blow but the results do not seem to be good.

The wave tank length is in x direction ranging from (0 to 30), height is in the z direction ranging from (-1.0 to 1.0) with still water at 0,and the gravity vector is defined as (0 0 -9.81).

Case-1
I got the free surface elevation plot near the inlet, shown in the image file (Free Surface Elevation Case-1) (x axis-time and y axiz-surface elevation).It does not seem to be correct.At time=0 the elevation is not zero. Plot of waves at time 15s is shown in image file (waves at time 15s case-1)

Case-2
I just changed the coordinates of the wave tank in z direction range to (0.0 to 2.0) with still water at 1.0.The free surface elevation plot near the inlet is not correct shown in the image file (Free Surface Elevation Case-1).Plot of waves at time 15s is shown in image file (waves at time 15s case-1)

In both the cases the results look wrong.Just changing the coordinates changed the results drastically.Kindly can somebody please advise me.Thank you.

Regards,
Gautami.

Free Surface Elevation Case-1.png

Free Surface Elevation Case-2.png

waves at time 15s case-1.jpg

waves at time 15s case-2.jpg

[nandiganavishal]

Hello All,

I am trying to solve a Laminar flow problem inside an axisymmetric cylindrical geometry with slip, where the slip length is specified by the user.

For this, I am trying to implement, the following boundary condition for velocity at the wall of a cylindrical geometry

dUz/dn = Uz/Ls; Ur = 0 everywhere in the domain

where Ls: Slip velocity
n = normal direction (radial direction in this case)
Uz = axial velocity
Ur = radial velocity


I would like to know if this BC can be implemented using GroovyBC/swak4Foam. Please let me know.

Thanks

Vishal

[gschaider]

Originally Posted by nandiganavishal

Hello All,

I am trying to solve a Laminar flow problem inside an axisymmetric cylindrical geometry with slip, where the slip length is specified by the user.

For this, I am trying to implement, the following boundary condition for velocity at the wall of a cylindrical geometry

dUz/dn = Uz/Ls; Ur = 0 everywhere in the domain

where Ls: Slip velocity
n = normal direction (radial direction in this case)
Uz = axial velocity
Ur = radial velocity


I would like to know if this BC can be implemented using GroovyBC/swak4Foam. Please let me know.

Thanks

Vishal

Not directly I'm afraid

[nandiganavishal]

Thanks, Bernhard.

Could you direct me how to implement this robin boundary condition.

There is a small correction in my previous post, Ls represents, user specified slip length and not slip velocity.

Thanks

Vishal

[gschaider]

Originally Posted by nandiganavishal

Thanks, Bernhard.

Could you direct me how to implement this robin boundary condition.

There is a small correction in my previous post, Ls represents, user specified slip length and not slip velocity.

Thanks

Vishal

Robin is called "mixed" in OF. groovyBC does that out of the box (see the Docu)

The problem with what you want is that you want Dirichlet in one direction and Neumann in the other. So probably directtionMixed might be what you want. There is even a groovyBC-Variant for that (groovyMixedBC I think. Haven't used it for some time)

See http://www.cfd-online.com/Forums/ope...-mixed-bc.html for a description of what directionMixed/mixed do. In the groovy-variant all described values can be set by expressions

[gkarlsen]

Originally Posted by gautami

Dear All,

I was able to run the wavetank simulation with groovyBC at the inlet(2nd order Stokes waves).Amplitude=0.1 and wavelength=5.The simulation does not blow but the results do not seem to be good.

The wave tank length is in x direction ranging from (0 to 30), height is in the z direction ranging from (-1.0 to 1.0) with still water at 0,and the gravity vector is defined as (0 0 -9.81).

Case-1
I got the free surface elevation plot near the inlet, shown in the image file (Free Surface Elevation Case-1) (x axis-time and y axiz-surface elevation).It does not seem to be correct.At time=0 the elevation is not zero. Plot of waves at time 15s is shown in image file (waves at time 15s case-1)

Are you sure you are plotting the right parameter? Your two last plots seem to make sense. I also plotted the surface elevation at an arbitrary point, and from what I can see it fluctuates around Z=0 (see attached).

Regards,
Geir