Evolution of flux-closure domain arrays in oxide multilayers with misfit strain