free_surface.h

Go to the documentation of this file.

/*
  Copyright (C) 2011 - 2018 by the authors of the ASPECT code.

  This file is part of ASPECT.

  ASPECT is free software; you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation; either version 2, or (at your option)
  any later version.

  ASPECT is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with ASPECT; see the file LICENSE.  If not see
  <http://www.gnu.org/licenses/>.
*/


#ifndef _aspect_free_surface_h
#define _aspect_free_surface_h

#include <aspect/simulator.h>
#include <aspect/simulator/assemblers/interface.h>

namespace aspect
{
  using namespace dealii;

  namespace Assemblers
  {
    template <int dim>
    class ApplyStabilization: public Assemblers::Interface<dim>,
      public SimulatorAccess<dim>
    {
      public:
        virtual ~ApplyStabilization () {};

        virtual
        void
        execute (internal::Assembly::Scratch::ScratchBase<dim>   &scratch,
                 internal::Assembly::CopyData::CopyDataBase<dim> &data) const;
    };
  }

  template<int dim>
  class FreeSurfaceHandler
  {
    public:
      FreeSurfaceHandler(Simulator<dim> &, ParameterHandler &prm);

      ~FreeSurfaceHandler();

      void execute();

      void set_assemblers();

      void setup_dofs();

      static
      void declare_parameters (ParameterHandler &prm);

      void parse_parameters (ParameterHandler &prm);


      double get_stabilization_term () const;

    private:
      void make_constraints ();

      void project_velocity_onto_boundary (LinearAlgebra::Vector &output);

      void compute_mesh_displacements ();

      void interpolate_mesh_velocity ();

      Simulator<dim> &sim;

      const FESystem<dim> free_surface_fe;

      DoFHandler<dim> free_surface_dof_handler;


      double free_surface_theta;

      LinearAlgebra::BlockVector mesh_velocity;

      LinearAlgebra::Vector mesh_displacements;

      LinearAlgebra::Vector fs_mesh_velocity;

      IndexSet mesh_locally_owned;

      IndexSet mesh_locally_relevant;

      ConstraintMatrix mesh_displacement_constraints;

      ConstraintMatrix mesh_vertex_constraints;

      struct SurfaceAdvection
      {
        enum Direction { normal, vertical };
      };

      typename SurfaceAdvection::Direction advection_direction;


      std::set<types::boundary_id> tangential_mesh_boundary_indicators;

      friend class Simulator<dim>;
      friend class SimulatorAccess<dim>;
  };
}


#endif