The formation and evolution of the massive disk galaxies in the nearby universe (z < 0.1) represents an important open question in the present-day galactic astrophysics research. One class of such galaxies are the giant low surface brightness galaxies, e.g, Malin1 ,UGC1378, UGC 1382 etc, that are observed to show a complex morphology- a central high surface brightness stellar disk (HSB) surrounded by an extended low surface brightness (LSB) disk. The extended LSB disk could form due to external accretion or by secular evolution processes. In this talk, I will discuss our exploration with IllustrisTNG50 simulation data to identify and study various properties of such double-exponential, massive disk galaxies in the stellar mass range >= 1011 solar mass. The structural properties of such galaxies are obtained by 2D GALFIT modeling. The radial scale length of the LSB disks are found to lie in the range of ~10-30 kpc, in agreement with observations. The specific star-formation properties of the above double-disk galaxies are studied to understand their distribution from blue star-forming to red quenched region. Finally, we study such galaxies in the Baryonic Tully-Fisher relation. Our theoretical exploration will be potentially useful in exploring such galaxies in the upcoming deep observed data such as from LSST.