Conversion to Polar CoordinatesConverting to cylindrical coordinates,
CoordinatesFromCartesian[{x,y,z},Cylindrical]
SetCoordinates[Cylindrical]
the wave equation becomes
waveCylindrical :=The vibrations of a drumhead satisfy the wave equation, where u is the displacement from the horizontal equilibrium position. If the initial conditions are independent of theta, then the solution will also be indepenent of theta. This simplifies the wave equation.D[u[r,theta,t],{t,2}] == Div[c^2 Grad[u[r,theta,t]]]
waveCylindrical
Clear[u]For a drum of radius one, the boundary condition at the outer edge would bewave := D[u[r,t],{t,2}] == Div[c^2 Grad[u[r,t]]]
wave
u(1,t) = 0, t > 0,
since the drumhead does not move at r = 1. We need another boundary condition, since the partial differential equation is second order in r. We will be able to determine an appropriate boudnary condition as we solve the differential equation.
We need two initial conditions, since the equation is second order in t. Let's assume the drumhead is initially displaced from equilibrium, but with no initial velocity.
u(r,0) = u0(r)
(d/dt) u(r,0) = 0