From blueprint to
permit-ready engineering

#!/usr/bin/env python3
"""
SYNC 3.0 Demand-Side Offtake Model
Hourly heat absorption simulation — coastal desal + district network
"""

import numpy as np
import pandas as pd
from dataclasses import dataclass
from typing import Dict, List, Tuple
from enum import Enum

class Season(Enum):
    WINTER = 0
    SPRING = 1
    SUMMER = 2
    FALL   = 3

class DemandSideOfftakeModel:
    """
    Hourly simulation engine.
    Tests whether demand curve absorbs campus rejection profile.
    """

    def calculate_rejection(self, P_it, T_h=323.0, T_c=298.0):
        """Calculate thermal rejection from IT load (reversed Carnot)"""
        cop_rev = T_c / (T_h - T_c)
        return P_it / cop_rev

    def run_simulation(self, campus_capacity_mw=50.0,
                        redundancy_factor=1.3, pipe_efficiency=0.92,
                        allow_curtailment=True):
        """Run full 8760-hour annual simulation."""
        results = []
        for hour in range(8760):
            # ... seasonal profile lookup, IT load, rejection, demand absorption
            P_it = min(self.it_load_24h[hour % 24] * redundancy_factor,
                        campus_capacity_mw)
            Q_rej = self.calculate_rejection(P_it) * pipe_efficiency
            Q_abs = Q_desal + Q_district + Q_ag + Q_other
            Q_net = Q_abs - Q_rej
            # curtailment if demand insufficient
            curtailment_mw = abs(Q_net) / 0.25 if Q_net < 0 else 0.0
            results.append({ ... })
        return pd.DataFrame(results), self._generate_statistics(df)

def calculate_diffuser_dilution(
    Q_brine_m3s: float = 0.15,
    S_brine_ppt: float = 70.0,
    S_ambient_ppt: float = 35.0,
    T_brine_C: float = 35.0,
    T_ambient_C: float = 18.0,
    U_current_ms: float = 0.25,
    n_ports: int = 8,
    d_port_m: float = 0.15
) -> dict:
    """
    Calculate near-field brine plume dilution.
    Tests compliance with CA Ocean Plan 2 ppt / 100m limit.
    """
    # Buoyancy and momentum flux per port
    Q_port = Q_brine_m3s / n_ports
    v_port = Q_port / (math.pi * (d_port_m/2)**2)
    M = Q_port * v_port
    B = g * S_ratio * Q_port

    # Simplified near-field dilution (Fischer et al., 1979 / CORMIX)
    D1 = (Q_brine_m3s / (M**0.5)) * ((g * S_ratio * Q_brine_m3s /
          (U_current_ms**3))**(1/3))
    D_total = D1 * (n_ports**(2/3))  # port multiplier

    delta_S = (S_brine_ppt - S_ambient_ppt) / D_total
    delta_T = (T_brine_C  - T_ambient_C)   / D_total

    return {
        'D_calculated': round(D_total, 1),
        'delta_S_ppt':  round(delta_S, 3),   # → 0.23 ppt
        'delta_T_C':    round(delta_T, 3),   # → 0.11°C (flag)
        'S_compliant':  delta_S <= 2.0,       # True ✓
        'T_compliant':  delta_T <= 0.1,       # False ⚠
    }