We consider the Boussinesq PDE perturbed by a time-dependent forcing. Even though there is no smoothing effect for arbitrary smooth initial data, we are able to apply the method of self-consistent bounds to deduce the existence of smooth classical periodic solutions in the vicinity of 0. The proof is non-perturbative and relies on construction of periodic isolating segments in the Galerkin projections.
We present a Lohner type algorithm for the computation of rigorous bounds for the solutions of ordinary differential equations and its derivatives with respect to the initial conditions up to an arbitrary order.
JavaScript is turned off in your web browser. Turn it on to take full advantage of this site, then refresh the page.