This study constructs a hybrid neural network model by integrating the physical constraints of the Bernoulli equation and Nikolaev's formula. The model is designed to explore and predict the variation pattern of the cold end temperature in a vortex tube. The input parameters include inlet pressure, inlet temperature, and cold mass fraction, with the cold end temperature as the output parameter. The network employs a multilayer feedforward model and the Levenberg-Marquardt learning algorithm, using a hyperbolic tangent function as the activation function. To evaluate the statistical validity of the developed model, the coefficient of determination (R²) and root mean square error (RMSE) are utilized, along with an analysis of the model's uncertainty and robustness. The hybrid model achieves an R² of 0.9936 and an RMSE of 0.3392, demonstrating strong performance in terms of uncertainty and robustness. These results indicate that the model accurately predicts the cold end temperature variation in the vortex tube. Furthermore, the findings reveal an optimal pressure range (0.49 MPa to 0.76 MPa) and cold mass fraction range (0.1 to 0.2) that minimize the cold end temperature.