CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   SteadyState solver for compressible flow (https://www.cfd-online.com/Forums/openfoam-solving/58316-steadystate-solver-compressible-flow.html)

matteo_gautero February 15, 2008 08:38

Hi to all, I have a flow with
 
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 February 15, 2008 17:09

check out rhoSimpleFoam. th
 
check out rhoSimpleFoam.

this will be what you are looking for.

Regards

Shaun.D

matteo_gautero February 18, 2008 09:00

Thanks Shaun. I have a problem
 
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 February 19, 2008 03:24

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

myting February 19, 2008 03:37

Hi, I also want to know how
 
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 February 19, 2008 03:56

look at these tutorials: rhoE
 
look at these tutorials:
rhoExplicitPorousSimpleFoam
rhoImplicitPorousSimpleFoam
they use rhoSimpleFoam solver

shaun February 19, 2008 06:37

you could also simple search t
 
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

myting February 19, 2008 06:59

Hi, thanks, I think, I will
 
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

matteo_gautero February 26, 2008 06:47

Hi to all, I set up a simple
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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

ttian February 26, 2008 07:02

Hi, can you give me your em
 
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 March 3, 2008 11:14

Hi, can anybody send me the
 
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 March 3, 2008 11:23

Hi Marco I think you want t
 
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 March 3, 2008 11:23

Hi Marco, when you save the
 
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 http://www.cfd-online.com/OpenFOAM_D...part/happy.gif

Regards

Francesco

sharmak March 3, 2008 12:51

Thanks Francesco and Cedric,
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#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 March 4, 2008 03:15

Hi, maybe you have to chang
 
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 March 4, 2008 03:43

Hi Marco, I think it could
 
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 March 4, 2008 04:04

Hi, with my old post I want
 
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 March 4, 2008 06:00

Hi Thomas, FoamX doesn't wo
 
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 March 4, 2008 06:05

Hi Francesco, maybe you are
 
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&)
#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 March 4, 2008 06:15

Hi Marco, can you give me y
 
Hi Marco,

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

Bye


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