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

buoyantboussinesqsimplefoam solver for laminar flow

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

Like Tree2Likes
  • 1 Post By sixwp
  • 1 Post By nimasam

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 10, 2011, 15:49
Exclamation buoyantboussinesqsimplefoam solver for laminar flow
  #1
New Member
 
alireza golzari
Join Date: Jun 2011
Posts: 3
Rep Power: 14
alireza.glz is on a distinguished road
Hi Dear foamers
I'm a new user of openfoam and I'm trying to simulate a pipe through which a fluid passes in laminar flow and is being heated.
for this problem I'm using buoyantboussinesqsimplefoam solver and my RASProperties file is as follows:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

RASModel laminar;

turbulence off;

printCoeffs on;
/////////////////////////////////////////////////

but openfoam does not run my case and following errors appear:

Create time

Create mesh for time = 0


Reading g
Reading thermophysical properties

Reading field T

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

Selecting RAS turbulence model laminar
Calculating field beta*(g.h)


Starting time loop

Time = 1



LHS and RHS of - have different dimensions
dimensions : [0 3 -2 0 0 0 0] - [1 0 -2 0 0 0 0]
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/alireza/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/alireza/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam:perator-(Foam::dimensionSet const&, Foam::dimensionSet const&) in "/home/alireza/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#3 Foam::tmp<Foam::GeometricField<Foam::typeOfSum<dou ble, double>::type, Foam::fvsPatchField, Foam::surfaceMesh> > Foam:perator-<double, double, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::tmp<Foam::GeometricField< double, Foam::fvsPatchField, Foam::surfaceMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > const&) in "/home/alireza/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantBoussinesqSimpleFoam"
#4 main in "/home/alireza/OpenFOAM/OpenFOAM-1.6/applications/bin/linuxGccDPOpt/buoyantBoussinesqSimpleFoam"
#5 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#6 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/i386/elf/start.S:122


From function operator-(const dimensionSet& ds1, const dimensionSet& ds2)
in file dimensionSet/dimensionSet.C at line 423.

FOAM aborting

Aborted
/////////////////////////////////////////////////////

unfortunately I don't understand the source of these errors.
could anybody here help me to correct my case?
Please give me your valuable advises, please!
alireza.glz is offline   Reply With Quote

Old   July 11, 2011, 05:39
Default
  #2
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 15
sixwp is on a distinguished road
Hi Alireza,

the error tells you that you have wrongly defined the boundary unit.

if you compare:
[0 3 -2 0 0 0 0] with [1 0 -2 0 0 0 0]
m3.s-2 with kg.s-2 => therefore a difference in kg/m3

My guess is:
I bet your p/p_rgh initial value are define with [1 -1 -2 0 0 0 0]
Change it to: [0 2 -2 0 0 0 0] and you'll be fine.

The reason: depending the solver, density can be included or not

Best regards,
DanAndrea87 likes this.
__________________
Sam
sixwp is offline   Reply With Quote

Old   July 12, 2011, 02:22
Default
  #3
New Member
 
alireza golzari
Join Date: Jun 2011
Posts: 3
Rep Power: 14
alireza.glz is on a distinguished road
Dear samuel
Thanks a lot
You are right
I made the changes and my case run
Thanks a milion for your help
alireza.glz is offline   Reply With Quote

Old   July 12, 2011, 03:57
Default
  #4
Member
 
Samuel ARNAUD
Join Date: Feb 2011
Location: Grenoble, FRANCE
Posts: 39
Rep Power: 15
sixwp is on a distinguished road
you're most welcome
__________________
Sam
sixwp is offline   Reply With Quote

Old   May 2, 2012, 08:03
Default
  #5
Member
 
anonymous
Join Date: Mar 2012
Posts: 45
Rep Power: 14
libia87 is on a distinguished road
Hi!

I have a similar simulation, and when I run the problem, an error ocurrs:

--> FOAM FATAL ERROR:
Invalid wall function specification
Patch type for patch top must be wall
Current patch type is patch



From function kappatJayatillekeWallFunctionFvPatchScalarField::c heckType()
in file derivedFvPatchFields/wallFunctions/kappatWallFunctions/kappatJayatillekeWallFunction/kappatJayatillekeWallFunctionFvPatchScalarField.C at line 56.

FOAM aborting


I change the geometry and everythings it's ok, but the error says that the surface "top" it's not a wall, and I don't know what I have to do,

some ideas??
libia87 is offline   Reply With Quote

Old   May 2, 2012, 09:06
Default
  #6
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
nimasam is on a distinguished road
Quote:
Invalid wall function specification
Patch type for patch top must be wall
Current patch type is patch
patch type is patch so change it into wall, for changing the patch type you need edit your blockMeshDict or if you import a file from gambit change your boundary file
look at constant/polyMesh

P.S. if you read the error carefully, the solution is in it !!!
libia87 likes this.
nimasam is offline   Reply With Quote

Old   May 3, 2012, 04:22
Default
  #7
Member
 
anonymous
Join Date: Mar 2012
Posts: 45
Rep Power: 14
libia87 is on a distinguished road
Quote:
Originally Posted by nimasam View Post
patch type is patch so change it into wall, for changing the patch type you need edit your blockMeshDict or if you import a file from gambit change your boundary file
look at constant/polyMesh

P.S. if you read the error carefully, the solution is in it !!!
Ok, thanks!! I'm new with openfoam and it's difficult to me understand all the warnings that the program gives me. Thank you again! my simulation runs!

Another question, If I run my simulation and then I change some parameters and I run it again, paraFoam gives me the same results as the first run. How can I restart the simulation to obtain the new results? thanks!! (I change the writeInterval, but in my directory appears the same folders as before).
libia87 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Solver for an incompressible, turbulent flow with heat transfer tH3f0rC3 OpenFOAM Running, Solving & CFD 9 June 17, 2019 07:12
How do I select solver options for external flow over an aircraft by fluent? hadieliasi FLUENT 5 May 2, 2011 04:54
pre-conditioning for low mach number compressible flow solver Shenren_CN Main CFD Forum 0 April 29, 2011 22:07
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 22:31


All times are GMT -4. The time now is 07:43.