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.

[shaun]

check out rhoSimpleFoam.

this will be what you are looking for.

Regards

Shaun.D

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.

[graser]

You just have to set up the case by hand without FoamX.

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.

[graser]

look at these tutorials:
rhoExplicitPorousSimpleFoam
rhoImplicitPorousSimpleFoam
they use rhoSimpleFoam solver

[shaun]

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

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

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

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

[sharmak]

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

[cedric_duprat]

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

[francesco_b]

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

[sharmak]

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

[tian]

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

[francesco_b]

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

[tian]

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

[sharmak]

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]

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

[tian]

Hi Marco,

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

Bye