rimseal.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	rimseal (node1, node2, nodem, nelem, lakon, kon, ipkon, nactdog, identity, ielprop, prop, iflag, v, xflow, f, nodef, idirf, df, cp, R, physcon, dvi, numf, set, mi, ttime, time, iaxial, co, vold)

Function/Subroutine Documentation

rimseal()

subroutine rimseal	(	integer	node1,
		integer	node2,
		integer	nodem,
		integer	nelem,
		character*8, dimension(*)	lakon,
		integer, dimension(*)	kon,
		integer, dimension(*)	ipkon,
		integer, dimension(0:3,*)	nactdog,
		logical	identity,
		integer, dimension(*)	ielprop,
		real*8, dimension(*)	prop,
		integer	iflag,
		real*8, dimension(0:mi(2),*)	v,
		real*8	xflow,
		real*8	f,
		integer, dimension(*)	nodef,
		integer, dimension(*)	idirf,
		real*8, dimension(*)	df,
		real*8	cp,
		real*8	R,
		real*8, dimension(*)	physcon,
		real*8	dvi,
		integer	numf,
		character*81, dimension(*)	set,
		integer, dimension(*)	mi,
		real*8	ttime,
		real*8	time,
		integer	iaxial,
		real*8, dimension(3,*)	co,
		real*8, dimension(0:mi(2),*)	vold
	)

!     
!     rimseal element
!
!     author: Yannick Muller
!     
      implicit none
!     
      logical identity
      character*8 lakon(*)
      character*81 set(*)
!     
      integer nelem,nactdog(0:3,*),node1,node2,nodem,numf,
     &     ielprop(*),nodef(*),idirf(*),index,iflag,mi(*),
     &     inv,ipkon(*),kon(*),kgas,nelem_in,nelem_out,iaxial,
     &     element0,node10,node20,node11,node21,node12,node22,node_cav,
     &     node_main,node_main2,node_in1,node_out1,node_in2,node_out2
!
      real*8 prop(*),v(0:mi(2),*),xflow,f,df(*),kappa,r,a,d,
     &     p1,p2,t1,t2,aeff,c1,c2,c3,cd,cp,physcon(*),p2p1,km1,dvi,
     &     kp1,kdkm1,tdkp1,km1dk,x,y,ca1,cb1,ca2,cb2,dt1,alambda,
     &     reynolds,pi,xflow_oil,s,tcav,pcav,pmin,pmax,ttime,time,
     &     tref,alpha1, alpha2, alpha3, gf,kf,mrtap_ref_ein,
     &     mrtap_ref_aus, m_ref_ein, m_ref_aus,maus_zu_mref,
     &     mein_zu_mref, a_aus, a_ein, a_ges,m_aus, m_ein, m_sperr,
     &     vold(0:mi(2),*),co(3,*)
!
      pi=4.d0*datan(1.d0)   

      if (iflag.eq.0) then
         identity=.true.
!     
         if(nactdog(2,node1).ne.0)then
            identity=.false.
         elseif(nactdog(2,node2).ne.0)then
            identity=.false.
         elseif(nactdog(1,nodem).ne.0)then
            identity=.false.
         endif
!     
      elseif (iflag.eq.1) then      
!     
         p1=v(2,node1)
         call rimseal_calc(p1)
      elseif (iflag.eq.2) then
!
      elseif (iflag.eq.3) then
!
      endif
      return