stokes_law.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/point.h>
#include <deal.II/base/vectorization.h>
 
#include <deal.II/lac/la_parallel_vector.h>
 
#include <deal.II/matrix_free/matrix_free.h>
 
#include <meltpooldg/compressible_flow/data_types.hpp>
#include <meltpooldg/compressible_flow/utils.hpp>
#include <meltpooldg/fluid_structure_interaction/fluid_structure_interaction_data.hpp>
#include <meltpooldg/particles/obstacle_field.hpp>
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType>
  struct StokesLawSphericalParticleForce
  {
    using VectorType = dealii::LinearAlgebra::distributed::Vector<number>;
 
    StokesLawSphericalParticleForce(const VectorType                     &solution,
                                    const MatrixFreeContext<dim, number> &matrix_free,
                                    const number                          dynamic_viscosity);
 
    void
    add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const;
 
  private:
    const MatrixFreeContext<dim, number> matrix_free;
 
    const VectorType &solution;
 
    const number dynamic_viscosity;
  };
 
  template <int dim, typename number, typename ObstacleType>
  struct StokesLawFluidForce final : public CompressibleFlow::ExternalFlowForce<dim, number>
  {
    using ConservedVariablesType = CompressibleFlow::ConservedVariablesType<dim, number>;
    using VectorType             = dealii::LinearAlgebra::distributed::Vector<number>;
 
  public:
    StokesLawFluidForce(const VectorType                         &solution,
                        ObstacleField<dim, number, ObstacleType> &obstacle_handler,
                        const MatrixFreeContext<dim, number>     &matrix_free,
                        const number                              dynamic_viscosity);
 
    ConservedVariablesType
    value(number                                                     time_step_size,
          unsigned int                                               cell_batch_id,
          const dealii::Point<dim, dealii::VectorizedArray<number>> &q_point,
          const ConservedVariablesType                              &w_q) override;
 
  private:
    const VectorType &solution;
 
    const number dynamic_viscosity;
 
    const MatrixFreeContext<dim, number> matrix_free;
 
    ObstacleField<dim, number, ObstacleType> &obstacle_handler;
  };
 
} // namespace MeltPoolDG