TURB_WlmUpdateBndlay.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: Model wall stresses based on turbulent BL assumption with mixing
!          length model for eddy viscosity, and pressure gradient as well as
!          usteady terms taken into account.
!
! Description: Wall stresses are derived from modeling terms in BL equations.
!              Viscous term has been modeled in routine TURB_FloWlmUpdateLoglay
!              based on mixing lenth eddy viscosity model which equivalent to
!              log layer assumption at zero pressure gradient with surface
!              roughness included. BL convective term is neglected. Pressure 
!              gradient and time derivative term are computed here and added 
!              to the viscous model. The method is inspired by paper of
!              Hoffman and Benocci, "Approximate wall BC for LES", 5th
!              Advances in Turbulence, Siena, Italy 1994. Surface roughness
!              were not considered in their model.
!
! Input: region  = data of current region.
!        patch   = current patch.
!
! Output: total and wall parallel components of wall stresses in body fitted 
!         coordinate.
!
! Notes: Unsteady term is not modeled yet. Its model term may be included 
!        in future.
!
!******************************************************************************
!
! $Id: TURB_WlmUpdateBndlay.F90,v 1.5 2008/12/06 08:44:43 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE turb_wlmupdatebndlay( region,patch )

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

#ifdef RFLO
#include "Indexing.h"
#endif

! ... parameters
  TYPE(t_region) :: region
  TYPE(t_patch)  :: patch

! ... loop variables
  INTEGER :: ijkval

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

  INTEGER :: ijbeg, ijend
  REAL(RFREAL)          :: wdist, dpdxi, dpdzt, utau, abvel
  REAL(RFREAL), POINTER :: vals(:,:)

#ifdef RFLO
  INTEGER :: n1, n2, ioff
#endif

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

  rcsidentstring = '$RCSfile: TURB_WlmUpdateBndlay.F90,v $ $Revision: 1.5 $'

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

! get pointers ----------------------------------------------------------------

  vals  => patch%valBola%vals 

! get dimensions

#ifdef RFLO
  n1    = abs(patch%l1end-patch%l1beg)
  n2    = abs(patch%l2end-patch%l2beg)
  ioff  = n1 + 1
  ijbeg = indij( 0, 0,ioff)
  ijend = indij(n1,n2,ioff)
#endif
#ifdef RFLU
  ijbeg = 1
  ijend = patch%nBFaces
#endif

! compute tau-wall and heat flux from total contributions

  DO ijkval=ijbeg,ijend

    wdist = vals(ijkval,wlm_vals_wdist)
    dpdxi = vals(ijkval,wlm_vals_dpdxi)
    dpdzt = vals(ijkval,wlm_vals_dpdzt)

    vals(ijkval,wlm_vals_tauuy) = vals(ijkval,wlm_vals_tauuy)-wdist*dpdxi     
    vals(ijkval,wlm_vals_tauwy) = vals(ijkval,wlm_vals_tauwy)-wdist*dpdzt

    vals(ijkval,wlm_vals_tauvx) = vals(ijkval,wlm_vals_tauuy)     
    vals(ijkval,wlm_vals_tauvz) = vals(ijkval,wlm_vals_tauwy)      

! - store total wall stress in heat-flux array for heat transfer comp. later
    vals(ijkval,wlm_vals_hflux) = &
          sqrt( vals(ijkval,wlm_vals_tauuy)*vals(ijkval,wlm_vals_tauuy) + & 
                vals(ijkval,wlm_vals_tauwy)*vals(ijkval,wlm_vals_tauwy) )

! - Note: utau should not be derived from here (effect of pressure gradient 
!   exist), as it is used to initiate iteration on utau at the next stage 
!   in TURB_FloWlmUpdateLoglay (model for BL viscous term), in which pressure
!   gradient is not (yet) taken into account. On the contrary, for post-
!   processing/visualisation utau has to be extracted separately from wall 
!   stress tauWall computed in this routine: utau=sqrt(tauWall/rho)
   
!checkprobe--------------------------------------------------------------------
!    utau = SQRT( vals(ijkVal,WLM_VALS_HFLUX)/vals(ijkVal,WLM_VALS_DENS) )
!    abVel= SQRT( vals(ijkVal,WLM_VALS_XIV)**2+vals(ijkVal,WLM_VALS_ZTV)**2 )
!    write(*,*) region%procId,patch%lbound,ijkVal,vals(ijkVal,WLM_VALS_TAUUY), &
!               vals(ijkVal,WLM_VALS_TAUWY),utau/abVel,dPdXi,dPdZt
!------------------------------------------------------------------------------   

  ENDDO       ! ijkVal

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

  CALL deregisterfunction( global )

END SUBROUTINE turb_wlmupdatebndlay

!******************************************************************************
!
! RCS Revision history:
! 
! $Log: TURB_WlmUpdateBndlay.F90,v $
! Revision 1.5  2008/12/06 08:44:43  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:55  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2004/03/24 03:37:03  wasistho
! prepared for RFLU
!
! Revision 1.2  2004/03/12 02:55:35  wasistho
! changed rocturb routine names
!
! Revision 1.1  2004/03/05 04:37:01  wasistho
! changed nomenclature
!
! Revision 1.3  2004/03/02 03:51:27  wasistho
! forgot colon after Id and Log
!
!
!******************************************************************************