30:T5d0,md`Property | Value ---|--- COP | ${res.COP.toFixed(3)} COP_reference | ${res.COP_ref.toFixed(3)} Improvement with IHX | ${improvementWithIHX.toFixed(2)} % Evaporator Capacity | ${mat.format(res.Q, { notation: "fixed", precision: 2 })} Evaporator Capacity ref | ${mat.format(res.Q_ref, { notation: "fixed", precision: 2 })} SubCooling with IHX | ${mat.format( mat.subtract(mat.to(res.state[3].T, "K"), mat.to(res.state[4].T, "K")), { notation: "fixed", precision: 1 } )} SuperHeating with IHX | ${mat.format( mat.subtract(mat.to(res.state[1].T, "K"), mat.to(res.state[0].T, "K")), { notation: "fixed", precision: 1 } )} Evaproation pressure | ${mat.format( mat.to(mat.unit(cycle[0].P, "Pa"), "psi"), { notation: "fixed", precision: 1 } )} Hot Side pressure | ${mat.format(mat.to(mat.unit(cycle[3].P, "Pa"), "psi"), { notation: "fixed", precision: 1 })} IHX Q | ${mat.format(mat.to(res.Q_IHX, "kW"), { notation: "fixed", precision: 2 })} IHX LMTD | ${mat.format(mat.to(res.IHX.LMTD, "K"), { notation: "fixed", precision: 2 })} IHX area | ${mat.format(res.IHX.A.to("m^2"), { notation: "fixed", precision: 2 })} IHX cube side | ${mat.format( mat.to(mat.pow(res.IHX.Volume, 1 / 3), "mm"), { notation: "fixed", precision: 2 } )} IHX K x A | ${mat.format( mat.dotMultiply(res.IHX.U, res.IHX.A).to("W/K"), { notation: "fixed", precision: 2 } )} `31:T620,function nelderMead1D(f, x0, opts = {}) { const step = opts.step ?? 1e-1; const tol = opts.tol ?? 1e-8; const maxIter = opts.maxIter ?? 200; const abs = Math.abs let x = [x0, x0 + (step || 1e-3)]; let fx = [evalSafe(f, x[0]), evalSafe(f, x[1])]; const order = () => { if (fx[0] <= fx[1]) return; [x[0], x[1]] = [x[1], x[0]]; [fx[0], fx[1]] = [fx[1], fx[0]]; }; order(); let iter = 0; while (iter++ < maxIter) { if (abs(x[1] - x[0]) <= tol || abs(fx[1] - fx[0]) <= tol) break; const c = x[0]; // centroid is the best point in 1D const xr = c + (c - x[1]); // reflection (alpha=1) const fr = evalSafe(f, xr); if (fr < fx[0]) { // expansion (gamma=2) const xe = c + 2 * (xr - c); const fe = evalSafe(f, xe); if (fe < fr) { x[1] = xe; fx[1] = fe; } else { x[1] = xr; fx[1] = fr; } } else if (fr < fx[1]) { // accept reflection x[1] = xr; fx[1] = fr; } else { // contraction (rho=0.5) or shrink (sigma=0.5) const xc = c + 0.5 * (x[1] - c); const fc = evalSafe(f, xc); if (fc < fx[1]) { x[1] = xc; fx[1] = fc; } else { x[1] = c + 0.5 * (x[1] - c); fx[1] = evalSafe(f, x[1]); } } order(); } return { x: x[0], fx: fx[0], iterations: iter - 1, success: iter <= maxIter }; function evalSafe(fn, xi) { try { const v = fn(xi); return Number.isFinite(v) && !Number.isNaN(v) ? v : Infinity; } catch (e) { return Infinity; } } }32:T5b9,cycleCalculation = mat.compile(` c = [{},{},{},{},{},{}]; # cycle with IHX # 6 to 1 Evaporation c[1].P = pLow; c[1].T = evap.T + evap.superHeating; c[1].D = prop('D', {'T|gas':c[1].T, P:c[1].P}); c[1].H = prop('H', {'T|gas':c[1].T, P:c[1].P}); c[1].S = prop('S', {'T|gas':c[1].T, P:c[1].P}); # 1 to 2 IHX Low c[2].P = c[1].P; # fix this c[4].T = cond.T - cond.subCooling; H_eta = prop('H', {T: c[4].T, P:c[2].P}); c[2].H = IHX.epsilon * (H_eta - c[1].H) + c[1].H; c[2].T = prop('T', c[2]); c[2].D = prop('D', c[2]); c[2].S = prop('S', c[2]); # 2 to 3 Compression of vapor c[3].P = pHigh + cond.P_drop; H_i = prop('H', { 'P': c[3].P, 'S': c[2].S }); c[3].H = (H_i - c[2].H) / etaS + c[2].H; c[3].T = prop('T', c[3]); c[3].D = prop('D', c[3]); c[3].S = prop('S', c[3]); # 3 to 4 Condensation or Gas Cooling c[4].P = pHigh; c4 = {P:c[4].P}; c4[isTransCritical ? "T": "T|liquid"] = c[4].T c[4].D = prop('D', c4); c[4].H = prop('H', c4); c[4].S = prop('S', c4); # 4 to 5 IHX high c[5].H = c[1].H - c[2].H + c[4].H; c[5].P = c[4].P; c[5].T = prop('T', c[5]); c[5].D = prop('D', c[5]); c[5].S = prop('S', c[5]); # 5 to 6 Expansion c[6].H = c[5].H; c[6].P = c[1].P + evap.P_drop; c[6].T = prop('T', c[6]); c[6].D = prop('D', c[6]); c[6].S = prop('S', c[6]); # Work, Energy and Performance W_comp = mDot*(c[3].H - c[2].H); # Q_h = mDot*(c[4].H - c[5].H); Q_c = mDot*(c[1].H - c[6].H); Q_IHX = mDot*(c[2].H-c[1].H); evap_COP = Q_c/W_comp; # cond_COP = Q_h/W_comp; `)33:T81c,IHX_calc = mat.compile(` # # IHX information Q_IHX = mDot * (c[2].H - c[1].H) deltaA = (c[4].T to K ) - (c[2].T to K); deltaB = (c[5].T to K) - (c[1].T to K); # LMTD IHX.LMTD = (deltaA - deltaB) / log(deltaA / deltaB) # Cell calcs IHX.cellVol = IHX.cellSize ^ 3; cellSizeToAreaFactor = 3.8424; # for schwarz D with a uniform cell size IHX.cellArea = cellSizeToAreaFactor * IHX.cellSize ^ 2; # Characteristic Length V_body / Surface Area # L = 0.260254 mm / 2; # characteristic length # Estimation of the volume of the hot or cold fluid # Half of the remaining space of the Volume of the fluids ( volume of the material as # the area of heat exchange * thickness) L = ((IHX.cellVol - IHX.cellArea* IHX.thickness) / 2) / IHX.cellArea # ## IHX of the cold side being heated # Reynolds rho = mean(c[1].D, c[2].D); # L = A_t*L / A_int transversal_fraction = (IHX.cellSize ^ 2 - 2*sqrt(2*IHX.cellSize^2)*IHX.thickness) / ( IHX.cellSize ^ 2 ) A = IHX.transversal_area / 2 * transversal_fraction; u = mDot / rho / A; # flow speed miu = mean(prop('V', sat_gas), prop('V', c[2])); Re = rho * u * L / miu # Prandtl number cold side k = mean(prop('L', sat_gas), prop('L', c[2])); Cp = mean(prop('C', sat_gas), prop('C', c[2])); Pr = Cp * miu / k # Nusselt number heating (cold side) Nu = 0.023 * Re ^ (4/5) * Pr ^ 0.3 # Convection heat transfer cold side h1 = Nu * k / L to W / (m^2 K) # ## IHX of the hot side being cooled # Reynolds rho = mean(c[4].D, c[5].D); u = mDot / rho / A; # flow speed miu = mean(prop('V', sat_liq), prop('V', c[5])); Re = rho * u * L / miu # Prandtl number cold side Cp = mean(prop('C', sat_liq), prop('C', c[5])); k = mean(prop('L', sat_liq), prop('L', c[5])); Pr = Cp * miu / k # Nusselt number heating (Dittus-Bolter equation) Nu = 0.023 * Re ^ (4/5) * Pr ^ 0.4 h2 = Nu * k / L to W / (m^2 K) IHX.U = 1/(1 / h1 + IHX.thickness/IHX.k + 1 / h2); # IHX.U = 1/(IHX.thickness/IHX.k); IHX.A = Q_IHX / (IHX.U * IHX.LMTD); # Recuprator's Volume IHX.Volume = IHX.A * IHX.cellVol / IHX.cellArea to mm ^ 3 # Side of a IHX Cube IHX.Volume ^ (1 / 3) `)34:T510,secondCycleCalculation = mat.compile(` prop(DesiredProperty, FluidState) = props(DesiredProperty, fluid, FluidState); # Without IHX c_ref = [{},{},{},{}]; # cycle without IHX c_ref[1].P = pLow; c_ref[1].T = evap.T + evap.superHeating; c_ref[1].D = prop('D', {'T|gas':c_ref[1].T, P:c_ref[1].P}); c_ref[1].H = prop('H', {'T|gas':c_ref[1].T, P:c_ref[1].P}); c_ref[1].S = prop('S', {'T|gas':c_ref[1].T, P:c_ref[1].P}); # 1 to 2 Compresison of vapor c_ref[2].P = pHigh_ref + cond.P_drop;; H_i = prop('H', { 'P': c_ref[2].P, 'S': c_ref[1].S }); c_ref[2].H = (H_i - c_ref[1].H) / etaS1 + c_ref[1].H; c_ref[2].T = prop('T', c_ref[2]); c_ref[2].D = prop('D', c_ref[2]); c_ref[2].S = prop('S', c_ref[2]); # 2 to 3 Condensation or Gas Cooling c_ref[3].T = cond.T-cond.subCooling; c_ref[3].P = pHigh_ref; c3 = {P:c_ref[3].P}; c3[isTransCritical ? "T": "T|liquid"] = c_ref[3].T c_ref[3].D = prop('D', c3); c_ref[3].H = prop('H', c3); c_ref[3].S = prop('S', c3); # 3 to 4 Expansion c_ref[4].H = c_ref[3].H; c_ref[4].P = c_ref[1].P + evap.P_drop; c_ref[4].T = prop('T', c_ref[4]); c_ref[4].D = prop('D', c_ref[4]); c_ref[4].S = prop('S', c_ref[4]); # Work, Energy and Performance W_comp_ref = mDot*(c_ref[2].H - c_ref[1].H); Q_c_ref = mDot*(c_ref[1].H - c_ref[4].H); evap_COP_ref = Q_c_ref/W_comp_ref; `)35:T650,fluids = [ "Acetone", "CO2", "Air", "Ammonia", "Argon", "Benzene", "iso-Butane", "n-Butane", "cis-2-Butene", "trans-2-Butene", "iso-Butene", "1-Butene", "Carbon dioxide", "Carbon monoxide", "Carbonyl sulfide", "n-Decane", "Deuterium", "ortho-Deuterium", "para-Deuterium", "Dimethyl carbonate", "Dimethyl ether", "n-Dodecane", "D4", "D5", "D6", "Ethane", "Ethanol", "Ethylene", "Ethyl benzene", "Fluorine", "Helium", "cyclo-Hexane", "iso-Hexane", "n-Hexane", "n-Heptane", "HFE143m", "normal-Hydrogen", "ortho-Hydrogen", "para-Hydrogen", "Hydrogen sulfide", "Krypton", "Methane", "Methanol", "Methyl linoleate", "Methyl linolenate", "Methyl oleate", "Methyl palmitate", "Methyl stearate", "MD2M", "MD3M", "MD4M", "MDM", "MM", "Neon", "Nitrogen", "Nitrous oxide", "n-Nonane", "n-Octane", "Oxygen", "cyclo-Pentane", "iso-Pentane", "n-Pentane", "neo-Pentane", "cyclo-Propane", "n-Propane", "Propylene", "Propyne", "R11", "R113", "R114", "R115", "R116", "R12", "R123", "R1233zd(E)", "R1234yf", "R1234ze(E)", "R1234ze(Z)", "R124", "R1243zf", "R125", "R13", "R134a", "R13I1", "R14", "R141b", "R142b", "R143a", "R152A", "R161", "R21", "R218", "R22", "R227EA", "R23", "R236EA", "R236FA", "R245ca", "R245fa", "R32", "R365MFC", "R40", "R404A", "R407C", "R41", "R410A", "R507A", "SulfurDioxide", "SulfurHexafluoride", "Toluene", "n-Undecane", "Water", "Xenon", "m-Xylene", "o-Xylene", "p-Xylene" ]36:[] 2e:{"title":"IHX (Internal Heat Exchanger)","cells":[{"id":0,"value":"# IHX (Internal Heat Exchanger)\n\nCalculation of a [vapor compression cycle](https://en.wikipedia.org/wiki/Vapor-compression_refrigeration) with an [internal heat exchanger](https://en.wikipedia.org/wiki/Subcooling)\n\nThis takes into account trans-critical examples. Thus the hot side represents either a condenser or gas cooler.\n\n## To do\n- Separate subcooling and superheating with and without IHX\n- Fix LMTD issue (I think it's a mathjs issue)\n- Include option to optimize subcooling (regardless if it's transcrítical)\n\n## Maybe\n\n- consider the compressor efficiency estimations according to AHRI Standard 540\n - [scroll compressors](https://doi.org/10.1016/j.applthermaleng.2022.119432)\n - reciprocating compressors: [Analysis of map-based models for reciprocating compressors and optimum\nselection of rating points](https://doi.org/10.1016/j.ijrefrig.2023.06.002)\n- filter refrigerants ","pinned":false,"hidden":false,"mode":"md"},{"id":2277,"value":"# IHX efficiency\n\nFormula\n\n${tex.block`\\eta_{IHX} = \\frac{H(at\\ T_{IHX\\ LP\\ out}\\ \\&\\ P_{IHX\\ LP\\ out}) - H(at\\ T_{IHX\\ LP\\ in}\\ \\&\\ P_{IHX\\ LP\\ in})}{H(at\\ T_{IHX\\ HP\\ in}\\ \\&\\ P_{IHX\\ LP\\ out}) - H(at\\ T_{IHX\\ LP\\ in}\\ \\&\\ P_{IHX\\ LP\\ in})}\n`}","pinned":false,"hidden":false,"mode":"md"},{"id":947,"value":"Plot.plot({\n y: { type: \"log\", label: \"P [Pa]\" },\n x: { label: \"H [J/kg]\" },\n title: \"PH diagram\",\n marks: [\n Plot.lineY(saturatedLiquid, { x: \"H\", y: \"P\", z: null, stroke: \"gray\" }),\n Plot.lineY(saturatedGas, { x: \"H\", y: \"P\", stroke: \"gray\" }),\n Plot.lineY([...cycle, cycle[0]], { x: \"H\", y: \"P\", stroke: \"green\" }),\n Plot.lineY([...cycle_ref, cycle_ref[0]], { x: \"H\", y: \"P\", stroke: \"blue\", marker:true, strokeDasharray:[2, 10]}),,\n Plot.dot(cycle, { x: \"H\", y: \"P\", fill: \"white\", stroke: \"green\", r: 10 }),\n Plot.text(cycle, { x: \"H\", y: \"P\", text: \"text\" })\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":893,"value":"viewof condT = Inputs.range([20, 50], {\n value: 40,\n step: 0.1,\n label: \"Hot Side Temperature [C]\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1471,"value":"viewof evapT = Inputs.range([10, -50], {\n value: Math.min(-20, condT - 1),\n step: 0.1,\n label: \"Evaporation Temperature [C]\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1560,"value":"viewof IHXeta = Inputs.range([0.0, 0.99], {\n value: 0.75,\n step: 0.01,\n max: 1,\n label: \"Efficiency of IHX\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1719,"value":"Plot.plot({\n title: \"COP Improvement with IHX\",\n x: { domain: [-10, 100], label: \"%\" },\n marks: [\n Plot.barX([{ x: improvementWithIHX }], { x: \"x\" }),\n Plot.text([{ x: improvementWithIHX }], {\n x: \"x\",\n text: (d) => `${d.x.toFixed(2)} %`,\n dx: 5,\n textAnchor: \"start\"\n })\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2185,"value":"Plot.plot({\n title: \"Q_cooling Improvement with IHX\",\n x: { domain: [-10, 100], label: \"%\" },\n marks: [\n Plot.barX([{ x: Q_improvementWithIHX }], { x: \"x\" }),\n Plot.text([{ x: Q_improvementWithIHX }], {\n x: \"x\",\n text: (d) => `${d.x.toFixed(2)} %`,\n dx: 5,\n textAnchor: \"start\"\n })\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2187,"value":"Plot.plot({\n title: \"Compressors W Deterioration with IHX\",\n x: { domain: [-10, 100], label: \"%\" },\n marks: [\n Plot.barX([{ x: W_deteriorationWithIHX }], { x: \"x\" }),\n Plot.text([{ x: W_deteriorationWithIHX }], {\n x: \"x\",\n text: (d) => `${d.x.toFixed(2)} %`,\n dx: 5,\n textAnchor: \"start\"\n })\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1511,"value":"viewof etaS = Inputs.range([0.5, 0.95], {\n value: 0.75,\n step: 0.01,\n label: \"Isentropic Efficiency of the Compressor\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2046,"value":"viewof etaS1 = Inputs.range([0.5, 0.95], {\n value: 0.75,\n step: 0.01,\n label: \"Isentropic Efficiency of the Compressor Without IHX\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2741,"value":"pressureSolution = subCoolingShouldBeOptimized\n ? nelderMead1D((p) => -get_P_high_in_trans_critical(p), pCrit / 1e6, {\n tol: 0.001\n })\n : null","pinned":true,"hidden":false,"mode":"ojs"},{"id":2919,"value":"pressureSolutionRef = subCoolingShouldBeOptimized\n ? nelderMead1D((p) => -get_P_high_in_trans_critical_ref(p), pCrit / 1e6, {\n tol: 0.001\n })\n : null","pinned":true,"hidden":false,"mode":"ojs"},{"id":2789,"value":"subCoolingShouldBeOptimized = isTransCritical || optimize_subc","pinned":true,"hidden":false,"mode":"ojs"},{"id":2370,"value":"condP = isTransCritical || optimize_subc\n ? pressureSolution.x\n : mat\n .unit(\n cp.PropsSI(\n \"P\",\n \"T|liquid\",\n mat.unit(condT - condSC, \"degC\").toNumber(\"K\"),\n \"Q\",\n 0,\n refrigerant\n ),\n \"Pa\"\n )\n .toNumber(\"MPa\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":2716,"value":"condPref = isTransCritical\n ? pressureSolutionRef.x\n : mat\n .unit(\n cp.PropsSI(\n \"P\",\n \"T|liquid\",\n mat.unit(condT - condSC, \"degC\").toNumber(\"K\"),\n \"Q\",\n 0,\n refrigerant\n ),\n \"Pa\"\n )\n .toNumber(\"MPa\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":1669,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":2002,"value":"viewof IHXcond = Inputs.range([1, 500], {value:50, step:0.1, label:'conductivty of IHX [W / (m K)]'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2008,"value":"viewof IHXcellSize = Inputs.range([5, 20], {value:10, step:0.1, label:'cell size of IHX [mm]'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2004,"value":"viewof IHXthickness = Inputs.range([0.1, 5], {value:1, step:0.1, label:'thickness of IHX [mm]'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1479,"value":"viewof evapPd = Inputs.range([0, 30000], {value:0, step:10, label:'Evaporator Pressure Drop [Pa]'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1485,"value":"viewof evapSH = Inputs.range([0,20], {value:5, step:0.1, label:'Evaporator Super Heating [K]'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2480,"value":"> **note:** due to the hot side temperature this cycle is **${isTransCritical ? \"trans-critical\" : \"sub-critical\"}** uses a **${typeOfHotSideHeatExchanger}** ","pinned":false,"hidden":false,"mode":"md"},{"id":2773,"value":"viewof optimize_subc = Inputs.toggle({\n label: \"Optimize sub_cooling\",\n disabled: true\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1503,"value":"viewof condSC = Inputs.range([0, 10], {\n value: isTransCritical ? 0 : 1,\n step: 0.1,\n label: \"Hot side Sub Cooling [K]\",\n disabled: subCoolingShouldBeOptimized\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":3016,"value":"viewof condSC_ref = Inputs.range([0, 10], {\n value: 2,\n step: 0.1,\n label: \"Hot side Sub Cooling [K] reference\",\n disabled: subCoolingShouldBeOptimized\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1495,"value":"viewof condPd = Inputs.range([0, 60000], {\n value: 0,\n step: 100,\n label: \"Hot Side Pressure Drop [Pa]\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1462,"value":"viewof mDot = Inputs.range([0.1, 300], {\n value: 88.1,\n step: 0.1,\n label: \"Mass Flow with IHX [kg/hr]\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2812,"value":"viewof transversal_area = Inputs.range([0, 10_000], {\n // L / A\n value: 1963,\n step: 0.1,\n label: \"Transversal area of heat exchanger\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1432,"value":"viewof refrigerant = Inputs.select(fluids, {\n value: \"R404A\",\n label: \"Refrigerant\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1845,"value":"viewof calcInputs = Inputs.textarea({spellcheck:false, label: \"Math Inputs\", value:\"2 + 2\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1847,"value":"{\n let results;\n try {\n results = mat.evaluate(calcInputs, scope);\n } catch (error) {\n return htl.html`

${error.toString()}

`;\n }\n return htl.html`

${formatCalc(results)}

`;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1395,"value":"state","pinned":true,"hidden":false,"mode":"ojs"},{"id":2514,"value":"# Appendix\n\nThe rest of the notebook are the source calculations used.","pinned":false,"hidden":false,"mode":"md"},{"id":2669,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":2603,"value":"function secant(f, x0, x1, tolerance = 1e-6, maxIter = 100) {\n const abs = Math.abs;\n for (let i = 0; i < maxIter; i++) {\n const f0 = f(x0);\n const f1 = f(x1);\n if (abs(f1) < tolerance) return x1;\n\n const x2 = x1 - f1 * ((x1 - x0) / (f1 - f0));\n x0 = x1;\n x1 = x2;\n }\n return x1;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2459,"value":"typeOfHotSideHeatExchanger = isTransCritical ? \"Gas Cooler\" : \"Condenser\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":1850,"value":"formatCalc = mat.typed(\"formatCalc\", {\n \"number | Unit | Matrix | Object | Matrix | Fraction | BigNumber | Complex\": (\n r\n ) => mat.format(r, 10),\n ResultSet: mat.typed.referToSelf(\n (self) => (r) => r.entries.map((res) => self(res)).join(\"\\n\")\n ),\n undefined: (r) => \"\",\n string: (r) => r,\n Help: (r) => r.toString()\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":965,"value":"tCrit = cp.PropsSI('Tcrit','',0,'',0, refrigerant)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2378,"value":"pCrit = cp.PropsSI(\"Pcrit\", \"\", 0, \"\", 0, refrigerant)","pinned":false,"hidden":false,"mode":"ojs"},{"id":2629,"value":"function get_P_high_in_trans_critical(p_high) {\n let scope = new Map([\n [\"fluid\", refrigerant],\n [\"isTranscritical\", isTransCritical],\n [\"pHigh\", mat.unit(p_high, \"MPa\")],\n // \"pHigh_ref\", mat.unit(p_high, \"MPa\")],\n [\"mDot\", mat.unit(mDot, \"kg/hour\")],\n [\n \"evap\",\n {\n T: mat.unit(evapT, \"degC\"),\n P_drop: mat.unit(evapPd, \"Pa\"),\n superHeating: mat.unit(evapSH, \"K\")\n }\n ],\n [\n \"cond\",\n {\n T: mat.unit(condT, \"degC\"),\n P_drop: mat.unit(condPd, \"Pa\"),\n subCooling: mat.unit(isTransCritical ? 0 : condSC, \"K\")\n }\n ],\n [\"etaS\", etaS],\n [\n \"IHX\",\n {\n epsilon: IHXeta\n }\n ],\n [\"isTransCritical\", isTransCritical]\n ]);\n preparations.evaluate(scope);\n cycleCalculation.evaluate(scope);\n return scope.get(\"evap_COP\");\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2729,"value":"function get_P_high_in_trans_critical_ref(p_high) {\n let scope = new Map([\n [\"fluid\", refrigerant],\n [\"isTranscritical\", isTransCritical],\n //[\"pHigh\", mat.unit(p_high, \"MPa\")],\n [\"pHigh_ref\", mat.unit(p_high, \"MPa\")],\n [\"mDot\", mat.unit(mDot, \"kg/hour\")],\n [\n \"evap\",\n {\n T: mat.unit(evapT, \"degC\"),\n P_drop: mat.unit(evapPd, \"Pa\"),\n superHeating: mat.unit(evapSH, \"K\")\n }\n ],\n [\n \"cond\",\n {\n T: mat.unit(condT, \"degC\"),\n P_drop: mat.unit(condPd, \"Pa\"),\n subCooling: mat.unit(isTransCritical ? 0 : condSC, \"K\")\n }\n ],\n //[\"etaS\", etaS],\n [\"etaS1\", etaS1],\n [\n \"IHX\",\n {\n epsilon: IHXeta\n }\n ],\n [\"isTransCritical\", isTransCritical]\n ]);\n preparations.evaluate(scope);\n secondCycleCalculation.evaluate(scope);\n return scope.get(\"evap_COP_ref\");\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2364,"value":"isTransCritical = mat.unit(condT, \"degC\").toNumber(\"K\") > tCrit","pinned":false,"hidden":false,"mode":"ojs"},{"id":2191,"value":"Q_improvementWithIHX = (res.Q.toNumber(\"W\") / res.Q_ref.toNumber(\"W\") - 1) * 100","pinned":false,"hidden":false,"mode":"ojs"},{"id":2189,"value":"W_deteriorationWithIHX = (res.W_comp.toNumber(\"W\") /\n res.W_comp_ref.toNumber(\"W\") -\n 1) *\n 100","pinned":false,"hidden":false,"mode":"ojs"},{"id":1608,"value":"cycle = res.state.map((e, i) => ({\n P: mat.number(e.P, \"Pa\"),\n H: mat.number(e.H, \"J/(kg)\"),\n T: mat.number(e.T, \"K\"),\n text: i + 1\n}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":1644,"value":"cycle_ref = res.state_ref.map((e, i) => ({\n P: mat.number(e.P, \"Pa\"),\n H: mat.number(e.H, \"J/(kg)\"),\n T: mat.number(e.T, \"K\"),\n text: i + 1\n}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":1732,"value":"improvementWithIHX = (res.COP / res.COP_ref - 1) * 100","pinned":false,"hidden":false,"mode":"ojs"},{"id":2341,"value":"tTriple = cp.PropsSI('Ttriple','',0,'',0,refrigerant)","pinned":false,"hidden":false,"mode":"ojs"},{"id":1603,"value":"temps = {\n const temps = mat.range(tTriple, tCrit, 5, true).toArray();\n if (temps.slice(-1) < tCrit) temps.push(tCrit);\n return temps;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1001,"value":"saturatedLiquid = temps.map( T => ({\n P:cp.PropsSI('P','T|liquid', T,'Q',0, refrigerant),\n H:cp.PropsSI('H','T|liquid', T,'Q',0, refrigerant),\n T:T\n}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":1010,"value":"saturatedGas = temps.map((T) => ({\n P: cp.PropsSI(\"P\", \"T|gas\", T, \"Q\", 1, refrigerant),\n H: cp.PropsSI(\"H\", \"T|gas\", T, \"Q\", 1, refrigerant),\n T\n}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":1399,"value":"state = res.state.map((e) => ({\n \"P [bar]\": e.P.toNumber(\"bar\"),\n \"T [C]\": e.T.toNumber(\"degC\"),\n \"H [J/kg]\": e.H.toNumber(\"J/kg\"),\n \"S [J/ (kg K)\": e.S.toNumber(\"J/kg/ K\"),\n \"D [kg/m^3]\": e.D.toNumber(\"kg/m^3\")\n}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":2165,"value":"scope = new Map([\n [\"fluid\", refrigerant],\n [\"isTransCritical\", isTransCritical],\n [\"shouldBeOptimized\", subCoolingShouldBeOptimized],\n [\"pHigh\", mat.unit(condP, \"MPa\")],\n [\"pHigh_ref\", mat.unit(condPref, \"MPa\")],\n [\"mDot\", mat.unit(mDot, \"kg/hour\")],\n [\n \"evap\",\n {\n T: mat.unit(evapT, \"degC\"),\n P_drop: mat.unit(evapPd, \"Pa\"),\n superHeating: mat.unit(evapSH, \"K\")\n }\n ],\n [\n \"cond\",\n {\n T: mat.unit(condT, \"degC\"),\n P_drop: mat.unit(condPd, \"Pa\"),\n subCooling: mat.unit(isTransCritical ? 0 : condSC, \"K\"),\n subCooling_ref: mat.unit(isTransCritical ? 0 : condSC_ref, \"K\")\n }\n ],\n [\"etaS\", etaS],\n [\"etaS1\", etaS1],\n [\n \"IHX\",\n {\n epsilon: IHXeta,\n thickness: mat.unit(IHXthickness, \"mm\"),\n cellSize: mat.unit(IHXcellSize, \"mm\"),\n k: mat.unit(IHXcond, \"W/(m K)\"),\n transversal_area: mat.unit(transversal_area, \"mm^2\")\n }\n ]\n])","pinned":false,"hidden":false,"mode":"ojs"},{"id":927,"value":"res = {\n preparations.evaluate(scope);\n cycleCalculation.evaluate(scope);\n IHX_calc.evaluate(scope);\n secondCycleCalculation.evaluate(scope);\n return {\n tlow: scope.get(\"evap\").T,\n COP: scope.get(\"evap_COP\"),\n COP_ref: scope.get(\"evap_COP_ref\"),\n state: scope.get(\"c\").toArray(),\n state_ref: scope.get(\"c_ref\").toArray(),\n Q: scope.get(\"Q_c\"),\n IHX: scope.get(\"IHX\"),\n Q_ref: scope.get(\"Q_c_ref\"),\n W_comp: scope.get(\"W_comp\"),\n W_comp_ref: scope.get(\"W_comp_ref\"),\n Q_IHX: scope.get(\"Q_IHX\")\n };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2997,"value":"preparations = mat.compile(`\n# Short function to get fluid properties\nprop(DesiredProperty, FluidState) = props(DesiredProperty, fluid, FluidState);\n\n# Define low and high pressure\nsat_gas = {'T|gas': evap.T, Q: 100%}\npLow = prop('P', sat_gas );\nsat_liq = {'T|liquid': cond.T, Q: 0%};\n`)","pinned":false,"hidden":false,"mode":"ojs"},{"id":916,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":1963,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":1636,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":898,"value":"$35","pinned":false,"hidden":false,"mode":"ojs"},{"id":57,"value":"# Import\n[mathjs](http://mathjs.org) & [coolprop](http://coolprop.org)","pinned":false,"hidden":false,"mode":"md"},{"id":2757,"value":"md`mathjs: ${mat.version}`","pinned":true,"hidden":false,"mode":"ojs"},{"id":1285,"value":"import {mat, cp} from \"@dvd101x/coolprop\"","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W36","stdlib":"1","titling":"auto"}