{"id":"9496458e-cc86-4c26-ad9e-a4268ec41dd4","arxiv_id":"2505.07433","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Highly resolved DNS indicates acoustic pulses precede and trigger baroclinic vorticity during onset of RTI and KHRTI, with new POD and enstrophy budget comparisons.","lead":"This paper uses massive direct numerical simulations to compare the onset of Rayleigh-Taylor and Kelvin-Helmholtz Rayleigh-Taylor instabilities, concluding that acoustic pressure pulses, not baroclinic vorticity, are the initial trigger in both cases. It also adds proper orthogonal decomposition and enstrophy budget analyses of the pressure field during onset.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Causal precedence of pressure over vorticity is inherited from the impulsive partition-removal initialization; without a smooth-start or convergence check the 'unequivocal' acoustic-trigger claim is not supported.","rationale":"The reader's weakest assumption identifies exactly the same initialization confound: impulsively removing the partition generates acoustic waves at t = 0, so the observed pressure-first, vorticity-after ordering may be predetermined by the setup rather than by the physics of RTI/KHRTI. I agree with that reading. The paper otherwise presents a careful and internally consistent DNS comparison, and the pressure and vorticity phenomenology is plausible; the issue is not internal inconsistency but external validity of the central causal claim. A single smooth-start test would settle whether the claim survives, and the absence of such a control is what makes the claim overreach. Because this can be tested directly and the paper's other contributions remain valuable, the conditional verdict remains appropriate; no change to the reader's verdict is needed.","tokens_in":19641,"tokens_out":4792,"duration_ms":49236,"concrete_test":"Run the RTI case again with the partition removed smoothly over a finite ramp duration, for example by linearly increasing the interface permeability from 0 to 1 over a time comparable to the acoustic crossing time of the domain (t_ramp ~ L/a), while keeping grid, time step, and all other parameters unchanged. If the O(10^-6) pressure pulses and their precedence over baroclinic vorticity are unchanged in amplitude and time lag, the acoustic-trigger claim is physical; if the signal amplitude scales with the ramp rate or the ordering disappears, the claimed trigger is an artifact of the impulsive initialization.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing concern is that the asserted temporal ordering—pressure first, baroclinic vorticity later—is built into the numerical protocol rather than demonstrated as physics. In Section 2, both the RTI and KHRTI cases are initialized by impulsively removing a partition at t = 0. For a compressible gas, this discontinuous release of the wall constraint is itself an acoustic source; it will generate pressure transients regardless of the instability mechanism. Section 3.1 then reports disturbance-pressure magnitudes of order 10^-6 and interprets the wave fronts in Figs. 2 and 3 as the trigger. Because the initial transient is the very signal being claimed as causal, the conclusion in Section 4 that 'acoustics trigger the instability, not the baroclinic vorticity' cannot be separated from the initialization artifact without a control. No convergence study, grid-refinement test, or alternative initial condition is presented, and the absence of explicit perturbations does not control for this, since the sharp density/temperature interface itself generates numerical pressure fluctuations at grid scale. The phrase 'demonstrated here unequivocally here for the first time' also conflicts with the paper's own citations [22,35], which already attribute the onset to acoustic excitation; this weakens the novelty of the causality claim. The concern is load-bearing because the entire central conclusion rests on this single temporal ordering.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents three-dimensional compressible Navier-Stokes direct numerical simulations of the onset of Rayleigh-Taylor instability (RTI, 4.19 billion grid points) and Kelvin-Helmholtz Rayleigh-Taylor instability (KHRTI, 480 million grid points), with setups patterned after experiments by Read and Akula et al. The authors analyze disturbance pressure fields, spectra, proper orthogonal decomposition, and a compressible enstrophy transport equation to compare the two instabilities. Their headline claim is that acoustic pulses trigger both instabilities, with baroclinic torque generation appearing as a downstream consequence, and that this is demonstrated 'unequivocally' for the first time. The paper also reports dominant viscous terms in the enstrophy budget and identifies regular and anomalous/shift POD modes.","tokens_in":19912,"tokens_out":6201,"duration_ms":57436,"significance":"If the causal claim were established, the paper would challenge the textbook attribution of RTI genesis to baroclinic torque and would provide a physically distinct picture of KHRTI onset. The computational effort is substantial, and the descriptive material—pressure-pulse morphology, spectral signatures, POD modes, and enstrophy budgets for two very large DNS runs—is potentially valuable as a reference database. The paper is explicit about its numerical methods and gives reproducible-looking parameter settings. However, the central claim is not supported by the evidence as presented: the simulations initialize the instability by impulsive partition removal, which is itself an acoustic source, and no convergence study, smooth-start control, or quantitative causality measure is provided. The finding is therefore currently a plausible interpretation of one numerical experiment per case rather than an unequivocal demonstration.","major_comments":[{"comment":"In both setups the instability is initiated by impulsive removal of the partition at t=0. For the compressible Navier-Stokes equations, this discontinuous removal is an acoustic source by construction, so the pressure pulses of order 10^-6 reported in Section 3.1 and their temporal precedence over baroclinic vorticity may reflect the initialization protocol rather than the physical onset mechanism of an unforced instability. The paper presents no smooth-start or alternative-initialization run, and no grid-refinement or time-step convergence study for the disturbance-pressure signal. Consequently, the sentence in Section 4 that acoustics trigger the instability 'has been demonstrated here unequivocally' overstates what the simulations can establish.","section":"Section 2 (initialization) and Section 3.1"},{"comment":"The causal ordering is inferred from visual inspection of contours at a small number of time instants: pressure pulses are visible in Figs. 2 and 3 at early times, and baroclinic vorticity features appear later in Figs. 11 and 12. No quantitative diagnostic connects the two: for example, a time-lagged correlation between p' and baroclinic torque, or a budget of the vorticity equation showing when T3 first exceeds a threshold, would be needed to demonstrate that the acoustic field causes the vorticity rather than accompanying it. Without such a measure, the claim 'acoustics trigger the instability' remains an interpretation of temporal ordering in a single run per case.","section":"Section 3.4, Figs. 11-12 and Section 4"},{"comment":"The paper's own literature review says that prior works [22] and [35] already attribute RTI and KHRTI genesis to acoustic excitation, yet Section 4 claims this is demonstrated 'for the first time.' This is internally inconsistent: either the prior works already made the causal claim, in which case the present novelty must be the specific comparative or quantitative evidence, or they did not, in which case the paper should specify what, exactly, is new. The wording should be corrected and the contribution re-scoped.","section":"Section 1 and Section 4"}],"minor_comments":[{"comment":"The word 'acosutic' appears in the discussion of pressure-pulse propagation and should be 'acoustic.'","section":"Section 3.1"},{"comment":"The terms 'unomodal' and 'nonomodal' appear to be typos for 'nonmodal'; please correct them throughout the spectrum discussion.","section":"Section 3.2"},{"comment":"The phrase 'demonstrated here unequivocally here for the first time' repeats 'here'; please revise.","section":"Section 4"},{"comment":"The scalar factors in Eq. (18), such as 2/rho^2 for T3 and 4/rho for T5, are not reflected in the term definitions given in the bullet list; please reconcile the notation so the terms match the equation exactly.","section":"Eq. (18) and bullet list"},{"comment":"The caption refers to 'line a', 'line b', and 'line c' while the legend uses solid, dash-dotted, and dashed lines; please align the labels.","section":"Figure 6"}],"recommendation":"major_revision","confidential_remarks":"The descriptive DNS content is likely sound and could find a useful home in this journal, but the strong causal claim is not supported by the current evidence. I would ask the authors to either provide a convergence study and at least one smooth-start or otherwise controlled initialization, or to substantially soften the causal claim and re-frame the paper as a comparative description of pressure and vorticity fields during onset. As written, the 'unequivocal' statement is the main risk to the manuscript's credibility."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline claim—that acoustic pulses trigger RTI and KHRTI before baroclinic vorticity—is not new to this paper: the authors have made it in [22] and [35], and they acknowledge as much in Section 3.1. What this paper adds is a side-by-side comparison of the onset stages of both instabilities using giant DNS, pressure-based POD, and a compressible enstrophy budget. Those analyses are new and, taken on their own, useful.\n\nThe DNS is genuinely large: 4.19 billion points for RTI, 480 million for KHRTI, with high-order dispersion-relation-preserving schemes and a non-overlapping parallel closure. The pressure fields at onset show a clean contrast—interface-normal pulses in RTI, radial upstream propagation in the shear plane for KHRTI—and the POD cleanly identifies regular pairs and an anomalous mode in each. The CETE budget, showing viscous terms dominant at onset and baroclinicity subdominant, is a useful benchmark, even if it depends on the bulk-viscosity regression model.\n\nNow the soft spots. The central causal ordering is not demonstrated by this data. Both setups remove the partition impulsively at t = 0, and that is an acoustic source by construction. So the observation that pressure pulses precede baroclinic vorticity is baked into the initialization; it would be surprising if it didn't happen. There is no smooth-start control, no convergence study, and the O(10^-6) pressure amplitudes are not checked against a grid-refinement set, so numerical artifacts cannot be ruled out. Section 4's 'unequivocally for the first time' also contradicts the paper's own citations and Section 3.1, which attribute the acoustic trigger to [22] and [35]. That internal contradiction isn't fatal, but it needs resolving.\n\nI would not reject the paper. The comparative POD and CETE analyses stand independently of the trigger claim, and the DNS database is valuable. A serious referee should ask for (1) a control run with a smooth partition release, (2) a grid-convergence check on the pressure signals, and (3) rewriting of the causal claim to separate what is new from what is inherited. With those, it could be a solid contribution. Without them, the 'unequivocal' claim should not appear in print.","headline":"Large careful DNS of RTI/KHRTI onset with useful new POD and enstrophy-budget comparisons, but the 'acoustic trigger' claim is inherited from the authors' prior work and is not made unequivocal by the impulsive initialization.","tokens_in":20417,"tokens_out":3332,"would_cite":true,"duration_ms":29368,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Acoustic pressure pulses, not baroclinic vorticity, initiate Rayleigh-Taylor and Kelvin-Helmholtz Rayleigh-Taylor instabilities.","keywords":["Rayleigh-Taylor instability","Kelvin-Helmholtz Rayleigh-Taylor instability","acoustic trigger","disturbance pressure field","compressible enstrophy transport equation","proper orthogonal decomposition","direct numerical simulation","baroclinic vorticity"],"falsifier":"Repeat the RTI or KHRTI simulation with a gradual, finite-time partition removal or an initially perturbed interface that produces no acoustic transient; if baroclinic vorticity appears at the interface before, or simultaneously with, any pressure pulse, the acoustic-trigger claim is wrong. A second check is a grid-convergence study: if the $10^{-6}$ pressure pulses shrink or change sign as the mesh is refined, they are numerical artifacts rather than physical signals.","tokens_in":19442,"feed_emoji":"🔊","tokens_out":11337,"duration_ms":90947,"temperature":0.7,"pith_summary":"This paper asks what actually starts the Rayleigh-Taylor instability (RTI) and the combined Kelvin-Helmholtz Rayleigh-Taylor instability (KHRTI) once a dense, cold layer of air is released above a light, warm layer. Using highly resolved compressible Navier-Stokes simulations—4.19 billion points for RTI and 480 million for KHRTI—the authors track disturbance pressure and vorticity from the first instants after the separating partition is removed. They conclude that the first signal is an acoustic pressure pulse, of order $10^{-6}$ relative to hydrostatic pressure, and that baroclinic vorticity, long treated as the trigger, is generated later as a consequence of that pulse. The distinction matters because it moves the causal origin of these instabilities from the vorticity production term to the compressible pressure field. The paper also documents how the pulse travels differently in the two cases and provides a benchmark dataset for onset-stage dynamics.","feed_headline":"Acoustics trigger Rayleigh-Taylor onset before vorticity forms","feed_subtitle":"Petascale simulations show pressure pulses arrive first, overturning the usual baroclinic-torque story.","key_machinery":"The object that carries the argument is the disturbance pressure field $p' = p - p_{\\mathrm{hydrostatic}}$, resolved to amplitudes of order $10^{-6}$ in two perpendicular planes, read together with the compressible enstrophy transport equation (CETE) budget. The machinery that makes the claim visible is the combination of a highly accurate dispersion-relation-preserving compact scheme, a non-overlapping parallel subdomain closure, and a non-zero bulk viscosity for air obtained by regression of measured acoustic attenuation data. The pressure field supplies the temporal and spatial ordering (pulse first, vorticity later), while the CETE budget shows that at onset the enstrophy growth is dominated by viscous terms—especially the bulk-viscosity term $T_4$—with the baroclinic term $T_3$ playing a secondary role.","core_discovery":"The paper's central claim is a causal ordering: in both RTI and KHRTI, the incipient mechanism of instability is acoustic pulse propagation, not baroclinic torque. After partition removal, compression and rarefaction fronts travel from the interface—normal to the interface in RTI, and radially upstream in the shear plane with interface-normal fronts in the perpendicular plane for KHRTI—and it is these travelling pressure gradients that misalign $\\nabla p$ and $\\nabla \\rho$, producing baroclinic vorticity at the side-wall/interface junctions. The authors state that this ordering has been demonstrated for the first time in their simulations. Supporting observations include: the first five POD modes of disturbance pressure capture 95.30% of the variance for RTI and 98.61% for KHRTI, each case showing two regular mode pairs and one anomalous shift mode; and the compressible enstrophy transport equation budget at onset is dominated by viscous terms, with the baroclinic term sub-dominant during the stage studied.","pith_inferences":["Inference: If the acoustic-trigger ordering holds, unforced experiments should show a measurable pressure transient arriving at the interface before any vorticity signature appears; placing fast pressure sensors near the interface in a laboratory RTI tank would test this directly.","Inference: The same causal ordering may extend to other impulsively started instabilities, such as Richtmyer-Meshkov flows, where a shock provides an explicit pressure pulse; in those settings the vortex-stretching and baroclinic terms in enstrophy budgets may likewise be consequences of the acoustic forcing rather than independent triggers.","Inference: Because the claim rests on the impulsive partition removal, a parameter study that gradually withdraws the partition over a finite time, or introduces the density mismatch without a pressure transient, would clarify how much of the acoustic trigger is initialization-dependent; this is a natural numerical extension of the present work.","Inference: The identification of an anomalous shift mode in the pressure POD suggests that reduced-order models of RTI/KHRTI onset could treat that mode as the acoustic transient degree of freedom and the regular pairs as the instability carriers, potentially allowing prediction of when the acoustic stage ends and vorticity-driven growth begins."],"forward_implications":["If the acoustic trigger is real, models that initialize RTI or KHRTI with prescribed vorticity perturbations omit the initiating mechanism and will mis-time the onset of mixing.","The different propagation geometry—interface-normal pulses for RTI, radial upstream pulses in the shear plane for KHRTI—implies that shear not only adds kinetic energy but redirects the acoustic wavefronts, which is why KHRTI spectra look chaotic over a wide wavenumber range.","The POD decomposition of disturbance pressure provides a low-dimensional description of onset: regular mode pairs track interfacial and propagating pressure signals, while the anomalous shift mode carries the transient adjustment; these are direct candidates for reduced-order models of instability onset.","The enstrophy budget result—viscous dominance with baroclinicity secondary at onset—recasts the usual emphasis on baroclinic torque for onset-stage enstrophy production, though the authors expect the baroclinic term to grow once coherent spikes, bubbles, and KH eddies form.","The pressure-probe time series show side-wall locations have stronger disturbance pressure than the center for both instabilities, identifying side-wall/interface junctions as the preferred sites where the acoustic pulse first converts into vorticity."],"supporting_citations":[{"why":"It supplies the earlier RTI acoustic-excitation result that this paper extends to KHRTI.","marker":"[22]"},{"why":"It provides the prior KHRTI simulation showing onset via pressure pulses and the upstream radial propagation direction used here.","marker":"[35]"},{"why":"It supplies the measured bulk-viscosity values for air whose regression makes the viscous term $T_4$ dominant in the enstrophy budget.","marker":"[39]"},{"why":"It supplies the non-overlapping parallel subdomain closure that allows the 4.19-billion-point RTI computation.","marker":"[57]"},{"why":"It defines the experimental RTI configuration whose geometry and parameters the RTI domain follows.","marker":"[2]"},{"why":"It defines the experimental KHRTI configuration that motivates the KHRTI domain size and velocity difference.","marker":"[32]"},{"why":"It derives the compressible enstrophy transport equation used for the budget analysis.","marker":"[54]"},{"why":"It provides the regular and anomalous shift-mode interpretation used to classify the POD modes.","marker":"[36]"}],"fun_headline_variants":["Acoustic pulses, not baroclinic torque, trigger Rayleigh-Taylor onset","Sound waves seed Rayleigh-Taylor and Kelvin-Helmholtz instabilities","First demonstration: acoustic trigger precedes vorticity in RTI and KHRTI","Petascale simulations show pressure pulses initiate fluid instabilities","Pressure fronts, not vortices, start Rayleigh-Taylor and Kelvin-Helmholtz onset"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole ordering depends on the assumption that instantly removing the partition—which necessarily fires acoustic waves into the fluid—is a faithful stand-in for how these instabilities begin in unforced experiments, so that the pressure-before-vorticity sequence is not an artifact of the starting procedure.","fun_headline_variants_meta":{"raw":{"variants":["Acoustic pulses, not baroclinic torque, trigger Rayleigh-Taylor onset","Sound waves seed Rayleigh-Taylor and Kelvin-Helmholtz instabilities","First demonstration: acoustic trigger precedes vorticity in RTI and KHRTI","Petascale simulations show pressure pulses initiate fluid instabilities","Pressure fronts, not vortices, start Rayleigh-Taylor and Kelvin-Helmholtz onset"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000333,"raw_usage":{"total_tokens":1869,"prompt_tokens":980,"completion_tokens":889,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":791}},"tokens_in":596,"tokens_out":889,"duration_ms":8187,"temperature":1.0,"reasoning_tokens":791,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:16:22.047495+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the RTI or KHRTI simulation with a gradual, finite-time partition removal or an initially perturbed interface that produces no acoustic transient; if baroclinic vorticity appears at the interface before, or simultaneously with, any pressure pulse, the acoustic-trigger claim is wrong. A second check is a grid-convergence study: if the $10^{-6}$ pressure pulses shrink or change sign as the mesh is refined, they are numerical artifacts rather than physical signals.","supporting_citations":[{"cited_title":"Sengupta, P","cited_arxiv_id":null,"evidence_quote":"It supplies the earlier RTI acoustic-excitation result that this paper extends to KHRTI."},{"cited_title":"Joshi, T","cited_arxiv_id":null,"evidence_quote":"It provides the prior KHRTI simulation showing onset via pressure pulses and the upstream radial propagation direction used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the measured bulk-viscosity values for air whose regression makes the viscous term $T_4$ dominant in the enstrophy budget."},{"cited_title":"Sundaram, A","cited_arxiv_id":null,"evidence_quote":"It supplies the non-overlapping parallel subdomain closure that allows the 4.19-billion-point RTI computation."},{"cited_title":"Read, Experimental investigation of turbulent mixing by Rayleigh- Taylor instability, Physica D Nonlinear Phenomena 12 (1-3) (1984) 45– 58","cited_arxiv_id":null,"evidence_quote":"It defines the experimental RTI configuration whose geometry and parameters the RTI domain follows."},{"cited_title":"Akula, P","cited_arxiv_id":null,"evidence_quote":"It defines the experimental KHRTI configuration that motivates the KHRTI domain size and velocity difference."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It derives the compressible enstrophy transport equation used for the budget analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the regular and anomalous shift-mode interpretation used to classify the POD modes."}],"review_version":1}