subsonic_inlet Subroutine

  1. bc_primitive.f90
  2. bc_primitive
  3. subsonic_inlet

private subroutine subsonic_inlet(face, bc, dims)

Subsonic inlet boundary condition. All the state variables's value expect pressure is fixed and pressure is copied from inside the domain

Arguments

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

Name of the face at which boundary condition is called
type(boundarytype), intent(in) :: bc
type(extent), intent(in) :: dims

Calls

proc~~subsonic_inlet~~CallsGraph proc~subsonic_inlet subsonic_inlet proc~copy3 copy3 proc~subsonic_inlet->proc~copy3 proc~fix fix proc~subsonic_inlet->proc~fix

Called by

proc~~subsonic_inlet~~CalledByGraph proc~subsonic_inlet subsonic_inlet proc~populate_ghost_primitive populate_ghost_primitive proc~populate_ghost_primitive->proc~subsonic_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

subsonic_inlet

Source Code

      subroutine subsonic_inlet(face, bc, dims)
        !< Subsonic inlet boundary condition. 
        !< All the state variables's value expect pressure
        !< is fixed and pressure is copied from inside the 
        !< domain
        implicit none
        type(extent), intent(in) :: dims
        type(boundarytype), intent(in) :: bc
        character(len=*), intent(in) :: face
        !< 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)
        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_tw, 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
        call copy3(pressure, "flat", face, bc, dims)
      end subroutine subsonic_inlet