How to extend SecondOrderODEProblem in DifferentialEquations.jl?

The Westervelt equation looks like a partial differential equation so you would need to discretize the spatial dimensions somehow (finite elements, finite volume, spectral methods etc).

Also Please read: make it easier to help you