ASPECT
diffusion.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2020 - 2022 by the authors of the ASPECT code.
3 
4  This file is part of ASPECT.
5 
6  ASPECT is free software; you can redistribute it and/or modify
7  it under the terms of the GNU General Public License as published by
8  the Free Software Foundation; either version 2, or (at your option)
9  any later version.
10 
11  ASPECT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU General Public License for more details.
15 
16  You should have received a copy of the GNU General Public License
17  along with ASPECT; see the file LICENSE. If not see
18  <http://www.gnu.org/licenses/>.
19 */
20 
21 
22 #ifndef _aspect_mesh_deformation_diffusion_h
23 #define _aspect_mesh_deformation_diffusion_h
24 
28 
29 
30 namespace aspect
31 {
32  using namespace dealii;
33 
34 
35  namespace MeshDeformation
36  {
42  template<int dim>
43  class Diffusion : public Interface<dim>, public SimulatorAccess<dim>
44  {
45  public:
46  Diffusion();
47 
52  void initialize() override;
53 
58  void update() override;
59 
60 
61 
68  void
69  compute_velocity_constraints_on_boundary(const DoFHandler<dim> &mesh_deformation_dof_handler,
70  AffineConstraints<double> &mesh_velocity_constraints,
71  const std::set<types::boundary_id> &boundary_id) const override;
72 
76  bool needs_surface_stabilization () const override;
77 
81  static
82  void declare_parameters (ParameterHandler &prm);
83 
87  void parse_parameters (ParameterHandler &prm) override;
88 
89  private:
95  void diffuse_boundary (const DoFHandler<dim> &free_surface_dof_handler,
96  const IndexSet &mesh_locally_owned,
97  const IndexSet &mesh_locally_relevant,
98  LinearAlgebra::Vector &output,
99  const std::set<types::boundary_id> &boundary_id) const;
100 
107  void check_diffusion_time_step (const DoFHandler<dim> &mesh_deformation_dof_handler,
108  const std::set<types::boundary_id> &boundary_ids) const;
109 
115  double diffusivity;
116 
121 
127 
133  };
134  }
135 }
136 
137 
138 #endif
::TrilinosWrappers::MPI::Vector Vector
Definition: global.h:223
void declare_parameters(ParameterHandler &prm)
unsigned int timesteps_between_diffusion
Definition: diffusion.h:120
std::set< types::boundary_id > additional_tangential_mesh_boundary_indicators
Definition: diffusion.h:132
Definition: compat.h:88