darcy_damping_operation.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/base/aligned_vector.h>
#include <deal.II/base/vectorization.h>
 
#include <deal.II/lac/generic_linear_algebra.h>
#include <deal.II/lac/la_parallel_vector.h>
#include <deal.II/lac/vector.h>
 
#include <meltpooldg/core/material.hpp>
#include <meltpooldg/core/scratch_data.hpp>
#include <meltpooldg/flow/darcy_damping_data.hpp>
#include <meltpooldg/post_processing/generic_data_out.hpp>
 
 
namespace MeltPoolDG::Flow
{
  template <int dim, typename number>
  class DarcyDampingOperation
  {
  public:
    using VectorType = dealii::LinearAlgebra::distributed::Vector<number>;
 
    explicit DarcyDampingOperation(const DarcyDampingData<number>      &data_in,
                                   const ScratchData<dim, dim, number> &scratch_data,
                                   const unsigned int flow_vel_hanging_nodes_dof_idx,
                                   const unsigned int flow_quad_idx);
 
    void
    reinit();
 
    void
    set_darcy_damping_at_q(const Material<number> &material,
                           const VectorType       &ls_as_heaviside,
                           const VectorType       &temperature,
                           const unsigned int      ls_hanging_nodes_dof_idx,
                           const unsigned int      heat_dof_idx);
    void
    set_darcy_damping_at_q(const Material<number> &material,
                           const VectorType       &ls_as_heaviside,
                           const VectorType       &temperature,
                           const VectorType       *interface_temperature,
                           const unsigned int      ls_hanging_nodes_dof_idx,
                           const unsigned int      heat_cut_dof_idx,
                           const unsigned int      heat_cont_no_bc_dof_idx);
 
    void
    assemble_rhs(VectorType &force_rhs, const VectorType &velocity_vec, const bool zero_out = true);
 
    void
    attach_output_vectors(GenericDataOut<dim, number> &data_out) const;
 
    dealii::VectorizedArray<number> &
    get_damping(const unsigned int cell, const unsigned int q);
 
    const dealii::VectorizedArray<number> &
    get_damping(const unsigned int cell, const unsigned int q) const;
 
  private:
    const number mushy_zone_morphology;
 
    const number avoid_div_zero_constant;
 
    const ScratchData<dim, dim, number> &scratch_data;
 
    const unsigned int flow_vel_hanging_nodes_dof_idx;
 
    const unsigned int flow_quad_idx;
 
    mutable VectorType damping;
 
    mutable dealii::AlignedVector<dealii::VectorizedArray<number>> damping_at_q;
 
    mutable dealii::Vector<number> damping_output;
 
    dealii::VectorizedArray<number>
    compute_darcy_damping_coefficient(const dealii::VectorizedArray<number> &solid_fraction) const;
 
    dealii::AlignedVector<dealii::VectorizedArray<number>> &
    get_damping_at_q();
  };
} // namespace MeltPoolDG::Flow