FCSys.Species.O2.Gas

O₂ gas

Information

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

Package Content

Name	Description
Fixed	Fixed properties

---

FCSys.Species.O2.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/(207.2e-7*U.Pa*U.s) and θ = U.m*U.K/(26.8e-3*U.W)) are based on data of gas at 1 atm and 300 K from Incropera and DeWitt [Incropera2002, pp. 920–921]. Table 1 lists the properties at other temperatures.

Table 1: Properties of O₂ gas at 1 atm [Incropera2002, pp. 920–921]

T /U.K	cp*U.kg*U.K /(U.J*Data.m)	η *U.Pa*U.s	θ*U.W /(U.m*U.K)
100	0.962e3	1/76.4e-7	1/9.25e-3
150	0.921e3	1/114.8e-7	1/13.8e-3
200	0.915e3	1/147.5e-7	1/18.3e-3
250	0.915e3	1/178.6e-7	1/22.6e-3
300	0.920e3	1/207.2e-7	1/26.8e-3
350	0.929e3	1/233.5e-7	1/29.6e-3
400	0.942e3	1/258.2e-7	1/33.0e-3
450	0.956e3	1/281.4e-7	1/36.3e-3
500	0.972e3	1/303.3e-7	1/41.2e-3
550	0.988e3	1/324.0e-7	1/44.1e-3
600	1.003e3	1/343.7e-7	1/47.3e-3
700	1.031e3	1/380.8e-7	1/52.8e-3
800	1.054e3	1/415.2e-7	1/58.9e-3
900	1.074e3	1/447.2e-7	1/64.9e-3
1000	1.090e3	1/477.0e-7	1/71.0e-3
1100	1.103e3	1/505.5e-7	1/75.8e-3
1200	1.115e3	1/532.5e-7	1/81.9e-3
1300	1.125e3	1/588.4e-7	1/87.1e-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.O2.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_O2	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]
Advanced
Pressure	p_stop	-Modelica.Constants.inf	Pressure below which the simulation should terminate [M/(L.T2)]

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.O2.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/(20.72e-6*U.Pa*U.s),
    redeclare parameter Q.ResistivityThermal theta=U.m*U.K/(26.8e-3*U.W),
    final tauprime,
    final N0,
    n_chem=1,
    p_IC=environment.p_O2);

  // See the documentation for a table of values.

  parameter Q.Pressure p_stop=-Modelica.Constants.inf 
    "Pressure below which the simulation should terminate";

equation 
  if p_stop > 0 then
    when p < p_stop then
      terminate("There is no more " + Data.formula + ".");
    end when;
  end if;

end Fixed;