applications/solvers/solidMechanics/deprecatedSolvers/materialModels/fvPatchFields/tractionDisplacement/tractionDisplacementFvPatchVectorField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation, either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "tractionDisplacementFvPatchVectorField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "volFields.H"
  29 #include "rheologyModel.H"
  30
  31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  32
  33 namespace Foam
  34 {
  35
  36 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  37
  38 tractionDisplacementFvPatchVectorField::
  39 tractionDisplacementFvPatchVectorField
  40 (
  41     const fvPatch& p,
  42     const DimensionedField<vector, volMesh>& iF
  43 )
  44 :
  45     fixedGradientFvPatchVectorField(p, iF),
  46     UName_("undefined"),
  47     rheologyName_("undefined"),
  48     traction_(p.size(), vector::zero),
  49     pressure_(p.size(), 0.0)
  50 {
  51     fvPatchVectorField::operator=(patchInternalField());
  52     gradient() = vector::zero;
  53 }
  54
  55
  56 tractionDisplacementFvPatchVectorField::
  57 tractionDisplacementFvPatchVectorField
  58 (
  59     const fvPatch& p,
  60     const DimensionedField<vector, volMesh>& iF,
  61     const dictionary& dict
  62 )
  63 :
  64     fixedGradientFvPatchVectorField(p, iF),
  65     UName_(dict.lookup("U")),
  66     rheologyName_(dict.lookup("rheology")),
  67     traction_("traction", dict, p.size()),
  68     pressure_("pressure", dict, p.size())
  69 {
  70     fvPatchVectorField::operator=(patchInternalField());
  71     gradient() = vector::zero;
  72     Info << "rf: " << rheologyName_ << endl;
  73 }
  74
  75
  76 tractionDisplacementFvPatchVectorField::
  77 tractionDisplacementFvPatchVectorField
  78 (
  79     const tractionDisplacementFvPatchVectorField& tdpvf,
  80     const fvPatch& p,
  81     const DimensionedField<vector, volMesh>& iF,
  82     const fvPatchFieldMapper& mapper
  83 )
  84 :
  85     fixedGradientFvPatchVectorField(tdpvf, p, iF, mapper),
  86     UName_(tdpvf.UName_),
  87     rheologyName_(tdpvf.rheologyName_),
  88     traction_(tdpvf.traction_, mapper),
  89     pressure_(tdpvf.pressure_, mapper)
  90 {}
  91
  92
  93 tractionDisplacementFvPatchVectorField::
  94 tractionDisplacementFvPatchVectorField
  95 (
  96     const tractionDisplacementFvPatchVectorField& tdpvf
  97 )
  98 :
  99     fixedGradientFvPatchVectorField(tdpvf),
 100     UName_(tdpvf.UName_),
 101     rheologyName_(tdpvf.rheologyName_),
 102     traction_(tdpvf.traction_),
 103     pressure_(tdpvf.pressure_)
 104 {}
 105
 106
 107 tractionDisplacementFvPatchVectorField::
 108 tractionDisplacementFvPatchVectorField
 109 (
 110     const tractionDisplacementFvPatchVectorField& tdpvf,
 111     const DimensionedField<vector, volMesh>& iF
 112 )
 113 :
 114     fixedGradientFvPatchVectorField(tdpvf, iF),
 115     UName_(tdpvf.UName_),
 116     rheologyName_(tdpvf.rheologyName_),
 117     traction_(tdpvf.traction_),
 118     pressure_(tdpvf.pressure_)
 119 {}
 120
 121
 122 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 123
 124 void tractionDisplacementFvPatchVectorField::autoMap
 125 (
 126     const fvPatchFieldMapper& m
 127 )
 128 {
 129     fixedGradientFvPatchVectorField::autoMap(m);
 130     traction_.autoMap(m);
 131     pressure_.autoMap(m);
 132 }
 133
 134
 135 // Reverse-map the given fvPatchField onto this fvPatchField
 136 void tractionDisplacementFvPatchVectorField::rmap
 137 (
 138     const fvPatchVectorField& ptf,
 139     const labelList& addr
 140 )
 141 {
 142     fixedGradientFvPatchVectorField::rmap(ptf, addr);
 143
 144     const tractionDisplacementFvPatchVectorField& dmptf =
 145         refCast<const tractionDisplacementFvPatchVectorField>(ptf);
 146
 147     traction_.rmap(dmptf.traction_, addr);
 148     pressure_.rmap(dmptf.pressure_, addr);
 149 }
 150
 151
 152 // Update the coefficients associated with the patch field
 153 void tractionDisplacementFvPatchVectorField::updateCoeffs()
 154 {
 155     if (updated())
 156     {
 157         return;
 158     }
 159
 160     // Looking up rheology
 161     const rheologyModel& rheology =
 162         this->db().objectRegistry::lookupObject<rheologyModel>(rheologyName_);
 163
 164     const scalarField mu = rheology.mu()().boundaryField()[patch().index()];
 165     const scalarField lambda =
 166         rheology.lambda()().boundaryField()[patch().index()];
 167
 168     vectorField n = patch().nf();
 169
 170     const fvPatchField<tensor>& gradU =
 171         patch().lookupPatchField<volTensorField, tensor>("grad(" +UName_ + ")");
 172
 173     gradient() =
 174     (
 175         (traction_ - (pressure_)*n)
 176       - (n & (mu*gradU.T() - (mu + lambda)*gradU))
 177       - n*lambda*tr(gradU)
 178     )/(2.0*mu + lambda);
 179
 180     fixedGradientFvPatchVectorField::updateCoeffs();
 181 }
 182
 183
 184 // Write
 185 void tractionDisplacementFvPatchVectorField::write(Ostream& os) const
 186 {
 187     fvPatchVectorField::write(os);
 188     os.writeKeyword("U") << UName_ << token::END_STATEMENT << nl;
 189     os.writeKeyword("rheology") << rheologyName_ << token::END_STATEMENT << nl;
 190     traction_.writeEntry("traction", os);
 191     pressure_.writeEntry("pressure", os);
 192     writeEntry("value", os);
 193 }
 194
 195
 196 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 197
 198 makePatchTypeField(fvPatchVectorField, tractionDisplacementFvPatchVectorField);
 199
 200 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 201
 202 } // End namespace Foam
 203
 204 // ************************************************************************* //