Given a partition of a system into two smaller systems, the energy $U$ is devided into $U_1$ and $U_2$, with $$U=\mathcal{P}(U_1,U_2):=U_1+U_2,$$ so that $U_2$ is given by $U-U_1$. Here the operation $\mathcal{P}$ denotes the partition rule.
What I've written above is basically the extensivity of energy and it somewhat enforces the idea of packages of energy, or particles. The idea of extensivity is closely related to additivity. A small energy flow from system one to system two is then given via
$$U_1\rightarrow \mathbb{fl}^-(U_1):= U_1-\epsilon, \\ U_2\rightarrow \mathbb{fl}^+(U_2):=U_2+\epsilon.$$
What we don't want to drop is energy conservation, so that
$$U=\mathcal{P}(U_1,U_2)=\mathcal{P}(\mathbb{fl}^-(U_1),\mathbb{fl}^+(U_2)),$$
where $\mathbb{fl}^\pm$ define what a flow of energy is supposed to mean, and due to energy conservation these functions in part determined by $\mathcal{P}$.
Are there any good arguments why in thermodynamics or in classical statistical physics, the implicit laws for this kind of divistion should be formulated with a plus sign (other than the fact that there are good theories in which this idea works fulfill)?
After all, there are some models for which the notions of extensivity or additivity are not so clear for all quantities (e.g. some of the entropy-definitions of the last decades).
The relations $$\mathcal{P}(U_2,U_1)=\mathcal{P}(U_1,U_2),$$ $$\mathcal{P}(\mathbb{fl}^-(U_1),\mathbb{fl}^+(U_2))=\mathcal{P}(\mathbb{fl}^+(U_1),\mathbb{fl}^-(U_2)),$$ also seem very natural - although this might not be totally necessary, as the values of $U_1$ and $U_2$ are elements of the reals and therefore ordered/distinguishable.
Also, I think the order of putting subsystems together can't matter
$$\mathcal{P}(U_1,\mathcal{P}(U_2,U_3))=\mathcal{P}(\mathcal{P}(U_1,U_2),U_3).$$
There should be more of these kind of restrictions.
So for example, a first ansatz would be that one could consider $$U=\mathcal{P}(U_1,U_2):=U_1U_2,$$ so that $$\mathbb{fl}^-(U_1)\equiv\tfrac{U_2}{\mathbb{fl}^+(U_2)}U_1.$$
One realization of that (which is constructed such that the change via $\mathbb{fl}^-$ is the same as in the additive case) would be
$$U_1\rightarrow \mathbb{fl}^-(U_1):= \left(1-\tfrac{\epsilon}{U_1}\right)U_1 =U_1-\epsilon, \\ U_2\rightarrow \mathbb{fl}^+(U_2):= \left(1-\tfrac{\epsilon}{U_1}\right)^{-1}U_2 =U_2+\left(\tfrac{U_2}{U_1}\right)\epsilon+ O\left((\tfrac{\epsilon}{U_1})^2\right).$$
