Initialize at openlb-1-3

1 files changed, 110 insertions, 0 deletions

diff --git a/src/particles/hlbm/hlbmDynamics2D.h b/src/particles/hlbm/hlbmDynamics2D.h
new file mode 100644
index 0000000..afc7bb5
--- /dev/null
+++ b/src/particles/hlbm/hlbmDynamics2D.h
@@ -0,0 +1,110 @@
+/*  DESCRIPTOR Boltzmann sample, written in C++, using the OpenLB
+ *  library
+ *
+ *  Copyright (C) 2006-2016 Thomas Henn, Fabian Klemens, Robin Trunk, Davide Dapelo
+ *  E-mail contact: info@openlb.net
+ *  The most recent release of OpenLB can be downloaded at
+ *  <http://www.openlb.net/>
+ *
+ *  This program is free software; you can redistribute it and/or
+ *  modify it under the terms of the GNU General Public License
+ *  as published by the Free Software Foundation; either version 2
+ *  of the License, or (at your option) any later version.
+ *
+ *  This program is distributed in the hope that it will be useful,
+ *  but WITHOUT ANY WARRANTY; without even the implied warranty of
+ *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ *  GNU General Public License for more details.
+ *
+ *  You should have received a copy of the GNU General Public
+ *  License along with this program; if not, write to the Free
+ *  Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
+ *  Boston, MA  02110-1301, USA.
+ */
+
+#ifndef PARTICLEDYNAMICS_2D_H
+#define PARTICLEDYNAMICS_2D_H
+
+#include "core/superLattice2D.h"
+#include "functors/analytical/indicator/smoothIndicatorBaseF2D.h"
+
+
+namespace olb {
+
+template<typename T, typename DESCRIPTOR>
+class ParticleDynamics2D {
+private:
+  SuperLattice2D<T, DESCRIPTOR>& _sLattice;
+  UnitConverter<T,DESCRIPTOR> const& _converter;
+  SuperGeometry2D<T>& _superGeometry;
+  std::shared_ptr<IndicatorF2D<T> > _indicatorF;
+  std::vector<SmoothIndicatorF2D<T,T,true>* > _vectorOfIndicator;
+  T _lengthX;
+  T _lengthY;
+  Vector<T,2> _accExt;
+  bool _escapeFromDomain=false;
+  bool _oldWallCollision=true;
+
+public:
+  ParticleDynamics2D(SuperLattice2D<T, DESCRIPTOR>& sLattice,
+                     UnitConverter<T,DESCRIPTOR> const& converter,
+                     SuperGeometry2D<T>& superGeometry,
+                     T lengthX, T lengthY, Vector<T,2> accExt = Vector<T,2> (0.,0.))
+    : _sLattice(sLattice),
+      _converter(converter),
+      _superGeometry(superGeometry),
+      _lengthX(lengthX),
+      _lengthY(lengthY),
+      _accExt(accExt)
+  {}
+
+  ParticleDynamics2D(SuperLattice2D<T, DESCRIPTOR>& sLattice,
+                     UnitConverter<T,DESCRIPTOR> const& converter,
+                     SuperGeometry2D<T>& superGeometry, std::shared_ptr<IndicatorF2D<T> > indicatorF,
+                     bool escapeFromDomain, bool oldWallCollision=false,
+                     Vector<T,2> accExt = Vector<T,2> (0.,0.))
+    : _sLattice(sLattice),
+      _converter(converter),
+      _superGeometry(superGeometry),
+      _indicatorF(indicatorF),
+      _accExt(accExt),
+      _escapeFromDomain(escapeFromDomain),
+      _oldWallCollision(oldWallCollision)
+  {}
+
+  void addCircle(Vector< T, 2> center, T radius, T density, T epsilon,
+                 Vector<S,2> vel = Vector<S,2> (0.,0.));
+
+  void addCuboid(Vector< T, 2> center, T xLength, T yLength, T density,
+                 T epsilon, T theta=0, Vector<S,2> vel = Vector<S,2> (0.,0.));
+
+  void addTriangle(Vector< T, 2> center, T radius, T density, T epsilon,
+                   T theta, Vector<S,2> vel = Vector<S,2> (0.,0.));
+
+  void addParticle(SmoothIndicatorF2D<T, T, true>& indicator);
+
+  void computeBoundaryForce(std::vector<SmoothIndicatorF2D<T,T,true>* >& indicator);
+
+  void addWallColl(SmoothIndicatorF2D<T,T,true>& indicator, T delta);
+
+  void verletIntegration(SmoothIndicatorF2D<T,T,true>& indicator);
+
+  void updateParticleDynamics(std::string name, SmoothIndicatorF2D<T,T,true>& indicator);
+
+  void checkAndRemoveEscaped();
+
+  void addParticleField(SmoothIndicatorF2D<T,T,true>& indicator);
+
+  void simulateTimestep(std::string name);
+
+  void print();
+
+  void load(std::string filename, T epsilon);
+
+  void save(std::string filename);
+
+};
+
+}
+
+#endif