[{"_1":2,"_3":-5,"_4":-5},"loaderData",{"_5":6,"_7":8,"_9":10},"actionData","errors","root",{"_294":295,"_296":297},"routes/_app.x.$handle.$name",{"_96":95,"_31":44,"_104":86,"_26":105,"_13":14,"_106":-5,"_11":107,"_108":109,"_110":111,"_112":113,"_114":115,"_20":-5,"_116":-5,"_117":-5,"_118":24,"_119":120,"_121":58,"_122":-5,"_123":-5,"_124":125},"routes/_app.x.$handle.$name.branch.$branch.version.$version.diff",{"_11":12,"_13":14,"_15":16,"_17":-5,"_18":-5},"project",{"_42":43,"_31":44,"_37":45,"_46":-5,"_47":24,"_48":49,"_50":49,"_51":52,"_53":-5,"_18":54,"_55":56,"_57":58,"_59":-5,"_60":-5,"_61":62,"_63":64,"_65":-5,"_66":-5,"_67":-5,"_68":-5,"_69":70,"_71":70,"_72":73,"_74":73,"_75":73,"_76":70,"_77":70,"_78":79,"_80":81,"_82":83,"_84":24},"version",275,"changes",[19],"aiChatInfo","userId",{"_20":21,"_22":-7,"_23":24,"_25":24,"_26":27,"_28":27},"file",{"_29":30,"_31":32,"_33":34,"_35":36,"_37":38,"_39":14,"_40":41,"_23":24},"previous","deleted",false,"deletedHere","path","src/simulation/state.js","latestPath","fileId","8f57c610-8b55-11f1-aa55-1607ee4eb77e","name","state.js","fileType","script","content","import { HEIGHT, WIDTH } from \"../browser/config.js\";\nimport { createGrid } from \"../browser/grid.js\";\nimport { buildLampMask, lampHalfWidth } from \"./geometry.js\";\n\nexport function createSimulationState(settings) {\n const cellCount = WIDTH * HEIGHT;\n const wax = new Float32Array(cellCount);\n const temperature = new Float32Array(cellCount);\n const grid = createGrid(WIDTH, HEIGHT, wax);\n\n const state = {\n particles: [],\n wax,\n temperature,\n velocityX: new Float32Array(cellCount),\n velocityY: new Float32Array(cellCount),\n divergence: new Float32Array(cellCount),\n pressure: new Float32Array(cellCount),\n pressureNext: new Float32Array(cellCount),\n vesselMask: null,\n tick: 0,\n time: 0,\n diagnostics: {},\n indexOf: grid.indexOf,\n };\n\n state.vesselMask = buildLampMask(state.indexOf);\n seedParticles(state, settings);\n return state;\n}\n\nexport function seedParticles(state, settings) {\n state.particles.length = 0;\n const center = (WIDTH - 1) / 2;\n\n for (\n let y = HEIGHT - 7.4;\n y <= HEIGHT - 2.1 && state.particles.length < settings.particleCount;\n y += 0.52\n ) {\n for (\n let x = center - 5.6;\n x <= center + 5.6 && state.particles.length < settings.particleCount;\n x += 0.52\n ) {\n if (Math.abs(x - center) >= lampHalfWidth(y) - 0.5) continue;\n const n = state.particles.length;\n state.particles.push({\n x: x + 0.06 * Math.sin(n * 1.7),\n y: y + 0.04 * Math.cos(n * 2.1),\n vx: 0,\n vy: 0,\n temperature: 0.2 + 0.035 * Math.sin(x * 1.9),\n });\n }\n }\n\n state.tick = 0;\n state.time = 0;\n}\n","createdAt",["D",1785333501824],"branchVersion","branchId","821943c5-8c6e-4b33-9dd7-51e1353405e8","id","6fa4f715-b236-46ca-b183-2820945d7d51","laval",["D",1784854669144],"deletedAt","legacyEnvVars","privacy","public","httpPrivacy","description","Learn fluid simulation through a lava lamp.","imageUrl","40e69d1a-3254-4858-920b-7b1449761575","pinnedFileId","b0d26520-86fa-11f1-9ca8-1607ee4eb77e","pinPreview",true,"pinMigrated","screenshotVersion","likeCount",0,"orgId","54bae3d1-8223-41e5-b22d-81eb28a263c7","tursoName","tursoHostname","tursoMeasuredAt","pinned","userPerm","r","valPerm","blobPerm","rw","emailPerm","sqlitePerm","projectPerm","telemetryPerm","user",{"_96":95,"_103":99},"org",{"_96":95,"_97":-5,"_98":99,"_18":54,"_100":101,"_102":24},"branch",{"_13":85,"_31":86,"_42":41,"_87":43,"_88":-5,"_89":-5,"_90":-5,"_91":-5,"_92":-5,"_93":44,"_94":95},"isOwned",291,"main","projectId","forkedBranchId","forkedBranchVersion","forkedProjectPrivacy","forkedProjectId","forkedProjectOrgId","projectName","userHandle","pchinjr","handle","displayName","avatar","https://images.clerk.dev/oauth_github/img_2OaSpooJk6yZ5KFqyjtIVA4z3SH.jpeg","tier","free","canRestrictVals","profileImageUrl","branchName","","adoptOrgId",{"_42":43,"_31":44,"_37":290,"_46":-5,"_47":24,"_48":49,"_50":49,"_51":52,"_53":-5,"_18":54,"_55":56,"_57":58,"_59":-5,"_60":-5,"_61":62,"_63":64,"_65":-5,"_66":-5,"_67":-5,"_68":-5,"_69":70,"_71":70,"_72":73,"_74":73,"_75":73,"_76":70,"_77":70,"_78":291,"_80":292,"_82":293,"_84":24},"tree",[166,180,181],"files",[166,167,168,169,170,171,172,173,174,175,176,177,178,179,180,181],"counts",{"_155":62,"_156":157,"_158":62,"_159":160,"_161":162},"contents",[126,127,128],"upstreamProject","upstreamBranch","emptyEnvVars","flags",[],"hasHttpFiles","grants","orgs","customDomains",{},{"_42":148,"_31":149,"_29":150,"_131":-5,"_40":41,"_39":151,"_133":58,"_37":152,"_13":62,"_33":153,"_35":-5,"_23":24,"_138":-5,"_26":154},{"_42":140,"_31":141,"_29":142,"_131":-5,"_40":41,"_39":143,"_133":24,"_37":144,"_13":145,"_33":20,"_35":146,"_23":24,"_138":-5,"_26":147},{"_42":129,"_31":130,"_29":56,"_131":-5,"_40":41,"_39":132,"_133":58,"_37":134,"_13":135,"_33":136,"_35":137,"_23":24,"_138":44,"_26":139},"c81208e2-cbe0-40c9-9d9e-48d5d170c77d","main.ts","parentId",258,"latest",["D",1785252584636],145,"http","import { documentHtml } from \"./src/ui/document.ts\";\nimport { styles } from \"./src/ui/styles.ts\";\n\nconst headers = {\n html: { \"content-type\": \"text/html; charset=utf-8\" },\n css: { \"content-type\": \"text/css; charset=utf-8\" },\n};\n\nexport default function (request: Request) {\n const pathname = new URL(request.url).pathname;\n\n if (pathname === \"/styles.css\") {\n return new Response(styles, { headers: headers.css });\n }\n\n return new Response(documentHtml, { headers: headers.html });\n}\n","label",[130],"db1c70c7-2a1c-43f2-88e8-1f9c216ba72d","README.md","d9aa0084-86fa-11f1-87e4-1607ee4eb77e",269,["D",1785252677616],45,"# laval\n\nA step-by-step fluid simulation learning project for TypeScript developers.\n\n## Current checkpoint: particle volume pressure and modular architecture\n\nThis project builds a complete educational lava-lamp simulation one concept at a time. It begins with a simple Eulerian fluid solver, then adds thermal, material, diagnostic, shaped-vessel, and hybrid-particle models while keeping each operation behind a clearly named TypeScript function.\n\nThe live lessons now cover:\n\n1. **The grid** — select a cell and inspect its `(x, y)` coordinate and flat-array index.\n2. **A scalar field** — paint an abstract dye amount from `0` to `1`.\n3. **Diffusion** — observe values move toward neighboring values.\n4. **A velocity field** — store a two-dimensional vector in every cell and render it as an arrow.\n5. **Advection** — transport dye through the velocity field by tracing backward and sampling the previous scalar field.\n6. **Advection plus diffusion** — separate movement of the pattern from spreading of the pattern.\n7. **Divergence** — measure local imbalance between flow entering and leaving each cell.\n8. **Pressure projection** — solve a pressure field and use its gradient to reduce divergence.\n9. **Boundary conditions** — prevent velocity from pointing through the outer walls while preserving motion parallel to them.\n10. **Temperature and buoyancy** — store heat in another scalar field and convert it into upward velocity.\n11. **Rising heat gap** — reveal why local buoyancy followed directly by backward advection can pin bottom-row heat.\n12. **Coupled thermal step** — project and constrain velocity before using it to advect temperature.\n13. **Continuous heating and cooling** — add temperature at the heater while removing a fraction throughout the field every step.\n14. **Conduction and wall cooling** — exchange temperature with neighboring cells and remove additional heat at the glass boundary.\n15. **Wax material field** — store and advect wax fraction independently from temperature.\n16. **Temperature-driven density** — combine wax fraction and temperature into a signed buoyancy force.\n17. **Temperature-driven phase** — derive a liquid-wax fraction and use the solid fraction to resist motion.\n18. **Wax interface** — render mixed cells fractionally and counter numerical smearing with an explicit sharpening pass.\n19. **Latent heat** — store melt progress so phase can change while temperature remains near the melting point.\n20. **Surface tension** — estimate interface normals and curvature, then apply a curvature force near the wax boundary.\n21. **Phase-dependent viscosity** — diffuse velocity using a blend rate derived from carrier fraction, wax fraction, and melt progress.\n22. **Stable complete timestep** — place every operation in one authoritative order and subdivide frame time into smaller solves.\n23. **Continuous simulator** — schedule the stable timestep at a fixed rate with presets, pause/step controls, and live diagnostics.\n24. **Uneven heating** — start from a broad bottom wax pool and introduce deterministic heater variation that seeds localized buoyant regions.\n\n## Boundary model\n\nThe outer cells act as solid walls. `enforceBoundaries(velocity)` removes the component of velocity normal to each wall:\n\n- left and right walls set horizontal velocity to zero;\n- top and bottom walls set vertical velocity to zero;\n- the tangential component remains, so fluid may slide along the wall.\n\nThis is a simple **free-slip, no-through-wall** boundary condition. The lesson highlights outward-pointing wall cells in red and reports a leakage count. One boundary step should reduce that count to zero.\n\n## Current heat model\n\nThe temperature field stores a normalized heat value from `0` to `1` in every cell. `applyBuoyancy(velocity, temperature)` adds upward velocity in proportion to local temperature and the selected buoyancy strength. The isolated heat lesson intentionally changes velocity without moving temperature.\n\n## The rising-heat gap\n\n`stepRisingHeat()` composes two operations already introduced separately:\n\n1. `applyBuoyancy(velocity, temperature)` changes vertical velocity.\n2. `advect(temperature, velocity)` moves the temperature field through that velocity.\n\nThis exposes a missing link: local buoyancy may leave neighboring destination cells stationary, so backward advection cannot pull heat into them.\n\n## Current coupled step\n\n`stepCoupledThermal()` uses the fuller ordering:\n\n1. Apply buoyancy to velocity.\n2. Pressure-project velocity so the local push becomes a divergence-reduced surrounding flow.\n3. Enforce no-through-wall boundaries.\n4. Advect temperature through the corrected velocity.\n\nProjection does not move temperature itself. It prepares the velocity field that advection uses.\n\n## Current heat-cycle step\n\nEach step now performs two scalar-field updates before the coupled flow calculation:\n\n1. Multiply every temperature by `1 - coolingRate`.\n2. Add `heaterStrength` to cells above the bottom-center heater, clamped to `1`.\n\nThe result depends on a competition between energy input and loss. A stronger heater accumulates heat faster; stronger cooling shortens the lifetime of transported heat.\n\n## Current conduction step\n\nThe temperature update now distinguishes three mechanisms:\n\n1. **Ambient cooling** scales down temperature throughout the volume.\n2. **Conduction** moves each cell toward the average of its direct neighbors, even when velocity is zero.\n3. **Wall cooling** removes additional temperature only from edge cells touching the glass.\n\nAdvection remains different from conduction: advection carries heat with bulk motion; conduction spreads heat through local neighbor exchange.\n\n## Current material-field step\n\nEvery cell now contains another scalar, `waxFraction`:\n\n- `0` means carrier liquid.\n- `1` means wax.\n- A value between `0` and `1` is a numerically mixed cell produced by sampling.\n\nTemperature and wax fraction are separate fields. Both are transported through the same velocity field, but temperature represents thermal state while wax fraction represents material identity.\n\n## Current density step\n\nWax buoyancy is derived rather than stored:\n\n```ts\ndensityAnomaly = waxFraction * (neutralTemperature - temperature);\nverticalVelocity += densityAnomaly * densityResponse;\n```\n\nHot wax produces a negative canvas-y change and rises. Cold wax produces a positive canvas-y change and sinks. Carrier-only cells have `waxFraction = 0`, so this wax-specific force is zero regardless of their temperature.\n\n## Current phase step\n\nWax identity remains in `waxFraction`; phase is derived separately from temperature:\n\n```ts\nliquidFraction = smoothMelt(temperature, meltingTemperature, transitionWidth);\nsolidFraction = waxFraction * (1 - liquidFraction);\nvelocity *= 1 - solidFraction * solidDrag;\n```\n\nThe smooth transition avoids an abrupt numerical switch at one exact temperature. Cold wax resists motion, while hot wax becomes increasingly mobile without ceasing to be wax.\n\n## Current interface step\n\nSemi-Lagrangian advection and bilinear sampling produce fractional wax cells at a blob boundary. Rendering now blends the carrier and phase colors using the actual `waxFraction`. An optional sharpening pass applies a smoothstep curve and blends toward it:\n\n```ts\nsharpened = wax * wax * (3 - 2 * wax);\nwax = mix(wax, sharpened, interfaceSharpness);\n```\n\nThis makes the boundary more legible, but it is not perfectly mass-conserving.\n\n## Current latent-heat step\n\nEach wax cell now stores melt progress separately from temperature. Excess temperature above the melting point can be absorbed into melt progress; cooling wax can release that stored phase energy back into temperature. The result is a visible temperature plateau during melting and solidification.\n\nThis is an enthalpy-style teaching abstraction with normalized values, not a calibrated joule-based energy model.\n\n## Current surface-force step\n\nThe wax-fraction gradient estimates an interface normal. The divergence of that normalized direction estimates curvature. A continuum-surface-force approximation then adds velocity near the interface:\n\n```ts\nnormal = normalize(gradient(waxFraction));\ncurvature = -divergence(normal);\nforce = surfaceTension * curvature * gradient(waxFraction);\n```\n\nThis encourages coherent blobs by resisting ragged, highly curved interfaces.\n\n## Current viscosity step\n\nViscosity is represented as neighbor-based velocity diffusion. Each cell derives an effective per-step blend from its material and phase:\n\n```ts\nwaxViscosity = mix(coldWaxViscosity, hotWaxViscosity, liquidFraction);\neffectiveViscosity = mix(carrierViscosity, waxViscosity, waxFraction);\nvelocity = diffuseVelocity(velocity, effectiveViscosity);\n```\n\nThis smooths local momentum differences rather than merely scaling velocity toward zero.\n\n## Current complete timestep\n\nOne substep now follows this order:\n\n1. Add heater energy and apply ambient cooling.\n2. Conduct temperature and cool the walls.\n3. Exchange latent heat.\n4. Apply density, phase drag, and surface-tension forces.\n5. Diffuse velocity by material-dependent viscosity.\n6. Pressure-project velocity and enforce walls.\n7. Advect temperature, wax fraction, and melt progress.\n8. Sharpen the wax interface.\n\nA frame divides its requested time by the substep count. Rate-like changes and advection distance use that smaller duration.\n\n## Complete simulator runtime\n\nThe final module runs the stable timestep through a fixed-rate accumulator. Playback speed is measured in simulation frames per second and remains independent of display refresh rate. Presets set related parameters together, while diagnostics expose wax amount, average temperature, melt fraction, and maximum speed.\n\nThe simulator can be paused for cell inspection or advanced one frame at a time.\n\n## Current symmetry-breaking step\n\nA perfectly uniform pool over a perfectly uniform heater can remain artificially symmetric. The new heater scales its local input with a deterministic spatial and slowly changing perturbation. The perturbation changes only temperature; density and buoyancy amplify those small thermal differences into distinct regions of motion.\n\nAt zero unevenness, the heater remains uniform for direct comparison.\n\n## Current pinch-off step\n\nA 28 × 18 grid cannot resolve a wax neck once it becomes thinner than one cell. This lesson therefore makes that missing sub-grid decision visible:\n\n1. Detect a partially filled cell connecting thicker wax on opposite sides.\n2. Require velocity strain so a quiet narrow shape is not cut automatically.\n3. Remove a strength-scaled fraction from the neck.\n4. Return the removed wax to the thicker neighboring regions.\n5. Measure the remaining wax error after correction.\n\nCyan outlines mark detected necks. Recombination is already represented by the material field: when advection brings fractional interfaces into contact, they become one connected wax region and surface tension smooths the joined boundary.\n\nThis is a conservative topology correction for teaching, not a resolved Rayleigh–Plateau calculation.\n\n## Current conservation audit\n\nThe simulator now records the sum of every wax fraction before a substep. After advection, pinch-off, and interface sharpening, it measures the numerical drift and optionally redistributes that difference through cells that already contain wax:\n\n```ts\ntargetMass = sum(waxBefore);\nrunPhysics();\nrawDrift = sum(waxAfter) - targetMass;\nredistribute(-rawDrift, waxCells);\n```\n\nAt 0% correction, the raw numerical drift remains visible. At 100%, the corrected drift should remain near zero. This preserves normalized wax volume globally while allowing wax to move freely between individual cells.\n\n## Current energy-budget step\n\nThermal energy is approximated as the sum of temperature and melt progress. The ledger measures the system at four points in each substep:\n\n```ts\nstart = sum(temperature + meltProgress);\nheaterInput = afterHeater - start;\ncoolingLoss = afterCooling - afterHeater;\nnumericalDrift = finish - afterCooling;\n```\n\nTemperature and latent storage can exchange energy without changing their combined total. The total is not held constant because the lamp is an open system: its heater deliberately adds energy and its cooling boundaries deliberately remove it.\n\n## Current closed-vessel step\n\nEvery cell represents a full volume split between two complementary fractions:\n\n```ts\ncarrierFraction = 1 - waxFraction;\nwaxFraction + carrierFraction === 1;\n```\n\nThe earlier paint gesture created wax and therefore increased the field total before conservation began. In the closed-vessel lesson, painting an empty location withdraws the same amount from existing wax, preferring the bottom reservoir. The gesture changes the distribution without changing the inventory.\n\n## Current cycle-diagnostics step\n\nA single frame cannot establish whether the lamp is cycling. The simulator now records three wax-weighted measurements over time:\n\n- Center height: the vertical center of mass of the wax field.\n- Rising fraction: wax in cells with upward velocity.\n- Falling fraction: wax in cells with downward velocity.\n\nA repeating rise and fall in center height is evidence of circulation. A flat trace indicates a steady arrangement even when small local velocities remain.\n\n## Current blob-tracking step\n\nA blob is derived rather than stored. Cells above a configurable wax-fraction threshold form a binary mask, and a four-neighbor flood fill labels each connected component. The diagnostics report component count, the largest component’s share of all wax, and apparent split or merge events when that count changes.\n\nThe threshold matters: a weak fractional bridge may connect two regions at 20% but separate them at 50%. This makes blob identity a measurement choice layered over the continuous wax field.\n\n## Current persistent-identity step\n\nConnected components are now matched to the previous frame by cell overlap. The strongest overlap inherits the previous identity and outline color; an unmatched component receives a new identity. Multiple children sharing one prior identity indicate a split, while one region overlapping multiple prior identities indicates a merge.\n\n## What this does not claim\n\nThe pressure solution is approximate and the collocated grid is intentionally simple. The boundary operation does not yet model no-slip wall friction, internal obstacles, or boundary-aware pressure coupling. Buoyancy is the first generated force, but temperature does not yet move, spread, cool, or interact with material density. There is currently no wax/carrier-liquid phase interface.\n\n## Build order\n\n1. Grid and scalar field — complete\n2. Diffusion — complete\n3. Velocity vector field — complete\n4. Advection — complete\n5. Advection plus diffusion — complete\n6. Divergence — complete\n7. Pressure projection — complete\n8. Boundary conditions — complete\n9. Temperature and buoyancy — complete\n10. Rising heat gap — complete\n11. Coupled thermal step — complete\n12. Continuous heating and cooling — complete\n13. Temperature diffusion and wall cooling — complete\n14. Wax material field — complete\n15. Temperature-driven density — complete\n16. Temperature-driven phase change and solid drag — complete\n17. Wax interface — complete\n18. Latent heat — complete\n19. Cohesion and surface tension — complete\n20. Phase-dependent viscosity — complete\n21. Stable complete timestep — complete\n22. Continuous simulator and calibration — complete\n23. Uneven heating and bottom wax pool — complete\n24. Pinch-off and recombination — complete\n25. Global wax conservation — complete\n26. Thermal energy budget — complete\n27. Closed-vessel material inventory — complete\n28. Cycle diagnostics — complete\n29. Connected-component blob tracking — complete\n30. Persistent blob identity — complete\n31. Motion-tolerant identity matching — complete\n32. Coordinated cycle calibration — complete\n33. Physical-unit interpretation — complete\n34. Final streamlined simulator — complete\n35. Shaped-vessel circulation — complete\n36. Hybrid particle wax — complete\n37. Particle volume pressure — complete\n\nEach checkpoint adds one concept, one visualization, and one falsifiable claim.\n\n\n## Final simulator\n\nThe final mode runs the same authoritative physics assembled throughout the course. Motion-tolerant identity matching searches nearby previous cells when exact overlap fails. Calibration presets adjust related coefficients together. The units lesson interprets cells-per-frame using selected length and time scales without claiming the normalized coefficients are measured material constants.\n\nThe completed simulation includes a closed wax inventory, heating and cooling, conduction, latent heat, temperature-dependent buoyancy and phase drag, surface tension, material-dependent viscosity, pressure projection, wall boundaries, advection, pinch-off, wax-volume correction, cycle diagnostics, and persistent blob tracking.\n\n\n## Shaped-vessel circulation\n\nModule 36 replaces the rectangular active domain with a tapered lamp mask. Falling wax is confined toward a narrow base, the heater is concentrated around the center, and temperature is removed along the stair-stepped glass boundary. The complete substep still performs phase exchange, buoyancy, surface tension, viscosity, pressure projection, advection, pinch-off, and wax-volume correction.\n\nThe geometry boundary is an educational cell mask rather than a cut-cell pressure solve, but it directly tests the hypothesis that a shaped base and localized heat source help fallen wax pool, re-form, and rise repeatedly.\n\n\n### Bottom-collection correction\n\nThe first shaped-vessel attempt only cancelled velocity at stair-stepped walls. That prevented leakage but did not reproduce the inward pressure response of a genuinely sloped boundary, so fallen wax still remained a horizontal layer.\n\nModule 36 now applies a visible, adjustable inward collection force to cool wax near the base. It also corrects the base-concentration mapping: higher concentration now means a narrower Gaussian heater. This is an explicit coarse-grid correction for missing cut-cell wall physics.\n\n\n## Hybrid particle module\n\nModule 36 now represents wax as a fixed set of Lagrangian particles rather than a fraction stored in Eulerian cells. Each particle stores position, velocity, and temperature. Temperature changes its buoyant acceleration, pairwise cohesion keeps nearby wax together, and analytic collisions constrain particles to the tapered vessel.\n\nParticle velocities are deposited onto the existing grid, pressure-projected there, and sampled back to the particles. A temporary wax field is rasterized only for diagnostics and rendering; it is no longer the authoritative material inventory. Conservation follows from the fixed particle count.\n\n\n### Particle mass and thermal breakup\n\nThe shaped particle module now defaults to 220 wax particles and exposes an 80–320 particle inventory control. Changing the inventory resets the example because particle count is the conserved mass.\n\nBreakup is not driven by mass alone. The focused heater now moves slowly across the base and varies within the particle body, creating temperature and buoyancy differences. Cohesion is stronger for cold particles and weaker for hot particles, allowing heated lobes to stretch away while the cold reservoir remains available for later reheating.\n\n\n## Particle volume pressure\n\nModule 37 estimates each particle's local density with a compact smoothing kernel. When neighborhood density exceeds the selected rest density, symmetric outward pressure is applied before cohesion. A hard-core term also prevents extreme overlap.\n\nParticle velocities are now deposited bilinearly across four grid cells and projected velocities are sampled back bilinearly. This removes the nearest-cell quantization that previously gave an entire collapsed clump one shared velocity sample.\n\n\n## Source architecture\n\nThe browser now executes genuine JavaScript project files directly. No simulation function is stored in, reconstructed from, or served through a source-code string.\n\n- `main.ts` is the HTTP composition layer for the page and stylesheet.\n- `src/ui/document.ts` contains the HTML document.\n- `src/ui/styles.ts` contains the stylesheet.\n- `src/browser/config.js` exports the grid dimensions.\n- `src/browser/grid.js` exports readable grid indexing and scalar access functions.\n- `src/browser/lessons.js` exports the lesson explanations and pseudocode.\n- `src/browser/app.js` is the formatted browser entry module.\n\n`app.js` is organized into named sections for field setup, thermal and material physics, shaped-vessel particle physics, diagnostics and runtime, rendering and input, and the lesson controller. It imports the smaller modules normally through browser ES-module imports.\n\nThe mutable simulation fields remain in one closure for now. That is intentional: moving them into an explicit state object would change every physics function at once and make a behavior-preserving review difficult. The important boundary of this refactor is now real and inspectable: functions are ordinary formatted JavaScript, lesson data and grid utilities are genuine modules, and the previous string assembler has been removed.\n\n`src/contracts.ts` remains the conceptual type map for the fields and simulation stages.\n",[141],"ec02ba1f-9c0f-4799-880e-882f8edb6681","src","c305d0ec-86fa-11f1-ace5-1607ee4eb77e",2,["D",1784854699691],"directory",[149],"pulls","branches",1,"forks","envVars",{"_163":62,"_165":62},"connections",{"_163":62,"_164":62},"count","unfilledCount","emptyCount",{"_42":148,"_31":149,"_29":150,"_131":-5,"_40":41,"_39":151,"_133":58,"_37":287,"_13":62,"_33":153,"_35":-5,"_23":24,"_138":-5,"_183":288,"_26":289},{"_42":282,"_31":257,"_29":251,"_131":150,"_40":41,"_39":283,"_133":58,"_37":284,"_13":62,"_33":153,"_35":-5,"_23":24,"_138":-5,"_183":285,"_26":286},{"_42":274,"_31":275,"_29":276,"_131":251,"_40":41,"_39":277,"_133":24,"_37":278,"_13":151,"_33":34,"_35":279,"_23":24,"_138":-5,"_183":280,"_26":281},{"_42":266,"_31":267,"_29":268,"_131":251,"_40":41,"_39":269,"_133":58,"_37":270,"_13":62,"_33":34,"_35":271,"_23":24,"_138":-5,"_183":272,"_26":273},{"_42":258,"_31":259,"_29":260,"_131":251,"_40":41,"_39":261,"_133":58,"_37":262,"_13":62,"_33":34,"_35":263,"_23":24,"_138":-5,"_183":264,"_26":265},{"_42":248,"_31":249,"_29":250,"_131":251,"_40":41,"_39":252,"_133":24,"_37":253,"_13":62,"_33":34,"_35":254,"_23":24,"_138":-5,"_183":255,"_26":256},{"_42":243,"_31":226,"_29":222,"_131":150,"_40":41,"_39":244,"_133":58,"_37":245,"_13":62,"_33":153,"_35":-5,"_23":24,"_138":-5,"_183":246,"_26":247},{"_42":235,"_31":236,"_29":237,"_131":222,"_40":41,"_39":238,"_133":58,"_37":239,"_13":62,"_33":34,"_35":240,"_23":24,"_138":-5,"_183":241,"_26":242},{"_42":227,"_31":228,"_29":229,"_131":222,"_40":41,"_39":230,"_133":58,"_37":231,"_13":62,"_33":34,"_35":232,"_23":24,"_138":-5,"_183":233,"_26":234},{"_42":221,"_31":32,"_29":30,"_131":222,"_40":41,"_39":14,"_133":58,"_37":223,"_13":62,"_33":34,"_35":36,"_23":24,"_138":-5,"_183":224,"_26":225},{"_42":216,"_31":207,"_29":201,"_131":150,"_40":41,"_39":217,"_133":58,"_37":218,"_13":62,"_33":153,"_35":-5,"_23":24,"_138":-5,"_183":219,"_26":220},{"_42":208,"_31":209,"_29":210,"_131":201,"_40":41,"_39":211,"_133":24,"_37":212,"_13":157,"_33":34,"_35":213,"_23":24,"_138":-5,"_183":214,"_26":215},{"_42":198,"_31":199,"_29":200,"_131":201,"_40":41,"_39":202,"_133":24,"_37":203,"_13":62,"_33":34,"_35":204,"_23":24,"_138":-5,"_183":205,"_26":206},{"_42":189,"_31":190,"_29":191,"_131":150,"_40":41,"_39":192,"_133":24,"_37":193,"_13":194,"_33":34,"_35":195,"_23":24,"_138":-5,"_183":196,"_26":197},{"_42":140,"_31":141,"_29":142,"_131":-5,"_40":41,"_39":143,"_133":24,"_37":186,"_13":145,"_33":20,"_35":146,"_23":24,"_138":-5,"_183":187,"_26":188},{"_42":129,"_31":130,"_29":56,"_131":-5,"_40":41,"_39":132,"_133":58,"_37":182,"_13":135,"_33":136,"_35":137,"_23":24,"_138":44,"_183":184,"_26":185},["D",1785252584636],"children",[],[130],["D",1785252677616],[],[141],"b048c151-ff04-4baa-843a-72fdd4180858","contracts.ts","db9d981a-86fa-11f1-aed7-1607ee4eb77e",229,["D",1785250348184],40,"/**\n * Named contracts for the fluid-simulation learning checkpoints.\n *\n * Each interface exposes the operation a lesson teaches while allowing the\n * numerical implementation to remain hidden until the learner is ready for it.\n */\n\nexport type CellCoordinate = Readonly<{\n x: number;\n y: number;\n}>;\n\nexport interface ScalarField {\n readonly width: number;\n readonly height: number;\n at(cell: CellCoordinate): number;\n set(cell: CellCoordinate, value: number): void;\n clear(): void;\n}\n\nexport type Vector2 = Readonly<{\n x: number;\n y: number;\n}>;\n\nexport interface VectorField {\n readonly width: number;\n readonly height: number;\n at(cell: CellCoordinate): Vector2;\n set(cell: CellCoordinate, value: Vector2): void;\n clear(): void;\n}\n\nexport interface ScalarFieldRenderer {\n render(field: ScalarField): void;\n}\n\nexport interface VectorFieldRenderer {\n render(field: VectorField): void;\n}\n\nexport interface FieldUpdate {\n step(current: ScalarField): ScalarField;\n}\n\nexport interface DiffusionUpdate extends FieldUpdate {\n readonly strength: number;\n}\n\nexport interface AdvectionUpdate {\n readonly distancePerStep: number;\n step(quantity: ScalarField, velocity: VectorField): ScalarField;\n}\n\nexport interface TransportAndDiffusionStep {\n step(quantity: ScalarField, velocity: VectorField): ScalarField;\n}\n\nexport interface DivergenceMeasurement {\n measure(velocity: VectorField): ScalarField;\n}\n\nexport interface PressureProjection {\n solvePressure(divergence: ScalarField): ScalarField;\n subtractPressureGradient(\n velocity: VectorField,\n pressure: ScalarField,\n ): VectorField;\n}\n\n/**\n * A free-slip, no-through-wall boundary operation.\n *\n * It removes velocity normal to the outer walls while preserving velocity\n * tangent to them. A future no-slip condition can additionally remove the\n * tangential component.\n */\nexport interface BoundaryCondition {\n enforce(velocity: VectorField): VectorField;\n countOutwardLeaks(velocity: VectorField): number;\n}\n\n/** Converts a temperature scalar field into a vertical velocity change. */\nexport interface BuoyancyForce {\n readonly strength: number;\n apply(velocity: VectorField, temperature: ScalarField): VectorField;\n}\n\n/** Applies buoyancy, then transports temperature through the new velocity. */\nexport interface ThermalFlowStep {\n step(velocity: VectorField, temperature: ScalarField): Readonly<{\n velocity: VectorField;\n temperature: ScalarField;\n }>;\n}\n\n/** Orders force, projection, walls, and temperature transport into one step. */\nexport interface CoupledThermalStep {\n step(\n velocity: VectorField,\n temperature: ScalarField,\n projection: PressureProjection,\n boundaries: BoundaryCondition,\n ): Readonly<{ velocity: VectorField; temperature: ScalarField }>;\n}\n\n/** Adds heat at the lamp base and removes a fraction of heat everywhere. */\nexport interface ThermalSourceAndSink {\n readonly heaterStrength: number;\n readonly coolingRate: number;\n apply(temperature: ScalarField): ScalarField;\n}\n\n/** Separates neighbor-based conduction from additional boundary heat loss. */\nexport interface ThermalConduction {\n readonly diffusionStrength: number;\n readonly wallCoolingRate: number;\n diffuse(temperature: ScalarField): ScalarField;\n coolWalls(temperature: ScalarField): ScalarField;\n}\n\n/** A scalar fraction identifying wax independently from thermal state. */\nexport interface MaterialField extends ScalarField {\n /** 0 is carrier liquid, 1 is wax, intermediate values are mixed cells. */\n at(cell: CellCoordinate): number;\n}\n\nexport interface MaterialTransport {\n advect(material: MaterialField, velocity: VectorField): MaterialField;\n}\n\n/** Derives a signed wax-density force relative to the carrier liquid. */\nexport interface TemperatureDrivenDensity {\n readonly neutralTemperature: number;\n readonly response: number;\n apply(\n velocity: VectorField,\n temperature: ScalarField,\n waxFraction: MaterialField,\n ): VectorField;\n}\n\n/** Derives phase state from temperature without replacing material identity. */\nexport interface TemperatureDrivenPhase {\n readonly meltingTemperature: number;\n readonly transitionWidth: number;\n readonly solidDrag: number;\n liquidFraction(temperature: number): number;\n applyDrag(\n velocity: VectorField,\n temperature: ScalarField,\n waxFraction: MaterialField,\n ): VectorField;\n}\n\n/** Controls fractional-cell sharpening after material advection. */\nexport interface MaterialInterface {\n readonly sharpness: number;\n sharpen(waxFraction: MaterialField): MaterialField;\n}\n\n/** Stores phase energy separately so melting can occur at near-constant temperature. */\nexport interface LatentHeatExchange {\n readonly coupling: number;\n exchange(\n temperature: ScalarField,\n waxFraction: MaterialField,\n meltProgress: ScalarField,\n ): Readonly<{ temperature: ScalarField; meltProgress: ScalarField }>;\n}\n\n/** Applies a curvature force near the wax–carrier interface. */\nexport interface SurfaceTensionForce {\n readonly strength: number;\n apply(velocity: VectorField, waxFraction: MaterialField): VectorField;\n}\n\n/** Diffuses velocity at a rate derived from material identity and phase. */\nexport interface PhaseDependentViscosity {\n readonly carrierViscosity: number;\n readonly hotWaxViscosity: number;\n readonly coldWaxViscosity: number;\n diffuse(\n velocity: VectorField,\n waxFraction: MaterialField,\n meltProgress: ScalarField,\n ): VectorField;\n}\n\n/** Runs the authoritative coupled order using one or more stable substeps. */\nexport interface CompleteTimestep {\n readonly frameTime: number;\n readonly substeps: number;\n step(): void;\n}\n\nexport interface SimulatorRuntime {\n readonly playing: boolean;\n readonly stepsPerSecond: number;\n play(): void;\n pause(): void;\n step(): void;\n reset(): void;\n}\n\n/** Adds deterministic spatial variation to a distributed bottom heater. */\nexport interface UnevenHeatSource {\n readonly unevenness: number;\n apply(temperature: ScalarField, simulationStep: number): ScalarField;\n}\n\n/** Makes a coarse-grid topology correction while preserving total wax. */\nexport interface PinchOffCorrection {\n readonly neckSensitivity: number;\n readonly pinchStrength: number;\n detect(\n waxFraction: MaterialField,\n velocity: VectorField,\n ): ReadonlyArray;\n apply(waxFraction: MaterialField, velocity: VectorField): Readonly<{\n waxFraction: MaterialField;\n detectedNecks: number;\n conservationError: number;\n }>;\n}\n\n/** Audits and corrects global wax-volume drift across a complete substep. */\nexport interface WaxConservationAudit {\n readonly correctionStrength: number;\n total(waxFraction: MaterialField): number;\n correct(waxFraction: MaterialField, targetTotal: number): Readonly<{\n waxFraction: MaterialField;\n rawDrift: number;\n correctedDrift: number;\n }>;\n}\n\n/** Separates intentional thermal transfers from unexplained numerical drift. */\nexport interface ThermalEnergyBudget {\n total(temperature: ScalarField, meltProgress: ScalarField): number;\n step(): Readonly<{\n heaterInput: number;\n coolingLoss: number;\n latentStorage: number;\n numericalDrift: number;\n }>;\n}\n\n/** Relocates a fixed material inventory instead of creating wax while editing. */\nexport interface ClosedVesselInventory {\n totalWax(waxFraction: MaterialField): number;\n carrierFractionAt(cell: CellCoordinate): number;\n relocate(\n waxFraction: MaterialField,\n target: CellCoordinate,\n ): Readonly<{ waxFraction: MaterialField; moved: number }>;\n}\n\n/** Reduces a sequence of full fields to signals that reveal circulation. */\nexport interface CycleDiagnostics {\n readonly centerHeightHistory: ReadonlyArray;\n record(waxFraction: MaterialField, velocity: VectorField): Readonly<{\n centerHeight: number;\n risingFraction: number;\n fallingFraction: number;\n }>;\n}\n\nexport interface TrackedBlob {\n readonly label: number;\n readonly waxAmount: number;\n readonly cells: ReadonlyArray;\n}\n\n/** Derives connected wax regions from a thresholded material field. */\nexport interface BlobTracker {\n readonly connectionThreshold: number;\n identify(waxFraction: MaterialField): ReadonlyArray;\n}\n\n/** Matches derived regions across adjacent frames using spatial overlap. */\nexport interface PersistentBlobTracker {\n match(previous: ReadonlyArray, current: ReadonlyArray): Readonly<{\n blobs: ReadonlyArray;\n splits: number;\n merges: number;\n }>;\n}\n\nexport interface WaxParticle {\n position: Readonly<{ x: number; y: number }>;\n velocity: Readonly<{ x: number; y: number }>;\n temperature: number;\n}\n\n/** Couples conserved Lagrangian wax particles to an Eulerian pressure grid. */\nexport interface HybridParticleLamp {\n readonly particles: ReadonlyArray;\n depositVelocity(): VectorField;\n projectCarrierFlow(): VectorField;\n sampleProjectedVelocity(velocity: VectorField): void;\n step(): void;\n}\n\n/** Preserves finite particle volume with neighbor-density pressure. */\nexport interface ParticleVolumeConstraint {\n readonly restDensity: number;\n readonly pressureStrength: number;\n estimateDensity(particles: ReadonlyArray): ReadonlyArray;\n applyPressure(particles: ReadonlyArray): void;\n}\n\nexport const roadmap = [\n \"grid-and-scalar-field\",\n \"diffusion\",\n \"velocity-field\",\n \"advection\",\n \"advection-and-diffusion\",\n \"divergence\",\n \"pressure-projection\",\n \"boundaries\",\n \"temperature-and-buoyancy\",\n \"rising-heat-gap\",\n \"coupled-thermal-step\",\n \"continuous-heating-and-cooling\",\n \"temperature-diffusion-and-wall-cooling\",\n \"wax-material-field\",\n \"temperature-driven-density\",\n \"phase-change-and-viscosity\",\n \"two-phase-flow\",\n \"stable-complete-timestep\",\n \"continuous-simulator\",\n \"uneven-heating\",\n \"pinch-off-and-recombination\",\n \"global-wax-conservation\",\n \"thermal-energy-budget\",\n \"closed-vessel-inventory\",\n \"cycle-diagnostics\",\n \"connected-component-blob-tracking\",\n \"persistent-blob-identity\",\n \"motion-tolerant-matching\",\n \"cycle-calibration\",\n \"physical-unit-interpretation\",\n \"final-lava-lamp-simulator\",\n \"shaped-vessel-circulation\",\n \"particle-volume-pressure\",\n] as const;\n",[],[149,190],"d779e42e-99ec-4a68-86e3-d67878e23ba3","styles.ts","4cb17186-8a96-11f1-a9d3-1607ee4eb77e","34766f7c-8a96-11f1-8334-1607ee4eb77e",232,["D",1785251356153],"export const styles = \" :root { color-scheme:dark; --bg:#0d0f12; --panel:#15181d; --line:#303640; --text:#edf0f4; --muted:#929aa6; --blue:#50a7ff; --yellow:#ffd166; }\\n * { box-sizing:border-box; }\\n body { margin:0; background:var(--bg); color:var(--text); font:15px/1.5 ui-monospace,SFMono-Regular,Menlo,monospace; }\\n button,input,select { font:inherit; }\\n select { width:100%; color:var(--text); background:#242a32; border:1px solid #414955; border-radius:7px; padding:7px; }\\n .final-only { display:none; } .final-only.visible { display:inline-block; }\\n .shell { width:min(1050px,calc(100% - 28px)); margin:auto; padding:32px 0 60px; }\\n header { display:flex; justify-content:space-between; gap:30px; align-items:flex-end; margin-bottom:26px; }\\n h1 { margin:0; font:750 52px/.9 system-ui,sans-serif; letter-spacing:-.06em; }\\n h1 span { color:var(--blue); }\\n header p { max-width:560px; margin:0; color:var(--muted); }\\n .eyebrow { color:var(--yellow); text-transform:uppercase; letter-spacing:.14em; font-size:11px; margin-bottom:8px; }\\n .lesson-nav { display:grid; grid-template-columns:repeat(auto-fit,minmax(135px,1fr)); border:1px solid var(--line); border-radius:12px; overflow:hidden; margin-bottom:16px; }\\n .lesson-button { padding:13px; color:var(--muted); background:#111419; border:0; cursor:pointer; text-align:left; }\\n .lesson-button+.lesson-button { border-left:1px solid var(--line); }\\n .lesson-button b { display:block; color:inherit; margin-bottom:2px; }\\n .lesson-button small { font-size:10px; }\\n .lesson-button.active { color:var(--text); background:#20252c; box-shadow:inset 0 -3px var(--blue); }\\n .workspace { display:grid; grid-template-columns:minmax(0,1fr) 290px; gap:16px; align-items:start; }\\n .card { background:var(--panel); border:1px solid var(--line); border-radius:14px; overflow:hidden; }\\n .canvas-head { display:flex; justify-content:space-between; align-items:center; gap:14px; padding:12px 14px; border-bottom:1px solid var(--line); }\\n .canvas-head span { color:var(--muted); font-size:11px; }\\n .actions { display:flex; gap:7px; }\\n button.control { border:1px solid #414955; color:var(--text); background:#242a32; border-radius:7px; padding:6px 10px; cursor:pointer; }\\n button.control.primary { background:var(--blue); color:#07111b; border-color:var(--blue); font-weight:700; }\\n button:disabled { opacity:.35; cursor:not-allowed; }\\n .canvas-wrap { padding:14px; }\\n canvas { display:block; width:100%; aspect-ratio:14/9; background:#0b1017; border:1px solid #343b45; image-rendering:pixelated; cursor:crosshair; }\\n canvas.mini-chart { width:100%; height:92px; aspect-ratio:auto; cursor:default; margin-top:8px; }\\n .readout { display:grid; grid-template-columns:repeat(4,1fr); border-top:1px solid var(--line); }\\n .readout div { padding:10px 12px; min-width:0; }\\n .readout div+div { border-left:1px solid var(--line); }\\n .readout b { display:block; color:var(--yellow); white-space:nowrap; }\\n .readout small { color:var(--muted); font-size:9px; text-transform:uppercase; letter-spacing:.08em; }\\n aside { padding:17px; }\\n aside h2 { font:700 17px system-ui,sans-serif; margin:0 0 8px; }\\n aside p { color:var(--muted); margin:0 0 18px; }\\n .code { padding:13px; background:#0d1117; border:1px solid #2b323b; border-radius:9px; color:#c9d1d9; white-space:pre-wrap; font-size:12px; margin-bottom:18px; }\\n .code em { color:var(--blue); font-style:normal; }\\n .concept { border-top:1px solid var(--line); padding-top:14px; margin-top:14px; }\\n .concept b { display:block; margin-bottom:5px; }\\n .concept span { color:var(--muted); font-size:12px; }\\n .slider { display:none; margin-bottom:18px; }\\n .slider.visible { display:block; }\\n .slider label { display:flex; justify-content:space-between; margin-bottom:6px; font-size:12px; }\\n input[type=range] { width:100%; accent-color:var(--blue); }\\n .roadmap { margin-top:16px; padding:18px; }\\n .roadmap h2 { font:700 18px system-ui,sans-serif; margin:0 0 14px; }\\n .steps { display:grid; grid-template-columns:repeat(4,1fr); gap:8px; }\\n .step { border:1px solid var(--line); border-radius:9px; padding:11px; color:var(--muted); min-height:82px; }\\n .step.current { border-color:#397fbd; background:#142131; color:var(--text); }\\n .step small { color:var(--blue); }\\n .step b { display:block; font:650 13px system-ui,sans-serif; margin-top:6px; }\\n .truth { margin-top:14px; padding:12px; border-left:3px solid var(--yellow); background:#1d1b15; color:#bdb7a7; font-size:12px; }\\n @media(max-width:780px){ header{align-items:flex-start;flex-direction:column}.workspace{grid-template-columns:1fr}.steps{grid-template-columns:repeat(2,1fr)} }\\n @media(max-width:500px){.lesson-nav{grid-template-columns:1fr}.lesson-button+.lesson-button{border-left:0;border-top:1px solid var(--line)}.readout{grid-template-columns:repeat(2,1fr)}.readout div:nth-child(3){border-left:0;border-top:1px solid var(--line)}.readout div:nth-child(4){border-top:1px solid var(--line)}}\";\n",[],[149,207,199],"ui","eb9ea1d5-2259-4f0e-a158-7d6faa610780","document.ts","52b29326-8a96-11f1-bb6b-1607ee4eb77e",250,["D",1785252326240],"export const documentHtml = \"\\n\\n\\n \\n \\n laval — fluid simulation fundamentals\\n \\n\\n\\n

\\n

fluid simulation fundamentals

laval

\\n

Before simulating a lava lamp, learn how a computer represents a fluid. A sustained thermal flow needs an energy source and an energy sink. The heater adds temperature while ambient cooling continuously removes it.

\\n

\\n\\n

\\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n \\n

\\n\\n

\\n

Click a cell to inspect its coordinates.

\\n

—selected cell

\\n

—array index

\\n

—cell value

\\n

0update steps

\\n

\\n\\n

\\n

section 1

\\n

A fluid needs an address

\\n

We divide a rectangular region into cells. Each cell has an x and y coordinate, just like a pixel.

\\n

neighbor blend per step25%

Each cell closes this percentage of the difference between itself and its neighbors.

\\n

distance per step0.75 cells

How far backward each cell looks along its velocity vector.

\\n

pressure iterations60

More iterations approximate the pressure solution more closely.

mean |divergence|: —

\\n

wall leakage

outward-pointing wall cells: —

\\n

buoyancy strength0.50

How much upward velocity one step adds in hotter cells.

\\n

heater per step0.12cooling per step2%

The heater adds temperature at the bottom; cooling removes this fraction from every cell.

\\n

temperature diffusion15%wall cooling8%

Conduction exchanges heat with neighbors; wall cooling removes heat only at the glass boundary.

\\n

neutral temperature0.50density response1.20

Wax above neutral gets an upward push; wax below neutral gets a downward push.

\\n

melting temperature0.55transition width0.10solid drag85%

Cold wax resists velocity; hot wax becomes mobile across a smooth melting range.

\\n

interface sharpening50%

Pushes mixed wax fractions toward carrier (0) or wax (1) after advection.

\\n

latent heat coupling70%

How strongly excess heat near the melt point changes melt progress instead of temperature.

\\n

surface tension0.35

Scales the curvature force applied only near the wax–carrier interface.

\\n

carrier viscosity5%hot-wax viscosity10%cold-wax viscosity45%

Each value is the fraction of the difference toward neighboring velocity closed per step.

\\n

time per frame1.00substeps4

max travel per substep: —

\\n

lamp presetsteps per second6

wax —avg temp —melt —max speed —

\\n

heater unevenness60%

Creates deterministic hot and cool patches across the bottom heater.

\\n

neck sensitivity55%pinch strength45%

detected neck cells: 0 · wax error: 0.000

Cyan outlines show narrow bridges. Removed wax is returned to nearby thicker wax.

\\n

mass correction100%

raw drift 0.000after correction 0.000

Measures wax before the substep and returns numerical gain or loss across mixed interface cells.

\\n

energy ledger this substep

heater +0.000cooling 0.000stored latent 0.000numerical drift 0.000

Thermal energy is temperature plus energy stored in melt progress. Heating and cooling are intentional transfers; drift is not.

\\n

closed material inventory

total wax 0.0last moved 0.00

Clicking an empty cell withdraws the same wax amount from existing cells, preferring the bottom reservoir.

\\n

wax motion

height 0%rising 0%falling 0%

The blue trace is the wax center of mass. Repeated peaks and valleys indicate a circulating cycle.

\\n

wax connection threshold35%

blobs 0largest 0%split events 0merge events 0

Colored outlines identify cells belonging to the same four-neighbor connected region.

\\n

persistent identity matching

next identity 1matched 0

A region keeps its outline color when most cells overlap a region from the previous frame.

\\n

identity search radius1 cell

When direct overlap fails, search this far around the new region for its previous identity.

\\n

cycle quality

height range 0%active blobs 0

A useful lamp has a changing center height and more than one connected wax region without becoming unstable.

\\n

cell width0.5 cmseconds per frame0.10 s

reported maximum speed: 0.0 cm/s

\\n

complete model

blobs 0cycle range 0%rising 0%falling 0%

This runs the closed, mass-corrected, heated, cooled, two-phase timestep with pinch-off and persistent blob tracking.

\\n

wax particle mass220base concentration75%particle cohesion70%glass cooling10%

particles 0 · blobs 0

\\n

target neighbor density3.0volume pressure90%

maximum particle density: 0.0

\\n

const field = createScalarField(28, 18);\\n\\nfield.at(x, y);

\\n

What to noticeA two-dimensional cell can live in a one-dimensional array.

\\n

Deliberately not includedThere is no motion, pressure, viscosity, or fluid physics yet.

\\n

\\n\\n

\\n

The path from this grid to a lava lamp

\\n

checkpoint 01Grid + scalar field

\\n

02Velocity vectors

\\n

03Advection

\\n

04Diffusion

\\n

05Divergence

\\n

06Pressure projection

\\n

07Boundaries

\\n

08Heat + buoyancy

\\n

09Heat transport gap

\\n

10Coupled thermal step

\\n

11Heating + cooling

\\n

12Conduction + wall cooling

\\n

13Wax material field

\\n

14Temperature-driven density

\\n

15Phase change + solid drag

\\n

16Wax interface

\\n

17Latent heat

\\n

18Cohesion + surface tension

\\n

19Phase-dependent viscosity

\\n

20Stable complete timestep

\\n

21Continuous calibrated lamp

\\n

22Uneven heating + bottom pool

\\n

23Pinch-off + recombination

\\n

24Wax conservation

\\n

25Energy budget

\\n

26Closed vessel

\\n

27Cycle diagnostics

\\n

28Blob components

\\n

29Persistent identity

\\n

30Motion matching

\\n

31Cycle calibration

\\n

32Physical units

\\n

33Final simulator

\\n

34Shaped lamp vessel

\\n

35Particle volume pressure

\\n

Current claim: this demonstrates a grid, a scalar field, and neighbor-based smoothing. It is not yet a fluid simulation.

\\n

\\n\\n\";\n",[],[149,207,209],"f1340c10-1007-45af-9382-df0fe0c3c541",230,["D",1785251315458],[177,178],[149,207],"128ec34c-61c8-4b69-a4e1-582b0af09489","747a0682-8b55-11f1-8efb-1607ee4eb77e",["D",1785333501824],[],[149,226,32],"simulation","1d87a5e4-2b2f-4b04-8da4-fcf63eba3723","settings.js","8aa47a82-8b55-11f1-84b0-1607ee4eb77e",273,["D",1785333493953],"export const defaultSettings = {\n particleCount: 220,\n heaterStrength: 0.12,\n coolingRate: 0.02,\n glassCoolingRate: 0.1,\n heaterFocus: 0.8,\n neutralTemperature: 0.5,\n densityResponse: 1.2,\n particleCohesion: 0.7,\n restDensity: 3,\n volumePressure: 0.9,\n pressureIterations: 60,\n timeStep: 1,\n substeps: 4,\n playbackRate: 6,\n};\n\nexport function copySettings(overrides = {}) {\n return { ...defaultSettings, ...overrides };\n}\n",[],[149,226,228],"da779fb8-0d28-455c-ad85-0dcb7040e179","geometry.js","8cce612e-8b55-11f1-be8d-1607ee4eb77e",274,["D",1785333497586],"import { HEIGHT, WIDTH } from \"../browser/config.js\";\n\nexport function lampHalfWidth(y) {\n const t = Math.max(0, Math.min(1, y / (HEIGHT - 1)));\n return 4.5 + 6.2 * Math.pow(Math.sin(Math.PI * t), 0.72);\n}\n\nexport function buildLampMask(indexOf) {\n const mask = new Uint8Array(WIDTH * HEIGHT);\n const center = (WIDTH - 1) / 2;\n\n for (let y = 0; y < HEIGHT; y++) {\n for (let x = 0; x < WIDTH; x++) {\n mask[indexOf(x, y)] =\n y > 0 && y < HEIGHT - 1 && Math.abs(x - center) <= lampHalfWidth(y) ? 1 : 0;\n }\n }\n\n return mask;\n}\n\nexport function collideWithLamp(particle) {\n const center = (WIDTH - 1) / 2;\n const minimumY = 1.35;\n const maximumY = HEIGHT - 2.05;\n\n particle.y = Math.max(minimumY, Math.min(maximumY, particle.y));\n if (particle.y === minimumY && particle.vy < 0) particle.vy *= -0.18;\n if (particle.y === maximumY && particle.vy > 0) particle.vy *= -0.12;\n\n const horizontalLimit = lampHalfWidth(particle.y) - 0.48;\n if (particle.x < center - horizontalLimit) {\n particle.x = center - horizontalLimit;\n particle.vx = Math.abs(particle.vx) * 0.2;\n }\n if (particle.x > center + horizontalLimit) {\n particle.x = center + horizontalLimit;\n particle.vx = -Math.abs(particle.vx) * 0.2;\n }\n}\n",[],[149,226,236],"b99a0b6b-59b7-4f82-8ffb-9b85d7bec8a1",272,["D",1785333456767],[173,174,175],[149,226],"c64b9493-5588-4f6b-9dcc-553cd6e29e41","lessons.js","801014a4-8a98-11f1-add9-1607ee4eb77e","18484c56-8a98-11f1-acc0-1607ee4eb77e",248,["D",1785252301329],"export const copy={\n grid:{number:\"section 1\",title:\"A fluid needs an address\",text:\"We divide a rectangular region into cells. Each cell has an x and y coordinate, just like a pixel.\",instruction:\"Click a cell to inspect its coordinates.\",code:\"const field = createScalarField(28, 18);\\n\\nfield.at(x, y);\",notice:\"A two-dimensional cell can live in a one-dimensional array.\",notYet:\"There is no stored fluid quantity and no motion yet.\"},\n scalar:{number:\"section 2\",title:\"One number in every cell\",text:\"A scalar field stores one value at each location. Here the value represents an abstract amount of blue dye.\",instruction:\"Click and drag to paint dye into cells.\",code:\"field.set(x, y, dyeAmount);\\n\\nrenderScalarField(field);\",notice:\"Color is only a picture of the numbers. The Float32Array is the model.\",notYet:\"The dye does not move or spread by itself.\"},\n diffusion:{number:\"section 3\",title:\"State changes in steps\",text:\"A simulation repeatedly derives a new field from the current field. The slider controls how far each cell moves toward its neighbor average during the next step.\",instruction:\"For a fair comparison: choose a percentage, press example, then press step.\",code:\"while (running) {\\n field = diffuseScalarField(field);\\n renderScalarField(field);\\n}\",notice:\"The slider affects future steps; it cannot recalculate steps already taken. Compare 5% and 100% from the same example state.\",notYet:\"This is diffusion only—not velocity, advection, pressure, or a complete fluid solver.\"},\n velocity:{number:\"checkpoint 2\",title:\"Two numbers in every cell\",text:\"A velocity field stores a horizontal component vx and a vertical component vy at each location. An arrow visualizes that pair.\",instruction:\"Drag from one cell toward another to assign a direction, or press example.\",code:\"velocity.set(x, y, { vx, vy });\\n\\nrenderVectorField(velocity);\",notice:\"Arrow direction shows the vector direction. Arrow length represents its magnitude.\",notYet:\"The arrows describe possible motion, but nothing follows them yet. Transport comes next with advection.\"},\n advection:{number:\"checkpoint 3\",title:\"Use velocity to transport dye\",text:\"Advection connects the fields. For every destination cell, we trace backward along its velocity and sample the dye value found there.\",instruction:\"Press example, then advance the dye one step at a time.\",code:\"dye = advectScalarField(\\n dye,\\n velocity\\n);\\n\\nrenderScalarField(dye);\",notice:\"The velocity field is not the dye. It tells the update where each new dye value should come from.\",notYet:\"This uses semi-Lagrangian sampling. It is stable and visual, but it can blur dye and does not conserve its total amount exactly.\"},\n transportDiffusion:{number:\"checkpoint 4\",title:\"Transport, then spread\",text:\"A fluid step can apply multiple updates in sequence. First advection transports the dye. Then diffusion exchanges dye with neighboring cells.\",instruction:\"Press example, choose both strengths, and advance one combined step at a time.\",code:\"dye = advect(dye, velocity);\\ndye = diffuse(dye, diffusion);\\n\\nrenderScalarField(dye);\",notice:\"Advection changes where the pattern goes. Diffusion changes how sharp or spread-out the pattern is.\",notYet:\"This diffuses the dye quantity, not the velocity field. We also have not enforced incompressibility or solved pressure.\"},\n divergence:{number:\"checkpoint 5\",title:\"Where flow expands or compresses\",text:\"Divergence measures the net flow leaving each cell. Positive values behave like sources; negative values behave like sinks; values near zero have balanced inflow and outflow.\",instruction:\"Press example. Red cells are sources, blue cells are sinks. Click a cell to inspect its divergence.\",code:\"divergence = measureDivergence(velocity);\\n\\n// target for incompressible flow\\ndivergence ≈ 0;\",notice:\"A velocity field can look smooth while still creating or destroying volume. Divergence gives us a number that exposes that problem.\",notYet:\"We only measure the error here. The next lesson will solve pressure and subtract its gradient to reduce divergence.\"},\n projection:{number:\"checkpoint 6\",title:\"Use pressure to correct velocity\",text:\"Projection converts divergence into a pressure field, then subtracts the pressure gradient from velocity. The correction pushes away from compression and resists expansion.\",instruction:\"Press example, note the red and blue error, then press step to project the velocity.\",code:\"divergence = measureDivergence(velocity);\\npressure = solvePressure(divergence);\\nvelocity -= gradient(pressure);\",notice:\"Pressure is not painted manually. It is solved from the divergence error, used to correct velocity, and recalculated on the next step.\",notYet:\"The iterative pressure solve is approximate, and wall boundary conditions are still simplified. The next lesson will make boundaries explicit.\"},\n boundaries:{number:\"checkpoint 7\",title:\"Fluid cannot cross a solid wall\",text:\"A wall removes the velocity component pointing through it. Motion parallel to the wall can remain, so the fluid slides along the surface instead of escaping.\",instruction:\"Press example. Red wall cells leak outward; press step to enforce no-through boundaries.\",code:\"velocity = enforceNoThroughWalls(velocity);\\n\\nnormalVelocityAtWall = 0;\",notice:\"At left and right walls we remove horizontal velocity. At top and bottom walls we remove vertical velocity.\",notYet:\"This is a free-slip wall: tangential motion remains. We have not yet implemented no-slip friction, obstacles, or coupled projection at the boundary.\"},\n temperature:{number:\"checkpoint 8\",title:\"Heat creates upward acceleration\",text:\"Temperature is another scalar field. Hotter fluid is less dense, so buoyancy adds upward velocity where temperature is high.\",instruction:\"Press example, choose a buoyancy strength, then press step. Watch the arrows grow upward inside the hot region.\",code:\"velocity = applyBuoyancy(\\n velocity,\\n temperature\\n);\",notice:\"Temperature does not move anything directly. It changes velocity; velocity can move dye later through advection.\",notYet:\"Temperature does not diffuse or advect yet. We also omit ambient cooling, gravity calibration, density, and wax phase changes.\"},\n thermalFlow:{number:\"checkpoint 9\",title:\"The local-force gap\",text:\"Buoyancy changes velocity only where heat already exists. Backward advection in neighboring stationary cells cannot discover that heat, so a bottom row can remain pinned.\",instruction:\"Paint the bottom row and step. This deliberately preserves the gap you noticed.\",code:\"velocity = applyBuoyancy(velocity, temperature);\\ntemperature = advect(temperature, velocity);\",notice:\"Composing two correct operations is not enough if the velocity field has not been made into a coherent fluid-wide flow.\",notYet:\"This version skips pressure projection and boundary enforcement between force and transport. The next module inserts them.\"},\n coupledThermal:{number:\"checkpoint 10\",title:\"Spread the push before transport\",text:\"Pressure projection turns the local buoyant push into a divergence-reduced velocity field. Wall enforcement keeps that corrected flow inside before it transports temperature.\",instruction:\"Paint a hot bottom row or press example, then step. Compare its motion with Module 11.\",code:\"velocity = applyBuoyancy(velocity, temperature);\\nvelocity = project(velocity);\\nvelocity = enforceWalls(velocity);\\ntemperature = advect(temperature, velocity);\",notice:\"Projection does not move heat. It reshapes velocity so neighboring cells participate; advection then moves the temperature.\",notYet:\"The heater is not continuous and temperature does not diffuse or cool. Those processes come next.\"},\n heatedFlow:{number:\"checkpoint 11\",title:\"Temperature needs a source and a sink\",text:\"A real lamp keeps receiving energy from its bulb while losing energy to its surroundings. Each step adds heat near the bottom and removes a percentage everywhere.\",instruction:\"Press clear, then step repeatedly. Adjust heater and cooling to see whether heat accumulates or fades.\",code:\"temperature = cool(temperature);\\ntemperature = addBottomHeat(temperature);\\nvelocity = applyBuoyancyAndProject(velocity, temperature);\\ntemperature = advect(temperature, velocity);\",notice:\"The temperature field is no longer just initial data. It is driven over time by competing input and loss rates.\",notYet:\"Cooling is uniform rather than stronger at the glass, and heat does not diffuse by conduction. We still have only one material.\"},\n conductiveFlow:{number:\"checkpoint 12\",title:\"Heat also travels without bulk motion\",text:\"Conduction exchanges temperature between neighboring cells. The glass boundary then removes extra heat only from edge cells, creating a spatially meaningful heat sink.\",instruction:\"Press clear and step repeatedly. Compare low and high conduction, then change wall cooling.\",code:\"temperature = addBottomHeat(temperature);\\ntemperature = diffuseTemperature(temperature);\\ntemperature = coolAtWalls(temperature);\\nvelocity = solveFlow(velocity, temperature);\\ntemperature = advect(temperature, velocity);\",notice:\"Advection carries heat with moving fluid. Conduction spreads heat between neighbors even when the local velocity is zero.\",notYet:\"The walls all use one cooling rate and the heater has an abstract strength. We still have one continuous material rather than wax in a carrier liquid.\"},\n waxField:{number:\"checkpoint 13\",title:\"Material identity is another scalar field\",text:\"Wax fraction stores how much of each cell is wax: 0 is carrier liquid, 1 is wax, and values between represent a numerically mixed cell.\",instruction:\"Press example or paint wax, then step. Pink shows wax while the underlying warmth still affects the shared velocity.\",code:\"wax = advect(wax, velocity);\\ntemperature = advect(temperature, velocity);\",notice:\"Wax fraction and temperature coexist at every cell. Advection can transport both fields through the same velocity without confusing what they mean.\",notYet:\"Wax does not yet have different density or viscosity, and it does not melt or solidify. This is material tracking, not phase-change physics.\"},\n densityWax:{number:\"checkpoint 14\",title:\"Wax buoyancy depends on relative density\",text:\"Temperature changes wax density relative to the carrier. Wax hotter than neutral gets an upward force; colder wax gets a downward force.\",instruction:\"Press example: the left blob is hot and low, the right blob is cold and high. Step and compare their directions.\",code:\"densityError = wax * (neutralTemperature - temperature);\\nvelocity.y += densityError * response;\",notice:\"The force combines temperature and material identity. Equally hot carrier liquid does not receive this wax-specific density force.\",notYet:\"Values are normalized rather than kg/m³, and the relationship is linear. Wax still has no temperature-dependent viscosity or solid state.\"},\n phaseWax:{number:\"checkpoint 15\",title:\"Temperature changes how freely wax moves\",text:\"Wax identity remains constant, but a smooth phase fraction is derived from temperature. Cold wax is mostly solid and damps local velocity; hot wax is mostly liquid and follows the flow.\",instruction:\"Press example, then step. Watch the initially cold wax resist motion until heating moves it through the melting range.\",code:\"liquid = smoothMelt(temperature, meltPoint, width);\\nvelocity *= 1 - wax * (1 - liquid) * solidDrag;\",notice:\"Material identity and phase state are different: wax remains wax while its liquid fraction changes continuously from 0 to 1.\",notYet:\"Solid drag is a teaching approximation, not a full viscosity or viscoelastic stress solve. Latent heat and surface tension are still absent.\"},\n interfaceWax:{number:\"checkpoint 16\",title:\"A wax boundary is made of fractional cells\",text:\"Bilinear advection creates cells containing partial wax. The renderer now respects that fraction, and an optional sharpening pass nudges mixed values toward 0 or 1.\",instruction:\"Press example and step. Compare 0% and 100% sharpening, and inspect the wax value at the edge.\",code:\"wax = advect(wax, velocity);\\nwax = sharpenInterface(wax, strength);\",notice:\"Faint edge cells are numerical mixtures, not fully solid wax. Wax fraction must affect both the simulation and its rendering.\",notYet:\"This smoothstep sharpening is not mass-conserving and is not a full VOF, level-set, or surface-tension method.\"},\n latentWax:{number:\"checkpoint 17\",title:\"Melting consumes energy\",text:\"Wax now stores melt progress separately from temperature. Near the melting point, added energy advances the phase instead of immediately making the wax hotter.\",instruction:\"Press example, select a wax cell, and step. Watch temperature linger near the melt point while melt progress rises.\",code:\"energy = exchangeLatentHeat(temperature, meltProgress);\\n// temperature plateaus while phase changes\",notice:\"A phase transition can store energy without a matching temperature increase. Solidifying releases that stored energy back into temperature.\",notYet:\"The latent capacity is normalized and approximate; we do not yet use joules, specific heat, or a fully conservative enthalpy formulation.\"},\n surfaceWax:{number:\"checkpoint 18\",title:\"The wax boundary creates a force\",text:\"A gradient estimates the interface normal. Changes in that normal estimate curvature. Surface tension applies a force proportional to curvature near the wax boundary.\",instruction:\"Press example and step. Compare surface tension at 0 and 1; cyan cells mark the mixed interface.\",code:\"normal = normalize(gradient(wax));\\ncurvature = -divergence(normal);\\nforce = surfaceTension * curvature * gradient(wax);\",notice:\"Surface tension acts at the interface, not throughout the wax. It resists highly curved, ragged boundaries and encourages coherent blobs.\",notYet:\"This coarse continuum-surface-force approximation can create numerical currents and is not quantitatively calibrated.\"},\n viscosityWax:{number:\"checkpoint 19\",title:\"Viscosity diffuses momentum\",text:\"Each cell moves its velocity toward the average of neighboring velocities. The blend strength depends on whether the cell is carrier, hot liquid wax, or cold wax.\",instruction:\"Press example and step. Compare low and high cold-wax viscosity, then inspect viscosity at hot and cold wax cells.\",code:\"viscosity = mixByMaterialAndPhase(carrier, hotWax, coldWax);\\nvelocity = diffuseVelocity(velocity, viscosity);\",notice:\"Diffusing dye spreads a scalar quantity. Viscosity uses similar neighbor exchange on velocity, smoothing differences in momentum.\",notYet:\"These are normalized per-step blend rates, not Pa·s. Variable-viscosity diffusion is approximate and solid drag remains a separate teaching force.\"},\n stableLamp:{number:\"checkpoint 20\",title:\"One frame can contain several smaller solves\",text:\"The complete physics order is fixed. A frame time is divided by the selected substep count, and rate-like updates are scaled by that smaller duration.\",instruction:\"Press example. Compare one substep with eight at the same time per frame, and watch maximum travel per substep.\",code:\"for (const dt of subdivide(frameTime)) {\\n addHeatAndCool(dt);\\n updatePhaseAndForces(dt);\\n diffuseAndProjectVelocity(dt);\\n enforceWalls();\\n advectFields(dt);\\n}\",notice:\"Substeps do not add new physics. They reduce how far information travels before fields are recomputed, improving numerical stability.\",notYet:\"This still uses normalized units and explicit approximations. The final lesson adds continuous time, playback controls, and lamp-oriented presets.\"},\n simulator:{number:\"checkpoint 21\",title:\"The complete lamp runs continuously\",text:\"A fixed-rate animation loop repeatedly executes the stable timestep. Presets choose coherent groups of parameters rather than changing one coefficient in isolation.\",instruction:\"Choose a preset, press example, then play. Pause at any time to inspect a cell or advance one frame manually.\",code:\"while (playing) {\\n accumulator += elapsedTime;\\n while (accumulator >= frameInterval)\\n stepCompleteLamp();\\n}\",notice:\"The simulator runtime does not add physics. It schedules the validated timestep repeatedly and exposes diagnostics that reveal drift or instability.\",notYet:\"Parameters remain normalized and visually calibrated. This is a complete educational browser simulator, not a quantitatively validated industrial multiphase CFD model.\"},\n unevenHeating:{number:\"checkpoint 22\",title:\"Perfect symmetry prevents separate blobs\",text:\"A broad wax pool now rests above a heater whose strength varies across space and slowly over time. Small thermal differences seed separate buoyant regions.\",instruction:\"Press example, compare unevenness at 0% and 100%, then play. Watch for localized hot spots instead of uniform lifting.\",code:\"localHeat = heater * (1 + spatialPerturbation);\\ntemperature.bottom += localHeat;\",notice:\"The perturbation does not create motion directly. It creates temperature differences that density and buoyancy amplify into different rising regions.\",notYet:\"Pinch-off and merging still depend on the coarse interface method. The next step focuses on necking, separation, and recombination.\"},\n pinchOff:{number:\"checkpoint 23\",title:\"Thin wax bridges can become separate blobs\",text:\"The coarse grid cannot resolve a neck thinner than one cell, so this lesson adds an explicit topology rule: detect a strained, partially filled bridge, remove wax from its neck, and return that wax to its thicker neighbors.\",instruction:\"Press example, then play. Cyan cells are neck candidates. Compare pinch strength at 0% and 100%, and watch separated blobs merge naturally when advection brings their interfaces together.\",code:\"neck = detectThinBridge(wax, velocity);\\nremoved = pinch(neck);\\nredistributeConservatively(removed);\",notice:\"Pinch-off is a change in connectivity. The algorithm severs only narrow, strained bridges and conserves the wax it moves. Recombination needs no opposite rule: overlapping interfaces already become one connected field.\",notYet:\"This is a transparent teaching approximation, not a resolved Rayleigh–Plateau instability. A finer adaptive grid or particle/interface method would produce breakup from the governing forces themselves.\"},\n conservation:{number:\"checkpoint 24\",title:\"Wax may move, but its total amount should not change\",text:\"Advection and interface sharpening operate cell by cell, so rounding and nonlinear sharpening can create or destroy wax numerically. This lesson records total wax before each substep and audits it afterward.\",instruction:\"Press example and play. Set mass correction to 0% to expose raw drift, then 100% and compare the two diagnostics.\",code:\"targetMass = sum(waxBefore);\\nrunPhysics();\\nerror = targetMass - sum(waxAfter);\\nredistribute(error, interfaceCells);\",notice:\"Conservation is a constraint on the whole field, not a property of one cell. Individual cell values may change while their sum remains constant.\",notYet:\"The correction preserves normalized wax volume, not physical mass derived from density and cell volume. Energy and momentum also deserve separate conservation audits.\"},\n energyBudget:{number:\"checkpoint 25\",title:\"Heat has sources, sinks, and storage\",text:\"Unlike wax volume, thermal energy is not expected to remain constant. The heater adds it, ambient and wall cooling remove it, and melting transfers it between temperature and latent storage.\",instruction:\"Press example and play. Change heater and cooling strength, then watch which ledger entries respond. Stored latent energy can rise even while temperature pauses.\",code:\"energy = sum(temperature + meltProgress);\\nchange = heaterInput + coolingLoss + numericalDrift;\",notice:\"A conservation law is an accounting identity: system change equals energy crossing the boundary plus numerical error. Constant energy is expected only in a closed system.\",notYet:\"The ledger uses normalized energy units and omits kinetic and gravitational potential energy. The next refinement would audit the complete mechanical-plus-thermal budget.\"},\n closedVessel:{number:\"checkpoint 26\",title:\"A closed lamp has a fixed material inventory\",text:\"Previously, painting set a cell to wax without removing wax anywhere else. Conservation then faithfully preserved that newly enlarged total. In this lesson, painting relocates wax from the existing field instead.\",instruction:\"Press example. Note the total wax, then click or drag through empty cells. Wax leaves other locations—starting with the bottom reservoir—as the selected cells fill.\",code:\"amount = withdraw(waxInventory);\\ntarget.wax += amount;\\n// carrier fraction is always 1 - wax fraction\",notice:\"Every cell is already full: wax fraction plus carrier-liquid fraction equals one. Adding wax to a cell means displacing carrier there and taking wax from somewhere else in the closed lamp.\",notYet:\"The relocation gesture is an educational editing constraint, not fluid motion. During play, advection and pressure remain responsible for physical transport.\"},\n cycleDiagnostics:{number:\"checkpoint 27\",title:\"A cycle is a pattern measured over time\",text:\"One frame cannot tell us whether wax is circulating. We track the wax-weighted center height and the fractions currently moving upward or downward, then retain a short history.\",instruction:\"Press example and play. Look for the height trace to rise and fall repeatedly. A flat trace means the lamp has approached a steady arrangement, even if local cells still move.\",code:\"height = weightedCenter(wax);\\nrising = sum(wax where velocity.y < 0);\\nfalling = sum(wax where velocity.y > 0);\",notice:\"Oscillation is a property of a sequence of states, not one state. Center of mass compresses the whole wax field into one useful—but incomplete—signal.\",notYet:\"The trace does not identify individual blobs. Object tracking would be needed to measure each blob’s lifetime, speed, split, and merge events.\"},\n particleVolume:{number:\"checkpoint 35\",title:\"Wax particles need pressure as well as cohesion\",text:\"A smoothing kernel estimates crowding. Density above the target creates outward pressure, preventing many particles from occupying one point.\",instruction:\"Press example and play. Raise volume pressure until wax spans several cells, then tune cohesion to keep that finite volume connected.\",code:\"density = kernelSum(neighbors);\\npressure = max(0, density - restDensity);\",notice:\"Pressure preserves particle volume; cohesion shapes it. Bilinear transfer couples each particle to four grid cells.\",notYet:\"This is an SPH-style teaching approximation, not an iterative PBF solve.\"},\n shapedLamp:{number:\"checkpoint 34\",title:\"The vessel shape completes the circulation path\",text:\"Wax is a configurable fixed particle mass moving through the projected carrier grid. Uneven heating creates internal density differences, while hot particles weaken their cohesion enough for a large cluster to neck and separate.\",instruction:\"Press example and play. Compare low and high base concentration. Watch fallen wax gather centrally, reheat locally, and launch again.\",code:\"inside = lampShape(cell);\\nheat = focusedBaseHeater(x);\\nenforceCurvedWalls(velocity);\",notice:\"More particles create a deformable wax body rather than a tiny point cluster. Cold particles cohere strongly at the base; unevenly heated particles rise at different rates and can pull the body apart.\",notYet:\"The particle-to-grid transfer is deliberately simple and the rendered surface is an overlap visualization rather than a calibrated SPH kernel.\"},\n finalLamp:{number:\"checkpoint 33\",title:\"The complete educational lava lamp\",text:\"The final mode runs one authoritative timestep: heat, cooling, phase exchange, buoyancy, surface force, viscosity, pressure projection, walls, transport, pinch-off, and conservation.\",instruction:\"Choose a preset, press example, then play. Use the compact diagnostics to judge whether wax is forming blobs and completing a rise-and-fall cycle.\",code:\"while (playing) stepLavaLamp();\",notice:\"The final picture is the consequence of interacting fields and constraints, not a scripted blob animation.\",notYet:\"This remains a coarse normalized teaching model, not a quantitatively validated industrial multiphase CFD solver.\"},\n physicalUnits:{number:\"checkpoint 32\",title:\"Normalized values need a scale to become units\",text:\"A velocity of one means one cell per simulation frame. Choose a cell width and frame duration to interpret that same value as centimeters per second.\",instruction:\"Play, then change cell width and seconds per frame. The motion is unchanged; only its physical interpretation changes.\",code:\"speedCmPerSecond = velocityCellsPerFrame * cellWidthCm / secondsPerFrame;\",notice:\"Units do not appear automatically. A nondimensional simulation becomes physical only after selecting reference scales.\",notYet:\"The coefficients are visually calibrated and cannot yet be claimed as measured wax viscosity, density, or surface tension.\"},\n calibration:{number:\"checkpoint 31\",title:\"A lamp is calibrated as a system\",text:\"Heating, cooling, viscosity, density response, and surface tension must produce a bounded repeating cycle together. No single best slider value exists in isolation.\",instruction:\"Compare the balanced, slow, and active presets. Watch height range and blob count instead of judging only instantaneous speed.\",code:\"preset = calibrateTogether({heat,cooling,viscosity,density,surface});\",notice:\"Calibration chooses a useful regime of behavior; it does not add new governing physics.\",notYet:\"These presets are visual targets rather than fits to measurements from a specific manufactured lamp.\"},\n motionMatch:{number:\"checkpoint 30\",title:\"Identity matching needs motion tolerance\",text:\"A fast blob may move into entirely new cells between frames. When direct overlap fails, the tracker searches nearby previous cells before declaring a new identity.\",instruction:\"Play and compare a search radius of zero with one or two cells. Larger radii preserve identity through faster motion but can confuse nearby blobs.\",code:\"match = overlap || nearbyPreviousIdentity(searchRadius);\",notice:\"Tracking is an inference problem layered over the simulation. Greater tolerance reduces lost identities but increases ambiguous matches.\",notYet:\"The matcher expands spatially rather than predicting a velocity-directed future centroid.\"},\n blobIdentity:{number:\"checkpoint 29\",title:\"Identity comes from overlap across frames\",text:\"Each connected region compares its cells with the previous frame. The strongest overlap inherits the old identity; unmatched regions receive a new one.\",instruction:\"Press example and play. Outline colors now persist as blobs travel. Splits share a parent identity and merges overlap multiple previous identities.\",code:\"identity = largestOverlap(currentBlob, previousBlobs);\\nif (!identity) identity = nextId++;\",notice:\"Connectivity answers what is one blob now. Temporal matching answers whether it is the same blob as before.\",notYet:\"Fast motion can defeat pure overlap matching. A production tracker would also compare centroid, volume, and velocity.\"},\n blobTracking:{number:\"checkpoint 28\",title:\"A blob is a connected region, not a special particle\",text:\"We threshold the wax fraction, then flood-fill neighboring wax cells. Every connected component receives a label, size, and colored outline.\",instruction:\"Press example and play. Change the connection threshold and watch weak fractional bridges either join regions or separate them. Split and merge counts respond when the component count changes.\",code:\"mask = wax > threshold;\\nblobs = connectedComponents(mask, fourNeighbors);\",notice:\"Blob identity is derived from the wax field. The simulation never stores a blob object; connected shapes emerge from neighboring cell values.\",notYet:\"A changing component count is only an apparent split or merge. Robust event tracking must match component overlap across frames so threshold flicker is not mistaken for topology change.\"}\n };\n\n",[],[149,257,249],"browser","99efe816-cc1f-47eb-826a-14d16620ce76","grid.js","ff57144c-8a98-11f1-9ff2-1607ee4eb77e",256,["D",1785252514857],"export function createGrid(width, height, field) {\n function indexOf(x, y) {\n return y * width + x;\n }\n\n function at(x, y) {\n return field[indexOf(x, y)];\n }\n\n function set(x, y, value) {\n field[indexOf(x, y)] = Math.max(0, Math.min(1, value));\n }\n\n return { indexOf, at, set };\n}\n",[],[149,257,259],"57520181-98db-49c7-9d74-d9f904cff695","config.js","fb6a58ee-8a98-11f1-950d-1607ee4eb77e",255,["D",1785252508284],"export const WIDTH = 28;\nexport const HEIGHT = 18;\nexport const CELL = 24;\n",[],[149,257,267],"f54bf197-8440-41a7-88f5-0079e92e2436","app.js","83f9cfba-8a98-11f1-b875-1607ee4eb77e",257,["D",1785252538414],"import { CELL, HEIGHT, WIDTH } from \"./config.js\";\nimport { createGrid } from \"./grid.js\";\nimport { copy } from \"./lessons.js\";\n\n(function () {\n // Mutable simulation fields and runtime state remain together while the physics\n // functions are migrated behind explicit interfaces in later passes.\n var canvas = document.getElementById(\"field\"),\n ctx = canvas.getContext(\"2d\");\n var field = new Float32Array(WIDTH * HEIGHT),\n next = new Float32Array(WIDTH * HEIGHT);\n var velocityX = new Float32Array(WIDTH * HEIGHT),\n velocityY = new Float32Array(WIDTH * HEIGHT),\n velocityNextX = new Float32Array(WIDTH * HEIGHT),\n velocityNextY = new Float32Array(WIDTH * HEIGHT);\n var divergence = new Float32Array(WIDTH * HEIGHT);\n var pressure = new Float32Array(WIDTH * HEIGHT),\n pressureNext = new Float32Array(WIDTH * HEIGHT);\n var temperature = new Float32Array(WIDTH * HEIGHT),\n temperatureNext = new Float32Array(WIDTH * HEIGHT);\n var wax = new Float32Array(WIDTH * HEIGHT),\n waxNext = new Float32Array(WIDTH * HEIGHT);\n var meltProgress = new Float32Array(WIDTH * HEIGHT),\n meltProgressNext = new Float32Array(WIDTH * HEIGHT);\n var surfaceNormalX = new Float32Array(WIDTH * HEIGHT),\n surfaceNormalY = new Float32Array(WIDTH * HEIGHT),\n surfaceCurvature = new Float32Array(WIDTH * HEIGHT),\n neckMask = new Float32Array(WIDTH * HEIGHT),\n blobLabels = new Int16Array(WIDTH * HEIGHT),\n previousBlobLabels = new Int16Array(WIDTH * HEIGHT),\n vesselMask = new Uint8Array(WIDTH * HEIGHT);\n var mode = \"grid\",\n selected = null,\n painting = false,\n vectorOrigin = null,\n tickCount = 0,\n stepScale = 1,\n running = false,\n animationAccumulator = 0,\n lastAnimationTime = 0,\n cycleHistory = [],\n previousBlobCount = 0,\n splitEvents = 0,\n mergeEvents = 0,\n nextBlobId = 1,\n waxParticles = [];\n var rate = document.getElementById(\"rate\"),\n distance = document.getElementById(\"distance\"),\n iterations = document.getElementById(\"iterations\"),\n buoyancy = document.getElementById(\"buoyancy\"),\n heater = document.getElementById(\"heater\"),\n cooling = document.getElementById(\"cooling\"),\n conduction = document.getElementById(\"conduction\"),\n wallCooling = document.getElementById(\"wallCooling\"),\n neutralTemp = document.getElementById(\"neutralTemp\"),\n densityResponse = document.getElementById(\"densityResponse\"),\n meltTemp = document.getElementById(\"meltTemp\"),\n transitionWidth = document.getElementById(\"transitionWidth\"),\n solidDrag = document.getElementById(\"solidDrag\"),\n interfaceSharpness = document.getElementById(\"interfaceSharpness\"),\n latentCoupling = document.getElementById(\"latentCoupling\"),\n surfaceTension = document.getElementById(\"surfaceTension\"),\n carrierViscosity = document.getElementById(\"carrierViscosity\"),\n hotViscosity = document.getElementById(\"hotViscosity\"),\n coldViscosity = document.getElementById(\"coldViscosity\"),\n timeStep = document.getElementById(\"timeStep\"),\n substeps = document.getElementById(\"substeps\"),\n playbackRate = document.getElementById(\"playbackRate\"),\n lampPreset = document.getElementById(\"lampPreset\"),\n heaterUnevenness = document.getElementById(\"heaterUnevenness\"),\n neckSensitivity = document.getElementById(\"neckSensitivity\"),\n pinchStrength = document.getElementById(\"pinchStrength\"),\n conservationStrength = document.getElementById(\"conservationStrength\"),\n blobThreshold = document.getElementById(\"blobThreshold\"),\n matchRadius = document.getElementById(\"matchRadius\"),\n cellSizeCm = document.getElementById(\"cellSizeCm\"),\n secondsPerFrame = document.getElementById(\"secondsPerFrame\"),\n baseFocus = document.getElementById(\"baseFocus\"),\n glassCooling = document.getElementById(\"glassCooling\"),\n collectionStrength = document.getElementById(\"collectionStrength\"),\n particleCount = document.getElementById(\"particleCount\"),\n restDensity = document.getElementById(\"restDensity\"),\n volumePressure = document.getElementById(\"volumePressure\");\n\n const { indexOf, at, set } = createGrid(WIDTH, HEIGHT, field);\n\n // FIELD SETUP AND LESSON SEEDING\n function clearCurrent() {\n if (mode === \"shapedLamp\" || mode === \"particleVolume\") waxParticles.length = 0;\n if (\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n )\n setRunning(false);\n if (\n mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ) {\n temperature.fill(0);\n temperatureNext.fill(0);\n wax.fill(0);\n waxNext.fill(0);\n meltProgress.fill(0);\n meltProgressNext.fill(0);\n velocityX.fill(0);\n velocityY.fill(0);\n divergence.fill(0);\n pressure.fill(0);\n pressureNext.fill(0);\n } else if (\n mode === \"velocity\" ||\n mode === \"divergence\" ||\n mode === \"projection\" ||\n mode === \"boundaries\"\n ) {\n velocityX.fill(0);\n velocityY.fill(0);\n divergence.fill(0);\n pressure.fill(0);\n pressureNext.fill(0);\n } else if (mode === \"advection\" || mode === \"transportDiffusion\") {\n field.fill(0);\n velocityX.fill(0);\n velocityY.fill(0);\n } else {\n field.fill(0);\n }\n tickCount = 0;\n selected = null;\n draw();\n updateReadout();\n }\n function seedCurrent() {\n clearCurrent();\n if (mode === \"shapedLamp\" || mode === \"particleVolume\") {\n buildLampMask();\n seedWaxParticles();\n } else if (\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ) {\n for (var y = HEIGHT - 4; y < HEIGHT - 1; y++)\n for (var x = 3; x < WIDTH - 3; x++) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.25;\n meltProgress[indexOf(x, y)] = 0;\n }\n } else if (mode === \"stableLamp\" || mode === \"simulator\") {\n for (var y = HEIGHT - 8; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 5; x <= Math.floor(WIDTH / 2) + 4; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 5);\n if (dx * dx + dy * dy < 22) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.62;\n meltProgress[indexOf(x, y)] = 0.65;\n }\n }\n } else if (mode === \"viscosityWax\") {\n for (var y = HEIGHT - 8; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 5; x <= Math.floor(WIDTH / 2) + 4; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 5);\n if (dx * dx + dy * dy < 22) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.68;\n meltProgress[indexOf(x, y)] = 0.75;\n }\n }\n } else if (mode === \"surfaceWax\") {\n for (var y = HEIGHT - 8; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 5; x <= Math.floor(WIDTH / 2) + 4; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 5);\n if (dx * dx + dy * dy < 22) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.7;\n meltProgress[indexOf(x, y)] = 0.8;\n }\n }\n } else if (mode === \"latentWax\") {\n for (var y = HEIGHT - 7; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 4; x <= Math.floor(WIDTH / 2) + 3; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 4);\n if (dx * dx + dy * dy < 16) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = +meltTemp.value - 0.12;\n }\n }\n } else if (mode === \"interfaceWax\") {\n for (var y = HEIGHT - 7; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 4; x <= Math.floor(WIDTH / 2) + 3; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 4);\n if (dx * dx + dy * dy < 16) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.65;\n }\n }\n } else if (mode === \"phaseWax\") {\n for (var y = HEIGHT - 7; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 4; x <= Math.floor(WIDTH / 2) + 3; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 4);\n if (dx * dx + dy * dy < 16) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.35;\n }\n }\n } else if (mode === \"densityWax\") {\n for (var y = HEIGHT - 6; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 5; x < Math.floor(WIDTH / 2); x++) {\n var dx = x - (WIDTH / 2 - 3),\n dy = y - (HEIGHT - 4);\n if (dx * dx + dy * dy < 10) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.85;\n }\n }\n for (var y = 2; y < 6; y++)\n for (var x = Math.floor(WIDTH / 2) + 1; x <= Math.floor(WIDTH / 2) + 5; x++) {\n var dx = x - (WIDTH / 2 + 3),\n dy = y - 4;\n if (dx * dx + dy * dy < 10) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.2;\n }\n }\n } else if (mode === \"waxField\") {\n for (var y = HEIGHT - 6; y < HEIGHT - 2; y++)\n for (var x = Math.floor(WIDTH / 2) - 4; x <= Math.floor(WIDTH / 2) + 3; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 4);\n if (dx * dx + dy * dy < 14) {\n wax[indexOf(x, y)] = 1;\n temperature[indexOf(x, y)] = 0.8;\n }\n }\n } else if (\n mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\"\n ) {\n for (var y = HEIGHT - 5; y < HEIGHT - 1; y++)\n for (var x = WIDTH / 2 - 4; x < WIDTH / 2 + 4; x++) {\n var dx = x - (WIDTH - 1) / 2,\n dy = y - (HEIGHT - 3);\n if (dx * dx + dy * dy < 13) temperature[indexOf(x, y)] = 1;\n }\n } else if (mode === \"heatedFlow\") {\n applyHeatSourceAndCooling();\n } else if (mode === \"boundaries\") {\n seedBoundaryExample();\n updateBoundaryMetric(null, countBoundaryLeaks());\n } else if (mode === \"divergence\" || mode === \"projection\") {\n seedDivergenceExample();\n pressure.fill(0);\n pressureNext.fill(0);\n if (mode === \"projection\") updateProjectionMetric(null, meanAbsoluteDivergence());\n } else if (mode === \"velocity\" || mode === \"advection\" || mode === \"transportDiffusion\") {\n var cx = (WIDTH - 1) / 2,\n cy = (HEIGHT - 1) / 2;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var dx = x - cx,\n dy = y - cy,\n length = Math.sqrt(dx * dx + dy * dy) || 1;\n velocityX[indexOf(x, y)] = -dy / length;\n velocityY[indexOf(x, y)] = dx / length;\n }\n if (mode === \"advection\" || mode === \"transportDiffusion\") {\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var dx = x - 19,\n dy = y - 9;\n if (dx * dx + dy * dy < 7) set(x, y, 1);\n }\n }\n } else {\n for (var y = 6; y < 12; y++)\n for (var x = 10; x < 18; x++) {\n var dx = x - 13.5,\n dy = y - 8.5;\n if (dx * dx + dy * dy < 12) set(x, y, 1);\n }\n }\n draw();\n updateReadout();\n }\n\n function seedBoundaryExample() {\n var cx = (WIDTH - 1) / 2,\n cy = (HEIGHT - 1) / 2;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var dx = x - cx,\n dy = y - cy,\n length = Math.sqrt(dx * dx + dy * dy) || 1;\n var radialX = dx / length,\n radialY = dy / length;\n var vx = 0.7 * radialX - 0.7 * radialY,\n vy = 0.7 * radialY + 0.7 * radialX;\n var speed = Math.sqrt(vx * vx + vy * vy);\n if (speed > 1) {\n vx /= speed;\n vy /= speed;\n }\n velocityX[indexOf(x, y)] = vx;\n velocityY[indexOf(x, y)] = vy;\n }\n }\n function isBoundaryLeak(x, y) {\n var i = indexOf(x, y);\n return (\n (x === 0 && velocityX[i] < 0) ||\n (x === WIDTH - 1 && velocityX[i] > 0) ||\n (y === 0 && velocityY[i] < 0) ||\n (y === HEIGHT - 1 && velocityY[i] > 0)\n );\n }\n function countBoundaryLeaks() {\n var count = 0;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) if (isBoundaryLeak(x, y)) count++;\n return count;\n }\n function updateBoundaryMetric(before, after) {\n document.getElementById(\"boundaryMetric\").textContent =\n before === null ? after : before + \" → \" + after;\n }\n function enforceBoundaries(commitUI) {\n var before = countBoundaryLeaks();\n for (var y = 0; y < HEIGHT; y++) {\n velocityX[indexOf(0, y)] = 0;\n velocityX[indexOf(WIDTH - 1, y)] = 0;\n }\n for (var x = 0; x < WIDTH; x++) {\n velocityY[indexOf(x, 0)] = 0;\n velocityY[indexOf(x, HEIGHT - 1)] = 0;\n }\n var after = countBoundaryLeaks();\n updateBoundaryMetric(before, after);\n if (commitUI !== false) {\n tickCount++;\n draw();\n updateReadout();\n }\n }\n\n function applyBuoyancy(commitUI) {\n var strength = +buoyancy.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n velocityY[i] -= temperature[i] * strength;\n velocityY[i] = Math.max(-1.5, velocityY[i]);\n }\n if (commitUI !== false) {\n tickCount++;\n draw();\n updateReadout();\n }\n }\n\n function sampleTemperature(x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n var x0 = Math.floor(x),\n y0 = Math.floor(y),\n x1 = Math.min(WIDTH - 1, x0 + 1),\n y1 = Math.min(HEIGHT - 1, y0 + 1);\n var tx = x - x0,\n ty = y - y0;\n var top = temperature[indexOf(x0, y0)] * (1 - tx) + temperature[indexOf(x1, y0)] * tx;\n var bottom = temperature[indexOf(x0, y1)] * (1 - tx) + temperature[indexOf(x1, y1)] * tx;\n return top * (1 - ty) + bottom * ty;\n }\n\n function advectTemperature() {\n var stepDistance = +distance.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n temperatureNext[i] = sampleTemperature(\n x - velocityX[i] * stepDistance,\n y - velocityY[i] * stepDistance,\n );\n }\n var swap = temperature;\n temperature = temperatureNext;\n temperatureNext = swap;\n }\n\n function stepRisingHeat() {\n applyBuoyancy(false);\n advectTemperature();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function stepCoupledThermal() {\n applyBuoyancy(false);\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function applyHeatSourceAndCooling() {\n var loss = 1 - +cooling.value,\n gain = +heater.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= loss;\n for (var y = HEIGHT - 3; y < HEIGHT - 1; y++)\n for (var x = Math.floor(WIDTH / 2) - 3; x <= Math.floor(WIDTH / 2) + 2; x++) {\n var i = indexOf(x, y);\n temperature[i] = Math.min(1, temperature[i] + gain);\n }\n }\n\n function stepHeatedFlow() {\n applyHeatSourceAndCooling();\n applyBuoyancy(false);\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n tickCount++;\n draw();\n updateReadout();\n }\n\n // THERMAL AND MATERIAL PHYSICS\n\n function addBottomHeat() {\n var gain = +heater.value * stepScale;\n for (var y = HEIGHT - 3; y < HEIGHT - 1; y++)\n for (var x = Math.floor(WIDTH / 2) - 3; x <= Math.floor(WIDTH / 2) + 2; x++) {\n var i = indexOf(x, y);\n temperature[i] = Math.min(1, temperature[i] + gain);\n }\n }\n\n function diffuseTemperature() {\n var strength = +conduction.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var center = temperature[indexOf(x, y)],\n sum = 0,\n count = 0;\n if (x > 0) {\n sum += temperature[indexOf(x - 1, y)];\n count++;\n }\n if (x < WIDTH - 1) {\n sum += temperature[indexOf(x + 1, y)];\n count++;\n }\n if (y > 0) {\n sum += temperature[indexOf(x, y - 1)];\n count++;\n }\n if (y < HEIGHT - 1) {\n sum += temperature[indexOf(x, y + 1)];\n count++;\n }\n temperatureNext[indexOf(x, y)] = center + strength * (sum / count - center);\n }\n var swap = temperature;\n temperature = temperatureNext;\n temperatureNext = swap;\n }\n\n function coolTemperatureAtWalls() {\n var keep = 1 - +wallCooling.value * stepScale;\n for (var x = 0; x < WIDTH; x++) {\n temperature[indexOf(x, 0)] *= keep;\n temperature[indexOf(x, HEIGHT - 1)] *= keep;\n }\n for (var y = 1; y < HEIGHT - 1; y++) {\n temperature[indexOf(0, y)] *= keep;\n temperature[indexOf(WIDTH - 1, y)] *= keep;\n }\n }\n\n function stepConductiveFlow() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyBuoyancy(false);\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function sampleWax(x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n var x0 = Math.floor(x),\n y0 = Math.floor(y),\n x1 = Math.min(WIDTH - 1, x0 + 1),\n y1 = Math.min(HEIGHT - 1, y0 + 1);\n var tx = x - x0,\n ty = y - y0;\n var top = wax[indexOf(x0, y0)] * (1 - tx) + wax[indexOf(x1, y0)] * tx;\n var bottom = wax[indexOf(x0, y1)] * (1 - tx) + wax[indexOf(x1, y1)] * tx;\n return top * (1 - ty) + bottom * ty;\n }\n\n function advectWax() {\n var stepDistance = +distance.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n waxNext[i] = sampleWax(x - velocityX[i] * stepDistance, y - velocityY[i] * stepDistance);\n }\n var swap = wax;\n wax = waxNext;\n waxNext = swap;\n }\n\n function stepWaxTransport() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyBuoyancy(false);\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function applyWaxDensityBuoyancy() {\n var neutral = +neutralTemp.value,\n response = +densityResponse.value * stepScale;\n for (var i = 0; i < wax.length; i++) {\n var densityAnomaly = wax[i] * (neutral - temperature[i]);\n velocityY[i] += densityAnomaly * response;\n velocityY[i] = Math.max(-1.5, Math.min(1.5, velocityY[i]));\n }\n }\n\n function stepDensityWax() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyWaxDensityBuoyancy();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function liquidWaxFraction(tempValue) {\n var half = Math.max(0.001, +transitionWidth.value),\n start = +meltTemp.value - half,\n end = +meltTemp.value + half;\n var t = Math.max(0, Math.min(1, (tempValue - start) / (end - start)));\n return t * t * (3 - 2 * t);\n }\n\n function applyPhaseDrag() {\n var drag = +solidDrag.value * stepScale;\n for (var i = 0; i < wax.length; i++) {\n var solidFraction = wax[i] * (1 - liquidWaxFraction(temperature[i]));\n var keep = 1 - solidFraction * drag;\n velocityX[i] *= keep;\n velocityY[i] *= keep;\n }\n }\n\n function stepPhaseWax() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyWaxDensityBuoyancy();\n applyPhaseDrag();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function sharpenWaxInterface() {\n var strength = +interfaceSharpness.value;\n for (var i = 0; i < wax.length; i++) {\n var w = Math.max(0, Math.min(1, wax[i]));\n var sharpened = w * w * (3 - 2 * w);\n wax[i] = w * (1 - strength) + sharpened * strength;\n }\n }\n\n function stepInterfaceWax() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyWaxDensityBuoyancy();\n applyPhaseDrag();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n sharpenWaxInterface();\n tickCount++;\n draw();\n updateReadout();\n }\n\n function applyLatentHeat() {\n var melt = +meltTemp.value,\n coupling = +latentCoupling.value * stepScale;\n for (var i = 0; i < wax.length; i++) {\n if (wax[i] <= 0) continue;\n if (temperature[i] > melt && meltProgress[i] < wax[i]) {\n var absorbed = Math.min((temperature[i] - melt) * coupling, wax[i] - meltProgress[i]);\n meltProgress[i] += absorbed;\n temperature[i] -= absorbed;\n } else if (temperature[i] < melt && meltProgress[i] > 0) {\n var released = Math.min((melt - temperature[i]) * coupling, meltProgress[i]);\n meltProgress[i] -= released;\n temperature[i] += released;\n }\n }\n }\n\n function applyLatentPhaseDrag() {\n var drag = +solidDrag.value * stepScale;\n for (var i = 0; i < wax.length; i++) {\n var liquid = wax[i] > 0 ? Math.max(0, Math.min(1, meltProgress[i] / wax[i])) : 0;\n var keep = 1 - wax[i] * (1 - liquid) * drag;\n velocityX[i] *= keep;\n velocityY[i] *= keep;\n }\n }\n\n function sampleMeltProgress(x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n var x0 = Math.floor(x),\n y0 = Math.floor(y),\n x1 = Math.min(WIDTH - 1, x0 + 1),\n y1 = Math.min(HEIGHT - 1, y0 + 1),\n tx = x - x0,\n ty = y - y0;\n var top = meltProgress[indexOf(x0, y0)] * (1 - tx) + meltProgress[indexOf(x1, y0)] * tx;\n var bottom = meltProgress[indexOf(x0, y1)] * (1 - tx) + meltProgress[indexOf(x1, y1)] * tx;\n return top * (1 - ty) + bottom * ty;\n }\n\n function advectMeltProgress() {\n var stepDistance = +distance.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n meltProgressNext[i] = sampleMeltProgress(\n x - velocityX[i] * stepDistance,\n y - velocityY[i] * stepDistance,\n );\n }\n var swap = meltProgress;\n meltProgress = meltProgressNext;\n meltProgressNext = swap;\n }\n\n function stepLatentWax() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n sharpenWaxInterface();\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n tickCount++;\n draw();\n updateReadout();\n }\n\n function waxAt(x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n return wax[indexOf(x, y)];\n }\n function normalAt(array, x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n return array[indexOf(x, y)];\n }\n\n function applySurfaceTension() {\n var sigma = +surfaceTension.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var gx = 0.5 * (waxAt(x + 1, y) - waxAt(x - 1, y)),\n gy = 0.5 * (waxAt(x, y + 1) - waxAt(x, y - 1)),\n length = Math.sqrt(gx * gx + gy * gy);\n var i = indexOf(x, y);\n surfaceNormalX[i] = length > 0.0001 ? gx / length : 0;\n surfaceNormalY[i] = length > 0.0001 ? gy / length : 0;\n }\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y),\n gx = 0.5 * (waxAt(x + 1, y) - waxAt(x - 1, y)),\n gy = 0.5 * (waxAt(x, y + 1) - waxAt(x, y - 1));\n var curvature =\n -0.5 *\n (normalAt(surfaceNormalX, x + 1, y) -\n normalAt(surfaceNormalX, x - 1, y) +\n (normalAt(surfaceNormalY, x, y + 1) - normalAt(surfaceNormalY, x, y - 1)));\n surfaceCurvature[i] = curvature;\n velocityX[i] += sigma * curvature * gx;\n velocityY[i] += sigma * curvature * gy;\n }\n }\n\n function stepSurfaceWax() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applyBottomCollection();\n applySurfaceTension();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n sharpenWaxInterface();\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n tickCount++;\n draw();\n updateReadout();\n }\n\n function viscosityAtIndex(i) {\n var w = wax[i],\n liquid = w > 0 ? Math.max(0, Math.min(1, meltProgress[i] / w)) : 0;\n var waxViscosity = +coldViscosity.value * (1 - liquid) + +hotViscosity.value * liquid;\n return +carrierViscosity.value * (1 - w) + waxViscosity * w;\n }\n\n function diffuseVelocityByViscosity() {\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y),\n sumX = 0,\n sumY = 0,\n count = 0;\n if (x > 0) {\n sumX += velocityX[indexOf(x - 1, y)];\n sumY += velocityY[indexOf(x - 1, y)];\n count++;\n }\n if (x < WIDTH - 1) {\n sumX += velocityX[indexOf(x + 1, y)];\n sumY += velocityY[indexOf(x + 1, y)];\n count++;\n }\n if (y > 0) {\n sumX += velocityX[indexOf(x, y - 1)];\n sumY += velocityY[indexOf(x, y - 1)];\n count++;\n }\n if (y < HEIGHT - 1) {\n sumX += velocityX[indexOf(x, y + 1)];\n sumY += velocityY[indexOf(x, y + 1)];\n count++;\n }\n var blend = viscosityAtIndex(i) * stepScale;\n velocityNextX[i] = velocityX[i] + blend * (sumX / count - velocityX[i]);\n velocityNextY[i] = velocityY[i] + blend * (sumY / count - velocityY[i]);\n }\n var swap = velocityX;\n velocityX = velocityNextX;\n velocityNextX = swap;\n swap = velocityY;\n velocityY = velocityNextY;\n velocityNextY = swap;\n }\n\n function stepViscosityWax() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applySurfaceTension();\n diffuseVelocityByViscosity();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n sharpenWaxInterface();\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n tickCount++;\n draw();\n updateReadout();\n }\n\n function runCompleteSubstep() {\n addBottomHeat();\n var ambientKeep = 1 - +cooling.value * stepScale;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applySurfaceTension();\n diffuseVelocityByViscosity();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n sharpenWaxInterface();\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n }\n\n function updateStabilityMetric() {\n var maxSpeed = 0;\n for (var i = 0; i < velocityX.length; i++)\n maxSpeed = Math.max(\n maxSpeed,\n Math.sqrt(velocityX[i] * velocityX[i] + velocityY[i] * velocityY[i]),\n );\n document.getElementById(\"stabilityMetric\").textContent =\n (maxSpeed * +distance.value * stepScale).toFixed(2) + \" cells\";\n }\n\n function stepStableLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runCompleteSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n draw();\n updateReadout();\n }\n\n function addUnevenBottomHeat() {\n var gain = +heater.value * stepScale,\n uneven = +heaterUnevenness.value;\n for (var y = HEIGHT - 3; y < HEIGHT - 1; y++)\n for (var x = 3; x < WIDTH - 3; x++) {\n var wave =\n 0.5 +\n 0.5 * Math.sin(x * 1.37 + tickCount * 0.19) +\n 0.25 * Math.sin(x * 0.53 - tickCount * 0.11);\n wave = Math.max(0, Math.min(1, wave));\n var localGain = gain * (1 - uneven + uneven * (0.35 + 1.3 * wave));\n var i = indexOf(x, y);\n temperature[i] = Math.min(1, temperature[i] + localGain);\n }\n }\n\n function runUnevenSubstep() {\n addUnevenBottomHeat();\n var ambientKeep = 1 - +cooling.value * stepScale;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applySurfaceTension();\n diffuseVelocityByViscosity();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n sharpenWaxInterface();\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n }\n\n function totalWax() {\n var total = 0;\n for (var i = 0; i < wax.length; i++) total += wax[i];\n return total;\n }\n function detectAndPinchNecks() {\n neckMask.fill(0);\n var before = totalWax(),\n removed = 0,\n count = 0,\n sensitivity = +neckSensitivity.value,\n strength = +pinchStrength.value;\n waxNext.set(wax);\n for (var y = 1; y < HEIGHT - 3; y++)\n for (var x = 1; x < WIDTH - 1; x++) {\n var i = indexOf(x, y),\n w = wax[i];\n if (w < 0.08 || w > 0.92) continue;\n var left = wax[indexOf(x - 1, y)],\n right = wax[indexOf(x + 1, y)],\n up = wax[indexOf(x, y - 1)],\n down = wax[indexOf(x, y + 1)];\n var vertical = up + down,\n horizontal = left + right;\n var verticalBridge = vertical > 1.05 && horizontal < 1.05 - 0.7 * sensitivity;\n var horizontalBridge = horizontal > 1.05 && vertical < 1.05 - 0.7 * sensitivity;\n var strain =\n Math.abs(velocityY[indexOf(x, y - 1)] - velocityY[indexOf(x, y + 1)]) +\n Math.abs(velocityX[indexOf(x - 1, y)] - velocityX[indexOf(x + 1, y)]);\n if (!(verticalBridge || horizontalBridge) || strain < 0.04 * (1 - sensitivity)) continue;\n var cut = Math.min(w, 0.18 * strength * stepScale * (0.6 + Math.min(2, strain)));\n if (cut <= 0) continue;\n neckMask[i] = 1;\n waxNext[i] -= cut;\n removed += cut;\n count++;\n var a = verticalBridge ? indexOf(x, y - 1) : indexOf(x - 1, y),\n b = verticalBridge ? indexOf(x, y + 1) : indexOf(x + 1, y);\n var capacityA = 1 - waxNext[a],\n capacityB = 1 - waxNext[b],\n capacity = capacityA + capacityB;\n if (capacity > 0) {\n var putA = Math.min(capacityA, (cut * capacityA) / capacity),\n putB = Math.min(capacityB, cut - putA);\n waxNext[a] += putA;\n waxNext[b] += putB;\n }\n }\n wax.set(waxNext);\n var after = totalWax(),\n error = before - after;\n if (Math.abs(error) > 1e-7) {\n for (var pass = 0; pass < 3 && Math.abs(error) > 1e-6; pass++) {\n var capacity = 0;\n for (var i = 0; i < wax.length; i++)\n capacity += error > 0 ? (wax[i] > 0.05 ? 1 - wax[i] : 0) : wax[i];\n if (capacity <= 0) break;\n for (var i = 0; i < wax.length; i++) {\n var room = error > 0 ? (wax[i] > 0.05 ? 1 - wax[i] : 0) : wax[i];\n wax[i] = Math.max(0, Math.min(1, wax[i] + (error * room) / capacity));\n }\n error = before - totalWax();\n }\n }\n document.getElementById(\"neckMetric\").textContent = count;\n document.getElementById(\"waxErrorMetric\").textContent = Math.abs(before - totalWax()).toFixed(\n 3,\n );\n }\n\n function runPinchSubstep() {\n addUnevenBottomHeat();\n var ambientKeep = 1 - +cooling.value * stepScale;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applySurfaceTension();\n diffuseVelocityByViscosity();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n detectAndPinchNecks();\n sharpenWaxInterface();\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n }\n\n function correctWaxMass(target, strength) {\n var raw = totalWax() - target,\n remaining = -raw * strength;\n for (var pass = 0; pass < 4 && Math.abs(remaining) > 1e-7; pass++) {\n var capacity = 0;\n for (var i = 0; i < wax.length; i++)\n capacity += remaining > 0 ? (wax[i] > 0.02 ? 1 - wax[i] : 0) : wax[i];\n if (capacity <= 0) break;\n for (var i = 0; i < wax.length; i++) {\n var room = remaining > 0 ? (wax[i] > 0.02 ? 1 - wax[i] : 0) : wax[i];\n wax[i] = Math.max(0, Math.min(1, wax[i] + (remaining * room) / capacity));\n }\n remaining = target - totalWax();\n }\n document.getElementById(\"rawDriftMetric\").textContent = (raw >= 0 ? \"+\" : \"\") + raw.toFixed(3);\n var corrected = totalWax() - target;\n document.getElementById(\"correctedDriftMetric\").textContent =\n (corrected >= 0 ? \"+\" : \"\") + corrected.toFixed(3);\n }\n\n function totalThermalEnergy() {\n var total = 0;\n for (var i = 0; i < temperature.length; i++) total += temperature[i] + meltProgress[i];\n return total;\n }\n function analyzeBlobs(countEvents, persistent) {\n blobLabels.fill(-1);\n var threshold = +blobThreshold.value,\n components = [],\n queue = [];\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var st = indexOf(x, y);\n if (wax[st] < threshold || blobLabels[st] >= 0) continue;\n var cells = [],\n ci = components.length;\n blobLabels[st] = ci;\n queue = [st];\n for (var q = 0; q < queue.length; q++) {\n var i = queue[q],\n cx = i % WIDTH,\n cy = Math.floor(i / WIDTH);\n cells.push(i);\n var ns = [\n [cx - 1, cy],\n [cx + 1, cy],\n [cx, cy - 1],\n [cx, cy + 1],\n ];\n for (var n = 0; n < 4; n++) {\n var nx = ns[n][0],\n ny = ns[n][1];\n if (nx < 0 || nx >= WIDTH || ny < 0 || ny >= HEIGHT) continue;\n var ni = indexOf(nx, ny);\n if (wax[ni] >= threshold && blobLabels[ni] < 0) {\n blobLabels[ni] = ci;\n queue.push(ni);\n }\n }\n }\n components.push(cells);\n }\n var ids = [],\n matched = 0,\n parents = {};\n for (var c = 0; c < components.length; c++) {\n var votes = {};\n for (var k = 0; k < components[c].length; k++) {\n var old = previousBlobLabels[components[c][k]];\n if (old >= 0) votes[old] = (votes[old] || 0) + 1;\n else if (persistent) {\n var rr = +matchRadius.value,\n cx = components[c][k] % WIDTH,\n cy = Math.floor(components[c][k] / WIDTH);\n for (var sy = Math.max(0, cy - rr); sy <= Math.min(HEIGHT - 1, cy + rr); sy++)\n for (var sx = Math.max(0, cx - rr); sx <= Math.min(WIDTH - 1, cx + rr); sx++) {\n var nearby = previousBlobLabels[indexOf(sx, sy)];\n if (nearby >= 0) votes[nearby] = (votes[nearby] || 0) + 0.15;\n }\n }\n }\n var best = -1,\n bestN = 0,\n keys = Object.keys(votes);\n for (var j = 0; j < keys.length; j++) {\n var id = +keys[j];\n if (votes[id] > bestN) {\n best = id;\n bestN = votes[id];\n }\n }\n if (persistent && best >= 0) {\n ids[c] = best;\n matched++;\n parents[best] = (parents[best] || 0) + 1;\n if (countEvents && keys.length > 1) mergeEvents++;\n } else ids[c] = nextBlobId++;\n }\n if (countEvents) for (var id in parents) if (parents[id] > 1) splitEvents += parents[id] - 1;\n var total = totalWax(),\n largest = 0;\n for (var c = 0; c < components.length; c++) {\n var size = 0;\n for (var k = 0; k < components[c].length; k++) {\n var i = components[c][k];\n size += wax[i];\n blobLabels[i] = ids[c];\n }\n largest = Math.max(largest, size);\n }\n previousBlobLabels.set(blobLabels);\n previousBlobCount = components.length;\n document.getElementById(\"blobCountMetric\").textContent = components.length;\n document.getElementById(\"largestBlobMetric\").textContent =\n Math.round((total ? largest / total : 0) * 100) + \"%\";\n document.getElementById(\"splitMetric\").textContent = splitEvents;\n document.getElementById(\"mergeMetric\").textContent = mergeEvents;\n document.getElementById(\"nextBlobMetric\").textContent = nextBlobId;\n document.getElementById(\"matchedBlobMetric\").textContent = matched;\n }\n function stepBlobTrackingLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runClosedVesselSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n recordCycleDiagnostics();\n analyzeBlobs(true);\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function stepBlobIdentityLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runClosedVesselSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n recordCycleDiagnostics();\n analyzeBlobs(true, true);\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function updateCourseDiagnostics() {\n var min = 1,\n max = 0;\n for (var i = 0; i < cycleHistory.length; i++) {\n min = Math.min(min, cycleHistory[i]);\n max = Math.max(max, cycleHistory[i]);\n }\n var range = Math.max(0, max - min),\n blobs = +document.getElementById(\"blobCountMetric\").textContent || 0;\n document.getElementById(\"cycleRangeMetric\").textContent = Math.round(range * 100) + \"%\";\n document.getElementById(\"activeBlobsMetric\").textContent = blobs;\n var maxSpeed = 0;\n for (var i = 0; i < velocityX.length; i++)\n maxSpeed = Math.max(\n maxSpeed,\n Math.sqrt(velocityX[i] * velocityX[i] + velocityY[i] * velocityY[i]),\n );\n document.getElementById(\"physicalSpeedMetric\").textContent =\n ((maxSpeed * +cellSizeCm.value) / +secondsPerFrame.value).toFixed(1) + \" cm/s\";\n document.getElementById(\"finalBlobsMetric\").textContent = blobs;\n document.getElementById(\"finalCycleMetric\").textContent = Math.round(range * 100) + \"%\";\n document.getElementById(\"finalRisingMetric\").textContent =\n document.getElementById(\"risingMetric\").textContent;\n document.getElementById(\"finalFallingMetric\").textContent =\n document.getElementById(\"fallingMetric\").textContent;\n }\n function lampHalfWidth(y) {\n var t = Math.max(0, Math.min(1, y / (HEIGHT - 1)));\n return 4.5 + 6.2 * Math.pow(Math.sin(Math.PI * t), 0.72);\n }\n // SHAPED VESSEL AND PARTICLE PHYSICS\n\n function buildLampMask() {\n var c = (WIDTH - 1) / 2;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++)\n vesselMask[indexOf(x, y)] =\n y > 0 && y < HEIGHT - 1 && Math.abs(x - c) <= lampHalfWidth(y) ? 1 : 0;\n }\n function seedWaxParticles() {\n waxParticles.length = 0;\n var c = (WIDTH - 1) / 2,\n target = +particleCount.value;\n for (var y = HEIGHT - 7.4; y <= HEIGHT - 2.1 && waxParticles.length < target; y += 0.52)\n for (var x = c - 5.6; x <= c + 5.6 && waxParticles.length < target; x += 0.52) {\n if (Math.abs(x - c) < lampHalfWidth(y) - 0.5) {\n var n = waxParticles.length;\n waxParticles.push({\n x: x + 0.06 * Math.sin(n * 1.7),\n y: y + 0.04 * Math.cos(n * 2.1),\n vx: 0,\n vy: 0,\n temp: 0.2 + 0.035 * Math.sin(x * 1.9),\n });\n }\n }\n rasterizeWaxParticles();\n }\n function rasterizeWaxParticles() {\n wax.fill(0);\n temperature.fill(0);\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n x = Math.max(0, Math.min(WIDTH - 1, Math.round(p.x))),\n y = Math.max(0, Math.min(HEIGHT - 1, Math.round(p.y))),\n i = indexOf(x, y);\n wax[i] = Math.min(1, wax[i] + 0.32);\n temperature[i] = Math.max(temperature[i], p.temp);\n }\n }\n function drawWaxParticles() {\n ctx.save();\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n hot = Math.max(0, Math.min(1, p.temp));\n ctx.beginPath();\n ctx.arc((p.x + 0.5) * CELL, (p.y + 0.5) * CELL, CELL * 0.42, 0, Math.PI * 2);\n ctx.fillStyle =\n \"rgba(\" +\n Math.round(155 + 90 * hot) +\n \",\" +\n Math.round(70 + 20 * hot) +\n \",\" +\n Math.round(205 - 45 * hot) +\n \",.42)\";\n ctx.fill();\n }\n ctx.restore();\n }\n function collideParticle(p) {\n var c = (WIDTH - 1) / 2,\n minY = 1.35,\n maxY = HEIGHT - 2.05;\n p.y = Math.max(minY, Math.min(maxY, p.y));\n if (p.y === minY && p.vy < 0) p.vy *= -0.18;\n if (p.y === maxY && p.vy > 0) p.vy *= -0.12;\n var limit = lampHalfWidth(p.y) - 0.48;\n if (p.x < c - limit) {\n p.x = c - limit;\n p.vx = Math.abs(p.vx) * 0.2;\n }\n if (p.x > c + limit) {\n p.x = c + limit;\n p.vx = -Math.abs(p.vx) * 0.2;\n }\n }\n function updateParticleThermal() {\n var c = (WIDTH - 1) / 2,\n width = 4 - 3 * +baseFocus.value,\n dt = stepScale,\n hotspot = 1.8 * Math.sin(tickCount * 0.045);\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n d = (p.x - (c + hotspot)) / width,\n uneven = 0.65 + 0.45 * (0.5 + 0.5 * Math.sin(p.x * 2.1 + n * 0.17));\n if (p.y > HEIGHT - 5.5) p.temp += +heater.value * dt * Math.exp(-d * d) * uneven;\n p.temp *= 1 - +cooling.value * dt;\n var wallGap = lampHalfWidth(p.y) - Math.abs(p.x - c);\n if (p.y < 4 || wallGap < 1.2) p.temp *= 1 - +glassCooling.value * dt;\n p.temp = Math.max(0, Math.min(1, p.temp));\n }\n }\n function applyParticleForces() {\n var dt = stepScale,\n c = +collectionStrength.value,\n neutral = +neutralTemp.value,\n response = +densityResponse.value;\n for (var a = 0; a < waxParticles.length; a++) {\n var p = waxParticles[a];\n p.vy += (neutral - p.temp) * response * 0.24 * dt;\n p.vx *= 0.992;\n p.vy *= 0.992;\n }\n for (var a = 0; a < waxParticles.length; a++)\n for (var q = a + 1; q < waxParticles.length; q++) {\n var p = waxParticles[a],\n o = waxParticles[q],\n dx = o.x - p.x,\n dy = o.y - p.y,\n d2 = dx * dx + dy * dy;\n if (d2 < 0.01 || d2 > 2.25) continue;\n var d = Math.sqrt(d2),\n phaseCohesion = 0.22 + 0.78 * (1 - (p.temp + o.temp) * 0.5),\n force = (d - 0.58) * 0.032 * c * phaseCohesion * dt,\n fx = (force * dx) / d,\n fy = (force * dy) / d;\n p.vx += fx;\n p.vy += fy;\n o.vx -= fx;\n o.vy -= fy;\n }\n }\n function applyParticleVolumePressure() {\n var count = waxParticles.length,\n h = 1.3,\n density = new Float32Array(count),\n maxD = 0,\n target = +restDensity.value,\n strength = +volumePressure.value * stepScale;\n for (var a = 0; a < count; a++) {\n density[a] = 1;\n for (var q = a + 1; q < count; q++) {\n var p = waxParticles[a],\n o = waxParticles[q],\n dx = o.x - p.x,\n dy = o.y - p.y,\n d = Math.sqrt(dx * dx + dy * dy);\n if (d < h) {\n var w = 1 - d / h;\n w *= w;\n density[a] += w;\n density[q] += w;\n }\n }\n }\n for (var a = 0; a < count; a++) maxD = Math.max(maxD, density[a]);\n for (var a = 0; a < count; a++)\n for (var q = a + 1; q < count; q++) {\n var p = waxParticles[a],\n o = waxParticles[q],\n dx = o.x - p.x,\n dy = o.y - p.y,\n d = Math.sqrt(dx * dx + dy * dy);\n if (d < 0.001 || d >= h) continue;\n var excess = Math.max(0, (density[a] + density[q]) * 0.5 - target),\n hard = Math.max(0, 0.48 - d) * 2,\n force = (excess * (1 - d / h) * 0.055 + hard * 0.12) * strength,\n fx = (force * dx) / d,\n fy = (force * dy) / d;\n p.vx -= fx;\n p.vy -= fy;\n o.vx += fx;\n o.vy += fy;\n }\n document.getElementById(\"particleDensityMetric\").textContent = maxD.toFixed(1);\n }\n function depositParticleVelocityBilinear() {\n velocityX.fill(0);\n velocityY.fill(0);\n var weights = new Float32Array(WIDTH * HEIGHT);\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n x0 = Math.floor(p.x),\n y0 = Math.floor(p.y),\n tx = p.x - x0,\n ty = p.y - y0;\n for (var oy = 0; oy <= 1; oy++)\n for (var ox = 0; ox <= 1; ox++) {\n var x = Math.max(0, Math.min(WIDTH - 1, x0 + ox)),\n y = Math.max(0, Math.min(HEIGHT - 1, y0 + oy)),\n w = (ox ? tx : 1 - tx) * (oy ? ty : 1 - ty),\n i = indexOf(x, y);\n velocityX[i] += p.vx * w;\n velocityY[i] += p.vy * w;\n weights[i] += w;\n }\n }\n for (var i = 0; i < weights.length; i++)\n if (weights[i]) {\n velocityX[i] /= weights[i];\n velocityY[i] /= weights[i];\n }\n }\n function sampleProjectedVelocityBilinear(p) {\n var x0 = Math.floor(p.x),\n y0 = Math.floor(p.y),\n tx = p.x - x0,\n ty = p.y - y0,\n vx = 0,\n vy = 0;\n for (var oy = 0; oy <= 1; oy++)\n for (var ox = 0; ox <= 1; ox++) {\n var x = Math.max(0, Math.min(WIDTH - 1, x0 + ox)),\n y = Math.max(0, Math.min(HEIGHT - 1, y0 + oy)),\n w = (ox ? tx : 1 - tx) * (oy ? ty : 1 - ty),\n i = indexOf(x, y);\n vx += velocityX[i] * w;\n vy += velocityY[i] * w;\n }\n return { vx: vx, vy: vy };\n }\n function projectParticleMotionBilinear() {\n depositParticleVelocityBilinear();\n projectVelocity(false);\n enforceBoundaries(false);\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n if (!vesselMask[i]) velocityX[i] = velocityY[i] = 0;\n }\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n v = sampleProjectedVelocityBilinear(p);\n p.vx = p.vx * 0.52 + v.vx * 0.48;\n p.vy = p.vy * 0.52 + v.vy * 0.48;\n }\n }\n function runParticleVolumeSubstep() {\n updateParticleThermal();\n applyParticleVolumePressure();\n applyParticleForces();\n projectParticleMotionBilinear();\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n];\n p.x += p.vx * stepScale;\n p.y += p.vy * stepScale;\n collideParticle(p);\n }\n rasterizeWaxParticles();\n }\n function stepParticleVolume() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runParticleVolumeSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n recordCycleDiagnostics();\n analyzeBlobs(true, true);\n updateCourseDiagnostics();\n document.getElementById(\"shapedWaxMetric\").textContent = waxParticles.length;\n document.getElementById(\"shapedBlobMetric\").textContent =\n document.getElementById(\"blobCountMetric\").textContent;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n function projectParticleMotion() {\n velocityX.fill(0);\n velocityY.fill(0);\n var counts = new Float32Array(WIDTH * HEIGHT);\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n x = Math.round(p.x),\n y = Math.round(p.y),\n i = indexOf(Math.max(0, Math.min(WIDTH - 1, x)), Math.max(0, Math.min(HEIGHT - 1, y)));\n velocityX[i] += p.vx;\n velocityY[i] += p.vy;\n counts[i]++;\n }\n for (var i = 0; i < counts.length; i++)\n if (counts[i]) {\n velocityX[i] /= counts[i];\n velocityY[i] /= counts[i];\n }\n projectVelocity(false);\n enforceBoundaries(false);\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n if (!vesselMask[i]) velocityX[i] = velocityY[i] = 0;\n }\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n],\n i = indexOf(\n Math.max(0, Math.min(WIDTH - 1, Math.round(p.x))),\n Math.max(0, Math.min(HEIGHT - 1, Math.round(p.y))),\n );\n p.vx = p.vx * 0.55 + velocityX[i] * 0.45;\n p.vy = p.vy * 0.55 + velocityY[i] * 0.45;\n }\n }\n function runShapedLampSubstep() {\n updateParticleThermal();\n applyParticleForces();\n projectParticleMotion();\n for (var n = 0; n < waxParticles.length; n++) {\n var p = waxParticles[n];\n p.x += p.vx * stepScale;\n p.y += p.vy * stepScale;\n collideParticle(p);\n }\n rasterizeWaxParticles();\n }\n function stepShapedLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runShapedLampSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n recordCycleDiagnostics();\n analyzeBlobs(true, true);\n updateCourseDiagnostics();\n document.getElementById(\"shapedWaxMetric\").textContent = waxParticles.length;\n document.getElementById(\"shapedBlobMetric\").textContent =\n document.getElementById(\"blobCountMetric\").textContent;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n function stepRemainingLesson() {\n stepBlobIdentityLamp();\n updateCourseDiagnostics();\n }\n\n // DIAGNOSTICS, CONSERVATION, AND RUNTIME\n\n function recordCycleDiagnostics() {\n var total = 0,\n heightSum = 0,\n rising = 0,\n falling = 0;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y),\n w = wax[i];\n total += w;\n heightSum += w * (1 - y / (HEIGHT - 1));\n if (velocityY[i] < -0.03) rising += w;\n else if (velocityY[i] > 0.03) falling += w;\n }\n var height = total > 0 ? heightSum / total : 0;\n cycleHistory.push(height);\n if (cycleHistory.length > 100) cycleHistory.shift();\n document.getElementById(\"waxHeightMetric\").textContent = Math.round(height * 100) + \"%\";\n document.getElementById(\"risingMetric\").textContent =\n Math.round((total ? rising / total : 0) * 100) + \"%\";\n document.getElementById(\"fallingMetric\").textContent =\n Math.round((total ? falling / total : 0) * 100) + \"%\";\n var chart = document.getElementById(\"cycleChart\"),\n c = chart.getContext(\"2d\"),\n w = chart.width,\n h = chart.height;\n c.clearRect(0, 0, w, h);\n c.fillStyle = \"#0b1017\";\n c.fillRect(0, 0, w, h);\n c.strokeStyle = \"#303640\";\n c.lineWidth = 1;\n for (var g = 1; g < 4; g++) {\n var gy = (g * h) / 4;\n c.beginPath();\n c.moveTo(0, gy);\n c.lineTo(w, gy);\n c.stroke();\n }\n if (cycleHistory.length > 1) {\n c.strokeStyle = \"#50a7ff\";\n c.lineWidth = 2;\n c.beginPath();\n for (var n = 0; n < cycleHistory.length; n++) {\n var px = (n * w) / 99,\n py = h - 4 - cycleHistory[n] * (h - 8);\n if (n === 0) c.moveTo(px, py);\n else c.lineTo(px, py);\n }\n c.stroke();\n }\n }\n function stepCycleDiagnosticsLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runClosedVesselSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n recordCycleDiagnostics();\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function runClosedVesselSubstep() {\n var targetWax = totalWax();\n runEnergyBudgetSubstep();\n correctWaxMass(targetWax, 1);\n document.getElementById(\"inventoryMetric\").textContent = totalWax().toFixed(1);\n }\n function stepClosedVesselLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runClosedVesselSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function runEnergyBudgetSubstep() {\n var start = totalThermalEnergy(),\n targetWax = totalWax();\n addUnevenBottomHeat();\n var afterHeater = totalThermalEnergy();\n var ambientKeep = 1 - +cooling.value * stepScale;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n var afterCooling = totalThermalEnergy();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applySurfaceTension();\n diffuseVelocityByViscosity();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n detectAndPinchNecks();\n sharpenWaxInterface();\n correctWaxMass(targetWax, 1);\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n var finish = totalThermalEnergy(),\n energyIn = afterHeater - start,\n energyOut = afterCooling - afterHeater,\n drift = finish - afterCooling;\n var latent = 0;\n for (var i = 0; i < meltProgress.length; i++) latent += meltProgress[i];\n document.getElementById(\"energyInMetric\").textContent = \"+\" + energyIn.toFixed(3);\n document.getElementById(\"energyOutMetric\").textContent = energyOut.toFixed(3);\n document.getElementById(\"latentEnergyMetric\").textContent = latent.toFixed(3);\n document.getElementById(\"energyDriftMetric\").textContent =\n (drift >= 0 ? \"+\" : \"\") + drift.toFixed(3);\n }\n\n function stepEnergyBudgetLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runEnergyBudgetSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function runConservationSubstep() {\n var target = totalWax();\n addUnevenBottomHeat();\n var ambientKeep = 1 - +cooling.value * stepScale;\n for (var i = 0; i < temperature.length; i++) temperature[i] *= ambientKeep;\n diffuseTemperature();\n coolTemperatureAtWalls();\n applyLatentHeat();\n applyWaxDensityBuoyancy();\n applyLatentPhaseDrag();\n applySurfaceTension();\n diffuseVelocityByViscosity();\n projectVelocity(false);\n enforceBoundaries(false);\n advectTemperature();\n advectWax();\n advectMeltProgress();\n detectAndPinchNecks();\n sharpenWaxInterface();\n correctWaxMass(target, +conservationStrength.value);\n for (var i = 0; i < meltProgress.length; i++)\n meltProgress[i] = Math.min(meltProgress[i], wax[i]);\n }\n\n function stepConservationLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runConservationSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function stepPinchLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runPinchSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function stepUnevenLamp() {\n var count = +substeps.value;\n stepScale = +timeStep.value / count;\n for (var n = 0; n < count; n++) runUnevenSubstep();\n updateStabilityMetric();\n stepScale = 1;\n tickCount++;\n draw();\n updateReadout();\n updateSimulatorDiagnostics();\n }\n\n function updateSimulatorDiagnostics() {\n var waxTotal = 0,\n tempTotal = 0,\n meltTotal = 0,\n maxSpeed = 0;\n for (var i = 0; i < wax.length; i++) {\n waxTotal += wax[i];\n tempTotal += temperature[i];\n meltTotal += meltProgress[i];\n maxSpeed = Math.max(\n maxSpeed,\n Math.sqrt(velocityX[i] * velocityX[i] + velocityY[i] * velocityY[i]),\n );\n }\n document.getElementById(\"waxMetric\").textContent = waxTotal.toFixed(1);\n document.getElementById(\"tempMetric\").textContent = (tempTotal / temperature.length).toFixed(2);\n document.getElementById(\"meltMetric\").textContent = (\n waxTotal > 0 ? meltTotal / waxTotal : 0\n ).toFixed(2);\n document.getElementById(\"speedMetric\").textContent = maxSpeed.toFixed(2);\n }\n\n function stepFinalLamp() {\n stepStableLamp();\n updateSimulatorDiagnostics();\n }\n\n function setRunning(value) {\n running = value;\n document.getElementById(\"play\").textContent = running ? \"pause\" : \"play\";\n if (!running) animationAccumulator = 0;\n }\n\n function applyLampPreset() {\n var preset = lampPreset.value;\n if (preset === \"slow\") {\n heater.value = 0.08;\n cooling.value = 0.03;\n densityResponse.value = 0.8;\n surfaceTension.value = 0.5;\n hotViscosity.value = 0.2;\n coldViscosity.value = 0.5;\n timeStep.value = 0.75;\n substeps.value = 6;\n } else if (preset === \"active\") {\n heater.value = 0.2;\n cooling.value = 0.02;\n densityResponse.value = 1.8;\n surfaceTension.value = 0.25;\n hotViscosity.value = 0.05;\n coldViscosity.value = 0.35;\n timeStep.value = 1.25;\n substeps.value = 6;\n } else {\n heater.value = 0.12;\n cooling.value = 0.02;\n densityResponse.value = 1.2;\n surfaceTension.value = 0.35;\n hotViscosity.value = 0.1;\n coldViscosity.value = 0.45;\n timeStep.value = 1;\n substeps.value = 4;\n }\n document.getElementById(\"heaterOut\").textContent = (+heater.value).toFixed(2);\n document.getElementById(\"coolingOut\").textContent = Math.round(+cooling.value * 100) + \"%\";\n document.getElementById(\"densityOut\").textContent = (+densityResponse.value).toFixed(2);\n document.getElementById(\"surfaceOut\").textContent = (+surfaceTension.value).toFixed(2);\n document.getElementById(\"hotViscosityOut\").textContent =\n Math.round(+hotViscosity.value * 100) + \"%\";\n document.getElementById(\"coldViscosityOut\").textContent =\n Math.round(+coldViscosity.value * 100) + \"%\";\n document.getElementById(\"timeStepOut\").textContent = (+timeStep.value).toFixed(2);\n document.getElementById(\"substepsOut\").textContent = substeps.value;\n }\n\n function animationLoop(now) {\n if (!lastAnimationTime) lastAnimationTime = now;\n var elapsed = Math.min(0.25, (now - lastAnimationTime) / 1000);\n lastAnimationTime = now;\n if (\n running &&\n (mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\")\n ) {\n animationAccumulator += elapsed * +playbackRate.value;\n var safety = 0;\n while (animationAccumulator >= 1 && safety < 4) {\n if (mode === \"particleVolume\") stepParticleVolume();\n else if (mode === \"shapedLamp\") stepShapedLamp();\n else if (\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\"\n )\n stepRemainingLesson();\n else if (mode === \"blobIdentity\") stepBlobIdentityLamp();\n else if (mode === \"blobTracking\") stepBlobTrackingLamp();\n else if (mode === \"cycleDiagnostics\") stepCycleDiagnosticsLamp();\n else if (mode === \"closedVessel\") stepClosedVesselLamp();\n else if (mode === \"energyBudget\") stepEnergyBudgetLamp();\n else if (mode === \"conservation\") stepConservationLamp();\n else if (mode === \"pinchOff\") stepPinchLamp();\n else if (mode === \"unevenHeating\") stepUnevenLamp();\n else stepFinalLamp();\n animationAccumulator -= 1;\n safety++;\n }\n }\n requestAnimationFrame(animationLoop);\n }\n\n // FIELD RENDERING, ADVECTION, AND POINTER INPUT\n\n function colorForLatentWax(tempValue, waxValue, progress) {\n var liquid = waxValue > 0 ? Math.max(0, Math.min(1, progress / waxValue)) : 0;\n var phaseR = 150 + 90 * liquid,\n phaseG = 95 - 20 * liquid,\n phaseB = 205 - 55 * liquid;\n var baseR = 11,\n baseG = 16,\n baseB = 23;\n if (tempValue > 0) {\n baseR = 40 + 215 * tempValue;\n baseG = 35 + 95 * tempValue;\n baseB = 50 - 25 * tempValue;\n }\n return (\n \"rgb(\" +\n Math.round(baseR * (1 - waxValue) + phaseR * waxValue) +\n \",\" +\n Math.round(baseG * (1 - waxValue) + phaseG * waxValue) +\n \",\" +\n Math.round(baseB * (1 - waxValue) + phaseB * waxValue) +\n \")\"\n );\n }\n\n function colorForPhaseWax(tempValue, waxValue) {\n var liquid = liquidWaxFraction(tempValue);\n var phaseR = 150 + 90 * liquid,\n phaseG = 95 - 20 * liquid,\n phaseB = 205 - 55 * liquid;\n var baseR = 11,\n baseG = 16,\n baseB = 23;\n if (tempValue > 0) {\n baseR = 40 + 215 * tempValue;\n baseG = 35 + 95 * tempValue;\n baseB = 50 - 25 * tempValue;\n }\n return (\n \"rgb(\" +\n Math.round(baseR * (1 - waxValue) + phaseR * waxValue) +\n \",\" +\n Math.round(baseG * (1 - waxValue) + phaseG * waxValue) +\n \",\" +\n Math.round(baseB * (1 - waxValue) + phaseB * waxValue) +\n \")\"\n );\n }\n\n function colorForWax(tempValue, waxValue) {\n var baseR = 11,\n baseG = 16,\n baseB = 23;\n if (tempValue > 0) {\n baseR = 40 + 215 * tempValue;\n baseG = 35 + 95 * tempValue;\n baseB = 50 - 25 * tempValue;\n }\n var waxR = 224,\n waxG = 80,\n waxB = 170;\n return (\n \"rgb(\" +\n Math.round(baseR * (1 - waxValue) + waxR * waxValue) +\n \",\" +\n Math.round(baseG * (1 - waxValue) + waxG * waxValue) +\n \",\" +\n Math.round(baseB * (1 - waxValue) + waxB * waxValue) +\n \")\"\n );\n }\n\n function colorForTemperature(value) {\n if (value <= 0) return \"#0b1017\";\n var r = Math.round(40 + 215 * value),\n g = Math.round(35 + 95 * value),\n b = Math.round(50 - 25 * value);\n return \"rgb(\" + r + \",\" + g + \",\" + b + \")\";\n }\n\n function seedDivergenceExample() {\n var source = { x: 8, y: 9 },\n sink = { x: 20, y: 9 };\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var sx = x - source.x,\n sy = y - source.y,\n sr = sx * sx + sy * sy + 3;\n var kx = x - sink.x,\n ky = y - sink.y,\n kr = kx * kx + ky * ky + 3;\n var vx = (5 * sx) / sr - (5 * kx) / kr,\n vy = (5 * sy) / sr - (5 * ky) / kr;\n var length = Math.sqrt(vx * vx + vy * vy);\n if (length > 1) {\n vx /= length;\n vy /= length;\n }\n velocityX[indexOf(x, y)] = vx;\n velocityY[indexOf(x, y)] = vy;\n }\n computeDivergence();\n }\n function velocityAt(array, x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n return array[indexOf(x, y)];\n }\n function computeDivergence() {\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var horizontal = (velocityAt(velocityX, x + 1, y) - velocityAt(velocityX, x - 1, y)) * 0.5;\n var vertical = (velocityAt(velocityY, x, y + 1) - velocityAt(velocityY, x, y - 1)) * 0.5;\n divergence[indexOf(x, y)] = horizontal + vertical;\n }\n }\n\n function scalarArrayAt(array, x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n return array[indexOf(x, y)];\n }\n function meanAbsoluteDivergence() {\n var total = 0;\n for (var i = 0; i < divergence.length; i++) total += Math.abs(divergence[i]);\n return total / divergence.length;\n }\n function updateProjectionMetric(before, after) {\n var metric = document.getElementById(\"projectionMetric\");\n if (!metric) return;\n metric.textContent =\n before === null ? after.toFixed(4) : before.toFixed(4) + \" → \" + after.toFixed(4);\n }\n function projectVelocity(commitUI) {\n computeDivergence();\n var before = meanAbsoluteDivergence();\n pressure.fill(0);\n pressureNext.fill(0);\n for (var iteration = 0; iteration < +iterations.value; iteration++) {\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n pressureNext[indexOf(x, y)] =\n (scalarArrayAt(pressure, x - 1, y) +\n scalarArrayAt(pressure, x + 1, y) +\n scalarArrayAt(pressure, x, y - 1) +\n scalarArrayAt(pressure, x, y + 1) -\n divergence[indexOf(x, y)]) /\n 4;\n }\n var swap = pressure;\n pressure = pressureNext;\n pressureNext = swap;\n }\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n velocityX[indexOf(x, y)] -=\n 0.5 * (scalarArrayAt(pressure, x + 1, y) - scalarArrayAt(pressure, x - 1, y));\n velocityY[indexOf(x, y)] -=\n 0.5 * (scalarArrayAt(pressure, x, y + 1) - scalarArrayAt(pressure, x, y - 1));\n }\n computeDivergence();\n var after = meanAbsoluteDivergence();\n updateProjectionMetric(before, after);\n if (commitUI !== false) {\n tickCount++;\n draw();\n updateReadout();\n }\n }\n\n // The math is intentionally kept behind this name for the first checkpoint.\n function diffuseScalarField(commitUI) {\n var strength = +rate.value;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var center = at(x, y),\n sum = 0,\n count = 0;\n if (x > 0) {\n sum += at(x - 1, y);\n count++;\n }\n if (x < WIDTH - 1) {\n sum += at(x + 1, y);\n count++;\n }\n if (y > 0) {\n sum += at(x, y - 1);\n count++;\n }\n if (y < HEIGHT - 1) {\n sum += at(x, y + 1);\n count++;\n }\n next[indexOf(x, y)] = center + strength * (sum / count - center);\n }\n var swap = field;\n field = next;\n next = swap;\n if (commitUI !== false) {\n tickCount++;\n draw();\n updateReadout();\n }\n }\n\n function sampleScalar(x, y) {\n x = Math.max(0, Math.min(WIDTH - 1, x));\n y = Math.max(0, Math.min(HEIGHT - 1, y));\n var x0 = Math.floor(x),\n y0 = Math.floor(y),\n x1 = Math.min(WIDTH - 1, x0 + 1),\n y1 = Math.min(HEIGHT - 1, y0 + 1);\n var tx = x - x0,\n ty = y - y0;\n var top = at(x0, y0) * (1 - tx) + at(x1, y0) * tx;\n var bottom = at(x0, y1) * (1 - tx) + at(x1, y1) * tx;\n return top * (1 - ty) + bottom * ty;\n }\n function advectScalarField(commitUI) {\n var stepDistance = +distance.value * stepScale;\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n var i = indexOf(x, y);\n next[i] = sampleScalar(x - velocityX[i] * stepDistance, y - velocityY[i] * stepDistance);\n }\n var swap = field;\n field = next;\n next = swap;\n if (commitUI !== false) {\n tickCount++;\n draw();\n updateReadout();\n }\n }\n function advectAndDiffuse() {\n advectScalarField(false);\n diffuseScalarField(false);\n tickCount++;\n draw();\n updateReadout();\n }\n\n function colorForDivergence(value) {\n var amount = Math.min(1, Math.abs(value) * 2.2);\n if (value > 0)\n return (\n \"rgb(\" +\n Math.round(22 + 145 * amount) +\n \",\" +\n Math.round(25 + 45 * amount) +\n \",\" +\n Math.round(34 + 25 * amount) +\n \")\"\n );\n if (value < 0)\n return (\n \"rgb(\" +\n Math.round(22 + 45 * amount) +\n \",\" +\n Math.round(25 + 35 * amount) +\n \",\" +\n Math.round(34 + 145 * amount) +\n \")\"\n );\n return \"#111821\";\n }\n function colorFor(value) {\n if (value <= 0) return \"#0b1017\";\n var r = Math.round(40 + value * 65),\n g = Math.round(105 + value * 110),\n b = Math.round(170 + value * 80);\n return \"rgb(\" + r + \",\" + g + \",\" + b + \")\";\n }\n function draw() {\n ctx.clearRect(0, 0, canvas.width, canvas.height);\n for (var y = 0; y < HEIGHT; y++)\n for (var x = 0; x < WIDTH; x++) {\n ctx.fillStyle =\n mode === \"shapedLamp\" && !vesselMask[indexOf(x, y)]\n ? \"#05070a\"\n : mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ? colorForLatentWax(\n temperature[indexOf(x, y)],\n wax[indexOf(x, y)],\n meltProgress[indexOf(x, y)],\n )\n : mode === \"phaseWax\" || mode === \"interfaceWax\"\n ? colorForPhaseWax(temperature[indexOf(x, y)], wax[indexOf(x, y)])\n : mode === \"waxField\" || mode === \"densityWax\"\n ? colorForWax(temperature[indexOf(x, y)], wax[indexOf(x, y)])\n : mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\"\n ? colorForTemperature(temperature[indexOf(x, y)])\n : mode === \"boundaries\" && isBoundaryLeak(x, y)\n ? \"#6b2929\"\n : mode === \"boundaries\" &&\n (x === 0 || x === WIDTH - 1 || y === 0 || y === HEIGHT - 1)\n ? \"#343b45\"\n : mode === \"divergence\" || mode === \"projection\"\n ? colorForDivergence(divergence[indexOf(x, y)])\n : mode === \"grid\" || mode === \"velocity\" || mode === \"boundaries\"\n ? \"#111821\"\n : colorFor(at(x, y));\n ctx.fillRect(x * CELL, y * CELL, CELL, CELL);\n ctx.strokeStyle =\n (mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\") &&\n blobLabels[indexOf(x, y)] >= 0\n ? [\"#50a7ff\", \"#ffd166\", \"#ef6f9a\", \"#75d39b\", \"#b28cff\", \"#ff9f5a\"][\n blobLabels[indexOf(x, y)] % 6\n ]\n : mode === \"pinchOff\" && neckMask[indexOf(x, y)]\n ? \"#7de3ff\"\n : mode === \"surfaceWax\" && wax[indexOf(x, y)] > 0.05 && wax[indexOf(x, y)] < 0.95\n ? \"#7de3ff\"\n : \"#28313c\";\n ctx.strokeRect(x * CELL + 0.5, y * CELL + 0.5, CELL - 1, CELL - 1);\n }\n if (mode === \"shapedLamp\" || mode === \"particleVolume\") drawWaxParticles();\n if (\n mode === \"velocity\" ||\n mode === \"advection\" ||\n mode === \"transportDiffusion\" ||\n mode === \"divergence\" ||\n mode === \"projection\" ||\n mode === \"boundaries\" ||\n mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ) {\n for (var y = 0; y < HEIGHT; y += 2)\n for (var x = 0; x < WIDTH; x += 2) {\n var vx = velocityX[indexOf(x, y)],\n vy = velocityY[indexOf(x, y)];\n drawArrow(x * CELL + CELL / 2, y * CELL + CELL / 2, vx * CELL * 0.72, vy * CELL * 0.72);\n }\n if (mode === \"boundaries\") {\n for (var x = 0; x < WIDTH; x++) {\n drawArrow(\n x * CELL + CELL / 2,\n CELL / 2,\n velocityX[indexOf(x, 0)] * CELL * 0.72,\n velocityY[indexOf(x, 0)] * CELL * 0.72,\n );\n drawArrow(\n x * CELL + CELL / 2,\n (HEIGHT - 1) * CELL + CELL / 2,\n velocityX[indexOf(x, HEIGHT - 1)] * CELL * 0.72,\n velocityY[indexOf(x, HEIGHT - 1)] * CELL * 0.72,\n );\n }\n for (var y = 1; y < HEIGHT - 1; y++) {\n drawArrow(\n CELL / 2,\n y * CELL + CELL / 2,\n velocityX[indexOf(0, y)] * CELL * 0.72,\n velocityY[indexOf(0, y)] * CELL * 0.72,\n );\n drawArrow(\n (WIDTH - 1) * CELL + CELL / 2,\n y * CELL + CELL / 2,\n velocityX[indexOf(WIDTH - 1, y)] * CELL * 0.72,\n velocityY[indexOf(WIDTH - 1, y)] * CELL * 0.72,\n );\n }\n }\n if (selected)\n drawArrow(\n selected.x * CELL + CELL / 2,\n selected.y * CELL + CELL / 2,\n velocityX[indexOf(selected.x, selected.y)] * CELL * 0.72,\n velocityY[indexOf(selected.x, selected.y)] * CELL * 0.72,\n true,\n );\n }\n if (selected) {\n ctx.strokeStyle = \"#ffd166\";\n ctx.lineWidth = 3;\n ctx.strokeRect(selected.x * CELL + 2, selected.y * CELL + 2, CELL - 4, CELL - 4);\n ctx.lineWidth = 1;\n }\n }\n function drawArrow(x, y, dx, dy, highlight) {\n var endX = x + dx,\n endY = y + dy,\n angle = Math.atan2(dy, dx),\n head = 4;\n ctx.strokeStyle = highlight ? \"#ffd166\" : \"#50a7ff\";\n ctx.fillStyle = ctx.strokeStyle;\n ctx.lineWidth = highlight ? 2.5 : 1.4;\n ctx.beginPath();\n ctx.moveTo(x, y);\n ctx.lineTo(endX, endY);\n ctx.stroke();\n if (Math.abs(dx) + Math.abs(dy) > 0.01) {\n ctx.beginPath();\n ctx.moveTo(endX, endY);\n ctx.lineTo(endX - head * Math.cos(angle - 0.55), endY - head * Math.sin(angle - 0.55));\n ctx.lineTo(endX - head * Math.cos(angle + 0.55), endY - head * Math.sin(angle + 0.55));\n ctx.closePath();\n ctx.fill();\n }\n ctx.lineWidth = 1;\n }\n function cellFromEvent(event) {\n var rect = canvas.getBoundingClientRect();\n return {\n x: Math.max(\n 0,\n Math.min(\n WIDTH - 1,\n Math.floor(((event.clientX - rect.left) * canvas.width) / rect.width / CELL),\n ),\n ),\n y: Math.max(\n 0,\n Math.min(\n HEIGHT - 1,\n Math.floor(((event.clientY - rect.top) * canvas.height) / rect.height / CELL),\n ),\n ),\n };\n }\n function relocateWax(cell) {\n if (mode === \"shapedLamp\" && !vesselMask[indexOf(cell.x, cell.y)]) return 0;\n var target = indexOf(cell.x, cell.y),\n needed = 1 - wax[target],\n moved = 0;\n if (needed <= 0) return 0;\n for (var y = HEIGHT - 2; y >= 1 && needed > 1e-6; y--)\n for (var x = 1; x < WIDTH - 1 && needed > 1e-6; x++) {\n var donor = indexOf(x, y);\n if (donor === target || wax[donor] <= 0) continue;\n var take = Math.min(needed, wax[donor]);\n var ratio = wax[donor] > 0 ? take / wax[donor] : 0;\n wax[donor] -= take;\n meltProgress[donor] *= 1 - ratio;\n wax[target] += take;\n needed -= take;\n moved += take;\n }\n document.getElementById(\"inventoryMetric\").textContent = totalWax().toFixed(1);\n document.getElementById(\"movedWaxMetric\").textContent = moved.toFixed(2);\n return moved;\n }\n\n function interact(event) {\n var cell = cellFromEvent(event);\n if (\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ) {\n selected = cell;\n if (painting || event.type === \"pointerdown\") relocateWax(cell);\n } else if (\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ) {\n selected = cell;\n if (painting || event.type === \"pointerdown\") wax[indexOf(cell.x, cell.y)] = 1;\n } else if (\n mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\"\n ) {\n selected = cell;\n if (painting || event.type === \"pointerdown\") temperature[indexOf(cell.x, cell.y)] = 1;\n } else if (\n mode === \"velocity\" ||\n mode === \"divergence\" ||\n mode === \"projection\" ||\n mode === \"boundaries\"\n ) {\n if (!vectorOrigin) vectorOrigin = cell;\n selected = vectorOrigin;\n var dx = cell.x - vectorOrigin.x,\n dy = cell.y - vectorOrigin.y,\n length = Math.sqrt(dx * dx + dy * dy);\n if (length > 0) {\n velocityX[indexOf(selected.x, selected.y)] = dx / length;\n velocityY[indexOf(selected.x, selected.y)] = dy / length;\n if (mode === \"divergence\" || mode === \"projection\") {\n computeDivergence();\n if (mode === \"projection\") updateProjectionMetric(null, meanAbsoluteDivergence());\n }\n if (mode === \"boundaries\") updateBoundaryMetric(null, countBoundaryLeaks());\n }\n } else {\n selected = cell;\n if (mode !== \"grid\" && (painting || event.type === \"pointerdown\")) set(cell.x, cell.y, 1);\n }\n draw();\n updateReadout();\n }\n function updateReadout() {\n document.getElementById(\"coordinate\").textContent = selected\n ? \"(\" + selected.x + \", \" + selected.y + \")\"\n : \"—\";\n document.getElementById(\"index\").textContent = selected ? indexOf(selected.x, selected.y) : \"—\";\n document.getElementById(\"value\").textContent = selected\n ? mode === \"divergence\" || mode === \"projection\"\n ? divergence[indexOf(selected.x, selected.y)].toFixed(3)\n : mode === \"velocity\" || mode === \"boundaries\"\n ? \"(\" +\n velocityX[indexOf(selected.x, selected.y)].toFixed(2) +\n \", \" +\n velocityY[indexOf(selected.x, selected.y)].toFixed(2) +\n \")\"\n : mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n ? \"wax \" +\n wax[indexOf(selected.x, selected.y)].toFixed(2) +\n \" / speed \" +\n Math.sqrt(\n velocityX[indexOf(selected.x, selected.y)] ** 2 +\n velocityY[indexOf(selected.x, selected.y)] ** 2,\n ).toFixed(2)\n : mode === \"viscosityWax\"\n ? \"wax \" +\n wax[indexOf(selected.x, selected.y)].toFixed(2) +\n \" / visc \" +\n viscosityAtIndex(indexOf(selected.x, selected.y)).toFixed(2)\n : mode === \"surfaceWax\"\n ? \"wax \" +\n wax[indexOf(selected.x, selected.y)].toFixed(2) +\n \" / curve \" +\n surfaceCurvature[indexOf(selected.x, selected.y)].toFixed(2)\n : mode === \"latentWax\"\n ? \"melt \" +\n (wax[indexOf(selected.x, selected.y)] > 0\n ? meltProgress[indexOf(selected.x, selected.y)] /\n wax[indexOf(selected.x, selected.y)]\n : 0\n ).toFixed(2) +\n \" / temp \" +\n temperature[indexOf(selected.x, selected.y)].toFixed(2)\n : mode === \"phaseWax\" || mode === \"interfaceWax\"\n ? \"wax \" +\n wax[indexOf(selected.x, selected.y)].toFixed(2) +\n \" / liquid \" +\n liquidWaxFraction(temperature[indexOf(selected.x, selected.y)]).toFixed(2) +\n \" / temp \" +\n temperature[indexOf(selected.x, selected.y)].toFixed(2)\n : mode === \"waxField\" || mode === \"densityWax\"\n ? \"wax \" +\n wax[indexOf(selected.x, selected.y)].toFixed(2) +\n \" / temp \" +\n temperature[indexOf(selected.x, selected.y)].toFixed(2)\n : mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\"\n ? temperature[indexOf(selected.x, selected.y)].toFixed(3) +\n \" / vy \" +\n velocityY[indexOf(selected.x, selected.y)].toFixed(2)\n : at(selected.x, selected.y).toFixed(3)\n : \"—\";\n document.getElementById(\"ticks\").textContent = tickCount;\n }\n\n // LESSON CONTROLLER AND DOM EVENTS\n\n function selectMode(nextMode) {\n if (\n (mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\") &&\n nextMode !== mode\n )\n setRunning(false);\n mode = nextMode;\n var c = copy[mode];\n document.querySelectorAll(\".lesson-button\").forEach(function (b) {\n b.classList.toggle(\"active\", b.dataset.mode === mode);\n });\n document.getElementById(\"sectionNumber\").textContent = c.number;\n document.getElementById(\"sectionTitle\").textContent = c.title;\n document.getElementById(\"sectionText\").textContent = c.text;\n document.getElementById(\"instruction\").textContent = c.instruction;\n document.getElementById(\"code\").innerHTML = c.code;\n document.getElementById(\"notice\").textContent = c.notice;\n document.getElementById(\"notYet\").textContent = c.notYet;\n document.getElementById(\"advance\").disabled = !(\n mode === \"diffusion\" ||\n mode === \"advection\" ||\n mode === \"transportDiffusion\" ||\n mode === \"projection\" ||\n mode === \"boundaries\" ||\n mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\"\n );\n document\n .getElementById(\"diffusionControl\")\n .classList.toggle(\"visible\", mode === \"diffusion\" || mode === \"transportDiffusion\");\n document\n .getElementById(\"advectionControl\")\n .classList.toggle(\n \"visible\",\n mode === \"advection\" ||\n mode === \"transportDiffusion\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document.getElementById(\"projectionControl\").classList.toggle(\"visible\", mode === \"projection\");\n document.getElementById(\"boundaryControl\").classList.toggle(\"visible\", mode === \"boundaries\");\n document\n .getElementById(\"buoyancyControl\")\n .classList.toggle(\n \"visible\",\n mode === \"temperature\" ||\n mode === \"thermalFlow\" ||\n mode === \"coupledThermal\" ||\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\",\n );\n document\n .getElementById(\"heatCycleControl\")\n .classList.toggle(\n \"visible\",\n mode === \"heatedFlow\" ||\n mode === \"conductiveFlow\" ||\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"conductionControl\")\n .classList.toggle(\n \"visible\",\n mode === \"conductiveFlow\" ||\n mode === \"waxField\" ||\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"densityControl\")\n .classList.toggle(\n \"visible\",\n mode === \"densityWax\" ||\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"phaseControl\")\n .classList.toggle(\n \"visible\",\n mode === \"phaseWax\" ||\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"interfaceControl\")\n .classList.toggle(\n \"visible\",\n mode === \"interfaceWax\" ||\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"latentControl\")\n .classList.toggle(\n \"visible\",\n mode === \"latentWax\" ||\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"surfaceControl\")\n .classList.toggle(\n \"visible\",\n mode === \"surfaceWax\" ||\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"viscosityControl\")\n .classList.toggle(\n \"visible\",\n mode === \"viscosityWax\" ||\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"stabilityControl\")\n .classList.toggle(\n \"visible\",\n mode === \"stableLamp\" ||\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"simulatorControl\")\n .classList.toggle(\n \"visible\",\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document\n .getElementById(\"unevenControl\")\n .classList.toggle(\"visible\", mode === \"unevenHeating\" || mode === \"pinchOff\");\n document\n .getElementById(\"pinchControl\")\n .classList.toggle(\"visible\", mode === \"pinchOff\" || mode === \"conservation\");\n document\n .getElementById(\"conservationControl\")\n .classList.toggle(\"visible\", mode === \"conservation\");\n document\n .getElementById(\"energyBudgetControl\")\n .classList.toggle(\"visible\", mode === \"energyBudget\");\n document\n .getElementById(\"closedVesselControl\")\n .classList.toggle(\"visible\", mode === \"closedVessel\" || mode === \"cycleDiagnostics\");\n document\n .getElementById(\"cycleDiagnosticsControl\")\n .classList.toggle(\"visible\", mode === \"cycleDiagnostics\" || mode === \"blobTracking\");\n document\n .getElementById(\"blobTrackingControl\")\n .classList.toggle(\"visible\", mode === \"blobTracking\" || mode === \"blobIdentity\");\n document\n .getElementById(\"blobIdentityControl\")\n .classList.toggle(\"visible\", mode === \"blobIdentity\" || mode === \"motionMatch\");\n document\n .getElementById(\"motionMatchControl\")\n .classList.toggle(\"visible\", mode === \"motionMatch\");\n document\n .getElementById(\"calibrationControl\")\n .classList.toggle(\"visible\", mode === \"calibration\");\n document\n .getElementById(\"physicalUnitsControl\")\n .classList.toggle(\"visible\", mode === \"physicalUnits\");\n document.getElementById(\"finalLampControl\").classList.toggle(\"visible\", mode === \"finalLamp\");\n document\n .getElementById(\"shapedLampControl\")\n .classList.toggle(\"visible\", mode === \"shapedLamp\" || mode === \"particleVolume\");\n document\n .getElementById(\"particleVolumeControl\")\n .classList.toggle(\"visible\", mode === \"particleVolume\");\n document\n .getElementById(\"play\")\n .classList.toggle(\n \"visible\",\n mode === \"simulator\" ||\n mode === \"unevenHeating\" ||\n mode === \"pinchOff\" ||\n mode === \"conservation\" ||\n mode === \"energyBudget\" ||\n mode === \"closedVessel\" ||\n mode === \"cycleDiagnostics\" ||\n mode === \"blobTracking\" ||\n mode === \"blobIdentity\" ||\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\" ||\n mode === \"shapedLamp\" ||\n mode === \"particleVolume\",\n );\n document.querySelectorAll(\".step\").forEach(function (s, i) {\n var current =\n mode === \"particleVolume\"\n ? i === 34\n : mode === \"shapedLamp\"\n ? i === 33\n : mode === \"finalLamp\"\n ? i === 32\n : mode === \"physicalUnits\"\n ? i === 31\n : mode === \"calibration\"\n ? i === 30\n : mode === \"motionMatch\"\n ? i === 29\n : mode === \"blobIdentity\"\n ? i === 28\n : mode === \"blobTracking\"\n ? i === 27\n : mode === \"cycleDiagnostics\"\n ? i === 26\n : mode === \"closedVessel\"\n ? i === 25\n : mode === \"energyBudget\"\n ? i === 24\n : mode === \"conservation\"\n ? i === 23\n : mode === \"pinchOff\"\n ? i === 22\n : mode === \"unevenHeating\"\n ? i === 21\n : mode === \"simulator\"\n ? i === 20\n : mode === \"stableLamp\"\n ? i === 19\n : mode === \"viscosityWax\"\n ? i === 18\n : mode === \"surfaceWax\"\n ? i === 17\n : mode === \"latentWax\"\n ? i === 16\n : mode === \"interfaceWax\"\n ? i === 15\n : mode === \"phaseWax\"\n ? i === 14\n : mode === \"densityWax\"\n ? i === 13\n : mode === \"waxField\"\n ? i === 12\n : mode === \"conductiveFlow\"\n ? i === 11\n : mode === \"heatedFlow\"\n ? i === 10\n : mode === \"coupledThermal\"\n ? i === 9\n : mode === \"thermalFlow\"\n ? i === 8\n : mode === \"temperature\"\n ? i === 7\n : mode === \"boundaries\"\n ? i === 6\n : mode === \"projection\"\n ? i === 5\n : mode === \"divergence\"\n ? i === 4\n : mode ===\n \"transportDiffusion\"\n ? i === 3\n : mode === \"advection\"\n ? i === 2\n : mode === \"velocity\"\n ? i === 1\n : i === 0;\n s.classList.toggle(\"current\", current);\n });\n document.getElementById(\"truth\").innerHTML =\n mode === \"particleVolume\"\n ? \"Current claim: neighbor-density pressure preserves finite wax volume while cohesion keeps it connected.\"\n : mode === \"shapedLamp\"\n ? \"Current claim: tapered walls and localized heating can return falling wax to a repeating collection-and-launch cycle.\"\n : mode === \"finalLamp\"\n ? \"Current claim: this is the complete educational lava-lamp simulator assembled from the preceding fields, forces, constraints, and diagnostics.\"\n : mode === \"physicalUnits\"\n ? \"Current claim: reference length and time scales convert normalized motion into an interpretable physical speed.\"\n : mode === \"calibration\"\n ? \"Current claim: useful lava-lamp behavior occupies a coordinated parameter regime rather than one magic coefficient.\"\n : mode === \"motionMatch\"\n ? \"Current claim: a spatial tolerance can preserve identity when a blob moves beyond direct cell overlap.\"\n : mode === \"blobIdentity\"\n ? \"Current claim: spatial overlap can preserve a derived blob identity across successive field states.\"\n : mode === \"blobTracking\"\n ? \"Current claim: blobs are connected components derived from the wax-fraction field at a chosen threshold.\"\n : mode === \"cycleDiagnostics\"\n ? \"Current claim: repeated changes in wax center height reveal a circulation cycle that no single simulation frame can establish.\"\n : mode === \"closedVessel\"\n ? \"Current claim: wax fraction and carrier fraction sum to one in every cell, while the closed lamp keeps one fixed total wax inventory.\"\n : mode === \"energyBudget\"\n ? \"Current claim: thermal energy changes because the lamp is an open system; its change can be separated into heater input, cooling loss, latent storage, and numerical drift.\"\n : mode === \"conservation\"\n ? \"Current claim: local field values may change freely while a global correction keeps their total wax constant.\"\n : mode === \"pinchOff\"\n ? \"Current claim: a conservative topology correction can teach pinch-off on a grid too coarse to resolve the physical instability directly.\"\n : mode === \"unevenHeating\"\n ? \"Current claim: small deterministic heater differences break symmetry and seed localized buoyant wax regions from a shared bottom pool.\"\n : mode === \"simulator\"\n ? \"Current claim: this continuously runs the complete educational lava-lamp timestep with coherent presets and live stability diagnostics.\"\n : mode === \"stableLamp\"\n ? \"Current claim: subdividing the same frame time recomputes coupled fields more frequently and reduces per-substep travel.\"\n : mode === \"viscosityWax\"\n ? \"Current claim: phase-dependent viscosity smooths neighboring velocity differences at a rate derived from material and melt progress.\"\n : mode === \"surfaceWax\"\n ? \"Current claim: an interface-normal and curvature estimate can create a surface force that encourages coherent wax blobs.\"\n : mode === \"latentWax\"\n ? \"Current claim: latent heat allows phase progress to change while temperature remains near the melting point.\"\n : mode === \"interfaceWax\"\n ? \"Current claim: fractional wax cells must render fractionally, and interface sharpening can counter numerical smearing at a conservation cost.\"\n : mode === \"phaseWax\"\n ? \"Current claim: wax identity can remain constant while a temperature-derived liquid fraction controls resistance to motion.\"\n : mode === \"densityWax\"\n ? \"Current claim: combining wax fraction with temperature produces a signed density force: hot wax rises and cold wax sinks relative to the carrier.\"\n : mode === \"waxField\"\n ? \"Current claim: material identity can be stored as a scalar fraction and transported independently from temperature through the same velocity field.\"\n : mode === \"conductiveFlow\"\n ? \"Current claim: conduction spreads temperature through neighbor exchange while wall cooling removes heat at the domain boundary.\"\n : mode === \"heatedFlow\"\n ? \"Current claim: sustained thermal motion depends on the balance between heat added each step and heat lost each step.\"\n : mode === \"coupledThermal\"\n ? \"Current claim: projection spreads a local buoyant push into a coherent velocity field before advection transports temperature.\"\n : mode === \"thermalFlow\"\n ? \"Current claim: buoyancy changes velocity, then advection uses that velocity to carry the temperature field upward.\"\n : mode === \"temperature\"\n ? \"Current claim: a temperature scalar field can add upward velocity through a buoyancy force.\"\n : mode === \"boundaries\"\n ? \"Current claim: no-through boundary enforcement removes outward wall-normal velocity while preserving tangential motion.\"\n : mode === \"projection\"\n ? \"Current claim: an iterative pressure solve reduces the measured divergence when its gradient is removed from velocity.\"\n : mode === \"divergence\"\n ? \"Current claim: this computes discrete divergence from the velocity field and visualizes sources and sinks.\"\n : mode === \"transportDiffusion\"\n ? \"Current claim: each step transports the dye, then diffuses it into neighboring cells.\"\n : mode === \"advection\"\n ? \"Current claim: this transports a scalar field through a fixed velocity field, one discrete step at a time.\"\n : mode === \"velocity\"\n ? \"Current claim: this stores and renders a velocity vector in each cell. It does not move dye yet.\"\n : \"Current claim: this demonstrates a grid, a scalar field, and neighbor-based smoothing. It is not yet a fluid simulation.\";\n if (mode === \"finalLamp\") {\n document.querySelectorAll(\".slider\").forEach(function (x) {\n x.classList.remove(\"visible\");\n });\n document.getElementById(\"simulatorControl\").classList.add(\"visible\");\n document.getElementById(\"finalLampControl\").classList.add(\"visible\");\n }\n if (mode === \"particleVolume\") {\n document.querySelectorAll(\".slider\").forEach(function (x) {\n x.classList.remove(\"visible\");\n });\n document.getElementById(\"simulatorControl\").classList.add(\"visible\");\n document.getElementById(\"heatCycleControl\").classList.add(\"visible\");\n document.getElementById(\"shapedLampControl\").classList.add(\"visible\");\n document.getElementById(\"particleVolumeControl\").classList.add(\"visible\");\n }\n if (mode === \"shapedLamp\") {\n document.querySelectorAll(\".slider\").forEach(function (x) {\n x.classList.remove(\"visible\");\n });\n document.getElementById(\"simulatorControl\").classList.add(\"visible\");\n document.getElementById(\"heatCycleControl\").classList.add(\"visible\");\n document.getElementById(\"shapedLampControl\").classList.add(\"visible\");\n }\n draw();\n updateReadout();\n }\n\n canvas.addEventListener(\"pointerdown\", function (e) {\n painting = true;\n canvas.setPointerCapture(e.pointerId);\n interact(e);\n });\n canvas.addEventListener(\"pointermove\", function (e) {\n if (painting) interact(e);\n });\n canvas.addEventListener(\"pointerup\", function () {\n painting = false;\n vectorOrigin = null;\n });\n document.querySelectorAll(\".lesson-button\").forEach(function (b) {\n b.addEventListener(\"click\", function () {\n selectMode(b.dataset.mode);\n });\n });\n document.getElementById(\"clear\").addEventListener(\"click\", clearCurrent);\n document.getElementById(\"seed\").addEventListener(\"click\", function () {\n seedCurrent();\n if (mode === \"blobTracking\") analyzeBlobs(false);\n });\n document.getElementById(\"advance\").addEventListener(\"click\", function () {\n if (mode === \"diffusion\") diffuseScalarField();\n else if (mode === \"advection\") advectScalarField();\n else if (mode === \"transportDiffusion\") advectAndDiffuse();\n else if (mode === \"projection\") projectVelocity();\n else if (mode === \"boundaries\") enforceBoundaries();\n else if (mode === \"temperature\") applyBuoyancy();\n else if (mode === \"thermalFlow\") stepRisingHeat();\n else if (mode === \"coupledThermal\") stepCoupledThermal();\n else if (mode === \"heatedFlow\") stepHeatedFlow();\n else if (mode === \"conductiveFlow\") stepConductiveFlow();\n else if (mode === \"waxField\") stepWaxTransport();\n else if (mode === \"densityWax\") stepDensityWax();\n else if (mode === \"phaseWax\") stepPhaseWax();\n else if (mode === \"interfaceWax\") stepInterfaceWax();\n else if (mode === \"latentWax\") stepLatentWax();\n else if (mode === \"surfaceWax\") stepSurfaceWax();\n else if (mode === \"viscosityWax\") stepViscosityWax();\n else if (mode === \"stableLamp\") stepStableLamp();\n else if (mode === \"simulator\") stepFinalLamp();\n else if (mode === \"unevenHeating\") stepUnevenLamp();\n else if (mode === \"pinchOff\") stepPinchLamp();\n else if (mode === \"conservation\") stepConservationLamp();\n else if (mode === \"energyBudget\") stepEnergyBudgetLamp();\n else if (mode === \"closedVessel\") stepClosedVesselLamp();\n else if (mode === \"cycleDiagnostics\") stepCycleDiagnosticsLamp();\n else if (mode === \"blobTracking\") stepBlobTrackingLamp();\n else if (mode === \"blobIdentity\") stepBlobIdentityLamp();\n else if (\n mode === \"motionMatch\" ||\n mode === \"calibration\" ||\n mode === \"physicalUnits\" ||\n mode === \"finalLamp\"\n )\n stepRemainingLesson();\n else if (mode === \"shapedLamp\") stepShapedLamp();\n else if (mode === \"particleVolume\") stepParticleVolume();\n });\n rate.addEventListener(\"input\", function () {\n document.getElementById(\"rateOut\").textContent = Math.round(+rate.value * 100) + \"%\";\n });\n distance.addEventListener(\"input\", function () {\n document.getElementById(\"distanceOut\").textContent = (+distance.value).toFixed(2) + \" cells\";\n });\n iterations.addEventListener(\"input\", function () {\n document.getElementById(\"iterationsOut\").textContent = iterations.value;\n });\n buoyancy.addEventListener(\"input\", function () {\n document.getElementById(\"buoyancyOut\").textContent = (+buoyancy.value).toFixed(2);\n });\n heater.addEventListener(\"input\", function () {\n document.getElementById(\"heaterOut\").textContent = (+heater.value).toFixed(2);\n });\n cooling.addEventListener(\"input\", function () {\n document.getElementById(\"coolingOut\").textContent = Math.round(+cooling.value * 100) + \"%\";\n });\n conduction.addEventListener(\"input\", function () {\n document.getElementById(\"conductionOut\").textContent =\n Math.round(+conduction.value * 100) + \"%\";\n });\n wallCooling.addEventListener(\"input\", function () {\n document.getElementById(\"wallCoolingOut\").textContent =\n Math.round(+wallCooling.value * 100) + \"%\";\n });\n neutralTemp.addEventListener(\"input\", function () {\n document.getElementById(\"neutralOut\").textContent = (+neutralTemp.value).toFixed(2);\n });\n densityResponse.addEventListener(\"input\", function () {\n document.getElementById(\"densityOut\").textContent = (+densityResponse.value).toFixed(2);\n });\n meltTemp.addEventListener(\"input\", function () {\n document.getElementById(\"meltOut\").textContent = (+meltTemp.value).toFixed(2);\n });\n transitionWidth.addEventListener(\"input\", function () {\n document.getElementById(\"transitionOut\").textContent = (+transitionWidth.value).toFixed(2);\n });\n solidDrag.addEventListener(\"input\", function () {\n document.getElementById(\"solidDragOut\").textContent = Math.round(+solidDrag.value * 100) + \"%\";\n });\n interfaceSharpness.addEventListener(\"input\", function () {\n document.getElementById(\"interfaceOut\").textContent =\n Math.round(+interfaceSharpness.value * 100) + \"%\";\n });\n latentCoupling.addEventListener(\"input\", function () {\n document.getElementById(\"latentOut\").textContent =\n Math.round(+latentCoupling.value * 100) + \"%\";\n });\n surfaceTension.addEventListener(\"input\", function () {\n document.getElementById(\"surfaceOut\").textContent = (+surfaceTension.value).toFixed(2);\n });\n carrierViscosity.addEventListener(\"input\", function () {\n document.getElementById(\"carrierViscosityOut\").textContent =\n Math.round(+carrierViscosity.value * 100) + \"%\";\n });\n hotViscosity.addEventListener(\"input\", function () {\n document.getElementById(\"hotViscosityOut\").textContent =\n Math.round(+hotViscosity.value * 100) + \"%\";\n });\n coldViscosity.addEventListener(\"input\", function () {\n document.getElementById(\"coldViscosityOut\").textContent =\n Math.round(+coldViscosity.value * 100) + \"%\";\n });\n timeStep.addEventListener(\"input\", function () {\n document.getElementById(\"timeStepOut\").textContent = (+timeStep.value).toFixed(2);\n });\n substeps.addEventListener(\"input\", function () {\n document.getElementById(\"substepsOut\").textContent = substeps.value;\n });\n playbackRate.addEventListener(\"input\", function () {\n document.getElementById(\"playbackOut\").textContent = playbackRate.value;\n });\n lampPreset.addEventListener(\"change\", function () {\n applyLampPreset();\n clearCurrent();\n seedCurrent();\n updateSimulatorDiagnostics();\n });\n document.getElementById(\"play\").addEventListener(\"click\", function () {\n setRunning(!running);\n });\n heaterUnevenness.addEventListener(\"input\", function () {\n document.getElementById(\"unevenOut\").textContent =\n Math.round(+heaterUnevenness.value * 100) + \"%\";\n });\n neckSensitivity.addEventListener(\"input\", function () {\n document.getElementById(\"neckOut\").textContent = Math.round(+neckSensitivity.value * 100) + \"%\";\n });\n pinchStrength.addEventListener(\"input\", function () {\n document.getElementById(\"pinchOut\").textContent = Math.round(+pinchStrength.value * 100) + \"%\";\n });\n conservationStrength.addEventListener(\"input\", function () {\n document.getElementById(\"conservationOut\").textContent =\n Math.round(+conservationStrength.value * 100) + \"%\";\n });\n restDensity.addEventListener(\"input\", function () {\n document.getElementById(\"restDensityOut\").textContent = (+restDensity.value).toFixed(1);\n });\n volumePressure.addEventListener(\"input\", function () {\n document.getElementById(\"volumePressureOut\").textContent =\n Math.round(+volumePressure.value * 100) + \"%\";\n });\n particleCount.addEventListener(\"input\", function () {\n document.getElementById(\"particleCountOut\").textContent = particleCount.value;\n });\n particleCount.addEventListener(\"change\", function () {\n setRunning(false);\n seedCurrent();\n recordCycleDiagnostics();\n analyzeBlobs(false, true);\n draw();\n });\n baseFocus.addEventListener(\"input\", function () {\n document.getElementById(\"baseFocusOut\").textContent = Math.round(+baseFocus.value * 100) + \"%\";\n });\n collectionStrength.addEventListener(\"input\", function () {\n document.getElementById(\"collectionOut\").textContent =\n Math.round(+collectionStrength.value * 100) + \"%\";\n });\n glassCooling.addEventListener(\"input\", function () {\n document.getElementById(\"glassCoolingOut\").textContent =\n Math.round(+glassCooling.value * 100) + \"%\";\n });\n matchRadius.addEventListener(\"input\", function () {\n document.getElementById(\"matchRadiusOut\").textContent =\n matchRadius.value + \" cell\" + (matchRadius.value === \"1\" ? \"\" : \"s\");\n });\n cellSizeCm.addEventListener(\"input\", function () {\n document.getElementById(\"cellSizeOut\").textContent = (+cellSizeCm.value).toFixed(1) + \" cm\";\n updateCourseDiagnostics();\n });\n secondsPerFrame.addEventListener(\"input\", function () {\n document.getElementById(\"secondsOut\").textContent = (+secondsPerFrame.value).toFixed(2) + \" s\";\n updateCourseDiagnostics();\n });\n blobThreshold.addEventListener(\"input\", function () {\n document.getElementById(\"blobThresholdOut\").textContent =\n Math.round(+blobThreshold.value * 100) + \"%\";\n analyzeBlobs(false);\n draw();\n });\n applyLampPreset();\n buildLampMask();\n selectMode(\"particleVolume\");\n seedCurrent();\n document.getElementById(\"inventoryMetric\").textContent = totalWax().toFixed(1);\n recordCycleDiagnostics();\n analyzeBlobs(false, true);\n updateCourseDiagnostics();\n updateSimulatorDiagnostics();\n requestAnimationFrame(animationLoop);\n})();\n\nexport {};\n",[],[149,257,275],"efea0db5-813a-4c4b-9b5e-662afae10362",244,["D",1785252127214],[168,169,170,171],[149,257],["D",1784854699691],[167,172,176,179],[149],["D",1784854669144],{"_96":95,"_103":99},{"_96":95,"_97":-5,"_98":99,"_18":54,"_100":101,"_102":24},{"_13":85,"_31":86,"_42":41,"_87":43,"_88":-5,"_89":-5,"_90":-5,"_91":-5,"_92":-5,"_93":44,"_94":95},"ctx",{"_78":-5,"_123":327,"_328":329,"_330":331,"_332":24,"_333":24,"_334":24,"_335":58,"_336":-5,"_337":-7,"_338":-7,"_339":62,"_340":24,"_341":24,"_342":-5,"_343":24,"_344":24,"_345":24,"_346":-5},"ENV",{"_298":299,"_300":301,"_302":86,"_303":304,"_305":306,"_307":308,"_309":310,"_311":312,"_313":314,"_315":316,"_317":318,"_319":320,"_321":322,"_323":324,"_325":326},"REACT_APP_SENTRY_DSN","https://586020d86f05489db63a554690acf1e9@o1398359.ingest.sentry.io/6724499","RENDER_GIT_COMMIT","fdb67791413f4386b4dbf18f354029a617f6160f","RENDER_GIT_BRANCH","POSTHOG_API_TOKEN","phc_nmoyuAkStEZuUp69rYlg1DZgrIt19Gx0Y4Hwf8tGgJP","POSTHOG_API_HOST","https://wart.val.town","WEB_API_TEMPLATE","https://{val}.web.val.run","APP_URL","https://www.val.town","MODULE_URL","https://esm.town","VANITY_WEB_API_TEMPLATE","https://{label}.val.run/","DENOLS_URL","https://lsp.val.town","NODE_ENV","production","CLOUDFLARE_TURNSTILE_SITE_KEY","0x4AAAAAACtaIZJH7s9rBMZ6","GOOGLE_CLIENT_ID","875989696899-jtsqkhlbdv5jrlsk8631p7ipjg3c6o62.apps.googleusercontent.com","STRIPE_PUBLISHABLE_KEY","pk_live_51N2zvWEpZTTLipIkVEwXrUG9AT0yeOOEYxszarZT1jpsZV58Ly6wlg3SltZpHN26CUalOBmUSY7adQRQiSkCZz2500qSXNkf8Y","API_URL","https://api.val.town",[],"preferencesKeybindings","default","preferencesTheme","light","preferencesFixers","preferencesCopilot","preferencesSelectionAi","preferencesCloseBrackets","preferencesChatPanelState","preferencesDprint","preferencesDenofmt","notificationCount","flagDenols","isBusinessOrgMember","businessOrgHandle","isBusinessOrgAdmin","isHigherTownieTier","navCollapsed","currentOrgId"]

Diff · v275