level_set_elem_base.h

Go to the documentation of this file.

/*
 * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
 * Copyright (C) 2013-2019  Manav Bhatia
 *
 * This library is free software; you can 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 2.1 of the License, or (at your option) any later version.
 *
 * This library 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
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */


#ifndef __mast__level_set_elem_base__
#define __mast__level_set_elem_base__

// MAST includes
#include "base/elem_base.h"

namespace MAST {
    
    // Forward declerations
    class ElementPropertyCardBase;
    class BoundaryConditionBase;
    class FEMOperatorMatrix;
    template <typename ValType> class FieldFunction;
    
    
    class LevelSetElementBase:
    public MAST::ElementBase
    {
    public:
        LevelSetElementBase(MAST::SystemInitialization&             sys,
                            MAST::AssemblyBase&                     assembly,
                            const MAST::GeomElem&                   elem);
        
        
        virtual ~LevelSetElementBase();
        
        
        void
        set_velocity_function(const MAST::FieldFunction<Real>& vel) {
            
            libmesh_assert(!_phi_vel);
            _phi_vel = &vel;
        }
        
        void set_propagation_mode(bool f) {
            _if_propagation = f;
        }

        
        void set_reference_solution_for_initialization(const RealVectorX& sol);
        
        
        virtual bool
        internal_residual (bool request_jacobian,
                           RealVectorX& f,
                           RealMatrixX& jac);
        
        
        virtual bool
        velocity_residual (bool request_jacobian,
                           RealVectorX& f,
                           RealMatrixX& jac_xdot,
                           RealMatrixX& jac);
        
        bool
        side_external_residual (bool request_jacobian,
                                RealVectorX& f,
                                RealMatrixX& jac,
                                std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
        
        bool
        volume_external_residual (bool request_jacobian,
                                  RealVectorX& f,
                                  RealMatrixX& jac,
                                  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
        
        virtual bool
        internal_residual_sensitivity (bool request_jacobian,
                                       RealVectorX& f,
                                       RealMatrixX& jac);
        virtual bool
        velocity_residual_sensitivity (bool request_jacobian,
                                       RealVectorX& f,
                                       RealMatrixX& jac);
        
        bool
        side_external_residual_sensitivity (bool request_jacobian,
                                            RealVectorX& f,
                                            RealMatrixX& jac,
                                            std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
        
        bool
        volume_external_residual_sensitivity (bool request_jacobian,
                                              RealVectorX& f,
                                              RealMatrixX& jac,
                                              std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
        
        Real
        volume();

        Real
        perimeter();

        Real
        perimeter_sensitivity();

        
        Real
        volume_boundary_velocity_on_side(unsigned int s);

    protected:

        void _velocity_and_source(const unsigned int qp,
                                  const libMesh::Point& p,
                                  const Real t,
                                  const MAST::FEMOperatorMatrix& Bmat,
                                  const std::vector<MAST::FEMOperatorMatrix>& dBmat,
                                  RealVectorX& vel,
                                  Real&        source);
        
        void _tau(const MAST::FEBase& fe,
                  unsigned int qp,
                  const MAST::FEMOperatorMatrix& Bmat,
                  const std::vector<MAST::FEMOperatorMatrix>& dBmat,
                  const RealVectorX& vel,
                  RealMatrixX& tau);
        
        
        
        
        void
        _dc_operator(const MAST::FEBase& fe,
                     const unsigned int qp,
                     const std::vector<MAST::FEMOperatorMatrix>& dB_mat,
                     const RealVectorX& vel,
                     Real& dc);
        
        void
        _calculate_dxidX (const MAST::FEBase& fe,
                          const unsigned int qp,
                          RealMatrixX& dxi_dX);
        
        void _initialize_fem_operators(const unsigned int qp,
                                       const MAST::FEBase& fe,
                                       MAST::FEMOperatorMatrix& Bmat,
                                       std::vector<MAST::FEMOperatorMatrix>& dBmat);
        
        
        const MAST::FieldFunction<Real>*  _phi_vel;

        bool  _if_propagation;
        
        RealVectorX    _ref_sol;
    };
    
}


#endif // __mast__level_set_elem_base__