  Simulations of high-density deuterium plasmas in a lower single-null magnetic
configuration based on a TCV discharge are presented. We evolve the dynamics of
three charged species (electrons, D$^{+}$ and D$_{2}^{+}$), interacting with
two neutrals species (D and D$_2$) through ionization, charge-exchange,
recombination and molecular dissociation processes. The plasma is modelled by
using the drift-reduced fluid Braginskii equations, while the neutral dynamics
is described by a kinetic model. To control the divertor conditions, a D$_2$
puffing is used and the effect of increasing the puffing strength is
investigated. The increase in fuelling leads to an increase of density in the
scrape-off layer and a decrease of the plasma temperature. At the same time,
the particle and heat fluxes to the divertor target decrease and the detachment
of the inner target is observed. The analysis of particle and transport balance
in the divertor volume shows that the decrease of the particle flux is caused
by a decrease of the local neutral ionization together with a decrease of the
parallel velocity, both caused by the lower plasma temperature. The relative
importance of the different collision terms is assessed, showing the crucial
role of molecular interactions, as they are responsible for increasing the
atomic neutral density and temperature, since most of the D neutrals are
produced by molecular activated recombination and D$_2$ dissociation. The
presence of strong electric fields in high-density plasmas is also shown,
revealing the role of the $E \times B$ drift in setting the asymmetry between
the divertor targets. Simulation results are in agreement with experimental
observations of increased density decay length, attributed to a decrease of
parallel transport, together with an increase of plasma blob size and radial
velocity.
