zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: Updated forces and forceCoeffs bin writing

Browse Source
This commit is contained in:
andy
2013-07-29 13:09:41 +01:00
parent ab782605c8
commit 1e711927a7
2 changed files with 36 additions and 49 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
src/postProcessing/functionObjects/forces/forceCoeffs/forceCoeffs.C
View File

@ -63,16 +63,14 @@ void Foam::forceCoeffs::writeFileHeader(const label i)
<< "# bins : " << nBin_ << nl
<< "# start : " << binMin_ << nl
<< "# delta : " << binDx_ << nl
<< "# direction : " << binDir_ << nl
<< "# Time";
<< "# direction : " << binDir_ << nl;
for (label j = 0; j < nBin_; j++)
{
const word jn = Foam::name(j);
file(i)
<< tab << "Cm" << jn << tab << "Cd" << jn << tab << "Cl" << jn;
}
file(i)
<< "# Time"
<< tab << "bin"
<< tab << "Cm"
<< tab << "Cd"
<< tab << "Cl";
}
else
{
@ -216,17 +214,16 @@ void Foam::forceCoeffs::write()
}
}
file(1)<< obr_.time().value();
forAll(coeffs[0], i)
{
file(1)
<< obr_.time().value()
<< tab << i
<< tab << coeffs[2][i]
<< tab << coeffs[1][i]
<< tab << coeffs[0][i];
<< tab << coeffs[0][i]
<< endl;
}
file(1) << endl;
}
if (log_)

60
src/postProcessing/functionObjects/forces/forces/forces.C
View File

@ -94,29 +94,19 @@ void Foam::forces::writeFileHeader(const label i)
<< "# bins : " << nBin_ << nl
<< "# start : " << binMin_ << nl
<< "# delta : " << binDx_ << nl
<< "# direction : " << binDir_ << nl
<< "# Time";
<< "# direction : " << binDir_ << nl;
for (label j = 0; j < nBin_; j++)
{
const word jn = Foam::name(j);
file(i)
<< "# Time"
<< tab << "bin"
<< tab << "forces(pressure,viscous,porous)"
<< tab << "moment(pressure,viscous,porous)";
file(i)
<< tab
<< "forces" << jn << "(pressure,viscous,porous) "
<< "moment" << jn << "(pressure,viscous,porous)";
}
if (localSystem_)
{
for (label j = 0; j < nBin_; j++)
{
const word jn = Foam::name(j);
file(i)
<< tab
<< "localForces" << jn << "(pressure,viscous,porous) "
<< "localMoments" << jn << "(pressure,viscous,porous)";
}
file(i)
<< tab << "localForces(pressure,viscous,porous)"
<< tab << "localMoment(pressure,viscous,porous)";
}
}
else
@ -435,27 +425,17 @@ void Foam::forces::writeBins()
f[0][i] += f[0][i-1];
f[1][i] += f[1][i-1];
f[2][i] += f[2][i-1];
m[0][i] += m[0][i-1];
m[1][i] += m[1][i-1];
m[2][i] += m[2][i-1];
}
}
file(1) << obr_.time().value();
forAll(f[0], i)
{
file(1)
<< tab
<< "(" << f[0][i] << "," << f[1][i] << "," << f[2][i] << ") "
<< "(" << m[0][i] << "," << m[1][i] << "," << m[2][i] << ")";
}
List<Field<vector> > lf(3);
List<Field<vector> > lm(3);
if (localSystem_)
{
List<Field<vector> > lf(3);
List<Field<vector> > lm(3);
lf[0] = coordSys_.localVector(force_[0]);
lf[1] = coordSys_.localVector(force_[1]);
lf[2] = coordSys_.localVector(force_[2]);
@ -475,17 +455,27 @@ void Foam::forces::writeBins()
lm[2][i] += lm[2][i-1];
}
}
}
forAll(lf[0], i)
forAll(f[0], i)
{
file(1)
<< obr_.time().value()
<< tab << i
<< tab << "(" << f[0][i] << "," << f[1][i] << "," << f[2][i] << ")"
<< tab << "(" << m[0][i] << "," << m[1][i] << "," << m[2][i] << ")";
if (localSystem_)
{
file(1)
<< tab
<< "(" << lf[0][i] << "," << lf[1][i] << "," << lf[2][i] << ") "
<< "(" << lf[0][i] << "," << lf[1][i] << "," << lf[2][i] << ")"
<< tab
<< "(" << lm[0][i] << "," << lm[1][i] << "," << lm[2][i] << ")";
}
}
file(1) << endl;
file(1) << endl;
}
}