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

Adding vectors for global pi and tau as well as a triplet of orientation vectors to visualise rotational motion to written fields for post-processing.

Browse Source
This commit is contained in:
graham
2009-01-22 17:20:48 +00:00
parent 2723998c7b
commit 1b9bef5dd3
1 changed files with 46 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

46
src/lagrangian/molecularDynamics/molecule/molecule/moleculeIO.C
View File

@ -169,6 +169,37 @@ void Foam::molecule::writeFields(const moleculeCloud& mC)
IOField<label> special(mC.fieldIOobject("special", IOobject::NO_READ), np); IOField<label> special(mC.fieldIOobject("special", IOobject::NO_READ), np);
IOField<label> id(mC.fieldIOobject("id", IOobject::NO_READ), np); IOField<label> id(mC.fieldIOobject("id", IOobject::NO_READ), np);
// Post processing fields
IOField<vector> piGlobal
(
mC.fieldIOobject("piGlobal", IOobject::NO_READ),
np
);
IOField<vector> tauGlobal
(
mC.fieldIOobject("tauGlobal", IOobject::NO_READ),
np
);
IOField<vector> orientation1
(
mC.fieldIOobject("orientation1", IOobject::NO_READ),
np
);
IOField<vector> orientation2
(
mC.fieldIOobject("orientation2", IOobject::NO_READ),
np
);
IOField<vector> orientation3
(
mC.fieldIOobject("orientation3", IOobject::NO_READ),
np
);
label i = 0; label i = 0;
forAllConstIter(moleculeCloud, mC, iter) forAllConstIter(moleculeCloud, mC, iter)
{ {
@ -182,6 +213,14 @@ void Foam::molecule::writeFields(const moleculeCloud& mC)
specialPosition[i] = mol.specialPosition_; specialPosition[i] = mol.specialPosition_;
special[i] = mol.special_; special[i] = mol.special_;
id[i] = mol.id_; id[i] = mol.id_;
piGlobal[i] = mol.Q_ & mol.pi_;
tauGlobal[i] = mol.Q_ & mol.tau_;
orientation1[i] = mol.Q_ & vector(1,0,0);
orientation2[i] = mol.Q_ & vector(0,1,0);
orientation3[i] = mol.Q_ & vector(0,0,1);
i++; i++;
} }
@ -193,6 +232,13 @@ void Foam::molecule::writeFields(const moleculeCloud& mC)
specialPosition.write(); specialPosition.write();
special.write(); special.write();
id.write(); id.write();
piGlobal.write();
tauGlobal.write();
orientation1.write();
orientation2.write();
orientation3.write();
} }