In the context of cosmology, a perfect fluid is generally described by an equation of state (EoS) of the following form: $$\tag{1} p = w \rho, $$ where $p$ is the fluid's pressure and $\rho$ is its energy density, while $w$ is a parameter (generally a simple constant). Non-linear relations are possible, for example the van der Wals EoS (the state parameter $w$ in (1) could be turned into a function of $\rho$): $$\tag{2} p(\rho) = \frac{w \rho}{1 - \alpha \rho} - \beta \rho^2. $$ A mixture of independant fluids could be defined by adding several densities and pressures: \begin{align} \rho &= \sum_i \rho_i, & p &= \sum_i p_i. \tag{3} \end{align} Now, I'm wondering if it's posible to mix the EoS of several superposed fluids by adding crossed terms. For example: \begin{align} p_1(\rho_1, \rho_2) &= w_1 \rho_1 + \alpha \rho_2^3, \tag{4} \\[2ex] p_2(\rho_1, \rho_2) &= w_2 \rho_2 + \beta \rho_1 + \gamma \rho_1^2. \tag{5} \end{align} I'll take radiation and matter and the following EoS as an example that puzzles me: \begin{align} p_{\text{mat}} &= \frac{\rho_{\text{mat}} \, \rho_{\text{rad}}}{3(\rho_{\text{mat}} + \rho_{\text{rad}})}, \tag{6} \\[1ex] p_{\text{rad}} &= \frac{\rho_{\text{rad}}^2}{3(\rho_{\text{mat}} + \rho_{\text{rad}})}. \tag{7} \end{align} Notice that \begin{align} p_{\text{mat}} + p_{\text{rad}} &= \frac{1}{3} \, \rho_{\text{rad}}, & \frac{p_{\text{rad}}}{p_{\text{mat}}} &= \frac{\rho_{\text{rad}}}{\rho_{\text{mat}}}. \tag{8} \end{align} Usually, radiation propagating in vacuum has the well known EoS $p_{\text{rad}} = \frac{1}{3} \, \rho_{\text{rad}}$, so (7) shouldn't be allowed. But the total energy-momentum tensor of the mixture has the following trace: $$\tag{9} \rho - 3 p = \rho_{\text{mat}} + \rho_{\text{rad}} - 3 p_{\text{mat}} - 3 p_{\text{rad}} = \rho_{\text{mat}}. $$ In this theoretical example, radiation is propagating in a matter fluid (it is superposed to it and is interacting with it!). Matter could react to the presence of radiation by absorbing or emitting some radiation, so produces a pressure when $\rho_{\text{rad}} \ne 0$, and radiation reacts to matter in a manner so its usual vacuum relation is modified.
So I have a few questions:
Since non-linear EoS are possible (relation (2) is an example), is it possible to mix several EoS as (4) and (5), for mixtures of interacting fluids, by adding crossed terms to their EoS? Since I never saw this before, I would like to see some specific examples from thermodynamics.
Are EoS (6) and (7) actually make any physical sense, especially for the radiation pressure which isn't in the usual $\frac{1}{3} \, \rho_{\text{rad}}$ form? How can we interpret or justify them in a better way than what I wrote?
This question is related to two other questions about the same subject:
Is it valid to add energy densities of *interacting* perfect fluids?