In this paper, we investigate the nonlinear characteristics of the flow in a two-dimensional and a three-dimensional counterflow model with symmetrical structures. Through numerical simulations, we obtain the velocity fields of the fluid flow in these models for different Re. The numerical results are analyzed to understand the nonlinear characteristics and differences between the two-dimensional and three-dimensional models. The findings indicate that, when Re varies, both the two-dimensional and three-dimensional models exhibit solution bifurcations and nonlinear phenomena such as symmetry breaking, self-sustained oscillations, and chaos. As Re increases, the two-dimensional counterflow model displays a unique solution, an asymmetric solution, and an oscillating solution. Specifically, when Re < 4320, both the laminar and turbulent models show a unique, symmetric, and steady-state velocity distribution. For 4652 < Re < 8639, the two-dimensional model solutions are not unique, presenting a pair of antisymmetric, asymmetric solutions that nevertheless remain steady-state. When Re > 8639, the solution becomes oscillatory and unsteady. The three-dimensional counterflow model exhibits a two-dimensional solution independent of the Z-axis. At Re = 4652, both the three-dimensional and two-dimensional models produce the same unique, symmetric, and steady-state velocity distribution with no three-dimensional flow (W = 0). At Re = 8639, the three-dimensional model solutions are not unique, showing a pair of antisymmetric, asymmetric solutions, while still being steady and time-independent. At Re = 87,627, the three-dimensional model solution becomes oscillatory and unsteady. By elucidating the flow characteristics and nonlinear features of both models, this study compares the differences between the two-dimensional and three-dimensional flows, thereby laying the groundwork for simplification of the problem and further theoretical research.