wpi.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	wpi (W, PI, Q, SQTT, kappa, RGAS)

Function/Subroutine Documentation

wpi()

subroutine wpi	(	real*8	W,
		real*8	PI,
		real*8	Q,
		real*8	SQTT,
		real*8	kappa,
		real*8	RGAS
	)

!
!-----------------------------------------------------------------------
!                                                                      |
!     Dieses Unterprogramm berechnet die Stroemungs-Geschwindigkeit    |
!     fuer das eingegebene Druckverhaeltnis PI.                        |
!                                                                      |
!     Eingabe-Groessen:                                                |
!       PI     = Druckverhaeltnis PS/PT                                |
!       Q      = reduzierter Durchsatz                                 |
!       SQTT   = SQRT (Totaltemperatur)                                |
!                                                                      |
!     Ausgabe-Groessen:                                                |
!       W      = Stroemungs-Geschwindigkeit                            |
!                                                                      |
!-----------------------------------------------------------------------
!
      IMPLICIT CHARACTER*1 (a-z)
!       INCLUDE 'comkapfk.inc'
      real*8    w, pi, q, sqtt,kappaq,kappa,rgas,pikrit,kappah,wkritf
!
!-----------------------------------------------------------------------
!
      kappaq = 1/kappa
      pikrit = (2./(kappa+1.)) ** (kappa/(kappa-1.))
!
      kappah = 2. * kappa / (kappa + 1.)
      wkritf = sqrt( kappah * rgas )
!
      IF (pi.GE.1.) THEN
!       Druckverhaeltnis groesser gleich 1
        w    = 0.
      ELSEIF (pi.GT.pikrit) THEN
!       Druckverhaeltnis unterkritisch
        IF (q.GT.0.) THEN
          w    = q * rgas * sqtt * pi**(-kappaq)
        ELSE
          w    = 0.
        ENDIF
      ELSEIF (pi.GT.0.) THEN
!       Druckverhaeltnis ueberkritisch
        w    = wkritf * sqtt
      ELSE
!       Druckverhaeltnis ungueltig
        w    = 1.e20
      ENDIF
!
      RETURN