41 Foam::forceCoeffs::forceCoeffs
46 const bool loadFromFiles
49 forces(name, obr, dict, loadFromFiles),
50 liftDir_(vector::zero),
51 dragDir_(vector::zero),
52 pitchAxis_(vector::zero),
76 dict.
lookup(
"liftDir") >> liftDir_;
77 dict.
lookup(
"dragDir") >> dragDir_;
78 dict.
lookup(
"pitchAxis") >> pitchAxis_;
81 dict.
lookup(
"magUInf") >> magUInf_;
84 dict.
lookup(
"lRef") >> lRef_;
85 dict.
lookup(
"Aref") >> Aref_;
92 if (forcesFilePtr_.valid())
95 <<
"# Time" <<
tab <<
"Cd" <<
tab <<
"Cl" <<
tab <<
"Cm" <<
endl;
121 scalar
pDyn = 0.5*rhoRef_*magUInf_*magUInf_;
126 scalar liftForce = totForce & liftDir_;
127 scalar dragForce = totForce & dragDir_;
128 scalar pitchMoment = totMoment & pitchAxis_;
130 scalar Cl = liftForce/(Aref_*
pDyn);
131 scalar Cd = dragForce/(Aref_*
pDyn);
132 scalar Cm = pitchMoment/(Aref_*lRef_*
pDyn);
137 << obr_.time().value() <<
tab
142 Info<<
"forceCoeffs output:" <<
nl
143 <<
" Cd = " << Cd <<
nl
144 <<
" Cl = " << Cl <<
nl
145 <<
" Cm = " << Cm <<
nl