external_force.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 EXTERNAL_FORCE_H
#define EXTERNAL_FORCE_H



#include "base_data_package.h"

namespace SPH {
    class ExternalForce
    {
    public:
        ExternalForce();
        virtual ~ExternalForce() {};
        virtual Vecd InducedAcceleration(Vecd& position) = 0;
    };

    class Gravity : public ExternalForce
    {
    protected:
        Vecd global_acceleration_;
        Vecd zero_potential_reference_;
    public:
        Gravity(Vecd gravity_vector, Vecd reference_position = Vecd(0));
        virtual ~Gravity() {};

        virtual Vecd InducedAcceleration(Vecd& position) override;
        Real getPotential(Vecd& position);
    };
}
#endif //EXTERNAL_FORCE_H