Supersonic outlet boundary condition. All the values of state variables are copied from inside the domain
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
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 |
subroutine supersonic_outlet(face, bc, dims)
!< Supersonic outlet boundary condition.
!< All the values of state variables are 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
call copy3(density , "flat", face, bc, dims)
call copy3(x_speed , "flat", face, bc, dims)
call copy3(y_speed , "flat", face, bc, dims)
call copy3(z_speed , "flat", face, bc, dims)
call copy3(pressure, "flat", face, bc, dims)
select case (turbulence)
case('none')
!do nothing
continue
case('sa', 'saBC')
call copy3(tv, "flat", face, bc, dims)
case('sst', 'sst2003')
call copy3(tk, "flat", face, bc, dims)
call copy3(tw, "flat", face, bc, dims)
case('kkl')
call copy3(tk, "flat", face, bc, dims)
call copy3(tkl, "flat", face, bc, dims)
case DEFAULT
Fatal_error
end select
select case(trim(transition))
case('lctm2015')
call copy3(tgm, "flat", face, bc, dims)
case DEFAULT
continue
end select
end subroutine supersonic_outlet