cl_GM_HelmholtzTransport_Hydrogen.hpp

Go to the documentation of this file.

/*
 * BELFEM -- The Berkeley Lab Finite Element Framework
 * Copyright (c) 2026, The Regents of the University of California,
 * through Lawrence Berkeley National Laboratory (subject to receipt of any required
 * approvals from the U.S. Dept. of Energy).  All rights reserved.
 *
 * Developers: Christian Messe, Gregory Giard
 *
 * See the top-level LICENSE file for the complete license and disclaimer.
 */
 
#ifndef BELFEM_CL_GM_HELMHOLTZTRANSPORT_HYDROGEN_HPP
#define BELFEM_CL_GM_HELMHOLTZTRANSPORT_HYDROGEN_HPP
 
#include "typedefs.hpp"
#include "cl_Vector.hpp"
#include "en_Helmholtz.hpp"
#include "cl_GM_HelmholtzTransport.hpp"
 
namespace belfem
{
    namespace gasmodels
    {
//----------------------------------------------------------------------------

        class HelmholtzTransport_Hydrogen : public HelmholtzTransport
        {
//----------------------------------------------------------------------------
        private:
//----------------------------------------------------------------------------

            const real mM = 2.01588 ;

            const real mSigma = 0.297 ;

            const real mEpsilonKb = 30.41 ;

            const real mTcrit = 33.145 ;

            const real mRhoScale = 90.909090909 ;

            Vector< real > mA ;

            Vector< real > mB ;

            Vector< real > mC ;

            real mConstEta0 ;

            real mConstBeta ;
 
//- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
//  thermal conductivity, Assael et al. 2011
//- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

            const HelmholtzModel mModel ;

            real mTcritEoS ;
            real mPcritEoS ;
            real mRhocritEoS ;

            Vector< real > mLambdaA1 ;
            Vector< real > mLambdaA2 ;

            Vector< real > mLambdaB1 ;
            Vector< real > mLambdaB2 ;

            const real mQd    = 1.0 / 4.0e-10 ;  
            const real mXi0   = 1.5e-10 ;        
            const real mGammaC = 0.052 ;         
            real       mTrefC ;                  
 
//----------------------------------------------------------------------------
        public:
//----------------------------------------------------------------------------

            HelmholtzTransport_Hydrogen(
                    Gas & aParent,
                    const HelmholtzModel aModel );
 
            ~HelmholtzTransport_Hydrogen() = default ;
 
//----------------------------------------------------------------------------

            real
            mu( const real T, const real p ) const ;
 
//----------------------------------------------------------------------------

            real
            lambda( const real T, const real p ) const ;
 
//----------------------------------------------------------------------------
        private:
//----------------------------------------------------------------------------

            real
            eta_0( const real T ) const ;
 
//----------------------------------------------------------------------------

            real
            eta_1( const real T ) const ;
 
//----------------------------------------------------------------------------

            real
            lambda_0( const real T ) const ;
 
//----------------------------------------------------------------------------

            real
            lambda_e( const real T, const real rho ) const ;
 
//----------------------------------------------------------------------------

            real
            lambda_c( const real T, const real p, const real rho ) const ;
 
//----------------------------------------------------------------------------

            real
            chi( const real T, const real v ) const ;
 
//----------------------------------------------------------------------------
        } ;
 
//----------------------------------------------------------------------------
    }
}
#endif //BELFEM_CL_GM_HELMHOLTZTRANSPORT_HYDROGEN_HPP