  We investigate the consequences of space-momentum correlations in quark
phase-space distributions for coalescence processes at the hadronization
transition. Thus far it has been proved difficult to reconcile such
correlations with the empirically observed constituent quark number scaling
(CQNS) at the Relativistic Heavy-Ion Collider (RHIC). To address this problem
we combine our earlier developed quark recombination model with quark
phase-space distributions computed from relativistic Langevin simulations in an
expanding Quark-Gluon Plasma (QGP). Hadronization is based on resonance
formation within a Boltzmann equation which recovers thermal equilibrium and
obeys energy conservation in the quark-coalescence process, while the fireball
background is adjusted to hydrodynamic simulations of semi-central Au-Au
collisions at RHIC. To facilitate the applicability of the Langevin process, we
focus on strange and charm quarks. Their interactions in the QGP are modeled
using leading-order perturbative QCD augmented by effective Lagrangians with
resonances which smoothly merge into hadronic states formed at T_c. The
interaction strength is adjusted to reproduce the empirical saturation value
for the quark-elliptic flow, v_{2,q}^{sat}~7-8%. The resulting phi and J/\psi
elliptic flow recover CQNS over a large range in transverse momentum (p_T)
within a few percent. As a function of transverse kinetic energy, both the
quark spectra from the Langevin simulations and the meson spectra generated via
resonance recombination recover CQNS from zero to at least 3 GeV.
