ConvectiveFluxes.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 convective fluxes, add them to dissipation in order
!          to obtain the residual.
!
! Description: none.
!
! Input: region = data of current region.
!
! Output: region%levels%mixt%rhs = convective fluxes + dissipation.
!
! Notes: none.
!
! ******************************************************************************
!
! $Id: ConvectiveFluxes.F90,v 1.10 2008/12/06 08:44:08 mtcampbe Exp $
!
! Copyright: (c) 2001-2006 by the University of Illinois
!
! ******************************************************************************

SUBROUTINE convectivefluxes( region )

  USE moddatatypes
  USE modglobal, ONLY: t_global
  USE moddatastruct, ONLY : t_region
  USE moderror
  USE modparameters

#ifdef RFLU
  USE rflu_modnscbc

  USE modinterfaces, ONLY: rflu_computefluxinv
#endif


#ifdef RFLO
  USE modinterfaces, ONLY : rflo_getdimensdummy, rflo_getcelloffset, &
        rflo_centralflux, rflo_roefluxfirst, rflo_roefluxsecond
#include "Indexing.h"
#endif

  IMPLICIT NONE

! ******************************************************************************
! Definitions and declarations
! ******************************************************************************

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

  TYPE(t_region), TARGET :: region

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

#ifdef RFLO
  INTEGER :: ibc, ic, idcbeg, idcend, iec, jdcbeg, jdcend, kdcbeg, kdcend
  INTEGER :: ilev, icoff, ijcoff
#endif
#ifdef RFLU
  INTEGER :: icg
#endif
  INTEGER :: spacediscr, spaceorder
  REAL(RFREAL), POINTER :: diss(:,:), rhs(:,:)
  TYPE(t_global), POINTER :: global
#ifdef RFLU
  TYPE(t_region), POINTER :: pregion
#endif

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

  global => region%global

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

! ******************************************************************************
! Get dimensions and pointers
! ******************************************************************************

#ifdef RFLO
  ilev = region%currLevel

  CALL rflo_getdimensdummy( region,ilev,idcbeg,idcend,jdcbeg,jdcend, &
                            kdcbeg,kdcend )
  CALL rflo_getcelloffset( region,ilev,icoff,ijcoff )
  ibc = indijk(idcbeg,jdcbeg,kdcbeg,icoff,ijcoff)
  iec = indijk(idcend,jdcend,kdcend,icoff,ijcoff)

  diss => region%levels(ilev)%mixt%diss
  rhs  => region%levels(ilev)%mixt%rhs
#endif
#ifdef RFLU
  pregion => region
  
  diss => region%mixt%diss
  rhs  => region%mixt%rhs  
#endif

  spacediscr = region%mixtInput%spaceDiscr
  spaceorder = region%mixtInput%spaceOrder

#ifdef RFLO
! ******************************************************************************
! Initialize residual (set = dissipation)
! ******************************************************************************

  DO ic=ibc,iec
    rhs(cv_mixt_dens,ic) = -diss(cv_mixt_dens,ic)
    rhs(cv_mixt_xmom,ic) = -diss(cv_mixt_xmom,ic)
    rhs(cv_mixt_ymom,ic) = -diss(cv_mixt_ymom,ic)
    rhs(cv_mixt_zmom,ic) = -diss(cv_mixt_zmom,ic)
    rhs(cv_mixt_ener,ic) = -diss(cv_mixt_ener,ic)
  ENDDO

! ******************************************************************************
! Compute fluxes
! ******************************************************************************

! ==============================================================================
! 2nd-order central scheme
! ==============================================================================

  IF (spacediscr==discr_cen_scal .AND. &
      (spaceorder==discr_order_1 .OR. spaceorder==discr_order_2)) THEN
    CALL rflo_centralflux( region )
  ENDIF

! ==============================================================================
! Roe upwind scheme 
! ==============================================================================

  IF (spacediscr == discr_upw_roe) THEN
    IF (spaceorder == discr_order_1) THEN
      CALL rflo_roefluxfirst( region )
    ELSE IF (spaceorder == discr_order_2) THEN
      CALL rflo_roefluxsecond( region )
    ENDIF
  ENDIF
#endif

#ifdef RFLU
! ******************************************************************************
! Initialize residual (set = dissipation). NOTE only need to do this for 
! unsteady flows because for steady flows this is now done by a dedicated 
! routine (RFLU_EMS_SetRhs).
! ******************************************************************************

  IF ( global%flowType == flow_unsteady ) THEN ! Unsteady
    DO icg = 1,region%grid%nCellsTot
      rhs(cv_mixt_dens,icg) = -diss(cv_mixt_dens,icg)
      rhs(cv_mixt_xmom,icg) = -diss(cv_mixt_xmom,icg)
      rhs(cv_mixt_ymom,icg) = -diss(cv_mixt_ymom,icg)
      rhs(cv_mixt_zmom,icg) = -diss(cv_mixt_zmom,icg)
      rhs(cv_mixt_ener,icg) = -diss(cv_mixt_ener,icg)
    END DO ! icg
  END IF ! global%flowType

! ******************************************************************************
! Compute fluxes
! ******************************************************************************

! ==============================================================================
! Check state of conserved variable vector
! ==============================================================================

  IF ( pregion%mixt%cvState /= cv_mixt_state_cons ) THEN 
    CALL errorstop(pregion%global,err_cv_state_invalid,__line__)
  END IF ! pRegion%mixt%cvState

! ==============================================================================
! Call flux functions 
! ==============================================================================

  CALL rflu_computefluxinv(pregion,flux_part_both)

! ******************************************************************************
! Compute Boundary Patch Rhs 
! ******************************************************************************

  IF ( rflu_nscbc_decidehavenscbc(pregion) .EQV. .true. ) THEN
    CALL rflu_nscbc_comprhs(pregion)
  END IF ! RFLU_NSCBC_DecideHaveNSCBC(pRegion)
#endif

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

  CALL deregisterfunction( region%global )

END SUBROUTINE convectivefluxes

! ******************************************************************************
!
! RCS Revision history:
!
! $Log: ConvectiveFluxes.F90,v $
! Revision 1.10  2008/12/06 08:44:08  mtcampbe
! Updated license.
!
! Revision 1.9  2008/11/19 22:17:22  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.8  2006/08/19 15:38:26  mparmar
! Added computing of fluxes for boundary variables
!
! Revision 1.7  2006/05/01 20:57:43  haselbac
! Converted to new flux routine
!
! Revision 1.6  2006/03/26 20:21:12  haselbac
! Rewrite bcos of of GL model
!
! Revision 1.5  2005/11/14 16:53:37  haselbac
! Added AUSM flux computation for pseudo-gas
!
! Revision 1.4  2005/11/10 01:55:16  haselbac
! Added new logic for AUSM scheme
!
! Revision 1.3  2005/07/14 21:39:28  haselbac
! Added AUSM flux function
!
! Revision 1.2  2005/05/16 20:37:56  haselbac
! Changed args to pRegion, rmvd OLES call, changed init of rhs, cosmetics
!
! Revision 1.1  2004/12/01 16:48:25  haselbac
! Initial revision after changing case
!
! Revision 1.17  2004/02/13 02:56:23  haselbac
! Added calls (commented out for now) to fast Roe routines
!
! Revision 1.16  2004/01/29 22:52:38  haselbac
! Rewrote for positive species update
!
! Revision 1.15  2003/12/04 03:22:57  haselbac
! Added second-order schemes
!
! Revision 1.14  2003/11/20 16:40:35  mdbrandy
! Backing out RocfluidMP changes from 11-17-03
!
! Revision 1.11  2003/05/29 17:28:42  jblazek
! Implemented Roe scheme.
!
! Revision 1.10  2003/05/16 22:05:40  haselbac
! Added HLLC flux
!
! Revision 1.9  2003/05/15 02:57:02  jblazek
! Inlined index function.
!
! Revision 1.8  2002/09/09 13:56:51  haselbac
! mixtInput under region, bug fix for iec, added state check
!
! Revision 1.7  2002/09/05 17:40:20  jblazek
! Variable global moved into regions().
!
! Revision 1.6  2002/07/29 17:13:08  jblazek
! Clean up after RFLU and TURB.
!
! Revision 1.5  2002/07/25 14:46:49  haselbac
! Added optimal LES flux routine
!
! Revision 1.4  2002/05/04 16:20:24  haselbac
! Added RFLU statements
!
! Revision 1.3  2002/02/21 23:25:05  jblazek
! Blocks renamed as regions.
!
! Revision 1.2  2002/01/28 23:55:22  jblazek
! Added flux computation (central scheme).
!
! Revision 1.1  2002/01/23 03:51:24  jblazek
! Added low-level time-stepping routines.
!
! ******************************************************************************