/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     |
    \\  /    A nd           | Copyright (C) 2009-2010 OpenCFD Ltd.
     \\/     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 2 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, write to the Free Software Foundation,
    Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA

Class
    Foam::kinematicSingleLayer

Description
    Kinematic form of single-cell layer surface film model

SourceFiles
    kinematicSingleLayer.C

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

#ifndef kinematicSingleLayer_H
#define kinematicSingleLayer_H

#include "surfaceFilmModel.H"
#include "fvMesh.H"
#include "volFields.H"
#include "surfaceFields.H"
#include "fvMatrices.H"

#include "injectionModelList.H"

#define DEBUGT(x) std::cout << tab.c_str()<<"["<< __FILE__ << ":" << __LINE__ << "] "<< #x " = " << x << std::endl;
#define TRACET(s) std::cout << tab.c_str()<<"["<< __FILE__ << ":" << __LINE__ << "] "<< #s << std::endl; s;
#define DEBUG(x) std::cout << "["<< __FILE__ << ":" << __LINE__ << "] "<< #x " = " << x << std::endl;
#define TRACE(s) std::cout << "["<< __FILE__ << ":" << __LINE__ << "] "<< #s << std::endl; s;
#define INFO Info << tab.c_str()


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

namespace Foam
{
namespace regionModels
{
namespace surfaceFilmModels
{

/*---------------------------------------------------------------------------*\
                   Class kinematicSingleLayer Declaration
\*---------------------------------------------------------------------------*/

class kinematicSingleLayer
:
    public surfaceFilmModel
{
private:

    // Private member functions

        //- Disallow default bitwise copy construct
        kinematicSingleLayer(const kinematicSingleLayer&);

        //- Disallow default bitwise assignment
        void operator=(const kinematicSingleLayer&);


protected:

    // Protected data

            string tab;
            void tabAdd();
            void tabSubtract();

        // Solution parameters

            //- Momentum predictor
            Switch momentumPredictor_;

            //- Number of outer correctors
            label nOuterCorr_;

            //- Number of PISO-like correctors
            label nCorr_;

            //- Number of non-orthogonal correctors
            label nNonOrthCorr_;

            //- Cumulative continuity error
            scalar cumulativeContErr_;


        // Model parameters

            //- Skin frition coefficient for film/primary region interface
            scalar Cf_;


        // Thermo properties

            // Fields

                //- Density / [kg/m3]
                volScalarField rho_;

                //- Dynamic viscosity / [Pa.s]
                volScalarField mu_;

                //- Surface tension / [m/s2]
                volScalarField sigma_;


        // Fields

            //- Film thickness / [m]
            volScalarField delta_;

            //- Velocity - mean / [m/s]
            volVectorField U_;

            //- Velocity - surface / [m/s]
            volVectorField Us_;

            //- Velocity - wall / [m/s]
            volVectorField Uw_;

            //- Film thickness*density (helper field) / [kg/m2]
            volScalarField deltaRho_;

            //- Mass flux (includes film thickness) / [kg.m/s]
            surfaceScalarField phi_;


            // Transfer fields

                //- Film mass available for transfer to the primary region
                volScalarField primaryMassTrans_;

                //- Film mass available for transfer to cloud
                volScalarField cloudMassTrans_;

                //- Parcel diameters originating from film to cloud
                volScalarField cloudDiameterTrans_;


        // Source term fields

            // Film region - registered to the film region mesh
            // Note: need boundary value mapped from primary region, and then
            // pushed into the patch internal field

                //- Momementum / [kg/m/s2]
                volVectorField USp_;

                //- Pressure / [Pa]
                volScalarField pSp_;

                //- Mass / [kg/m2/s]
                volScalarField rhoSp_;

                //- Impingment Mass / [kg/m2/s]
                volScalarField massImp_;


            // Primary region - registered to the primary region mesh
            // Internal use only - not read-in

                //- Momementum / [kg/m/s2]
                volVectorField USpPrimary_;

                //- Pressure / [Pa]
                volScalarField pSpPrimary_;

                //- Mass / [kg/m2/s]
                volScalarField rhoSpPrimary_;


        // Fields mapped from primary region - registered to the film region
        // Note: need both boundary AND patch internal fields to be mapped

            //- Velocity / [m/s]
            volVectorField UPrimary_;

            //- Pressure / [Pa]
            volScalarField pPrimary_;

            //- Density / [kg/m3]
            volScalarField rhoPrimary_;

            //- Viscosity / [Pa.s]
            volScalarField muPrimary_;


        // Sub-models

            //- Available mass for transfer via sub-models
            scalarField availableMass_;

            //- Cloud injection
            injectionModelList injection_;


       // Checks

           //- Cumulative mass added via sources [kg]
           scalar addedMassTotal_;


    // Protected member functions

        //- Read control parameters from dictionary
        virtual bool read();

        //- Correct the thermo fields
        virtual void correctThermoFields();

        //- Reset source term fields
        virtual void resetPrimaryRegionSourceTerms();

        //- Transfer thermo fields from the primary region to the film region
        virtual void transferPrimaryRegionThermoFields();

        //- Transfer source fields from the primary region to the film region
        virtual void transferPrimaryRegionSourceFields();

        // Explicit pressure source contribution
        virtual tmp<volScalarField> pu();

        // Implicit pressure source coefficient
        virtual tmp<volScalarField> pp();

        //- Update the film sub-models
        virtual void updateSubmodels();

        //- Continuity check
        virtual void continuityCheck();

        //- Update film surface velocities
        virtual void updateSurfaceVelocities();

        //- Return the stress term for the momentum equation
        virtual tmp<fvVectorMatrix> tau(volVectorField& dU) const;

        //- Constrain a film region master/slave boundaries of a field to a
        //  given value
        template<class Type>
        void constrainFilmField
        (
            Type& field,
            const typename Type::cmptType& value
        );


        // Equations

            //- Solve continuity equation
            virtual void solveContinuity();

            //- Solve for film velocity
            virtual tmp<fvVectorMatrix> solveMomentum
            (
                const volScalarField& pu,
                const volScalarField& pp
            );

            //- Solve coupled velocity-thickness equations
            virtual void solveThickness
            (
                const volScalarField& pu,
                const volScalarField& pp,
                const fvVectorMatrix& UEqn
            );


public:

    //- Runtime type information
    TypeName("kinematicSingleLayer");


    // Constructors

        //- Construct from components
        kinematicSingleLayer
        (
            const word& modelType,
            const fvMesh& mesh,
            const dimensionedVector& g,
            const bool readFields = true
        );


    //- Destructor
    virtual ~kinematicSingleLayer();


    // Member Functions

        // Solution parameters

            //- Courant number evaluation
            virtual scalar CourantNumber() const;

            //- Return the momentum predictor
            inline const Switch& momentumPredictor() const;

            //- Return the number of outer correctors
            inline label nOuterCorr() const;

            //- Return the number of PISO correctors
            inline label nCorr() const;

            //- Return the number of non-orthogonal correctors
            inline label nNonOrthCorr() const;


        // Model parameters

            //- Return the skin friction coefficient
            inline scalar Cf() const;


        // Thermo properties

            //- Return const access to the dynamic viscosity / [Pa.s]
            inline const volScalarField& mu() const;

            //- Return const access to the surface tension / [m/s2]
            inline const volScalarField& sigma() const;


        // Fields

            //- Return const access to the film thickness / [m]
            inline const volScalarField& delta() const;

            //- Return the film velocity [m/s]
            virtual const volVectorField& U() const;

            //- Return the film surface velocity [m/s]
            virtual const volVectorField& Us() const;

            //- Return the film wall velocity [m/s]
            virtual const volVectorField& Uw() const;

            //- Return the film flux [kg.m/s]
            virtual const surfaceScalarField& phi() const;

            //- Return the film density [kg/m3]
            virtual const volScalarField& rho() const;

            //- Return the film mean temperature [K]
            virtual const volScalarField& T() const;

            //- Return the film surface temperature [K]
            virtual const volScalarField& Ts() const;

            //- Return the film wall temperature [K]
            virtual const volScalarField& Tw() const;

            //- Return the wall heat flux
            virtual tmp<DimensionedField<scalar, volMesh> > qWall() const;

            //- Return the incident radiation
            //qrad virtual const volScalarField& qRad() const;
            virtual tmp<DimensionedField<scalar, volMesh> > qRad() const;

            //- Return the film specific heat capacity [J/kg/K]
            virtual const volScalarField& Cp() const;

            //- Return the film thermal conductivity [W/m/K]
            virtual const volScalarField& kappa() const;

            //- Return the film wettness (0=dry, 1=wet)
            virtual const volScalarField& omega() const;

            // Transfer fields - to the primary region

                //- Return mass transfer source - Eulerian phase only
                virtual tmp<volScalarField> primaryMassTrans() const;

                //- Return the film mass available for transfer to cloud
                virtual const volScalarField& cloudMassTrans() const;

                //- Return the parcel diameters originating from film to cloud
                virtual const volScalarField& cloudDiameterTrans() const;


        // External helper functions

            //- External hook to add sources to the film
            virtual void addSources
            (
                const label patchI,            // patchI on primary region
                const label faceI,             // faceI of patchI
                const scalar massSource,       // [kg]
                const vector& momentumSource,  // [kg.m/s] (tangential momentum)
                const scalar pressureSource,   // [kg.m/s] (normal momentum)
                const scalar energySource = 0  // [J]
            );


         // Source fields (read/write access)

            // Primary region

                //- Momementum / [kg/m/s2]
                inline volVectorField& USpPrimary();

                //- Pressure / [Pa]
                inline volScalarField& pSpPrimary();

                //- Mass / [kg/m2/s]
                inline volScalarField& rhoSpPrimary();


            // Film region

                //- Momentum / [kg/m/s2]
                inline volVectorField& USp();

                //- Pressure / [Pa]
                inline volScalarField& pSp();

                //- Mass / [kg/m2/s]
                inline volScalarField& rhoSp();

                //- Momentum / [kg/m/s2]
                inline const volVectorField& USp() const;

                //- Pressure / [Pa]
                inline const volScalarField& pSp() const;

                //- Mass / [kg/m2/s]
                inline const volScalarField& rhoSp() const;


        // Fields mapped from primary region

            //- Velocity / [m/s]
            inline const volVectorField& UPrimary() const;

            //- Pressure / [Pa]
            inline const volScalarField& pPrimary() const;

            //- Density / [kg/m3]
            inline const volScalarField& rhoPrimary() const;

            //- Viscosity / [Pa.s]
            inline const volScalarField& muPrimary() const;


        // Sub-models

            //- Injection
            inline injectionModelList& injection();


        // Helper functions

            //- Return the current film mass
            inline tmp<volScalarField> mass() const;

            //- Return the net film mass available over the next integration
            inline tmp<volScalarField> netMass() const;

            //- Return the gravity normal-to-patch component contribution
            inline tmp<volScalarField> gNorm() const;

            //- Return the gravity normal-to-patch component contribution
            //  Clipped so that only non-zero if g & nHat_ < 0
            inline tmp<volScalarField> gNormClipped() const;

            //- Return the gravity tangential component contributions
            inline tmp<volVectorField> gTan() const;


        // Evolution

            //- Pre-evolve film hook
            virtual void preEvolveRegion();

            //- Evolve the film equations
            virtual void evolveRegion();


        // Source fields

            // Mapped into primary region

                //- Return total mass source - Eulerian phase only
                virtual tmp<DimensionedField<scalar, volMesh> > Srho() const;

                //- Return mass source for specie i - Eulerian phase only
                virtual tmp<DimensionedField<scalar, volMesh> > Srho
                (
                    const label i
                ) const;

                //- Return enthalpy source - Eulerian phase only
                virtual tmp<DimensionedField<scalar, volMesh> > Sh() const;


        // I-O

            //- Provide some feedback
            virtual void info() const;
            void slottedCircle();
};


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

} // End namespace surfaceFilmModels
} // End namespace regionModels
} // End namespace Foam

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

#ifdef NoRepository
#   include "kinematicSingleLayerTemplates.C"
#endif

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

#include "kinematicSingleLayerI.H"

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

#endif

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