CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

BuoyantFoam problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 4, 2007, 04:34
Default Hi all, I have made a geome
  #1
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
Hi all,

I have made a geometry that is basically a cube with an inlet of fresh air (298K) and an outlet of wasted hot air, walls at 303 K, ground at 298K and ceiling at 308K. I fix the inlet velocity as 1 m/s and gravity in environmental properties. All I want is to simulate the profile for the air between inlet and outlet.

Here there are 2 pictures of the case:



When I run the case I obtain:

#0 Foam::error::printStack(Foam:stream&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::compressible::turbulenceModels::kEpsilon::kE psilon(Foam::GeometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libcompressibleTurbulenceModels.so"
#4 Foam::compressible::turbulenceModel::adddictionary ConstructorToTable<foam::compressible::turbulencem odels::kepsilon>::New(Foam::GeometricField<double, > const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libcompressibleTurbulenceModels.so"
#5 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libcompressibleTurbulenceModels.so"
#6 main in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantFoam"
#7 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantFoam"

Is there anything that I forget?

Thank you!!
Miquel
miquelvinya is offline   Reply With Quote

Old   October 4, 2007, 05:31
Default I believe there is something w
  #2
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
I believe there is something wrong with the internal field that you have initialised in k epsilon file... Make sure that they are not initialised to zero
prashant24983 is offline   Reply With Quote

Old   October 4, 2007, 07:24
Default Thanks Prashant, I tryed yo
  #3
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
Thanks Prashant,

I tryed your option because I forgot k and epsilon values for the internal field, but the problem persists with the same message...

#0 Foam::error::printStack(Foam:stream&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::compressible::turbulenceModels::kEpsilon::kE psilon(Foam::GeometricField<do uble,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libcompressibleTurbulenc eModels.so"
#4 Foam::compressible::turbulenceModel::adddictionary ConstructorToTable<foam::compr essible::turbulencemodels::kepsilon>::New(Foam::Ge ometricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libcompressibleTurbulenc eModels.so"
#5 Foam::compressible::turbulenceModel::New(Foam::Geo metricField<double,> const&, Foam::GeometricField<foam::vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double,> const&, Foam::basicThermo&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libcompressibleTurbulenc eModels.so"
#6 main in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantFoam "
#7 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantFoam "

Wha can I try now? Thanks in advance,
Miquel
miquelvinya is offline   Reply With Quote

Old   October 4, 2007, 08:25
Default Check the R file in the 0 dire
  #4
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
Check the R file in the 0 directory, just like in the buoyantFoam tutorial set the internalField likewise

dimensions [0 2 -2 0 0 0 0];
internalField uniform (0 0 0 0 0 0 0 0 0)

followed by the patches, I believe the zeroGradient patch should work.

If the problem persists I believe the problem is with the implementation of boundary conditions.

Can you post the case or mail it at prashantojha@gmail.com
prashant24983 is offline   Reply With Quote

Old   October 5, 2007, 05:13
Default I would specially thank to Pra
  #5
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
I would specially thank to Prashant for helping me.

I changed intitial conditions such as k, epsilon and pressure. With this there is no zerodivision problem. And then I reduced timestep to 0.05 s.

Here there are the results:


=)
Miquel
miquelvinya is offline   Reply With Quote

Old   October 5, 2007, 10:51
Default The next step was to create a
  #6
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
The next step was to create a more complex case based in the above described. When I try to run the case the solver shows me this message:

Create time

Create mesh for time = 0


Reading environmentalProperties
Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
#0 Foam::error::printStack(Foam:stream&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xffffe420]
#3 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::calculate() in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMo dels.so"
#4 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::hThermo(Foam::fvMesh const&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMo dels.so"
#5 Foam::basicThermo::addfvMeshConstructorToTable<foa m::hthermo<foam::puremixture<f oam::consttransport<foam::speciethermo<foam::hcons tthermo<foam::perfectgas> > > > > >::New(Foam::fvMesh const&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMo dels.so"
#6 Foam::basicThermo::New(Foam::fvMesh const&) in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMo dels.so"
#7 main in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantFoam "
#8 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/mvinya/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/buoyantFoam "

What can I do now? Thank you!!
Miquel
miquelvinya is offline   Reply With Quote

Old   October 6, 2007, 01:18
Default are you trying to add aother s
  #7
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
are you trying to add aother specie in the previous case? What complexities are you talking about?
prashant24983 is offline   Reply With Quote

Old   October 8, 2007, 03:15
Default Hello! The geometry is more
  #8
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
Hello!

The geometry is more complex and there are surfaces up to T = 473 K. Is the buoyantFoam thermodynamic package prepared for these temperatures? Must I change the solver?

Thank you!
Miquel
miquelvinya is offline   Reply With Quote

Old   October 12, 2007, 08:37
Default Excuse me. I want to know w
  #9
New Member
 
Giro
Join Date: Mar 2009
Posts: 13
Rep Power: 8
girogirozakk is on a distinguished road
Excuse me.

I want to know what type of boundary aer there?
Now , I try same-type calculation.
But , I cannot understand what type of boundary "Inlet Hole" & "Outlet Hole".

I tried "freestreem" or "atmosphere" , but cannot calc.
So , I cannot understand to set U,p ,etc..

Please teach me the setting.


Sorry,I'm not good English-speaker , my sentence is Hard to read...




Thank you!
Giro.
girogirozakk is offline   Reply With Quote

Old   October 15, 2007, 09:15
Default Hello, The values used for
  #10
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
Hello,

The values used for my first buoyantFoam simulation were:

INLET

- Type: inlet
- U_y = 0.1 m/s
- k = 0.1
- epsilon = 0.01
- P = 1.0E5 Pa
- T = 295 K

OUTLET

- Type: outlet
- U = 0.0 m/s
- k = 0.1
- epsilon = 0.01
- P = 1.0E5 Pa

WALLS

- Type: fixedTemperatureWallFunction
- U = 0.0 m/s
- k = 0.1
- epsilon = 0.01
- P = 1.0E5 Pa
- T = 303 K

INTERNAL FIELD

- U = 0.0 m/s
- k = 0.1
- epsilon = 0.01
- P = 1.0E5 Pa
- T = 298 K

My simulaion has 90.000 tetrahedra and I fixed dt=0.001 s.

My new simulation crashes. Probably is because wall temperaturea are above 373K.

Is air thermodynamic package prepared for properties at that temperature (hThermo...)??

If not, why my crash report shows me a sigFpe message??



Thanks for your attention,
Miquel
miquelvinya is offline   Reply With Quote

Old   October 15, 2007, 09:28
Default Please check your thermoPhysic
  #11
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
Please check your thermoPhysicalProperties file.

Try this for Air instead of previous form... Hope this helps.

mixture air 1 28.9 1000 2.544e+06 1.8e-05 0.7;
prashant24983 is offline   Reply With Quote

Old   October 15, 2007, 22:18
Default Thanks ! I could get some res
  #12
New Member
 
Giro
Join Date: Mar 2009
Posts: 13
Rep Power: 8
girogirozakk is on a distinguished road
Thanks !
I could get some result for my calc.

Giro
girogirozakk is offline   Reply With Quote

Old   October 16, 2007, 02:28
Default Thanks Prashant, The case i
  #13
New Member
 
Miquel Vinya
Join Date: Mar 2009
Location: Tarragona
Posts: 13
Rep Power: 8
miquelvinya is on a distinguished road
Thanks Prashant,

The case is correctly running.

Regards,
Miquel
miquelvinya is offline   Reply With Quote

Old   November 6, 2008, 07:17
Default Hi, i'm trying to use the s
  #14
ep4
Member
 
Pattyn Eric
Join Date: Mar 2009
Posts: 61
Rep Power: 8
ep4 is on a distinguished road
Hi,

i'm trying to use the solver buoyantFoam for natural convection and having some difficulties: temperatures increasing untill the simulation crashes.

I would like to make a simple test. I have a flat plate case wich gives satisfying results with turbFoam. I would like to reproduce these results with buoyantFoam, setting g= (0 0 0) in the envireonnemental properties (no buoyancy effects) and using constant temperature boundaries, equal to the initial intern field (no heat transfer phenomenons). However, it doesn't work. I have also tried with adiabatic boundaries. For each case, i obtain after several time steps temperature of order of 10000000 K !!

It makes me think to stability problems. Using the Euler or backwrd time schemes which are implicit (User Guide U-112), i don't understand, with my theoritical knowledge of CFD, why i encounter problems of stability.
I have similar problems with buoyantSimpleFoam.

Could someone say me what i'm doing wrong?

Thank you

Eric
ep4 is offline   Reply With Quote

Old   November 6, 2008, 07:49
Default Use "fixedFluxBuoyantPressure"
  #15
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
Use "fixedFluxBuoyantPressure" for the flat plate boundary at "p". I think the problem is in getting the boundary conditions right. Can you paste your "p" boundary condition here?
prashant24983 is offline   Reply With Quote

Old   November 7, 2008, 04:07
Default Hi everybody, i have change
  #16
ep4
Member
 
Pattyn Eric
Join Date: Mar 2009
Posts: 61
Rep Power: 8
ep4 is on a distinguished road
Hi everybody,

i have changed the intern initial velocity field by (10 0 0) in place of (0 0 0). It is more stable, i mean, the really high temperatures appear later...

Any idea to help me?

Thank you
ep4 is offline   Reply With Quote

Old   November 7, 2008, 04:58
Default Hello Pattyn, I have reali
  #17
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
Hello Pattyn,

I have realized that the code for buoyantFoam and buoyantSimpleFoam has changed considerably from version 1.4.1 to 1.5. Even though, I am not very clear about the declaration of pd and p in version 1.5, I can help you with some basics based on how I understand the buoyantFoam from version 1.4.1.

pd is dynamic pressure. In the simple form is 0.5*rho*V^2. In your declaration of pd you provide a '0' value to fluid coming in and going out. Which I must say is not right. For inlet and outlet provide a zeroGradient boundary condition.

Also make sure that you are not providing a fixedValue boundary condition for inlet and outlet patches in 'T'

Most times your solution may blow up because of the courant number.Check the courant number as your iterations progresses.

I have few questions, What is Up? is it wall? if its not then your pd should look something like this.

internalField uniform 0;

boundaryField
{
WallDown
{
type fixedFluxBuoyantPressure;
value uniform 0;
}

Inlet
{
type zeroGradient;
}

Outlet
{
type zeroGradient;
}

Up
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
prashant24983 is offline   Reply With Quote

Old   November 7, 2008, 05:12
Default If your boundary conditions ar
  #18
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
If your boundary conditions are not right, this leads to wrong calculation of the pEqn which eventually results in irrational p values, and irrational velocities.. so your maximum courant number eventually will make you worry!
prashant24983 is offline   Reply With Quote

Old   November 7, 2008, 08:33
Default Hello Eric, Its really dif
  #19
Member
 
Prashant Ojha
Join Date: Mar 2009
Posts: 38
Rep Power: 8
prashant24983 is on a distinguished road
Hello Eric,

Its really difficult to say what might be the exact trouble!

If you are using turbulence do you observe bounding errors? Just switch off turbulence and see if the error remains.

If possible send your case to me prashantojha@gmail.com

Though, I can look at it only tomorrow.

Regards,
prashant24983 is offline   Reply With Quote

Old   November 7, 2008, 09:03
Default Thank you for your help. I'm g
  #20
ep4
Member
 
Pattyn Eric
Join Date: Mar 2009
Posts: 61
Rep Power: 8
ep4 is on a distinguished road
Thank you for your help. I'm going to search a solution this afternoon and if it doesn't work i will send the case to you.

I have the impression that the problem could come from the T file given that the case runs well with turbFoam.

It appears me logical tu use a fixed temperature at the inlet and the plate, a zeroGradient condition at the boundary Up (infront of the plate).
And i'm wondering if it would be possible no to impose anything at the Outlet. It is not logical for me to impose a temperature or a gradient at this boundary.

Any comments?

Thank you
ep4 is offline   Reply With Quote

Reply

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
BuoyantFoam and total pressure bc problem is a bug ariorus OpenFOAM Bugs 3 January 28, 2008 05:20
TotalPressure and buoyantFoam ariorus OpenFOAM Running, Solving & CFD 1 January 22, 2008 09:41
ThermoPhysicalProperties in buoyantFoam prashant24983 OpenFOAM Running, Solving & CFD 0 October 6, 2007 09:40
BuoyantFoam braennstroem OpenFOAM Running, Solving & CFD 22 September 19, 2007 16:55
Release 13 buoyantFoam braennstroem OpenFOAM 0 March 30, 2006 02:43


All times are GMT -4. The time now is 04:50.