Studying the physical characteristics of the melt region in polymer-based laser powder bed fusion (PBF-LB/P) offers good comprehension that is essential for quality control, process optimization, understanding of material behavior, and preserving process consistency. Among the few numerical modelling investigations of PBF-LB/P, this study focuses on 3D multiphysics simulation of melt pool morphologies, beginning with powder bed preparation incorporating the actual particle size distribution using the discrete element method (DEM). A thermal fluid dynamic numerical model is then developed based on the finite volume method (FVM) using Flow-3D Weld to simulate PBF-LB while using polyamide 11 (PA11) as the primary material. The cooling after sintering plays a crucial role for the formation of the melt region due to the high viscosity of the polymer material effectively resulting in creeping flow. Here, we propose a model for both sintering as well as the subsequent cooling step, a two-step model capable of simulating the finished structure accurately with relatively low computational resources. The model's accuracy is verified through a mesh dependency analysis and its predictive power via validation against internal single-track experiments providing dimension data on the melt region width and depth. The simulation results allow for a detailed exploration of the effects of laser settings on melt region morphologies and the occurrence of defects. The model's scope is further extended to multi-track simulation, revealing the connection between surface roughness and laser settings. In general, this study provides significant contributions to the comprehension of PBF-LB/P methodologies which in turn can be used for further process improvement.