compute_fluxes Subroutine

  1. scheme.f90
  2. scheme
  3. compute_fluxes

public subroutine compute_fluxes(F, G, H, Ifaces, Jfaces, Kfaces, scheme, flow, bc, dims)

Arguments

Type IntentOptional AttributesName
real(kind=wp), intent(inout), dimension(:, :, :, :):: F

Store fluxes throught the I faces
real(kind=wp), intent(inout), dimension(:, :, :, :):: G

Store fluxes throught the J faces
real(kind=wp), intent(inout), dimension(:, :, :, :):: H

Store fluxes throught the K faces
type(facetype), intent(in), dimension(-2:dims%imx+3,-2:dims%jmx+2,-2:dims%kmx+2):: Ifaces

Store face quantites for I faces

type(facetype), intent(in), dimension(-2:dims%imx+2,-2:dims%jmx+3,-2:dims%kmx+2):: Jfaces

Store face quantites for J faces

type(facetype), intent(in), dimension(-2:dims%imx+2,-2:dims%jmx+2,-2:dims%kmx+3):: Kfaces

Store face quantites for K faces

type(schemetype), intent(in) :: scheme

finite-volume Schemes
type(flowtype), intent(in) :: flow

Information about fluid flow: freestream-speed, ref-viscosity,etc.
type(boundarytype), intent(in) :: bc

boundary conditions and fixed values
type(extent), intent(in) :: dims

extent of the 3D domain

Calls

proc~~compute_fluxes~7~~CallsGraph proc~compute_fluxes~7 compute_fluxes compute_fluxes_van_leer compute_fluxes_van_leer proc~compute_fluxes~7->compute_fluxes_van_leer compute_fluxes_ldfss0 compute_fluxes_ldfss0 proc~compute_fluxes~7->compute_fluxes_ldfss0 compute_fluxes_ausm compute_fluxes_ausm proc~compute_fluxes~7->compute_fluxes_ausm compute_fluxes_ausmp compute_fluxes_ausmp proc~compute_fluxes~7->compute_fluxes_ausmp compute_fluxes_slau compute_fluxes_slau proc~compute_fluxes~7->compute_fluxes_slau compute_fluxes_ausmup compute_fluxes_ausmup proc~compute_fluxes~7->compute_fluxes_ausmup

Called by

proc~~compute_fluxes~7~~CalledByGraph proc~compute_fluxes~7 compute_fluxes proc~get_total_conservative_residue get_total_conservative_Residue proc~get_total_conservative_residue->proc~compute_fluxes~7 proc~get_next_solution get_next_solution proc~get_next_solution->proc~get_total_conservative_residue proc~iterate_one_more_time_step iterate_one_more_time_step proc~iterate_one_more_time_step->proc~get_next_solution program~main main program~main->proc~iterate_one_more_time_step

Contents

Source Code

compute_fluxes

Source Code

        subroutine compute_fluxes(F,G,H, Ifaces, Jfaces, Kfaces, scheme, flow, bc, dims)
        
            implicit none
            type(schemetype), intent(in) :: scheme
            !< finite-volume Schemes
            type(flowtype), intent(in) :: flow
            !< Information about fluid flow: freestream-speed, ref-viscosity,etc.
            type(extent), intent(in) :: dims
            !< extent of the 3D domain
            real(wp), dimension(:, :, :, :), intent(inout) :: F
            !< Store fluxes throught the I faces
            real(wp), dimension(:, :, :, :), intent(inout) :: G
            !< Store fluxes throught the J faces
            real(wp), dimension(:, :, :, :), intent(inout) :: H
            !< Store fluxes throught the K faces
            type(facetype), dimension(-2:dims%imx+3,-2:dims%jmx+2,-2:dims%kmx+2), intent(in) :: Ifaces
            !< Store face quantites for I faces 
            type(facetype), dimension(-2:dims%imx+2,-2:dims%jmx+3,-2:dims%kmx+2), intent(in) :: Jfaces
            !< Store face quantites for J faces 
            type(facetype), dimension(-2:dims%imx+2,-2:dims%jmx+2,-2:dims%kmx+3), intent(in) :: Kfaces
            !< Store face quantites for K faces 
            type(boundarytype), intent(in) :: bc
            !< boundary conditions and fixed values

            select case (scheme%scheme_name)
                case ("van_leer")
                  call compute_fluxes_van_leer(F,G,H,x_qp_left,x_qp_right,y_qp_left,y_qp_right,z_qp_left,z_qp_right,Ifaces,Jfaces,Kfaces,flow,bc,dims)
                case ("ausm")
                  call compute_fluxes_ausm(F,G,H,x_qp_left,x_qp_right,y_qp_left,y_qp_right,z_qp_left,z_qp_right,Ifaces,Jfaces,Kfaces,flow,bc,dims)
                case ("ausmP")
                  call compute_fluxes_ausmP(F,G,H,x_qp_left,x_qp_right,y_qp_left,y_qp_right,z_qp_left,z_qp_right,Ifaces,Jfaces,Kfaces,flow,bc,dims)
                case ("ausmUP")
                  call compute_fluxes_ausmUP(F,G,H,x_qp_left,x_qp_right,y_qp_left,y_qp_right,z_qp_left,z_qp_right,Ifaces,Jfaces,Kfaces,flow,bc,dims)
                case ("slau")
                  call compute_fluxes_slau(F,G,H,x_qp_left,x_qp_right,y_qp_left,y_qp_right,z_qp_left,z_qp_right,Ifaces,Jfaces,Kfaces,flow,bc,dims)
                case ("ldfss0")
                  call compute_fluxes_ldfss0(F,G,H,x_qp_left,x_qp_right,y_qp_left,y_qp_right,z_qp_left,z_qp_right,Ifaces,Jfaces,Kfaces,flow,bc,dims)
                case default
                    Fatal_error
            end select
            
        end subroutine compute_fluxes