ConvectiveFluxesMP.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: calculate convective fluxes for RocfluidMP framework.
!
! Description: none.
!
! Input: region = data of current region.
!
! Output: regions%levels%mixt = new solution after one time step.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: ConvectiveFluxesMP.F90,v 1.5 2008/12/06 08:44:08 mtcampbe Exp $
!
! Copyright: (c) 2003-2004 by the University of Illinois
!
!******************************************************************************

SUBROUTINE convectivefluxesmp( region )

  USE moddatatypes
  USE moddatastruct, ONLY: t_region
  USE modbndpatch, ONLY: t_patch
  USE modglobal, ONLY: t_global
  USE moderror
  USE modparameters
  USE modinterfaces, ONLY: convectivefluxes

#ifdef PEUL
  USE modinterfaceseulerian, ONLY : peul_convectivefluxes
#endif
#ifdef SPEC
#ifdef RFLU
  USE rflu_modconvertcv, ONLY: rflu_scalarconvertcvcons2prim, &
                               rflu_scalarconvertcvprim2cons
  USE modinterfaces, ONLY: rflu_scalarfirst, &
                           rflu_scalarsecond, &
                           rflu_scalarfirstpatch, &
                           rflu_scalarsecondpatch, &
                           rflu_scalarinitrhs
#endif
#endif
#ifdef TURB
  USE modinterfacesturbulence, ONLY : turb_ransconvectivefluxes
#endif

  IMPLICIT NONE

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

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

  TYPE(t_region) :: region

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

  INTEGER :: ipatch,spacediscr,spaceorder
  TYPE(t_global), POINTER :: global

#ifdef RFLU
  TYPE(t_patch), POINTER :: ppatch
  TYPE(t_region), POINTER :: pregion
#endif

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

  global => region%global

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

! *****************************************************************************
! Get parameters
! *****************************************************************************

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

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

! =============================================================================
! Mixture
! =============================================================================

  CALL convectivefluxes( region )

#ifdef PEUL
! =============================================================================
! Smoke
! =============================================================================

  IF ( global%peulUsed ) CALL peul_convectivefluxes(region)
#endif

#ifdef SPEC
! =============================================================================
! Species
! =============================================================================

  IF ( global%specUsed .EQV. .true. ) THEN
#ifdef RFLU
    pregion => region%pRegion

    CALL rflu_scalarconvertcvcons2prim(pregion,pregion%spec%cv, &
                                       pregion%spec%cvState)

    CALL rflu_scalarinitrhs(pregion,pregion%specInput%nSpecies, &
                            pregion%spec%diss,pregion%spec%rhs)

    SELECT CASE ( spaceorder ) 
      CASE ( 1 ) 
        CALL rflu_scalarfirst(pregion,pregion%specInput%nSpecies, &
                              pregion%spec%cv,pregion%spec%rhs)

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

          CALL rflu_scalarfirstpatch(pregion,ppatch, &
                                     pregion%specInput%nSpecies, &
                                     pregion%spec%cv,ppatch%spec, &
                                     pregion%spec%rhs)
        END DO ! iPatch                                     
      CASE ( 2 )     
        CALL rflu_scalarsecond(pregion,pregion%specInput%nSpecies, &
                               pregion%spec%cv,pregion%spec%gradCell, &
                               pregion%spec%rhs)

        DO ipatch = 1,pregion%grid%nPatches
          ppatch => pregion%patches(ipatch)
          
          SELECT CASE ( ppatch%spaceOrder ) 
            CASE ( 1 ) 
              CALL rflu_scalarfirstpatch(pregion,ppatch, &
                                         pregion%specInput%nSpecies, &
                                         pregion%spec%cv,ppatch%spec, &
                                         pregion%spec%rhs)            
            CASE ( 2 )           
              CALL rflu_scalarsecondpatch(pregion,ppatch, &
                                          pregion%specInput%nSpecies, &
                                          pregion%spec%cv, &
                                          pregion%spec%gradCell, &
                                          ppatch%spec,pregion%spec%rhs)
            CASE default
              CALL errorstop(global,err_reached_default,__line__)
          END SELECT ! pPatch%spaceOrder
        END DO ! iPatch        
      CASE default 
        CALL errorstop(global,err_reached_default,__line__)
    END SELECT ! spaceOrder

    CALL rflu_scalarconvertcvprim2cons(pregion,pregion%spec%cv, &
                                       pregion%spec%cvState)
#endif
  END IF ! global%specUsed
#endif

#ifdef TURB
! =============================================================================
! Turbulence
! =============================================================================

  IF ((region%mixtInput%flowModel == flow_navst) .AND. &
      (region%mixtInput%turbModel /= turb_model_none)) THEN
    CALL turb_ransconvectivefluxes( region )
  ENDIF
#endif

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

  CALL deregisterfunction( global )

END SUBROUTINE convectivefluxesmp

!******************************************************************************
!
! RCS Revision history:
!
! $Log: ConvectiveFluxesMP.F90,v $
! Revision 1.5  2008/12/06 08:44:08  mtcampbe
! Updated license.
!
! Revision 1.4  2008/11/19 22:17:22  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.3  2006/08/19 15:38:28  mparmar
! Renamed patch variables
!
! Revision 1.2  2006/04/15 16:53:04  haselbac
! Added option of 1st order bfluxes with 2nd order scheme
!
! Revision 1.1  2004/12/01 16:48:27  haselbac
! Initial revision after changing case
!
! Revision 1.12  2004/07/28 15:29:18  jferry
! created global variable for spec use
!
! Revision 1.11  2004/03/05 22:09:00  jferry
! created global variables for peul, plag, and inrt use
!
! Revision 1.10  2004/01/29 22:52:40  haselbac
! Added second-order for species
!
! Revision 1.9  2003/11/25 21:01:38  haselbac
! Added rocspecies support
!
! Revision 1.8  2003/11/20 16:40:35  mdbrandy
! Backing out RocfluidMP changes from 11-17-03
!
! Revision 1.5  2003/10/03 20:11:55  wasistho
! initial installation of turbModel SA and DES
!
! Revision 1.4  2003/10/01 23:52:09  jblazek
! Corrected bug in moving noslip wall BC and grid speeds.
!
! Revision 1.3  2003/09/26 21:45:34  fnajjar
! Modified ModInterfaces calls to proper physical modules
!
! Revision 1.2  2003/05/15 02:57:02  jblazek
! Inlined index function.
!
! Revision 1.1  2003/03/28 19:45:46  fnajjar
! Initial import for RocfluidMP
!
!******************************************************************************