Valence double parton distribution functions of the nucleon are evaluated in the framework of a simple model, where the conservation of the longitudinal momentum is taken into account. The leading-order DGLAP QCD evolution from the low quark-model scale to higher renormalization scales is carried out via the Mellin moments of the distributions. Results of the valence quark correlation function show that in general the double distributions cannot be approximated as a product of the single-particle distributions.Theoretical interest in double parton distributions (dPDFs) in hadrons has been recently renewed (see the reviews [1,2] and the references therein). The activity is largely triggered by the experimental program at the LHC, where the double parton scattering (DPS) may bring substantial contribution in certain production processes [3,4]. Phenomenological aspects of the inclusion of the multi-parton effects have been widely pursued since the early days of the parton model [5], while many formal features of the multiparton evolution [6,7] follow from a pioneering study on fragmentation functions [8]. Elements of the QCD formulation are presented in [9][10][11][12], while a 2 GPD interpretation is provided in [13]. Positivity bounds for dPDFs were recently discussed in [14].Dynamic modeling of dPDFs has been up to now very little explored. In fact, the only calculations which nonperturbatively predict the valence dPDFs of the nucleon have been carried out in the framework of the MIT bag model [15] and in the constituent quark model [16]. In both calculations, the approximate implementation of relativity has the side effect of producing distributions with unphysical support, i.e., extending outside the interval 0 ≤ x 1 , x 2 , x 1 + x 2 ≤ 1. An important aspect of the dPDFs is that they allow to detect tight diquark correlations, in which case the dPDF becomes a function of the sum x 1 + x 2 . In the (x 1 , x 2 ) distribution plot the fingerprints are the straight-line structures inclined by 45 • .The purpose of this talk is to explore the valence dPDF of the proton in a very simple model which includes the exact momentum-conservation constraints and complies to the Lorentz invariance. Satisfaction of these requirements results in proper theoretical features, such as the correct support in the x-variables or the quarknumber and momentum sum rules [17]. In the considered model, the correlations between partons are due Talk presended by WB at the Light Cone 2013 Conference, Skiathos, Greece.