no_slip Subroutine

  1. bc_primitive.f90
  2. bc_primitive
  3. no_slip

private subroutine no_slip(face, bc, dims)

No-slip wall boundary condition. All the component of velocity throught face is zero

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~~no_slip~~CallsGraph proc~no_slip no_slip proc~set_omega_at_wall set_omega_at_wall proc~no_slip->proc~set_omega_at_wall proc~copy3 copy3 proc~no_slip->proc~copy3

Called by

proc~~no_slip~~CalledByGraph proc~no_slip no_slip proc~wall wall proc~wall->proc~no_slip proc~populate_ghost_primitive populate_ghost_primitive proc~populate_ghost_primitive->proc~wall 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

no_slip

Source Code

      subroutine no_slip(face, bc, dims)
        !< No-slip wall boundary condition. All the 
        !< component of velocity throught face is zero
        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(x_speed, "anti", face, bc, dims)
        call copy3(y_speed, "anti", face, bc, dims)
        call copy3(z_speed, "anti", face, bc, dims)
        select case(turbulence)
          case("none")
            !do nothing
            continue
          case('sa', 'saBC')
            call copy3(tv  , "anti", face, bc, dims)
          case("sst", 'sst2003')
            call copy3(tk  , "anti", face, bc, dims)
            call set_omega_at_wall(face)
          case("kkl")
            call copy3(tk  , "anti", face, bc, dims)
            call copy3(tkl , "anti", 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 no_slip