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

Pressure Driven Supersonic Flow in Converging-Diverging Nozzle

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

Like Tree4Likes
  • 2 Post By jaason
  • 2 Post By cdm

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 4, 2013, 16:12
Default Pressure Driven Supersonic Flow in Converging-Diverging Nozzle
  #1
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 12
cdm is on a distinguished road
After extensively reading through the forums and running a number of test cases, I still can't find suitable boundary conditions for my geometry. Hopefully someone here will have had previous experience and may be able to provide some input.

I'm trying to simulate a supersonic converging-diverging nozzle for comparison to experimental results. The known conditions are: stagnation p and T at inlet of 300kPa and 1200K, wall temp roughly 600K, Mach number of 2.03 at outlet, and based on a gas dynamics analytical calculation for design conditions, the outlet static pressure should be roughly 37kPa (isentropic).

I've tried a number of boundary conditions as below, implemented in OF 2.1.1, and using rhoCentralFoam with kOmegaSST. For all cases, the temperature was defined as follows:

Code:
    INLET
    {
        type            totalTemperature;
        gamma           1.34;
        T0              uniform 1200;
        value           uniform 1200;
    }
    OUTLET
    {
        type            zeroGradient;
    }
    WALL
    {
        type            fixedValue;
        value           uniform 600;
    }
Case 1:
This case has been the most promising; however, the flow settles with an outlet static pressure around 31 to 36kPa, below the expected 37kPa. Maybe this is an effect of a turbulent boundary layer or something else? Further reducing lInf to 0.001 causes the simulation to crash with a high pressure gradient near the outlet.

Pressure, P
Code:
    INLET
    {
        type            totalPressure;
        p0              uniform 300000;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1.34;
        value           uniform 300000;
    }
    OUTLET
    {
        type            waveTransmissive;
        value           uniform 37000;
        field           p;
        gamma           1.34;
        phi             phi;
        rho             rho;
        psi             psi;
        lInf            0.1;
        fieldInf        37000;
    }
    WALL
    {
        type            zeroGradient;
    }
Velocity, U
Code:
    INLET
    {
        type            zeroGradient;
    }
    OUTLET
    {
        type            zeroGradient;
    }
    WALL
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
Case 2:
This case implemented the desired static pressure at the outlet of 37kPa. The velocity was defined as in Case 1. The solution crashes after a few thousand iterations with a large pressure gradient at the outlet (similar to the small lInf in Case 1).

Pressure, P
Code:
    INLET
    {
        type            totalPressure;
        p0              uniform 300000;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1.34;
        value           uniform 300000;
    }
    OUTLET
    {
        type            fixedValue;
        value           uniform 37000;
    }
    WALL
    {
        type            zeroGradient;
    }
Case 3:
This case is the same as Case 2, except for a fluxCorrectedVelocity boundary implemented at the outlet. Again, there is a large pressure gradient (even larger this time) at the exit.

Velocity, U
Code:
    OUTLET
    {
        type            fluxCorrectedVelocity;
        phi             phi;
        rho             rho;
    }
Case 4:
This case implements the fluxCorrectedVelocity outlet velocity boundary condition (as in Case 3) with a waveTransmissive outlet pressure boundary condition (similar to Case 1, but lInf = 0.01). The simulation does not crash, but the pressure is far too high (between roughly 200 and 400 kPa) and fluctuates wildly with no apparent trend to settle.


Conclusions:
Does anyone have any suggestions for further boundary conditions that I could explore? Or any comments on the above cases? Case 1 is the only usable result so far, but I'm not confident that waveTransmissive is beneficial in this case since I don't have any shocks in the developed flow and I want the flow to relax to a given static pressure at the outlet, not artificially at a distance from the outlet. But when I specify the desired static pressure, a pressure gradient develops at the outlet and crashes the simulation.

Thanks for reading, and I appreciate any help.

Last edited by cdm; June 4, 2013 at 17:30. Reason: Case 4 BC correction.
cdm is offline   Reply With Quote

Old   June 5, 2013, 09:00
Default nozzle
  #2
Member
 
Join Date: Nov 2012
Posts: 83
Rep Power: 13
Henning86 is on a distinguished road
im using this BC:

pressureInletOutletVelocity for U

are you just simulating the nozzle without a wake?
Henning86 is offline   Reply With Quote

Old   June 5, 2013, 13:49
Default
  #3
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 12
cdm is on a distinguished road
I tried using that boundary condition on the outlet, but it didn't help when specifying a static pressure. What pressure BC do you use with it?

There are no shocks in my system, and I'm not concerned with the wake. What I'm trying to investigate is temperature distribution in a constant area downstream of the throat.
cdm is offline   Reply With Quote

Old   July 19, 2013, 11:32
Default
  #4
New Member
 
Jason Pearl
Join Date: Jul 2013
Location: Burlington
Posts: 10
Rep Power: 12
jaason is on a distinguished road
I'm working on a similar problem supersonic flow in a MicroNozzle.

for your temperature field have you tried using totalTemperature BC for the inlet with inletOutletTotalTemperature at the exit? might be something to look into. Though it might be difficult to use when the walls are not adiabatic.

Also for the velocity I had better luck with pressureDirectedInletOutletVelocity and pressureNormalInletOutletVelocity (just be aware of patch orientation for the latter)

for pressure I am using totalPressure and waveTransmissive.

Also what is the exact function of Linf. Does it weight the impact of the farfields effect on the BC? where larger would Linf decrease the effect of the farfield on the BC?

cheers
jbliu2015 and Isalm Shahboun like this.
jaason is offline   Reply With Quote

Old   July 21, 2013, 17:24
Default
  #5
cdm
Member
 
Join Date: May 2013
Location: Canada
Posts: 32
Rep Power: 12
cdm is on a distinguished road
jaason, thanks for your input. I'll have a look at your suggested boundary conditions. Right now I've been using (inlet / outlet, respectively) totalPressure / waveTransmissive for Pressure, totalTemperature / zeroGradient for Temperature, and zeroGradient / inletOutlet (with inletValue uniform (0 0 0) for the outlet).

I believe what you say is roughly how lInf is expected to work. It relaxes the BC pressure such that it would realise a farfield pressure pInf at a distance lInf from the BC patch. If lInf = 0 it would essentially be forcing a fixedValue. However, for large lInf, I think I read that sometimes there is trouble with convergence and having large oscillations in the outlet pressure. So far it appears to work fine for my supersonic flow nozzle with over-expanded conditions.
jbliu2015 and PeterPete like this.
cdm is offline   Reply With Quote

Old   November 14, 2017, 07:30
Default
  #6
Member
 
Bal Krishan
Join Date: Jul 2017
Posts: 43
Rep Power: 8
balkrishan is on a distinguished road
Guys, I need your help..!

Do anybody have reference(paper or any thing) in which nozzle contour points for CD nozzle are given, if yes please share.

@balkrishan010@gmail.com
balkrishan is offline   Reply With Quote

Old   May 8, 2018, 15:18
Default
  #7
New Member
 
Join Date: Apr 2018
Location: Ohio, USA
Posts: 15
Rep Power: 7
shang2 is on a distinguished road
Hi! I was wondering if you ever got this to work.

Thanks!
shang2 is offline   Reply With Quote

Old   August 28, 2019, 02:37
Default CFD in Pelton turbine nozzle
  #8
New Member
 
Bishal Pandey
Join Date: Aug 2019
Posts: 1
Rep Power: 0
Bishal Pandey is on a distinguished road
I am having a problem keeping the total head constant. I am simulating a nozzle used in a Pelton turbine at different nozzle position. it works if I use total mass flow inlet and static pressure at the outlet. but when I move to total pressure inlet, I get "floating point error", I have tried refining the mesh also but it does not help. Have anyone has any idea about such a problem? what might be the problem?
Bishal Pandey is offline   Reply With Quote

Reply

Tags
converging-diverging, nozzle, pressure driven, supersonic

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
velocity and pressure driven flow in openfoam jr33 OpenFOAM Running, Solving & CFD 1 August 7, 2013 12:00
Supersonic nozzle exit conditions Betty Main CFD Forum 13 May 24, 2012 01:26
UDF for transient pressure inlet, mass flow check at nozzle exit and volumetric heat kokoory FLUENT 0 August 17, 2011 03:07
Urgent: How to model a stationary sphere in a pressure driven flow using Ansys CFX? farhan OpenFOAM Running, Solving & CFD 1 April 14, 2009 15:34
flow simulation across a small fan jane luo Main CFD Forum 15 April 12, 2004 18:49


All times are GMT -4. The time now is 22:39.