src/topoChangerFvMesh/linearValveFvMesh/linearValveFvMesh.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "linearValveFvMesh.H"
  27 #include "Time.H"
  28 #include "slidingInterface.H"
  29 #include "mapPolyMesh.H"
  30 #include "polyTopoChange.H"
  31 #include "addToRunTimeSelectionTable.H"
  32
  33 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37     defineTypeNameAndDebug(linearValveFvMesh, 0);
  38
  39     addToRunTimeSelectionTable(topoChangerFvMesh, linearValveFvMesh, IOobject);
  40 }
  41
  42
  43 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  44
  45 void Foam::linearValveFvMesh::addZonesAndModifiers()
  46 {
  47     // Add zones and modifiers for motion action
  48
  49     if
  50     (
  51         pointZones().size()
  52      || faceZones().size()
  53      || cellZones().size()
  54      || topoChanger_.size()
  55     )
  56     {
  57         Info<< "void linearValveFvMesh::addZonesAndModifiers() : "
  58             << "Zones and modifiers already present.  Skipping."
  59             << endl;
  60
  61         return;
  62     }
  63
  64     Info<< "Time = " << time().timeName() << endl
  65         << "Adding zones and modifiers to the mesh" << endl;
  66
  67     // Add zones
  68     List<pointZone*> pz(1);
  69
  70     // Add an empty zone for cut points
  71
  72     pz[0] = new pointZone
  73     (
  74         "cutPointZone",
  75         labelList(0),
  76         0,
  77         pointZones()
  78     );
  79
  80
  81     // Do face zones for slider
  82
  83     List<faceZone*> fz(3);
  84
  85     // Inner slider
  86     const word innerSliderName(motionDict_.subDict("slider").lookup("inside"));
  87     const polyPatch& innerSlider = boundaryMesh()[innerSliderName];
  88
  89     labelList isf(innerSlider.size());
  90
  91     forAll(isf, i)
  92     {
  93         isf[i] = innerSlider.start() + i;
  94     }
  95
  96     fz[0] = new faceZone
  97     (
  98         "insideSliderZone",
  99         isf,
 100         boolList(innerSlider.size(), false),
 101         0,
 102         faceZones()
 103     );
 104
 105     // Outer slider
 106     const word outerSliderName(motionDict_.subDict("slider").lookup("outside"));
 107     const polyPatch& outerSlider = boundaryMesh()[outerSliderName];
 108
 109     labelList osf(outerSlider.size());
 110
 111     forAll(osf, i)
 112     {
 113         osf[i] = outerSlider.start() + i;
 114     }
 115
 116     fz[1] = new faceZone
 117     (
 118         "outsideSliderZone",
 119         osf,
 120         boolList(outerSlider.size(), false),
 121         1,
 122         faceZones()
 123     );
 124
 125     // Add empty zone for cut faces
 126     fz[2] = new faceZone
 127     (
 128         "cutFaceZone",
 129         labelList(0),
 130         boolList(0, false),
 131         2,
 132         faceZones()
 133     );
 134
 135     List<cellZone*> cz(0);
 136
 137     Info<< "Adding point, face and cell zones" << endl;
 138     addZones(pz, fz, cz);
 139
 140     // Add a topology modifier
 141     Info<< "Adding topology modifiers" << endl;
 142     topoChanger_.setSize(1);
 143     topoChanger_.set
 144     (
 145         0,
 146         new slidingInterface
 147         (
 148             "mixerSlider",
 149             0,
 150             topoChanger_,
 151             outerSliderName + "Zone",
 152             innerSliderName + "Zone",
 153             "cutPointZone",
 154             "cutFaceZone",
 155             outerSliderName,
 156             innerSliderName,
 157             slidingInterface::INTEGRAL,
 158             true                          // Attach-detach action
 159         )
 160     );
 161     topoChanger_.writeOpt() = IOobject::AUTO_WRITE;
 162
 163     // Write mesh
 164     write();
 165 }
 166
 167
 168 void Foam::linearValveFvMesh::makeSlidersDead()
 169 {
 170     const polyTopoChanger& topoChanges = topoChanger_;
 171
 172     // Enable layering
 173     forAll(topoChanges, modI)
 174     {
 175         if (isA<slidingInterface>(topoChanges[modI]))
 176         {
 177             topoChanges[modI].disable();
 178         }
 179         else
 180         {
 181             FatalErrorIn("void Foam::linearValveFvMesh::makeSlidersDead()")
 182                 << "Don't know what to do with mesh modifier "
 183                 << modI << " of type " << topoChanges[modI].type()
 184                 << abort(FatalError);
 185         }
 186     }
 187 }
 188
 189
 190 void Foam::linearValveFvMesh::makeSlidersLive()
 191 {
 192     const polyTopoChanger& topoChanges = topoChanger_;
 193
 194     // Enable sliding interface
 195     forAll(topoChanges, modI)
 196     {
 197         if (isA<slidingInterface>(topoChanges[modI]))
 198         {
 199             topoChanges[modI].enable();
 200         }
 201         else
 202         {
 203             FatalErrorIn("void Foam::linearValveFvMesh::makeLayersLive()")
 204                 << "Don't know what to do with mesh modifier "
 205                 << modI << " of type " << topoChanges[modI].type()
 206                 << abort(FatalError);
 207         }
 208     }
 209 }
 210
 211
 212 bool Foam::linearValveFvMesh::attached() const
 213 {
 214     const polyTopoChanger& topoChanges = topoChanger_;
 215
 216     bool result = false;
 217
 218     forAll(topoChanges, modI)
 219     {
 220         if (isA<slidingInterface>(topoChanges[modI]))
 221         {
 222             result =
 223                 result
 224              || refCast<const slidingInterface>(topoChanges[modI]).attached();
 225         }
 226     }
 227
 228     // Check thal all sliders are in sync (debug only)
 229     forAll(topoChanges, modI)
 230     {
 231         if (isA<slidingInterface>(topoChanges[modI]))
 232         {
 233             if
 234             (
 235                 result
 236              != refCast<const slidingInterface>(topoChanges[modI]).attached()
 237             )
 238             {
 239                 FatalErrorIn("bool linearValveFvMesh::attached() const")
 240                     << "Slider " << modI
 241                     << " named " << topoChanges[modI].name()
 242                     << " out of sync: Should be" << result
 243                     << abort(FatalError);
 244             }
 245         }
 246     }
 247
 248     if (result)
 249     {
 250         Info<< "linearValveFvMesh: attached!" << endl;
 251     }
 252     else
 253     {
 254         Info<< "linearValveFvMesh: detached!" << endl;
 255     }
 256
 257     return result;
 258 }
 259
 260
 261 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 262
 263 // Construct from components
 264 Foam::linearValveFvMesh::linearValveFvMesh(const IOobject& io)
 265 :
 266     topoChangerFvMesh(io),
 267     motionDict_
 268     (
 269         IOdictionary
 270         (
 271             IOobject
 272             (
 273                 "dynamicMeshDict",
 274                 time().constant(),
 275                 *this,
 276                 IOobject::MUST_READ_IF_MODIFIED,
 277                 IOobject::NO_WRITE,
 278                 false
 279             )
 280         ).subDict(typeName + "Coeffs")
 281     ),
 282     msPtr_(motionSolver::New(*this))
 283 {
 284     addZonesAndModifiers();
 285 }
 286
 287
 288 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 289
 290 Foam::linearValveFvMesh::~linearValveFvMesh()
 291 {}
 292
 293
 294 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 295
 296 void Foam::linearValveFvMesh::update()
 297 {
 298     // Detaching the interface
 299     if (attached())
 300     {
 301         Info<< "Decoupling sliding interfaces" << endl;
 302         makeSlidersLive();
 303
 304         // Changing topology by hand
 305         resetMorph();
 306         setMorphTimeIndex(3*time().timeIndex());
 307         updateMesh();
 308
 309         msPtr_->updateMesh();
 310     }
 311     else
 312     {
 313         Info<< "Sliding interfaces decoupled" << endl;
 314     }
 315
 316     // Perform mesh motion
 317     makeSlidersDead();
 318
 319     // Changing topology by hand
 320     resetMorph();
 321     setMorphTimeIndex(3*time().timeIndex() + 1);
 322     updateMesh();
 323
 324     msPtr_->updateMesh();
 325
 326     if (topoChangeMap.valid())
 327     {
 328         if (topoChangeMap().hasMotionPoints())
 329         {
 330             Info<< "Topology change; executing pre-motion" << endl;
 331             movePoints(topoChangeMap().preMotionPoints());
 332         }
 333     }
 334
 335     // Solve for motion
 336     msPtr_->solve();
 337
 338     movePoints(msPtr_->curPoints());
 339
 340     // Attach the interface
 341     Info<< "Coupling sliding interfaces" << endl;
 342     makeSlidersLive();
 343     resetMorph();
 344     setMorphTimeIndex(3*time().timeIndex() + 2);
 345     updateMesh();
 346
 347     Info<< "Moving points post slider attach" << endl;
 348
 349     msPtr_->updateMesh();
 350
 351     Info<< "Sliding interfaces coupled: " << attached() << endl;
 352 }
 353
 354
 355 // ************************************************************************* //