CircularPeriodicBoundary Class Reference

Defines a pair of periodic walls. The particles are in {x: position_left<=normal*x <position_right}, with normal being the outward unit normal vector of the right wall. If a particle moves outside these boundaries, it will be shifted. More...

#include <CircularPeriodicBoundary.h>

Inheritance diagram for CircularPeriodicBoundary:

Public Member Functions
	CircularPeriodicBoundary ()
	CircularPeriodicBoundary (double innerRadius)
CircularPeriodicBoundary *	copy () const
	BaseBoundary copy method. More...
void	rotateParticle (BaseParticle *P, double angle)
int	createPeriodicParticles (BaseParticle *P, ParticleHandler &pH)
bool	checkBoundaryAfterParticleMoved (BaseParticle *P, ParticleHandler &pH)
void	read (std::istream &is)
	reads wall More...
void	print (std::ostream &os) const
	outputs wall More...
Public Member Functions inherited from BaseBoundary
virtual	~BaseBoundary ()
virtual int	createPeriodicParticles (BaseParticle *P UNUSED, ParticleHandler &pH UNUSED)
virtual bool	checkBoundaryAfterParticleMoved (BaseParticle *P UNUSED, ParticleHandler &pH UNUSED)
virtual void	checkBoundaryActionsBeforeTimeStep (ParticleHandler &pH UNUSED, WallHandler &wH UNUSED, RNG &random UNUSED)
void	set_Index (int index)
void	set_Id (int id)
void	setHandler (BoundaryHandler *handler)
virtual void	read (std::istream &is UNUSED)=0
	reads boundary More...
virtual void	print (std::ostream &os UNUSED) const =0
	outputs boundary More...
virtual void	moveInHandler (int newPos)

Private Attributes
double	innerRadius

Detailed Description

Defines a pair of periodic walls. The particles are in {x: position_left<=normal*x <position_right}, with normal being the outward unit normal vector of the right wall. If a particle moves outside these boundaries, it will be shifted.

Definition at line 32 of file CircularPeriodicBoundary.h.

Constructor & Destructor Documentation

CircularPeriodicBoundary::CircularPeriodicBoundary ( )

inline

Definition at line 35 of file CircularPeriodicBoundary.h.

References innerRadius.

Referenced by copy().

                               : BaseBoundary()
    {
        innerRadius=1.0;
        #ifdef CONSTUCTOR_OUTPUT
            std::cerr << "CircularPeriodicBoundary() finished" << std::endl;
        #endif              
    }

CircularPeriodicBoundary::CircularPeriodicBoundary ( double  innerRadius )

inline

Definition at line 43 of file CircularPeriodicBoundary.h.

References innerRadius.

                                                 : BaseBoundary()
    {
        this->innerRadius=innerRadius;
        #ifdef CONSTUCTOR_OUTPUT
            std::cerr << "CircularPeriodicBoundary(double innerRadius) finished" << std::endl;
        #endif              
    }   

Member Function Documentation

bool CircularPeriodicBoundary::checkBoundaryAfterParticleMoved ( BaseParticle *  P, ParticleHandler &  pH  )

inline

Todo:

{TW: Dinant, please confirm that i and oldI should be integer}

Definition at line 122 of file CircularPeriodicBoundary.h.

References ParticleHandler::addObject(), BaseParticle::copy(), BaseParticle::get_Displacement(), BaseParticle::get_Index(), BaseParticle::get_Position(), innerRadius, constants::pi, R, BaseHandler< T >::removeObject(), rotateParticle(), BaseParticle::set_Displacement(), Vec3D::X, and Vec3D::Y.

    {
        double R=sqrt(pow(P->get_Position().X,2)+pow(P->get_Position().Y,2));
        double alpha=atan2(P->get_Position().Y,P->get_Position().X);
        int i=log(R/innerRadius)/log(2.0)+1;
        double pieSize=2.0/pow(2.0,i)*constants::pi;
        
        double oldR=sqrt(pow(P->get_Position().X-P->get_Displacement().X,2)+pow(P->get_Position().Y-P->get_Displacement().Y,2));
        int oldI=log(oldR/innerRadius)/log(2.0)+1.0;
        
        if(i>0&&i>oldI) //Particle moves outward so it may have to be deleted
        {
            //std::cout<<"Particle="<<P->get_Index()<<" moving outward with alpha="<<alpha<<" and pieSize="<<pieSize<<" ";
            if(alpha<0||alpha>pieSize)
            {
                if(alpha>2.0*pieSize)
                {
                    //std::cout<<"and it is rotated into the pie"<<std::endl;
                    rotateParticle(P,-pieSize);
                }
                else
                {
                    //Delete particle
                    //std::cout<<"and it should be deleted"<<std::endl;
                    pH.removeObject(P->get_Index());
                    return true;
                }
            }
            //else
                //std::cout<<"and nothing happens"<<std::endl;
        }
        else if (i>=0&&i<oldI) //Particle moves inward so it has to be coppied
        {
            //std::cout<<"Particle="<<P->get_Index()<<" moving inward and is thus coppied with alpha="<<alpha<<" and pieSize="<<pieSize<<std::endl;
            //std::cout<<"i="<<i<<" oldI="<<oldI<<" R="<<R<<" oldR="<<oldR<<std::endl;
            //std::cout<<"Position="<<P->get_Position()<<" Displacement="<<P->get_Displacement()<<std::endl;
            BaseParticle* F0=P->copy();
            F0->set_Displacement(Vec3D(0.0,0.0,0.0));
            if(alpha<0)
            {
                rotateParticle(P,pieSize);
                rotateParticle(F0,0.5*pieSize);
            }
            else if(alpha<0.5*pieSize)
            {
                rotateParticle(F0,0.5*pieSize);
            }
            else if(alpha<pieSize)
            {
                rotateParticle(F0,-0.5*pieSize);
            }
            else
            {
                rotateParticle(P,-pieSize);
                rotateParticle(F0,-0.5*pieSize);
            }
            pH.addObject(F0);
        }
        else if(i>0&&alpha<0)
        {
            //std::cout<<"Particle="<<P->get_Index()<<" i="<<i<<" R="<<R<<" alpha="<<alpha<<" positive rotated pieSize="<<pieSize<<std::endl;           
            rotateParticle(P,pieSize);
        }
        else if(i>0&&alpha>pieSize)
        {
            //std::cout<<"Particle="<<P->get_Index()<<" i="<<i<<" R="<<R<<" alpha="<<alpha<<" negative rotated pieSize="<<pieSize<<std::endl;           
            rotateParticle(P,-pieSize);
        }
        return false;
    }   

CircularPeriodicBoundary* CircularPeriodicBoundary::copy ( ) const

inlinevirtual

BaseBoundary copy method.

It calls the copy constructor of this BaseBoundary, useful for polymorphism todo{Does this work correctly?}

Implements BaseBoundary.

Definition at line 51 of file CircularPeriodicBoundary.h.

References CircularPeriodicBoundary().

    {
        #ifdef CONSTUCTOR_OUTPUT
            std::cerr << "virtual CircularPeriodicBoundary* copy() const finished" << std::endl;
        #endif              
        return new CircularPeriodicBoundary(*this);
    }

int CircularPeriodicBoundary::createPeriodicParticles ( BaseParticle *  P, ParticleHandler &  pH  )

inline

Definition at line 73 of file CircularPeriodicBoundary.h.

References ParticleHandler::addObject(), BaseParticle::copy(), BaseParticle::get_InteractionRadius(), BaseParticle::get_PeriodicFromParticle(), BaseParticle::get_Position(), ParticleHandler::getLargestParticle(), innerRadius, constants::pi, R, rotateParticle(), BaseParticle::set_periodicFromParticle(), Vec3D::X, and Vec3D::Y.

    {
        double R=sqrt(pow(P->get_Position().X,2)+pow(P->get_Position().Y,2));
        double alpha=atan2(P->get_Position().Y,P->get_Position().X);
        int i=log(R/innerRadius)/log(2.0)+1;
        double pieSize=2.0/pow(2.0,i)*constants::pi;
        int C=0;
        //std::cout<<"R="<<R<<" alpha="<<alpha<<" i="<<i<<" pieSize="<<pieSize<<std::endl;
        
        //Check if the particle is close to it's inner Radius or is is close to zero alpha (small y)
        if(i>0&&(R-(P->get_InteractionRadius()+pH.getLargestParticle()->get_InteractionRadius())<pow(2.0,i-1)*innerRadius||P->get_Position().Y<(P->get_InteractionRadius()+pH.getLargestParticle()->get_InteractionRadius())))
        {
            //std::cout<<"Going to shift because "<<R-P->get_Radius()<<"<"<<pow(2,i-1)*innerRadius<<" or "<<P->get_Position().Y<<"<"<<P->get_Radius()<<std::endl;
            //std::cout<<*P<<" has been shifted"<<std::endl;

            BaseParticle* F0=P->copy();
            rotateParticle(F0,pieSize);

            //If Particle is Mdouble shifted, get correct original particle         
            BaseParticle* From=P;
            while(From->get_PeriodicFromParticle()!=NULL)
                From=From->get_PeriodicFromParticle();      
            F0->set_periodicFromParticle(From);

            //std::cout<<*F0<<" is the result"<<std::endl;
            pH.addObject(F0);
            C++;
        }
        //Check here only for i>0 becuase for i=1 they both give the same particle
        if(i>1&&R*R*(1-pow(cos(alpha-pieSize),2))<P->get_InteractionRadius()+pH.getLargestParticle()->get_InteractionRadius())
        {
            //std::cout<<*P<<" has been shifted back"<<std::endl;

            BaseParticle* F0=P->copy();
            rotateParticle(F0,-pieSize);

            //If Particle is Mdouble shifted, get correct original particle         
            BaseParticle* From=P;
            while(From->get_PeriodicFromParticle()!=NULL)
                From=From->get_PeriodicFromParticle();      
            F0->set_periodicFromParticle(From);

            //std::cout<<*F0<<" is the result"<<std::endl;
            pH.addObject(F0);
            C++;
        }
        return C;
    }

void CircularPeriodicBoundary::print ( std::ostream &  os ) const

inline

outputs wall

Definition at line 201 of file CircularPeriodicBoundary.h.

References innerRadius.

                                     {
        os << "CircularPeriodicBoundary innerRadius "<<innerRadius;
    }

void CircularPeriodicBoundary::read ( std::istream &  is )

inline

reads wall

Definition at line 195 of file CircularPeriodicBoundary.h.

References innerRadius.

                               { 
    std::string dummy;
        is >>  dummy >> innerRadius;
    }

void CircularPeriodicBoundary::rotateParticle ( BaseParticle *  P, double  angle  )

inline

todo{Do we need to update rotations and rotational velocitys?}

Definition at line 59 of file CircularPeriodicBoundary.h.

References BaseParticle::get_Position(), BaseParticle::get_Velocity(), R, BaseParticle::set_Position(), BaseParticle::set_Velocity(), Vec3D::X, Vec3D::Y, and Vec3D::Z.

Referenced by checkBoundaryAfterParticleMoved(), and createPeriodicParticles().

    {
        double R=sqrt(pow(P->get_Position().X,2)+pow(P->get_Position().Y,2));
        double alphaPos=atan2(P->get_Position().Y,P->get_Position().X);
        double V=sqrt(pow(P->get_Velocity().X,2)+pow(P->get_Velocity().Y,2));
        double alphaVel=atan2(P->get_Velocity().Y,P->get_Velocity().X);
        alphaPos+=angle;
        alphaVel+=angle;
        
        P->set_Position(Vec3D(cos(alphaPos)*R,sin(alphaPos*R),P->get_Position().Z));
        P->set_Velocity(Vec3D(cos(alphaVel)*V,sin(alphaVel*V),P->get_Velocity().Z));
    }

Member Data Documentation

double CircularPeriodicBoundary::innerRadius

private

Definition at line 208 of file CircularPeriodicBoundary.h.

Referenced by checkBoundaryAfterParticleMoved(), CircularPeriodicBoundary(), createPeriodicParticles(), print(), and read().

---

The documentation for this class was generated from the following file:

• CircularPeriodicBoundary.h