mass_conservation_check_process.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Simon Wenczowski
//
//
 
#ifndef KRATOS_MASS_CONSERVATION_CHECK_PROCESS_H
#define KRATOS_MASS_CONSERVATION_CHECK_PROCESS_H
 
// System includes
#include <string>
 
// External includes
 
// Project includes
#include "processes/process.h"
#include "custom_elements/two_fluid_navier_stokes.h"
#include "modified_shape_functions/tetrahedra_3d_4_modified_shape_functions.h"
#include "modified_shape_functions/triangle_2d_3_modified_shape_functions.h"
 
// Application includes
 
 
namespace Kratos
{
 
 
 
 
 
 
 
class KRATOS_API(FLUID_DYNAMICS_APPLICATION) MassConservationCheckProcess : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(MassConservationCheckProcess);
 
    typedef Node NodeType;
    typedef Geometry<NodeType> GeometryType;
 
 
    MassConservationCheckProcess(
        ModelPart& rModelPart,
        const bool PerformCorrections,
        const int CorrectionFreq,
        const bool WriteToLogFile,
        const std::string& LogFileName);
 
    MassConservationCheckProcess(
        ModelPart& rModelPart,
        Parameters& rParameters);
 
    ~MassConservationCheckProcess() override {}
 
 
 
 
    double ComputePositiveVolume();
 
    double ComputeNegativeVolume();
 
    double ComputeInterfaceArea();
 
    double ComputeFlowOverBoundary( const Kratos::Flags boundaryFlag );
 
    std::string Initialize();
 
    std::string ExecuteInTimeStep();
 
    // ///@}
    // ///@name Inquiry
    // ///@{
 
    // ///@}
    // ///@name Input and output
    // ///@{
 
    std::string Info() const override
    {
        std::stringstream buffer;
        buffer << "MassConservationCheckProcess";
        return buffer.str();
    }
 
    void PrintInfo(std::ostream& rOStream) const override {rOStream << "MassConservationCheckProcess";}
 
    void PrintData(std::ostream& rOStream) const override {}
 
 
 
 
private:
 
 
 
    // Reference to the model part
    const ModelPart& mrModelPart;
 
    // Process parameters
    int mCorrectionFreq = 1;
    bool mWriteToLogFile = true;
    bool mPerformCorrections = true;
    std::string mLogFileName = "mass_conservation.log";
 
    // Inital volume with negative distance field ("water" volume)
    double mInitialNegativeVolume = -1.0;
    // Inital volume with positive distance field ("air" volume)
    double mInitialPositiveVolume = -1.0;
 
    // Balance parameter resulting from an integration of the net inflow into the domain over time
    // The initial value is the "mInitialNegativeVolume" meaning "water" is considered here
    double mTheoreticalNegativeVolume = -1.0;
 
    // Net inflow into the domain (please consider that inflow at the outlet and outflow at the inlet are possible)
    double mQNet0 = 0.0;      // for the current time step (t)
    double mQNet1 = 0.0;      // for the past time step (t - 1)
    double mQNet2 = 0.0;      // for the past time step (t - 2)
 
 
 
    void ComputeVolumesAndInterface( double& positiveVolume, double& negativeVolume, double& interfaceArea );
 
    void CalculateNormal2D( array_1d<double,3>& An, const Geometry<Node >& pGeometry );
 
    void CalculateNormal3D( array_1d<double,3>& An, const Geometry<Node >& pGeometry );
 
    void ShiftDistanceField( double deltaDist );
 
    Triangle2D3<Node>::Pointer GenerateAuxTriangle( const Geometry<Node >& rGeom );
 
    void GenerateAuxLine(   const Geometry<Node >& rGeom,
                            const Vector& distance,
                            Line3D2<IndexedPoint>::Pointer& p_aux_line,
                            array_1d<double, 3>& aux_velocity1,
                            array_1d<double, 3>& aux_velocity2 );
 
 
 
 
 
 
    // DistanceModificationProcess() = delete;
 
    // DistanceModificationProcess& operator=(DistanceModificationProcess const& rOther) = delete;
 
    // DistanceModificationProcess(DistanceModificationProcess const& rOther) = delete;
 
 
}; // Class DistanceModificationProcess
 
 
 
 
 
};  // namespace Kratos.
 
#endif // KRATOS_MASS_CONSERVATION_CHECK_PROCESS_H