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

Help running an invicid compressible nozzle flow using openFOAM

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 7, 2010, 17:32
Post Help running an invicid compressible nozzle flow using openFOAM
  #1
New Member
 
Radoslav Bozinoski
Join Date: Feb 2010
Location: Davis, CA
Posts: 4
Rep Power: 7
rbozinoski is on a distinguished road
Hi Foamers,

I am very new to OpenFOAM and have only been playing around with it for about a month(very unsuccessfully I might add). I have heard great things about this cfd toolkit and wanted to start getting my feet wet with some of its capabilities. I've gone through some of the tutorials and have read a lot of journal articles where OpenFOAM was used and am convinced that OpenFOAM will be a great tool for my future research aspirations. Unfortunately, my first few attempts at using it for cases other than what's packaged with the tutorials have not gone very well. So what I finally did was focus all my efforts on learning OpenFOAM by first using the following Fluent tutorial.

http://courses.cit.cornell.edu/fluent/nozzle/step4.htm

This tutorial outlines a simple inviscid nozzle flow problem and I wanted to reproduce the results using OpenFOAM. I converted the mesh using 'fluentMeshToFoam' and then tried to first generate an initial flow field using potentialFoam and then launch an inviscid case using rhoSonicFoam. I was unsuccessful with both. I've attached a tarball of my case files and was hoping someone could help me understand what I am doing wrong and how to properly launch this case. I will admit that I don't have a clear understanding of how to properly implement the boundary conditions or how to get more information regarding the current capabilities in OpenFOAM.

Once again, thank you all in advance and any assistance would be greatly appreciated.
Attached Files
File Type: gz nozzle-fluent-OpenFOAM.tar.gz (65.5 KB, 205 views)

Last edited by rbozinoski; February 8, 2010 at 01:31.
rbozinoski is offline   Reply With Quote

Old   February 10, 2010, 10:20
Default
  #2
New Member
 
Radoslav Bozinoski
Join Date: Feb 2010
Location: Davis, CA
Posts: 4
Rep Power: 7
rbozinoski is on a distinguished road
Can anyone please offer any advice or help. Still very much lost with OpenFOAM.

Thank you in advance.
rbozinoski is offline   Reply With Quote

Old   February 11, 2010, 11:36
Default
  #3
Member
 
Carlos Xisto
Join Date: Nov 2009
Location: Covilhã, Portugal
Posts: 53
Rep Power: 7
xisto is on a distinguished road
Send a message via MSN to xisto
Hi Radoslav

What is the problem? you need to be more specific.

I have seen your case (rhoSonicFoam) and you make some mistakes.

In the /0/p the internalField of pressure must be the value that you specify at the outlet.

Try also to lower your deltaT.

Best luck

Cheers,

CX
xisto is offline   Reply With Quote

Old   February 12, 2010, 01:08
Default
  #4
New Member
 
Radoslav Bozinoski
Join Date: Feb 2010
Location: Davis, CA
Posts: 4
Rep Power: 7
rbozinoski is on a distinguished road
Hello Carlos,

Thank you so much for replying. I took your advice and set the internal pressure field to the exit pressure and lowered my time step but I still get the following error:

**********************

Create time

Create mesh for time = 0

Reading thermodynamicProperties

Reading field p

Reading field T

Reading field U


Starting time loop

Time = 1e-06


Max Courant Number = 0.00504052



request for surfaceScalarField phi from objectRegistry region0 failed
available objects of type surfaceScalarField are

4
(
phiv
limitedLinearLimiter(rho)
differenceFactors_
weightingFactors
)
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam:perator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSonicFoam"
#3 Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const& Foam:bjectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvsPatchField, Foam::surfaceMesh> >(Foam::word const&) const in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 Foam::totalPressureFvPatchScalarField::updateCoeff s(Foam::Field<Foam::Vector<double> > const&) in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensionSet const&) in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Fo am::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libfiniteVolume.so"
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSonicFoam"
#8 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSonicFoam"
#9 main in "/home/rbozinoski/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSonicFoam"
#10 __libc_start_main in "/lib/libc.so.6"
#11 _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/x86_64/elf/start.S:116


From function objectRegistry::lookupObject<Type>(const word&) const
in file /home/rbozinoski/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 140.

FOAM aborting

Aborted

******************************

I tracked down the problem to the total pressure/temperature boundary conditions. It seems to be complaining about the variable phi. I have tried using other bc at the inlet but all eventually diverge. I was wondering if you are having the same issues. Also, I am not sure if my control/thermodynamicProperties file is correct. I am simply trying to run this case inviscid and compressible.

I have attached my latest rhoSonic files. You will have to use the fluentMeshToFoam utility to generate the OpenFOAM mesh to run the simulation.

Any further advice will be greatly appreciated.

Thank you again for your help.
Attached Files
File Type: gz rhoSonic.tar.gz (65.2 KB, 70 views)
rbozinoski is offline   Reply With Quote

Old   February 12, 2010, 06:22
Default
  #5
Member
 
Carlos Xisto
Join Date: Nov 2009
Location: Covilhã, Portugal
Posts: 53
Rep Power: 7
xisto is on a distinguished road
Send a message via MSN to xisto
I never used a density based solver in OpenFoam, just a Pressured based.
If you want to try with the rhoPisoFoam check the boundary conditions in this post (the 0 file).

http://www.cfd-online.com/Forums/ope...tml#post244986

Another thing it seems that you have a pressure realation very high resulting in Ma=2.9, in the other and the the velocity (internalField and intial value) is to low 30 m/s. This is why it crashes (imagine you crashing against a wall).

My advice is to increase the velocity in the internal field and at the inlet to 2.9 Ma.

Good Luck

CX
xisto is offline   Reply With Quote

Old   February 12, 2010, 11:55
Default
  #6
New Member
 
Radoslav Bozinoski
Join Date: Feb 2010
Location: Davis, CA
Posts: 4
Rep Power: 7
rbozinoski is on a distinguished road
Carlos,

Yeah I think you are correct regarding the velocity initialization. I'll have to give that a try to see what happens. However, I am pretty sure that I will still get the same error I posted earlier regarding the total pressure and temperature bc's.

Thanks for pointing me to your other post. I think I will switch to the rhoPisoFoam solver and try and reproduce that subsonice case first and then the transonic case. I noticed that you posted your initialization files and was wondering if you would mind also posting the grid and any other files I would need to run this simulation. I just wanted to start from a working case and then build up to a transonic from there.

Thank you again for your help,

-RB
rbozinoski is offline   Reply With Quote

Old   February 12, 2010, 12:24
Default
  #7
Member
 
Carlos Xisto
Join Date: Nov 2009
Location: Covilhã, Portugal
Posts: 53
Rep Power: 7
xisto is on a distinguished road
Send a message via MSN to xisto
My email is carlos_xisto@hotmail.com

Just send me a request and I will mail the files to you

Cheers

CX
xisto is offline   Reply With Quote

Old   January 14, 2013, 03:16
Default
  #8
Member
 
Baris (Heewa)
Join Date: Jan 2013
Location: Japan
Posts: 98
Rep Power: 4
shipman is on a distinguished road
Dear Carlos,

Hi,
I am new user in openFoam

I try to run the cavitatingFoam for the throttle. For pressure at the inlet BC set as totalPressure so by changing the velocity i try to get new results but give following errors. COuld you help to me what does it mean and how can i correct it?

Thanks.

Phase-change corrected p max-min : 6.37897e+07 400
max(U) 14881.4
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/opt/openfoam211/platforms/linux64GccDPOpt/bin/cavitatingFoam"
Floating point exception (core dumped)
shipman is offline   Reply With Quote

Old   April 5, 2014, 04:13
Default Rocket Nozzle flow simulation in OpenFOAM
  #9
New Member
 
Ranjan
Join Date: Apr 2014
Posts: 7
Rep Power: 3
Ranjan is on a distinguished road
Dear Foamers

Can anyone suggest a solver I can use for simulating supersonic flow in a rocket nozzle?

Thanks
Ranjan is offline   Reply With Quote

Old   April 5, 2014, 12:29
Default rhoCentralFoam is OK
  #10
New Member
 
何刚
Join Date: Jan 2011
Posts: 25
Rep Power: 6
hg2lf is on a distinguished road
Quote:
Originally Posted by Ranjan View Post
Dear Foamers

Can anyone suggest a solver I can use for simulating supersonic flow in a rocket nozzle?

Thanks


you need to specify all variables except one at inlet (f.g. fixed values for temperature (static), and extrapolate all variables from outlet (zero gradient).
You can use totalPressure, total temperature, zeroGradient(U) for inlet of the nozzle, and specify static pressure and temperature.
hg2lf is offline   Reply With Quote

Reply

Tags
compressible, invicid, nozzle flow, openfoam

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
Incompressible and compressible flow. Confused student. Main CFD Forum 22 January 30, 2014 07:32
compressible flow calculation error using rhoSimpleFoam solver student4326 OpenFOAM 6 February 10, 2012 10:36
Statically Compiling OpenFOAM Issues herzfeldd OpenFOAM Installation 21 January 6, 2009 10:38
Axisymmetirc Compressible Flow in a C/c and Nozzle Nauman FLUENT 0 December 21, 2004 08:40
compressible channel flow.. R.D.Prabhu Main CFD Forum 0 July 17, 1998 17:23


All times are GMT -4. The time now is 20:46.