** ecker-main-input_1C_charge.inp *heading new element qec3d8 - coupled temp disp electrical electrochemistry DOFS active are : 9, 11, 32, 33 *preprint, model=yes, history=yes, echo=yes *physical constants, STEFAN BOLTZMANN=1.38064852e-23, absolute zero = 0., faraday constant = 96485.3329, universal gas constant = 8.314, boltzmann=1.38064852e-23, avogadro number= 6.02214086e+23, elementary charge = 1.602176634e-19 *parameter totalanodeWidth = 85.e-3 totalcathodeWidth = 85.e-3 totalanodeLength = 101.e-3 totalcathodeLength = 101.e-3 totalanodeThickness = 73.7e-6 totalcathodeThickness = 54.5e-6 totalseparatorThickness = 19.0e-6 anodeWidth = totalanodeWidth/2.0 anodeLength = totalanodeLength/2.0 areaJellyRollSingle = totalanodeLength * totalanodeWidth TotalAreaJellyRoll = 48.0 * areaJellyRollSingle capacityNominal = 7.5 * 3600. spCapacityPerUnitAreaJellyRoll = capacityNominal / TotalAreaJellyRoll chargeRate = 1.0 time = 3600./chargeRate Simtime = 3600./chargeRate + 200. Dtmax = 10.0 initDt = 5.0 deCurrCathode = spCapacityPerUnitAreaJellyRoll / time numNodesPerFace = 4 deCurrCathode = deCurrCathode / (numNodesPerFace * 1.0) * areaJellyRollSingle cutoffVoltage = 4.2 ** ** ANODE ** *node, nset= nalla 1 , -0.0425 , -0.0505 , 0.000000E+00 2 , 0.0425 , -0.0505 , 0.000000E+00 3 , 0.0425 , 0.0505 , 0.000000E+00 4 , -0.0425 , 0.0505 , 0.000000E+00 5 , -0.0425 , -0.0505 , 2.948000E-06 6 , 0.0425 , -0.0505 , 2.948000E-06 7 , 0.0425 , 0.0505 , 2.948000E-06 8 , -0.0425 , 0.0505 , 2.948000E-06 9 , -0.0425 , -0.0505 , 5.896000E-06 10 , 0.0425 , -0.0505 , 5.896000E-06 11 , 0.0425 , 0.0505 , 5.896000E-06 12 , -0.0425 , 0.0505 , 5.896000E-06 13 , -0.0425 , -0.0505 , 8.844000E-06 14 , 0.0425 , -0.0505 , 8.844000E-06 15 , 0.0425 , 0.0505 , 8.844000E-06 16 , -0.0425 , 0.0505 , 8.844000E-06 17 , -0.0425 , -0.0505 , 1.179200E-05 18 , 0.0425 , -0.0505 , 1.179200E-05 19 , 0.0425 , 0.0505 , 1.179200E-05 20 , -0.0425 , 0.0505 , 1.179200E-05 21 , -0.0425 , -0.0505 , 1.474000E-05 22 , 0.0425 , -0.0505 , 1.474000E-05 23 , 0.0425 , 0.0505 , 1.474000E-05 24 , -0.0425 , 0.0505 , 1.474000E-05 25 , -0.0425 , -0.0505 , 1.768800E-05 26 , 0.0425 , -0.0505 , 1.768800E-05 27 , 0.0425 , 0.0505 , 1.768800E-05 28 , -0.0425 , 0.0505 , 1.768800E-05 29 , -0.0425 , -0.0505 , 2.063600E-05 30 , 0.0425 , -0.0505 , 2.063600E-05 31 , 0.0425 , 0.0505 , 2.063600E-05 32 , -0.0425 , 0.0505 , 2.063600E-05 33 , -0.0425 , -0.0505 , 2.358400E-05 34 , 0.0425 , -0.0505 , 2.358400E-05 35 , 0.0425 , 0.0505 , 2.358400E-05 36 , -0.0425 , 0.0505 , 2.358400E-05 37 , -0.0425 , -0.0505 , 2.653200E-05 38 , 0.0425 , -0.0505 , 2.653200E-05 39 , 0.0425 , 0.0505 , 2.653200E-05 40 , -0.0425 , 0.0505 , 2.653200E-05 41 , -0.0425 , -0.0505 , 2.948000E-05 42 , 0.0425 , -0.0505 , 2.948000E-05 43 , 0.0425 , 0.0505 , 2.948000E-05 44 , -0.0425 , 0.0505 , 2.948000E-05 45 , -0.0425 , -0.0505 , 3.242800E-05 46 , 0.0425 , -0.0505 , 3.242800E-05 47 , 0.0425 , 0.0505 , 3.242800E-05 48 , -0.0425 , 0.0505 , 3.242800E-05 49 , -0.0425 , -0.0505 , 3.537600E-05 50 , 0.0425 , -0.0505 , 3.537600E-05 51 , 0.0425 , 0.0505 , 3.537600E-05 52 , -0.0425 , 0.0505 , 3.537600E-05 53 , -0.0425 , -0.0505 , 3.832400E-05 54 , 0.0425 , -0.0505 , 3.832400E-05 55 , 0.0425 , 0.0505 , 3.832400E-05 56 , -0.0425 , 0.0505 , 3.832400E-05 57 , -0.0425 , -0.0505 , 4.127200E-05 58 , 0.0425 , -0.0505 , 4.127200E-05 59 , 0.0425 , 0.0505 , 4.127200E-05 60 , -0.0425 , 0.0505 , 4.127200E-05 61 , -0.0425 , -0.0505 , 4.422000E-05 62 , 0.0425 , -0.0505 , 4.422000E-05 63 , 0.0425 , 0.0505 , 4.422000E-05 64 , -0.0425 , 0.0505 , 4.422000E-05 65 , -0.0425 , -0.0505 , 4.716800E-05 66 , 0.0425 , -0.0505 , 4.716800E-05 67 , 0.0425 , 0.0505 , 4.716800E-05 68 , -0.0425 , 0.0505 , 4.716800E-05 69 , -0.0425 , -0.0505 , 5.011600E-05 70 , 0.0425 , -0.0505 , 5.011600E-05 71 , 0.0425 , 0.0505 , 5.011600E-05 72 , -0.0425 , 0.0505 , 5.011600E-05 73 , -0.0425 , -0.0505 , 5.306400E-05 74 , 0.0425 , -0.0505 , 5.306400E-05 75 , 0.0425 , 0.0505 , 5.306400E-05 76 , -0.0425 , 0.0505 , 5.306400E-05 77 , -0.0425 , -0.0505 , 5.601200E-05 78 , 0.0425 , -0.0505 , 5.601200E-05 79 , 0.0425 , 0.0505 , 5.601200E-05 80 , -0.0425 , 0.0505 , 5.601200E-05 81 , -0.0425 , -0.0505 , 5.896000E-05 82 , 0.0425 , -0.0505 , 5.896000E-05 83 , 0.0425 , 0.0505 , 5.896000E-05 84 , -0.0425 , 0.0505 , 5.896000E-05 85 , -0.0425 , -0.0505 , 6.190800E-05 86 , 0.0425 , -0.0505 , 6.190800E-05 87 , 0.0425 , 0.0505 , 6.190800E-05 88 , -0.0425 , 0.0505 , 6.190800E-05 89 , -0.0425 , -0.0505 , 6.485600E-05 90 , 0.0425 , -0.0505 , 6.485600E-05 91 , 0.0425 , 0.0505 , 6.485600E-05 92 , -0.0425 , 0.0505 , 6.485600E-05 93 , -0.0425 , -0.0505 , 6.780400E-05 94 , 0.0425 , -0.0505 , 6.780400E-05 95 , 0.0425 , 0.0505 , 6.780400E-05 96 , -0.0425 , 0.0505 , 6.780400E-05 97 , -0.0425 , -0.0505 , 7.075200E-05 98 , 0.0425 , -0.0505 , 7.075200E-05 99 , 0.0425 , 0.0505 , 7.075200E-05 100 , -0.0425 , 0.0505 , 7.075200E-05 101 , -0.0425 , -0.0505 , 7.370000E-05 102 , 0.0425 , -0.0505 , 7.370000E-05 103 , 0.0425 , 0.0505 , 7.370000E-05 104 , -0.0425 , 0.0505 , 7.370000E-05 *nset, nset=nAnodeStart 1, 2, 3, 4 *nset, nset=nAnodeEnd, generate 101, 104, 1 *element, type = qec3d8, elset= ealla 1, 1, 2, 3, 4, 5, 6, 7, 8 2, 5, 6, 7, 8, 9, 10, 11, 12 3, 9, 10, 11, 12, 13, 14, 15, 16 4, 13, 14, 15, 16, 17, 18, 19, 20 5, 17, 18, 19, 20, 21, 22, 23, 24 6, 21, 22, 23, 24, 25, 26, 27, 28 7, 25, 26, 27, 28, 29, 30, 31, 32 8, 29, 30, 31, 32, 33, 34, 35, 36 9, 33, 34, 35, 36, 37, 38, 39, 40 10, 37, 38, 39, 40, 41, 42, 43, 44 11, 41, 42, 43, 44, 45, 46, 47, 48 12, 45, 46, 47, 48, 49, 50, 51, 52 13, 49, 50, 51, 52, 53, 54, 55, 56 14, 53, 54, 55, 56, 57, 58, 59, 60 15, 57, 58, 59, 60, 61, 62, 63, 64 16, 61, 62, 63, 64, 65, 66, 67, 68 17, 65, 66, 67, 68, 69, 70, 71, 72 18, 69, 70, 71, 72, 73, 74, 75, 76 19, 73, 74, 75, 76, 77, 78, 79, 80 20, 77, 78, 79, 80, 81, 82, 83, 84 21, 81, 82, 83, 84, 85, 86, 87, 88 22, 85, 86, 87, 88, 89, 90, 91, 92 23, 89, 90, 91, 92, 93, 94, 95, 96 24, 93, 94, 95, 96, 97, 98, 99, 100 25, 97, 98, 99, 100, 101, 102, 103, 104 *elset, elset=eoutput 1, 5, 10, 15, 20 , 25 *solid section, elset=ealla, material=ABQ_EChemPET_matAnode 1 *initial conditions, type=temperature nalla, 298.1500 *initial conditions, type=solid electric potential nalla, 0. *initial conditions, type=fluid electric potential nalla, -0.6126932965 *initial conditions, type=ion concentration nalla, 1000. ** ** SEPARATOR ** *node 105 , -0.0425 , -0.0505 , 7.542727E-05 106 , 0.0425 , -0.0505 , 7.542727E-05 107 , 0.0425 , 0.0505 , 7.542727E-05 108 , -0.0425 , 0.0505 , 7.542727E-05 109 , -0.0425 , -0.0505 , 7.715455E-05 110 , 0.0425 , -0.0505 , 7.715455E-05 111 , 0.0425 , 0.0505 , 7.715455E-05 112 , -0.0425 , 0.0505 , 7.715455E-05 113 , -0.0425 , -0.0505 , 7.888182E-05 114 , 0.0425 , -0.0505 , 7.888182E-05 115 , 0.0425 , 0.0505 , 7.888182E-05 116 , -0.0425 , 0.0505 , 7.888182E-05 117 , -0.0425 , -0.0505 , 8.060909E-05 118 , 0.0425 , -0.0505 , 8.060909E-05 119 , 0.0425 , 0.0505 , 8.060909E-05 120 , -0.0425 , 0.0505 , 8.060909E-05 121 , -0.0425 , -0.0505 , 8.233636E-05 122 , 0.0425 , -0.0505 , 8.233636E-05 123 , 0.0425 , 0.0505 , 8.233636E-05 124 , -0.0425 , 0.0505 , 8.233636E-05 125 , -0.0425 , -0.0505 , 8.406364E-05 126 , 0.0425 , -0.0505 , 8.406364E-05 127 , 0.0425 , 0.0505 , 8.406364E-05 128 , -0.0425 , 0.0505 , 8.406364E-05 129 , -0.0425 , -0.0505 , 8.579091E-05 130 , 0.0425 , -0.0505 , 8.579091E-05 131 , 0.0425 , 0.0505 , 8.579091E-05 132 , -0.0425 , 0.0505 , 8.579091E-05 133 , -0.0425 , -0.0505 , 8.751818E-05 134 , 0.0425 , -0.0505 , 8.751818E-05 135 , 0.0425 , 0.0505 , 8.751818E-05 136 , -0.0425 , 0.0505 , 8.751818E-05 137 , -0.0425 , -0.0505 , 8.924545E-05 138 , 0.0425 , -0.0505 , 8.924545E-05 139 , 0.0425 , 0.0505 , 8.924545E-05 140 , -0.0425 , 0.0505 , 8.924545E-05 141 , -0.0425 , -0.0505 , 9.097273E-05 142 , 0.0425 , -0.0505 , 9.097273E-05 143 , 0.0425 , 0.0505 , 9.097273E-05 144 , -0.0425 , 0.0505 , 9.097273E-05 145 , -0.0425 , -0.0505 , 9.270000E-05 146 , 0.0425 , -0.0505 , 9.270000E-05 147 , 0.0425 , 0.0505 , 9.270000E-05 148 , -0.0425 , 0.0505 , 9.270000E-05 *nset, generate, nset=nalls 105, 144, 1 *element, type = qec3d8, elset= ealls 101, 101, 102, 103, 104, 105, 106, 107, 108 102, 105, 106, 107, 108, 109, 110, 111, 112 103, 109, 110, 111, 112, 113, 114, 115, 116 104, 113, 114, 115, 116, 117, 118, 119, 120 105, 117, 118, 119, 120, 121, 122, 123, 124 106, 121, 122, 123, 124, 125, 126, 127, 128 107, 125, 126, 127, 128, 129, 130, 131, 132 108, 129, 130, 131, 132, 133, 134, 135, 136 109, 133, 134, 135, 136, 137, 138, 139, 140 110, 137, 138, 139, 140, 141, 142, 143, 144 111, 141, 142, 143, 144, 145, 146, 147, 148 *solid section, elset=ealls, material=ABQ_EChemPET_matSeparator 1 *initial conditions, type=ion concentration nalls, 1000. *initial conditions, type=fluid electric potential nalls, -0.6126932965 *initial conditions, type=temperature nalls, 298.15 ** ** CATHODE ** *node, nset=nallc 149 , -0.0425 , -0.0505 , 9.488000E-05 150 , 0.0425 , -0.0505 , 9.488000E-05 151 , 0.0425 , 0.0505 , 9.488000E-05 152 , -0.0425 , 0.0505 , 9.488000E-05 153 , -0.0425 , -0.0505 , 9.706000E-05 154 , 0.0425 , -0.0505 , 9.706000E-05 155 , 0.0425 , 0.0505 , 9.706000E-05 156 , -0.0425 , 0.0505 , 9.706000E-05 157 , -0.0425 , -0.0505 , 9.924000E-05 158 , 0.0425 , -0.0505 , 9.924000E-05 159 , 0.0425 , 0.0505 , 9.924000E-05 160 , -0.0425 , 0.0505 , 9.924000E-05 161 , -0.0425 , -0.0505 , 1.014200E-04 162 , 0.0425 , -0.0505 , 1.014200E-04 163 , 0.0425 , 0.0505 , 1.014200E-04 164 , -0.0425 , 0.0505 , 1.014200E-04 165 , -0.0425 , -0.0505 , 1.036000E-04 166 , 0.0425 , -0.0505 , 1.036000E-04 167 , 0.0425 , 0.0505 , 1.036000E-04 168 , -0.0425 , 0.0505 , 1.036000E-04 169 , -0.0425 , -0.0505 , 1.057800E-04 170 , 0.0425 , -0.0505 , 1.057800E-04 171 , 0.0425 , 0.0505 , 1.057800E-04 172 , -0.0425 , 0.0505 , 1.057800E-04 173 , -0.0425 , -0.0505 , 1.079600E-04 174 , 0.0425 , -0.0505 , 1.079600E-04 175 , 0.0425 , 0.0505 , 1.079600E-04 176 , -0.0425 , 0.0505 , 1.079600E-04 177 , -0.0425 , -0.0505 , 1.101400E-04 178 , 0.0425 , -0.0505 , 1.101400E-04 179 , 0.0425 , 0.0505 , 1.101400E-04 180 , -0.0425 , 0.0505 , 1.101400E-04 181 , -0.0425 , -0.0505 , 1.123200E-04 182 , 0.0425 , -0.0505 , 1.123200E-04 183 , 0.0425 , 0.0505 , 1.123200E-04 184 , -0.0425 , 0.0505 , 1.123200E-04 185 , -0.0425 , -0.0505 , 1.145000E-04 186 , 0.0425 , -0.0505 , 1.145000E-04 187 , 0.0425 , 0.0505 , 1.145000E-04 188 , -0.0425 , 0.0505 , 1.145000E-04 189 , -0.0425 , -0.0505 , 1.166800E-04 190 , 0.0425 , -0.0505 , 1.166800E-04 191 , 0.0425 , 0.0505 , 1.166800E-04 192 , -0.0425 , 0.0505 , 1.166800E-04 193 , -0.0425 , -0.0505 , 1.188600E-04 194 , 0.0425 , -0.0505 , 1.188600E-04 195 , 0.0425 , 0.0505 , 1.188600E-04 196 , -0.0425 , 0.0505 , 1.188600E-04 197 , -0.0425 , -0.0505 , 1.210400E-04 198 , 0.0425 , -0.0505 , 1.210400E-04 199 , 0.0425 , 0.0505 , 1.210400E-04 200 , -0.0425 , 0.0505 , 1.210400E-04 201 , -0.0425 , -0.0505 , 1.232200E-04 202 , 0.0425 , -0.0505 , 1.232200E-04 203 , 0.0425 , 0.0505 , 1.232200E-04 204 , -0.0425 , 0.0505 , 1.232200E-04 205 , -0.0425 , -0.0505 , 1.254000E-04 206 , 0.0425 , -0.0505 , 1.254000E-04 207 , 0.0425 , 0.0505 , 1.254000E-04 208 , -0.0425 , 0.0505 , 1.254000E-04 209 , -0.0425 , -0.0505 , 1.275800E-04 210 , 0.0425 , -0.0505 , 1.275800E-04 211 , 0.0425 , 0.0505 , 1.275800E-04 212 , -0.0425 , 0.0505 , 1.275800E-04 213 , -0.0425 , -0.0505 , 1.297600E-04 214 , 0.0425 , -0.0505 , 1.297600E-04 215 , 0.0425 , 0.0505 , 1.297600E-04 216 , -0.0425 , 0.0505 , 1.297600E-04 217 , -0.0425 , -0.0505 , 1.319400E-04 218 , 0.0425 , -0.0505 , 1.319400E-04 219 , 0.0425 , 0.0505 , 1.319400E-04 220 , -0.0425 , 0.0505 , 1.319400E-04 221 , -0.0425 , -0.0505 , 1.341200E-04 222 , 0.0425 , -0.0505 , 1.341200E-04 223 , 0.0425 , 0.0505 , 1.341200E-04 224 , -0.0425 , 0.0505 , 1.341200E-04 225 , -0.0425 , -0.0505 , 1.363000E-04 226 , 0.0425 , -0.0505 , 1.363000E-04 227 , 0.0425 , 0.0505 , 1.363000E-04 228 , -0.0425 , 0.0505 , 1.363000E-04 229 , -0.0425 , -0.0505 , 1.384800E-04 230 , 0.0425 , -0.0505 , 1.384800E-04 231 , 0.0425 , 0.0505 , 1.384800E-04 232 , -0.0425 , 0.0505 , 1.384800E-04 233 , -0.0425 , -0.0505 , 1.406600E-04 234 , 0.0425 , -0.0505 , 1.406600E-04 235 , 0.0425 , 0.0505 , 1.406600E-04 236 , -0.0425 , 0.0505 , 1.406600E-04 237 , -0.0425 , -0.0505 , 1.428400E-04 238 , 0.0425 , -0.0505 , 1.428400E-04 239 , 0.0425 , 0.0505 , 1.428400E-04 240 , -0.0425 , 0.0505 , 1.428400E-04 241 , -0.0425 , -0.0505 , 1.450200E-04 242 , 0.0425 , -0.0505 , 1.450200E-04 243 , 0.0425 , 0.0505 , 1.450200E-04 244 , -0.0425 , 0.0505 , 1.450200E-04 245 , -0.0425 , -0.0505 , 1.472000E-04 246 , 0.0425 , -0.0505 , 1.472000E-04 247 , 0.0425 , 0.0505 , 1.472000E-04 248 , -0.0425 , 0.0505 , 1.472000E-04 *nset, nset=nallc, generate 145, 248, 1 *nset, nset=nCathodeStart 145, 146, 147, 148 *nset, nset=nCathodeEnd 245, 246, 247, 248 *element, type = qec3d8, elset= eallc 201, 145, 146, 147, 148, 149, 150, 151, 152 202, 149, 150, 151, 152, 153, 154, 155, 156 203, 153, 154, 155, 156, 157, 158, 159, 160 204, 157, 158, 159, 160, 161, 162, 163, 164 205, 161, 162, 163, 164, 165, 166, 167, 168 206, 165, 166, 167, 168, 169, 170, 171, 172 207, 169, 170, 171, 172, 173, 174, 175, 176 208, 173, 174, 175, 176, 177, 178, 179, 180 209, 177, 178, 179, 180, 181, 182, 183, 184 210, 181, 182, 183, 184, 185, 186, 187, 188 211, 185, 186, 187, 188, 189, 190, 191, 192 212, 189, 190, 191, 192, 193, 194, 195, 196 213, 193, 194, 195, 196, 197, 198, 199, 200 214, 197, 198, 199, 200, 201, 202, 203, 204 215, 201, 202, 203, 204, 205, 206, 207, 208 216, 205, 206, 207, 208, 209, 210, 211, 212 217, 209, 210, 211, 212, 213, 214, 215, 216 218, 213, 214, 215, 216, 217, 218, 219, 220 219, 217, 218, 219, 220, 221, 222, 223, 224 220, 221, 222, 223, 224, 225, 226, 227, 228 221, 225, 226, 227, 228, 229, 230, 231, 232 222, 229, 230, 231, 232, 233, 234, 235, 236 223, 233, 234, 235, 236, 237, 238, 239, 240 224, 237, 238, 239, 240, 241, 242, 243, 244 225, 241, 242, 243, 244, 245, 246, 247, 248 *elset, elset=eoutput 201, 205, 210, 215, 220 , 225 *elset, elset=eSensorCsurf 1, *nset, nset=nSensorPot 248 *solid section, elset=eallc, material=ABQ_EChemPET_matCathode 1 *initial conditions, type=temperature nallc, 298.1500 *initial conditions, type=solid electric potential nallc, 3.00784265 *initial conditions, type=fluid electric potential nallc, -0.6126932965 *initial conditions, type=ion concentration nallc, 1000. *elset, elset = eall ealla, ealls, eallc *include, input = ABQ_EchemPET_Types.inp *include, input = ecker-anode_propParam-charge.inp *include, input = ecker-cathode_propParam-charge.inp *include, input = ecker-separator_propParam-charge.inp ** ************************** *Amplitude, name=CUTOFFDETECTION, definition=USER, variables=0, properties=1 ************************** *step , name=charge, inc = 250000000, unsymm=yes *Coupled thermal-ElectroChemical , , , *amplitude, name=amp1 0, 1, 1,1 *boundary, op = new *boundary, op = new, amplitude=CUTOFFDETECTION nAnodeStart, 9, 9, 0.0 *amplitude, name=rampcurr 0,-1.,10000,-1. *CECURRENT, op = new, amp=rampcurr nCathodeEnd, 9, - *Film, amplitude=amp1 EALL, ffs, 298.15, 0.10 *output,field,freq=5 *node output epot, epote, nnce, nt, recur, rfl,rfle recure, rflce *element output hfl, mfle, ecde, conce, elecpote, elecpot, TEMP csurfavg, *output,history *element output, elset =eoutput csurfavg, elecpot, elecpote, conce, temp *output, history, freq=1, sensor, name=VOLTAGECUTOFF *node output, nset =nSensorPot ePot, *monitor, node=245, dof=9 *end step