In recent years, the demand to provide best-estimate predictions with confident bounds is increasing for the nuclear reactor performance and safety analysis. The Organisation for Economic Co-operation and Development Nuclear Energy Agency has been developing an international benchmark of the light water reactor (LWR) uncertainty analysis in modeling (UAM) for the examination of uncertainty quantification and propagation methodologies with various modeling and simulation code systems. The objective of the present work is to develop an uncertainty propagation mechanism based on the stochastic sampling method by taking into account the uncertainties of both basic nuclear data and fuel modeling parameters in the simulation of pressurized water reactors (PWRs) that can be incorporated in the conventional LWR simulation approach. More specifically, the Three Mile Island Unit 1–related exercises from the LWR-UAM benchmark have been modeled using the coupled TRACE/PARCS code system in the three-dimensional core representation. The input uncertainties of the neutronics simulation include few-group cross sections and kinetics parameters generated using the Sampler/Polaris sequence of SCALE 6.2.1. Several heat transfer–related variables for the fuel modeling were considered as sources of input uncertainty of the thermal-hydraulics simulations, including the thermal conductivity of fuel and cladding, fuel heat capacity, and the gap conductance. Dakota was used to sample input parameters of the coupled code system and to perform the uncertainty analysis. Two types of simulations were conducted: steady-state calculation at hot full-power condition and transient scenario initiated by the spatially asymmetric rod ejection accident. Quantities of interest for the steady-state calculation, including core multiplication factor and power peaking factors, were calculated with associated uncertainties. For transient calculations, best-estimate plus uncertainty results of the time evolution of core reactivity, core power, and peak fuel temperature were generated and analyzed. The Wilks’ formula was used to determine the necessary sample size to achieve a 95% confidence of the 95% limit of output quantities of interest. Although the uncertainty propagation and quantification method presented in this paper was developed for PWRs, it could be in general applicable to the multiphysics uncertainty quantification of other types of LWR cores.