src/navalFoam/navalFoam.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright held by original author
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 Application
  26     interDyMFoam
  27
  28 Description
  29     Solver for 2 incompressible fluids capturing the interface with
  30     dynamic mesh motion.
  31
  32     Turbulence modelling is generic, i.e.  laminar, RAS or LES may be selected.
  33
  34 Author
  35     Hrvoje Jasak, Wikki Ltd.  All rights reserved.
  36
  37 \*---------------------------------------------------------------------------*/
  38
  39 #include "fvCFD.H"
  40 #include "MULES.H"
  41 #include "dynamicFvMesh.H"
  42 #include "interfaceProperties.H"
  43 #include "subCycle.H"
  44 #include "twoPhaseMixture.H"
  45 #include "turbulenceModel.H"
  46 #include "numericalBeach.H"
  47
  48 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  49
  50 int main(int argc, char *argv[])
  51 {
  52
  53 #   include "setRootCase.H"
  54 #   include "createTime.H"
  55 #   include "createDynamicFvMesh.H"
  56 #   include "readGravitationalAcceleration.H"
  57 #   include "initContinuityErrs.H"
  58 #   include "initTotalVolume.H"
  59 #   include "createFields.H"
  60 #   include "readControls.H"
  61 #   include "correctPhi.H"
  62
  63 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  64
  65     // Execute function objects
  66     runTime.functionObjects().execute();
  67
  68     Info<< "\nStarting time loop\n" << endl;
  69
  70     while (runTime.run())
  71     {
  72 #       include "readControls.H"
  73 #       include "checkTotalVolume.H"
  74 #       include "CourantNo.H"
  75 #       include "setDeltaT.H"
  76
  77         // Make the fluxes absolute
  78         fvc::makeAbsolute(phi, U);
  79
  80         runTime++;
  81
  82         Info<< "Time = " << runTime.timeName() << nl << endl;
  83
  84         if (runTime.value() > startMotionTime)
  85         {
  86             bool meshChanged = mesh.update();
  87
  88 #           include "volContinuity.H"
  89
  90             // Mesh motion update
  91             if (correctPhi && meshChanged)
  92             {
  93 #               include "correctPhi.H"
  94             }
  95
  96             // Make the fluxes relative to the mesh motion
  97             fvc::makeRelative(phi, U);
  98
  99             if (meshChanged)
 100             {
 101 #               include "CourantNo.H"
 102             }
 103         }
 104
 105         // Pressure-velocity corrector
 106         int oCorr = 0;
 107         do
 108         {
 109             gh = (g & mesh.C());
 110             ghf = (g & mesh.Cf());
 111
 112             twoPhaseProperties.correct();
 113
 114 #           include "UEqn.H"
 115
 116             // --- PISO loop
 117             for (int corr=0; corr < nCorr; corr++)
 118             {
 119 #               include "pdEqn.H"
 120             }
 121
 122 #           include "limitU.H"
 123
 124             if (mesh.moving())
 125             {
 126                 // Make the fluxes relative
 127                 phi -= fvc::meshPhi(rho, U);
 128                 rhoPhi = phi*fvc::interpolate(rho);
 129             }
 130
 131 #           include "movingMeshRhoUContinuityErrs.H"
 132
 133             // Calculate static pressure
 134             if (oCorr == nOuterCorr - 1)
 135             {
 136                 // Calculate static pressure
 137                 Info << "Calculating static pressure" << endl;
 138 #               include "pEqn.H"
 139             }
 140
 141 #           include "alphaEqn.H"
 142
 143             turbulence->correct();
 144         } while (++oCorr < nOuterCorr);
 145
 146         runTime.write();
 147
 148         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
 149             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
 150             << nl << endl;
 151     }
 152
 153     Info<< "End\n" << endl;
 154
 155     return(0);
 156 }
 157
 158
 159 // ************************************************************************* //