dem_time_integrators.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/timer.h>
 
#include <meltpooldg/particles/dem_util.hpp>
#include <meltpooldg/particles/particle_accessor.hpp>
#include <meltpooldg/particles/particle_iterator.hpp>
#include <meltpooldg/time_integration/solution_history.hpp>
 
#include <concepts>
#include <ranges>
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType, typename ContainerType>
    requires std::ranges::range<ContainerType> and
             std::same_as<std::ranges::range_value_t<ContainerType>,
                          DEMParticleAccessor<dim, number>>
  void
  symplectic_euler_advance_time_step(const number time_step, ContainerType particle_range)
  {
    for (DEMParticleAccessor<dim, number> &particle : particle_range)
      {
        for (unsigned d = 0; d < dim; ++d)
          {
            particle.linear_velocity(d) += time_step * particle.linear_acceleration(d);
            particle.get_location()[d] += time_step * particle.linear_velocity(d);
            particle.linear_acceleration(d) =
              particle.force(d) / particle.get_property(ObstacleType::Properties::mass);
          }
 
        for (unsigned d = 0; d < axial_dim<dim>; ++d)
          {
            particle.angular_velocity(d) += time_step * particle.angular_acceleration(d);
            particle.angular_acceleration(d) =
              particle.torque(d) /
              particle.get_property(ObstacleType::Properties::moment_of_inertia);
          }
      }
  }
} // namespace MeltPoolDG