ASPECT
diffusion.h
Go to the documentation of this file.
1 /*
2  Copyright (C) 2020 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 
66  void
67  compute_velocity_constraints_on_boundary(const DoFHandler<dim> &mesh_deformation_dof_handler,
68  AffineConstraints<double> &mesh_velocity_constraints,
69  const std::set<types::boundary_id> &boundary_id) const override;
70 
74  static
76 
80  void parse_parameters (ParameterHandler &prm) override;
81 
82  private:
88  void diffuse_boundary (const DoFHandler<dim> &free_surface_dof_handler,
89  const IndexSet &mesh_locally_owned,
90  const IndexSet &mesh_locally_relevant,
91  LinearAlgebra::Vector &output,
92  const std::set<types::boundary_id> &boundary_id) const;
93 
100  void check_diffusion_time_step (const DoFHandler<dim> &mesh_deformation_dof_handler,
101  const std::set<types::boundary_id> &boundary_ids) const;
102 
108  double diffusivity;
109 
114 
120 
126  };
127  }
128 }
129 
130 
131 #endif
void declare_parameters(ParameterHandler &prm)
unsigned int timesteps_between_diffusion
Definition: diffusion.h:113
std::set< types::boundary_id > additional_tangential_mesh_boundary_indicators
Definition: diffusion.h:125
unsigned int boundary_id