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

Subsonic Inflow/Outflow Boundarty Conditions

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 29, 2006, 16:11
Default Subsonic Inflow/Outflow Boundarty Conditions
  #1
Ruben
Guest
 
Posts: n/a
Hello,

I'm testing my code for solving the compressible inviscid Euler equations with very simple examples. Noy I'm testing the subsonic case.

I know that in the inflow total temperature, total pressure and the angle/s must to be specified. But, I'd like to know how to "translate" these conditions to the conservatives quantities (density, momentum and energy).

For the outflow I read that the most common quantity to be imposed is the pressure.

Many thanks!
  Reply With Quote

Old   May 29, 2006, 21:52
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #2
khaiching
Guest
 
Posts: n/a
I presume that you may try to use the isentropic relations??
  Reply With Quote

Old   May 30, 2006, 01:32
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #3
Ruben
Guest
 
Posts: n/a
I found the isentropic relations here

http://www.grc.nasa.gov/WWW/K-12/airplane/isentrop.html

but:

- how can I determine the Mach to compute the density, pressure and temperature. I pressume that I need the Riemann invariants but I don't know...

- what is the total density? Is it a given data?

  Reply With Quote

Old   May 30, 2006, 08:56
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #4
os
Guest
 
Posts: n/a
dear ruben,

with fixed total pressure and temperature the total density is also known (gas equation). To get the Mach number the most simple possibility is to extrapolate the pressure or the density from the inner flowfield. Together with the prescribed total state you can determine the Machnumber and from the Machnumber and then prescribed flow angle all conservation variables.

  Reply With Quote

Old   May 30, 2006, 10:02
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #5
Ruben
Guest
 
Posts: n/a
Ok, thanks and excuse me for my trivial question about density...

About the Mach: I don't understand when the authors say "I impose a inlet Mach number of ..." Can I impose a inlet mach number?

  Reply With Quote

Old   May 30, 2006, 19:36
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #6
Peter Rozovski
Guest
 
Posts: n/a
When you solve the compressible inviscid Euler equations at the inlet subsonic boundary you cannot set all values because 1) you have the characteristics in 1D case (or characteristics surface in 2D/3D cases) running to the inlet boundary from inside calculation region and 2) you don't know the inlet speed (or inlet Mach number). So you have to use total temperature T0, total pressure P0, angle(s) in 2D(3D) cases and Riemann invariant condition from characterictics coming from inside region from previous time step. Such conditions are called "non-reflecting" conditions because they allow disturbances of flow parameters to leave calculation region and so increase speed of convergence of numerical solution.

Above means that you write system of non-linear equations for every inlet grid point. For 2D case we have

P/P0 = Fp isoentopic (cappa) T/T0 = Ft isoentopic (cappa) v/u = angle( given funtion) f(u,v,p,t, parameters at the nearest points at previous time layer) = 0 (Riemann invariant for corresponding characteristic)

and solve it by Newton method, for example.

Hope, above answers your question.

  Reply With Quote

Old   May 31, 2006, 04:58
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #7
Ruben
Guest
 
Posts: n/a
Thanks Peter!
  Reply With Quote

Old   May 31, 2006, 12:13
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #8
Peter Rozovski
Guest
 
Posts: n/a
Ruben,

One correction about reflecting and non-reflecting boundary conditions (in my previous post I was non-correct). Non-reflecting boundary conditions mean that you set Riemann invariants along characteristics coming from outside calculation region. In 1D subsonic flow, for example, we have two such characteristics - C+ and C0. So we have to set 2 conditions like h+u2/2 = H0 = cpT0 (along C0) and u+2a/(cappa-1)= 2a0/(cappa-1) (along C+). I don't have textbook in front of me so better to verify these relations. If you simulate case when cappa is not constant (cp = cp(T)), then you can use similar approach but formulae would be slightly complicated (for Newton iterations it doesn't matter).

  Reply With Quote

Old   June 1, 2006, 03:59
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #9
Ruben
Guest
 
Posts: n/a
My boundary conditions are working!!

Thanks a lot another time Peter.
  Reply With Quote

Old   June 1, 2006, 04:00
Default Re: Subsonic Inflow/Outflow Boundarty Conditions
  #10
os
Guest
 
Posts: n/a
It is not possible to impose the Mach number at inlet. The precsibed total state at inlet does not fix the Mach number. The Machnumber depends also from the exit pressure, i.e. the boundary condition at outlet. It is like a test rig. At inlet you have a vessel with a total state and guide vanes for the inlet angle. At the end of the test rig the flow is blown for example in an exhaust pipe with given outlet pressure.
  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
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
subsonic inlet-exit boundary conditions Cfd Student Main CFD Forum 6 July 26, 2005 10:16
A problem about setting boundary conditions lyang Main CFD Forum 0 September 19, 1999 18:29


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