pressure_face_inverse_density Function

public function pressure_face_inverse_density(transport, owner, nb) result(beta)

Face-centered inverse density coefficient for variable-density projection.

Arguments

Type IntentOptional Attributes Name
type(transport_properties_t), intent(in) :: transport
integer, intent(in) :: owner
integer, intent(in) :: nb

Return Value real(kind=rk)


Source Code

   function pressure_face_inverse_density(transport, owner, nb) result(beta)
      type(transport_properties_t), intent(in) :: transport
      integer, intent(in) :: owner
      integer, intent(in) :: nb
      real(rk) :: beta

      if (.not. allocated(transport%rho)) then
         call fatal_error('numerics', 'transport rho must be allocated for variable-density projection')
      end if

      if (nb > 0) then
         beta = half * (one / max(transport%rho(owner), tiny_safe) + &
                        one / max(transport%rho(nb), tiny_safe))
      else
         beta = one / max(transport%rho(owner), tiny_safe)
      end if
   end function pressure_face_inverse_density