ViscousFluxes.F90

Go to the documentation of this file.

! *********************************************************************
! * 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: Compute viscous fluxes and add them to dissipation.
!
! Description: Distinction made between laminar and turbulent (total)
!              viscous fluxes; laminar flux is computed by generic routine
!              ViscousFluxEddy if TurbModel is inactive, 
!              else turbulent flux routine is called; upon output 
!              viscous fluxes added to numerical dissipation fluxes
!
! Input: 
!   region      Data of current region.
!
! Output: None.
!
! Notes: None.
!
! ******************************************************************************
!
! $Id: ViscousFluxes.F90,v 1.14 2008/12/06 08:44:10 mtcampbe Exp $
!
! Copyright: (c) 2001-2006 by the University of Illinois
!
! ******************************************************************************

SUBROUTINE viscousfluxes(region)

  USE moderror
  USE modparameters
  USE moddatatypes
  USE modbndpatch, ONLY: t_patch
  USE modglobal, ONLY: t_global
  USE moddatastruct, ONLY: t_region
  USE modbndpatch, ONLY: t_patch
  
#ifdef RFLU
  USE rflu_modconvertcv, ONLY: rflu_convertcvcons2prim, & 
                               rflu_convertcvprim2cons  
  USE rflu_moddifferentiationbfaces
  USE rflu_moddifferentiationfaces
  USE rflu_modviscousflux  
  USE modinterfaces, ONLY : rflu_decideneedbgradface 
#endif  
  
#ifdef RFLO
  USE modinterfaces, ONLY : rflo_calcgradvector,rflo_viscousflux
#endif
#ifdef TURB
  USE modinterfacesturbulence, ONLY : turb_coviscousfluxes
#endif  
  IMPLICIT NONE

! ******************************************************************************
! Declarations and definitions
! ******************************************************************************  

! ==============================================================================
! Arguments
! ==============================================================================

  TYPE(t_region) :: region

! ==============================================================================
! Locals
! ==============================================================================

  INTEGER :: ilev, turbmodel
#ifdef RFLU
  INTEGER :: ipatch
  INTEGER :: varinfo(cv_mixt_xvel:cv_mixt_temp)
#endif
#ifdef RFLO
  REAL(RFREAL), POINTER :: tv(:,:)
#endif
#ifdef RFLU
  TYPE(t_patch), POINTER :: ppatch
  TYPE(t_region), POINTER :: pregion
#endif
  TYPE(t_global), POINTER :: global

! ******************************************************************************
! Start
! ******************************************************************************

  global => region%global

  CALL registerfunction(global,'ViscousFluxes',&
  'ViscousFluxes.F90')

! ******************************************************************************
! Define input parameters and gridlevel
! ******************************************************************************
 
  turbmodel = region%mixtInput%turbModel 
 
#ifdef RFLO
  ilev =  region%currLevel
  tv   => region%levels(ilev)%mixt%tv
#endif

! ******************************************************************************
! Add viscous residual from laminar or laminar+turbulent stress contributions
! ******************************************************************************

#ifdef RFLO
  CALL rflo_calcgradvector( region,dv_mixt_uvel,dv_mixt_temp, &
                                   gr_mixt_ux,  gr_mixt_tz, &
                                   region%levels(ilev)%mixt%dv, &
                                   region%levels(ilev)%mixt%gradi, &
                                   region%levels(ilev)%mixt%gradj, &
                                   region%levels(ilev)%mixt%gradk )
                                   
  IF (turbmodel == turb_model_none) THEN    ! only laminar
    CALL rflo_viscousflux( region,tv_mixt_muel,tv_mixt_tcol,tv )
  ELSE
#ifdef TURB
    CALL turb_coviscousfluxes( region )       ! laminar + turbulent
#endif
  ENDIF
#endif

#ifdef RFLU
  pregion => region%pRegion

  CALL rflu_convertcvcons2prim(pregion,cv_mixt_state_duvwt)
  
  varinfo(cv_mixt_xvel) = v_mixt_xvel
  varinfo(cv_mixt_yvel) = v_mixt_yvel
  varinfo(cv_mixt_zvel) = v_mixt_zvel
  varinfo(cv_mixt_temp) = v_mixt_temp
  
  CALL rflu_computegradfaceswrapper(pregion,cv_mixt_xvel,cv_mixt_temp, & 
                                    grf_mixt_xvel,grf_mixt_temp, &
                                    pregion%mixt%cv,pregion%mixt%gradFace)
                              
  IF ( pregion%grid%nFacesConstr > 0 ) THEN                              
    CALL rflu_computegradfacesconstr(pregion,cv_mixt_xvel,cv_mixt_temp, & 
                                     grf_mixt_xvel,grf_mixt_temp,varinfo, &
                                     pregion%mixt%cv,pregion%mixt%gradFace)                              
  END IF ! pRegion%grid%nFacesConstr      

  DO ipatch = 1,pregion%grid%nPatches
    ppatch => pregion%patches(ipatch)

    IF ( rflu_decideneedbgradface(pregion,ppatch) .EQV. .true. ) THEN
    
      CALL rflu_computegradbfaceswrapper(pregion,ppatch,cv_mixt_xvel, &
                                         cv_mixt_temp,grbf_mixt_xvel, &
                                         grbf_mixt_temp,pregion%mixt%cv, &
                                         ppatch%mixt%gradFace)
    
      IF ( ppatch%cReconst /= constr_none ) THEN                                        
        CALL rflu_computebfgradconstrwrapper(pregion,ppatch,cv_mixt_xvel, &
                                             cv_mixt_temp,grbf_mixt_xvel, &
                                             grbf_mixt_temp,varinfo, &
                                             pregion%mixt%cv, &
                                             ppatch%mixt%gradFace)                           
      END IF ! pPatch%cReconst                                           
    END IF ! RFLU_DecideNeedBGradFace
  END DO ! iPatch
                                                                   
  IF ( turbmodel == turb_model_none ) THEN ! Only laminar                                    
    CALL rflu_enforceheatflux(pregion,pregion%mixt%tv,tv_mixt_tcol)  
    CALL rflu_viscousfluxes(pregion,pregion%mixt%tv,tv_mixt_muel,tv_mixt_tcol)
    CALL rflu_viscousfluxespatches(pregion,pregion%mixt%tv,tv_mixt_muel, &
                                   tv_mixt_tcol)
  ELSE 
#ifdef TURB
    CALL turb_coviscousfluxes( pregion )   ! laminar + turbulent
#endif
  END IF ! turbModel                                                                
                                                              
  CALL rflu_convertcvprim2cons(pregion,cv_mixt_state_cons)
#endif  

! ******************************************************************************
! End
! ****************************************************************************** 

  CALL deregisterfunction(region%global)

END SUBROUTINE viscousfluxes

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: ViscousFluxes.F90,v $
! Revision 1.14  2008/12/06 08:44:10  mtcampbe
! Updated license.
!
! Revision 1.13  2008/11/19 22:17:24  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.12  2006/08/19 15:38:37  mparmar
! Moved region%mixt%bGradFace to patch%mixt%gradFace
!
! Revision 1.11  2006/04/15 16:53:35  haselbac
! Added IF on cReconst flag on patch
!
! Revision 1.10  2006/04/07 15:19:15  haselbac
! Removed tabs
!
! Revision 1.9  2006/04/07 14:39:22  haselbac
! Changed calls to bface grad routines, now inside patch loop
!
! Revision 1.8  2005/10/27 18:53:54  haselbac
! Adapted to new cell grad routines
!
! Revision 1.7  2005/10/14 14:02:47  haselbac
! Added call to RFLU_EnforceHeatFlux
!
! Revision 1.6  2005/10/05 13:49:57  haselbac
! Adapted to new face grads, cosmetics
!
! Revision 1.5  2005/05/16 20:39:51  haselbac
! Now USE RFLU_ModViscousFlux
!
! Revision 1.4  2005/03/09 14:52:04  haselbac
! Cosmetics only
!
! Revision 1.3  2005/01/05 01:45:44  haselbac
! Added better comments
!
! Revision 1.2  2005/01/05 01:43:31  haselbac
! Added init of global coeffs, now also defined for all regions
!
! Revision 1.1  2004/12/01 16:52:17  haselbac
! Initial revision after changing case
!
! Revision 1.30  2004/08/02 23:13:14  wasistho
! mv libfloflu/viscousFluxEddy(Patch) to rocflo/RFLO_viscousFlux(Patch)
!
! Revision 1.29  2004/07/03 01:26:59  wasistho
! activate turbulent viscous fluxes in Rocflu
!
! Revision 1.28  2004/03/13 03:09:52  wasistho
! TURB_CoViscousFluxesFlo to TURB_CoViscousFluxes
!
! Revision 1.27  2004/03/08 23:32:41  wasistho
! changed turb file/routine name
!
! Revision 1.26  2004/01/29 22:52:50  haselbac
! Added ONLY to USE module statement
!
! Revision 1.25  2003/12/04 03:23:10  haselbac
! Added RFLU specific routines
!
! Revision 1.24  2003/11/20 16:40:36  mdbrandy
! Backing out RocfluidMP changes from 11-17-03
!
! Revision 1.21  2003/10/03 20:49:45  haselbac
! Renamed ip variable to iPatch
!
! Revision 1.20  2003/08/28 20:32:33  wasistho
! excluced ModInterfacesTurbulence,Radiation,Periodic from ModInterfaces
!
! Revision 1.19  2003/05/15 02:57:02  jblazek
! Inlined index function.
!
! Revision 1.18  2003/04/10 01:22:41  jblazek
! Got rid of pRegion in ViscousFluxesMP.
!
! Revision 1.17  2002/12/12 02:54:40  wasistho
! Added TURB_ViscousFluxes interface
!
! Revision 1.16  2002/09/09 14:09:08  haselbac
! mixtInput now under regions, added gradient routines
!
! Revision 1.15  2002/09/05 17:40:20  jblazek
! Variable global moved into regions().
!
! Revision 1.13  2002/08/27 23:58:20  wasistho
! TURB_TotalVisFluxes is set back to TURB_ViscousFluxes
!
! Revision 1.12  2002/08/16 16:51:49  wasistho
! TURB_ViscousFluxes renamed to TURB_TotalVisFluxes
!
! Revision 1.11  2002/08/15 19:48:05  jblazek
! Implemented grid deformation capability.
!
! Revision 1.10  2002/08/01 01:34:04  wasistho
! Included RFLU capability and changed CalcGradVelT
!
! Revision 1.9  2002/07/29 17:13:08  jblazek
! Clean up after RFLU and TURB.
!
! Revision 1.8  2002/07/27 08:13:16  wasistho
! prepared for rocturb implementation
!
! Revision 1.7  2002/07/23 01:08:33  wasistho
! Turbulent part of viscous flux moved to rocturb
!
! Revision 1.6  2002/07/22 22:59:11  jblazek
! Some more clean up.
!
! Revision 1.5  2002/07/19 23:38:24  wasistho
! made compliant with CODING RULE
!
! Revision 1.4  2002/06/14 21:17:01  wasistho
! Added #ifdef RFLO until RFLU ready
!
! Revision 1.3  2002/05/21 01:51:56  wasistho
! add viscous terms
!
! Revision 1.1  2002/01/23 03:51:24  jblazek
! Added low-level time-stepping routines.
!
! ******************************************************************************