REVIEW 2 major objections 4 minor 43 references
Time-harmonic waves in Korteweg and nematic-Korteweg fluids
T0 review · 2 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The paper derives the Helmholtz–Korteweg equation for time-harmonic waves in Korteweg fluids, adds a nematic variant with orientational coupling, and predicts evanescent-wave penetration depth and orientation-dependent scattering.
desk verdict The paper's framework is sound but its headline penetration-depth prediction is algebraically wrong, so the advertised experimentally-verifiable claims do not hold as stated. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the Helmholtz–Korteweg operator $-\omega^2 - c_0^2\Delta + \rho_0^2 u_1\Delta^2$, a fourth-order linear operator whose nematic counterpart carries the extra term $\rho_0^2 u_2\nabla\cdot\nabla[\mathbf{n}\cdot H_S\mathbf{n}]$. It turns the time-harmonic ansatz into a quartic dispersion relation $-\omega^2 + c_0^2 k^2 + \rho_0^2(u_1+u_2\cos^2\xi)k^4=0$, whose roots include real wavenumbers for propagating waves and purely imaginary wavenumbers for evanescent waves. The characteristic time $\tau_1=4\rho_0\sqrt{u_1}/c_0^2$ sets the crossover: for $\omega\tau_1\ll1$ the equation behaves like the classical Helmholtz equation, while for $\omega\tau_1\gg1$ the fourth-order term dominates and yields the $\delta\propto\omega^{-1/2}$ penetration law. Boundary conditions enter as paired conditions on $S$ and $\Delta S$ with nematic corrections, and the circular obstacle is treated by a boundary-layer split into a Helmholtz solution plus an $O(\ell^2)$ correction.
What would settle it
Measure the penetration depth of an evanescent acoustic wave in an aligned nematic liquid crystal as a function of frequency and of the angle between the propagation direction and the director; the paper predicts $\delta\propto\omega^{-1/2}$ for $\omega\tau_1\gg1$ and a depth that is largest along the director. A depth that falls off exponentially with frequency, or an orientation dependence opposite to the prediction, would refute the central claim.
Extended reading notes
Core claim
The authors claim that a small time-harmonic condensation $S(\mathbf{x})e^{-i\omega t}$ in a Korteweg fluid obeys $-\omega^2 S - c_0^2\Delta S + \rho_0^2 u_1 \Delta^2 S = 0$, the Helmholtz–Korteweg equation (13), and that adding the nematic orientational stress changes the plane-wave dispersion relation to $-\omega^2 + c_0^2 k^2 + (\rho_0^2 u_1 + \rho_0^2 u_2 \cos^2\xi)k^4 = 0$ (51), where $\xi$ is the angle between the propagation direction and the nematic director. These dispersion relations match the earlier variational nematoacoustic result in the inviscid regime. The extension is to allow purely imaginary wavenumbers, giving an evanescent wave whose penetration depth is $\delta = c_0/(2\omega)\bigl[(-1+\sqrt{1+\tau_1^2\omega^2})/(\tau_1^2\omega^2)\bigr]^{-1/2}$, asymptotic to $\delta\approx c_0\sqrt{\tau_1\omega}/(2\omega)$ for $\omega\tau_1\gg1$. In the nematic case the effective time scale $\tau_2(\xi)$ makes both sound speed and penetration depth anisotropic, with maxima along the director. The authors further claim that a circular sound-soft obstacle scatters with reduced amplitude when the incident wave propagates orthogonally to the director, a boundary-layer effect captured as an $O(\ell^2)$ correction to classical Mie scattering.
Load-bearing premise
The load-bearing assumption is that the liquid-crystal alignment direction is uniform over the acoustic length scale, so all terms containing spatial gradients of that direction are dropped from the nematic wave equation; if the alignment bends or twists on that scale, the predicted dispersion and orientation-dependent penetration depths change.
Editorial extensions
If this is right
- At high frequencies ($\omega\tau_1\gg1$), evanescent sound in a Korteweg fluid penetrates a distance that shrinks as $\omega^{-1/2}$ rather than vanishing exponentially with frequency.
- In a nematic-Korteweg fluid the speed of sound is largest when propagation is parallel to the director, and evanescent penetration is deepest in that same direction, so rotating the director reorients both acoustic quantities.
- For a circular sound-soft obstacle, the scattered amplitude is predicted to be reduced when the incident wave is orthogonal to the director, an anisotropy that does not appear at a plane boundary for sound-soft conditions.
- The fourth-order boundary value problems split into paired conditions on $S$ and $\Delta S$, so standard Helmholtz solvers can be adapted to the Korteweg case with modest changes.
- As the Korteweg length $\ell$ tends to zero, the scattering solution is a singular perturbation of classical Helmholtz scattering, with an $O(\ell^2)$ boundary-layer correction, so ordinary acoustics predictions acquire small prescribed corrections.
Reading between the lines
- Extending the paper's results, a strongly anchored nematic cell, where the director bends near the walls, should show a spatially varying sound speed that the current $\nabla\mathbf{n}=0$ equation does not capture.
- The $\delta\propto\omega^{-1/2}$ scaling is distinctive enough that swept-frequency ultrasonic reflectometry on a single aligned cell could separate it from ordinary dissipative skin-depth behavior; the paper does not report such a measurement.
- Extending the circular-scattering result, an external field that rotates the director could act as an acoustic switch, changing the scattered amplitude by changing the effective wavenumber; the paper notes tunability but does not quantify switching dynamics.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives a fourth-order linear PDE, called the Helmholtz–Korteweg equation, for time-harmonic acoustic disturbances in Korteweg fluids, and a nematic variant that includes an orientational stress contribution. It then analyzes plane-wave dispersion, evanescent-wave penetration depth, reflection at plane interfaces, and scattering by a circular obstacle, claiming new experimentally testable predictions for nematic-Korteweg fluids. The dispersion relations (36) and (51) are correctly derived and match Virga's earlier results. The central new quantitative claim, however, is the penetration-depth formula (49) and its asymptotic form (50).
Significance. If the derivation were correct, the paper would offer a compact fourth-order framework for time-harmonic acoustics in Korteweg and nematic-Korteweg fluids, with concrete predictions for evanescent waves and scattering that could guide experiments. The paper correctly reproduces Virga's dispersion relation and provides a plausible boundary-layer treatment for the circular obstacle, and it uses numerical simulations to illustrate the orientation dependence of scattering. These strengths are, however, outweighed by a concrete algebraic error in the penetration-depth calculation, which invalidates the paper's principal advertised quantitative prediction. The additional assumption that the director field is undistorted on acoustic scales, although inherited from prior work, is not adequately justified for the new scattering geometries considered.
major comments (2)
- [Section III, Eqs. (47)–(50)] Equation (48) does not solve the quadratic (47). Substituting α = 4(c/c0)²[-1+√(1+τ1²ω²)]/(τ1²ω²) into (47) leaves a residual -1 + 4(r-1)(r-2)/(τ1²ω²) with r = √(1+τ1²ω²), which is not identically zero. The correct positive root is α = 2(c/c0)²[1+√(1+τ1²ω²)]/(τ1²ω²), leading to δ = c0 τ1/√[2(1+√(1+τ1²ω²))]. This corrected δ tends to the finite value c0τ1/2 as ωτ1→0 instead of diverging, and its high-frequency asymptote is δ ∼ c0√(τ1/(2ω)), which differs from the paper's (50) by a factor √2. Because the abstract and Section III present (49)–(50) as new, experimentally verifiable predictions, this error directly undermines the central quantitative claim of the paper.
- [Section II, before Eq. (25)] The assumption ∇n = 0 removes from the nematic Helmholtz–Korteweg equation (24) all director-gradient terms, including ∇n∇S and (∇S·n)(∇·n). This assumption is inherited from Virga's planar-wave analysis, but the present paper applies the resulting equation to new settings, notably the circular-obstacle scattering problem in Section V and the director-discontinuity reflection in Figure 2, where director gradients are not obviously negligible. The authors do not justify that the director is effectively constant on acoustic wavelengths in these geometries, so the nematic dispersion relation (51) and the orientation-dependent scattering predictions are not fully established.
minor comments (4)
- [Section I, Eq. (10)] Equation (10) contains a spurious minus sign on the right-hand side: from the general wave equation (9), substituting the time-harmonic ansatz gives -ρ0ω²S = ∇·(∇·σ), so the second minus sign in (10) is an error. The final Helmholtz–Korteweg equation (13) is nevertheless correct, but the intermediate formula is misleading.
- [Section I, after Eq. (12)] The text says 'divide by ρ0 e^{-ωt}', but the exponential in the ansatz (7) is e^{-iωt}; this should be corrected.
- [Throughout] There are several typographical and formatting issues in the equations, such as the spacing in τ1²ω² in (38) and the missing ℜ operator in (12). These should be cleaned up.
- [Section V, Eq. (64)] The condition γ ≫ 1 ≫ ℓ mixes a dimensionless parameter γ with a parameter ℓ that presumably carries units of length. The statement should be made dimensionless or clarified.
Circularity Check
No significant circularity; derivation is self-contained and Virga's dispersion relation is used only as a consistency check.
full rationale
The derivation is self-contained. The Helmholtz–Korteweg equation (13) is obtained by linearizing the Euler–Korteweg system, substituting the time-harmonic ansatz, and dropping O(epsilon^2) terms; no target result is used as an input. The nematic equation (26) is obtained the same way from the constitutive relation (3), with the gradient-free director assumption stated explicitly rather than smuggled in. The plane-wave dispersion relations (36), (38), and (51) are consequences of substituting the plane-wave ansatz into the derived equations, and Virga's result is used only as an a posteriori consistency check. The evanescent penetration depth (49) is computed from the dispersion relation, not fitted to data; parameters such as u1, u2, c0, and rho0 are material constants from the constitutive theory. The only self-citations ([26], and the Firedrake manual [33] in a footnote) are peripheral: [26] is a remark about compressible Leslie–Ericksen equations and [33] is a software citation; neither supplies a load-bearing premise. The algebraic objection to eq. (48) raised by the reviewer concerns correctness of the formula, not circularity, and does not change this assessment.
Assumptions & free parameters
assumptions (7)
- domain assumption Small perturbation linearization: the condensation s and velocity v are O(ε) with ε<<1, and terms O(ε²) are neglected.
- domain assumption Time-harmonic ansatz s = Re(S(x)e^{-iωt}) and analytic continuation of the resulting equation to the complex plane.
- domain assumption Undistorted director field at the acoustic scale: ∇n = 0.
- ad hoc to paper Boundary conditions for the fourth-order equation are obtained by postulating that the classical Helmholtz sound-soft/sound-hard conditions also apply, giving homogeneous Dirichlet/Neumann conditions on both S and ΔS (or their normal derivatives).
- domain assumption The wave number decomposition k = k(R)+ik(I) with k(R)=ω/c and k(I) = -ω/c√α, α≥0, inspired by total internal reflection.
- domain assumption Physical regime for scattering: ρ0²u1 ≈ ℓ², ρ0²u2 ≈ γ^{-1}ℓ², with γ≫1≫ℓ.
- domain assumption The identity ∇·((n⊗n)∇S) = n·H_S n, which uses ∇·n = 0.
Cite this review
Pith. "Pith review of Time-harmonic waves in Korteweg and nematic-Korteweg fluids." pith.science (2026). https://pith.science/paper/NJDBNKWG
@misc{pith2026241113354,
author = {Pith},
title = {Pith review of: Time-harmonic waves in Korteweg and nematic-Korteweg fluids},
year = {2026},
howpublished = {\url{https://pith.science/paper/NJDBNKWG}},
note = {Machine review of arXiv:2411.13354}
}
read the original abstract
We derive the Helmholtz--Korteweg equation, which models acoustic waves in Korteweg fluids. We further derive a nematic variant of the Helmholtz-Korteweg equation, which incorporates an additional orientational term in the stress tensor. Its dispersion relation coincides with that arising in Virga's analysis of the Euler-Korteweg equations, which we extend to consider imaginary wave numbers and the effect of boundary conditions. In particular, our extensions allow us to analyze the effect of nematic orientation on the penetration depth of evanescent plane waves, and on the scattering of sound waves by obstacles. Furthermore, we make new, experimentally-verifiable predictions for the effect of boundary conditions for a modification of the Mullen-L\"uthi-Stephen experiment, and for the scattering of acoustic waves in nematic-Korteweg fluids by a circular obstacle.
Figures
Reference graph
Works this paper leans on
-
[26]
S. Benzoni-Gavage, R. Danchin, and S. Descombes, On the well-posedness for the Euler-Korteweg model in sev- eral space dimensions, Indiana University Mathematics Journal 56, 1499 (2007)
work page 2007
-
[1]
As previously discussed, (1) can be de- rived within the hyperelasticity framework, in particular, the elastic energy density associated with (1) is given by W (ρ, ∇ρ) = 1 2 c2 0ρ2 + 1 2 u1∥∇ρ∥2. (14) This energy shows that as u1 increases, a larger penalty is imposed on the gradient of the density field, and in the limit this imposes a constraint of spat...
-
[2]
From (2) we know that the acoustic pres- sure is given by c2 0ρ0S(⃗ x) − ρ3 0u1∆S(⃗ x) = 0
Sound-soft boundary conditions, which correspond to imposing that the acoustic pressure vanishes along Γ. From (2) we know that the acoustic pres- sure is given by c2 0ρ0S(⃗ x) − ρ3 0u1∆S(⃗ x) = 0. (17) Thus, assuming the sound-soft boundary condi- tions for the classical Helmholtz equations also ap- ply, we obtain that the sound-soft boundary condi- tion...
-
[3]
Sound-hard boundary conditions, which corre- spond to imposing that the normal component of the fluid velocity ⃗ ν· ⃗ vvanishes along Γ. Using (8) and assuming the fluid velocity is also time- harmonic, we get −iωρ0(⃗V · ⃗ ν) = c2 0∂⃗ νS(⃗ x) − ρ2 0u1∂⃗ ν∆S(⃗ x). (18) Thus, assuming that sound-hard boundary condi- tions for the classical Helmholtz equatio...
-
[4]
Impedance boundary conditions, which correspond to imposing that the normal component of the fluid velocity ⃗ ν· ⃗ vis proportional to the excess pressure along Γ. Imposing the acoustic pressure (2) to be proportional to ⃗ ν· ⃗ v, computed from (8) using the assumption that the fluid velocity is also time- harmonic, we obtain c2 0ρ0S(⃗ x) − ρ3 0u1∆S(⃗ x) ...
-
[5]
Sound-soft boundary conditions will change since excess pressure is now defined, from (4), as c2 0ρ0S(⃗ x) − ρ3 0u1∆S(⃗ x) − u2ρ3 0 ⃗ n· HS⃗ n = 0. (29) Sound-soft boundary conditions thus correspond to imposing homogeneous Dirichlet boundary condi- tions on S(⃗ x) and ∆S(⃗ x) = − u2 u1 ⃗ n· HS⃗ n . (30)
-
[6]
Sound-hard boundary conditions also change since the normal component of the fluid velocity⃗ ν· ⃗ vnow satisfies the equation −iωρ0(⃗ v· ⃗ ν) = c2 0∂⃗ νS(⃗ x) − ρ2 0u1∂⃗ ν∆S(⃗ x) (31) − ρ2 0u2∂⃗ ν ⃗ n· HS⃗ n . (32) Sound-hard boundary conditions thus correspond to imposing homogeneous Neumann boundary con- ditions on S(⃗ x) and ∂⃗ ν∆S(⃗ x) = − u2 u1 ∂⃗ ν ...
-
[7]
Some computation shows that the impedance boundary conditions for the nematic Helmholtz– Korteweg equation are equivalent to imposing Robin boundary conditions on S(⃗ x) and ∂⃗ ν∆S(⃗ x) = iζ∆S(⃗ x) + iζ u2 u1 ⃗ n· HS⃗ n − u2 u1 ∂⃗ ν ⃗ n· HS⃗ n . (34) III. PLANE W A VES We wish to build intuition about the Helmholtz– Korteweg and nematic Helmholtz–Korteweg...
Show all 43 references
-
[8]
For sound-soft boundary conditions along the real axis we have s0eikd1x1 + s0Aeikd1x1 = 0, thus we need to impose A = −1
-
[9]
For sound-hard boundary conditions along the real axis we have s0ikd2eikd1x1 − s0ikd2Aeikd1x1 = 0, thus we need to impose A = 1
-
[10]
Mathematical modelling of multicomponent systems
For impedance boundary conditions along the real axis we have iA(d2k − ζ)eikd1x1 = −i(d2k + ζ)eikd1x1 , (59) thus we need to impose A = k d2+ζk −1 d2−ζk −1 . First we focus our attention on the case where k ∈ R. Using our previous results, we know that the wave- number k of th...
-
[11]
D. J. Korteweg, Sur la forme que prennent les ´ equations du mouvements des fluides si l’on tient compte des forces capillaires caus´ ees par des variations de densit´ e con- sid´ erables mais connues et sur la th´ eorie de la capillarit´ e dans l’hypoth` ese d’une variation c...
1901
-
[12]
Truesdell, W
C. Truesdell, W. Noll, and S. S. Antman, The non-linear field theories of mechanics (Springer, Berlin, 2004)
2004
-
[13]
Lauro, A note on a Korteweg fluid and the hydro- dynamic form of the logarithmic Schr¨ odinger equation, Geophysical and Astrophysical Fluid Dynamics 102, 373 (2008)
G. Lauro, A note on a Korteweg fluid and the hydro- dynamic form of the logarithmic Schr¨ odinger equation, Geophysical and Astrophysical Fluid Dynamics 102, 373 (2008)
2008
-
[14]
Benzoni-Gavage, Planar traveling waves in capil- lary fluids, Differential and Integral Equations 26, 439 (2013)
S. Benzoni-Gavage, Planar traveling waves in capil- lary fluids, Differential and Integral Equations 26, 439 (2013)
2013
-
[15]
E. G. Virga, Variational theory for nematoacoustics, Physical Review E 80, 031705 (2009)
2009
-
[16]
G. G. Capriz, Continua with microstructure , Springer Tracts in Natural Philosophy, Vol. 35 (Springer-Verlag, New York, 1989)
1989
-
[17]
R. A. Toupin, Elastic materials with couple-stresses, Archive for Rational Mechanics and Analysis 11, 385 (1962)
1962
-
[18]
R. A. Toupin, Theories of elasticity with couple-stress, Archive for Rational Mechanics and Analysis 17, 85 (1964)
1964
-
[19]
Fried and M
E. Fried and M. E. Gurtin, Tractions, balances, and boundary conditions for nonsimple materials with appli- cation to liquid flow at small-length scales, Archive for Rational Mechanics and Analysis 182, 513 (2006)
2006
-
[20]
Giovangigli, Kinetic derivation of diffuse-interface fluid models, Physical Review E 102 (2020)
V. Giovangigli, Kinetic derivation of diffuse-interface fluid models, Physical Review E 102 (2020)
2020
-
[21]
J. E. Dunn and J. Serrin, On the thermomechanics of 10 interstitial working, Archive for Rational Mechanics and Analysis 88, 95 (1985)
1985
-
[22]
D. M. Anderson, G. B. McFadden, and A. A. Wheeler, Diffuse-interface methods in fluid mechanics, Annual Re- view of Fluid Mechanics 30, 139 (1998)
1998
-
[23]
M. M. Mehrabadi, S. C. Cowin, and M. Massoudi, Con- servation laws and constitutive relations for density- gradient-dependent viscous fluids, Continuum Mechanics and Thermodynamics 17, 183 (2005)
2005
-
[24]
Benzoni-Gavage, S
S. Benzoni-Gavage, S. Descombes, D. Jamet, and L. Mazet, Structure of Korteweg models and stability of diffuse interfaces, Interfaces and Free Boundaries 7, 371 (2005)
2005
-
[25]
Mathematical Fluid Dynamics
S. Benzoni-Gavage, Propagating phase boundaries and capillary fluids (2010), notes from the international sum- mer school on “Mathematical Fluid Dynamics”, held at Levico Terme (Trento)
2010
-
[27]
Murata and Y
M. Murata and Y. Shibata, The global well-posedness for the compressible fluid model of Korteweg type, SIAM Journal on Mathematical Analysis 52, 6313 (2020)
2020
-
[28]
K. Tsuda, On the existence and stability of time peri- odic solution to the compressible Navier–Stokes equation on the whole space, Archive for Rational Mechanics and Analysis 219, 637 (2016)
2016
-
[29]
Giesselmann, C
J. Giesselmann, C. Lattanzio, and A. E. Tzavaras, Rela- tive energy for the Korteweg theory and related Hamil- tonian flows in gas dynamics, Archive for Rational Me- chanics and Analysis 223, 1427 (2017)
2017
-
[30]
A. M. Sonnet and E. G. Virga, Dissipative ordered fluids : theories for liquid crystals (Springer, New York, 2012)
2012
-
[31]
Benzoni-Gavage and D
S. Benzoni-Gavage and D. Chiron, Long wave asymptotics for the Euler–Korteweg system, Revista Matem´ atica Iberoamericana34, 245 (2018)
2018
-
[32]
J. L. Dion and A. D. Jacob, A new hypothesis on ul- trasonic interaction with nematic liquid crystal, Applied Physics Letters 31, 490 (1977)
1977
-
[33]
J. L. Dion, The orienting action of ultrasound on liquid crystals related to the theorem of minimum entropy pro- duction, Journal of Applied Physics 50, 2965 (1979)
1979
-
[34]
M. E. Mullen, B. L¨ uthi, and M. J. Stephen, Sound veloc- ity in a nematic liquid crystal, Physical Review Letters 28, 799 (1972)
1972
-
[35]
J. V. Selinger, M. S. Spector, V. A. Greanya, B. T. Wes- lowski, D. K. Shenoy, and R. Shashidhar, Acoustic re- alignment of nematic liquid crystals, Physical Review E 66, 051708 (2002)
2002
-
[36]
P. E. Farrell, G. Russo, and U. Zerbinati, Kinetic deriva- tion of a compressible Leslie–Ericksen equation for rari- fied calamitic gases (2023), arXiv:2312.15210 [math-ph]
2023 arXiv
-
[37]
The PDE analysis for these equations, and a numerical analysis of their discretisations, will be presented elsewhere
Since the Helmholtz–Korteweg equation is fourth order, we employ C 1-conforming Argyris finite elements [32] in Firedrake [33]. The PDE analysis for these equations, and a numerical analysis of their discretisations, will be presented elsewhere
-
[38]
E. G. Virga, Variational theories for liquid crystals , Applied Mathematics and Mathematical Computation, Vol. 8 (CRC Press, Taylor & Francis Group, 1994)
1994
-
[39]
M. I. Vishik and L. A. Lyusternik, The solution of some perturbation problems for matrices and selfadjoint or non-selfadjoint differential equations, Russian Mathe- matical Surveys 15, 1 (1960)
1960
-
[40]
Advanced numerical meth- ods for PDEs
A. Moiola, Scattering of time-harmonic acoustic waves: Helmholtz equation, boundary integral equations and bem, Lecture notes for the “Advanced numerical meth- ods for PDEs” class, University of Pavia, Department of Mathematics (2024)
2024
-
[41]
A. F. Oskooi, L. Zhang, Y. Avniel, and S. G. Johnson, The failure of perfectly matched layers, and towards their redemption by adiabatic absorbers, Optics Express 16, 11376 (2008)
2008
-
[42]
J. H. Argyris, I. Fried, and D. W. Scharpf, The TUBA family of plate elements for the matrix displacement method, The Aeronautical Journal 72, 701 (1968)
1968
-
[43]
D. A. Ham, P. H. J. Kelly, L. Mitchell, C. J. Cotter, R. C. Kirby, K. Sagiyama, N. Bouziani, S. Vorderwuel- becke, T. J. Gregory, J. Betteridge, D. R. Shapero, R. W. Nixon-Hill, C. J. Ward, P. E. Farrell, P. D. Brubeck, I. Marsden, T. H. Gibson, M. Homolya, T. Sun, A. T. T. Mc...
2023
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.