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).

You display the results on the following plots:
  • 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:
    Figure 1. EQUATION_DISPLAY
    Σ i , I N H ^ i ( T I N ) Y i m ˙ Σ i , O U T H ^ i ( T O U T ) Y i m ˙ + Σ i , Re a c t i o n H ^ i ( T R ) m ˙ i + Q O + Q E = 0
    (5283)
    where T I N and T O U T are the temperatures at the inlet and outlet respectively. m ˙ is the mass flux and Y i is the mass fraction of species i . The enthalpy of formation, H ^ i of species i , that is associated with the mass that is consumed by the reaction, is converted to internal electrochemical heat ( Q E ), internal Ohmic heat ( Q O ), and external electric power ( P ):
    Figure 2. EQUATION_DISPLAY
    Σ i , Re a c t i o n H ^ i ( T R ) m ˙ i + Q O + Q E + P = 0
    (5284)
    combining Eqn. (5283) and Eqn. (5284), gives:
    Figure 3. EQUATION_DISPLAY
    P = Σ i , I N H ^ i ( T I N ) Y i m ˙ Σ i , O U T H ^ i ( T O U T ) Y i m ˙
    (5285)
    The electrical power that is available to the external circuit is calculated by the anode-cathode potential difference, U , and the total current that is available to the external circuit, I :
    Figure 4. EQUATION_DISPLAY
    P = I U
    (5286)
    Simcenter STAR-CCM+ calculates the ohmic heating by accounting for electric potential losses throughout every electrically or ionically conductive substrate:
    Figure 5. EQUATION_DISPLAY
    Q O = σ | ϕ | 2 d V
    (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, i, using:
    Figure 6. EQUATION_DISPLAY
    QE=i(η+TEeqoxT|P)
    (5288)
    where η is the overpotential and EeqoxT|P 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.

To create the reports:
  1. Right-click the Reports node and select New > User > Surface Integral.
  2. Rename the Surface Integral 1 node to Current_Anode.
  3. Select Current_Anode and set the following properties:
    Property Setting
    Field Function Boundary Specific Electric Current
    Parts Regions > currentcollector-anode > Boundaries > contact
  4. In the same way, create 11 more Surface Integral reports and name them according to the table in the next step.
  5. Set the following properties for the reports:
    Report Field Function Parts
    Current_AnodeTPB Boundary Specific Electric Current > Boundary Specific Electric Current of conductor-anode Regions > anode > Boundaries > mem-Interface [anode-TPB]
    Current_Cathode Boundary Specific Electric Current Regions > currentcollector-cathode > Boundaries > contact
    Current_CathodeTPB Boundary Specific Electric Current > Boundary Specific Electric Current of conductor-cathode Regions > cathode > Boundaries > mem-Interface [cathode-TPB]
    Electric_Power ElectricPowerFlux
    • Regions > currentcollector-anode > Boundaries > contact
    • Regions > currentcollector-cathode > Boundaries > contact
    Enthalpy_Gain Boundary Heat Flux
    • Regions > air > Boundaries > inlet-air
    • Regions > air > Boundaries > outlet-air
    • Regions > fuel > Boundaries > inlet-hydrogen
    • Regions > fuel > Boundaries > outlet-hydrogen
    Mass_Gain_H2 Mass Flux H2
    • Regions > fuel > Boundaries > inlet-hydrogen
    • Regions > fuel > Boundaries > outlet-hydrogen
    Mass_Gain_H2O Mass Flux H2O
    • Regions > fuel > Boundaries > inlet-hydrogen
    • Regions > fuel > Boundaries > outlet-hydrogen
    Mass_Gain_O2 Mass Flux O2
    • Regions > air > Boundaries > inlet-air
    • Regions > air > Boundaries > outlet-air
    Mass_Reaction_H2 Boundary Species Electrochemical Reaction Flux > Boundary Species Electrochemical Reaction Flux of H2 Regions > anode > Boundaries > mem-Interface [anode-TPB]
    Mass_Reaction_H2O Boundary Species Electrochemical Reaction Flux > Boundary Species Electrochemical Reaction Flux of H2O Regions > anode > Boundaries > mem-Interface [anode-TPB]
    Mass_Reaction_O2 Boundary Species Electrochemical Reaction Flux > Boundary Species Electrochemical Reaction Flux of O2 Regions > cathode > Boundaries > mem-Interface [cathode-TPB]
  6. In the same way as above, create a Volume Integral report and set the following properties:
    Property Setting
    Field Function Boundary Specific Electric Current > Boundary Specific Electric Current of conductor-anode
    Parts
    • Regions > anode
    • Regions > cathode
Create monitors and plots for the reports:
  1. To create reports and a plot for current conservation:
    1. Multi-select the Current_Anode, Current_AnodeTPB, Current_Cathode, and Current_CathodeTPB nodes.
    2. Right-click one of the selected nodes and select Create Monitor and Plot from Report.
    3. In the Create Plot From Reports dialog, click Single Plot.
    4. Rename the Reports Plot node to Current Conservation.
  2. 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.
  3. Expand the Plots node and set the following properties:
    Node Property Setting
    Current Conservation Title Current Conservation
    Energy Conservation Title Energy Conservation


    Axes


    Left Axis
    Minimum 0.0
    Maximum 0.01
    Mass Conservation Title Mass Conservation


    Axes


    Left Axis
    Minimum -8.0E-10
    Maximum 8.0E-10
  4. Group specific monitors according to the Normalization Option.
    1. Right-click the Monitors node and select New Group.
    2. Repeat the above step to create two further groups.
    3. Re-name the groups:
      • Normalization Option: Auto
      • Normalization Option: Manual
      • Normalization Option: Off
    4. Using the multi-select option, drag and drop the monitors into the groups as specified:
      Group Monitor
      Normalization Option: Auto Continuity
      Electric Potential
      Energy
      H2
      H2O
      N2
      O2
      X-momentum
      Y-momentum
      Z-momentum
      Normalization Option: Manual Current_Anode Monitor
      Current_CathodeTPB Monitor
      Enthalpy_Gain Monitor
      Mass_Gain_H2 Monitor
      Mass_Gain_H2O Monitor
      Mass_Gain_O2 Monitor
      Mass_Reaction_H2 Monitor
      Mass_Reaction_H2O Monitor
      Mass_Reaction_O2 Monitor
      Normalization Option: Off Current_AnodeTPB Monitor
      Current_Cathode Monitor
      Electric Power Monitor
  5. 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
  6. 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.

  1. 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
  2. Save the simulation.