During the last decade, the relation between activity cycle periods with stellar parameters has received special attention. The construction of reliable registries of activity reveals that solar type stars exhibit activity cycles with periods from few years to decades and, in same cases, long and short activity cycles coexist suggesting that two dynamos could operate in these stars. In particular, Epsilon Eridani is an active young K2V star (0.8 Gyr), which exhibits a short and long-term chromospheric cycles of near 3 and 13-yr periods. Additionally, between 1985 and 1992, the star went through a broad activity minimum, similar to the solar Maunder Minimum-state. Motivated by these results, we found in Epsilon Eridani a great opportunity to test the dynamo theory. Based on the model developed in Sraibman & Minotti (2019), in this work we built a non linear axisymmetric dynamo for Epsilon Eridani. The time series of the simulated magnetic field components near the surface integrated in all the stellar disc exhibits both the long and short-activity cycles with periods similar to the ones detected from observations and also time intervals of low activity which could be associated to the broad Minimun. The short activity cycle associated to the magnetic reversal could be explained by the differential rotation, while the long cycle is associated to the meridional mass flows induced by the Lorentz force. In this way, we show that a single non-linear dynamo model derived from first principles with accurate stellar parameters could reproduce coexisting activity cycles.

A dynamo model is presented, based on a previously introduced kinematic model, in which the reaction of the magnetic field on the mass flow through the Lorentz force is included. Given the base mass flow corresponding to the case with no magnetic field, and assuming that the modification of this flow due to the Lorentz force can be treated as a perturbation, a complete model of the large-scale magnetic field dynamics can be obtained. The input needed consists in the large-scale meridional and zonal flows, the small-scale magnetic diffusivity, and a constant parameter entering the expression of the $\alpha$-effect. When applied to a solar-like star, the model shows a realistic dynamics of the magnetic field, including cycle duration, consistent field amplitudes with the correct parity, progression of the zonal magnetic field towards the equator, and motion toward the poles of the radial field at high latitudes. Also, the radial and zonal components show a correct phase relation, and, at the surface level, the magnetic helicity is predominantly negative in the northern hemisphere and positive in the southern hemisphere.