brinkman_penalization.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/exception_macros.h>
#include <deal.II/base/exceptions.h>
#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/brinkman_penalization_data.hpp>
#include <meltpooldg/fluid_structure_interaction/fluid_structure_interaction_util.hpp>
#include <meltpooldg/particles/cell_batch_particle_cache.hpp>
#include <meltpooldg/particles/obstacle_field.hpp>
#include <meltpooldg/utilities/matrix_free_util.hpp>
 
 
namespace MeltPoolDG
{
  template <int dim, typename number, typename ObstacleType>
  struct BrinkmanObstacleForce
  {
    using VectorType = dealii::LinearAlgebra::distributed::Vector<number>;
 
    BrinkmanObstacleForce(const VectorType                       &solution,
                          const MatrixFreeContext<dim, number>   &matrix_free,
                          const BrinkmanPenalizationData<number> &data);
 
    void
    add_load_to_obstacles(ObstacleField<dim, number, ObstacleType> &obstacle_field) const;
 
  private:
    const BrinkmanPenalizationData<number> brinkman_penalization_data;
 
    const MatrixFreeContext<dim, number> matrix_free;
 
    const VectorType &solution;
  };
 
 
  template <int dim, typename number, typename ObstacleType>
  struct BrinkmanPenalizationResidualContribution final
    : public CompressibleFlow::ExternalFlowForce<dim, number>
  {
    using ConservedVariablesType = CompressibleFlow::ConservedVariablesType<dim, number>;
 
  public:
    BrinkmanPenalizationResidualContribution(
      const BrinkmanPenalizationData<number> &brinkman_penalization_data,
      const std::shared_ptr<MatrixFreeCellBatchParticleCache<dim, number, ObstacleType>>
        &cell_cache);
 
    ConservedVariablesType
    value(number                                                     time_step_size,
          const unsigned int                                         cell_batch_id,
          const dealii::Point<dim, dealii::VectorizedArray<number>> &q_point,
          const ConservedVariablesType                              &w_q) override;
 
  private:
    const BrinkmanPenalizationData<number> brinkman_penalization_data;
 
    const std::shared_ptr<MatrixFreeCellBatchParticleCache<dim, number, ObstacleType>> cell_cache;
  };
 
 
  template <int dim, typename number, typename ObstacleType>
  struct BrinkmanPenalizationJacobianContribution final
    : public CompressibleFlow::ExternalFlowForceJacobian<dim, number>
  {
    using ConservedVariablesType = CompressibleFlow::ConservedVariablesType<dim, number>;
 
  public:
    BrinkmanPenalizationJacobianContribution(
      const BrinkmanPenalizationData<number> &brinkman_penalization_data,
      const std::shared_ptr<MatrixFreeCellBatchParticleCache<dim, number, ObstacleType>>
        &cell_cache);
 
    ConservedVariablesType
    value(number                                                     time_step_size,
          unsigned int                                               cell_batch_id,
          const dealii::Point<dim, dealii::VectorizedArray<number>> &q_point,
          const ConservedVariablesType                              &w_q,
          const ConservedVariablesType                              &delta_w_q) override;
 
  private:
    const BrinkmanPenalizationData<number> brinkman_penalization_data;
 
    const std::shared_ptr<MatrixFreeCellBatchParticleCache<dim, number, ObstacleType>> cell_cache;
  };
} // namespace MeltPoolDG