GIT: interfaceTracking: relocate library

src/dynamicFaMesh/interfaceTrackingFvMesh/functionObjects/pointHistory/pointHistory.C

@@ -0,0 +1,215 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | www.openfoam.com
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+    Copyright (C) 2019 Zeljko Tukovic, FSB Zagreb.
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*----------------------------------------------------------------------------*/
+
+#include "pointHistory.H"
+#include "addToRunTimeSelectionTable.H"
+#include "IOmanip.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    defineTypeNameAndDebug(pointHistory, 0);
+
+    addToRunTimeSelectionTable
+    (
+        functionObject,
+        pointHistory,
+        dictionary
+    );
+}
+
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+bool Foam::pointHistory::writeData()
+{
+    const fvMesh& mesh =
+        time_.lookupObject<fvMesh>(polyMesh::defaultRegion);
+
+    vector position(Zero);
+
+    if (processor_ == Pstream::myProcNo())
+    {
+        position = mesh.points()[historyPointID_];
+    }
+
+    reduce(position, sumOp<vector>());
+
+    if (Pstream::master())
+    {
+        historyFilePtr_() << setprecision(12);
+
+        historyFilePtr_()
+            << time_.time().value() << tab
+            << position.x() << tab
+            << position.y() << tab
+            << position.z() << endl;
+    }
+
+    return true;
+}
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::pointHistory::pointHistory
+(
+    const word& name,
+    const Time& runTime,
+    const dictionary& dict
+)
+:
+    functionObject(name),
+    name_(name),
+    time_(runTime),
+    regionName_(polyMesh::defaultRegion),
+    historyPointID_(-1),
+    refHistoryPoint_(dict.lookup("refHistoryPoint")),
+    processor_(-1),
+    fileName_(dict.get<word>("fileName")),
+    historyFilePtr_()
+{
+    Info<< "Creating " << this->name() << " function object." << endl;
+
+    dict.readIfPresent("region", regionName_);
+    dict.readIfPresent("historyPointID", historyPointID_);
+
+    const fvMesh& mesh =
+        time_.lookupObject<fvMesh>(regionName_);
+
+    const vectorField& points = mesh.points();
+
+    List<scalar> minDist(Pstream::nProcs(), GREAT);
+
+    if (historyPointID_ == -1)
+    {
+        forAll(points, pointI)
+        {
+            scalar dist = mag(refHistoryPoint_ - points[pointI]);
+
+            if (dist < minDist[Pstream::myProcNo()])
+            {
+                minDist[Pstream::myProcNo()] = dist;
+                historyPointID_ = pointI;
+            }
+        }
+    }
+
+    Pstream::gatherList(minDist);
+    Pstream::scatterList(minDist);
+
+    processor_ = -1;
+    scalar min = GREAT;
+
+    forAll(minDist, procI)
+    {
+        if (minDist[procI] < min)
+        {
+            min = minDist[procI];
+            processor_ = procI;
+        }
+    }
+
+    if (processor_ == Pstream::myProcNo())
+    {
+        Pout<< "History point ID: " << historyPointID_ << nl
+            << "History point coordinates: "
+            << points[historyPointID_] << nl
+            << "Reference point coordinates: " << refHistoryPoint_
+            << endl;
+    }
+
+    // Create history file if not already created
+    if (historyFilePtr_.empty())
+    {
+        // File update
+        if (Pstream::master())
+        {
+            fileName historyDir;
+
+            word startTimeName =
+                time_.timeName(mesh.time().startTime().value());
+
+            if (Pstream::parRun())
+            {
+                // Put in undecomposed case (Note: gives problems for
+                // distributed data running)
+                historyDir = time_.path()/".."/"history"/startTimeName;
+            }
+            else
+            {
+                historyDir = time_.path()/"history"/startTimeName;
+            }
+
+            // Create directory if does not exist.
+            mkDir(historyDir);
+
+
+            // Open new file at start up
+
+            // OStringStream FileName;
+            // FileName() << "point_" << historyPointID_ << ".dat";
+
+            historyFilePtr_.reset
+            (
+                new OFstream(historyDir/fileName_)
+            );
+
+            // Add headers to output data
+            if (historyFilePtr_.valid())
+            {
+                historyFilePtr_()
+                    << "# Time" << tab << "X" << tab
+                        << "Y" << tab << "Z";
+
+                historyFilePtr_() << endl;
+            }
+        }
+    }
+
+    // Write start time data
+    writeData();
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::pointHistory::execute()
+{
+    return writeData();
+}
+
+
+bool Foam::pointHistory::read(const dictionary& dict)
+{
+    dict.readIfPresent("region", regionName_);
+
+    return true;
+}
+
+
+// ************************************************************************* //

src/dynamicFaMesh/interfaceTrackingFvMesh/functionObjects/pointHistory/pointHistory.H

@@ -0,0 +1,135 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | www.openfoam.com
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+    Copyright (C) 2019 Zeljko Tukovic, FSB Zagreb.
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::pointHistory
+
+Description
+
+SourceFiles
+    pointHistory.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef pointHistory_H
+#define pointHistory_H
+
+#include "functionObject.H"
+#include "dictionary.H"
+#include "fvMesh.H"
+#include "OFstream.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                        Class pointHistory Declaration
+\*---------------------------------------------------------------------------*/
+
+class pointHistory
+:
+    public functionObject
+{
+    // Private Data
+
+        //- Name
+        const word name_;
+
+        //- Reference to main object registry
+        const Time& time_;
+
+        //- Region name
+        word regionName_;
+
+        //- History point ID
+        label historyPointID_;
+
+        //- History point
+        vector refHistoryPoint_;
+
+        //- Processor of history point
+        label processor_;
+
+        //- Output file name
+        word fileName_;
+
+        //- Output file stream
+        autoPtr<OFstream> historyFilePtr_;
+
+    // Private Member Functions
+
+        //- Write data
+        bool writeData();
+
+        //- No copy construct
+        pointHistory(const pointHistory&) = delete;
+
+        //- No copy assignment
+        void operator=(const pointHistory&) = delete;
+
+public:
+
+    //- Runtime type information
+    TypeName("pointHistory");
+
+
+    // Constructors
+
+        //- Construct from components
+        pointHistory
+        (
+            const word& name,
+            const Time& runTime,
+            const dictionary&
+        );
+
+
+    // Member Functions
+
+        //- execute is called at each ++ or += of the time-loop
+        virtual bool execute();
+
+        //- Read and set the function object if its data has changed
+        virtual bool read(const dictionary& dict);
+
+        //- No-op
+        virtual bool write()
+        {
+            return true;
+        }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //

src/dynamicFaMesh/interfaceTrackingFvMesh/functionObjects/writeFreeSurface/writeFreeSurface.C

@@ -0,0 +1,114 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | www.openfoam.com
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+    Copyright (C) 2019 Zeljko Tukovic, FSB Zagreb.
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+\*----------------------------------------------------------------------------*/
+
+#include "writeFreeSurface.H"
+#include "addToRunTimeSelectionTable.H"
+#include "IOmanip.H"
+#include "interfaceTrackingFvMesh.H"
+
+// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
+
+namespace Foam
+{
+    defineTypeNameAndDebug(writeFreeSurface, 0);
+
+    addToRunTimeSelectionTable
+    (
+        functionObject,
+        writeFreeSurface,
+        dictionary
+    );
+}
+
+// * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
+
+bool Foam::writeFreeSurface::writeData()
+{
+    if (time_.outputTime())
+    {
+        const fvMesh& mesh =
+            time_.lookupObject<fvMesh>(polyMesh::defaultRegion);
+
+        interfaceTrackingFvMesh& itm =
+            refCast<interfaceTrackingFvMesh>
+            (
+                const_cast<dynamicFvMesh&>
+                (
+                    mesh.lookupObject<dynamicFvMesh>("fvSolution")
+                )
+            );
+
+        itm.writeVTKControlPoints();
+    }
+
+    return true;
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
+
+Foam::writeFreeSurface::writeFreeSurface
+(
+    const word& name,
+    const Time& runTime,
+    const dictionary& dict
+)
+:
+    functionObject(name),
+    name_(name),
+    time_(runTime),
+    regionName_(polyMesh::defaultRegion)
+{
+    Info<< "Creating " << this->name() << " function object." << endl;
+
+    dict.readIfPresent("region", regionName_);
+}
+
+
+// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
+
+bool Foam::writeFreeSurface::start()
+{
+    return writeData();
+}
+
+
+bool Foam::writeFreeSurface::execute()
+{
+    return writeData();
+}
+
+
+bool Foam::writeFreeSurface::read(const dictionary& dict)
+{
+    dict.readIfPresent("region", regionName_);
+
+    return true;
+}
+
+
+// ************************************************************************* //

src/dynamicFaMesh/interfaceTrackingFvMesh/functionObjects/writeFreeSurface/writeFreeSurface.H

@@ -0,0 +1,125 @@
+/*---------------------------------------------------------------------------*\
+  =========                 |
+  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
+   \\    /   O peration     |
+    \\  /    A nd           | www.openfoam.com
+     \\/     M anipulation  |
+-------------------------------------------------------------------------------
+    Copyright (C) 2019 Zeljko Tukovic, FSB Zagreb.
+-------------------------------------------------------------------------------
+License
+    This file is part of OpenFOAM.
+
+    OpenFOAM is free software: you can redistribute it and/or modify it
+    under the terms of the GNU General Public License as published by
+    the Free Software Foundation, either version 3 of the License, or
+    (at your option) any later version.
+
+    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
+    for more details.
+
+    You should have received a copy of the GNU General Public License
+    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
+
+Class
+    Foam::writeFreeSurface
+
+Description
+    A function object to write the control points on the free surface
+
+SourceFiles
+    writeFreeSurface.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef writeFreeSurface_H
+#define writeFreeSurface_H
+
+#include "functionObject.H"
+#include "dictionary.H"
+#include "fvMesh.H"
+#include "OFstream.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+/*---------------------------------------------------------------------------*\
+                        Class writeFreeSurface Declaration
+\*---------------------------------------------------------------------------*/
+
+class writeFreeSurface
+:
+    public functionObject
+{
+    // Private Data
+
+        //- Name
+        const word name_;
+
+        //- Reference to main object registry
+        const Time& time_;
+
+        //- Region name
+        word regionName_;
+
+
+    // Private Member Functions
+
+        //- Write data
+        bool writeData();
+
+        //- No copy construct
+        writeFreeSurface(const writeFreeSurface&) = delete;
+
+        //- No copy assignment
+        void operator=(const writeFreeSurface&) = delete;
+
+public:
+
+    //- Runtime type information
+    TypeName("writeFreeSurface");
+
+
+    // Constructors
+
+        //- Construct from components
+        writeFreeSurface
+        (
+            const word& name,
+            const Time& runTime,
+            const dictionary& dict
+        );
+
+
+    // Member Functions
+
+        //- start is called at the start of the time-loop
+        virtual bool start();
+
+        //- execute is called at each ++ or += of the time-loop
+        virtual bool execute();
+
+        //- Read and set the function object if its data has changed
+        virtual bool read(const dictionary& dict);
+
+        //- No-op
+        virtual bool write()
+        {
+            return true;
+        }
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //