| import assert from 'node:assert/strict'; |
| import { climates, defaults, baseline, configure, validateDesign, solarAt, beamShadeFraction, |
| thermalStep, ventilationAt, simulateDay, summarize, comparisonText, ROOM } from '../build/model.mjs'; |
| let count = 0; |
| function test(name, run) { run(); count++; console.log('PASS', name); } |
| const near = (a,b,tolerance=1e-8) => assert.ok(Math.abs(a-b) <= tolerance, `${a} differs from ${b} by more than ${tolerance}`); |
| test('Exact thermal integration and energy conservation', () => { |
| const input = { inside: 20, outside: 35, conductance: 90, capacity: 1e6, gains: 200, seconds: 600 }; |
| const result = thermalStep(input); |
| const equilibrium = 35 + 200/90; |
| near(result.temperature, equilibrium + (20-equilibrium)*Math.exp(-90*600/1e6)); |
| near(result.residualJ, 0, 1e-6); |
| near(thermalStep({...input, inside:25, outside:25, gains:0}).temperature, 25); |
| near(thermalStep({...input, conductance:0}).temperature, 20.12); |
| }); |
| test('Heat flow is signed and thermal storage slows changes', () => { |
| const input = { inside:25, outside:35, conductance:50, capacity:1e6, gains:0, seconds:600 }; |
| assert.ok(thermalStep({...input, conductance:150}).temperature > thermalStep(input).temperature); |
| assert.ok(thermalStep({...input, outside:15, conductance:150}).temperature < thermalStep({...input,outside:15}).temperature); |
| assert.ok(thermalStep({...input,capacity:5e6}).temperature < thermalStep(input).temperature); |
| }); |
| test('Night solar is zero for every facade and climate', () => { |
| for (const climate of Object.values(climates)) for (let orientation=0;orientation<=315;orientation+=45) { |
| for (const hour of [0,1,3,21,23]) { |
| const s = solarAt(climate,{...defaults,orientation},hour); |
| assert.equal(s.watts,0); assert.equal(s.beam,0); assert.equal(s.diffuse,0); |
| } |
| for(let hour=0;hour<24;hour+=0.125) { |
| const s=solarAt(climate,{...defaults,orientation},hour); |
| assert.ok(Number.isFinite(s.watts) && s.watts>=0); |
| } |
| } |
| }); |
| test('Horizontal overhang blocks more direct beam at high sun and greater depth', () => { |
| assert.ok(beamShadeFraction(.9,68,0)>beamShadeFraction(.9,30,0)); |
| assert.ok(beamShadeFraction(1.8,30,0)>beamShadeFraction(.9,30,0)); |
| assert.equal(beamShadeFraction(0,68,0),0); |
| assert.equal(beamShadeFraction(2.4,-3,0),0); |
| for (const climate of Object.values(climates)) for(let hour=0;hour<24;hour+=.25) { |
| assert.ok(solarAt(climate,{...defaults,shade:2.4},hour).watts <= solarAt(climate,{...defaults,shade:0},hour).watts + 1e-8); |
| } |
| }); |
| test('Purge boundaries, proportional ramp and actual infiltration', () => { |
| for(const [outdoor,inside] of [[20,25],[24.5,25],[26,26.5],[35,25]]) { |
| const v=ventilationAt(climates.humid,defaults,outdoor,inside); |
| assert.equal(v.purge,false); near(v.ach,ROOM.infiltrationACH); |
| } |
| const partial=ventilationAt(climates.humid,defaults,24,25.25); |
| const full=ventilationAt(climates.humid,defaults,24,25.5); |
| near(partial.openingFraction,.5); near(full.openingFraction,1); |
| assert.ok(full.ach>partial.ach); |
| }); |
| test('Complete configuration validation is atomic and rejects coercion', () => { |
| const original={design:{...defaults},climate:'humid',hour:16,lens:'heat'}; |
| const frozen=JSON.stringify(original); |
| for(const input of [ |
| {shade:1.5,lens:'invalid'}, {glazing:20,climate:'nowhere'}, {mass:50,hour:24}, |
| {shade:null},{mass:true},{glazing:'34'},{hour:NaN},{unknown:1},{orientation:135.1},[], |
| ]) { |
| assert.throws(()=>configure(original,input)); |
| assert.equal(JSON.stringify(original),frozen); |
| } |
| const changed=configure(original,{shade:1.5,climate:'arid',hour:23,lens:'air'}); |
| assert.equal(changed.design.shade,1.5); assert.equal(changed.hour,23); |
| assert.equal(JSON.stringify(original),frozen); |
| assert.deepEqual(validateDesign(defaults),defaults); |
| }); |
| test('Occupied duration is exactly twelve hours; gain integral has correct units', () => { |
| const s=simulateDay(climates.humid,defaults); |
| const result=summarize(s); |
| assert.equal(s.steps.filter(x=>x.occupied).length,72); |
| assert.equal(result.occupiedHours,12); |
| const constant={...s,steps:s.steps.map(row=>({...row,solarW:100}))}; |
| near(summarize(constant).solarKWh,2.4); |
| assert.equal(s.hourly.length,24); |
| }); |
| test('Periodic convergence is independent of starting temperature', () => { |
| for(const climate of Object.values(climates)) for(const design of [defaults,baseline,{...defaults,mass:100,ventilation:5}]) { |
| const cold=simulateDay(climate,design,{initialTemperature:5}); |
| const hot=simulateDay(climate,design,{initialTemperature:50}); |
| assert.ok(cold.converged && hot.converged); |
| near(summarize(cold).peak.temperature,summarize(hot).peak.temperature,.02); |
| for(let i=0;i<cold.steps.length;i++) near(cold.steps[i].temperature,hot.steps[i].temperature,.02); |
| } |
| const unsettled=simulateDay(climates.humid,defaults,{maximumDays:1,initialTemperature:5}); |
| assert.equal(unsettled.converged,false); assert.equal(unsettled.warmupDays,1); |
| }); |
| test('Ten-minute steps agree with five-minute steps within declared educational tolerance', () => { |
| for(const climate of Object.values(climates)) for(const design of [defaults,baseline,{...defaults,mass:0,ventilation:100}]) { |
| const ten=summarize(simulateDay(climate,design)); |
| const five=summarize(simulateDay(climate,design,{minutes:5})); |
| near(ten.peak.temperature,five.peak.temperature,.2); |
| near(ten.bandHours,five.bandHours,.5); |
| near(ten.meanACH,five.meanACH,.2); |
| near(ten.purgeHours,five.purgeHours,.5); |
| } |
| }); |
| test('Design grid remains finite, deterministic, energy-balanced and converged', () => { |
| let cases=0; |
| for(const climate of Object.values(climates)) for(const orientation of [0,90,180,270]) |
| for(const glazing of [12,72]) for(const shade of [0,2.4]) |
| for(const mass of [0,100]) for(const ventilation of [5,100]) { |
| const design={orientation,glazing,shade,mass,ventilation}; |
| const s=simulateDay(climate,design); |
| assert.ok(s.converged,JSON.stringify({climate:climate.name,design})); |
| assert.ok(s.steps.every(r=>Number.isFinite(r.temperature) && r.ach>=.35 && r.solarW>=0)); |
| assert.ok(s.steps.every(r=>Math.abs(r.residualJ)<1e-5)); |
| assert.ok(summarize(s).bandPercent>=0 && summarize(s).bandPercent<=100); |
| cases++; |
| } |
| assert.equal(cases,192); |
| assert.deepEqual(simulateDay(climates.humid,defaults),simulateDay(climates.humid,defaults)); |
| }); |
| test('Worse schemes are reported as higher, never as zero improvement', () => { |
| const reference={peak:{temperature:30}}; |
| assert.match(comparisonText({peak:{temperature:32}},reference),/2.0°C higher/); |
| assert.match(comparisonText({peak:{temperature:28}},reference),/2.0°C lower/); |
| }); |
| console.log(`\n${count} test groups passed, including 192 parameter-grid cases.`); |
|
|
|
|
