Running: ./testmodel.py --libraries=/home/hudson/saved_omc/libraries/.openmodelica/libraries/ --ompython_omhome=/usr Chemical_Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.conf.json loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/ModelicaServices 4.0.0+maint.om/package.mo", uses=false) loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/Complex 4.0.0+maint.om/package.mo", uses=false) loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/Modelica 4.0.0+maint.om/package.mo", uses=false) loadFile("/home/hudson/saved_omc/libraries/.openmodelica/libraries/Chemical 1.4.1-master/package.mo", uses=false) Using package Chemical with version 1.4.1 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/Chemical 1.4.1-master/package.mo) Using package Modelica with version 4.0.0 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/Modelica 4.0.0+maint.om/package.mo) Using package Complex with version 4.0.0 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/Complex 4.0.0+maint.om/package.mo) Using package ModelicaServices with version 4.0.0 (/home/hudson/saved_omc/libraries/.openmodelica/libraries/ModelicaServices 4.0.0+maint.om/package.mo) Running command: translateModel(Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential,tolerance=1e-06,outputFormat="empty",numberOfIntervals=5000,variableFilter="",fileNamePrefix="Chemical_Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential") translateModel(Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential,tolerance=1e-06,outputFormat="empty",numberOfIntervals=5000,variableFilter="",fileNamePrefix="Chemical_Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential") Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/ModelicaServices 4.0.0+maint.om/package.mo): time 0.0009818/0.0009818, allocations: 107.1 kB / 16.37 MB, free: 6.363 MB / 14.72 MB Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/Complex 4.0.0+maint.om/package.mo): time 0.0009496/0.0009496, allocations: 191.5 kB / 17.31 MB, free: 5.91 MB / 14.72 MB Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/Modelica 4.0.0+maint.om/package.mo): time 1.251/1.251, allocations: 222.9 MB / 241 MB, free: 15.19 MB / 206.1 MB Notification: Performance of loadFile(/home/hudson/saved_omc/libraries/.openmodelica/libraries/Chemical 1.4.1-master/package.mo): time 0.07995/0.07995, allocations: 15.31 MB / 306.5 MB, free: 15.46 MB / 254.1 MB Notification: Performance of FrontEnd - Absyn->SCode: time 2.058e-05/2.06e-05, allocations: 2.281 kB / 371.7 MB, free: 50.76 MB / 302.1 MB [/home/hudson/saved_omc/libraries/.openmodelica/libraries/Chemical 1.4.1-master/package.mo:4740:7-4740:66:writable] Warning: An inner declaration for outer component system could not be found and was automatically generated. [/home/hudson/saved_omc/libraries/.openmodelica/libraries/Chemical 1.4.1-master/package.mo:4740:7-4740:66:writable] Notification: The diagnostics message for the missing inner is: Your model is using an outer "system" component but an inner "system" component is not defined. For simulation drag Modelica.Fluid.System into your model to specify system properties. Notification: Performance of NFInst.instantiate(Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential): time 0.02742/0.02746, allocations: 28.01 MB / 399.7 MB, free: 34.48 MB / 302.1 MB Notification: Performance of NFInst.instExpressions: time 0.02146/0.04894, allocations: 28.27 MB / 428 MB, free: 10.75 MB / 302.1 MB Notification: Performance of NFInst.updateImplicitVariability: time 0.001924/0.05092, allocations: 23.88 kB / 428 MB, free: 10.72 MB / 302.1 MB Notification: Performance of NFTyping.typeComponents: time 0.001432/0.05236, allocations: 0.6711 MB / 428.7 MB, free: 10.07 MB / 302.1 MB Notification: Performance of NFTyping.typeBindings: time 0.003498/0.05587, allocations: 1.311 MB / 430 MB, free: 8.758 MB / 302.1 MB Notification: Performance of NFTyping.typeClassSections: time 0.00909/0.06498, allocations: 4.262 MB / 434.2 MB, free: 4.469 MB / 302.1 MB Notification: Performance of NFFlatten.flatten: time 0.00252/0.06752, allocations: 2.449 MB / 436.7 MB, free: 2.016 MB / 302.1 MB Notification: Performance of NFFlatten.resolveConnections: time 0.0009453/0.06849, allocations: 1.134 MB / 437.8 MB, free: 0.8672 MB / 302.1 MB Notification: Performance of NFEvalConstants.evaluate: time 0.0008185/0.06931, allocations: 1.052 MB / 438.9 MB, free: 15.81 MB / 318.1 MB Notification: Performance of NFSimplifyModel.simplify: time 0.0009926/0.07032, allocations: 1.082 MB / 439.9 MB, free: 14.73 MB / 318.1 MB Notification: Performance of NFPackage.collectConstants: time 0.0002135/0.07054, allocations: 168 kB / 440.1 MB, free: 14.56 MB / 318.1 MB Notification: Performance of NFFlatten.collectFunctions: time 0.006766/0.07731, allocations: 2.676 MB / 442.8 MB, free: 11.88 MB / 318.1 MB Notification: Performance of combineBinaries: time 0.002061/0.0794, allocations: 2.926 MB / 445.7 MB, free: 8.926 MB / 318.1 MB Notification: Performance of replaceArrayConstructors: time 0.0009578/0.08036, allocations: 1.753 MB / 447.5 MB, free: 7.156 MB / 318.1 MB Notification: Performance of NFVerifyModel.verify: time 0.000401/0.08077, allocations: 291.2 kB / 447.7 MB, free: 6.871 MB / 318.1 MB Notification: Performance of FrontEnd: time 0.0002732/0.08105, allocations: 83.8 kB / 447.8 MB, free: 6.789 MB / 318.1 MB Notification: Model statistics after passing the front-end and creating the data structures used by the back-end: * Number of equations: 376 (362) * Number of variables: 376 (355) Notification: Performance of Bindings: time 0.006735/0.08779, allocations: 8.777 MB / 456.6 MB, free: 13.74 MB / 334.1 MB Notification: Performance of FunctionAlias: time 0.0009797/0.08878, allocations: 1.173 MB / 457.8 MB, free: 12.56 MB / 334.1 MB Notification: Performance of Early Inline: time 0.008278/0.09707, allocations: 8.631 MB / 466.4 MB, free: 3.867 MB / 334.1 MB Notification: Performance of simplify1: time 0.0004667/0.09756, allocations: 0.4916 MB / 466.9 MB, free: 3.375 MB / 334.1 MB Notification: Performance of Alias: time 0.008961/0.1065, allocations: 8.51 MB / 475.4 MB, free: 9.91 MB / 350.1 MB Notification: Performance of simplify2: time 0.0004263/0.107, allocations: 323.7 kB / 475.7 MB, free: 9.594 MB / 350.1 MB Notification: Performance of Events: time 0.0003191/0.1073, allocations: 305.2 kB / 476 MB, free: 9.289 MB / 350.1 MB Notification: Performance of Detect States: time 0.0006832/0.108, allocations: 0.7912 MB / 476.8 MB, free: 8.48 MB / 350.1 MB Notification: Performance of Partitioning: time 0.00116/0.1092, allocations: 1.248 MB / 478.1 MB, free: 6.926 MB / 350.1 MB Error: Internal error NBSlice.fillDependencyArray failed because number of flattened indices 1 for dependency electrodeReaction1.h_mix could not be devided by the body size 3 without rest. Error: Internal error NBAdjacency.Matrix.createPseudo failed for: [ARRY] (3) electrodeReaction1.products.h_outflow = electrodeReaction1.h_mix * fill(1.0, 3) ($RES_SIM_116) Error: Internal error NBAdjacency.Matrix.create failed to create adjacency matrix for system: System Variables (81/102) *************************** (1) [ALGB] (1) protected Real electrone1.molarVolume (min = 0.0) (2) [ALGB] (1) protected Real electrone1.u0 (3) [ALGB] (1) Real PbSO4_.a (4) [ALGB] (1) Real $FUN_13 (5) [ALGB] (1) Real HSO4.port_a.u (6) [ALGB] (1) Real $FUN_12 (7) [ALGB] (3) stream Real[3] electrodeReaction1.products.h_outflow (nominal = {2e4 for $i1 in 1:3}) (8) [ALGB] (1) Real $FUN_11 (9) [ALGB] (1) Real $FUN_10 (10) [ALGB] (1) flow Real HSO4.port_a.q (11) [ALGB] (1) protected Real electrodeReaction1.du (12) [ALGB] (1) Real H2O.x (min = 0.0, max = 1.0) (13) [ALGB] (1) protected Real electrone.pressure (14) [ALGB] (2) flow Real[2] electrodeReaction.products.q (15) [ALGB] (1) protected Real electrone.u0 (16) [ALGB] (2) Real[2] electrodeReaction.products.u (17) [ALGB] (1) Real PbSO4.port_a.u (18) [ALGB] (1) Real H.port_a.u (19) [ALGB] (2) stream Real[2] electrodeReaction1.substrates.h_outflow (nominal = {2e4 for $i1 in 1:2}) (20) [ALGB] (1) protected Real electrone1.pressure (21) [ALGB] (1) flow Real PbSO4.port_a.q (22) [ALGB] (1) flow Real H.port_a.q (23) [ALGB] (1) Real electrone1.solution.v (24) [ALGB] (1) Real $FUN_9 (25) [ALGB] (1) Real electrodeReaction.rr (start = 0.0) (26) [ALGB] (1) Real Pb.a (27) [ALGB] (1) Real $FUN_8 (28) [ALGB] (1) protected Real electrone1.molarEntropyPure (29) [ALGB] (1) Real $FUN_7 (30) [ALGB] (1) Real $FUN_6 (31) [ALGB] (1) Real H2O.a (32) [ALGB] (1) Real $FUN_5 (33) [ALGB] (1) Real $FUN_4 (34) [ALGB] (1) Real electrodeReaction1.h_mix (nominal = 2e4) (35) [ALGB] (1) Real $FUN_3 (36) [ALGB] (1) Real $FUN_2 (37) [ALGB] (1) protected Real electrone.molarVolume (min = 0.0) (38) [ALGB] (1) Real H.a (39) [ALGB] (1) Real $FUN_1 (40) [ALGB] (1) stream Real electrone.port_a.h_outflow (nominal = 2e4) (41) [ALGB] (1) Real electrodeReaction.h_mix (nominal = 2e4) (42) [DER-] (1) Real $DER.electrone1.temperature (43) [ALGB] (1) protected Real electrone.molarEntropyPure (44) [ALGB] (1) flow Real electrone.port_a.q (45) [ALGB] (1) Real electrone.solution.v (46) [ALGB] (1) Real electrone.port_a.u (47) [ALGB] (4) stream Real[4] electrodeReaction.substrates.h_outflow (nominal = {2e4 for $i1 in 1:4}) (48) [ALGB] (1) stream Real electrone1.port_a.h_outflow (nominal = 2e4) (49) [ALGB] (1) Real electrodeReaction1.rr (start = 0.0) (50) [ALGB] (3) Real[3] electrodeReaction1.products.u (51) [ALGB] (1) flow Real electrone1.port_a.q (52) [DER-] (1) Real $DER.electrone.temperature (53) [ALGB] (1) protected Real H2O.SelfClustering_K = exp(-(-(81.6348 + 32.845554 * H2O.Temperature)) / (H2O.Temperature * 8.31446261815324)) (min = 0.0, max = 1.0) (54) [ALGB] (2) stream Real[2] electrodeReaction.products.h_outflow (nominal = {2e4 for $i1 in 1:2}) (55) [ALGB] (3) flow Real[3] electrodeReaction1.products.q (56) [ALGB] (1) Real electrone1.port_a.u (57) [ALGB] (1) flow Real Pb.port_a.q (58) [ALGB] (1) Real Pb.port_a.u (59) [ALGB] (1) Real voltageSensor.v (60) [ALGB] (1) Real electrone.a (61) [ALGB] (1) Real PbSO4_.port_a.u (62) [ALGB] (1) flow Real PbSO4_.port_a.q (63) [ALGB] (1) flow Real H2O.port_a.q (64) [ALGB] (1) protected Real electrone.uPure (65) [ALGB] (1) protected Real electrodeReaction.du (66) [ALGB] (1) Real H2O.port_a.u (67) [ALGB] (1) protected Real electrone1.uPure (68) [ALGB] (1) Real PbO2.port_a.u (69) [ALGB] (1) Real HSO4.a (70) [DISC] (1) Boolean $SEV_3 (71) [DISC] (1) Boolean $SEV_2 (72) [DISC] (1) Boolean $SEV_1 (73) [ALGB] (1) flow Real PbO2.port_a.q (74) [ALGB] (2) flow Real[2] electrodeReaction1.substrates.q (75) [DISC] (1) Boolean $SEV_0 (76) [ALGB] (4) flow Real[4] electrodeReaction.substrates.q (77) [ALGB] (2) Real[2] electrodeReaction1.substrates.u (78) [ALGB] (4) Real[4] electrodeReaction.substrates.u (79) [ALGB] (1) Real electrone1.a (80) [ALGB] (1) Real PbO2.a (81) [ALGB] (1) Real PbSO4.a System Equations (88/102) *************************** (1) [SCAL] (1) electrodeReaction.du = electrodeReaction.p * electrodeReaction.products.u - electrodeReaction.s * electrodeReaction.substrates.u ($RES_SIM_132) (2) [SCAL] (1) voltageSensor.v = electrone1.solution.v - electrone.solution.v ($RES_SIM_133) (3) [SCAL] (1) HSO4.port_a.u = (8.31446261815324 * HSO4.Temperature * $FUN_2 + ((-885750.0) - HSO4.Temperature * ((-445.68170384034886) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.HSO4.stateOfMatter.molarVolumePure(Chemical.Substances.HydrogenSulfate_aqueous(0.097, -1.0, -752870.0, -885750.0, 1.0, 0.0, {"http://www.mhhe.com/physsci/chemistry/chang7/ssg/graphics/chang7/pdf/cng7pa08.pdf"}, 0.0, 0.0, false, 0.0, 0.0, 1800.0), HSO4.Temperature, HSO4.Pressure, HSO4.ElectricPotential, HSO4.MoleFractionBasedIonicStrength) / HSO4.Temperature) * ((-1e5) + HSO4.Pressure)))) - 96485.33212331001 * HSO4.ElectricPotential ($RES_SIM_171) (4) [SCAL] (1) $FUN_13 = log(PbSO4.a) ($RES_$AUX_379) (5) [SCAL] (1) H.port_a.u = 8.31446261815324 * H.Temperature * $FUN_4 + H.Temperature * ((Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.H.stateOfMatter.molarVolumePure(Chemical.Substances.Proton_aqueous(0.001007, 1.0, 0.0, 0.0, 1.0, 0.0, {"http://www.vias.org/genchem/standard_enthalpies_table.html"}, 0.0, 0.0, false, 0.0, 0.0, 1000.0), H.Temperature, H.Pressure, H.ElectricPotential, H.MoleFractionBasedIonicStrength) / H.Temperature) * ((-1e5) + H.Pressure)) + 96485.33212331001 * H.ElectricPotential ($RES_SIM_137) (6) [SCAL] (1) electrone.u0 = (electrone.substanceData.DfH + ((-298.15) + electrone.temperature) * electrone.substanceData.Cp) - electrone.temperature * ((electrone.substanceData.Cp * log(0.0033540164346805303 * electrone.temperature) + 0.0033540164346805303 * (electrone.substanceData.DfH - electrone.substanceData.DfG)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone.stateOfMatter.molarVolumePure(electrone.substanceData, electrone.temperature, electrone.pressure, electrone.solution.v, 0.0) / electrone.temperature) * ((-1e5) + electrone.pressure)) ($RES_SIM_90) (7) [SCAL] (1) HSO4.a = 1.0 ($RES_SIM_172) (8) [SCAL] (1) H.a = 1.0 ($RES_SIM_138) (9) [SCAL] (1) electrone.molarEntropyPure = (electrone.substanceData.Cp * log(0.0033540164346805303 * electrone.temperature) + 0.0033540164346805303 * (electrone.substanceData.DfH - electrone.substanceData.DfG)) - ((electrone.substanceData.MolarWeight / electrone.substanceData.density) / electrone.temperature) * ((-1e5) + electrone.pressure) ($RES_SIM_91) (10) [SCAL] (1) electrone1.port_a.u = 8.31446261815324 * electrone1.temperature * $FUN_10 + ((electrone1.substanceData.DfH + ((-298.15) + electrone1.temperature) * electrone1.substanceData.Cp) - electrone1.temperature * ((electrone1.substanceData.Cp * log(0.0033540164346805303 * electrone1.temperature) + 0.0033540164346805303 * (electrone1.substanceData.DfH - electrone1.substanceData.DfG)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone1.stateOfMatter.molarVolumePure(electrone1.substanceData, electrone1.temperature, electrone1.pressure, electrone1.solution.v, 0.0) / electrone1.temperature) * ((-1e5) + electrone1.pressure))) + 96485.33212331001 * electrone1.substanceData.z * electrone1.solution.v ($RES_SIM_56) (11) [SCAL] (1) H2O.SelfClustering_K = exp((0.12027235504272604 * (81.6348 + 32.845554 * H2O.Temperature)) / H2O.Temperature) ($RES_$AUX_377) (12) [SCAL] (1) electrone.port_a.h_outflow = (electrone.substanceData.DfH + ((-298.15) + electrone.temperature) * electrone.substanceData.Cp) + 96485.33212331001 * electrone.substanceData.z * electrone.solution.v ($RES_SIM_92) (13) [SCAL] (1) electrone1.a = electrone1.substanceData.gamma ($RES_SIM_57) (14) [SCAL] (1) $DER.electrone1.temperature = 0.0 ($RES_SIM_213) (15) [SCAL] (1) electrone1.pressure = cathode.BasePressure ($RES_SIM_216) (16) [SCAL] (1) electrodeReaction1.substrates[1].q + Pb.port_a.q = 0.0 ($RES_SIM_290) (17) [SCAL] (1) H2O.port_a.u = 8.31446261815324 * H2O.Temperature * $FUN_12 + (((-285830.0) + 75.3 * ((-298.15) + H2O.Temperature)) - H2O.Temperature * (((-196.54536307227906) + 75.3 * log(0.0033540164346805303 * H2O.Temperature)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.H2O.stateOfMatter.molarVolumePure(Chemical.Substances.Water_liquid(0.018015, 0.0, -227230.0, -285830.0, 1.0, 75.3, {"http://www.vias.org/genchem/standard_enthalpies_table.html"}, 0.0, 0.0, true, -81.6348, 32.845554, 1000.0), H2O.Temperature, H2O.Pressure, H2O.ElectricPotential, H2O.MoleFractionBasedIonicStrength) / H2O.Temperature) * ((-1e5) + H2O.Pressure))) ($RES_SIM_22) (18) [SCAL] (1) H2O.a = H2O.x ($RES_SIM_23) (19) [SCAL] (1) electrone1.molarVolume = electrone1.substanceData.MolarWeight / electrone1.substanceData.density + Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone1.stateOfMatter.molarVolumePure(electrone1.substanceData, electrone1.temperature, electrone1.pressure, electrone1.solution.v, 0.0) * log(Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone1.stateOfMatter.activityCoefficient(electrone1.substanceData, electrone1.temperature, electrone1.pressure, electrone1.solution.v, 0.0)) ($RES_SIM_60) (20) [SCAL] (1) electrone1.uPure = (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone1.stateOfMatter.molarEnthalpyElectroneutral(electrone1.substanceData, electrone1.temperature, electrone1.pressure, electrone1.solution.v, 0.0) - electrone1.temperature * Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone1.stateOfMatter.molarEntropyPure(electrone1.substanceData, electrone1.temperature, electrone1.pressure, electrone1.solution.v, 0.0)) + 96485.33212331001 * electrone1.substanceData.z * electrone1.solution.v ($RES_SIM_61) (21) [SCAL] (1) electrone1.u0 = (electrone1.substanceData.DfH + ((-298.15) + electrone1.temperature) * electrone1.substanceData.Cp) - electrone1.temperature * ((electrone1.substanceData.Cp * log(0.0033540164346805303 * electrone1.temperature) + 0.0033540164346805303 * (electrone1.substanceData.DfH - electrone1.substanceData.DfG)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone1.stateOfMatter.molarVolumePure(electrone1.substanceData, electrone1.temperature, electrone1.pressure, electrone1.solution.v, 0.0) / electrone1.temperature) * ((-1e5) + electrone1.pressure)) ($RES_SIM_62) (22) [SCAL] (1) 96485.33212331001 * electrone.substanceData.z * electrone.port_a.q = 0.0 ($RES_SIM_109) (23) [SCAL] (1) electrone1.molarEntropyPure = (electrone1.substanceData.Cp * log(0.0033540164346805303 * electrone1.temperature) + 0.0033540164346805303 * (electrone1.substanceData.DfH - electrone1.substanceData.DfG)) - ((electrone1.substanceData.MolarWeight / electrone1.substanceData.density) / electrone1.temperature) * ((-1e5) + electrone1.pressure) ($RES_SIM_63) (24) [SCAL] (1) electrone1.port_a.h_outflow = (electrone1.substanceData.DfH + ((-298.15) + electrone1.temperature) * electrone1.substanceData.Cp) + 96485.33212331001 * electrone1.substanceData.z * electrone1.solution.v ($RES_SIM_64) (25) [SCAL] (1) Pb.port_a.u = 8.31446261815324 * Pb.Temperature * $FUN_1 + (26.4 * ((-298.15) + Pb.Temperature) - Pb.Temperature * (26.4 * log(0.0033540164346805303 * Pb.Temperature) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.Pb.stateOfMatter.molarVolumePure(Chemical.Substances.Lead_solid(0.2072, 0.0, 0.0, 0.0, 1.0, 26.4, {"http://www.vias.org/genchem/standard_enthalpies_table.html, http://www.update.uu.se/~jolkkonen/pdf/CRC_TD.pdf"}, 0.0, 0.0, false, 0.0, 0.0, 1000.0), Pb.Temperature, Pb.Pressure, Pb.ElectricPotential, Pb.MoleFractionBasedIonicStrength) / Pb.Temperature) * ((-1e5) + Pb.Pressure))) ($RES_SIM_188) (26) [SCAL] (1) Pb.a = 1.0 ($RES_SIM_189) (27) [SCAL] (1) PbO2.port_a.q + electrodeReaction.substrates[1].q = 0.0 ($RES_SIM_262) (28) [SCAL] (1) electrone1.port_a.q + electrodeReaction.substrates[4].q = 0.0 ($RES_SIM_264) (29) [-IF-] (1)if $SEV_0 then (29) [----] [SCAL] (1) (electrodeReaction1.h_mix * electrodeReaction1.products.q) / fill(1.0, 3) + electrodeReaction1.substrates.q * {26.4 * ((-298.15) + Pb.Temperature), ($OMC$PositiveMax(-electrodeReaction.substrates[2].q, 1e-7) * electrodeReaction.substrates[2].h_outflow + $OMC$PositiveMax(-HSO4.port_a.q, 1e-7) * ((-885750.0) - 96485.33212331001 * HSO4.ElectricPotential)) / ($OMC$PositiveMax(-electrodeReaction.substrates[2].q, 1e-7) + $OMC$PositiveMax(-HSO4.port_a.q, 1e-7))} = 0.0 ($RES_SIM_113) (29) [----] elseif $SEV_1 then (29) [----] [SCAL] (1) (electrodeReaction1.h_mix * electrodeReaction1.substrates.q) / fill(1.0, 2) + electrodeReaction1.products.q * {(-918400.0) + 103.2 * ((-298.15) + PbSO4_.Temperature), ($OMC$PositiveMax(-electrodeReaction.substrates[3].q, 1e-7) * electrodeReaction.substrates[3].h_outflow + $OMC$PositiveMax(-H.port_a.q, 1e-7) * (96485.33212331001 * H.ElectricPotential)) / ($OMC$PositiveMax(-electrodeReaction.substrates[3].q, 1e-7) + $OMC$PositiveMax(-H.port_a.q, 1e-7)), electrone.port_a.h_outflow} = 0.0 ($RES_SIM_114) (29) [----] else (29) [----] [SCAL] (1) electrodeReaction1.h_mix = 0.0 ($RES_SIM_115) (29) [----] end if; (30) [ARRY] (3) electrodeReaction1.products.h_outflow = electrodeReaction1.h_mix * fill(1.0, 3) ($RES_SIM_116) (31) [ARRY] (2) electrodeReaction1.substrates.h_outflow = electrodeReaction1.h_mix * fill(1.0, 2) ($RES_SIM_117) (32) [ARRY] (3) electrodeReaction1.rr * electrodeReaction1.p = -electrodeReaction1.products.q ($RES_SIM_118) (33) [SCAL] (1) PbSO4_.port_a.u = 8.31446261815324 * PbSO4_.Temperature * $FUN_3 + (((-918400.0) + 103.2 * ((-298.15) + PbSO4_.Temperature)) - PbSO4_.Temperature * (((-359.55056179775283) + 103.2 * log(0.0033540164346805303 * PbSO4_.Temperature)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.PbSO4_.stateOfMatter.molarVolumePure(Chemical.Substances.LeadSulfate_solid(0.30326, 0.0, -811200.0, -918400.0, 1.0, 103.2, {"http://www.mhhe.com/physsci/chemistry/chang7/ssg/graphics/chang7/pdf/cng7pa08.pdf, http://www.update.uu.se/~jolkkonen/pdf/CRC_TD.pdf"}, 0.0, 0.0, false, 0.0, 0.0, 1000.0), PbSO4_.Temperature, PbSO4_.Pressure, PbSO4_.ElectricPotential, PbSO4_.MoleFractionBasedIonicStrength) / PbSO4_.Temperature) * ((-1e5) + PbSO4_.Pressure))) ($RES_SIM_154) (34) [ARRY] (2) electrodeReaction1.rr * electrodeReaction1.s = electrodeReaction1.substrates.q ($RES_SIM_119) (35) [SCAL] (1) H2O.x = 1/(1.0 + H2O.SelfClustering_K) ($RES_SIM_37) (36) [SCAL] (1) PbSO4_.a = 1.0 ($RES_SIM_155) (37) [SCAL] (1) PbO2.port_a.u = 8.31446261815324 * PbO2.Temperature * $FUN_11 + (((-276600.0) + 64.6 * ((-298.15) + PbO2.Temperature)) - PbO2.Temperature * (((-193.19134663759854) + 64.6 * log(0.0033540164346805303 * PbO2.Temperature)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.PbO2.stateOfMatter.molarVolumePure(Chemical.Substances.LeadDioxide_solid(0.2391988, 0.0, -219000.0, -276600.0, 1.0, 64.6, {"http://www.vias.org/genchem/standard_enthalpies_table.html, http://www.update.uu.se/~jolkkonen/pdf/CRC_TD.pdf"}, 0.0, 0.0, false, 0.0, 0.0, 1000.0), PbO2.Temperature, PbO2.Pressure, PbO2.ElectricPotential, PbO2.MoleFractionBasedIonicStrength) / PbO2.Temperature) * ((-1e5) + PbO2.Pressure))) ($RES_SIM_39) (38) [SCAL] (1) $FUN_1 = log(Pb.a) ($RES_$AUX_391) (39) [SCAL] (1) $FUN_2 = log(HSO4.a) ($RES_$AUX_390) (40) [SCAL] (1) electrone.port_a.u = electrodeReaction1.products[3].u ($RES_SIM_311) (41) [SCAL] (1) electrone1.port_a.u = electrodeReaction.substrates[4].u ($RES_SIM_312) (42) [SCAL] (1) electrone.port_a.q + electrodeReaction1.products[3].q = 0.0 ($RES_SIM_275) (43) [SCAL] (1) H.port_a.u = electrodeReaction.substrates[3].u ($RES_SIM_313) (44) [SCAL] (1) H.port_a.u = electrodeReaction1.products[2].u ($RES_SIM_314) (45) [SCAL] (1) PbO2.port_a.u = electrodeReaction.substrates[1].u ($RES_SIM_315) (46) [SCAL] (1) electrodeReaction.products[2].u = H2O.port_a.u ($RES_SIM_316) (47) [SCAL] (1) electrodeReaction.products[1].u = PbSO4.port_a.u ($RES_SIM_317) (48) [SCAL] (1) PbSO4_.port_a.u = electrodeReaction1.products[1].u ($RES_SIM_318) (49) [SCAL] (1) HSO4.port_a.u = electrodeReaction1.substrates[2].u ($RES_SIM_319) (50) [SCAL] (1) electrodeReaction1.rr = -electrodeReaction1.du * electrodeReaction1.KC * $FUN_8 ($RES_SIM_120) (51) [SCAL] (1) electrodeReaction1.du = electrodeReaction1.p * electrodeReaction1.products.u - electrodeReaction1.s * electrodeReaction1.substrates.u ($RES_SIM_121) (52) [SCAL] (1) PbO2.a = 1.0 ($RES_SIM_40) (53) [-IF-] (1)if $SEV_2 then (53) [----] [SCAL] (1) (electrodeReaction.h_mix * electrodeReaction.products.q) / fill(1.0, 2) + electrodeReaction.substrates.q * {(-276600.0) + 64.6 * ((-298.15) + PbO2.Temperature), ($OMC$PositiveMax(-electrodeReaction1.substrates[2].q, 1e-7) * electrodeReaction1.substrates[2].h_outflow + $OMC$PositiveMax(-HSO4.port_a.q, 1e-7) * ((-885750.0) - 96485.33212331001 * HSO4.ElectricPotential)) / ($OMC$PositiveMax(-electrodeReaction1.substrates[2].q, 1e-7) + $OMC$PositiveMax(-HSO4.port_a.q, 1e-7)), ($OMC$PositiveMax(-electrodeReaction1.products[2].q, 1e-7) * electrodeReaction1.products[2].h_outflow + $OMC$PositiveMax(-H.port_a.q, 1e-7) * (96485.33212331001 * H.ElectricPotential)) / ($OMC$PositiveMax(-electrodeReaction1.products[2].q, 1e-7) + $OMC$PositiveMax(-H.port_a.q, 1e-7)), electrone1.port_a.h_outflow} = 0.0 ($RES_SIM_124) (53) [----] elseif $SEV_3 then (53) [----] [SCAL] (1) (electrodeReaction.h_mix * electrodeReaction.substrates.q) / fill(1.0, 4) + electrodeReaction.products.q * {(-918400.0) + 103.2 * ((-298.15) + PbSO4.Temperature), (-285830.0) + 75.3 * ((-298.15) + H2O.Temperature)} = 0.0 ($RES_SIM_125) (53) [----] else (53) [----] [SCAL] (1) electrodeReaction.h_mix = 0.0 ($RES_SIM_126) (53) [----] end if; (54) [SCAL] (1) $FUN_3 = log(PbSO4_.a) ($RES_$AUX_389) (55) [ARRY] (2) electrodeReaction.products.h_outflow = electrodeReaction.h_mix * fill(1.0, 2) ($RES_SIM_127) (56) [SCAL] (1) $FUN_4 = log(H.a) ($RES_$AUX_388) (57) [SCAL] (1) $SEV_0 = electrodeReaction1.rr > 0.0 ($RES_EVT_392) (58) [SCAL] (1) 96485.33212331001 * electrone1.substanceData.z * electrone1.port_a.q = 0.0 ($RES_SIM_81) (59) [ARRY] (4) electrodeReaction.substrates.h_outflow = electrodeReaction.h_mix * fill(1.0, 4) ($RES_SIM_128) (60) [SCAL] (1) $FUN_5 = abs(electrodeReaction.du) ($RES_$AUX_387) (61) [SCAL] (1) $SEV_1 = electrodeReaction1.rr < 0.0 ($RES_EVT_393) (62) [ARRY] (2) electrodeReaction.rr * electrodeReaction.p = -electrodeReaction.products.q ($RES_SIM_129) (63) [SCAL] (1) $FUN_6 = exp(-electrodeReaction.kE * $FUN_5) ($RES_$AUX_386) (64) [SCAL] (1) $SEV_2 = electrodeReaction.rr > 0.0 ($RES_EVT_394) (65) [SCAL] (1) $FUN_7 = abs(electrodeReaction1.du) ($RES_$AUX_385) (66) [SCAL] (1) $SEV_3 = electrodeReaction.rr < 0.0 ($RES_EVT_395) (67) [SCAL] (1) electrone.port_a.u = 8.31446261815324 * electrone.temperature * $FUN_9 + ((electrone.substanceData.DfH + ((-298.15) + electrone.temperature) * electrone.substanceData.Cp) - electrone.temperature * ((electrone.substanceData.Cp * log(0.0033540164346805303 * electrone.temperature) + 0.0033540164346805303 * (electrone.substanceData.DfH - electrone.substanceData.DfG)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone.stateOfMatter.molarVolumePure(electrone.substanceData, electrone.temperature, electrone.pressure, electrone.solution.v, 0.0) / electrone.temperature) * ((-1e5) + electrone.pressure))) + 96485.33212331001 * electrone.substanceData.z * electrone.solution.v ($RES_SIM_84) (68) [SCAL] (1) $FUN_8 = exp(-electrodeReaction1.kE * $FUN_7) ($RES_$AUX_384) (69) [SCAL] (1) electrone.a = electrone.substanceData.gamma ($RES_SIM_85) (70) [SCAL] (1) $FUN_9 = log(electrone.a) ($RES_$AUX_383) (71) [SCAL] (1) $FUN_10 = log(electrone1.a) ($RES_$AUX_382) (72) [SCAL] (1) $DER.electrone.temperature = 0.0 ($RES_SIM_242) (73) [SCAL] (1) $FUN_11 = log(PbO2.a) ($RES_$AUX_381) (74) [SCAL] (1) electrone.molarVolume = electrone.substanceData.MolarWeight / electrone.substanceData.density + Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone.stateOfMatter.molarVolumePure(electrone.substanceData, electrone.temperature, electrone.pressure, electrone.solution.v, 0.0) * log(Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone.stateOfMatter.activityCoefficient(electrone.substanceData, electrone.temperature, electrone.pressure, electrone.solution.v, 0.0)) ($RES_SIM_88) (75) [SCAL] (1) $FUN_12 = log(H2O.a) ($RES_$AUX_380) (76) [SCAL] (1) electrone.uPure = (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone.stateOfMatter.molarEnthalpyElectroneutral(electrone.substanceData, electrone.temperature, electrone.pressure, electrone.solution.v, 0.0) - electrone.temperature * Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.electrone.stateOfMatter.molarEntropyPure(electrone.substanceData, electrone.temperature, electrone.pressure, electrone.solution.v, 0.0)) + 96485.33212331001 * electrone.substanceData.z * electrone.solution.v ($RES_SIM_89) (77) [SCAL] (1) electrone.pressure = anode.BasePressure ($RES_SIM_245) (78) [SCAL] (1) HSO4.port_a.u = electrodeReaction.substrates[2].u ($RES_SIM_320) (79) [SCAL] (1) Pb.port_a.u = electrodeReaction1.substrates[1].u ($RES_SIM_321) (80) [SCAL] (1) electrodeReaction.substrates[2].q + electrodeReaction1.substrates[2].q + HSO4.port_a.q = 0.0 ($RES_SIM_285) (81) [SCAL] (1) H2O.port_a.q + electrodeReaction.products[2].q = 0.0 ($RES_SIM_286) (82) [SCAL] (1) PbSO4.a = 1.0 ($RES_SIM_6) (83) [SCAL] (1) PbSO4.port_a.q + electrodeReaction.products[1].q = 0.0 ($RES_SIM_287) (84) [SCAL] (1) PbSO4.port_a.u = 8.31446261815324 * PbSO4.Temperature * $FUN_13 + (((-918400.0) + 103.2 * ((-298.15) + PbSO4.Temperature)) - PbSO4.Temperature * (((-359.55056179775283) + 103.2 * log(0.0033540164346805303 * PbSO4.Temperature)) - (Chemical.Examples.CheckSubstancesData.StandardLeadAcidPotential.PbSO4.stateOfMatter.molarVolumePure(Chemical.Substances.LeadSulfate_solid(0.30326, 0.0, -811200.0, -918400.0, 1.0, 103.2, {"http://www.mhhe.com/physsci/chemistry/chang7/ssg/graphics/chang7/pdf/cng7pa08.pdf, http://www.update.uu.se/~jolkkonen/pdf/CRC_TD.pdf"}, 0.0, 0.0, false, 0.0, 0.0, 1000.0), PbSO4.Temperature, PbSO4.Pressure, PbSO4.ElectricPotential, PbSO4.MoleFractionBasedIonicStrength) / PbSO4.Temperature) * ((-1e5) + PbSO4.Pressure))) ($RES_SIM_5) (85) [SCAL] (1) electrodeReaction.substrates[3].q + electrodeReaction1.products[2].q + H.port_a.q = 0.0 ($RES_SIM_288) (86) [SCAL] (1) electrodeReaction1.products[1].q + PbSO4_.port_a.q = 0.0 ($RES_SIM_289) (87) [ARRY] (4) electrodeReaction.rr * electrodeReaction.s = electrodeReaction.substrates.q ($RES_SIM_130) (88) [SCAL] (1) electrodeReaction.rr = -electrodeReaction.du * electrodeReaction.KC * $FUN_6 ($RES_SIM_131)