import math import json from datetime import datetime, timedelta from astral import LocationInfo from astral.sun import sun, elevation, azimuth import pytz class BifacialSolarSimulation: def __init__(self): # Panel specifications (Hanersun N-TOPCon 700W) self.panel_power = 700 # watts self.panel_efficiency = 0.225 # 22.5% self.bifaciality_factor = 0.80 # Average of 75-85% self.panel_area = 2.384 * 1.303 # m² # Location: Lago Brown, Chile self.location = LocationInfo("Lago Brown", "Chile", "America/Santiago", -47.3919449, -72.3174973) self.timezone = pytz.timezone('America/Santiago') # Direct Normal Irradiance - constant when sun is shining (ideal conditions) self.dni_clear_sky = 900 # W/m² for clear sky conditions # Ground albedo (reflection factor) - clay soil year-round self.ground_albedo = 0.15 def get_solar_position(self, date, hour, minute=0): """Get precise solar position using astronomical calculations""" dt = datetime(date.year, date.month, date.day, hour, minute, 0) dt = self.timezone.localize(dt) # Get solar elevation and azimuth solar_elevation = elevation(self.location.observer, dt) solar_azimuth = azimuth(self.location.observer, dt) return solar_elevation, solar_azimuth def calculate_angle_of_incidence(self, solar_elevation, solar_azimuth, panel_tilt, panel_azimuth): """Calculate angle of incidence between sun and panel surface""" if solar_elevation <= 0: return 0 # Convert to radians se_rad = math.radians(solar_elevation) sa_rad = math.radians(solar_azimuth) pt_rad = math.radians(panel_tilt) pa_rad = math.radians(panel_azimuth) # Standard angle of incidence formula cos_incidence = ( math.sin(se_rad) * math.cos(pt_rad) + math.cos(se_rad) * math.sin(pt_rad) * math.cos(sa_rad - pa_rad) ) return cos_incidence def calculate_bifacial_irradiance(self, solar_elevation, solar_azimuth, panel_tilt, panel_azimuth): """Calculate irradiance with proper bifacial physics""" if solar_elevation <= 0: return 0, 0 # Standard angle of incidence for front side front_cos = self.calculate_angle_of_incidence(solar_elevation, solar_azimuth, panel_tilt, panel_azimuth) # Back side angle of incidence back_azimuth = (panel_azimuth + 180) % 360 back_cos = self.calculate_angle_of_incidence(solar_elevation, solar_azimuth, panel_tilt, back_azimuth) # Determine illumination scenario azimuth_to_front = abs(solar_azimuth - panel_azimuth) if azimuth_to_front > 180: azimuth_to_front = 360 - azimuth_to_front azimuth_to_back = abs(solar_azimuth - back_azimuth) if azimuth_to_back > 180: azimuth_to_back = 360 - azimuth_to_back # Determine primary illumination direction if azimuth_to_front < 90: # Front illumination (±90° from front) front_irr = self.dni_clear_sky * max(0, front_cos) back_irr = 0 elif azimuth_to_back < 90: # Back illumination - SPECIAL BIFACIAL PHYSICS front_irr = 0 if back_cos > 0: # For back illumination, vertical panels have significant advantage at low sun if solar_elevation < 30: # Low sun - bifacial geometry critical if panel_tilt == 90: # Vertical panel - best bifacial geometry bifacial_geometry_factor = 1.2 # 20% bonus for optimal geometry else: # Tilted panels lose efficiency for low sun back-illumination # More tilted = further from vertical = worse bifacial geometry vertical_deviation = abs(90 - panel_tilt) / 90.0 bifacial_geometry_factor = 1.0 - 0.4 * vertical_deviation # Stronger penalty else: bifacial_geometry_factor = 1.0 # High sun - less geometry effect back_irr = self.dni_clear_sky * back_cos * self.bifaciality_factor * bifacial_geometry_factor else: back_irr = 0 else: # Side illumination - minimal direct light front_irr = self.dni_clear_sky * max(0, front_cos) * 0.5 # Reduced efficiency back_irr = self.dni_clear_sky * max(0, back_cos) * self.bifaciality_factor * 0.5 return front_irr, back_irr def calculate_panel_irradiance(self, date, hour, minute, panel_tilt, panel_azimuth, panel_type='tilted'): """Calculate irradiance using simplified but physically accurate approach""" solar_elev, solar_az = self.get_solar_position(date, hour, minute) if solar_elev <= 0: return 0, 0 # No sun if panel_type == 'ew_vertical': # East-West vertical panel east_cos = self.calculate_angle_of_incidence(solar_elev, solar_az, panel_tilt, 90) west_cos = self.calculate_angle_of_incidence(solar_elev, solar_az, panel_tilt, 270) if east_cos > west_cos and east_cos > 0: front_direct = self.dni_clear_sky * east_cos back_direct = self.dni_clear_sky * west_cos * self.bifaciality_factor if west_cos > 0 else 0 else: front_direct = self.dni_clear_sky * west_cos if west_cos > 0 else 0 back_direct = self.dni_clear_sky * east_cos * self.bifaciality_factor if east_cos > 0 else 0 else: # North-South panel with proper bifacial physics front_direct, back_direct = self.calculate_bifacial_irradiance( solar_elev, solar_az, panel_tilt, panel_azimuth) # Diffuse irradiance (sky radiation) diffuse_factor = 0.15 sky_diffuse = self.dni_clear_sky * diffuse_factor front_diffuse = sky_diffuse * (1 + math.cos(math.radians(panel_tilt))) / 2 back_diffuse = sky_diffuse * (1 - math.cos(math.radians(panel_tilt))) / 2 # Ground reflected irradiance ground_reflected = self.dni_clear_sky * math.sin(math.radians(solar_elev)) * self.ground_albedo front_ground = ground_reflected * (1 - math.cos(math.radians(panel_tilt))) / 2 back_ground = ground_reflected * (1 + math.cos(math.radians(panel_tilt))) / 2 # Apply bifacial efficiency to back side indirect components angle_from_vertical = abs(90 - panel_tilt) bifacial_efficiency = 1.0 - (angle_from_vertical / 90.0) back_indirect_factor = self.bifaciality_factor * (0.4 + 0.6 * bifacial_efficiency) # Total irradiance front_irradiance = front_direct + front_diffuse + front_ground back_irradiance = back_direct + (back_diffuse + back_ground) * back_indirect_factor return front_irradiance, back_irradiance def calculate_monthly_production(self, panel_tilt, panel_azimuth, panel_type='tilted'): """Calculate monthly energy production for given panel configuration""" monthly_production = {} for month in range(1, 13): # Representative date for month (15th day) date = datetime(2024, month, 15) daily_production = 0 half_hourly_data = [] # Calculate for each half hour of the day (48 intervals) for half_hour in range(48): hour = half_hour // 2 minute = 0 if half_hour % 2 == 0 else 30 front_irr, back_irr = self.calculate_panel_irradiance(date, hour, minute, panel_tilt, panel_azimuth, panel_type) if front_irr > 0 or back_irr > 0: total_irradiance = front_irr + back_irr # Power calculation: Irradiance (W/m²) × Panel area × Panel efficiency power = (total_irradiance / 1000) * self.panel_power # Energy for 30 minutes (0.5 hour) daily_production += power * 0.5 # Calculate bifacial contribution percentage bifacial_contribution = (back_irr / total_irradiance * 100) if total_irradiance > 0 else 0 time_str = f"{hour:02d}:{minute:02d}" half_hourly_data.append({ 'time': time_str, 'solar_elevation': round(self.get_solar_position(date, hour, minute)[0], 1), 'front_irradiance': round(front_irr, 2), 'back_irradiance': round(back_irr, 2), 'total_irradiance': round(total_irradiance, 2), 'power': round(power, 1), 'bifacial_contribution': round(bifacial_contribution, 1) }) else: time_str = f"{hour:02d}:{minute:02d}" half_hourly_data.append({ 'time': time_str, 'solar_elevation': 0, 'front_irradiance': 0, 'back_irradiance': 0, 'total_irradiance': 0, 'power': 0, 'bifacial_contribution': 0 }) # Days in month days_in_month = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31][month-1] monthly_energy = daily_production * days_in_month / 1000 # kWh monthly_production[month] = { 'energy_kwh': round(monthly_energy, 1), 'daily_production': round(daily_production, 1), 'half_hourly_data': half_hourly_data } return monthly_production def run_simulation(self): """Run simulation for all configurations""" configurations = [ {'name': 'Vertical North-South', 'tilt': 90, 'azimuth': 0, 'type': 'ns_vertical'}, # Vertical, facing North {'name': 'Vertical East-West', 'tilt': 90, 'azimuth': 90, 'type': 'ew_vertical'}, # Vertical, E-W oriented {'name': 'Tilt 20°', 'tilt': 20, 'azimuth': 0, 'type': 'tilted'}, # Tilted facing North {'name': 'Tilt 30°', 'tilt': 30, 'azimuth': 0, 'type': 'tilted'}, # Tilted facing North {'name': 'Tilt 45°', 'tilt': 45, 'azimuth': 0, 'type': 'tilted'}, # Tilted facing North {'name': 'Tilt 60°', 'tilt': 60, 'azimuth': 0, 'type': 'tilted'}, # Tilted facing North {'name': 'Tilt 70°', 'tilt': 70, 'azimuth': 0, 'type': 'tilted'} # Tilted facing North ] results = {} for config in configurations: print(f" - {config['name']}") production = self.calculate_monthly_production(config['tilt'], config['azimuth'], config['type']) results[config['name']] = production return results def generate_html_report(self, results): """Generate comprehensive HTML report""" # Calculate annual totals annual_totals = {} for config_name, data in results.items(): annual_total = sum(data[month]['energy_kwh'] for month in range(1, 13)) annual_totals[config_name] = annual_total # Find best configuration best_config = max(annual_totals, key=annual_totals.get) best_value = annual_totals[best_config] months_czech = ['Leden', 'Únor', 'Březen', 'Duben', 'Květen', 'Červen', 'Červenec', 'Srpen', 'Září', 'Říjen', 'Listopad', 'Prosinec'] months_english = ['January', 'February', 'March', 'April', 'May', 'June', 'July', 'August', 'September', 'October', 'November', 'December'] html_content = f""" Analýza výkonu bifaciálního panelu - Lago Brown, Chile

🌞 Analýza výkonu bifaciálního panelu

Hanersun N-TOPCon 700W • Lago Brown, Chile • Astronomická simulace

Specifikace panelu a lokace

Panel Hanersun N-TOPCon

Výkon: 700 W
Účinnost: 22.5%
Bifacialita: 75-85%
Rozměry: 2384 × 1303 mm
Hmotnost: 38.7 kg

Lokace Lago Brown

Zeměpisná šířka: -47.39°
Zeměpisná délka: -72.32°
Nadmořská výška: 165 m
Klimatická oblast: Patagonská Chile
Albedo terénu: 0.15 (jíl)

📊 Roční souhrn výroby energie

""" # Add summary cards for config_name in ['Vertical North-South', 'Vertical East-West', 'Tilt 20°', 'Tilt 30°', 'Tilt 45°', 'Tilt 60°', 'Tilt 70°']: annual_value = int(annual_totals[config_name]) html_content += f"""

{config_name}

{annual_value}
kWh/rok
""" html_content += f"""
📈 Graf porovnání ročních výnosů by byl zde
Nejlepší výkon: {best_config} ({int(best_value)} kWh/rok)

📅 Měsíční výroba energie

""" # Monthly data table config_order = ['Vertical North-South', 'Vertical East-West', 'Tilt 20°', 'Tilt 30°', 'Tilt 45°', 'Tilt 60°', 'Tilt 70°'] for month in range(1, 13): month_name = months_english[month-1] row_class = 'highlight-row' if month in [5, 6, 7, 8] else '' html_content += f"" for config in config_order: value = results[config][month]['energy_kwh'] col_class = 'highlight-column' if config in ['Vertical North-South', 'Tilt 70°'] else '' if row_class and col_class: col_class += ' highlight-row' elif row_class: col_class = 'highlight-row' html_content += f"" html_content += "" # Annual total row html_content += "" for config in config_order: annual_value = annual_totals[config] col_class = 'highlight-column' if config in ['Vertical North-South', 'Tilt 70°'] else '' html_content += f"" html_content += "" html_content += """
MonthVertical North-SouthVertical East-WestTilt 20°Tilt 30°Tilt 45°Tilt 60°Tilt 70°
{month_name}{value}
ANNUAL TOTAL{annual_value:.1f}

💡 Poznámky k měsíční výrobě:

  • Hodnoty jsou v kWh za měsíc
  • Simulace používá přesné astronomické výpočty pro pozici slunce
  • Počítá s ideálními podmínkami: jasná obloha, bez mraků, bez zastínění
  • Bifaciální efekt je zahrnut ve všech konfiguracích
  • Panely směřují na sever (jižní polokoule)

⏰ Půlhodinová výroba (15. den v měsíci)

""" # Month tabs for i, month_name in enumerate(months_czech): active_class = 'active' if i == 0 else '' html_content += f"" html_content += "
" # Half-hourly data for each month for month in range(1, 13): month_name = months_czech[month-1] month_name_en = months_english[month-1] active_class = 'active' if month == 1 else '' html_content += f"""

Půlhodinová výroba - {month_name_en} (15. den)

""" # Find times with production half_hourly_data = {} for config in config_order: half_hourly_data[config] = results[config][month]['half_hourly_data'] # Get all daylight times (when sun is above horizon) all_times = [] date = datetime(2024, month, 15) for half_hour in range(48): hour = half_hour // 2 minute = 0 if half_hour % 2 == 0 else 30 solar_elev, solar_az = self.get_solar_position(date, hour, minute) if solar_elev > 0: # Sun is above horizon time_str = f"{hour:02d}:{minute:02d}" all_times.append((time_str, solar_elev, solar_az)) if all_times: for time_str, solar_elev, solar_az in all_times: html_content += f"" for config in config_order: power = 0 bifacial_contrib = 0 for time_data in half_hourly_data[config]: if time_data['time'] == time_str: power = time_data['power'] bifacial_contrib = time_data['bifacial_contribution'] break if bifacial_contrib > 0: html_content += f"" else: html_content += f"" html_content += "" else: html_content += "" html_content += "
Time (Elevation/Azimuth)Vertical North-SouthVertical East-WestTilt 20°Tilt 30°Tilt 45°Tilt 60°Tilt 70°
{time_str}
({solar_elev:.1f}°/{solar_az:.0f}°)		{power:.1f} W
({bifacial_contrib:.1f}%)		{power:.1f} W
Žádná výroba v tomto měsíci
" html_content += """

🎨 Legenda barevného označení:

9.7° - Elevace slunce (výška nad horizontem)
39° - Azimut slunce (směr od severu)
7.6% - Bifaciální příspěvek (% ze zadní strany panelu)

⚡ Poznámky k půlhodinové výrobě:

  • Hodnoty jsou ve wattech (W) pro každých 30 minut
  • Simulováno pro 15. den každého měsíce
  • Používá přesné astronomické výpočty pozice slunce
  • Půlhodinová data zahrnují bifaciální efekt zadní strany panelu
  • Zobrazeny pouze časy s výrobou energie
  • Elevace: čím vyšší, tím silnější slunce
  • Azimut: 0°=sever, 90°=východ, 180°=jih, 270°=západ
  • Bifaciální %: podíl energie ze zadní strany panelu

Astronomická simulace vytvořena pro panel Hanersun N-TOPCon 700W v lokalitě Lago Brown, Chile
Datum: 12.08.2025
Poznámka: Simulace používá přesné astronomické výpočty a předpokládá ideální podmínky - jasná obloha bez mraků

""" with open('/home/www/claudev.cz/orb3/solar_panel_report.html', 'w', encoding='utf-8') as f: f.write(html_content) def main(): simulation = BifacialSolarSimulation() print("🌞 Spouštím astronomickou simulaci bifaciálního panelu pro Lago Brown, Chile...") print(f"Panel: Hanersun N-TOPCon 700W Bifacial") print(f"Lokace: {simulation.location.latitude}°, {simulation.location.longitude}°") print("Konfigurace: Vertikální N-S, Vertikální E-W, Sklony 20°-70°") print() # Run simulation print("Running astronomical simulation for all configurations...") results = simulation.run_simulation() # Calculate annual totals annual_totals = {} for config_name, data in results.items(): annual_total = sum(data[month]['energy_kwh'] for month in range(1, 13)) annual_totals[config_name] = annual_total print() print("📊 Souhrn roční výroby energie:") for config, total in annual_totals.items(): print(f" {config}: {total:.1f} kWh/rok") # Generate HTML report simulation.generate_html_report(results) print() print("✅ HTML report byl vytvořen: solar_panel_report.html") # Find best configuration best_config = max(annual_totals, key=annual_totals.get) best_value = annual_totals[best_config] print(f"🎯 Nejlepší konfigurace: {best_config} ({best_value:.1f} kWh/rok)") return results, annual_totals if __name__ == "__main__": results, annual_totals = main()