interCondensatingEvaporatingFoam.C 4.24 KB
Newer Older
1 2 3 4
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
5
    \\  /    A nd           | Copyright (C) 2016-2019 OpenCFD Ltd.
6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
     \\/     M anipulation  |
-------------------------------------------------------------------------------
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/>.

Application
    interCondensatingEvaporatingFoam

Group
    grpMultiphaseSolvers

Description
31
    Solver for two incompressible, non-isothermal immiscible fluids with
32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58
    phase-change (evaporation-condensation) between a fluid and its vapour.
    Uses a VOF (volume of fluid) phase-fraction based interface capturing
    approach.

    The momentum, energy and other fluid properties are of the "mixture" and a
    single momentum equation is solved.

    Turbulence modelling is generic, i.e. laminar, RAS or LES may be selected.

\*---------------------------------------------------------------------------*/

#include "fvCFD.H"
#include "CMULES.H"
#include "subCycle.H"
#include "interfaceProperties.H"
#include "twoPhaseMixtureEThermo.H"
#include "temperaturePhaseChangeTwoPhaseMixture.H"
#include "turbulentTransportModel.H"
#include "turbulenceModel.H"
#include "pimpleControl.H"
#include "fvOptions.H"
#include "fixedFluxPressureFvPatchScalarField.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
59 60 61 62 63 64 65 66 67
    argList::addNote
    (
        "Solver for two incompressible, non-isothermal immiscible fluids with"
        " phase-change,"
        " using VOF phase-fraction based interface capturing.\n"
        "With optional mesh motion and mesh topology changes including"
        " adaptive re-meshing."
    );

68 69
    #include "postProcess.H"

70
    #include "addCheckCaseOptions.H"
71
    #include "setRootCaseLists.H"
72 73
    #include "createTime.H"
    #include "createMesh.H"
74
    #include "createControl.H"
75 76 77 78 79 80
    #include "createFields.H"
    #include "createFvOptions.H"
    #include "createTimeControls.H"
    #include "CourantNo.H"
    #include "setInitialDeltaT.H"

81 82
    volScalarField& T = thermo->T();

83 84
    turbulence->validate();

85 86 87 88 89 90
    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

    Info<< "\nStarting time loop\n" << endl;

    while (runTime.run())
    {
91
        #include "createTimeControls.H"
92 93 94
        #include "CourantNo.H"
        #include "setDeltaT.H"

95
        ++runTime;
96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112

        Info<< "Time = " << runTime.timeName() << nl << endl;

        // --- Pressure-velocity PIMPLE corrector loop
        while (pimple.loop())
        {
            #include "alphaControls.H"

            surfaceScalarField rhoPhi
            (
                IOobject
                (
                    "rhoPhi",
                    runTime.timeName(),
                    mesh
                ),
                mesh,
113
                dimensionedScalar(dimMass/dimTime, Zero)
114 115 116 117
            );

            mixture->correct();

118 119 120 121
            #include "alphaEqnSubCycle.H"

            solve(fvm::ddt(rho) + fvc::div(rhoPhi));

122 123 124 125 126 127 128 129 130 131 132 133 134 135 136
            #include "UEqn.H"
            #include "eEqn.H"

            // --- Pressure corrector loop
            while (pimple.correct())
            {
                #include "pEqn.H"
            }

            if (pimple.turbCorr())
            {
                turbulence->correct();
            }
        }

137 138
        rho = alpha1*rho1 + alpha2*rho2;

139 140
        runTime.write();

141
        runTime.printExecutionTime(Info);
142 143 144 145 146 147 148 149 150
    }

    Info<< "End\n" << endl;

    return 0;
}


// ************************************************************************* //