Rocstar-legacy/BcondOutflowPerf__GL_8F90_source.html

 ! *********************************************************************

 ! * Rocstar Simulation Suite                                          *

 ! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *

 ! *                                                                   *

 ! * Illinois Rocstar LLC                                              *

 ! * Champaign, IL                                                     *

 ! * www.illinoisrocstar.com                                           *

 ! * sales@illinoisrocstar.com                                         *

 ! *                                                                   *

 ! * License: See LICENSE file in top level of distribution package or *

 ! * http://opensource.org/licenses/NCSA                               *

 ! *********************************************************************

 ! *********************************************************************

 ! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *

 ! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *

 ! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *

 ! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *

 ! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *

 ! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *

 ! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *

 ! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *

 ! *********************************************************************

 ! *****************************************************************************

 !

 ! Purpose: Set outflow boundary condition for one cell for gas-liquid model.

 !

 ! Description:

 !

 ! Input:

 !   bcOpt    boundary treatment: subsonic, supersonic, or mixed

 !   betaP    compressibility at const. pressure

 !   betaT    compressibility at const. temperature

 !   cvg      specific heat of gas at constant volume (boundary cell)

 !   cvl      specific heat of liquid at constant volume (boundary cell)

 !   cvv      specific heat of vapor at constant volume (boundary cell)

 !   pin      static pressure at boundary cell

 !   Po       reference pressure

 !   ro       reference density

 !   To       reference Temperature

 !   pout     given static outlet pressure

 !   Rg       gas constant

 !   Rv       vapor constant

 !   rho      density at boundary cell

 !   rhoe     density * total energy at boundary cell

 !   rhou/v/w density * velocity components at boundary cell

 !   rhogpg   density of gas * volume fraction of gas at boundary cell

 !   rhovpv   density of vapor * volume fraction of vapor at boundary cell

 !   sx/y/zn  components of ortho-normalized face vector (outward facing)

 !

 ! Output:

 !   rhob      density at boundary

 !   rhou/v/wb density * velocity components at boundary

 !   rhoeb     density * total energy at boundary

 !   rhogpgb   density of gas * volume fraction of gas at boundary

 !   rhovpvb   density of vapor * volume fraction of vapor at boundary

 !

 ! Notes: None.

 !

 ! *****************************************************************************

 !

 ! $Id: BcondOutflowPerf_GL.F90,v 1.3 2008/12/06 08:44:08 mtcampbe Exp $

 !

 ! Copyright: (c) 2006 by the University of Illinois

 !

 ! *****************************************************************************


 SUBROUTINE bcondoutflowperf_gl(bcOpt,ro,Po,To,betaP,betaT,cvl,cvv,cvg,Rg,Rv, &

                                 pout,sxn,syn,szn,rho,rhou,rhov,rhow,rhoe, &

                                 rhogpg,rhovpv,pin,rhob,rhoub,rhovb,rhowb, &

                                 rhoeb,rhogpgb,rhovpvb)


   USE moddatatypes

   USE modparameters

   USE modinterfaces, ONLY: mixtgasliq_c, &

                            mixtliq_c2_bp, &

                            mixtperf_c2_grt, &

                            mixtliq_d_dobpppobttto, &

                            mixtperf_d_prt, &

                            mixtperf_t_cveovm2


   IMPLICIT NONE


 ! ******************************************************************************

 ! Declarations and definitions

 ! ******************************************************************************


 ! ==============================================================================

 ! Arguments

 ! ==============================================================================


   INTEGER :: bcopt

   REAL(RFREAL), INTENT(IN) :: betap,betat,cvg,cvl,cvv,pin,po,pout,rg,rho,rhoe, &

                               rhogpg,rhou,rhov,rhovpv,rhow,ro,rv,sxn,syn,szn,to

   REAL(RFREAL), INTENT(OUT) :: rhob,rhoeb,rhogpgb,rhoub,rhovb,rhovpvb,rhowb


 ! ==============================================================================

 ! Locals

 ! ==============================================================================


   REAL(RFREAL) :: bg2,bl2,bv2,cg2,cl2,cm,ct2,cv2,cvm,deltp,e,fg,fv,ic2p,ic2pb, &

                   mach,rhg,rhgb,rhl,rhlb,rholpl,rhv,rhvb,rrhoc,t,tb,u,ub,v,vb, &

                   vel2,vfg,vfgb,vfl,vflb,vfv,vfvb,vnd,w,wb


 ! ******************************************************************************

 ! Start

 ! ******************************************************************************


   u    = rhou/rho

   v    = rhov/rho

   w    = rhow/rho

   vel2 = (u*u + v*v + w*w)


   rholpl = rho - rhovpv - rhogpg

   cvm    = (rholpl*cvl + rhovpv*cvv + rhogpg*cvg)/rho

   e      = rhoe/rho

   t      = mixtperf_t_cveovm2(cvm,e,vel2)


   cl2 = mixtliq_c2_bp(betap)

   cv2 = mixtperf_c2_grt(1.0_rfreal,rv,t)

   cg2 = mixtperf_c2_grt(1.0_rfreal,rg,t)


   rhl = mixtliq_d_dobpppobttto(ro,betap,betat,pin,po,t,to)

   rhv = mixtperf_d_prt(pin,rv,t)

   rhg = mixtperf_d_prt(pin,rg,t)


   bl2 = -betat/betap

   bv2 = rhv*rv

   bg2 = rhg*rg


   vfl = rholpl/rhl

   vfv = rhovpv/rhv

   vfg = rhogpg/rhg

   cm  = mixtgasliq_c(cvm,rho,pin,rhl,rhv,rhg,vfl,vfv,vfg,cl2,cv2,cg2, &

                      bl2,bv2,bg2)

   mach = sqrt(vel2)/cm


 ! *****************************************************************************

 ! Subsonic outflow

 ! *****************************************************************************


   IF ( mach < 1.0_rfreal ) THEN

     rrhoc = 1.0_rfreal/(rho*cm)

     deltp = pin - pout

     ub    = u + sxn*deltp*rrhoc

     vb    = v + syn*deltp*rrhoc

     wb    = w + szn*deltp*rrhoc


     tb = t - (deltp*pin)/(rho*rho*cvm*cm*cm)

     fg = (vfg*(rhg*cg2 - rho*cm*cm +(bg2*pin)/(rho*cvm)))/(rhg*cg2*rho*cm*cm)

     fv = (vfv*(rhv*cv2 - rho*cm*cm +(bv2*pin)/(rho*cvm)))/(rhv*cv2*rho*cm*cm)


     vfgb = vfg -  deltp*fg

     vfvb = vfv -  deltp*fv

     vflb = 1.0_rfreal - vfgb - vfvb


     ic2p  = vfl/cl2 + vfg/cg2 + vfv/cv2

     ic2pb = (bl2*vfl)/cl2 + (bg2*vfg)/cg2 + (bv2*vfv)/cv2


     rhob = rho - (rhv-rhl)*(vfl-vflb) - (rhg-rhl)*(vfg-vfgb) &

             - deltp*ic2p + ic2pb*(t-tb)

     rhlb = rhl - deltp/cl2 + (bl2*(t-tb))/cl2

     rhvb = rhv - deltp/cv2 + (bv2*(t-tb))/cv2

     rhgb = rhg - deltp/cg2 + (bg2*(t-tb))/cg2


 ! =============================================================================

 !   Special treatment to prevent "deltp" from changing the sign of velocity

 !   components. This may happen for very small u, v, w.

 ! =============================================================================


     vnd = ub*sxn + vb*syn + wb*szn


     IF ( vnd < 0.0_rfreal ) THEN ! inflow at outflow boundary

       ub = sign(1.0_rfreal,u)*max(abs(ub),abs(u))

       vb = sign(1.0_rfreal,v)*max(abs(vb),abs(v))

       wb = sign(1.0_rfreal,w)*max(abs(wb),abs(w))

     END IF ! vnd


     rhoub = rhob*ub

     rhovb = rhob*vb

     rhowb = rhob*wb

     rhoeb = (rhlb*vflb*cvl +rhgb*vfgb*cvg +rhvb*vfvb*cvv)*tb + &

               0.5_rfreal*rhob*(ub*ub + vb*vb + wb*wb)

     rhogpgb = rhgb*vfgb

     rhovpvb = rhvb*vfvb


 ! *****************************************************************************

 ! Supersonic flow

 ! *****************************************************************************


   ELSE

     rhob    = rho

     rhoub   = rhou

     rhovb   = rhov

     rhowb   = rhow

     rhoeb   = rhoe

     rhogpgb = rhogpg

     rhovpvb = rhovpv

   END IF ! mach


 ! *****************************************************************************

 ! End

 ! *****************************************************************************


 END SUBROUTINE bcondoutflowperf_gl


 ! *****************************************************************************

 !

 ! RCS Revision history:

 !

 ! $Log: BcondOutflowPerf_GL.F90,v $

 ! Revision 1.3  2008/12/06 08:44:08  mtcampbe

 ! Updated license.

 !

 ! Revision 1.2  2008/11/19 22:17:22  mtcampbe

 ! Added Illinois Open Source License/Copyright

 !

 ! Revision 1.1  2006/03/26 20:20:45  haselbac

 ! Initial revision

 !

 ! *****************************************************************************


modparameters
Definition: ModParameters.F90:42

sign
static SURF_BEGIN_NAMESPACE double sign(double x)
Definition: compute_curvature.C:31

max
Vector_n max(const Array_n_const &v1, const Array_n_const &v2)
Definition: Vector_n.h:354

mixtperf_d_prt
real(rfreal) function mixtperf_d_prt(P, R, T)
Definition: MixtPerf_D.F90:71

sqrt
double sqrt(double d)
Definition: double.h:73

v
*********************************************************************Illinois Open Source License ****University of Illinois NCSA **Open Source License University of Illinois All rights reserved ****Developed free of to any person **obtaining a copy of this software and associated documentation to deal with the Software without including without limitation the rights to and or **sell copies of the and to permit persons to whom the **Software is furnished to do subject to the following this list of conditions and the following disclaimers ****Redistributions in binary form must reproduce the above **copyright this list of conditions and the following **disclaimers in the documentation and or other materials **provided with the distribution ****Neither the names of the Center for Simulation of Advanced the University of nor the names of its **contributors may be used to endorse or promote products derived **from this Software without specific prior written permission ****THE SOFTWARE IS PROVIDED AS WITHOUT WARRANTY OF ANY **EXPRESS OR INCLUDING BUT NOT LIMITED TO THE WARRANTIES **OF FITNESS FOR A PARTICULAR PURPOSE AND **NONINFRINGEMENT IN NO EVENT SHALL THE CONTRIBUTORS OR **COPYRIGHT HOLDERS BE LIABLE FOR ANY DAMAGES OR OTHER WHETHER IN AN ACTION OF TORT OR **ARISING OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE **USE OR OTHER DEALINGS WITH THE SOFTWARE v
Definition: roccomf90.h:20

mixtperf_c2_grt
real(rfreal) function mixtperf_c2_grt(G, R, T)
Definition: MixtPerf_C.F90:101

mixtliq_c2_bp
real(rfreal) function mixtliq_c2_bp(Bp)
Definition: MixtLiq_C.F90:54

block::pin
blockLoc pin(const blockLoc &l) const
Definition: split.cpp:77

modinterfaces
Definition: ModInterfaces.F90:39

mixtliq_d_dobpppobttto
real(rfreal) function mixtliq_d_dobpppobttto(Dz, Bp, Bt, P, Po, T, To)
Definition: MixtLiq_D.F90:40

mixtgasliq_c
real(rfreal) function mixtgasliq_c(Cvm, D, P, Dl, Dv, Dg, VFl, VFv, VFg, Cl2, Cv2, Cg2, Bl2, Bv2, Bg2)
Definition: MixtGasLiq_C.F90:39

mixtperf_t_cveovm2
real(rfreal) function mixtperf_t_cveovm2(Cv, Eo, Vm2)
Definition: MixtPerf_T.F90:70

bcondoutflowperf_gl
subroutine bcondoutflowperf_gl(bcOpt, ro, Po, To, betaP, betaT, cvl, cvv, cvg, Rg, Rv, pout, sxn, syn, szn, rho, rhou, rhov, rhow, rhoe, rhogpg, rhovpv, pin, rhob, rhoub, rhovb, rhowb, rhoeb, rhogpgb, rhovpvb)
Definition: BcondOutflowPerf_GL.F90:67

moddatatypes
Definition: ModDataTypes.F90:40