In this study, a coupled model integrating flow, temperature, phase separation, fibre alignment, and wall-slip has been developed to elucidate the complex behaviour observed during high moisture extrusion (HME) fibre formation. By departing from previous high-resolution approaches, the model uses a mean-field simplification to conveniently address wall-slip, thus avoiding the numerical intractability associated with resolving microscopic phases through solving the full Cahn-Hilliard equations. The critical simulation parameters are justified through prior studies and microscopy data and may to a certain extent be quantifiable from dead-stop experiments. The model can capture key qualitative features of HME, including the spatial distribution of fibres in the cooling die and their orientation, as observed in microscopy. Moreover, the model explains a potential delicate interplay between die cooling, phase separation/syneresis and protein melt flow characteristics. The study identifies extensional and pre-cooling die orientation of fibres as promising avenues for future model refinement