B
@B_chavarria
Paramétrico CR (Temperaturas)
import matplotlib.pyplot as plt
# Datos proporcionados
rangos = [20, 18, 16, 14, 12, 10, 8, 6, 4]
temp_pv = [90.32, 88.38, 86.44, 84.49, 82.55, 80.59, 78.64, 76.67, 74.7]
temp_absorbedor = [88.15, 86.43, 84.71, 83.00, 81.28, 79.55, 77.82, 76.09, 74.34]
temp_sustrato = [88.06, 86.34, 84.62, 82.90, 81.18, 79.46, 77.73, 75.99, 74.25]
temp_htf_salida = [87.42, 85.76, 84.10, 82.45, 80.79, 79.12, 77.46, 75.78, 74.1]
# Configuración del gráfico
Paramétrico CR (Temperaturas)
import matplotlib.pyplot as plt
# Datos proporcionados
rangos = [20, 18, 16, 14, 12, 10, 8, 6, 4]
temp_pv = [90.32, 88.38, 86.44, 84.49, 82.55, 80.59, 78.64, 76.67, 74.7]
temp_absorbedor = [88.15, 86.43, 84.71, 83, 81.28, 79.55, 77.82, 76.09, 74.34]
temp_sustrato = [88.06, 86.34, 84.62, 82.9, 81.18, 79.46, 77.73, 75.99, 74.25]
temp_htf_salida = [87.42, 85.76, 84.1, 82.45, 80.79, 79.12, 77.46, 75.78, 74.1]
# Configuración del gráfico
Paramétrico Ef Optica (Temperaturas)
import matplotlib.pyplot as plt
# Datos proporcionados
rangos = [0.9, 0.85, 0.8, 0.75, 0.7, 0.65, 0.6]
temp_pv = [81.44, 80.84, 80.23, 79.63, 79.02, 78.42, 77.82]
temp_absorbedor = [80.3, 79.76, 79.23, 78.7, 78.16, 77.63, 77.1]
temp_sustrato = [80.21, 79.67, 79.14, 78.6, 78.07, 77.54, 77.01]
temp_htf_salida = [79.85, 79.33, 78.81, 78.3, 77.79, 77.27, 76.76]
# Configuración del gráfico
Paramétrico Longitud del SRC-PVT (Temperaturas)
import matplotlib.pyplot as plt
# Datos proporcionados
rangos = [3, 6, 9, 12, 15, 18, 21, 24, 27, 30]
temp_pv = [83.82, 81.28, 80.68, 80.51, 80.51, 80.58, 80.69, 80.83, 80.99, 81.15]
temp_absorbedor = [80.19, 79.49, 79.51, 79.66, 79.84, 80.04, 80.24, 80.45, 80.66, 80.87]
temp_sustrato = [80.1, 79.4, 79.42, 79.56, 79.75, 79.95, 80.15, 80.36, 80.57, 80.77]
temp_htf_salida = [78.41, 78.76, 79.04, 79.29, 79.53, 79.76, 79.99, 80.21, 80.43, 80.65]
# Configuración del gráfico
Paramétrico Diametro (Temperaturas)
import matplotlib.pyplot as plt
# Datos proporcionados
rangos = [0.03, 0.025, 0.02, 0.015, 0.01, 0.005]
temp_pv = [80.59, 80.58, 80.61, 80.72, 81.14, 83.24]
temp_absorbedor = [79.55, 79.54, 79.56, 79.68, 80.1, 82.2]
temp_sustrato = [79.46, 79.45, 79.47, 79.59, 80.01, 82.1]
temp_htf_salida = [79.12, 79.12, 79.12, 79.12, 79.13, 79.14]
# Configuración del gráfico
Paramétrico AP (Temperaturas)
import matplotlib.pyplot as plt
# Datos proporcionados
rangos = [10, 20, 30, 40, 50, 60]
temp_pv = [78.96, 87.09, 95.17, 103.3, 111.4, 119.5]
temp_absorbedor = [78.15, 85.25, 92.45, 99.55, 106.75, 113.95]
temp_sustrato = [78.02, 85.19, 92.34, 99.5, 106.7, 113.9]
temp_htf_salida = [77.73, 84.66, 91.55, 98.45, 105.4, 112.3]
# Configuración del gráfico
Paramétrico Longitud SRC-PVT (Ef Elec-Term)
import matplotlib.pyplot as plt
# Datos
rangos = [3, 6, 9, 12, 15, 18, 21, 24, 27, 30]
eficiencia_electrica = [0.2656, 0.2562, 0.2504, 0.2461, 0.2428, 0.24, 0.2376, 0.2356, 0.2337, 0.232]
eficiencia_termica = [0.5505, 0.5728, 0.5912, 0.6078, 0.6236, 0.6388, 0.6536, 0.6681, 0.6824, 0.6965]
# Configuración del gráfico
fig, ax1 = plt.subplots(figsize=(10.67, 8)) # Controla el ancho del gráfico
Paramétrico AP (Ef Elec-Term)
import matplotlib.pyplot as plt
# Datos
rangos = [10, 20, 30, 40, 50, 60]
eficiencia_electrica = [0.2469, 0.2537, 0.2566, 0.258, 0.2586, 0.2587]
eficiencia_termica = [0.6072, 0.5753, 0.564, 0.5584, 0.5553, 0.5535]
# Configuración del gráfico
fig, ax1 = plt.subplots(figsize=(10.67, 8)) # Controla el ancho del gráfico
Paramétrico Diametro (Ef Elec-Term)
import matplotlib.pyplot as plt
# Nuevos datos
rangos = [0.03, 0.025, 0.02, 0.015, 0.01, 0.005]
eficiencia_electrica = [0.2488, 0.2488, 0.2488, 0.2488, 0.2486, 0.2476]
eficiencia_termica = [0.5969, 0.5969, 0.5969, 0.5969, 0.597, 0.5977]
# Configuración del gráfico
fig, ax1 = plt.subplots(figsize=(10.67, 8)) # Controla el ancho del gráfico
Paramétrico Ef Optica (Ef Elec-Term)
import matplotlib.pyplot as plt
# Nuevos datos
rangos = [0.9, 0.85, 0.8, 0.75, 0.7, 0.65, 0.6]
eficiencia_electrica = [0.2694, 0.2547, 0.24, 0.2252, 0.2105, 0.1956, 0.1808]
eficiencia_termica = [0.6441, 0.6103, 0.5766, 0.5429, 0.5093, 0.4757, 0.4421]
# Configuración del gráfico
fig, ax1 = plt.subplots(figsize=(10.67, 8)) # Controla el ancho del gráfico
Temp HTF de salida - CueMex
import matplotlib.pyplot as plt
# Datos proporcionados
meses = list(range(1, 13))
temp_htf_cuernavaca = [71.55, 87.84, 97.06, 86.09, 92.03, 90.56, 101.9, 88.03, 54.34, 47.6, 78.65, 71.97]
temp_htf_mexicali = [63.23, 82.52, 98.38, 109.7, 112.4, 115.5, 111.8, 106.4, 100.6, 81.51, 65.32, 58.37]
# Configuración del gráfico
plt.figure(figsize=(10.67, 8)) # Tamaño en pulgadas para una escala 4:3 y resolución 800x600
plt.plot(meses, temp_htf_cuernavaca, marker='o', color='blue', markersize=8, lines
Eficiencia Term Cue-Mex
import matplotlib.pyplot as plt
# Datos
meses = list(range(1, 13))
eficiencia_termica_cuernavaca = [0.5958, 0.5975, 0.603, 0.6028, 0.5992, 0.6011, 0.6018, 0.5975, 0.5918, 0.5921, 0.5925, 0.5908]
eficiencia_termica_mexicali = [0.5768, 0.5893, 0.5988, 0.6046, 0.6063, 0.6073, 0.6048, 0.6019, 0.5987, 0.5883, 0.5775, 0.5718]
# Configuración del gráfico
plt.figure(figsize=(10.67, 8)) # Tamaño en pulgadas para una escala 4:3 y resolución 800x600 (10.67 x 8 pulgadas)
Eficiencia Elec Cue-Mex
import matplotlib.pyplot as plt
# Datos
meses = list(range(1, 13))
eficiencia_cuernavaca = [0.2527, 0.245, 0.2408, 0.2459, 0.2431, 0.2438, 0.2385, 0.2449, 0.2494, 0.2546, 0.2493, 0.2525]
eficiencia_mexicali = [0.2565, 0.2475, 0.2401, 0.2349, 0.2336, 0.2322, 0.2339, 0.2364, 0.2392, 0.248, 0.2555, 0.2588]
# Configuración del gráfico
plt.figure(figsize=(10.67, 8)) # Tamaño en pulgadas para una escala 4:3 y resolución 800x600 (10.67 x 8 pulgadas)
Energía Term Cue-Mex
import matplotlib.pyplot as plt
# Datos
meses = list(range(1, 13))
energia_termica_cuernavaca = [6743088, 8911746, 10492632, 9158608, 11529427, 10117026, 12718343, 10610885, 4769532, 3553407, 7477019, 6937986]
energia_termica_mexicali = [5897151, 7940142, 10748241, 12355472, 14317028, 13886040, 13873084, 12288815, 10436601, 8480472, 6226605, 5177217]
# Convertir Wh a MWh
energia_termica_cuernavaca = [x / 1_000_000 for x in energia_termica_cuernavaca]
energia_termica_mexicali = [x / 1_000_000
Energía Elec Cue-Mex
import matplotlib.pyplot as plt
# Datos
meses = list(range(1, 13))
energia_elec_cuernavaca = [2859672, 3654675, 4189464, 3736312, 4677218, 4103054, 5039067, 4350135, 2013362, 1539120, 3146569, 2964654]
energia_elec_mexicali = [2622541, 3334548, 4310278, 4800336, 5516492, 5308891, 5365416, 4827144, 4168927, 3575184, 2754990, 2342763]
# Convertir Wh a MWh
energia_elec_cuernavaca = [x / 1_000_000 for x in energia_elec_cuernavaca]
energia_elec_mexicali = [x / 1_000_000 for x in energia_elec_mexic
Matriz singular
program solve_temperatures
implicit none
! Declaración de variables
real :: Tem_a, Tem_sky, Tem_f_in, Tem_pv, Tem_conc
real :: Tem_abs, Tem_sub, Tem_f_out, P_ele_pv, q_f
real :: tol
integer :: max_iter
! Variables de las ecuaciones internas
Gráfica de la parábola
import numpy as np
import matplotlib.pyplot as plt
# Definir la ecuación de la parábola en su forma estándar y el rango de valores de x
a = 1 / (4 * 0.6) # Calculamos 'a' basado en la distancia desde el vértice hasta el punto focal
x = np.linspace(-1, 1, 400)
y = a * x**2
# Graficar la parábola
plt.figure(figsize=(8, 6))
DM Cuer
import matplotlib.pyplot as plt
# Datos proporcionados
meses = list(range(1, 13))
temp_ambiente = [24.65, 23.85, 29.65, 27.85, 23.65, 26.75, 27.25, 23.85, 23.25, 24.25, 25.65, 23.35]
temp_cielo = [16.65, 15.85, 21.65, 19.85, 15.65, 18.75, 19.25, 15.85, 15.25, 16.25, 17.65, 15.35]
radiacion_total = [812.6, 942.5, 991.1, 953.5, 989.3, 971.3, 1020, 944.1, 827.1, 748.4, 805.5, 777.4]
radiacion_directa = [548.4, 746.5, 779.6, 673.5, 795.8, 740.1, 865.3, 749.3, 650.7, 507.8, 624.1, 574]
# Configura
DM Mex
import matplotlib.pyplot as plt
# Datos proporcionados
meses = list(range(1, 13))
temp_ambiente = [18.75, 22.85, 29.15, 33.55, 32.45, 37.75, 35.75, 35.25, 33.35, 24.05, 18.05, 18.55]
temp_cielo = [10.75, 14.85, 21.15, 25.55, 24.45, 29.75, 27.75, 27.25, 25.35, 16.05, 10.05, 10.55]
radiacion_total = [629.3, 813.5, 934.2, 1009, 1061, 1019, 990.8, 932.8, 887.5, 780.2, 652.5, 560.2]
radiacion_directa = [537.6, 706.1, 806.3, 877.9, 919.6, 892.5, 876.8, 823.9, 783.1, 680.6, 570.5, 485.2]
# Configura
Pot Cuer
import matplotlib.pyplot as plt
# Datos
meses = list(range(1, 13))
potencia_cuernavaca = [8313, 10975, 11262, 9936, 11606, 10826, 12381, 11012, 9737, 7758, 9337, 8694]
potencia_mexicali = [8273, 10486, 11618, 12372, 12889, 12433, 12306, 11688, 11237, 10128, 8746, 7533]
# Dividir los valores por 1000 para convertir de kW a MW
potencia_cuernavaca = [p/1000 for p in potencia_cuernavaca]
potencia_mexicali = [p/1000 for p in potencia_mexicali]