nlisomoisturemat.h

Go to the documentation of this file.

/*
 *
 *                 #####    #####   ######  ######  ###   ###
 *               ##   ##  ##   ##  ##      ##      ## ### ##
 *              ##   ##  ##   ##  ####    ####    ##  #  ##
 *             ##   ##  ##   ##  ##      ##      ##     ##
 *            ##   ##  ##   ##  ##      ##      ##     ##
 *            #####    #####   ##      ######  ##     ##
 *
 *
 *             OOFEM : Object Oriented Finite Element Code
 *
 *               Copyright (C) 1993 - 2025   Borek Patzak
 *
 *
 *
 *       Czech Technical University, Faculty of Civil Engineering,
 *   Department of Structural Mechanics, 166 29 Prague, Czech Republic
 *
 *  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 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
 *  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 nlisomoisturemat_h
#define nlisomoisturemat_h
 
#include "tm/Materials/isomoisturemat.h"
#include "floatarray.h"
#include "scalarfunction.h"
#include "function.h"


#define _IFT_NlIsoMoistureMaterial_Name "nlisomoisturemat"
#define _IFT_NlIsoMoistureMaterial_isothermtype "isothermtype"
#define _IFT_NlIsoMoistureMaterial_permeabilitytype "permeabilitytype"
#define _IFT_NlIsoMoistureMaterial_capillarytransporttype "capillarytransporttype"
#define _IFT_NlIsoMoistureMaterial_rhodry "rhodry"
#define _IFT_NlIsoMoistureMaterial_capa "capa"
#define _IFT_NlIsoMoistureMaterial_hx "hx"
#define _IFT_NlIsoMoistureMaterial_dx "dx"
#define _IFT_NlIsoMoistureMaterial_iso_offset "isooffset"
#define _IFT_NlIsoMoistureMaterial_iso_h "iso_h"
#define _IFT_NlIsoMoistureMaterial_iso_wh "iso_w(h)"
#define _IFT_NlIsoMoistureMaterial_dd "dd"
#define _IFT_NlIsoMoistureMaterial_wf "wf"
#define _IFT_NlIsoMoistureMaterial_b "b"
#define _IFT_NlIsoMoistureMaterial_uh "uh"
#define _IFT_NlIsoMoistureMaterial_a "a"
#define _IFT_NlIsoMoistureMaterial_nn "nn"
#define _IFT_NlIsoMoistureMaterial_c "c"
#define _IFT_NlIsoMoistureMaterial_k "k"
#define _IFT_NlIsoMoistureMaterial_vm "vm"
#define _IFT_NlIsoMoistureMaterial_perm_h "perm_h"
#define _IFT_NlIsoMoistureMaterial_perm_ch "perm_c(h)"
#define _IFT_NlIsoMoistureMaterial_hc "hc"
#define _IFT_NlIsoMoistureMaterial_alpha0 "alpha0"
#define _IFT_NlIsoMoistureMaterial_c1 "c1"
#define _IFT_NlIsoMoistureMaterial_n "n"
#define _IFT_NlIsoMoistureMaterial_alphah "alphah"
#define _IFT_NlIsoMoistureMaterial_betah "betah"
#define _IFT_NlIsoMoistureMaterial_gammah "gammah"
#define _IFT_NlIsoMoistureMaterial_rhoh2o "rhoh2o"
#define _IFT_NlIsoMoistureMaterial_capperm_h "capperm_h"
#define _IFT_NlIsoMoistureMaterial_capperm_dwh "capperm_dw(h)"
#define _IFT_NlIsoMoistureMaterial_capperm_wv "capperm_wv"
#define _IFT_NlIsoMoistureMaterial_capperm_dwwv "capperm_dw(wv)"
#define _IFT_NlIsoMoistureMaterial_abs "abs"
#define _IFT_NlIsoMoistureMaterial_pl "pl"
#define _IFT_NlIsoMoistureMaterial_mu "mu"
#define _IFT_NlIsoMoistureMaterial_timescale "timescale"
#define _IFT_NlIsoMoistureMaterial_wn "wn"
#define _IFT_NlIsoMoistureMaterial_alpha "alpha"
#define _IFT_NlIsoMoistureMaterial_capil_coef "capil_coef"
#define _IFT_NlIsoMoistureMaterial_t "t"
#define _IFT_NlIsoMoistureMaterial_ttf "ttf"
#define _IFT_NlIsoMoistureMaterial_vg_b "vg_b"
#define _IFT_NlIsoMoistureMaterial_vg_m "vg_m"
 
namespace oofem {
class NlIsoMoistureMaterial : public IsotropicMoistureTransferMaterial
{
protected:
    enum isothermType { linear, multilinear, Ricken, Kuenzel, Hansen, BSB, bilinear, vanGenuchten } Isotherm;

    double rhodry = 0.;

    double moistureCapacity = 0.;

    FloatArray iso_h;
    FloatArray iso_wh;

    double dd = 0.;

    double wf = 0., b = 0.;

    double uh = 0., A = 0., nn = 0.;

    double c = 0., k = 0., Vm = 0.;

    double hx = 0., dx = 0.;
    double iso_offset = 0.;
    double c1 = 0., c2 = 0., capa2 = 0.;

    double vG_b = 0., vG_m = 0.;

    double wn = 0.;

    ScalarFunction alpha;
 
    enum permeabilityType { multilin, Bazant, Xi, KunzelPerm } Permeability;

    FloatArray perm_h;
    FloatArray perm_ch;

    double C1 = 0., n = 0., alpha0 = 0., hC = 0.;

    double alphah = 0., betah = 0., gammah = 0.;
 
    enum capillaryTransportType { Multilin_h, Multilin_wV, KunzelCT } CapillaryTransport;
    double Abs = 0.;

    double mu = 0.;

    double PL = 101325.;

    double timeScale = 1.;

    double capillary_transport_coef = 1000.;

    double T = 0.;
    int T_TF = 0;
 

    FloatArray capPerm_h;
    FloatArray capPerm_Dwh;
    FloatArray capPerm_wV;
    FloatArray capPerm_DwwV;
    double rhoH2O = 0.;
 
public:
    NlIsoMoistureMaterial(int n, Domain *d) : IsotropicMoistureTransferMaterial(n, d) { }
 
    void initializeFrom(InputRecord &ir) override;

    double giveMoistureCapacity(GaussPoint *gp, TimeStep *tStep) const override;
    double giveMoistureContent(double humidity) const override;
    double givePermeability(GaussPoint *gp, TimeStep *tStep) const override;
    double computeCapTranspCoeff(double humidity) const;

    double computeVaporDiffusionCoeff(GaussPoint *gp, TimeStep *tStep) const;
    double computeSaturationWaterVaporPressure(GaussPoint *gp, TimeStep *tStep) const;
    double computeTemperatureEffectOnViscosity(GaussPoint *gp, TimeStep *tStep) const;
    double giveTemperature(GaussPoint *gp, TimeStep *tStep) const;
 
    const char *giveInputRecordName() const override { return _IFT_NlIsoMoistureMaterial_Name; }
    const char *giveClassName() const override { return "NlIsoMoistureMaterial"; }
 
    double giveHumidity(GaussPoint *gp, ValueModeType mode) const override;
 
 
    bool hasInternalSource() const override;
    void computeInternalSourceVector(FloatArray &val, GaussPoint *gp, TimeStep *tStep, ValueModeType mode) const override;
};
} // end namespace oofem
#endif // nlisomoisturemat_h