fixed_mesh_ale_utilities.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                   Kratos default license: kratos/license.txt
//
//  Main author:     Ruben Zorrilla
//
//
 
#if !defined(KRATOS_FIXED_MESH_ALE_UTILITIES_H_INCLUDED )
#define  KRATOS_FIXED_MESH_ALE_UTILITIES_H_INCLUDED
 
// System includes
 
 
// External includes
 
 
// Project includes
#include "containers/model.h"
#include "includes/define.h"
#include "factories/linear_solver_factory.h"
#include "linear_solvers/linear_solver.h"
#include "processes/calculate_embedded_nodal_variable_from_skin_process.h"
#include "solving_strategies/builder_and_solvers/residualbased_block_builder_and_solver.h"
#include "solving_strategies/schemes/residualbased_incrementalupdate_static_scheme.h"
#include "solving_strategies/strategies/residualbased_linear_strategy.h"
#include "solving_strategies/strategies/implicit_solving_strategy.h"
#include "spaces/ublas_space.h"
 
// Application includes
 
 
namespace Kratos
{
 
 
 
 
 
 
 
 
 
 
class KRATOS_API(MESH_MOVING_APPLICATION) FixedMeshALEUtilities
{
public:
 
    typedef Element::NodesArrayType NodesArrayType;
    typedef UblasSpace<double, Matrix, Vector> LocalSpaceType;
    typedef UblasSpace<double, CompressedMatrix, Vector> SparseSpaceType;
    typedef LinearSolver<SparseSpaceType, LocalSpaceType> LinearSolverType;
    typedef LinearSolverFactory<SparseSpaceType, LocalSpaceType> LinearSolverFactoryType;
    typedef ResidualBasedIncrementalUpdateStaticScheme<SparseSpaceType, LocalSpaceType> SchemeType;
    typedef ResidualBasedLinearStrategy<SparseSpaceType, LocalSpaceType, LinearSolverType> StrategyType;
    typedef ResidualBasedBlockBuilderAndSolver<SparseSpaceType, LocalSpaceType, LinearSolverType> BuilderAndSolverType;
    typedef CalculateEmbeddedNodalVariableFromSkinProcess<array_1d<double, 3>, SparseSpaceType, LocalSpaceType, LinearSolverType> EmbeddedNodalVariableProcessArrayType;
 
    typedef typename SchemeType::Pointer SchemePointerType;
    typedef typename StrategyType::Pointer StrategyPointerType;
    typedef typename BuilderAndSolverType::Pointer BuilderAndSolverPointerType;
 
    KRATOS_CLASS_POINTER_DEFINITION(FixedMeshALEUtilities);
 
 
    FixedMeshALEUtilities(
        Model &rModel,
        Parameters &rParameters);
 
    virtual ~FixedMeshALEUtilities() = default;
 
 
 
 
    virtual void Initialize(ModelPart &rOriginModelPart);
 
    virtual void SetVirtualMeshValuesFromOriginMesh();
 
    virtual void ComputeMeshMovement(const double DeltaTime);
 
    virtual void UndoMeshMovement();
 
    template <unsigned int TDim>
    void ProjectVirtualValues(
        ModelPart &rOriginModelPart,
        const unsigned int BufferSize = 3);
 
 
 
 
 
 
    std::string Info() const;
 
    void PrintInfo(std::ostream& rOStream) const;
 
    void PrintData(std::ostream& rOStream) const;
 
 
 
private:
 
 
 
    ModelPart &mrVirtualModelPart;
    ModelPart &mrStructureModelPart;
    ModelPart *mpOriginModelPart = nullptr;
 
    double mSearchTolerance;
    unsigned int mSearchMaxResults;
 
    Parameters mEmbeddedNodalVariableSettings;
 
    std::vector<const Variable<double>*> mScalarVariablesList;
    std::vector<const Variable<array_1d<double,3>>*> mArrayVariablesList;
 
    LinearSolverType::Pointer mpLinearSolver = nullptr;
    StrategyPointerType mpMeshMovingStrategy = nullptr;
 
 
 
 
    Parameters GetDefaultParameters();
 
    void SetLinearSolverPointer(const Parameters &rLinearSolverSettings);
 
    virtual void FillVirtualModelPart(ModelPart &rOriginModelPart);
 
    virtual void CreateVirtualModelPartElements(const ModelPart &rOriginModelPart);
 
    void SetMeshMovingStrategy();
 
    void InitializeVirtualMeshValues();
 
    void InitializeMeshDisplacementFixity();
 
    void SetMeshDisplacementFixityFromOriginModelPart();
 
    void SetEmbeddedNodalMeshDisplacement();
 
    void SolveMeshMovementStrategy(const double DeltaTime);
 
    void RevertMeshDisplacementFixity();
 
 
 
 
 
 
    FixedMeshALEUtilities& operator=(FixedMeshALEUtilities const& rOther) = delete;
 
    FixedMeshALEUtilities(FixedMeshALEUtilities const& rOther) = delete;
 
}; // Class FixedMeshALEUtilities
 
 
inline std::ostream& operator << (
    std::ostream& rOStream,
    const FixedMeshALEUtilities& rThis);
 
}  // namespace Kratos.
 
#endif // KRATOS_FIXED_MESH_ALE_UTILITIES_H_INCLUDED  defined