save_lagrangian_surface_process_p.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main authors:    Author Julio Marti.
//
 
//  this process saves the boundary of a Lagrangian part for EMBEDDED technique
 
#if !defined(KRATOS_SAVE_LAGRANGIAN_SURFACE_PROCESS_P_INCLUDED )
#define  KRATOS_SAVE_LAGRANGIAN_SURFACE_PROCESS_P_INCLUDED
 
 
 
// System includes
#include <string>
#include <iostream>
#include <algorithm>
 
// External includes
 
 
// Project includes
#include "includes/define.h"
#include "processes/process.h"
#include "includes/node.h"
#include "includes/condition.h"
#include "includes/model_part.h"
#include "includes/deprecated_variables.h"
 
 
 
namespace Kratos
{
 
 
 
 
 
 
 
 
class SaveLagrangianSurfaceProcess_p
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(SaveLagrangianSurfaceProcess_p);
 
 
    SaveLagrangianSurfaceProcess_p()
    {
    }
 
    virtual ~SaveLagrangianSurfaceProcess_p()
    {
    }
 
 
 
//      void operator()()
//      {
//          SaveStructure();
//      }
 
 
 
    //FLUID_MODEL_PART is the whole Lagrangian fluid domain, surface model_part is its boundary
    void  SaveSurfaceConditions_p(ModelPart& lagrangian_model_part, ModelPart& surface_model_part)
    {
        KRATOS_TRY
        surface_model_part.Elements().clear();
        surface_model_part.Conditions().clear();
        surface_model_part.Nodes().clear();
 
    KRATOS_WATCH("SaveSurfaceConditions");
 
        //ModelPart::IndexType default_index = 0;
 
    for(ModelPart::NodesContainerType::iterator in = lagrangian_model_part.NodesBegin() ;
                in != lagrangian_model_part.NodesEnd() ; ++in)
        {
        if (in->FastGetSolutionStepValue(IS_BOUNDARY)==1.0 && in->GetValue(NEIGHBOUR_CONDITIONS).size()!=0)
            surface_model_part.Nodes().push_back(*(in.base()));
 
        }
    int id = 1;
    Properties::Pointer properties = lagrangian_model_part.GetMesh().pGetProperties(1);
        for(ModelPart::ConditionsContainerType::iterator ic = lagrangian_model_part.ConditionsBegin() ; ic != lagrangian_model_part.ConditionsEnd() ; ++ic)
        {
        //surface_model_part.Conditions().push_back(*(ic.base()));
        Geometry< Node >& geom = ic->GetGeometry();
 
            std::vector<std::size_t> NodeIds(3);
 
        for(int i=0; i<3; i++) NodeIds[i]= geom[i].Id();
 
        surface_model_part.CreateNewElement("Element3D3N",id, NodeIds,0);//,default_index);
 
        id = id + 1;
        }
 
 
 
 
 
        KRATOS_CATCH("")
    }
 
 
 
 
 
 
 
 
    virtual std::string Info() const override
    {
        return "SaveLagrangianSurfaceProcess_p";
    }
 
    virtual void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "SaveLagrangianSurfaceProcess_p";
    }
 
    virtual void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    //ModelPart& mr_fluid_model_part;
    //ModelPart& mr_structure_model_part;
 
 
 
 
 
 
 
 
 
 
//      SaveLagrangianSurfaceProcess& operator=(SaveLagrangianSurfaceProcess const& rOther);
 
//      SaveLagrangianSurfaceProcess(SaveLagrangianSurfaceProcess const& rOther);
 
 
 
}; // Class SaveLagrangianSurfaceProcess
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  SaveLagrangianSurfaceProcess_p& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const SaveLagrangianSurfaceProcess_p& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_SAVE_LAGRANGIAN_SURFACE_PROCESS_INCLUDED  defined