deformedGeom.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2016 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 Application
25  deformedGeom
26 
27 Description
28  Deforms a polyMesh using a displacement field U and a scaling factor
29  supplied as an argument.
30 
31 \*---------------------------------------------------------------------------*/
32 
33 #include "argList.H"
34 #include "fvMesh.H"
35 #include "pointFields.H"
36 #include "IStringStream.H"
37 #include "volPointInterpolation.H"
38 
39 using namespace Foam;
40 
41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
42 
43 int main(int argc, char *argv[])
44 {
45  argList::validArgs.append("scaling factor");
46 
47  #include "setRootCase.H"
48 
49  const scalar scaleFactor = args.argRead<scalar>(1);
50 
51  #include "createTime.H"
52  #include "createMesh.H"
53 
54  volPointInterpolation pInterp(mesh);
55 
56  // Get times list
57  instantList Times = runTime.times();
58 
59  pointField zeroPoints(mesh.points());
60 
61  // skip "constant" time
62  for (label timeI = 1; timeI < Times.size(); ++timeI)
63  {
64  runTime.setTime(Times[timeI], timeI);
65 
66  Info<< "Time = " << runTime.timeName() << endl;
67 
68  IOobject Uheader
69  (
70  "U",
71  runTime.timeName(),
72  mesh,
74  );
75 
76  // Check U exists
77  if (Uheader.headerOk())
78  {
79  Info<< " Reading U" << endl;
80  volVectorField U(Uheader, mesh);
81 
82  pointField newPoints
83  (
84  zeroPoints
85  + scaleFactor*pInterp.interpolate(U)().primitiveField()
86  );
87 
88  mesh.polyMesh::movePoints(newPoints);
89  mesh.write();
90  }
91  else
92  {
93  Info<< " No U" << endl;
94  }
95 
96  Info<< endl;
97  }
98 
99  Info<< "End\n" << endl;
100 
101  return 0;
102 }
103 
104 
105 // ************************************************************************* //
U
Definition: pEqn.H:83
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
void size(const label)
Override size to be inconsistent with allocated storage.
Definition: ListI.H:76
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:253
static SLList< string > validArgs
A list of valid (mandatory) arguments.
Definition: argList.H:154
virtual const pointField & points() const
Return raw points.
Definition: polyMesh.C:979
dynamicFvMesh & mesh
Interpolate from cell centres to points (vertices) using inverse distance weighting.
T argRead(const label index) const
Read a value from the argument at index.
Definition: argListI.H:177
messageStream Info
Foam::argList args(argc, argv)
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:91
virtual bool write() const
Write mesh using IO settings from time.
Definition: fvMesh.C:870
Namespace for OpenFOAM.