@B_chavarria

import matplotlib.pyplot as plt # Datos meses = list(range(1, 13)) potencia_termica_cuernavaca = [19603, 26761, 28206, 24357, 28611, 26693, 31247, 26862, 23107, 18040, 22186, 20346] potencia_termica_mexicali = [18603, 24969, 28971, 31844, 33451, 32520, 31819, 29755, 28131, 24024, 19767, 16647] # Dividir los valores por 1000 para convertir de kW a MW potencia_termica_cuernavaca = [p/1000 for p in potencia_termica_cuernavaca] potencia_termica_mexicali = [p/1000 for p in potencia_termica_mexical

import matplotlib.pyplot as plt # Datos meses = list(range(1, 13)) temp_pv = [72.26, 88.82, 98.1, 86.97, 93.09, 91.54, 103.1, 89.01, 79.36, 67.97, 79.46, 72.71] temp_absorbedor = [71.63, 87.92, 97.14, 86.17, 92.12, 90.65, 102.0, 88.11, 78.58, 67.39, 78.72, 72.04] temp_sustrato = [71.55, 87.84, 97.06, 86.09, 92.04, 90.56, 101.9, 88.03, 78.51, 67.33, 78.65, 71.97] temp_htf_salida = [71.55, 87.84, 97.06, 86.09, 92.04, 90.56, 101.9, 88.03, 78.51, 67.33, 78.65, 71.97] temp_htf_entrada = [24.65, 23.

import matplotlib.pyplot as plt # Datos meses = list(range(1, 13)) temp_pv = [63.92, 83.44, 99.45, 110.9, 113.7, 116.8, 113.0, 107.5, 101.6, 82.4, 66.05, 58.98] temp_absorbedor = [63.29, 82.59, 98.46, 109.8, 112.5, 115.6, 111.9, 106.5, 100.6, 81.58, 65.38, 58.41] temp_sustrato = [63.23, 82.52, 98.38, 109.7, 112.4, 115.5, 111.8, 106.4, 100.6, 81.51, 65.32, 58.37] temp_htf_salida = [63.23, 82.52, 98.38, 109.7, 112.4, 115.5, 111.8, 106.4, 100.6, 81.51, 65.32, 58.37] temp_htf_entrada = [18.75, 22.

program solve_tem_conc implicit none real :: G_tot, A_conc, alpha_conc, sigma, epsilon_pv, A_pv, Tem_pv real :: h_conc, Tem_a, epsilon_conc, Tem_sky, Tem_conc real :: tol integer :: max_iter ! Asignar valores a las variables conocidas G_tot = 578 A_conc = 94

program main implicit none ! -------- Valores de entrada -------- " Real :: G_tot = 578 !Radiación total incidente" Real :: G_b = 374 !Radiación directa" Real :: Tem_a = 282 !25 + 273.15, Temperatura ambiente" Real :: Tem_sky = 282 !25 + 273.15, Temperatura del cielo"

import numpy as np # Parámetros del problema G_tot = 1000 # Irradiación solar total (W/m^2) A_conc = 10 # Área del concentrador (m^2) alpha_conc = 0.8 # Absorbancia del concentrador sigma = 5.67e-8 # Constante de Stefan-Boltzmann (W/m^2 K^4) epsilon_pv = 0.9 # Emisividad de las células fotovoltaicas A_pv = 5 # Área de las células fotovoltaicas (m^2) h_conc = 10 # Coeficiente de transferencia de calor por convección (W/m^2 K)

import numpy as np from scipy.optimize import fsolve # Parámetros G_tot = 1000 # Irradiancia total en W/m² G_b = 800 # Irradiancia directa en W/m² A_conc = 13.77 # Área del concentrador en m² A_pv = 1.2 # Área del PV en m² A_abs = 0.6 # Área del absorbedor en m² alpha_conc = 0.03 # Absorptividad del concentrador

import matplotlib.pyplot as plt # Datos relacion_concentracion = [20, 18, 16, 14, 12, 10, 8, 6, 4] temp_pv = [84.09, 82.71, 81.33, 79.96, 78.59, 77.22, 75.85, 74.49, 73.14] temp_absorbedor = [82.54, 81.33, 80.12, 78.91, 77.7, 76.5, 75.3, 74.1, 72.91] temp_sustrato = [82.47, 81.26, 80.05, 78.84, 77.63, 76.43, 75.23, 74.03, 72.83] temp_htf_salida = [82.46, 81.24, 80.03, 78.83, 77.62, 76.42, 75.22, 74.02, 72.83] # Crear figura y ejes

import matplotlib.pyplot as plt import numpy as np # Datos eficiencia_optica = [0.9, 0.85, 0.8, 0.75, 0.7, 0.65, 0.6] eficiencia_electrica = [0.2916, 0.2754, 0.2592, 0.243, 0.2268, 0.2106, 0.1944] eficiencia_termica = [0.4881, 0.4627, 0.4372, 0.4118, 0.3863, 0.3608, 0.3354] # Crear figura y ejes fig, ax = plt.subplots()