Slip wall boundary condition. Maintain flow tangency
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
character(len=*), | intent(in) | :: | face | |||
type(facetype), | intent(in), | dimension(-2:dims%imx+3,-2:dims%jmx+2,-2:dims%kmx+2) | :: | Ifaces | Input varaible which stores I faces' area and unit normal |
|
type(facetype), | intent(in), | dimension(-2:dims%imx+2,-2:dims%jmx+3,-2:dims%kmx+2) | :: | Jfaces | Input varaible which stores J faces' area and unit normal |
|
type(facetype), | intent(in), | dimension(-2:dims%imx+2,-2:dims%jmx+2,-2:dims%kmx+3) | :: | Kfaces | Input varaible which stores K faces' area and unit normal Name of the face at which boundary condition is called |
|
type(boundarytype), | intent(in) | :: | bc | |||
type(extent), | intent(in) | :: | dims |
subroutine slip_wall(face, Ifaces, Jfaces, Kfaces, bc, dims)
!< Slip wall boundary condition.
!< Maintain flow tangency
implicit none
type(extent), intent(in) :: dims
type(boundarytype), intent(in) :: bc
character(len=*), intent(in) :: face
type(facetype), dimension(-2:dims%imx+3,-2:dims%jmx+2,-2:dims%kmx+2), intent(in) :: Ifaces
!< Input varaible which stores I faces' area and unit normal
type(facetype), dimension(-2:dims%imx+2,-2:dims%jmx+3,-2:dims%kmx+2), intent(in) :: Jfaces
!< Input varaible which stores J faces' area and unit normal
type(facetype), dimension(-2:dims%imx+2,-2:dims%jmx+2,-2:dims%kmx+3), intent(in) :: Kfaces
!< Input varaible which stores K faces' area and unit normal
!< Name of the face at which boundary condition is called
call copy3(density , "symm", face, bc, dims)
call copy3(pressure, "symm", face, bc, dims)
select case (turbulence)
case('none')
!do nothing
continue
case('sa', 'saBC')
call copy3(tv, "symm", face, bc, dims)
case('sst', 'sst2003')
call copy3(tk, "symm", face, bc, dims)
call copy3(tw, "symm", face, bc, dims)
case('kkl')
call copy3(tk, "symm", face, bc, dims)
call copy3(tkl, "symm", 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
call flow_tangency(qp, face, Ifaces, Jfaces, Kfaces, dims)
end subroutine slip_wall