We study certain multilayer
free-boundary problems, in which the layer interfaces constitute a nested family of
convex closed surfaces, each characterized by a Bernoulli joining condition between
the potentials in the neighboring layers. In this context, we develop convex
variational methods based on a family of convexity-preserving free-boundary
perturbation operators, and we apply these methods in the study of the existence of
convex solutions.