linear_visco_elasticity.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (C) 2022 by Luca Heltai
//
// This file is part of the FSI-suite platform, based on the deal.II library.
//
// The FSI-suite platform is free software; you can use it, redistribute it,
// and/or modify it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 3.0 of the License,
// or (at your option) any later version. The full text of the license can be
// found in the file LICENSE at the top level of the FSI-suite platform
// distribution.
//
// ---------------------------------------------------------------------
#ifndef pdes_linear_visco_elasticity_h
#define pdes_linear_visco_elasticity_h
 
#include "parsed_tools/constants.h"
#include "parsed_tools/mapping_eulerian.h"
#include "pdes/linear_problem.h"
 
namespace PDEs
{
  using namespace dealii;
 
  template <int dim, int spacedim = dim, class LacType = LAC::LAdealii>
  class LinearViscoElasticity : public LinearProblem<dim, spacedim, LacType>
  {
  public:
    LinearViscoElasticity();
 
    virtual ~LinearViscoElasticity() = default;
 
    using ScratchData =
      typename LinearProblem<dim, spacedim, LacType>::ScratchData;
 
    using CopyData = typename LinearProblem<dim, spacedim, LacType>::CopyData;
 
    using VectorType =
      typename LinearProblem<dim, spacedim, LacType>::VectorType;
 
    void
    postprocess();
 
  protected:
    virtual void
    assemble_system() override;
 
    virtual void
    solve() override;
 
    const FEValuesExtractors::Vector displacement;
 
    const FEValuesExtractors::Vector velocity;
 
    ParsedTools::Constants constants_0;
 
    ParsedTools::Constants constants_1;
 
    std::set<types::material_id> material_ids_0;
 
    std::set<types::material_id> material_ids_1;
 
    ParsedTools::MappingEulerian<dim, spacedim> eulerian_mapping;
 
    double dt;
 
    double current_time;
 
    unsigned int current_cycle;
 
    typename LacType::BlockVector current_displacement;
    typename LacType::BlockVector current_displacement_locally_relevant;
  };
 
  namespace MPI
  {
    template <int dim, int spacedim = dim>
    using LinearViscoElasticity =
      PDEs::LinearViscoElasticity<dim, spacedim, LAC::LATrilinos>;
  }
 
  namespace Serial
  {
    template <int dim, int spacedim = dim>
    using LinearViscoElasticity =
      PDEs::LinearViscoElasticity<dim, spacedim, LAC::LAdealii>;
  }
} // namespace PDEs
#endif