A number of misunderstandings about modeling the apparent accelerated expansion of the Universe, and about the `weak singularity' are clarified: 1. Of the five definitions of the deceleration parameter given by Hirata and Seljak (HS), only $q_1$ is a correct invariant measure of acceleration/deceleration of expansion. The $q_3$ and $q_4$ are unrelated to acceleration in an inhomogeneous model. 2. The averaging over directions involved in the definition of $q_4$ does not correspond to what is done in observational astronomy. 3. HS's equation (38) connecting $q_4$ to the flow invariants gives self-contradictory results when applied at the centre of symmetry of the Lema\^{\i}tre-Tolman (L-T) model. The intermediate equation (31) that determines $q_{3'}$ is correct, but approximate, so it cannot be used for determining the sign of the deceleration parameter. Even so, at the centre of symmetry of the L-T model, it puts no limitation on the sign of $q_{3'}(0)$. 4. The `weak singularity' of Vanderveld {\it et al.} is a conical profile of mass density at the centre - a perfectly acceptable configuration. 5. The so-called `critical point' in the equations of the `inverse problem' for a central observer in an L-T model is a manifestation of the apparent horizon - a common property of the past light cones in zero-lambda L-T models, perfectly manageable if the equations are correctly integrated.

We develop models of void formation starting from a small initial fluctuation at recombination and growing to a realistic present day density profile in agreement with observations of voids. The model construction is an extension of previously developed algorithms for finding a Lemaitre-Tolman metric that evolves between two profiles of either density or velocity specified at two times. Of the 4 profiles of concern -- those of density and velocity at recombination and at the present day -- two can be specified and the other two follow from the derived model. We find that, in order to reproduce the present-day void density profiles, the initial velocity profile is more important than the initial density profile. Extrapolation of current CMB observations to the scales relevant to proto-voids is very uncertain. Even so, we find that it is very difficult to make both the initial density and velocity fluctuation amplitudes small enough, and still obtain a realistic void by today.

Structure formation within the Lemaitre-Tolman model is investigated in a general manner. We seek models such that the initial density perturbation within a homogeneous background has a smaller mass than the structure into which it will develop, and the perturbation then accretes more mass during evolution. This is a generalisation of the approach taken by Bonnor in 1956. It is proved that any two spherically symmetric density profiles specified on any two constant time slices can be joined by a Lemaitre-Tolman evolution, and exact implicit formulae for the arbitrary functions that determine the resulting L-T model are obtained. Examples of the process are investigated numerically.