Abstract: Successful fusion pilot plant (FPP) design hinges upon the ability to predict and control exhaust conditions to maximize the lifetime of plasma-facing components. This work describes a pathway toward high-fidelity, first-principles simulations with predictive capabilities for plasma particle fueling and detachment. It presents the coupling of a continuum full-f gyrokinetic turbulence model with atomic neutral models, using the Gkeyll code. To investigate how neutral interactions and plasma shaping fundamentally affect edge turbulent transport, we carry out simulations of DIII-D inner-wall-limited (IWL) plasmas. We specifically consider negative triangularity (NT) plasmas, which exhibit robust confinement properties without the presence of disruptive edge localized modes (ELMs). Results demonstrate good agreement with experimental data, with neutral interactions providing important particle fueling and heat loss channels. Our analysis reveals key differences in the shear flow and turbulent fluctuations that contribute to improved confinement properties in NT.