programmer's documentation
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
Functions/Subroutines
vitens.f90 File Reference

This function computes the equivalent viscosity at faces for a 3x3 symetric tensor. More...

Functions/Subroutines

subroutine vitens (w1, iwarnp, weighf, weighb, viscf, viscb)
 
subroutine symmetric_matrix_product (sout, s1, s2)
 

Detailed Description

This function computes the equivalent viscosity at faces for a 3x3 symetric tensor.

Function/Subroutine Documentation

subroutine symmetric_matrix_product ( double precision, dimension(6)  sout,
double precision, dimension(6)  s1,
double precision, dimension(6)  s2 
)
Parameters
[out]soutsout = s1 * s2
[in]s1symmetric matrix
[in]s2symmetric matrix
subroutine vitens ( double precision, dimension(6,ncelet), target  w1,
integer  iwarnp,
double precision, dimension(2,nfac)  weighf,
double precision, dimension(nfabor)  weighb,
double precision, dimension(nfac)  viscf,
double precision, dimension(nfabor)  viscb 
)
Parameters
[in]w1cell viscosity symmetric tensor
[in]iwarnpverbosity
[out]weighfinner face weight between cells i and j $ \frac{\vect{IF} \cdot \tens{K}_\celli} {\norm{\tens{K}_\celli \cdot \vect{S}}^2} $ and $ \frac{\vect{FJ} \cdot \tens{K}_\cellj} {\norm{\tens{K}_\cellj \cdot \vect{S}}^2} $
[out]weighbboundary face weight $ \frac{\vect{IF} \cdot \tens{K}_\celli} {\norm{\tens{K}_\celli \cdot \vect{S}}^2} $
[out]viscfinner face viscosity (times surface divided distance)
[out]viscbinner face viscosity (surface, must be consistent with flux BCs)