supersonic_inlet Subroutine

private subroutine supersonic_inlet(face, bc)

Supersonic inlet boundary condition All the values of state variables are fixed

Arguments

Type IntentOptional AttributesName
character(len=*), intent(in) :: face
type(boundarytype), intent(in) :: bc

Name of the face at which boundary condition is called


Calls

proc~~supersonic_inlet~~CallsGraph proc~supersonic_inlet supersonic_inlet proc~fix fix proc~supersonic_inlet->proc~fix

Called by

proc~~supersonic_inlet~~CalledByGraph proc~supersonic_inlet supersonic_inlet proc~populate_ghost_primitive populate_ghost_primitive proc~populate_ghost_primitive->proc~supersonic_inlet proc~get_total_conservative_residue get_total_conservative_Residue proc~get_total_conservative_residue->proc~populate_ghost_primitive 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


Source Code

      subroutine supersonic_inlet(face, bc)
        !< Supersonic inlet boundary condition
        !< All the values of state variables are fixed
        implicit none
        character(len=*), intent(in) :: face
        type(boundarytype), intent(in) :: bc
        !< Name of the face at which boundary condition is called
        if(current_iter<=2)then
        call fix(density , bc%fixed_density , face)
        call fix(x_speed , bc%fixed_x_speed , face)
        call fix(y_speed , bc%fixed_y_speed , face)
        call fix(z_speed , bc%fixed_z_speed , face)
        call fix(pressure, bc%fixed_pressure, face)
        select case (turbulence)
          case('none')
            !do nothing
            continue
          case('sa', 'saBC')
            call fix(tv, bc%fixed_tv, face)
          case('sst', 'sst2003')
            !call check_if_value_fixed(bc, "sst")
            call fix(tk, bc%fixed_tk, face)
            call fix(tw, bc%fixed_tw, face)
          case('kkl')
            !call check_if_value_fixed(bc, "kkl")
            call fix(tk, bc%fixed_tk, face)
            call fix(tkl, bc%fixed_tkl, face)
          case DEFAULT
            Fatal_error
        end select
        select case(trim(transition))
          case('lctm2015')
            !call check_if_value_fixed(bc, "lctm2015")
            call fix(tgm, bc%fixed_tgm, face)
          case DEFAULT
            continue
        end select
        end if
      end subroutine supersonic_inlet