Dilatational Dissipation in Fluid Dynamics

In the world of high-speed fluid dynamics, where air can be compressed like a spring, a subtle yet powerful phenomenon governs the flow of energy: dilatational dissipation. While often overlooked in low-speed scenarios, understanding this process is critical for accurately predicting the behavior of supersonic jets, re-entry vehicles, and even cosmic phenomena. Many conventional engineering models, developed for incompressible fluids, fail to account for this energy loss mechanism, leading to potentially dangerous design flaws. This article demystifies dilatational dissipation by breaking it down into its core components. The first part, "Principles and Mechanisms," will delve into the fundamental physics, distinguishing the "squeezing" motion of a fluid from its "swirling" motion and explaining how each contributes to energy loss. Subsequently, "Applications and Interdisciplinary Connections" will demonstrate the profound impact of this concept, from improving computational models for aerospace engineering to explaining the dynamics of astrophysical shock waves.

To truly grasp the nature of high-speed flight, we must look under the hood of a moving fluid and understand its inner workings. Imagine a substance that can simultaneously swirl like a galaxy and compress like a spring. This dual character is the key. The story of dilatational dissipation is the story of how a fluid pays the price for being compressed, a price measured in lost energy and generated heat.

Let’s begin with a simple picture. Any complex motion of a fluid can be thought of as a combination of two fundamental types of movement. The first is a swirling, tumbling, rotating motion—think of the beautiful eddies that form when you stir cream into your coffee. This is a motion that preserves volume; the fluid elements are sheared and spun around, but not squashed. In the language of physics, we call this the ​​solenoidal​​ component of the flow, because it is mathematically divergence-free.

The second type of motion is a squeezing and stretching, an expansion and compression—imagine pushing the plunger of a syringe filled with air. This motion changes the volume of the fluid elements. We call this the ​​dilatational​​ component, as it is associated with changes in density. This component is mathematically curl-free, meaning it lacks the local spinning quality of the solenoidal part.

The famous ​​Helmholtz decomposition​​ theorem tells us that any velocity field can be rigorously and uniquely split into these two parts: a solenoidal (swirling) soul and a dilatational (squeezing) soul. For a fluid we treat as incompressible, like water in your kitchen sink, the story is simple. The fluid is like a bag of un-squashable marbles; it can only swirl. Its dilatational soul is dormant. But for a compressible fluid, like the air rushing over a supersonic aircraft's wing, both souls are wide awake and interacting in a complex dance.

Motion is not free. When you stir your coffee, the swirling doesn't last forever; it gradually dies down. The energy of motion, or kinetic energy, is converted into heat, warming the coffee ever so slightly. This process is called ​​viscous dissipation​​, and it is the result of internal friction within the fluid.

If fluid motion has two souls, then it must pay two different prices for this motion. Viscous dissipation, it turns out, can also be split perfectly into two parts, one for each type of motion. The total rate of energy dissipation per unit volume, which we can call $\Phi$, is the sum of these two contributions:

Let's not be intimidated by the symbols. This equation tells a beautiful physical story.

The first term, $2\mu S_{ij}S_{ij}$, is the ​​solenoidal dissipation​​. It represents the energy lost due to the friction of swirling and shearing motions. The quantity $S_{ij}$ is the strain-rate tensor, which measures how fast the fluid elements are being deformed, and $\mu$ is the familiar ​​dynamic viscosity​​—the property that measures a fluid's resistance to shear, like the difference between honey and water. This is the only form of dissipation that exists in an incompressible fluid.

The second term, $\zeta \theta^2$, is the star of our show: the ​​dilatational dissipation​​. This is the energy lost due to the friction of volume changes. The quantity $\theta$ is the dilatation, representing the rate at which the fluid is expanding or compressing at a point. And $\zeta$ (zeta) is a less famous but equally important fluid property called the ​​bulk viscosity​​. It measures the fluid's internal resistance to being compressed or expanded. While shear viscosity $\mu$ is about resistance to changes in shape, bulk viscosity $\zeta$ is about resistance to changes in size. This term is zero unless the fluid is actively being compressed or expanded ($\theta \neq 0$).

When a flow becomes turbulent, it is a chaotic maelstrom of eddies of all sizes. In a compressible flow, this chaos includes not just swirling vortices but also fluctuating pockets of compression and expansion, like a frantic, disorganized version of the syringe plunger. In this violent environment, we must distinguish between two very different processes that involve dilatation.

First, we have the turbulent version of dilatational dissipation, often denoted $\varepsilon_d$. This is a truly dissipative, irreversible process. It is the component of viscous friction that specifically targets compressive motions, turning their kinetic energy directly and permanently into heat. Like all forms of friction, it is a one-way street dictated by the second law of thermodynamics. It is always a sink of kinetic energy.

Second, we have a more subtle and fascinating process called the ​​pressure-dilatation correlation​​, $\Pi_d = \overline{p' \theta'}$. This term represents the work done by fluctuating pressure fields ($p'$) on the fluctuating dilatation ($\theta'$). Unlike dissipation, this process is, in principle, reversible. Imagine a pocket of high pressure compressing a fluid element; this converts kinetic energy into internal energy. However, that parcel might then move into a low-pressure region and expand, getting a "kick" that converts internal energy back into kinetic energy.

To put it in an analogy, dilatational dissipation is like a leak in a bucket of water—the water (energy) is lost for good. Pressure-dilatation is like water sloshing back and forth between two connected buckets, one representing kinetic energy and the other internal energy. Depending on how the pressure and dilatation fluctuations are correlated, it can be a net source or a net sink of turbulent kinetic energy in the long run. In many high-speed flows featuring shock-like structures, high pressure is correlated with strong compression, making pressure-dilatation a significant net sink of kinetic energy.

What determines the balance between the swirling and squeezing motions? What decides how much energy is channeled into the dilatational modes to be dissipated? The conductor of this symphony is a dimensionless quantity called the ​​turbulent Mach number​​, $M_t$. It is defined as the ratio of the characteristic speed of the turbulent fluctuations to the local speed of sound: $M_t = \sqrt{2k}/a$, where $k$ is the turbulent kinetic energy and $a$ is the speed of sound.

Do not confuse this with the flight Mach number of an aircraft. The turbulent Mach number describes the compressibility of the turbulence itself.

At low $M_t$, the turbulence behaves as if it's incompressible. The swirls dominate, and the squeezes are negligible. But as $M_t$ increases, a remarkable thing happens. The nonlinear dynamics of the flow, primarily through the pressure field, begin to transfer energy from the powerful solenoidal modes into the dilatational modes. The chaotic swirls literally start generating sound waves and compressive fluctuations.

As a result, the fraction of kinetic energy residing in the dilatational component, $k^d$, grows. Consequently, the dilatational dissipation, $\varepsilon_d$, becomes an increasingly important sink for turbulent energy. Rigorous theoretical analysis and high-fidelity computer simulations reveal a beautifully simple scaling law: the dilatational dissipation rate, relative to the solenoidal rate, is proportional to the square of the turbulent Mach number.

This quadratic scaling is a cornerstone of modern turbulence theory. It tells us that as the internal Mach number of the turbulence doubles, the importance of dilatational dissipation quadruples.

This entire discussion is not merely academic; it has profound practical consequences for engineering. To design a supersonic jet or a reentry capsule, engineers rely on ​​Computational Fluid Dynamics (CFD)​​, which uses turbulence models to predict things like drag and heating.

The workhorse models of turbulence, such as the standard ​​k-ε model​​, were developed and calibrated for incompressible flows. Their "DNA" only contains information about solenoidal dissipation. They are structurally blind to the existence of dilatational dissipation and pressure-dilatation.

When these models are used for a high-speed flow with significant $M_t$, they sense that energy is disappearing faster than they expect. Not knowing about the separate channel of dilatational dissipation, they make a critical error: they blame the only mechanism they know, solenoidal dissipation, and artificially inflate its value. This is like a mechanic fixing a car's flat tire by over-revving the engine. The result is a fundamentally wrong prediction of the turbulence structure, often leading to an over-prediction of turbulent mixing and dissipation, which can translate into inaccurate designs.

The solution, born from the physics we've just discussed, is to introduce a ​​compressibility correction​​. Engineers explicitly add a new term to their models to account for dilatational dissipation. One of the most famous corrections, proposed by Sarkar, takes the form:

Here, $\varepsilon_s$ is the standard solenoidal dissipation from the baseline model, and the second term is the model for $\varepsilon_d$, perfectly capturing the $M_t^2$ scaling we discovered. This is a beautiful example of fundamental physics being translated directly into an engineering tool that helps us build safer and more efficient high-speed vehicles.

To complete our journey, we must ask one final, crucial question: is compressibility always important for the structure of turbulence whenever the flight speed is high? The surprising answer is no.

We must again distinguish between the flight Mach number, $M_{\infty}$, and the turbulent Mach number, $M_t$. An aircraft might be flying at Mach 5, but deep inside the boundary layer clinging to its wing, the turbulent eddies might be fluctuating at speeds much lower than the local speed of sound, resulting in a small $M_t$.

This insight was formalized by Mark Morkovin in his famous ​​Morkovin's Hypothesis​​. It states that if the turbulent Mach number $M_t$ is small (typically less than about 0.3), the direct effects of compressibility on the turbulence structure are negligible. In this regime, the dilatational dissipation and pressure-dilatation terms are of order $M_t^2$ and can be safely ignored. The turbulence, in its heart, behaves as if it were incompressible.

So, what is the dominant effect of compressibility in these high-speed flows? It is the large variation in the mean fluid properties—especially density and viscosity—due to aerodynamic heating. To handle this, modelers use a technique called ​​Favre (or density-weighted) averaging​​, which elegantly simplifies the equations. The incredible consequence of Morkovin's hypothesis is that once we use this averaging technique, we can apply our trusted incompressible turbulence models to a vast range of hypersonic flows, like predicting the heat transfer to a vehicle's skin.

This brings our story full circle. We began by separating motion into swirls and squeezes. We learned that each has a price, a form of dissipation. We saw how the turbulent Mach number conducts the symphony between them, and how engineers must teach their models to hear this music. And finally, we learned that sometimes, even in the roar of a hypersonic wind, the turbulence itself whispers in an incompressible tune, and the squeezes fade into silence.

Having explored the principles and mechanisms of dilatational dissipation, we now embark on a journey to see this fascinating concept in action. Like a master key, it unlocks our understanding of phenomena across a breathtaking range of scales, from the whisper-thin air screaming over a hypersonic wing to the silent, colossal shock waves that ripple through galaxies. We will see that this is not merely an esoteric correction term for specialists but a beautiful and unifying thread woven into the very fabric of fluid dynamics.

Perhaps the most immediate and visceral application of dilatational dissipation is in the realm of high-speed flight. When an aircraft punches through the air at supersonic speeds, the simple, smooth patterns of low-speed flow give way to a maelstrom of shock waves, expansions, and violent turbulence. To engineer vehicles that can withstand these conditions, we must be able to predict them, and it is here that our story begins.

Imagine a supersonic mixing layer, the turbulent battleground where two streams of air at different speeds clash. In an incompressible, low-speed flow, this layer grows by greedily entraining the surrounding air, a process driven by turbulent eddies. But as the Mach number increases, something remarkable happens: the mixing layer grows much more slowly. The turbulence seems to become sluggish, less effective. Why? A key reason is dilatational dissipation. The rapid compressions and expansions within the turbulent eddies themselves begin to drain energy from the turbulence, acting as a powerful brake. In our computer models, including a term for dilatational dissipation allows us to capture this very real physical suppression of mixing, which is crucial for designing everything from supersonic combustors to noise-reduction systems on jet engines.

The story gets even more critical when we consider the flow over a complete aircraft. Think of the flow separating from a surface, like air breaking away behind a steeply angled flap or a backward-facing step—a common feature in aerospace design. The point where this flow reattaches to the surface downstream dictates pressures, forces, and structural loads. Early turbulence models, blind to compressibility's braking effect, would often predict that this reattachment happens too soon. By accounting for dilatational dissipation, our simulations give a more realistic, longer reattachment length, better matching what we observe in experiments and allowing for safer, more efficient designs.

Now, let's raise the stakes to the brutal environment of atmospheric re-entry. A space capsule or hypersonic vehicle returning to Earth is subjected to extreme heating. Predicting this aerothermal load is a matter of life and death. Here, the role of dilatational dissipation becomes paramount. In these flows, powerful shock waves form, and the turbulence interacting with them is intensely compressed. Without a proper model for dilatational dissipation, simulations can wildly over-predict the levels of turbulence downstream of the shock. This "excess" turbulence, in turn, would lead to a catastrophic over-prediction of the heat transferred to the vehicle's surface. By correctly modeling how compression drains energy from the turbulence, dilatational dissipation terms effectively "tame" the simulated turbulence, bringing predictions of skin friction and heat flux into much better alignment with reality. It helps us avoid terrifying (and unphysical) "heat spikes" in our simulations and design thermal protection systems that work.

How, then, do we capture this intricate physics within a computer simulation? This is where we peer into the mind of the computational physicist. It’s an art form guided by rigorous principles.

The first principle is the most fundamental of all: conservation of energy. When dilatational dissipation removes energy from the turbulent velocity fluctuations, that energy doesn't simply vanish. It must go somewhere. Where? It is converted directly into thermal energy, heating the fluid. Therefore, any "sink" term we add to the equation for turbulent kinetic energy, $k$, to model this dissipation must appear as an equal and opposite "source" term in the equation for the fluid's internal energy. This beautiful and simple balance ensures our simulations respect the First Law of Thermodynamics, connecting the abstract world of turbulence statistics to the tangible property of temperature.

But the physics is richer still. Turbulence is not just a blob of random energy; it has structure and directionality, a property known as anisotropy. Imagine turbulence being squeezed by a shock wave. It doesn't just lose energy; it changes shape. The fluctuations might be flattened in one direction and elongated in another. To capture this, physicists employ more advanced tools like Reynolds Stress Models, which track the evolution of each component of the turbulent stress tensor, $R_{ij}$. Within this sophisticated framework, compressibility manifests in two distinct ways. First, there is ​​dilatational dissipation​​, an isotropic process that drains energy from all components, representing the "death" of turbulence into heat. Second, there is the ​​pressure-dilatation​​ term, which describes how the mean compression or expansion of the fluid directly deforms the shape of the turbulent eddies, transferring energy between components. Distinguishing these two effects is essential for accurately predicting complex phenomena like shock-boundary layer interactions, where the entire structure of turbulence is violently rearranged.

The creation of these model terms is itself a journey of discovery. Often, we don't have a perfect theory from first principles. Instead, we use physical reasoning. For instance, to model how compressibility affects the turbulence scale equation (the equation for $\varepsilon$), we can deduce the form of the correction term by imposing constraints: it must have the right physical dimensions, it must vanish in the incompressible limit, and its effect should saturate in extreme compressions. This line of reasoning leads to elegant mathematical forms that can be tested and calibrated against experiments and more detailed simulations, allowing us to systematically build more accurate and robust models.

This intimate knowledge of the physics can even guide the computational process itself. In a technique called goal-oriented mesh adaptation, we can tell the computer to focus its efforts on the regions where the most important physics is happening. If we want to accurately calculate the total dilatational dissipation in a flow, what better guide could there be than the dilatational dissipation itself? We can use its magnitude as an indicator, telling the simulation to automatically place more grid points in regions of strong compression, ensuring we capture these critical events with high fidelity. The physics, in a sense, tells the computer how to look at it.

The influence of dilatational dissipation extends far beyond our terrestrial machines, reaching into the silent depths of the cosmos and the subtle propagation of sound.

Consider a sound wave traveling through a turbulent patch of air. The wave is a dance of compression and expansion. As it passes through the turbulence, the turbulent eddies are repeatedly squeezed and stretched. This process activates dilatational dissipation, causing the eddies to extract energy from the sound wave, converting it into heat. The result is that the sound wave is damped; its energy is quite literally "eaten" by the turbulence. From a macroscopic perspective, the turbulence gives the fluid an effective "bulk viscosity," a property that quantifies its resistance to volume change. This phenomenon is a cornerstone of aeroacoustics, explaining how turbulence helps to muffle sound in everything from jet exhausts to the atmosphere itself.

Now, let us turn our gaze from the audible to the astronomical. The interstellar medium is not a tranquil void; it is a violent place, crisscrossed by colossal shock waves from supernova explosions and stellar winds. A shock front is an incredibly thin region where the bulk kinetic energy of a gas is violently and irreversibly converted into thermal energy, heating the gas from thousands to millions of degrees. The engine of this transformation is viscosity. Inside the shock, the fluid is compressed at an astronomical rate. This immense velocity divergence, $\nabla \cdot \boldsymbol{u}$, drives a powerful viscous dissipation that acts as the primary heating mechanism. This process, in its essence, is a dramatic and magnificent display of dilatational dissipation, providing the thermal energy that makes supernova remnants glow for millennia.

After this grand tour, it is natural to ask: where does this effect, which has proven so versatile, ultimately come from? To find the answer, we must return to the very foundations of fluid mechanics—the relationship between stress and strain in a fluid.

For a simple Newtonian fluid, we relate the viscous stress to the rate of fluid deformation. This relationship contains two constants: the shear viscosity, $\mu$, which governs resistance to shape change, and the bulk viscosity, $\zeta$, which governs resistance to volume change. The dilatational dissipation we have been discussing is directly proportional to this bulk viscosity, $\zeta$.

However, in the 19th century, Sir George Stokes made a convenient proposition. Based on arguments for a monatomic gas, he hypothesized that these two viscosities were related, such that $\zeta=0$. For over a century, ​​Stokes's hypothesis​​ has been a cornerstone of fluid dynamics, built into the standard Navier-Stokes equations taught to every student. By this assumption, dilatational dissipation is defined out of existence from the very start for a vast number of problems!

So why have we spent this chapter discussing it? Because in the complex world of turbulence and high-energy plasmas, Stokes's hypothesis breaks down. A turbulent flow, when viewed macroscopically, does not behave like a simple gas. It develops an effective bulk viscosity, even if the underlying fluid molecules obey Stokes's rule. The collective action of countless eddies interacting with a compressive strain gives rise to a resistance to volume change that the mean flow feels. The same is true for plasmas inside a shock wave.

Dilatational dissipation, therefore, is the ghost in the machine. It is the physical manifestation of the failure of a convenient, simplifying assumption. It reminds us that the richness of the universe is often found in the places where our simplest models break down, and that by exploring these frontiers, we uncover a deeper and more unified understanding of the world.