## FCSys.Species.H2.Gas

H2 gas

### Information

Extends from Modelica.Icons.Package (Icon for standard packages).

### Package Content

Name | Description
--- | ---
Fixed | Fixed properties

## FCSys.Species.H2.Gas.Fixed

Fixed properties

### Information

Assumptions:

1. The generalized resistivities (η, θ) are fixed (e.g., independent of thermodynamic state).
2. Ideal gas

The default resistivities (η = `1/(89.6e-7*U.Pa*U.s)` and θ = `U.m*U.K/(183e-3*U.W)`) are based on data of H2 gas at 1 atm and 300 K from Incropera and DeWitt [Incropera2002, pp. 919–920]. Table 1 lists the properties at other temperatures.

Table 1: Properties of H2 gas at 1 atm [Incropera2002, pp. 919–920].
T `/U.K` | cp`*U.kg*U.K /(U.J*Data.m)` | η`*U.Pa*U.s` | θ`*U.W /(U.m*U.K)`
--- | --- | --- | ---
100 | 11.23e3 | 1/42.1e-7 | 1/67.0e-3
150 | 12.60e3 | 1/56.0e-7 | 1/101e-3
200 | 13.54e3 | 1/68.1e-7 | 1/131e-3
250 | 14.06e3 | 1/78.9e-7 | 1/157e-3
300 | 14.31e3 | 1/89.6e-7 | 1/183e-3
350 | 14.43e3 | 1/98.8e-7 | 1/204e-3
400 | 14.48e3 | 1/108.2e-7 | 1/226e-3
450 | 14.50e3 | 1/117.2e-7 | 1/247e-3
500 | 14.52e3 | 1/126.4e-7 | 1/266e-3
550 | 14.53e3 | 1/134.3e-7 | 1/285e-3
600 | 14.55e3 | 1/142.4e-7 | 1/305e-3
700 | 14.61e3 | 1/157.8e-7 | 1/342e-3
800 | 14.70e3 | 1/172.4e-7 | 1/378e-3
900 | 14.83e3 | 1/186.5e-7 | 1/412e-3
1000 | 14.99e3 | 1/201.3e-7 | 1/448e-3
1100 | 15.17e3 | 1/213.0e-7 | 1/488e-3
1200 | 15.37e3 | 1/226.2e-7 | 1/528e-3
1300 | 15.59e3 | 1/238.5e-7 | 1/568e-3
1400 | 15.81e3 | 1/250.7e-7 | 1/610e-3
1500 | 16.02e3 | 1/262.7e-7 | 1/655e-3
1600 | 16.28e3 | 1/273.7e-7 | 1/697e-3
1700 | 16.58e3 | 1/284.9e-7 | 1/742e-3
1800 | 16.96e3 | 1/296.1e-7 | 1/786e-3
1900 | 17.49e3 | 1/307.2e-7 | 1/835e-3
2000 | 18.25e3 | 1/318.2e-7 | 1/878e-3

The specific heat capacity is not fixed because it would affect the chemical potential and result in an incorrect cell potential.

For more information, please see the Species model.

Extends from Fluid (Base model for a fluid species).

### Parameters

Type | Name | Default | Description
--- | --- | --- | ---
Integer | n_inter | 0 | Number of exchange connections with other phases
Material properties
replaceable package Data | FCSys.Characteristics.H2.Gas… | Characteristic data
Independence factors
NumberAbsolute | k_intra_Phi[n_intra, n_trans] | ones(n_intra, n_trans) | For translational exchange among species within the phase [1]
NumberAbsolute | k_intra_Q[n_intra] | ones(n_intra) | For thermal exchange among species within the phase [1]
Initialization
Velocity | phi.start[n_trans] | 0 | Velocity [L/T]
Current | I.start[n_trans] | 0 | Current [N/T]
Velocity | phi_boundaries.start[n_trans, Side] | 0 | Normal velocities at the boundaries [L/T]
Force | f.start[n_trans] | 0 | Total normal translational force on pairs of boundaries [L.M/T2]
Force | minusDeltaf.start[n_trans] | 0 | Dynamic and nonequilibrium compression forces [L.M/T2]
Geometry
Length | kL[:] | L[cartTrans] | Effective transport length [L]
Initialization
Init | initMaterial | Init.pressure | Method of initializing the material state
Init | initEnergy | Init.temperature | Method of initializing the thermal state
Amount | N_IC |  | Initial amount of material [N]
Density | rho_IC |  | Initial density [N/L3]
Volume | V_IC |  | Initial volume [L3]
PressureAbsolute | p_IC | environment.p_dry | Initial pressure [M/(L.T2)]
TemperatureAbsolute | T_IC |  | Initial temperature [L2.M/(N.T2)]
Potential | h_IC |  | Initial specific enthalpy [L2.M/(N.T2)]
Potential | g_IC |  | Initial Gibbs potential [L2.M/(N.T2)]
Assumptions
Integer | n_trans | 1 | Number of transport axes
Integer | n_chem | 1 | Number of reaction and phase change processes
Formulation of the conservation equations
ConsThermo | consMaterial | ConsThermo.dynamic | Material
Boolean | consRot | false | Conserve rotational momentum
ConsTrans | consTransX | ConsTrans.dynamic | X-axis translational momentum
ConsTrans | consTransY | ConsTrans.dynamic | Y-axis translational momentum
ConsTrans | consTransZ | ConsTrans.dynamic | Z-axis translational momentum
ConsThermo | consEnergy | ConsThermo.dynamic | Energy
Axes with upstream discretization
Boolean | upstreamX | true | X
Boolean | upstreamY | true | Y
Boolean | upstreamZ | true | Z
Flow conditions
Boolean | approxVelocity | true | Calculate normal boundary velocities assuming uniform density
NumberAbsolute | Nu_Phi[Axis] | {4,4,4} | Translational Nusselt numbers [1]
NumberAbsolute | Nu_Q | 1 | Thermal Nusselt number [1]

### Connectors

Type | Name | Description
--- | --- | ---
Intra | intra[n_intra] | Connectors to exchange translational momentum and energy within the phase
Inter | inter[n_inter] | Connectors to exchange translational momentum and energy with all other species
Dalton | dalton | Connector for additivity of pressure
Boundary | boundaries[n_trans, Side] | Connectors for transport
Chemical | chemical[n_chem] | Connector for reactions and phase change
Material properties
replaceable package Data | Characteristic data

### Modelica definition

    model Fixed "Fixed properties"
      extends Fluid(
        redeclare replaceable package Data = FCSys.Characteristics.H2.Gas (b_v=[1],
              n_v={-1,0}),
        redeclare parameter Q.Mobility mu=Data.mu(),
        redeclare parameter Q.TimeAbsolute nu=Data.nu(),
        redeclare parameter Q.Continuity zeta=Data.zeta(),
        redeclare parameter Q.Fluidity eta=1/(8.96e-6*U.Pa*U.s),
        redeclare parameter Q.ResistivityThermal theta=U.m*U.K/(0.183*U.W),
        final tauprime,
        final N0,
        n_chem=1,
        p_IC=environment.p_dry);

      // See the documentation for a table of values.

    end Fixed;