The total flux (\$\Phi\$) through an solenoidal inductor of length \$l\$ and \$N\$ turns is proportional to the current through the inductor and the inductance \$L\$ of the inductor according to
$$\Phi =L \cdot I $$ $$\Rightarrow \Phi =\frac{\pi r^2\mu_0 N^2 }{l}\cdot I \tag{1}$$
Clearly, in this case, if we double the amount of turns (whilst holding the length and current constant) we increase the total flux by a factor of 4. This is because the inductance of a solenoid depends quadratically on \$N\$.
Now I have just learnt about the concept of reluctance (denoted \$R\$) and magnetomotive force (denoted \$m.m.f\$) where \$m.m.f\$ is defined to be \$m.m.f \equiv N\cdot I\$. These two quantities are related to each other by the total flux via the formula
$$m.m.f = R\cdot \Phi \tag{2}$$
Now applying the above formula (eq 2) to the case of the solenoidal inductor, if we hold the current constant but double the number of turns (N) of the inductor, the \$m.m.f\$ only doubles (since \$m.m.f \equiv N\cdot I\$). According to equation (1), the total flux must quadruple since \$\Phi \propto N^2\$. Thus in order for equation (2) to remain valid, the reluctance must necessarily half. That is, in order for equations 1 and 2 to be simultaneously valid, the reluctance must be inversely proportional to the number of turns so that \$R\propto \frac{1}{N}\$. However, everywhere I look, I find that reluctance is always given by equations that are independent of \$N\$. For instance, wikipedia states that \$R=\frac{l}{\mu A}\$.
How can this be? How can equations 1 and 2 be simultaneously true whilst reluctance is independent of \$N\$?