particle_generator_network.h

Go to the documentation of this file.

/* -------------------------------------------------------------------------*
 *                              SPHinXsys                                   *
 * --------------------------------------------------------------------------*
 * SPHinXsys (pronunciation: s'finksis) is an acronym from Smoothed Particle    *
 * Hydrodynamics for industrial compleX systems. It provides C++ APIs for   *
 * physical accurate simulation and aims to model coupled industrial dynamic *
 * systems including fluid, solid, multi-body dynamics and beyond with SPH  *
 * (smoothed particle hydrodynamics), a meshless computational method using *
 * particle discretization.                                                 *
 *                                                                          *
 * SPHinXsys is partially funded by German Research Foundation              *
 * (Deutsche Forschungsgemeinschaft) DFG HU1527/6-1, HU1527/10-1                *
 * and HU1527/12-1.                                                         *
 *                                                                           *
 * Portions copyright (c) 2017-2020 Technical University of Munich and      *
 * the authors' affiliations.                                               *
 *                                                                           *
 * Licensed under the Apache License, Version 2.0 (the "License"); you may   *
 * not use this file except in compliance with the License. You may obtain a *
 * copy of the License at http://www.apache.org/licenses/LICENSE-2.0.        *
 *                                                                           *
 * --------------------------------------------------------------------------*/
#ifndef PARTICLE_GENERATOR_NETWORK_H
#define PARTICLE_GENERATOR_NETWORK_H

#include "sph_data_containers.h"
#include "base_particle_generator.h"
#include "complex_body.h"

namespace SPH
{
    class ParticleGeneratorNetwork : public ParticleGenerator
    {
    public:
        ParticleGeneratorNetwork(SPHBody &sph_body, const Vecd &starting_pnt, const Vecd &second_pnt, int iterator, Real grad_factor);
        virtual ~ParticleGeneratorNetwork(){};

        virtual void initializeGeometricVariables() override;

    protected:
        Vecd starting_pnt_;                                 
        Vecd second_pnt_;                                   
        size_t n_it_;                                       
        bool fascicles_;                                    
        size_t segments_in_branch_;                         
        Real segment_length_;                               
        Real angle_ = 0.3;                                  
        Real repulsivity_ = 0.175;                          
        Real grad_factor_;                                  
        std::vector<Real> fascicle_angles_ = {-1.25, 0.75}; 
        Real fascicle_ratio_ = 15.0;                        
        SPHBody &sph_body_;
        Shape &body_shape_;
        BaseCellLinkedList *cell_linked_list_;
        TreeBody *tree_;
        Vecd getGradientFromNearestPoints(Vecd pt, Real delta);
        bool createABranchIfValid(size_t parent_id, Real angle, Real repulsivity, size_t number_segments);
        Vecd createATentativeNewBranchPoint(Vecd init_point, Vecd dir);
        bool isCollision(const Vecd &new_point, const ListData &nearest_neighbor, size_t parent_id);
        virtual bool extraCheck(const Vecd &new_point) { return false; };

        void growAParticleOnBranch(TreeBody::Branch *branch, const Vecd &new_point, const Vecd &end_direction);
    };
}
#endif // PARTICLE_GENERATOR_NETWORK_H