basset_force_tools.h

Go to the documentation of this file.

//
//   Project Name:        Kratos
//   Last Modified by:    $Author: G.Casas (gcasas@cimmne.upc.edu) $
//   Date:                $Date: 2016-5-13 08:56:42 $
//   Revision:            $Revision: 1.5 $
//
//
//README::::look to the key word "VERSION" if you want to find all the points where you have to change something so that you can pass from a kdtree to a bin data search structure;
 
#if !defined(KRATOS_BASSET_FORCE_TOOLS)
#define KRATOS_BASSET_FORCE_TOOLS
 
// /* External includes */
#ifdef _OPENMP
#include <omp.h>
#endif
 
// System includes
#include <vector>
 
// Project includes
#include "includes/model_part.h"
#include "utilities/timer.h"
#include "utilities/openmp_utils.h"
#include "processes/find_nodal_neighbours_process.h"
 
//Database includes
#include "custom_utilities/search/discrete_particle_configure.h"
#include "includes/define.h"
#include "custom_elements/discrete_element.h"
#include "custom_elements/swimming_particle.h"
#include "custom_constitutive/history_force_laws/boussinesq_basset_history_force_law.h"
#include "includes/define.h"
#include "custom_utilities/AuxiliaryFunctions.h"
#include "custom_elements/spheric_particle.h"
#include "swimming_DEM_application.h"
#include "utilities/geometry_utilities.h"
 
namespace Kratos
{
class KRATOS_API(SWIMMING_DEM_APPLICATION) BassetForceTools
{
public:
 
typedef ModelPart::ElementsContainerType::iterator  ElementIterator;
typedef ModelPart::NodesContainerType::iterator     NodeIterator;
typedef ModelPart::NodesContainerType               NodesArrayType;
 
KRATOS_CLASS_POINTER_DEFINITION(BassetForceTools);
 
BassetForceTools(Parameters& r_parameters)
    : mFirstTimeAppending(true), mNumberOfQuadratureStepsInWindow(0), mTimeWindow(0.0)
{
    Parameters default_parameters( R"(
            {
                "do_use_mae": false,
                "m": 10,
                "window_time_interval": 0.1,
                "type": 2
            }
            )" );
 
    r_parameters.ValidateAndAssignDefaults(default_parameters);
    mBassetForceType = r_parameters["type"].GetInt();
 
}
 
virtual ~BassetForceTools(){}
 
 
void FillDaitcheVectors(const int N, const int order, const int n_steps_per_quad_step);
void FillHinsbergVectors(ModelPart& r_model_part, const int m, const int n_quad_delta_times_window);
void AppendIntegrands(ModelPart& r_model_part);
void AppendIntegrandsImplicit(ModelPart& r_model_part);
void AppendIntegrandsWindow(ModelPart& r_model_part);
 
//**************************************************************************************************************************************************
//**************************************************************************************************************************************************
 
private:
 
bool mFirstTimeAppending;
int mBassetForceType;
int mNumberOfQuadratureStepsInWindow;
int mNumberOfExponentials;
double mTimeWindow;
 
//**************************************************************************************************************************************************
//**************************************************************************************************************************************************
 
DenseVector<unsigned int> mElementsPartition;
 
 
DenseVector<unsigned int>& GetElementPartition()
{
    return mElementsPartition;
}
 
ElementIterator GetElementPartitionBegin(ModelPart& r_model_part, unsigned int k)
{
    return r_model_part.GetCommunicator().LocalMesh().Elements().ptr_begin() + mElementsPartition[k];
}
 
ElementIterator GetElementPartitionEnd(ModelPart& r_model_part, unsigned int k)
{
    return r_model_part.GetCommunicator().LocalMesh().Elements().ptr_begin() + mElementsPartition[k + 1];
}
 
//**************************************************************************************************************************************************
//**************************************************************************************************************************************************
 
}; // Class BassetForceTools
} // namespace Kratos.
 
#endif // KRATOS_CREATE_AND_DESTROY  defined