This study investigates the mixing behavior in a Y-junction mixer for supercritical water hydrolysis using large eddy simulation with a discrete phase model. Yield changes were simulated using a two-step reaction model with first-order kinetics, based on particles’ temporal temperature data. Effective mixing produced closely matched mass and particle flow temperature distributions, both exhibiting bell-shaped profiles near the mixed temperature. Although variations in flow rate within ±25 % and changes in the inlet temperatures of supercritical water from 350 °C to 430 °C and subcritical water from 100 °C to 170 °C did not significantly affect the overall mixing performance, they did alter the mixed temperature and, subsequently, yield changes. Additionally, backflow occurred when Richardson number for the subcritical inlet reached approximately 7. In effective mixing, simulated yields were approximately 15 % lower than the ideal theoretical yields, calculated using the reaction rate constant at the mixed temperature.