49 const List<point>& lPoints,
50 const List<point>& rPoints,
59 const vector si(lPoints[i] - lPoints.rcValue(i));
60 const point ci(0.5*(lPoints[i] + lPoints.rcValue(i)));
65 const vector sj(rPoints[j] - rPoints.rcValue(j));
66 const point cj(0.5*(rPoints[j] + rPoints.rcValue(j)));
87 const List<List<vector>>& compactFineSf,
88 const List<List<point>>& compactFineCf,
89 const DynamicList<List<point>>& compactPoints,
90 const DynamicList<label>& compactPatchId
97 forAll(visibleFaceFaces, facei)
101 Pout<<
"facei:" << facei <<
"/" << visibleFaceFaces.size()
105 const labelList& visibleFaces = visibleFaceFaces[facei];
109 const vector& Ai = compactSf[facei];
110 const scalar magAi =
mag(Ai);
112 forAll(visibleFaces, visibleFacei)
114 const label sloti = visibleFaces[visibleFacei];
115 const List<point>& lPoints = compactPoints[facei];
116 const List<point>& rPoints = compactPoints[sloti];
118 const scalar Fij2LI = calculateFij(lPoints, rPoints, alpha_);
120 Fij[facei][visibleFacei] = Fij2LI/magAi;
133 const dictionary&
dict 137 alpha_(
dict.getOrDefault(
"alpha", 0.21))
Different types of constants.
dimensionedScalar log(const dimensionedScalar &ds)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
void resize(const label len)
Adjust allocated size of list.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Ostream & endl(Ostream &os)
Add newline and flush stream.
List< labelList > labelListList
List of labelList.
Macros for easy insertion into run-time selection tables.
#define forAll(list, i)
Loop across all elements in list.
vectorField pointField
pointField is a vectorField.
virtual void calculate()
Calculate the view factors.
List< scalarList > scalarListList
List of scalarList.
constexpr scalar pi(M_PI)
viewFactor2LI(const fvMesh &mesh, const dictionary &dict)
Constructor.
int debug
Static debugging option.
defineTypeNameAndDebug(combustionModel, 0)
vector point
Point is a vector.
static scalar calculateFij(const List< point > &lPoints, const List< point > &rPoints, const scalar alpha)
Calculate view factor using the double-area integral.
Field< vector > vectorField
Specialisation of Field<T> for vector.
List< label > labelList
A List of labels.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
prefixOSstream Pout
OSstream wrapped stdout (std::cout) with parallel prefix.
addToRunTimeSelectionTable(functionObject, pointHistory, dictionary)
static constexpr const zero Zero
Global zero (0)