bessel_neumann_wave_equation
Graham Dennis
Solve the wave equation on a disk of radius 1 utilising cylindrical symmetry by using the Bessel function transform. The boundary condition at r=1 is d/dr u = 0
const real T = 10.0; const real mass = 1e-3; const real length = 1.0; const real mu = mass/length;
t
u uDot
u = exp(-100.0*(x-0.25)*(x-0.25)); uDot = 0.0;
100 100
L
L = -T*kx*kx/mu;
main
du_dt = uDot; duDot_dt = L[u];
amp
main
amp = u.Re();
amp
main
amp = u.Re();