@B_chavarria

program main_ET_PTC_HTF implicit none ! ======= Constantes físicas y de geometría ======= Double Precision :: sigma, G_GHI Double Precision :: A_PTC, alpha_PTC, epsilon_PTC Double Precision :: A_PV, epsilon_PV, R_conv_PTC, R_conv_PV Double Precision :: A_abs, alpha_abs, epsilon_abs, R_conv_abs Double Precision :: R_cond_PV, R_cond_sub, R_cond_abs Double Precision :: m_dot_HTF, C_p_HTF, h_HTF

program main_ET_PTC_HTF implicit none ! ======= Constantes físicas y de geometría ======= Double Precision :: sigma, G_GHI Double Precision :: A_PTC, alpha_PTC, epsilon_PTC Double Precision :: A_PV, epsilon_PV, R_conv_PTC, R_conv_PV Double Precision :: A_abs, alpha_abs, epsilon_abs, R_conv_abs Double Precision :: R_cond_PV, R_cond_sub, R_cond_abs Double Precision :: m_dot_HTF, C_p_HTF, h_HTF

program main_ET_PTC_HTF implicit none Double Precision :: T_PTC integer :: iter1, iter2 ! ====== Paso 1: Calcular T_PTC con Newton-Raphson 1 variable ====== call calc_T_PTC(T_PTC, iter1) print *, "-------------------------------------------" print *, "Temperatura PTC calculada:", T_PTC, "K", T_PTC - 273.15d0, "°C" print *, "Iteraciones método 1:", iter1

program Newton3Vars_HTF implicit none integer, parameter :: dp = selected_real_kind(15,300) ! variables incógnitas real(dp) :: T_abs, T_sub, T_out ! datos y constantes real(dp) :: T_a, T_sky, T_PV, T_in_HTF, T_PTC real(dp) :: q_dot_sol_PV, P_PV real(dp) :: G_GHI, A_abs, alpha_abs, epsilon_abs, A_PV, epsilon_PV real(dp) :: sigma, R_conv_abs, R_cond_abs, R_cond_sub, R_cond_PV, R_conv_PV

program Newton2Vars_HTF implicit none real(8) :: T_abs, T_sub real(8) :: T_a, T_sky, T_PV, T_in_HTF, T_PTC real(8) :: G_GHI, A_abs, alpha_abs, epsilon_abs real(8) :: sigma, R_conv_abs, R_cond_abs, R_cond_sub, R_cond_PV real(8) :: m_dot_HTF, C_p_HTF, h_HTF real(8) :: f1, f2, J11, J12, J21, J22 real(8) :: dT_abs, dT_sub, detJ real(8) :: tol, UA_HTF, NTU, eff_HTF

program Newton2Vars implicit none real(8) :: T_abs, T_sub, T_a, T_sky real(8) :: G_GHI, A_abs, alpha_abs, epsilon_abs real(8) :: sigma, R_conv_abs, R_cond_abs, R_cond_sub, R_conv_sub real(8) :: f1, f2, J11, J12, J21, J22 real(8) :: dT_abs, dT_sub, detJ real(8) :: tol integer :: iter, max_iter

program ET_PTC ! Declaración de variables real(8) :: T_PTC, T_PVi, T_a, T_sky real(8) :: G_GHI, A_PTC, alpha_PTC, epsilon_PV, epsilon_PTC real(8) :: sigma, A_PV, R_conv_PTC real(8) :: f_T_PTC, df_T_PTC, T_PTCe real(8) :: tol integer :: iter, max_iter

program newton_ptc implicit none ! Declaración de variables real(8) :: Tptc, Tpv, Ta, Tsky real(8) :: Gghi, Aptc, alpha_ptc, eps_pv, eps_ptc real(8) :: sigma, Apv, Rconv real(8) :: f, df, Tnew real(8) :: tol integer :: iter, max_iter

import numpy as np # Constantes físicas fundamentales q = 1.602e-19 # Carga elemental del electrón (Coulombs) k = 1.381e-23 # Constante de Boltzmann (Joules/Kelvin) # Función para ajustar parámetros según condiciones ambientales def ajustar_parametros(Jsc_ref, Voc_ref, G, T_cell, T_ref=25): """ Ajusta los parámetros fotovoltaicos a las condiciones reales de operación

import numpy as np # Constantes físicas fundamentales q = 1.602e-19 # Carga elemental del electrón (Coulombs) k = 1.381e-23 # Constante de Boltzmann (Joules/Kelvin) # Función para ajustar parámetros según condiciones ambientales def ajustar_parametros(Jsc_ref, Voc_ref, G, T_cell, T_ref=25): """ Ajusta los parámetros fotovoltaicos a las condiciones reales de operación

import numpy as np import matplotlib.pyplot as plt # Parámetros L = 0.9 # Longitud de los lados del triángulo (m) - Ampliado a 0.5 d_circulo = 0.45 # Diámetro del círculo (m) r_circulo = d_circulo / 2 # Radio del círculo (m) # Cálculo de la altura del triángulo h_triangulo = (np.sqrt(3) / 2) * L # Altura total del triángulo

import matplotlib.pyplot as plt # Datos de flujo másico, eficiencia eléctrica y térmica mass_flow = [0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18] electric_efficiency = [0.1607, 0.1715, 0.1752, 0.177, 0.1781, 0.1788, 0.1793, 0.1797, 0.18] thermal_efficiency = [0.5625, 0.5574, 0.5554, 0.5544, 0.5537, 0.5533, 0.553, 0.5527, 0.5526] # Tamaños de letra personalizables title_fontsize = 22 axis_label_fontsize = 20

import matplotlib.pyplot as plt # Datos de flujo másico, potencia eléctrica y térmica mass_flow = [0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18] electric_power = [1.542, 1.647, 1.682, 1.699, 1.709, 1.716, 1.721, 1.725, 1.728] thermal_power = [5.4, 5.351, 5.332, 5.322, 5.316, 5.312, 5.309, 5.306, 5.305] # Tamaños de letra personalizables title_fontsize = 22 axis_label_fontsize = 20

import matplotlib.pyplot as plt # Nuevos datos de diámetro y temperatura diameters = [0.02, 0.04, 0.06, 0.08, 0.1, 0.12, 0.14, 0.16, 0.18] pv_temp = [136.4, 103.5, 92.64, 87.18, 83.91, 81.71, 80.14, 78.96, 78.04] absorber_temp = [136.1, 103.3, 92.37, 86.92, 83.64, 81.45, 79.88, 78.7, 77.78] htf_inlet_temp = [134.5, 102, 91.23, 85.89, 82.7, 80.57, 79.06, 77.92, 77.04] # Configuración del tamaño de la figura fig, ax = plt.subplots(figsize=(12, 8))

import matplotlib.pyplot as plt import numpy as np # Datos de diámetro, eficiencia eléctrica y térmica diameter = [0.015, 0.03, 0.045, 0.06, 0.075, 0.09] electric_efficiency = [0.1813, 0.1734, 0.1684, 0.1646, 0.1614, 0.1585] thermal_efficiency = [0.5457, 0.5888, 0.6286, 0.6668, 0.7041, 0.7405] # Tamaños de letra personalizables title_fontsize = 22

import matplotlib.pyplot as plt # Datos de diámetro, potencia eléctrica y térmica diameters = [0.015, 0.03, 0.045, 0.06, 0.075, 0.09] electric_power = [1.74, 1.665, 1.617, 1.58, 1.549, 1.522] thermal_power = [5.239, 5.652, 6.034, 6.401, 6.759, 7.109] # Tamaños de letra personalizables title_fontsize = 22 axis_label_fontsize = 20

import matplotlib.pyplot as plt # Nuevos datos de diámetro y temperatura diameters = [0.015, 0.03, 0.045, 0.06, 0.075, 0.09] pv_temp = [78.8, 79.98, 81.61, 83.46, 85.46, 87.59] absorber_temp = [78.49, 79.84, 81.52, 83.4, 85.42, 87.55] htf_inlet_temp = [78.34, 79, 79.61, 80.19, 80.76, 81.32] # Configuración del tamaño de la figura fig, ax = plt.subplots(figsize=(12, 8))

import matplotlib.pyplot as plt # Nuevos datos de ancho, eficiencia eléctrica y térmica width = [0.03, 0.06, 0.09, 0.12] electric_efficiency = [0.1946, 0.1872, 0.1828, 0.1795] thermal_efficiency = [0.5061, 0.5245, 0.5394, 0.5529] # Tamaños de letra personalizables title_fontsize = 22 axis_label_fontsize = 20

import matplotlib.pyplot as plt # Datos de ancho, potencia eléctrica y térmica widths = [0.03, 0.06, 0.09, 0.12] electric_power = [1.868, 1.797, 1.755, 1.723] thermal_power = [4.859, 5.036, 5.178, 5.307] # Tamaños de letra personalizables title_fontsize = 22 axis_label_fontsize = 20

import matplotlib.pyplot as plt # Nuevos datos de ancho y temperatura widths = [0.03, 0.06, 0.09, 0.12] pv_temp = [78.91, 78.75, 78.92, 79.51] absorber_temp = [77.83, 78.21, 78.56, 79.25] htf_inlet_temp = [77.74, 78.02, 78.24, 78.45] # Configuración del tamaño de la figura fig, ax = plt.subplots(figsize=(12, 8))