rayleigh-flowlisted

# Rayleigh Flow (aerodynamics/high-speed/rayleigh-flow) Use when you must compute the state change of a perfect gas heated or cooled in a constant-area frictionless duct, the Rayleigh flow that closes the pair with the Fanno friction duct of the same wave. This leaf implements the closed-form Rayleigh-line station relations against the thermal-choking sonic state (the star state of the same-mass-flow duct, reached at M = 1), the maximum heat addition that thermally chokes the duct on either inlet branch, the exit Mach number after a prescribed heat addition recovered from the quadratic total-temperature balance, and the second-law entropy rise, in pure Python stdlib. It pairs with the same-pack isentropic-flow-relations (frictionless duct WITHOUT heat: total temperature frozen, isentropic throat A*) and the same-pack fanno-flow (adiabatic duct WITH wall friction) as the two non-isentropic constant-area duct mechanisms of the wave; heat in, no friction here. ## Domain quick reference Air at gamma = 1.4, R = 287.0 J/(kg K), cp = gamma*R/(gamma - 1) = 1004.5 J/(kg K), constant specific heat only. The star denotes the state of the same-mass-flow constant-area frictionless duct at M = 1, i.e. thermal choking, and f = 1 + gamma*M^2: - Static ratios: T/T* = M^2*(1+gamma)^2/f^2, p/p* = (1+gamma)/f, rho/rho* = f/(M^2*(1+gamma)) = (p/p*)/(T/T*) exactly. T/T* peaks at (1+gamma)^2/(4*gamma) = 36/35 = 1.02857142857 at M = 1/sqrt(gamma). - Total ratios: T0/T0* = (T/T*)*(1 + (gamma-