Supersonic inlet boundary condition All the values of state variables are fixed
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
character(len=*), | intent(in) | :: | face | |||
type(boundarytype), | intent(in) | :: | bc | Name of the face at which boundary condition is called |
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