Basic heat and wave equation solver for node and the browser
A simple solver of the heat/diffusion equation and wave equation meant for driving visualizations not for performing exacting scientific analysis. The forward Euler step is employed for simplicity and speed, as well as further simplifying the discrete equations with the convolution operator.
It solves the pde on an nxm square grid with n and m being defined by the
io = require'socket.io'listen6543engine = require'pdeEngine'field = enginedt: 0.1gamma: 0.02eqn: "wave"fieldsetResolution80 100iosocketson'connection'socketon'addSource'fieldaddSourcedatarow datacol magsetIntervalvar coeffs = fieldupdatedt: 0.1eqn: "diffusion"alpha: 0.5// Render coeffs //50
Of course you can employ any means of adding sources.
addSource method specifies the location of the source at that particular timestep with magnitude
mag. Right now it emplaces a hard-coded Gaussian source centred at
(row, col), with boundary cases already implemented.
The following are possible configuration options. If configs are omitted defaults will be used. The defaults are given after the
dt = specdt || 0.1dx = specdx || 1gamma = specgamma || 0.02 // wave decay factorvel = specvel || 2 // wave velocityalpha = specalpha || 1 // diffusion paramtereqn = speceqn || 'wave'
This is an early but fully working physics engine. Future releases may optimize algorithms or add new methods/configs, but this basic API should remain stable.
Tests will be coming I have them already for my own purposes. Just need to refactor and include.