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

sixDoFRigidBodyDisplacementPointPatchVectorField: Added optional under-relaxation

Browse Source
This commit is contained in:
Henry
2012-08-24 15:15:24 +01:00
parent 1ee6288d39
commit 20b5d43cd6
2 changed files with 19 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

21
src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyDisplacement/sixDoFRigidBodyDisplacementPointPatchVectorField.C
View File

@ -52,7 +52,8 @@ sixDoFRigidBodyDisplacementPointPatchVectorField
rhoInf_(1.0), rhoInf_(1.0),
rhoName_("rho"), rhoName_("rho"),
lookupGravity_(-1), lookupGravity_(-1),
g_(vector::zero) g_(vector::zero),
relaxationFactor_(1)
{} {}
@ -69,7 +70,8 @@ sixDoFRigidBodyDisplacementPointPatchVectorField
rhoInf_(1.0), rhoInf_(1.0),
rhoName_(dict.lookupOrDefault<word>("rhoName", "rho")), rhoName_(dict.lookupOrDefault<word>("rhoName", "rho")),
lookupGravity_(-1), lookupGravity_(-1),
g_(vector::zero) g_(vector::zero),
relaxationFactor_(dict.lookupOrDefault<scalar>("relaxationFactor", 1))
{ {
if (rhoName_ == "rhoInf") if (rhoName_ == "rhoInf")
{ {
@ -112,7 +114,8 @@ sixDoFRigidBodyDisplacementPointPatchVectorField
rhoInf_(ptf.rhoInf_), rhoInf_(ptf.rhoInf_),
rhoName_(ptf.rhoName_), rhoName_(ptf.rhoName_),
lookupGravity_(ptf.lookupGravity_), lookupGravity_(ptf.lookupGravity_),
g_(ptf.g_) g_(ptf.g_),
relaxationFactor_(ptf.relaxationFactor_)
{} {}
@ -129,7 +132,8 @@ sixDoFRigidBodyDisplacementPointPatchVectorField
rhoInf_(ptf.rhoInf_), rhoInf_(ptf.rhoInf_),
rhoName_(ptf.rhoName_), rhoName_(ptf.rhoName_),
lookupGravity_(ptf.lookupGravity_), lookupGravity_(ptf.lookupGravity_),
g_(ptf.g_) g_(ptf.g_),
relaxationFactor_(ptf.relaxationFactor_)
{} {}
@ -235,7 +239,11 @@ void sixDoFRigidBodyDisplacementPointPatchVectorField::updateCoeffs()
Field<vector>::operator= Field<vector>::operator=
( (
motion_.currentPosition(initialPoints_) - initialPoints_ relaxationFactor_
*(
motion_.currentPosition(initialPoints_)
- initialPoints_
)
); );
fixedValuePointPatchField<vector>::updateCoeffs(); fixedValuePointPatchField<vector>::updateCoeffs();
@ -258,6 +266,9 @@ void sixDoFRigidBodyDisplacementPointPatchVectorField::write(Ostream& os) const
os.writeKeyword("g") << g_ << token::END_STATEMENT << nl; os.writeKeyword("g") << g_ << token::END_STATEMENT << nl;
} }
os.writeKeyword("relaxationFactor")
<< relaxationFactor_ << token::END_STATEMENT << nl;
motion_.write(os); motion_.write(os);
initialPoints_.writeEntry("initialPoints", os); initialPoints_.writeEntry("initialPoints", os);

3
src/postProcessing/functionObjects/forces/pointPatchFields/derived/sixDoFRigidBodyDisplacement/sixDoFRigidBodyDisplacementPointPatchVectorField.H
View File

@ -84,6 +84,9 @@ class sixDoFRigidBodyDisplacementPointPatchVectorField
//- Gravity vector to store when not available from the db //- Gravity vector to store when not available from the db
vector g_; vector g_;
//- Optional under-relaxation factor for the motion
scalar relaxationFactor_;
public: public: