cutdg_operator.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/fe/fe_system.h>
 
#include <meltpooldg/compressible_flow/data_types.hpp>
#include <meltpooldg/compressible_flow/explicit_time_integration_utils.hpp>
#include <meltpooldg/compressible_flow/kernels.hpp>
#include <meltpooldg/compressible_flow/operation_scratch_data.hpp>
#include <meltpooldg/utilities/dg_generic_convection_diffusion_worker.hpp>
 
namespace MeltPoolDG::CompressibleFlow
{
  template <int dim, typename number, bool is_viscous = true>
  class CutDGOperator
  {
  public:
    using VectorType            = dealii::LinearAlgebra::distributed::Vector<number>;
    using MappingInfoType       = CutUtil::MappingInfoType<dim, number>;
    using MappingInfoVectorType = CutUtil::MappingInfoVectorType<dim, number>;
 
    using ConservedVariables         = ConservedVariablesType<dim, number>;
    using ConservedVariablesGradient = ConservedVariablesGradientType<dim, number>;
 
    using ConvectionDiffusionOperator =
      Utils::DGConvectionDiffusionOperator<dim,
                                           number,
                                           ConvectiveFlux<dim, number>,
                                           DiffusiveFlux<dim, number>>;
 
    using ConvectionOperator =
      Utils::DGConvectionOperator<dim, number, ConvectiveFlux<dim, number>>;
 
    explicit CutDGOperator(OperationScratchData<dim, number> &flow_scratch_data,
                           const MappingInfoType             &mapping_info_surface_in,
                           const MappingInfoVectorType       &mapping_info_cells_in,
                           const MappingInfoVectorType       &mapping_info_faces_in);
 
    void
    set_inflow_field_unfitted_boundary(std::shared_ptr<dealii::Function<dim>> &inflow_function);
 
    void
    set_unfitted_object_velocity(std::shared_ptr<dealii::Function<dim>> &velocity_function);
 
    template <CellEvaluatorType<dim, dim + 2, number, dealii::VectorizedArray<number>> Integrator>
    void
    do_cell_integral_rhs(
      Integrator                                                    &phi,
      const OperationScratchData<dim, number>                       &flow_scratch_data,
      const dealii::Tensor<1, dim, dealii::VectorizedArray<number>> *constant_body_force,
      const unsigned int                                             q) const;
 
    void
    do_surface_integral_rhs(
      dealii::FEPointEvaluation<dim + 2, dim, dim, dealii::VectorizedArray<number>> &phi,
      const dealii::VectorizedArray<number>   &interior_penalty_parameter,
      const OperationScratchData<dim, number> &flow_scratch_data,
      const unsigned int                       q) const;
 
    template <FaceEvaluatorType<dim, dim + 2, number, dealii::VectorizedArray<number>> Integrator>
    void
    do_face_integral_rhs(Integrator                              &phi_m,
                         Integrator                              &phi_p,
                         const dealii::VectorizedArray<number>   &interior_penalty_parameter,
                         const OperationScratchData<dim, number> &flow_scratch_data,
                         const unsigned int                       q) const;
 
    template <FaceEvaluatorType<dim, dim + 2, number, dealii::VectorizedArray<number>> Integrator>
    void
    do_boundary_face_integral_rhs(Integrator                            &phi,
                                  const dealii::VectorizedArray<number> &interior_penalty_parameter,
                                  const OperationScratchData<dim, number> &flow_scratch_data,
                                  const auto                               boundary_id,
                                  const unsigned int                       q) const;
 
    void
    local_apply_cell_rhs(const dealii::MatrixFree<dim, number>       &matrix_free,
                         VectorType                                  &dst,
                         const VectorType                            &src,
                         const std::pair<unsigned int, unsigned int> &cell_range) const;
 
    void
    local_apply_face_rhs(const dealii::MatrixFree<dim, number>       &matrix_free,
                         VectorType                                  &dst,
                         const VectorType                            &src,
                         const std::pair<unsigned int, unsigned int> &face_range) const;
 
    void
    local_apply_boundary_face_rhs(const dealii::MatrixFree<dim, number>       &matrix_free,
                                  VectorType                                  &dst,
                                  const VectorType                            &src,
                                  const std::pair<unsigned int, unsigned int> &face_range) const;
 
    void
    local_apply_cell_lhs(const dealii::MatrixFree<dim, number>       &matrix_free,
                         VectorType                                  &dst,
                         const VectorType                            &src,
                         const std::pair<unsigned int, unsigned int> &cell_range) const;
 
    void
    local_apply_face_lhs(const dealii::MatrixFree<dim, number>       &matrix_free,
                         VectorType                                  &dst,
                         const VectorType                            &src,
                         const std::pair<unsigned int, unsigned int> &face_range) const;
 
    void
    local_apply_boundary_face_lhs(const dealii::MatrixFree<dim, number>       &matrix_free,
                                  VectorType                                  &dst,
                                  const VectorType                            &src,
                                  const std::pair<unsigned int, unsigned int> &face_range) const;
 
    void
    create_rhs(const number     &time,
               const number     &time_step,
               VectorType       &dst,
               const VectorType &src) const;
 
    void
    vmult(VectorType &dst, const VectorType &src) const;
 
  private:
    OperationScratchData<dim, number> &flow_scratch_data;
 
    const MappingInfoType &mapping_info_surface;
 
    const MappingInfoVectorType &mapping_info_cells;
 
    const MappingInfoVectorType &mapping_info_faces;
 
    dealii::FESystem<dim> fe_point_temp;
 
    const unsigned int n_dofs_per_cell;
 
    mutable number current_time_step = 0.;
 
    mutable number current_inv_time_step = 0.;
 
    std::shared_ptr<dealii::Function<dim>> unfitted_object_velocity;
 
    std::shared_ptr<dealii::Function<dim>> unfitted_inflow;
 
    void
    get_adjacent_face_values_at_unfitted_boundary(
      const dealii::Point<dim, dealii::VectorizedArray<number>> &q_point,
      const ConservedVariables                                  &w_m,
      ConservedVariables                                        &w_p,
      const ConservedVariablesGradient                          &grad_w_m,
      ConservedVariablesGradient                                &grad_w_p) const;
  };
} // namespace MeltPoolDG::CompressibleFlow