CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

groovyBC

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree2Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   June 17, 2014, 23:49
Default Wave Tank Simulation
  #21
Member
 
Gautami Erukulla
Join Date: Mar 2009
Posts: 43
Rep Power: 9
gautami is on a distinguished road
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.
gautami is offline   Reply With Quote

Old   June 18, 2014, 09:53
Default
  #22
Member
 
zohreh imani nejad
Join Date: Feb 2014
Location: gonabad
Posts: 55
Rep Power: 4
imani is on a distinguished road
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
imani is offline   Reply With Quote

Old   June 18, 2014, 10:03
Default
  #23
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,189
Blog Entries: 1
Rep Power: 15
nimasam is on a distinguished road
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:
Quote:
( true/false value ? true : false )
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog (http://openfoam.blogfa.com/)
Training Course on OpenFOAM at (http://www.isme.ir/)
nimasam is offline   Reply With Quote

Old   June 18, 2014, 10:07
Default
  #24
Member
 
zohreh imani nejad
Join Date: Feb 2014
Location: gonabad
Posts: 55
Rep Power: 4
imani is on a distinguished road
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"
imani is offline   Reply With Quote

Old   June 18, 2014, 12:09
Default
  #25
Member
 
Gautami Erukulla
Join Date: Mar 2009
Posts: 43
Rep Power: 9
gautami is on a distinguished road
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.
gautami is offline   Reply With Quote

Old   June 18, 2014, 12:17
Default
  #26
Member
 
zohreh imani nejad
Join Date: Feb 2014
Location: gonabad
Posts: 55
Rep Power: 4
imani is on a distinguished road
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?
imani is offline   Reply With Quote

Old   June 18, 2014, 13:09
Default
  #27
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,189
Blog Entries: 1
Rep Power: 15
nimasam is on a distinguished road
well you should write something like this:
Quote:
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)",
HesamAUT likes this.
__________________
Telegram channel (https://telegram.me/openfoam4Iranian)
My Weblog (http://openfoam.blogfa.com/)
Training Course on OpenFOAM at (http://www.isme.ir/)
nimasam is offline   Reply With Quote

Old   June 18, 2014, 13:24
Default
  #28
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,926
Rep Power: 41
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by gautami View Post
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.
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   June 18, 2014, 13:44
Default
  #29
Member
 
Gautami Erukulla
Join Date: Mar 2009
Posts: 43
Rep Power: 9
gautami is on a distinguished road
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
gautami is offline   Reply With Quote

Old   June 18, 2014, 15:30
Default
  #30
New Member
 
Geir Karlsen
Join Date: Nov 2013
Location: Norway
Posts: 14
Rep Power: 4
gkarlsen is on a distinguished road
Quote:
Originally Posted by gautami View Post
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)

groovyWaveTank
gkarlsen is offline   Reply With Quote

Old   June 18, 2014, 16:27
Default
  #31
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,926
Rep Power: 41
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by gautami View Post
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
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   June 18, 2014, 22:35
Default
  #32
Member
 
Gautami Erukulla
Join Date: Mar 2009
Posts: 43
Rep Power: 9
gautami is on a distinguished road
Dear Gschaider,

I totally understand.

Thank you.

Regards,
Gautami.
gautami is offline   Reply With Quote

Old   June 19, 2014, 00:54
Default
  #33
Member
 
zohreh imani nejad
Join Date: Feb 2014
Location: gonabad
Posts: 55
Rep Power: 4
imani is on a distinguished road
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

Last edited by imani; June 19, 2014 at 04:12.
imani is offline   Reply With Quote

Old   June 19, 2014, 11:45
Default
  #34
Member
 
Gautami Erukulla
Join Date: Mar 2009
Posts: 43
Rep Power: 9
gautami is on a distinguished road
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.
gkarlsen likes this.
gautami is offline   Reply With Quote

Old   June 19, 2014, 14:50
Default
  #35
Member
 
Gautami Erukulla
Join Date: Mar 2009
Posts: 43
Rep Power: 9
gautami is on a distinguished road
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
gautami is offline   Reply With Quote

Old   July 2, 2014, 14:12
Default Boundary condition implementation
  #36
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
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
nandiganavishal is offline   Reply With Quote

Old   July 2, 2014, 18:43
Default
  #37
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,926
Rep Power: 41
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by nandiganavishal View Post
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
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   July 2, 2014, 21:39
Default
  #38
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
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
nandiganavishal is offline   Reply With Quote

Old   July 3, 2014, 15:11
Default
  #39
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,926
Rep Power: 41
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by nandiganavishal View Post
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 DirectionMixed mixed bc for a description of what directionMixed/mixed do. In the groovy-variant all described values can be set by expressions
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   July 7, 2014, 17:17
Default
  #40
New Member
 
Geir Karlsen
Join Date: Nov 2013
Location: Norway
Posts: 14
Rep Power: 4
gkarlsen is on a distinguished road
Quote:
Originally Posted by gautami View Post
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
Attached Images
File Type: png SurfaceElevation.png (21.8 KB, 10 views)
gkarlsen is offline   Reply With Quote

Reply

Tags
groovybc, swak4foam

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
groovyBC issue - k and epsilon sagnikmazumdar OpenFOAM Running, Solving & CFD 24 March 1, 2015 08:16
groovyBC and funkySetFields married and got a kid named swak4Foam gschaider OpenFOAM 164 January 13, 2015 03:52
groovyBC for oscillatory flow liybzd OpenFOAM 1 December 15, 2013 08:08
switch between existing BCs with groovyBC?! michielm OpenFOAM Pre-Processing 3 September 4, 2013 08:01
groovyBC and Eqn.setReference() benk OpenFOAM 3 June 2, 2011 08:49


All times are GMT -4. The time now is 23:00.