dragcoeff.f

Go to the documentation of this file.
00001 C:\opentelemac\v7p0\sources\telemac2d\dragcoeff.f
00002 !
00042                      SUBROUTINE DRAGCOEFF
00043 !                    ********************
00044 !
00045      & (V, D, VK, CW)
00046 !
00047 !***********************************************************************
00048 ! TELEMAC2D   V6P1                                   21/08/2010
00049 !***********************************************************************
00050 !
00051 !
00052 !
00053 !
00054 !
00055 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00056 !| CW             |<--| DRAG COEFFICIENT BEHIND A CYLINDER
00057 !| D              |-->| DIAMETER
00058 !| V              |-->| VELOCITY UPSTREAM
00059 !| VK             |-->| LAMINAR VISCOSITY
00060 !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
00061 !
00062       IMPLICIT NONE
00063 !
00064 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00065 !
00066       DOUBLE PRECISION, INTENT(IN)  :: V, D, VK
00067       DOUBLE PRECISION, INTENT(OUT) :: CW
00068 !
00069 !+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
00070 !
00071       DOUBLE PRECISION              :: RE
00072 !
00073 !=======================================================================!
00074 !=======================================================================!
00075 !
00076       RE = V * D / VK
00077 !
00078       IF (RE.LE.800.0D0) THEN
00079         CW = 3.07D0 / RE**(0.168D0)
00080       ELSEIF(RE.LE.6000.D0) THEN
00081         CW = 1.D0
00082       ELSEIF(RE.LE.11000.0D0) THEN
00083         CW = 1.0D0+0.2D0*(RE-6000.D0)/5000.D0
00084       ELSE
00085         CW = 1.2D0
00086       ENDIF
00087 !
00088 !=======================================================================!
00089 !=======================================================================!
00090 !
00091       RETURN
00092       END

Generated on Fri Aug 31 2013 18:12:58 by S.E.Bourban (HRW) using doxygen 1.7.0