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

Transonic jet simulation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 24, 2006, 14:40
Default Hi everybody, I am a beginn
  #1
New Member
 
Magnus Åberg
Join Date: Mar 2009
Posts: 11
Rep Power: 8
magnus is on a distinguished road
Hi everybody,

I am a beginner OpenFOAM user with a wish of doing an investigation of the acoustics generated by a transonic jet through a nozzle. For this I tried rhoSonicFoam with initial conditions for the velocity field given by an incompressible solution of the same case.
I then start the compressible case using a really low Courant number. After some 6000 iterations the calculation still looks good but then blows up very quickly.
This is the case both when i use higher order vanLeer but also when I tried the upwind scheme.

As I am a beginner user maybe this is not the right solver that I have choosen.
Does anybody have ideas/suggestions of what I can do?

Thanks for your kind answers!
/Magnus
magnus is offline   Reply With Quote

Old   January 25, 2006, 04:20
Default It is usually something physic
  #2
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
It is usually something physical happening that causes this blowing up after such a long time of stable running. Check the solution just before it blows up. Maybe there is something happening (inflow?) at the outlet? Is the initial shock (from starting from an incompatible solution) reaching the outlet. Or is the solution finally going supersonic?
mattijs is offline   Reply With Quote

Old   January 25, 2006, 04:38
Default IC and BC should always be an
  #3
New Member
 
Shiuh-Hwa Shyu (Steven)
Join Date: Mar 2009
Location: Taiwan
Posts: 17
Rep Power: 8
sshyu is on a distinguished road
IC and BC should always be an issue. A bad IC might eventially kills your calculation. But sometimes, the code might get survived and give you final solution. This would depend on the physical problem. As for BC, it would either give you a wrong solution, or just blow up. In my experience, the FDM tends to survive for quite a long period of time before it blows up, or gives you wrong solution. As for the spectral-like method, it would diverge quite quickly. This is a good thing though.
sshyu is offline   Reply With Quote

Old   January 25, 2006, 11:17
Default Many Thanks for your input!
  #4
New Member
 
Magnus Åberg
Join Date: Mar 2009
Posts: 11
Rep Power: 8
magnus is on a distinguished road
Many Thanks for your input!

I restarted my calculations, using first order schemes, with zero for pressure and velocity in the interior. I can see that the problem comes from the inlet patch
I looked throught my boundary conditions and I found that maybe I misunderstand the meaning of pressureInletVelocity. What is the value i need to specify there?

I want to specify totalpressure and totaltemperature what is the correct velocity BC for this setup?

/Magnus
magnus is offline   Reply With Quote

Old   January 25, 2006, 11:19
Default Many Thanks for your input!
  #5
New Member
 
Magnus Åberg
Join Date: Mar 2009
Posts: 11
Rep Power: 8
magnus is on a distinguished road
Many Thanks for your input!

I restarted my calculations, using first order schemes, with zero for pressure and velocity in the interior. I can see that the problem comes from the inlet patch
I looked throught my boundary conditions and I found that maybe I misunderstand the meaning of pressureInletVelocity. What is the value i need to specify there?

I want to specify totalpressure and totaltemperature what is the correct velocity BC for this setup?

/Magnus
magnus is offline   Reply With Quote

Old   February 1, 2006, 12:54
Default Hi, I realize that what i w
  #6
New Member
 
Magnus Åberg
Join Date: Mar 2009
Posts: 11
Rep Power: 8
magnus is on a distinguished road
Hi,

I realize that what i wanted to do is not implemented.
So I ask instead if somebody can give me a suggestion of which solver to use for a transonic turbulent jet with moving shocks. The inlet is subsonic and I whish to use LES for closure.

Thank you
/Magnus
magnus is offline   Reply With Quote

Old   February 1, 2006, 13:23
Default coodles with "transonic true
  #7
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
coodles with
"transonic true;"
added to the PISO subdictionary.
No way it will capture shocks properly though.
eugene is offline   Reply With Quote

Old   February 11, 2006, 15:57
Default Could I modify sonicFoam (merg
  #8
New Member
 
Thomas Hauser
Join Date: Mar 2009
Posts: 1
Rep Power: 0
thauser is on a distinguished road
Could I modify sonicFoam (merge with coodles) to create a LES solver for compressible flow with shocks? Or are there better ways?
thauser is offline   Reply With Quote

Old   June 14, 2006, 02:43
Default I have the same question as wh
  #9
Member
 
Dihao Tang
Join Date: Mar 2009
Posts: 78
Rep Power: 8
tangd is on a distinguished road
I have the same question as what Magnus asked:
"I want to specify totalpressure and totaltemperature what is the correct velocity BC for this setup?"
I also use coodles, hope somebody can give me the answer, thanks!
tangd 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
AVL Fire transonic simulation Ainil Main CFD Forum 1 August 24, 2010 15:53
transonic simulation with OpenFoam Shyam Main CFD Forum 0 June 12, 2008 02:39
Simulation of transonic flow over NACA0012 airfoil MSc Student CD-adapco 2 August 9, 2006 13:49
Compressible transonic airfoil RAE2822 simulation Stefano CD-adapco 9 June 21, 2006 10:47
Transonic DNS cfd-freak Main CFD Forum 1 April 16, 2003 12:39


All times are GMT -4. The time now is 13:15.