Constitutive Modelling of Korteweg Fluids Using Liu's Method

This paper employs Liu's method of multipliers to derive thermodynamically consistent constitutive models for Korteweg fluids, successfully incorporating temperature-dependent capillary effects and cross-coupling between mechanical and thermal phenomena while establishing a generalized Gibbs relation.

The Gist

Imagine you are trying to bake a perfect cake. You have a recipe (the laws of physics) that tells you how ingredients like flour and sugar (mass and energy) should behave. But sometimes, when you mix things, strange things happen at the edges—like the batter clinging to the bowl or forming a skin. In fluid dynamics, this "edge behavior" is called capillarity.

This paper is about creating a new, more accurate "recipe" for a special type of fluid called a Korteweg fluid. These are fluids where the "skin" or surface tension isn't just a sharp line; it's a fuzzy, gradual transition zone, like the mist between a cloud and the sky, rather than a hard wall.

Here is the breakdown of what the authors did, using simple analogies:

1. The Problem: The "Rulebook" vs. Reality

Scientists have long tried to write down the rules for how these fuzzy fluids move. The standard rulebook (thermodynamics) says that energy can't be created out of nothing and that disorder (entropy) must always increase or stay the same.

However, previous attempts to write the rules for these fuzzy fluids often felt like they were "cheating" the rulebook. They had to make special, arbitrary assumptions to make the math work. The authors of this paper wanted to see if they could derive the rules strictly from the fundamental laws without cheating, using a specific mathematical tool called Liu's Method.

Think of Liu's Method as a very strict referee in a game. The referee says, "You must follow these balance laws (mass, momentum, energy), and you must also follow the rule that entropy cannot decrease. If your proposed rules for the fluid break this, they are invalid."

2. The New Approach: A Better Way to Count

The authors applied this "strict referee" method to Korteweg fluids. They made a few clever changes to how they looked at the problem:

• The "Fuzzy" Ingredient: They realized that to describe the fuzzy edge, you can't just look at how much "stuff" (density) is there. You also have to look at how quickly that "stuff" is changing from one spot to the next (the density gradient). It's like not just counting the number of people in a room, but also measuring how crowded it is at the door versus the back.
• The "Spin" Factor: When fluids move, they can spin (like a tornado) or stretch (like taffy). Previous studies often ignored the spinning part to make the math easier. The authors kept the spinning part in their calculations. Surprisingly, this made the math simpler and revealed a hidden "helper" variable (a multiplier) that was previously hard to find.
• The Entropy Detective: Instead of guessing how "disorder" (entropy) flows, they let the strict referee (Liu's method) tell them exactly what the flow must look like. It turned out the flow of disorder is directly linked to the flow of heat and the movement of the fuzzy edges.

3. The Big Discoveries

By letting the math do the heavy lifting without forcing it, they found three major things:

• The Stress is Reversible: They confirmed that the special "Korteweg stresses" (the forces caused by the fuzzy edges) are like a perfect spring. If you push them, they push back. They exist even when the fluid is perfectly still (equilibrium). This confirms they are a fundamental part of the fluid's nature, not just a side effect of movement.
• Temperature Matters: They found that the "strength" of the fuzzy edge (the capillary coefficient) can change depending on the temperature. This is like saying the "stickiness" of the mist changes if you heat it up. This connects their work to recent microscopic theories (kinetic theory) that suggest this should happen.
• A New "Gibbs Relation": In thermodynamics, there is a famous equation (Gibbs relation) that links energy, heat, and pressure. The authors derived a new, expanded version of this equation. Their version includes a term for the "fuzziness" of the edge. It's like adding a new chapter to the rulebook that explains how the edge contributes to the fluid's total energy.

4. What This Means (According to the Paper)

The paper doesn't claim this will immediately cure diseases or build new engines. Instead, it claims to have fixed the theoretical foundation.

• Consistency: They proved that the rules for these fluids are perfectly consistent with the laws of thermodynamics.
• Flexibility: They showed that there are actually two slightly different ways to write the rules for how these fluids conduct heat and move (Case 1 and Case 2 in the paper), but both lead to the same physical outcome.
• The "Holographic" Property: They noted that for these fluids, the complex internal forces can sometimes be described as if they were coming from a single, simple "potential" (like a hill that the fluid rolls down). This connects fluid dynamics to deeper physics concepts, including quantum mechanics.

Summary

Think of this paper as a group of architects who went back to the blueprints of a complex building (Korteweg fluids). Previous architects had to use duct tape and guesswork to make the roof fit. These authors used a laser-level (Liu's method) and found that if you look at the building from a slightly different angle (keeping the "spin" and the "fuzzy edges" in mind), the roof fits perfectly on its own, and the building follows all the laws of physics naturally. They didn't just fix the roof; they also discovered a new, hidden room (the generalized Gibbs relation) that explains how the building's edges store energy.

Technical Summary

1. Problem Statement

The paper addresses the constitutive modelling of Korteweg fluids, which are fluids exhibiting capillary effects (surface tension) and described by a diffuse interface where the mass density gradient acts as an order parameter. The central challenge is ensuring thermodynamic consistency: deriving constitutive relations (stress tensor, heat flux, entropy) that are compatible with the Second Law of Thermodynamics (entropy balance law).

Previous approaches, such as the Coleman-Noll procedure or variational methods based on free energy functionals, often require specific assumptions about the structure of internal energy or entropy flux. The authors aim to resolve these issues using Liu's Method of Multipliers, a rigorous framework that derives constitutive constraints directly from the balance laws and the entropy inequality without ad hoc assumptions about the fluxes.

2. Methodology

The authors employ Liu's Method of Multipliers, a procedure that treats the entropy balance law as the primary equation and the conservation laws (mass, momentum, energy, and density gradient) as constraints.

Key Methodological Steps:

• Constitutive State Space: The state space is expanded to include the mass density gradient ($\nabla \rho$) and its higher derivatives, as well as the velocity gradient ($\mathbf{D}$), temperature ($\theta$), and their spatial derivatives.
• Specific Entropy Structure: A crucial assumption is that specific entropy $s$ depends on specific internal energy $e$ and the remaining state variables independently. This allows the recovery of classical thermodynamic relations while accommodating non-equilibrium contributions.
• Velocity Gradient Decomposition: Unlike standard approaches that often discard the skew-symmetric part of the velocity gradient (vorticity), the authors retain it. This decomposition into symmetric ($\mathbf{D}$) and skew-symmetric ($\mathbf{W}$) parts simplifies the derivation of multipliers related to capillary effects.
• Two-Step Derivation:
  1. Equilibrium Analysis: The equilibrium stress tensor and heat flux are derived by imposing the conditions that entropy production vanishes and is minimized in equilibrium.
  2. Non-Equilibrium Analysis: General constitutive relations for non-equilibrium fluxes are derived from the residual inequality, ensuring non-negative entropy production.

3. Key Contributions and Results

A. Derivation of Multipliers

The authors successfully determine the Lagrange multipliers associated with the balance laws:

• $\Lambda_e = 1/\theta$ (inverse temperature).
• $\Lambda_{\nabla \rho} = \frac{1}{\theta} \alpha(\rho, \theta) \nabla \rho$.
  • Significance: The multiplier $\alpha(\rho, \theta)$ is identified as the Korteweg coefficient. Crucially, the authors allow $\alpha$ to depend on temperature, a feature supported by recent kinetic-theory results (Enskog–Vlasov equation) but often neglected in macroscopic models.

B. Korteweg Stress Tensor

The equilibrium stress tensor ($t^{eq}_{ij}$) is derived as:
$$t^{eq}_{ij} = -\left[ p + \rho \frac{\partial \alpha}{\partial \rho} |\nabla \rho|^2 + \rho \alpha \Delta \rho \right] \delta_{ij} + \alpha \frac{\partial \rho}{\partial x_i} \frac{\partial \rho}{\partial x_j}$$

• Novelty: The derivation emphasizes that Korteweg stresses are equilibrium (reversible) stresses.
• Non-Uniqueness: The paper discusses the non-uniqueness of the stress tensor representation (e.g., presence or absence of the Hessian term $\nabla \otimes \nabla \rho$). It argues that different forms are divergence-equivalent and physically indistinguishable in the bulk, though they may differ at boundaries.

C. Entropy Flux and Cross-Coupling

The entropy flux $\varphi_j$ is determined from the residual inequality rather than assumed:
$$\varphi_j = \frac{1}{\theta} \left( q_j + \rho \alpha \nabla \rho \cdot \nabla \cdot \mathbf{v} \right)$$

• Cross-Coupling: The temperature dependence of $\alpha$ leads to a cross-coupling term in the non-equilibrium heat flux ($q^*$). Specifically, a temperature gradient can drive a mass flux and vice versa, a phenomenon absent if $\alpha$ is temperature-independent.
• Two cases for non-equilibrium fluxes are presented (entropic vs. energetic representations), both ensuring non-negative entropy production.

D. Generalized Gibbs Relation

A major theoretical outcome is the derivation of a generalized Gibbs relation as a direct consequence of the procedure (not an assumption):
$$\theta ds = de - \frac{p}{\rho^2} d\rho + \frac{\alpha(\rho, \theta)}{\rho} (\nabla \rho \cdot d(\nabla \rho))$$
This relation explicitly includes the differential of the density gradient, formally integrating capillary effects into the fundamental thermodynamic equation.

4. Significance and Implications

1. Thermodynamic Rigor: The study demonstrates that Liu's method can handle complex, weakly nonlocal fluids (Korteweg fluids) without requiring the specific internal energy to be pre-defined with capillary terms. The stress structure emerges naturally from the entropy inequality.
2. Temperature-Dependent Capillarity: By allowing the Korteweg coefficient $\alpha$ to depend on temperature, the model aligns with microscopic kinetic theory. This introduces thermo-mechanical coupling in the non-equilibrium fluxes, suggesting that temperature gradients can induce capillary-driven flows.
3. Unified Framework: The derivation of the generalized Gibbs relation and the holographic property of the stress tensor (where momentum balance reduces to a potential form) connects macroscopic fluid dynamics with quantum-mechanical fluid equations and variational approaches.
4. Future Applications: The methodology is presented as a robust foundation for modeling mixtures of Korteweg fluids and higher-grade fluids, where classical approaches often struggle with closure assumptions.

Conclusion

The paper provides a comprehensive and thermodynamically consistent framework for Korteweg fluids. By leveraging Liu's method and retaining specific structural details (like the skew-symmetric velocity gradient and temperature-dependent coefficients), the authors derive constitutive relations that recover classical results while extending them to include cross-coupling effects and a generalized Gibbs relation. This work bridges the gap between macroscopic continuum mechanics and microscopic kinetic theory for fluids with diffuse interfaces.