This function calculates the sigma for the planar viscosity spongelayer and multiply with targetState 'viscosity'. This function is currectly used to define viscosity sponge in musubi.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(tem_spongeLayer_plane_type) | :: | me |
Spacetime function to evaluate |
|||
| real(kind=rk), | intent(in) | :: | coord(n,3) |
barycentric Ids of an elements. 1st index goes over number of elements and 2nd index goes over x,y,z coordinates |
||
| integer, | intent(in) | :: | n |
Number of arrays to return |
return value
function viscSpongelayer_plane_for_coord(me, coord, n) & & result(res) ! -------------------------------------------------------------------------- !> Spacetime function to evaluate type(tem_spongeLayer_plane_type) :: me !> Number of arrays to return integer, intent(in) :: n !> barycentric Ids of an elements. !! 1st index goes over number of elements and !! 2nd index goes over x,y,z coordinates real(kind=rk), intent( in ) :: coord(n,3) !> return value real(kind=rk) :: res(n) ! -------------------------------------------------------------------------- integer :: i real(kind=rk) :: sigma, origin(3), normal(3), vec(3), proj_len ! -------------------------------------------------------------------------- origin(:) = me%origin normal(:) = me%normal do i = 1,n vec(:) = coord(i,:) - origin(:) proj_len = (vec(1)*normal(1)+ vec(2)*normal(2)+vec(3)*normal(3))/ & (normal(1)**2 + normal(2)**2 + normal(3)**2) sigma = 1.0_rk + (me%dampFactor-1.0_rk)*((proj_len)**me%dampExponent) if (proj_len > 0 .and. proj_len < 1) then res(i) = sigma * me%targetState(1) else if (proj_len > 1) then res(i) = me%dampFactor * me%targetState(1) else res(i) = me%targetState(1) end if end do end function viscSpongelayer_plane_for_coord