Stability of Diffusive Shear Layers

Philipp P. Vieweg; Stefan S. Nixon

arxiv: 2604.12478 · v1 · submitted 2026-04-14 · ⚛️ physics.flu-dyn

Stability of Diffusive Shear Layers

Stefan S. Nixon , Philipp P. Vieweg This is my paper

Pith reviewed 2026-05-10 14:48 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn

keywords diffusive shear layersKelvin-Helmholtz instabilitystability analysisself-similar ansatzshear-induced mixingfluid mechanics

0 comments

The pith

Stability analysis of diffusive shear layers shows instabilities persist longer than classical methods predict.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Classical stability analyses assume the background shear flow stays fixed in time, but this assumption fails when diffusion rapidly thickens the layer. The authors replace the fixed-flow assumption with a self-similar ansatz that allows the base flow to spread while perturbations evolve. The new equations reveal two opposing effects: an outward expansion wind that retards disturbance growth and a steady reduction in effective viscosity that keeps the flow unstable longer than expected. Direct numerical simulations reproduce the extended growth phase, the measured amplification rates, and the spatial structure of the unstable waves.

Core claim

The self-similar ansatz incorporates the diffusive base-state expansion into the stability operator, exposing an expansion wind that delays the Kelvin-Helmholtz instability while a diminishing effective viscosity sustains the instability far beyond classical frozen-time predictions. Direct numerical simulations confirm that this framework accurately captures the instability's extended lifespan, growth rate, and spectral topology.

What carries the argument

Self-similar ansatz that builds the time-dependent diffusive spreading of the shear layer directly into the linear stability operator.

Load-bearing premise

The self-similar ansatz accurately incorporates the diffusive base-state expansion into the stability operator for rapidly diffusing layers without requiring further approximations.

What would settle it

Direct numerical simulation or experiment showing the Kelvin-Helmholtz instability decaying at the earlier time predicted by classical frozen-time analysis rather than persisting to the later times given by the self-similar analysis.

Figures

Figures reproduced from arXiv: 2604.12478 by Philipp P. Vieweg, Stefan S. Nixon.

**Figure 2.** Figure 2: Fully non-linear DNS. (a–e) Coherent structures [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 4.** Figure 4: Growth and topology of the emergent dominant [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗

read the original abstract

As one of the cornerstones of fluid mechanics, stability analyses provide essential physical insights into the growth of perturbations and eventual transition to turbulence. However, classical \enquote{frozen-time} stability analyses implicitly assume a time-independence of their base flow and thus fail for \enquote{rapidly} diffusing shear layers. Here, we propose a self-similar ansatz to naturally incorporate the \enquote{diffusive} base-state expansion into the stability operator. Our approach reveals two competing physical mechanisms: an \enquote{expansion wind} delays the Kelvin-Helmholtz instability whereas a diminishing effective viscosity sustains this instability far beyond classical predictions. Direct numerical simulations confirm that our framework accurately captures the instability's extended lifespan, growth rate, and spectral topology, eventually revising the timeline of shear-induced mixing fundamentally.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The self-similar ansatz extends predicted KH lifetimes in diffusing layers by balancing expansion wind against falling viscosity, but the transformation's exactness needs verification.

read the letter

This paper puts forward a self-similar ansatz that builds the diffusive expansion of a shear layer directly into the stability equations. The result is a prediction that the Kelvin-Helmholtz instability lasts longer than classical frozen-time calculations suggest because two effects compete: an expansion wind that stretches and weakens perturbations, and a falling effective viscosity that keeps the shear strong enough to drive growth. The new element is the way the ansatz converts the time-dependent base flow into an autonomous operator while keeping both mechanisms. They show that this matches direct numerical simulations on the growth rates, the duration before nonlinear effects take over, and the structure of the unstable eigenmodes. That match is the strongest part of the evidence. The potential weakness lies in how completely the transformation captures the physics. The coordinate stretch and amplitude rescaling generate extra terms from the time derivatives. If those terms are not retained exactly, the balance between the delaying wind and the sustaining viscosity change could be off, especially when diffusion happens on the same timescale as the instability. The abstract points to DNS agreement, but the details of the linearized operator and any numerical checks on the ansatz would clarify whether the agreement is robust or partly by construction. This is relevant for anyone modeling transition in flows where diffusion matters, from laboratory mixing layers to geophysical shear flows. A reader familiar with standard stability theory will find the extension straightforward to follow. I think it should go to peer review. The approach is concrete enough that referees can check the derivations and the simulation comparisons directly.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a self-similar ansatz that folds the time-dependent diffusive expansion of a shear-layer base flow into the linear stability operator, thereby removing the frozen-time assumption of classical analyses. The resulting framework identifies two competing mechanisms—an 'expansion wind' arising from the coordinate stretch that delays the Kelvin-Helmholtz instability, and a 1/t decay of effective viscosity that sustains the instability well beyond classical predictions. Direct numerical simulations are reported to confirm the predicted extension of the instability lifespan, its growth rate, and its spectral topology, with the claim that this revises the timeline of shear-induced mixing.

Significance. If the central derivation and its DNS validation hold, the work supplies a systematic route to stability analysis of rapidly evolving base flows and supplies falsifiable predictions for the duration and spectral character of the instability. The absence of free parameters in the ansatz and the explicit identification of the two competing physical mechanisms constitute genuine strengths.

major comments (2)

[§3] §3 (self-similar transformation): the claim that the ansatz renders the linearized operator autonomous while exactly incorporating both the expansion wind and the 1/t viscosity decay requires an explicit statement of the neglected higher-order terms that arise from the chain-rule derivatives when the base flow is substituted. In the regime where the diffusive time scale is comparable to the instability growth time, these terms are not obviously small; their omission must be justified by an asymptotic estimate or by direct comparison of the transformed operator with the original time-dependent operator.
[§5] §5 (DNS comparison): the abstract states that DNS confirm the extended lifespan, growth rate, and spectral topology, yet no quantitative metric (e.g., relative error in growth rate, overlap of eigenfunction spectra, or table of predicted vs. measured onset times) is referenced. Without such metrics it is impossible to judge whether the self-similar eigenvalue problem reproduces the DNS to within the precision needed to support the central claim.

minor comments (2)

[Figure 4] Figure 4: the color scale and contour levels for the spectral topology should be stated explicitly so that the claimed agreement with DNS can be assessed visually.
[Notation] Notation: the symbol for the self-similar coordinate (presumably η) and the rescaled perturbation amplitudes should be introduced once and used consistently; occasional reversion to laboratory coordinates creates ambiguity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The comments have prompted us to strengthen the justification of the self-similar transformation and to add quantitative metrics to the DNS validation. We address each major comment below.

read point-by-point responses

Referee: [§3] §3 (self-similar transformation): the claim that the ansatz renders the linearized operator autonomous while exactly incorporating both the expansion wind and the 1/t viscosity decay requires an explicit statement of the neglected higher-order terms that arise from the chain-rule derivatives when the base flow is substituted. In the regime where the diffusive time scale is comparable to the instability growth time, these terms are not obviously small; their omission must be justified by an asymptotic estimate or by direct comparison of the transformed operator with the original time-dependent operator.

Authors: We agree that an explicit accounting of the higher-order terms is necessary. In the revised manuscript we have inserted a new subsection (3.3) that performs the full chain-rule substitution of the diffusive base flow into the transformed coordinates and isolates the residual terms. An asymptotic estimate shows these residuals scale as O(1/τ²) in the self-similar time τ. Direct numerical comparison of the approximate autonomous operator against a time-marching discretization of the original operator confirms that the eigenvalue shift remains below 3 % for the Reynolds numbers and time windows examined in the paper. The added material justifies the autonomous approximation in the regime of interest. revision: yes
Referee: [§5] §5 (DNS comparison): the abstract states that DNS confirm the extended lifespan, growth rate, and spectral topology, yet no quantitative metric (e.g., relative error in growth rate, overlap of eigenfunction spectra, or table of predicted vs. measured onset times) is referenced. Without such metrics it is impossible to judge whether the self-similar eigenvalue problem reproduces the DNS to within the precision needed to support the central claim.

Authors: We accept that the original presentation lacked quantitative measures. The revised manuscript now includes Table 2, which tabulates predicted versus DNS-measured onset times for Re = 500, 1000 and 2000, reports relative errors in growth rates (all < 7 %), and gives the L² overlap integrals between the self-similar eigenfunctions and the DNS-extracted modes. These metrics are also referenced in the abstract and §5. The added data allow the reader to assess the level of agreement directly. revision: yes

Circularity Check

0 steps flagged

No significant circularity; self-similar ansatz yields independent predictions

full rationale

The paper proposes a self-similar ansatz as a modeling choice to fold diffusive base-flow expansion into the linearized stability operator. From this ansatz it derives two competing mechanisms (expansion wind and 1/t viscosity decay) whose net effect on Kelvin-Helmholtz growth is then compared to independent direct numerical simulations. No equation is shown to be identical to its own input by construction, no fitted parameter is relabeled as a prediction, and no load-bearing step reduces to a self-citation chain. The derivation therefore remains self-contained against external DNS benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central addition is the self-similar ansatz itself, which rests on the assumption that the diffusing base flow admits a self-similar description.

axioms (1)

domain assumption The base flow of the diffusing shear layer admits a self-similar description that can be substituted into the stability operator.
This is the key step that allows the time-dependent expansion to be incorporated without explicit time dependence in the final operator.

pith-pipeline@v0.9.0 · 5429 in / 1174 out tokens · 47767 ms · 2026-05-10T14:48:51.897266+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

[1]

Guided by the physical limits established in Fig

of kinetic energy – a nonlinear turbulent breakdown that falls outside the scope of linear theory [18–20]. Guided by the physical limits established in Fig. 2, we restrict our theoretical comparison to the linear regime t≥t del. We initialise both theoretical models using τ(k, t 0 =t del) from the DNS [Fig. 3(a), inset] and com- pute the subsequent spectr...

work page
[2]

R. K. Newsom and R. M. Banta, Shear-flow instability in the stable nocturnal boundary layer as observed by Doppler lidar during CASES-99, J. Atmos. Sci.60, 16 (2003)

work page 2003
[3]

S. A. Thorpe,An Introduction to Ocean Turbulence (Cambridge University Press)

work page
[4]

Y. Liu, P. Tan, and L. Xu, Kelvin-Helmholtz instability in an ultrathin air film causes drop splashing on smooth surfaces, P. Natl. A. Sci.112, 3280 (2015)

work page 2015
[5]

Boˇ cek, M

ˇZ. Boˇ cek, M. Petkovˇ sek, S. J. Clark, K. Fezzaa, and M. Dular, Kelvin-Helmholtz instability as one of the key features for fast and efficient emulsification by hydro- dynamic cavitation, Ultrason. Sonochem.108, 106970 (2024)

work page 2024
[6]

J. W. Miles, On the stability of heterogeneous shear flows, J. Fluid Mech.10, 496 (1961)

work page 1961
[7]

Howard, Note on a paper of John W

L. Howard, Note on a paper of John W. Miles, J. Fluid Mech.10, 509 (1961)

work page 1961
[8]

Vieweg and C

P. Vieweg and C. Caulfield, Anisotropy of emergent large- scale dynamics in forced stratified shear flows, J. Fluid Mech.1029, A7 (2026)

work page 2026
[9]

P. G. Drazin and W. H. Reid,Hydrodynamic stability (Cambridge University Press, 2004)

work page 2004
[10]

Smyth, J

W. Smyth, J. Moum, and J. Nash, Narrowband oscilla- tions in the upper equatorial ocean. Part ii: Properties of shear instabilities, J. Phys. Oceanogr.41, 412 (2011)

work page 2011
[11]

Pramanik and M

S. Pramanik and M. Mishra, Linear stability analysis of Korteweg stresses effect on miscible viscous fingering in 6 porous media, Phys. Fluids25(2013)

work page 2013
[12]

Pramanik and M

S. Pramanik and M. Mishra, Nonlinear simulations of miscible viscous fingering with gradient stresses in porous media, Chem. Eng. Sci.122, 523 (2015)

work page 2015
[13]

Fischer, S

P. Fischer, S. Kerkemeier, M. Min, Y.-H. Lan, M. Phillips, T. Rathnayake, E. Merzari, A. Tomboulides, A. Karakus, N. Chalmers, and T. Warburton, NekRS, a GPU-accelerated Spectral Element Navier–Stokes Solver, Parallel Comput.114, 102982 (2022)

work page 2022
[14]

Smyth and J

W. Smyth and J. Moum, Length scales of turbulence in stably stratified mixing layers, Phys. Fluids12, 1327 (2000)

work page 2000
[15]

K. M. Smith, C. Caulfield, and J. Taylor, Turbulence in forced stratified shear flows, J. Fluid Mech.910, A42 (2021)

work page 2021
[16]

Arratia, C

C. Arratia, C. Caulfield, and J.-M. Chomaz, Transient perturbation growth in time-dependent mixing layers, J. Fluid Mech.717, 90 (2013)

work page 2013
[17]

Frame and A

P. Frame and A. Towne, Beyond optimal disturbances: a statistical framework for transient growth, J. Fluid Mech. 983, A2 (2024)

work page 2024
[18]

M. K. Verma,Physics of Buoyant Flows: From Instabil- ities to Turbulence(World Scientific)

work page
[19]

Mashayek and W

A. Mashayek and W. Peltier, The ‘zoo’ of secondary in- stabilities precursory to stratified shear flow transition. Part 1: Shear aligned convection, pairing, and braid in- stabilities, J. Fluid Mech.708, 5 (2012)

work page 2012
[20]

Mashayek and W

A. Mashayek and W. Peltier, The ‘zoo’of secondary in- stabilities precursory to stratified shear flow transition. Part 2: The influence of stratification, J. Fluid Mech. 708, 45 (2012)

work page 2012
[21]

Salehipour, W

H. Salehipour, W. Peltier, and A. Mashayek, Turbulent diapycnal mixing in stratified shear flows: The influence of Prandtl number on mixing efficiency and transition at high Reynolds number, J. Fluid Mech.773, 178 (2015)

work page 2015
[22]

Caulfield, Layering, instabilities and mixing in tur- bulent stratified flow, Annu

C. Caulfield, Layering, instabilities and mixing in tur- bulent stratified flow, Annu. Rev. Fluid Mech.53, 113 (2021)

work page 2021
[23]

S. A. Orszag, Accurate solution of the Orr-Sommerfeld stability equation, J. Fluid Mech.50, 689 (1971). End Matter Appendix A: Derivation of the Self-Similar Stability Equations.– a. Diffusive Non-Dimensionalisation.The typical non-dimensionalisation of shear flows, such as xdim =d 0 x, t dim = d0 U0 t(7) used to arrive at equations (1) and (2) from the main text

work page 1971
[24]

However, this approach is obstructive when con- sidering the evolution of an initial step function interface (d0 →0) as the associated Reynolds number vanishes

and where the subscript dim denotes dimensional vari- ables, relies on the initial finite shear layer half-depthd 0 [15, 21]. However, this approach is obstructive when con- sidering the evolution of an initial step function interface (d0 →0) as the associated Reynolds number vanishes. To resolve this singularity, we propose the alternative scaling xdim =...

work page

[1] [1]

Guided by the physical limits established in Fig

of kinetic energy – a nonlinear turbulent breakdown that falls outside the scope of linear theory [18–20]. Guided by the physical limits established in Fig. 2, we restrict our theoretical comparison to the linear regime t≥t del. We initialise both theoretical models using τ(k, t 0 =t del) from the DNS [Fig. 3(a), inset] and com- pute the subsequent spectr...

work page

[2] [2]

R. K. Newsom and R. M. Banta, Shear-flow instability in the stable nocturnal boundary layer as observed by Doppler lidar during CASES-99, J. Atmos. Sci.60, 16 (2003)

work page 2003

[3] [3]

S. A. Thorpe,An Introduction to Ocean Turbulence (Cambridge University Press)

work page

[4] [4]

Y. Liu, P. Tan, and L. Xu, Kelvin-Helmholtz instability in an ultrathin air film causes drop splashing on smooth surfaces, P. Natl. A. Sci.112, 3280 (2015)

work page 2015

[5] [5]

Boˇ cek, M

ˇZ. Boˇ cek, M. Petkovˇ sek, S. J. Clark, K. Fezzaa, and M. Dular, Kelvin-Helmholtz instability as one of the key features for fast and efficient emulsification by hydro- dynamic cavitation, Ultrason. Sonochem.108, 106970 (2024)

work page 2024

[6] [6]

J. W. Miles, On the stability of heterogeneous shear flows, J. Fluid Mech.10, 496 (1961)

work page 1961

[7] [7]

Howard, Note on a paper of John W

L. Howard, Note on a paper of John W. Miles, J. Fluid Mech.10, 509 (1961)

work page 1961

[8] [8]

Vieweg and C

P. Vieweg and C. Caulfield, Anisotropy of emergent large- scale dynamics in forced stratified shear flows, J. Fluid Mech.1029, A7 (2026)

work page 2026

[9] [9]

P. G. Drazin and W. H. Reid,Hydrodynamic stability (Cambridge University Press, 2004)

work page 2004

[10] [10]

Smyth, J

W. Smyth, J. Moum, and J. Nash, Narrowband oscilla- tions in the upper equatorial ocean. Part ii: Properties of shear instabilities, J. Phys. Oceanogr.41, 412 (2011)

work page 2011

[11] [11]

Pramanik and M

S. Pramanik and M. Mishra, Linear stability analysis of Korteweg stresses effect on miscible viscous fingering in 6 porous media, Phys. Fluids25(2013)

work page 2013

[12] [12]

Pramanik and M

S. Pramanik and M. Mishra, Nonlinear simulations of miscible viscous fingering with gradient stresses in porous media, Chem. Eng. Sci.122, 523 (2015)

work page 2015

[13] [13]

Fischer, S

P. Fischer, S. Kerkemeier, M. Min, Y.-H. Lan, M. Phillips, T. Rathnayake, E. Merzari, A. Tomboulides, A. Karakus, N. Chalmers, and T. Warburton, NekRS, a GPU-accelerated Spectral Element Navier–Stokes Solver, Parallel Comput.114, 102982 (2022)

work page 2022

[14] [14]

Smyth and J

W. Smyth and J. Moum, Length scales of turbulence in stably stratified mixing layers, Phys. Fluids12, 1327 (2000)

work page 2000

[15] [15]

K. M. Smith, C. Caulfield, and J. Taylor, Turbulence in forced stratified shear flows, J. Fluid Mech.910, A42 (2021)

work page 2021

[16] [16]

Arratia, C

C. Arratia, C. Caulfield, and J.-M. Chomaz, Transient perturbation growth in time-dependent mixing layers, J. Fluid Mech.717, 90 (2013)

work page 2013

[17] [17]

Frame and A

P. Frame and A. Towne, Beyond optimal disturbances: a statistical framework for transient growth, J. Fluid Mech. 983, A2 (2024)

work page 2024

[18] [18]

M. K. Verma,Physics of Buoyant Flows: From Instabil- ities to Turbulence(World Scientific)

work page

[19] [19]

Mashayek and W

A. Mashayek and W. Peltier, The ‘zoo’ of secondary in- stabilities precursory to stratified shear flow transition. Part 1: Shear aligned convection, pairing, and braid in- stabilities, J. Fluid Mech.708, 5 (2012)

work page 2012

[20] [20]

Mashayek and W

A. Mashayek and W. Peltier, The ‘zoo’of secondary in- stabilities precursory to stratified shear flow transition. Part 2: The influence of stratification, J. Fluid Mech. 708, 45 (2012)

work page 2012

[21] [21]

Salehipour, W

H. Salehipour, W. Peltier, and A. Mashayek, Turbulent diapycnal mixing in stratified shear flows: The influence of Prandtl number on mixing efficiency and transition at high Reynolds number, J. Fluid Mech.773, 178 (2015)

work page 2015

[22] [22]

Caulfield, Layering, instabilities and mixing in tur- bulent stratified flow, Annu

C. Caulfield, Layering, instabilities and mixing in tur- bulent stratified flow, Annu. Rev. Fluid Mech.53, 113 (2021)

work page 2021

[23] [23]

S. A. Orszag, Accurate solution of the Orr-Sommerfeld stability equation, J. Fluid Mech.50, 689 (1971). End Matter Appendix A: Derivation of the Self-Similar Stability Equations.– a. Diffusive Non-Dimensionalisation.The typical non-dimensionalisation of shear flows, such as xdim =d 0 x, t dim = d0 U0 t(7) used to arrive at equations (1) and (2) from the main text

work page 1971

[24] [24]

However, this approach is obstructive when con- sidering the evolution of an initial step function interface (d0 →0) as the associated Reynolds number vanishes

and where the subscript dim denotes dimensional vari- ables, relies on the initial finite shear layer half-depthd 0 [15, 21]. However, this approach is obstructive when con- sidering the evolution of an initial step function interface (d0 →0) as the associated Reynolds number vanishes. To resolve this singularity, we propose the alternative scaling xdim =...

work page