concentration_model.hh

Go to the documentation of this file.

/*!
 *
 * Copyright (C) 2007 Technical University of Liberec.  All rights reserved.
 *
 * Please make a following refer to Flow123d on your project site if you use the program for any purpose,
 * especially for academic research:
 * Flow123d, Research Centre: Advanced Remedial Technologies, Technical University of Liberec, Czech Republic
 *
 * This program is free software; you can redistribute it and/or modify it under the terms
 * of the GNU General Public License version 3 as published by the Free Software Foundation.
 *
 * This program 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 this program; if not,
 * write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 021110-1307, USA.
 *
 *
 * $Id$
 * $Revision$
 * $LastChangedBy$
 * $LastChangedDate$
 *
 * @file
 * @brief Discontinuous Galerkin method for equation of transport with dispersion.
 *  @author Jan Stebel
 */

#ifndef CONC_TRANS_MODEL_HH_
#define CONC_TRANS_MODEL_HH_

#include "advection_diffusion_model.hh"
#include "fields/multi_field.hh"




class ConcentrationTransportModel : public AdvectionDiffusionModel {
public:

    class ModelEqData : public TransportBase::TransportEqData {
    public:

        /// Boundary conditions (Dirichlet) for concentrations.
        BCField<3, FieldValue<3>::Vector> bc_conc;
        /// Initial concentrations.
        Field<3, FieldValue<3>::Vector> init_conc;
        /// Longitudal dispersivity (for each substance).
        Field<3, FieldValue<3>::Vector> disp_l;
        /// Transversal dispersivity (for each substance).
        Field<3, FieldValue<3>::Vector> disp_t;
        /// Molecular diffusivity (for each substance).
        Field<3, FieldValue<3>::Vector> diff_m;

        Field<3, FieldValue<3>::Scalar > rock_density;      ///< Rock matrix density.
        Field<3, FieldValue<3>::Vector > sorption_mult;     ///< Coefficient of linear sorption.


        MultiField<3, FieldValue<3>::Scalar> output_field;



        ModelEqData();

        static constexpr const char * name() { return "SoluteTransport"; }

        static string default_output_field() { return "conc"; }

        static IT::Selection &get_output_selection_input_type(const string &implementation, const string &description);

    };

protected:

    /// Derived class should implement getter for ModelEqData instance.
    virtual ModelEqData &data() = 0;

    /**
     * Create input type that can be passed to the derived class.
     * @param implementation String characterizing the numerical method, e.g. DG, FEM, FVM.
     * @param description    Comment used to describe the record key.
     * @return
     */
    static IT::Record &get_input_type(const string &implementation, const string &description);

    /**
     * Formula to calculate the dispersivity tensor.
     * @param velocity  Fluid velocity.
     * @param Dm        Molecular diffusivity.
     * @param alphaL    Longitudal dispersivity.
     * @param alphaT    Transversal dispersivity.
     * @param porosity  Porosity.
     * @param cross_cut Cross-section.
     * @param K         Dispersivity tensor (output).
     */
    void calculate_dispersivity_tensor(const arma::vec3 &velocity,
            double Dm,
            double alphaL,
            double alphaT,
            double porosity,
            double cross_cut,
            arma::mat33 &K);

    /// Indicator of change in advection vector field.
    bool flux_changed;

    /// Pointer to list of substances (needed e.g. for access to molar masses).
    SubstanceList *substances_;

    /// Density of liquid (a global constant).
    double solvent_density_;


public:

    ConcentrationTransportModel();

    static string balance_prefix() { return "mass"; }

    UnitSI balance_units();

    void init_from_input(const Input::Record &in_rec, SubstanceList &substances) override;

    void compute_mass_matrix_coefficient(const std::vector<arma::vec3 > &point_list,
            const ElementAccessor<3> &ele_acc,
            std::vector<double> &mm_coef) override;


    void compute_retardation_coefficient(const std::vector<arma::vec3 > &point_list,
            const ElementAccessor<3> &ele_acc,
            std::vector<std::vector<double> > &ret_coef) override;


    void compute_advection_diffusion_coefficients(const std::vector<arma::vec3 > &point_list,
            const std::vector<arma::vec3> &velocity,
            const ElementAccessor<3> &ele_acc,
            std::vector<std::vector<arma::vec3> > &ad_coef,
            std::vector<std::vector<arma::mat33> > &dif_coef) override;

    void compute_init_cond(const std::vector<arma::vec3> &point_list,
            const ElementAccessor<3> &ele_acc,
            std::vector< arma::vec > &init_values) override;

    void compute_dirichlet_bc(const std::vector<arma::vec3> &point_list,
            const ElementAccessor<3> &ele_acc,
            std::vector< arma::vec > &bc_values) override;

    void compute_source_coefficients(const std::vector<arma::vec3> &point_list,
                const ElementAccessor<3> &ele_acc,
                std::vector<arma::vec> &sources_conc,
                std::vector<arma::vec> &sources_density,
                std::vector<arma::vec> &sources_sigma) override;

    void compute_sources_sigma(const std::vector<arma::vec3> &point_list,
                const ElementAccessor<3> &ele_acc,
                std::vector<arma::vec> &sources_sigma) override;

    ~ConcentrationTransportModel() override;
};








#endif /* CONC_TRANS_MODEL_HH_ */