This article comprehensively examines entropy production in a peristaltic flow of Bingham-Papanastasiou fluid in an asymmetric wavy wall channel. The viscosity of the fluid is selected to depend on the temperature in addition to the shear rate. Furthermore, viscous dissipation is considered. The equations governing momentum and energy are adjusted to incorporate the existence of a nonlinear correlation between stress and strain rate. The energy dissipation caused by viscosity in the energy equation exhibits variation between Newtonian and non-Newtonian fluids due to their distinct rheological behaviors and response to shear forces. Newton’s law of viscosity states that the viscous dissipation term for Newtonian fluids follows a simple linear relationship between shear stress and shear rate. In contrast, non-Newtonian fluids exhibit more complex rheological behaviors, with the correlation between shear stress and shear rate described by various constitutive equations, depending on the specific type of non-Newtonian behavior. Consequently, the energy equation for non-Newtonian fluid flow problems becomes more intricate and nonlinear compared to the energy equation for Newtonian fluids, resulting in a more challenging analysis. Modeling of the equations is traditionally obtained using all steps involved in analyzing peristaltic flow. Using MATLAB's bvp5c solver, the obtained coupled set of differential equations is numerically solved. The research explores how the Bingham number and stress growth exponent parameters affect different dimensionless variables, including velocity, streamlines, temperature, axial pressure gradient, entropy production, and the Bejan number.