Dissipative Colored Fluid Dynamics

• Dissipative colored fluid refers to a class of models combining explicit dissipation with a 'color' variable that can represent tracers, noise, gauge charge, phase labels, or pigment.
• Key methodologies include inertial corrections in ocean eddy flows, colored noise in dissipative particle dynamics, and non-Abelian reductions in relativistic hydrodynamics.
• These models provide practical insights for interpreting phenomena across scales—from coherent Lagrangian eddies and computational fluid interfaces to multi-fluid dynamics in neutron-star physics.

Searching arXiv for papers and usages of “dissipative colored fluid”. “Dissipative colored fluid” is not a single canonical object but a family of technically distinct constructions whose common feature is the coexistence of a fluid description, an explicit dissipation mechanism, and a notion of “color.” In the arXiv literature assembled here, “color” may denote buoyant colored matter or tracers in geophysical flow, time-correlated stochastic forcing in dissipative particle dynamics, non-Abelian color charge in relativistic hydrodynamics, pigment or sediment content in computational fluid models, phase labels in ternary diffuse-interface mixtures, or the QCD origin of color–flavor-locked matter in neutron-star cores (Beron-Vera et al., 2014, Borreguero et al., 2019, Torrente-Lujan, 22 May 2026, Jaffer, 2017, Li et al., 2024, Bhopalam et al., 2024, Haskell et al., 2012).

1. Terminological scope

A recurring source of confusion is that “colored” does not have a uniform meaning across these works. The term attaches either to the advected material, to the noise process, to internal gauge charge, or to phase labels in a multiphase formulation.

Usage	Meaning of “color”	Dissipative mechanism
Ocean eddies	floating matter, drifters, \emph{Sargassum}	Stokes drag and Coriolis deflection
DPD with colored noise	time-correlated random forcing	dissipative DPD friction with KP noise
Non-Abelian hydrodynamics	adjoint color charge	first-order shear, bulk-like, vector-dissipative terms
Marbling and ink simulation	pigment or sediment content	viscosity and particle–fluid drag
Ternary phase-field FSI	three phase “colors” $c_i$	viscous and Cahn–Hilliard dissipation
Neutron-star multi-fluids	CFL microphysical origin	shear, bulk, drag, reactions, mutual friction

The literature therefore does not support a single universal definition. A precise reading requires identifying what carries the color variable and which irreversible processes enter the constitutive or dynamical equations.

2. Inertial colored matter near coherent ocean eddies

In rotating two-dimensional incompressible ocean surface flow, the relevant conservative reference model is the geostrophic velocity

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$

with vorticity

$$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$

Within this setting, coherent Lagrangian eddies are materially coherent vortices detected by the Haller–Beron-Vera geodesic eddy framework: they are time-evolving closed material curves that resist filamentation and stretching, and in the cited applications they are obtained from 90-day forward integrations of altimetry-derived velocities, typically with $\lambda=1$, indicating exceptional coherence (Beron-Vera et al., 2014).

The central paradox addressed in this formulation is that materially coherent eddies generate conservative-looking fluid motion, whereas observed drifting buoys, floats, and floating vegetation often exhibit dissipative-looking clustering, filamentation, and apparent convergence or expulsion near eddies. The resolution is to abandon neutral-fluid advection for finite-size inertial particle dynamics. In the reduced slow-manifold regime,

$\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$

where $\delta=\rho/\rho^{\mathrm p}$ and $\tau=\frac{2a^2}{9\nu\delta}$. The induced inertial velocity has nonzero divergence,

$$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$

so the effective dynamics is compressible even though the underlying fluid is incompressible (Beron-Vera et al., 2014).

The flux criterion across a coherent eddy boundary $\gamma$ is

$$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$

Its sign gives a polarity–buoyancy selection rule. Cyclonic eddies, for which $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$0, attract positively buoyant particles $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$1 and repel negatively buoyant particles $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$2. Anticyclonic eddies, for which $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$3, attract negatively buoyant particles and repel positively buoyant particles. The paper states the result explicitly: anticyclonic coherent Lagrangian eddies attract negatively buoyant finite-size particles and repel positively buoyant particles, while cyclonic coherent Lagrangian eddies attract positively buoyant finite-size particles and repel negatively buoyant particles (Beron-Vera et al., 2014).

This framework is linked to several observations. Two RAFOS floats inside an anticyclonic eddy behaved differently after slight vertical motion: the slightly lighter float spiraled outward and left the eddy, whereas the slightly heavier float remained trapped and spiraled inward. Surface drifters near a Loop Current ring and an Argos-tracked drifter inside the Agulhas ring “Astrid” dispersed away from anticyclonic coherent rings, consistent with repulsion of positively buoyant surface instruments. By contrast, MERIS ocean color showed spiral accumulation of floating \emph{Sargassum} inside a cyclonic Gulf Stream ring, aligned with backward-time attracting light inertial LCS computed from the inertial equation (Beron-Vera et al., 2014).

A common misconception is to interpret these patterns as evidence of fluid dissipation. The paper argues instead that dye or an ideal passive scalar would remain conservative-looking with respect to material eddy boundaries, whereas drifters, debris, or \emph{Sargassum} behave as inertial particles. In this sense, the “colored fluid” appearance is a property of the advected matter, not of the incompressible carrier flow (Beron-Vera et al., 2014).

3. Dissipative particle dynamics with colored noise

In mesoscopic simulation, the phrase denotes a DPD system in which the stochastic forcing is time-correlated rather than white. Standard DPD decomposes the pairwise force as

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$4

with dissipative and random forces

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$5

and classical fluctuation–dissipation conditions

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$6

The colored-noise variant replaces the classical pairwise Wiener noise by a single-particle Kangaroo-process force with algebraic correlation (Borreguero et al., 2019).

The Kangaroo-process autocorrelation is

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$7

and in two dimensions the stochastic force on particle $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$8 is

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$9

with $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$0 uniformly distributed in $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$1. In the simulations reported for colored-noise DPD, the conservative force is set to zero, so the equation of motion becomes

$$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$2

This is a structural change relative to classical DPD because the random forcing is no longer pairwise, finite-ranged, or delta-correlated (Borreguero et al., 2019).

The microscopic statistics differ substantially from white-noise DPD. The velocity pdf remains Gaussian, but the acceleration pdf becomes bimodal; in cooling regimes such as $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$3, $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$4, it develops pronounced non-Gaussian tails. The VACF is fit by

$$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$5

with reported exponents $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$6, $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$7, and $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$8 for C-DPD cases and $$\omega=-\nabla\cdot v^\perp=\frac{g}{f}\,\nabla^2\eta.$$9 for a representative classical DPD case. The diffusion coefficient

$\lambda=1$0

scales with friction qualitatively similarly in C-DPD and standard DPD, so altered microstatistics do not imply qualitatively altered large-scale diffusion for stationary parameter choices (Borreguero et al., 2019).

The paper emphasizes that the classical pairwise fluctuation–dissipation theorem does not strictly apply in this setting, yet equilibrium-like stationary states are still observed. An empirical calibration between noise amplitude and dissipation is reported: $\lambda=1$1 For the setup tested, the system reaches stationary temperature plateaus, and distinct initial conditions relax to the same stationary velocity and acceleration statistics. Reported numerical parameters include a $\lambda=1$2 two-dimensional box, $\lambda=1$3, density $\lambda=1$4, $\lambda=1$5, $\lambda=1$6, and a baseline colored-noise equilibrium case with $\lambda=1$7, $\lambda=1$8, $\lambda=1$9 (Borreguero et al., 2019).

Within this literature, “dissipative colored fluid” therefore means a mesoscopic DPD fluid in which dissipation is balanced by algebraically correlated noise rather than white noise. The color attribute refers to the temporal color of the stochastic process, not to pigment or gauge charge (Borreguero et al., 2019).

4. Non-Abelian dissipative colored hydrodynamics

A distinct relativistic usage constructs a $\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$0-dimensional dissipative colored fluid by Scherk–Schwarz reduction of a neutral viscous conformal fluid in $\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$1 dimensions on an $\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$2-dimensional unimodular group manifold. The higher-dimensional off-diagonal stress components descend to non-Abelian color currents, while the parent shear tensor produces three classes of first-order structures in the reduced theory: shear, bulk-like, and vector-dissipative terms (Torrente-Lujan, 22 May 2026).

The internal geometry is encoded by left-invariant one-forms $\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$3 satisfying

$\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$4

with unimodularity condition

$\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$5

The metric ansatz is

$\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$6

so the gauge fields $\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$7 arise from off-diagonal metric components. In the reduced theory the perfect color current is

$\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$8

and the sound speed, at fixed $\dot x=v^{\mathrm p}=v+\tau(\delta-1)f\,v^\perp,$9 and scalar moduli, is

$\delta=\rho/\rho^{\mathrm p}$0

For a conformal parent fluid with $\delta=\rho/\rho^{\mathrm p}$1, the reduced equation of state becomes

$\delta=\rho/\rho^{\mathrm p}$2

so conformality is generically broken by reduction (Torrente-Lujan, 22 May 2026).

The transport map is one of the main results: $\delta=\rho/\rho^{\mathrm p}$3 Here $\delta=\rho/\rho^{\mathrm p}$4 is the reduced shear viscosity, $\delta=\rho/\rho^{\mathrm p}$5 multiplies the isotropic expansion term and acts as an emergent bulk viscosity, and $\delta=\rho/\rho^{\mathrm p}$6 controls a transverse vector-dissipative or color-conductivity-like response. The dissipative correction to the stress tensor is written as

$\delta=\rho/\rho^{\mathrm p}$7

with

$\delta=\rho/\rho^{\mathrm p}$8

No independent first-order color-diffusion coefficient appears beyond the parent viscosity $\delta=\rho/\rho^{\mathrm p}$9; the viscous correction to the color current is fixed by the reduction map (Torrente-Lujan, 22 May 2026).

Two subtleties are stressed. First, the reduction induces a hydrodynamic-frame issue: even if the parent theory is in Landau frame, the reduced stress tensor is generally not. Second, the internal rapidity $\tau=\frac{2a^2}{9\nu\delta}$0 may be treated as constant, as a prescribed slowly varying background, or as an additional hydrodynamic variable, but the last option requires an extra equation of motion or an extended first law (Torrente-Lujan, 22 May 2026).

The second law descends cleanly only under unimodular reduction. The parent entropy production

$\tau=\frac{2a^2}{9\nu\delta}$1

reduces to

$\tau=\frac{2a^2}{9\nu\delta}$2

For non-unimodular groups, uncontrolled source terms appear and positivity is no longer automatic. The paper therefore treats the construction as a toy model for non-Abelian dissipative hydrodynamics with potential phenomenological relevance to quark–gluon plasma (Torrente-Lujan, 22 May 2026).

5. Computational and interfacial colored fluids

In computational and applied settings, “colored fluid” usually refers to pigment-bearing or phase-labeled flow, and dissipation enters through viscosity, drag, or phase-field diffusion.

In paint marbling, slow motions of a thin colored paint layer floating on a liquid bath are modeled as a two-dimensional incompressible Newtonian fluid. The relevant Reynolds number is

$\tau=\frac{2a^2}{9\nu\delta}$3

and for $\tau=\frac{2a^2}{9\nu\delta}$4, $\tau=\frac{2a^2}{9\nu\delta}$5, and $\tau=\frac{2a^2}{9\nu\delta}$6, the paper reports $\tau=\frac{2a^2}{9\nu\delta}$7. A characteristic decay length

$\tau=\frac{2a^2}{9\nu\delta}$8

controls the exponential decay of the Oseen-inspired short-stroke velocity field,

$\tau=\frac{2a^2}{9\nu\delta}$9

The corresponding stream function is

$$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$0

Closed-form homeomorphisms are then composed to generate marbling patterns. In this model, colors are passively advected tracers; no advection–diffusion equation is solved, and viscous dissipation shapes the velocity field rather than directly blurring color edges (Jaffer, 2017).

A more recent graphics formulation treats ink as a particle-laden flow using particle flow maps. The fluid obeys

$$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$1

while sediment particles satisfy

$$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$2

The visible color is carried by the pigment or sediment phase, with rendering based on the sediment volume fraction $$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$3; the paper states that no separate scalar advection–diffusion equation is solved and that in the reported results $$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$4. Dissipative forces are incorporated into flow maps through a covector path integral, and fluid–sediment coupling is enforced by the variable-coefficient Poisson problem

$$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$5

The reported phenomena include bulging and breakup of suspension drop tails, torus formation, torus disintegration, coalescence of sedimenting drops, vortex bulbs, viscous tails, fractal branching, and hierarchical structures (Li et al., 2024).

A third usage appears in diffuse-interface fluid–structure interaction with three immiscible incompressible fluids and an elastic solid. The phase or “color” variables are $$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$6, $$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$7, with $$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$8. The fluid equations are

$$\nabla\cdot v^{\mathrm p}=-\tau(\delta-1)f\,\omega,$$9

with

$\gamma$0

and the phase fields evolve by

$\gamma$1

The total energy combines fluid kinetic energy, Ginzburg–Landau mixing energy, fluid–solid surface energy, and solid kinetic plus elastic energy. Under static wetting and in the absence of line work, the paper proves the dissipation law

$\gamma$2

Here the colored-fluid attribute is literal: the three fluids are represented by color or phase fields, and dissipation arises from viscosity and compositional diffusion (Bhopalam et al., 2024).

These computational and interfacial models share a practical theme. Color is not an abstract label but a rendered, phased, or sedimented constituent, and the irreversible terms determine both physical transport and visual morphology.

6. Dissipative colored matter in neutron-star hydrodynamics

In neutron-star theory, the relevant formalism is Newtonian multi-fluid hydrodynamics with entrainment, reactions, viscosity, and other dissipative couplings. The paper develops the general framework for $\gamma$3 coupled constituents with densities $\gamma$4, mass densities $\gamma$5, velocities $\gamma$6, chemical potentials $\gamma$7, and relative velocities $\gamma$8. The energy functional is

$\gamma$9

and entrainment enters through

$$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$0

Dissipation is organized by Onsager symmetry, using affinities $$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$1, relative velocities, and velocity gradients as thermodynamic forces, with reaction rates, drag forces, and viscous stresses as fluxes (Haskell et al., 2012).

One application is to deconfined quark matter in the color–flavor-locked phase with a neutral $$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$2 condensate. The three-fluid model consists of the CFL condensate $$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$3, the neutral kaon condensate $$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$4, and the entropy/phonon gas $$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$5. The paper explicitly notes that although color enters through CFL pairing, the low-energy hydrodynamics is effectively color-neutral because color and flavor are locked. The momenta are

$$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$6

$$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$7

$$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$8

Because $$\Phi_\gamma=\tau(1-\delta)\,f\int_{U_\gamma}\omega\,\mathrm d^2x.$$9 in the co-moving background, the steady entropy Euler equation implies $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$00 in equilibrium (Haskell et al., 2012).

In the superfluid truncation relevant to CFL–$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$01 matter, the admissible dissipative sector contains a standard shear viscosity and six bulk viscosity coefficients: $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$02 Bulk viscosity is driven by weak processes involving kaons and is important when the weak rate resonates with the oscillation frequency. Shear viscosity is associated mainly with the normal entropy component, with phonon–phonon scattering cited as dominant at low temperature. If rotating superfluids form quantized vortices, the theory also admits the standard mutual-friction force

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$03

of which only the $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$04 term is dissipative (Haskell et al., 2012).

The entropy production is positive-definite and quadratic in the Onsager forces and fluxes, while the mode-damping rate is written as

$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$05

The formalism is then used to discuss the damping of $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$06- and $$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$07-modes in both hyperon matter and CFL–$$v(x,t)=\frac{g}{f}\,\nabla^\perp \eta(x,t), \qquad \dot x=v(x,t),$$08 matter. In this context, “dissipative colored fluid” refers not to a hydrodynamic color current but to a QCD-origin multi-fluid whose macroscopic irreversibility is governed by reactions, viscous terms, inter-fluid drag, entrainment, and mutual friction (Haskell et al., 2012).

The astrophysical usage therefore sits at the opposite end of the spectrum from pigment or phase-field models: color originates in microscopic strong-interaction physics, yet the effective hydrodynamics is formulated in terms of coupled neutral or effectively color-neutral fluids.