TURB_fluLesBLij.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 the scale similarity turbulent stress tensor, which
!          is the same as the Leonard stress term in the dynamic LES
!          models. SS-tauij uses bar-filter while Leonard stress uses
!          test-filter.
!
! Description: The turbulent stress tensor follows the SS model and can
!              be seen as Leonard tensor if the outer bar filter is replaced
!              by test filter. This tensor is defined by
!              l_{ij}=bar(bar[rho*ui]bar[rho*uj]/bar[rho])
!                    -bar(bar[rho*ui])bar(bar[rho*uj])/bar(bar[rho])
!              where bar denotes the bar-filter with filter-width defined by
!              nDel(DIRI:DIRK). On entry we have the solution components 
!              at cell faces and perform the filtering from face to face. 
!              We return the resulting tensor at all cell faces including  
!              dummies. Also 1/bar(bar[rho]) is available stored ffVar(1,:).
!
! Input: region = data of current region
!        nDel   = three components filter width paramater
!
! Output: lij = resulting stress tensor of scale similarity or Leonard term
!
! Notes: This routine relevant only for physical boundaries of FLU regions.
!  1. Removed CALL to RFLU_InterpCells2FacesPatches because it no longer exists 
!     This routine might become broken because of this
!
!******************************************************************************
!
! $Id: TURB_fluLesBLij.F90,v 1.5 2008/12/06 08:44:44 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE turb_flulesblij( region,nDel,lij )

  USE moddatatypes
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modbndpatch, ONLY   : t_patch
  USE moderror
  USE turb_modparameters
  IMPLICIT NONE

! ... parameters
  TYPE(t_region), POINTER :: region
  INTEGER                 :: ndel(diri:dirk)
  REAL(RFREAL), POINTER   :: lij(:,:)

! ... loop variables
  INTEGER :: ijkn, ipatch

! ... local variables
  CHARACTER(CHRLEN)       :: rcsidentstring
  TYPE(t_global), POINTER :: global
  TYPE(t_patch), POINTER  :: patch

  INTEGER :: ibn,ien,idbeg,idend,npatches,nbfaces
  REAL(RFREAL), POINTER :: bfvar(:,:),bffvar(:,:)
  REAL(RFREAL), POINTER :: tens(:,:),tensbar(:,:)
  REAL(RFREAL), POINTER :: cv(:,:)

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

  rcsidentstring = '$RCSfile: TURB_fluLesBLij.F90,v $'

  global => region%global
  CALL registerfunction( global,'TURB_FluLesBLij',&
  'TURB_fluLesBLij.F90' )

! get indices and pointers --------------------------------------------------

  npatches = region%grid%nPatches
  nbfaces  = 0

  DO ipatch = 1,npatches
    patch   => region%patches(ipatch)
    nbfaces =  nbfaces + patch%nBTris + patch%nBQuads
  END DO ! iPatch
  ibn = 1
  ien = nbfaces

  cv     => region%mixt%cv
  bfvar  => region%turb%bfVar
  bffvar => region%turb%bffVar

! allocate temporary arrays

  ALLOCATE( tens(e11:e33,ibn:ien),tensbar(e11:e33,ibn:ien) )

! ff-filtering bRhoF and bRhoUiF (i=1,3) stored in fVar to get 
! bar[bRhoF] and bar[bRhoUiF] stored in fbVar

  idbeg = cv_turb_dens
  idend = cv_turb_zmom

! for efficiency we store 1/bRhoF in bRhoF and 1/bbRhoF in bbRhoF;
! next, build components of scale-similarity (Bardina) turbulent stress 
! tensor lij-component: store bRhoUi*bRhoUj/bRho in tens

  DO ijkn = ibn,ien
    bfvar(cv_turb_dens,ijkn)=1._rfreal/bfvar(cv_turb_dens,ijkn)
    bffvar(cv_turb_dens,ijkn)=1._rfreal/bffvar(cv_turb_dens,ijkn)

    tens(e11,ijkn) = bfvar(cv_turb_xmom,ijkn)*bfvar(cv_turb_xmom,ijkn)* &
                     bfvar(cv_turb_dens,ijkn)
    tens(e12,ijkn) = bfvar(cv_turb_xmom,ijkn)*bfvar(cv_turb_ymom,ijkn)* &
                     bfvar(cv_turb_dens,ijkn)
    tens(e13,ijkn) = bfvar(cv_turb_xmom,ijkn)*bfvar(cv_turb_zmom,ijkn)* &
                     bfvar(cv_turb_dens,ijkn)
    tens(e22,ijkn) = bfvar(cv_turb_ymom,ijkn)*bfvar(cv_turb_ymom,ijkn)* &
                     bfvar(cv_turb_dens,ijkn)
    tens(e23,ijkn) = bfvar(cv_turb_ymom,ijkn)*bfvar(cv_turb_zmom,ijkn)* &
                     bfvar(cv_turb_dens,ijkn)
    tens(e33,ijkn) = bfvar(cv_turb_zmom,ijkn)*bfvar(cv_turb_zmom,ijkn)* &
                     bfvar(cv_turb_dens,ijkn)
  ENDDO

! ff-filter bRhoUi*bRhoUj/bRho stored in tens to get 
! bar[bRhoUi*bRhoUj/bRho] in tensBar

  idbeg = e11
  idend = e33

!  CALL RFLU_InterpFaces2FacesPatches( region,tens,tensBar )

! combine and find lij=bar[bRhoUi*bRhoUj/bRho]-bbRhoUi*bbRhoUj/bbRho

  DO ijkn=ibn,ien
    lij(e11,ijkn)=tensbar(e11,ijkn)-bffvar(cv_turb_dens,ijkn)* &
                  bffvar(cv_turb_xmom,ijkn)*bffvar(cv_turb_xmom,ijkn)

    lij(e12,ijkn)=tensbar(e12,ijkn)-bffvar(cv_turb_dens,ijkn)* &
                  bffvar(cv_turb_xmom,ijkn)*bffvar(cv_turb_ymom,ijkn)

    lij(e13,ijkn)=tensbar(e13,ijkn)-bffvar(cv_turb_dens,ijkn)* &
                  bffvar(cv_turb_xmom,ijkn)*bffvar(cv_turb_zmom,ijkn)

    lij(e22,ijkn)=tensbar(e22,ijkn)-bffvar(cv_turb_dens,ijkn)* &
                  bffvar(cv_turb_ymom,ijkn)*bffvar(cv_turb_ymom,ijkn)

    lij(e23,ijkn)=tensbar(e23,ijkn)-bffvar(cv_turb_dens,ijkn)* &
                  bffvar(cv_turb_ymom,ijkn)*bffvar(cv_turb_zmom,ijkn)

    lij(e33,ijkn)=tensbar(e33,ijkn)-bffvar(cv_turb_dens,ijkn)* &
                  bffvar(cv_turb_zmom,ijkn)*bffvar(cv_turb_zmom,ijkn)
  ENDDO

! deallocate temporary arrays

  DEALLOCATE( tens,tensbar )

! finalize --------------------------------------------------------------------

  CALL deregisterfunction( global )

END SUBROUTINE turb_flulesblij

!******************************************************************************
!
! RCS Revision history:
!
! $Log: TURB_fluLesBLij.F90,v $
! Revision 1.5  2008/12/06 08:44:44  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:56  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2006/08/19 15:41:03  mparmar
! Removed call to RFLU_InterpCells2FacesPatches
!
! Revision 1.2  2004/08/07 01:07:18  wasistho
! bugfixed, defined pointer cv
!
! Revision 1.1  2004/05/28 02:03:58  wasistho
! update unstructured grid LES
!
!
!
!******************************************************************************