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

Stability Problem with sonicFoam for Nozzle Flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 19, 2009, 05:25
Default Stability Problem with sonicFoam for Nozzle Flow
  #1
Member
 
Julian Krick
Join Date: May 2009
Location: Guelph
Posts: 88
Rep Power: 7
Julian K. is on a distinguished road
Hi,

I have a stability problem in my flow simulation. Here's the configuration:

I'd like to simulate the flow air through a convergent-divergent nozzle. The geometry starts at the inlet and stops at the outlet of the nozzle. The flow should be critical in the throat. For this I use a pressure difference of dp=300mbar, whereas at the inlet I apply atmospheric pressure. The the calculation I use the sonicFoam solver with 2nd order, bounded interpolation schemes (limitedLinear with phi=1/0.5 and MUSCL) for the divergence terms.


I set the BCs as follows:
/0/p:
Code:
    inlet
    {
        type            totalPressure;
        p0               uniform 1.01325e+05;
        U                U;
        phi              phi;
        rho              none;
        psi              none;
        gamma        1.4; 
        value           uniform 1.01325e+05;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 0.71325e+05;
    }
/0/U:
Code:
    inlet
    {
        type            pressureInletVelocity;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
What happens is the following:
The simulation starts fine. The p-, U- and T-fields develop as they should. Shocks occur, which cause a separation of the boundary layer resulting in a creation vortices, which are eventually transported towards the outlet.
However, at some point the solution becomes unstable. At the outlet there is a drastic increase in velocity and pressure (of order > 10^30). Also negative pressure occurs near the outlet. Then the simulation crashes.

I suppose the reason for this effect is the BC at the outlet. What I could do is to fix the pressure, velocity and temperature at the inlet and have 'zeroGradients' at the outlet. However, for my analysis it is necessary to induce the flow via a pressure difference between inlet and outlet.
__________________
grid generation: ICEM CFD 13.0
solver: CFX 13.0
Julian K. is offline   Reply With Quote

Old   June 19, 2009, 08:42
Default
  #2
Member
 
Julian Krick
Join Date: May 2009
Location: Guelph
Posts: 88
Rep Power: 7
Julian K. is on a distinguished road
As suggested by prapanj in this thread (Pressure Inlet Velocity) I used the 'waveTransmissive' BC for the pressure at the outlet. Now, the simulation seem to be stabilized, however it's still running so that, yet, I cannot tell if it works. I hope so .
__________________
grid generation: ICEM CFD 13.0
solver: CFX 13.0
Julian K. is offline   Reply With Quote

Old   June 29, 2009, 19:46
Default
  #3
New Member
 
Giro
Join Date: Mar 2009
Posts: 13
Rep Power: 7
girogirozakk is on a distinguished road
Hi

I try to solve nozzle-flow and my sample-files are opened on my URL.
http://giropenfoam.web.fc2.com/sample/index_sample.html

This is how you can help?
girogirozakk 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
Can 'shock waves' occur in viscous fluid flows? diaw Main CFD Forum 105 November 20, 2009 04:19
Problem in cavity flow problem saad Main CFD Forum 4 November 1, 2007 08:45
Need help with BC's for a internal flow problem ravi FLUENT 2 March 16, 2005 00:45
transient compressible flow problem (urgent plz) jehanzeb FLUENT 5 August 3, 2004 08:04
Periodic flow boundary condition problem sudha FLUENT 3 April 28, 2004 08:40


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