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

SteadyState solver for compressible flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 15, 2008, 09:38
Default Hi to all, I have a flow with
  #1
matteo_gautero
Guest
 
Posts: n/a
Hi to all,
I have a flow with Mach=0.3. I'm needing a steady state solver for compressible flow, but in the user guide there aren't any solver for steady state case. Can anyone help me?

Matteo.
  Reply With Quote

Old   February 15, 2008, 18:09
Default check out rhoSimpleFoam. th
  #2
Member
 
Shaun Darmody
Join Date: Mar 2009
Posts: 36
Rep Power: 8
shaun is on a distinguished road
check out rhoSimpleFoam.

this will be what you are looking for.

Regards

Shaun.D
shaun is offline   Reply With Quote

Old   February 18, 2008, 10:00
Default Thanks Shaun. I have a problem
  #3
matteo_gautero
Guest
 
Posts: n/a
Thanks Shaun. I have a problem. This solver is not present in FoamX and in the tutorials folder. I have to compile rhoSimpleFoam or I've done some errors in the installation procedure?

Regards,
Matteo.
  Reply With Quote

Old   February 19, 2008, 04:24
Default You just have to set up the ca
  #4
New Member
 
Maximilian Graser
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 18
Rep Power: 8
graser is on a distinguished road
You just have to set up the case by hand without FoamX.
graser is offline   Reply With Quote

Old   February 19, 2008, 04:37
Default Hi, I also want to know how
  #5
myting
Guest
 
Posts: n/a
Hi,

I also want to know how I can set up the case by hand. Can somebody give me an example for this solver please.
  Reply With Quote

Old   February 19, 2008, 04:56
Default look at these tutorials: rhoE
  #6
New Member
 
Maximilian Graser
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 18
Rep Power: 8
graser is on a distinguished road
look at these tutorials:
rhoExplicitPorousSimpleFoam
rhoImplicitPorousSimpleFoam
they use rhoSimpleFoam solver
graser is offline   Reply With Quote

Old   February 19, 2008, 07:37
Default you could also simple search t
  #7
Member
 
Shaun Darmody
Join Date: Mar 2009
Posts: 36
Rep Power: 8
shaun is on a distinguished road
you could also simple search the forum for "rhoSimpleFoam" and you will find some references to case setup.

I believe I can post the case setup tomorrow when at work.

Shaun.D
shaun is offline   Reply With Quote

Old   February 19, 2008, 07:59
Default Hi, thanks, I think, I will
  #8
myting
Guest
 
Posts: n/a
Hi,

thanks, I think, I will try to change the case "rhoImplicitPorousSimpleFoam" solver to "rhoSimpleFoam".


In the forum I found this information:

"I have tried to create a FoamX configuration for rhoSimpleFoam. Just untar it under $FOAM_SOLVERS/rhoSimpleFoam which will give you the directory FoamX with all the files needed for your first steps with rhoSimpleFoam."

I also will try it out today. I hope it is working under version 1.4.1 :-)

Thanks a lot
  Reply With Quote

Old   February 26, 2008, 07:47
Default Hi to all, I set up a simple
  #9
matteo_gautero
Guest
 
Posts: n/a
Hi to all,
I set up a simple case to understand what rhoSimpleFoam. In the forum I found a discussion where they said that to set up the case I had work with rhoTurbFoam. So I set up the case with rhoTurbFoam and the I modified the fvSolution and fvSchemes dictionary with the corresponding file present in simpleFoam. When I run the case I receive this error message:

Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
#0 Foam::error::printStack(Foam:stream&) in "/home/gauta/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/gauta/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/gauta/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMod els.so"
#4 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::hThermo(Foam::fvMesh const&) in "/home/gauta/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMod els.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/gauta/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMod els.so"
#6 Foam::basicThermo::New(Foam::fvMesh const&) in "/home/gauta/OpenFOAM/OpenFOAM-1.4.1/lib/linuxGccDPOpt/libbasicThermophysicalMod els.so"
#7 main in "/home/gauta/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/rhoSimpleFoa m"
#8 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#9 Foam::regIOobject::readIfModified() in "/home/gauta/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linuxGccDPOpt/rhoSimpleFoa m"
Floating point exception (core dumped)

Can anyone help me?
Thanks
  Reply With Quote

Old   February 26, 2008, 08:02
Default Hi, can you give me your em
  #10
ttian
Guest
 
Posts: n/a
Hi,

can you give me your email? I can send you a file.
Just untar it under $FOAM_SOLVERS/compressible/rhoSimpleFoam which will give you the directory FoamX with all the files needed for your first steps with rhoSimpleFoam. I found it in this forum.

Bye
  Reply With Quote

Old   March 3, 2008, 12:14
Default Hi, can anybody send me the
  #11
New Member
 
Marco Zardetto
Join Date: Mar 2009
Location: Italy
Posts: 6
Rep Power: 8
sharmak is on a distinguished road
Hi,

can anybody send me the file that I have to untar in the $FOAM_SOLVERS/compressible/rhoSimpleFoam directory in order to use the rhoSimpleFoam solver in FoamX?
The file I found in the forum (rhoSimpleFoam_FoamX.tar.gz)is a .unk file... I don't know how to use it.
Or anybody can send me an example to compile correctly the files (fvSolution, fvSchemes...) with this solver?

Thanks
Bye
sharmak is offline   Reply With Quote

Old   March 3, 2008, 12:23
Default Hi Marco I think you want t
  #12
Senior Member
 
Cedric DUPRAT
Join Date: Mar 2009
Location: Belgium
Posts: 179
Rep Power: 8
cedric_duprat is on a distinguished road
Hi Marco

I think you want to do something like that:
1- mv rhoSimpleFoam_FoamX.unk rhoSimpleFoam_FoamX.tar.gz
2- gunzip rhoSimpleFoam_FoamX.tar.gz
3- tar -xvf rhoSimpleFoam_FoamX.tar
and, I think you should get what you want
:o)

Cedric
cedric_duprat is offline   Reply With Quote

Old   March 3, 2008, 12:23
Default Hi Marco, when you save the
  #13
Member
 
Francesco Boschetto
Join Date: Mar 2009
Location: Italy
Posts: 56
Rep Power: 8
francesco_b is on a distinguished road
Hi Marco,

when you save the file, simply rename it rhoSimpleFoam_FoamX.tar.gz and you'll download the correct file. Then extract it and start to play

Regards

Francesco
francesco_b is offline   Reply With Quote

Old   March 3, 2008, 13:51
Default Thanks Francesco and Cedric,
  #14
New Member
 
Marco Zardetto
Join Date: Mar 2009
Location: Italy
Posts: 6
Rep Power: 8
sharmak is on a distinguished road
Thanks Francesco and Cedric,
you are right. Now the solver is in FoamX but it gives error creating a new case! Probably there are some errors in the file.

I'm sorry but I'm new to OpenFoam...

Anyhow I've to study a compressible air flow in steady state. I started from rhoturbFoam solver and than I modified manually some files. When the simulation starts it plots this errors:

Exec : rhoSimpleFoam . v990_prova
Date : Mar 03 2008
Time : 18:22:51
Host : computer-marco
PID : 6983
Root : /home/simulation/OpenFOAM/simulation-1.4/run/tutorials/rhoTurbFoam
Case : v990_prova
Nprocs : 1
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hThermo<puremixture<consttransport<speciethermo<hc onstthermo<perfectgas>>>>>
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
#2 Uninterpreted: [0xffffe420]
#3 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::calculate()
#4 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::hThermo(Foam::fvMesh const&)
#5 Foam::basicThermo::addfvMeshConstructorToTable<foa m::hthermo<foam::puremixture<f oam::consttransport<foam::speciethermo<foam::hcons tthermo<foam::perfectgas> > > > > >::New(Foam::fvMesh const&)
#6 Foam::basicThermo::New(Foam::fvMesh const&)
#7 main
#8 __libc_start_main
#9 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122
Floating point exception (core dumped)

What does it mean? How can I find the errors?

Thanks
Bye
sharmak is offline   Reply With Quote

Old   March 4, 2008, 04:15
Default Hi, maybe you have to chang
  #15
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 103
Rep Power: 8
tian is on a distinguished road
Hi,

maybe you have to change the fvScheme in "FoamX/defaults/system/" folder?

laplacian((rho*(1|A(U))),p) Gauss linear corrected

I am not sure it is the reason...

bye
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   March 4, 2008, 04:43
Default Hi Marco, I think it could
  #16
Member
 
Francesco Boschetto
Join Date: Mar 2009
Location: Italy
Posts: 56
Rep Power: 8
francesco_b is on a distinguished road
Hi Marco,

I think it could be a BC problem, can you tell us which is your set of BC? Which changes did you make in the code?

Regards

Francesco
francesco_b is offline   Reply With Quote

Old   March 4, 2008, 05:04
Default Hi, with my old post I want
  #17
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 103
Rep Power: 8
tian is on a distinguished road
Hi,

with my old post I want to tell about you have to change the file fvScheme that you can open it in FoamX without error.

Bye
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian is offline   Reply With Quote

Old   March 4, 2008, 07:00
Default Hi Thomas, FoamX doesn't wo
  #18
New Member
 
Marco Zardetto
Join Date: Mar 2009
Location: Italy
Posts: 6
Rep Power: 8
sharmak is on a distinguished road
Hi Thomas,

FoamX doesn't work in spite of the change. Maybe other changes are necessary and I tried to do them, but nothing.... The error, when I create the new case, is this:

class org.omg.CORBA.UNKNOWN

Thanks anyway
sharmak is offline   Reply With Quote

Old   March 4, 2008, 07:05
Default Hi Francesco, maybe you are
  #19
New Member
 
Marco Zardetto
Join Date: Mar 2009
Location: Italy
Posts: 6
Rep Power: 8
sharmak is on a distinguished road
Hi Francesco,

maybe you are right, I modified the BC and the simulation starts but now it stops here:

Starting time loop

Time = 0.05

DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 9.71047e-06, No Iterations 5
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.14451e-06, No Iterations 5
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 2.265e-06, No Iterations 6
#0 Foam::error::printStack(Foam:stream&)
#1 Foam::sigFpe::sigFpeHandler(int)
...
#14 __gxx_personality_v0 at /usr/src/packages/BUILD/glibc-2.3/csu/../sysdeps/i386/elf/start.S:122
Floating point exception (core dumped)

This is my problem, I think it is not difficult to solve but I'm not able:

I have an air duct with an inlet section and an outlet section. The motion is produced by a delta p between the two sections. Flow is compressible, turbulent and steady.

These are some of the changes I made to modify rhoTurbFoam solver:
application rhoSimpleFoam (in contolDict file)
ddtSchemes steadyState (in fvSchemes file)
SIMPLE (in fvSolution file)

Honestly I don't know how to put all the BC at time 0 (R, k, epsilon?). Do you think the mistake is here? Or in the solution for p?

These are the BC I put:
U=0 on the wall, zeroGradient on inlet and outlet
p=fixedValue on inlet and outlet, zeroGradient on the wall
k, epsilon, R, T zeroGradient (and internalField ?)

Please can you help me?
Thanks
sharmak is offline   Reply With Quote

Old   March 4, 2008, 07:15
Default Hi Marco, can you give me y
  #20
Senior Member
 
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 103
Rep Power: 8
tian is on a distinguished road
Hi Marco,

can you give me your email? I will send you my files.

Bye
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance.
tian 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
Steadystate pressurebased compressible solver cosimobianchini OpenFOAM Running, Solving & CFD 1 July 19, 2010 14:45
Steadystate versus Transient solver pda OpenFOAM Running, Solving & CFD 1 July 11, 2007 08:36
Steadystate Vs Transient solver amitshah OpenFOAM Running, Solving & CFD 1 August 23, 2006 02:54
Steadystate Euler solver jelmer OpenFOAM Running, Solving & CFD 1 June 19, 2006 08:24
External Flow-compressible flow solver-lift/drag Tom Brown Main CFD Forum 7 December 29, 2000 14:41


All times are GMT -4. The time now is 19:54.