add_wall_process.h

Go to the documentation of this file.

//
//   Project Name:        Kratos
//   Last Modified by:    $Author: anonymous $
//   Date:                $Date: 2007-11-06 12:34:26 $
//   Revision:            $Revision: 1.4 $
//
//  this process save structural elements in a separate list
 
#if !defined(KRATOS_ADD_WALL_PROCESS_INCLUDED )
#define  KRATOS_ADD_WALL_PROCESS_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/element.h"
#include "includes/model_part.h"
//#include "custom_elements/updated_lagrangian_fluid.h"
//#include "custom_elements/updated_lagrangian_fluid3D.h"
//#include "custom_elements/updated_lagrangian_fluid_inc.h"
//#include "custom_elements/updated_lagrangian_fluid3D_inc.h"
 
 
namespace Kratos
{
 
 
 
 
 
 
 
 
class AddWallProcess
    : public Process
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION(AddWallProcess);
 
 
    AddWallProcess()
    //ModelPart& fluid_model_part, ModelPart& structure_model_part, ModelPart& combined_model_part)
    //: mr_fluid_model_part(fluid_model_part), mr_structure_model_part(structure_model_part), mr_combined_model_part(combined_model_part)
    {
    //KRATOS_WATCH(" INSIDE ADD WALL NODES CONSTRUCTOR") 
    }
 
    ~AddWallProcess() override
    {
    }
 
 
 
    //  void operator()()
    //  {
    //      MergeParts();
    //  }
 
 
 
    void AddWall(ModelPart& fluid_model_part, ModelPart& wall_model_part)
    {
        KRATOS_TRY  
 
    //COMPUTING AVERAGE NODAL_H in FLUID
    double av_mesh_size=0.0;
    for(ModelPart::NodesContainerType::iterator in = fluid_model_part.NodesBegin() ;
                in != fluid_model_part.NodesEnd() ; ++in)
        {
    av_mesh_size+=in->FastGetSolutionStepValue(NODAL_H);
    }
 
    unsigned int n_fluid_nodes=fluid_model_part.Nodes().size();
    av_mesh_size/=n_fluid_nodes;
 
    unsigned wall_node_id=n_fluid_nodes+1;
 
    if (av_mesh_size==0)
        KRATOS_THROW_ERROR(std::logic_error,"your wall nodes will have NODAL_H=0","");
        for(ModelPart::NodesContainerType::iterator in = wall_model_part.NodesBegin() ;
                in != wall_model_part.NodesEnd() ; ++in)
        {
        //set the NODAL_H in the wall node to the average size and then add this node to the fluid model part
        in->FastGetSolutionStepValue(NODAL_H)=av_mesh_size;
        in->SetId(wall_node_id);
        wall_node_id++;
        //ADDING ALL THE SECOND MOULD WALL NODES TO THE MODEL PART
        fluid_model_part.AddNode(*(in.base()),0);
        
        
        }
    
        //fluid_model_part.AddNodes(wall_model_part.NodesBegin(), wall_model_part.NodesEnd());  
        //sorting and renumbering the fluid elements
        unsigned int id=1;
        for(ModelPart::NodesContainerType::iterator in = fluid_model_part.NodesBegin() ;
                in != fluid_model_part.NodesEnd() ; ++in)
        {
            in->SetId(id);
            id++;
        }
 
 
    
        KRATOS_CATCH("")
    }
 
 
 
 
 
 
 
    std::string Info() const override
    {
        return "AddWallProcess";
    }
 
    void PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "AddWallProcess";
    }
 
    void PrintData(std::ostream& rOStream) const override
    {
    }
 
 
 
 
 
protected:
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
private:
 
 
    //ModelPart& mr_fluid_model_part;
    //ModelPart& mr_structure_model_part;
    //ModelPart& mr_combined_model_part;
 
 
 
 
 
 
 
 
 
 
//      AddWallProcess& operator=(AddWallProcess const& rOther);
 
//      AddWallProcess(AddWallProcess const& rOther);
 
 
 
}; // Class AddWallProcess
 
 
 
 
 
 
inline std::istream& operator >> (std::istream& rIStream,
                                  AddWallProcess& rThis);
 
inline std::ostream& operator << (std::ostream& rOStream,
                                  const AddWallProcess& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
 
}  // namespace Kratos.
 
#endif // KRATOS_ADD_WALL_PROCESS_INCLUDED  defined