Checking the Conservation of Current, Energy, and Mass
In this part of the tutorial, you set up reports and generate monitors and plots. These plots check the balance of current, energy, and mass, respectively, throughout the solid oxide fuel cell (SOFC).
- Current Conservation
The current conservation plot compares the electric current on the anode and cathode current collector boundaries to the current on the anode and cathode triple phase boundaries (TPB).
- Energy Conservation
The energy that is generated from the SOFC is not entirely converted to external electrical power. Some energy is lost as heat at the outlet. The energy conservation plot balances the external power that is generated Eqn. (5286) with the flow boundary enthalpy losses Eqn. (5285).
The non-isothermal energy conservation equation is:(5283)where and are the temperatures at the inlet and outlet respectively. is the mass flux and is the mass fraction of species . The enthalpy of formation, of species , that is associated with the mass that is consumed by the reaction, is converted to internal electrochemical heat ( ), internal Ohmic heat ( ), and external electric power ( ):(5284)combining Eqn. (5283) and Eqn. (5284), gives:(5285)The electrical power that is available to the external circuit is calculated by the anode-cathode potential difference, , and the total current that is available to the external circuit, :(5286)Simcenter STAR-CCM+ calculates the ohmic heating by accounting for electric potential losses throughout every electrically or ionically conductive substrate:(5287)where is the electrical conductivity and is the electric potential gradient.Simcenter STAR-CCM+ calculates the electrochemical heat that is released from a reaction, , using:(5288)where is the overpotential and is the equilibrium potential temperature derivative. - Mass Conservation
The mass conservation plot compares the mass that is gained by the gas streams to the mass that is produced by electrochemical reactions at the TPB.
- Right-click the Reports node and select .
- Rename the Surface Integral 1 node to Current_Anode.
-
Select
Current_Anode and set the following properties:
Property Setting Field Function Boundary Specific Electric Current Parts - In the same way, create 11 more Surface Integral reports and name them according to the table in the next step.
-
Set the following properties for the reports:
Report Field Function Parts Current_AnodeTPB Current_Cathode Boundary Specific Electric Current Current_CathodeTPB Electric_Power ElectricPowerFlux Enthalpy_Gain Boundary Heat Flux Mass_Gain_H2 Mass Flux H2 Mass_Gain_H2O Mass Flux H2O Mass_Gain_O2 Mass Flux O2 Mass_Reaction_H2 Mass_Reaction_H2O Mass_Reaction_O2 -
In the same way as above, create a Volume Integral report
and set the following properties:
Property Setting Field Function Parts
-
To create reports and a plot for current conservation:
- Multi-select the Current_Anode, Current_AnodeTPB, Current_Cathode, and Current_CathodeTPB nodes.
- Right-click one of the selected nodes and select Create Monitor and Plot from Report.
- In the Create Plot From Reports dialog, click Single Plot.
- Rename the Reports Plot node to Current Conservation.
-
In the same way, create monitors and a plot for each of the following groups of reports:
- For Electric_Power and Enthalpy_Gain, create a plot named Energy Conservation.
- For Mass_Gain_H2, Mass_Gain_H2O, Mass_Gain_O2, Mass_Reaction_H2, Mass_Reaction_H2O, and Mass_Reaction_O2, create a plot named Mass Conservation.
-
Expand the Plots node and set the following
properties:
Node Property Setting Current Conservation Title Current Conservation Energy Conservation Title Energy Conservation
Axes
Left AxisMinimum 0.0 Maximum 0.01 Mass Conservation Title Mass Conservation
Axes
Left AxisMinimum -8.0E-10 Maximum 8.0E-10 -
Group specific monitors according to the Normalization Option.
-
Expand the Monitors node and set Normalization Option as follows:
Node Normalization Option Normalization Option: Auto Auto Normalization Option: Manual Manual Normalization Option: Off Off -
Expand the following nodes, multi-select each of the
Manual Normalization sub-nodes, and set
Normalization Factor to
-1.0:
- Current_Anode Monitor
- Current_CathodeTPB Monitor
- Enthalpy_Gain Monitor
- Mass_Gain_H2 Monitor
- Mass_Gain_H2O Monitor
- Mass_Gain_O2 Monitor
Setting the normalization factor to -1 inverts the monitor values.
The normalization factors are inverse of molecular weights—they are used to convert between molar flux and mass flux when checking conservation.
-
For the following monitor nodes, set
Normalization Factor of the
Manual Normalization sub-nodes as follows:
- Mass_Reaction_H2 Monitor: -0.496032
- Mass_Reaction_H2O Monitor: -0.0555093
- Mass_Reaction_O2 Monitor: -0.03125
- Save the simulation.