Page 522 - Design and Operation of Heat Exchangers and their Networks
P. 522

Appendix  505


                     fprintf("thm =%f, tcm =%f, 1-thm =%f, e_h =%f\n", ...
                         thm, tcm, 1 - thm, epsilon_h);

                case 46
                     % analytical
                     nua2b = special_function_nu(aA, 2 ∗ bA);
                     nu1 = special_function_nu(aA, bA);
                     epsilon_h =1-2 ∗ nua2b + (2 - nu1) ∗ nu1;
                     % numerical
                     th(:, :) = 1;
                     tc(:, :) = 0;
                     ntuh = NTUh(1);
                     ntuc = NTUc(1);
                     [th, tc, thm, tcm] = crossflow_unmixed_unmixed ...
                          (ntuh, ntuc, nx, ny, th, tc);
                     fprintf("thm_i =%f, tcm_i =%f\n", thm, tcm);

                     ntuh = NTUh(2);
                     ntuc = NTUc(2);
                     th(:, :) = thm;
                     for i=1:1:nx1
                          tc(i, 1) = tc(i, ny1);
                     end
                     [th, tc, thm, tcm] = crossflow_unmixed_unmixed ...
                          (ntuh, ntuc, nx, ny, th, tc);
                     fprintf("thm =%f, tcm =%f, 1-thm =%f, e_h =%f\n", ...
                         thm, tcm, 1 - thm, epsilon_h);

                case 55
                     % analytical
                     mua = special_function_mu(aA, bA);
                     mub = special_function_mu(bA, aA);
                     epsilon_h = mua + mub / Rh(1);
                     % numerical
                     th(:, :) = 1;
                     tc(:, :) = 0;
   517   518   519   520   521   522   523   524   525   526   527