{"id":"5166a7a1-786c-4b31-8bec-b471cf6ff805","arxiv_id":"2411.09021","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"This paper derives and validates a Keller-Miksis based model with time-delayed acoustic interactions between bubbles, showing finite-speed effects matter in dense clusters.","lead":"A new low-dimensional model simulates cavitation bubbles interacting through pressure waves that travel at finite speed, instead of the usual instant-action approximation. It predicts measurably different behavior in dense bubble clusters, where delayed acoustic feedback changes how bubbles expand and collapse.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reported QA-vs-IC differences may be controlled by the undocumented, case-dependent pruning coefficient C of Eq. (A2); no sensitivity analysis is given.","rationale":"The reader's weakest_assumption targets the linear superposition of spherical potentials and the Appendix A averaging of retarded invariants over emission nodes inside the receiving bubble. For the tested configurations, however, those effects are secondary: R_i is much smaller than the inter-bubble spacing, so the averaging window is a small fraction of the acoustic period, and the clusters are dilute enough that multiple scattering is unlikely to dominate. The more concrete and potentially fatal gap is the pruning coefficient C in Eq. (A2). The paper itself admits C is case-dependent and chosen by trial and error, yet the reported QA-vs-IC differences are the main evidence for the central claim. If C controls the effective interaction radius, the validation does not establish the claimed finite-speed effects. Because the code is openly available, a sensitivity study can settle this directly. The reader already issued CONDITIONAL, and this concern supports that verdict without changing it, so VERDICT remains UNCHANGED.","tokens_in":19170,"tokens_out":18338,"duration_ms":215651,"concrete_test":"Rerun the 51x51 bubble screen (Sec. VI B) and the spherical cluster cases (Secs. VI C-E) with pruning disabled (C = 0) and with C reduced by factors of 10 and 100 from the values used in the paper, reporting the actual C values. If the radial-amplitude patterns and ambient-pressure traces change by amounts comparable to the reported quasi-acoustic vs incompressible differences, the headline result is controlled by the pruning coefficient; if the results converge over this range, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central validations (bubble screen in Sec. VI B, cavitation onset in Secs. VI D and VI E) are produced by the numerical implementation of Eq. (27), whose interaction sums depend on the node-pruning criterion in Appendix A, Eq. (A2): an emission node is deleted when its individual pressure amplitude falls below C p_inf(t)/N. The paper states that C is case-dependent and 'determined mainly by trial and error', but it reports no C values and no sensitivity study. Because the criterion is applied node-by-node before summation, it effectively sets an interaction range: in the 51x51 bubble screen with D = 400 R0, an order-of-magnitude estimate at C = 1 gives a pruning radius comparable to the nearest-neighbor distance, so the nearest-neighbor interaction itself is borderline discarded. The dramatic quasi-acoustic vs incompressible differences in Fig. 6 could then reflect the retained interaction radius set by C rather than finite propagation speed. The same concern applies to the cluster and onset-of-cavitation cases. The underlying equations may be sound, but the reported numerical demonstrations are not yet shown to be independent of this free parameter.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a low-dimensional model for the radial dynamics, acoustic emissions, and acoustic interactions of cavitation bubbles and bubble clusters, based on a quasi-acoustic assumption that retains first-order Mach number corrections and finite propagation speed. The model couples a Keller-Miksis equation for each bubble with a Lagrangian wave-tracking scheme for the emitted acoustic field and computes inter-bubble interaction pressures from superposed retarded potentials. The claims are tested against literature cases: frequency response of small polydisperse clusters, resonance patterns in a bubble screen, asymmetric collapse of a spherical cluster, and pressure-drop-induced cavitation onset in two-bubble and cluster configurations. The central claim is that finite-speed (quasi-acoustic) interactions materially change predicted dynamics in dense systems compared with incompressible-interaction models.","tokens_in":19373,"tokens_out":4073,"duration_ms":40594,"significance":"If the central claim holds, the paper offers an efficient computational tool that is more faithful than instantaneous-interaction RP-type models for dense bubble clusters, and it is accompanied by a public implementation (APECSS v1.7) and data repository. The derivation of the Keller-Miksis equation from the quasi-acoustic assumption is internally consistent, and the emission-node expressions reduce correctly to the bubble-wall boundary conditions. The independent literature comparisons (Haghi et al., Fan et al., Ida, Maeda and Colonius) are appropriate and no target constants are fitted to reproduce those results. However, the numerical demonstrations that underpin the headline results depend on an undocumented, case-dependent pruning coefficient C in Eq. (A2), and the validation against Fan et al. is only qualitative; these issues must be resolved before the significance of the claimed time-delay effects can be assessed.","major_comments":[{"comment":"The pruning criterion for emission nodes contains the free coefficient C, which the authors state is case-dependent and 'determined mainly by trial and error'. No C values are reported for any of the figures, and no sensitivity study is given. Because the central QA-vs-IC differences in Fig. 6 and the cavitation-onset results in Figs. 10-17 are produced by the numerical implementation of Eq. (27) using this pruning, the observed differences could in part be controlled by C rather than by finite propagation speed. This is a load-bearing issue. The authors should report the actual C values used for each test case and provide a sensitivity analysis over a range of C, demonstrating that the qualitative and quantitative conclusions are unchanged.","section":"Appendix A, Eq. (A2); Sections VI B, VI D, VI E"},{"comment":"The interaction model represents the retarded invariants phi_j and g_j of a neighbor by the mean of all emission nodes located inside the receiving bubble. This is an ad-hoc approximation: it averages over a volume that may contain only a few nodes in sparse cases or many nodes with high gradients in dense cases, and it neglects the finite size and curvature of the receiving bubble. The paper provides no validation of this approximation. Since the authors themselves note that other weighting schemes are possible, the sensitivity of the reported cluster dynamics to this averaging choice should be assessed, for instance by comparing with a higher-order interpolation or with direct evaluation of the retarded potentials at the bubble center.","section":"Appendix A, interaction-averaging approximation"},{"comment":"The comparison with Fan, Li, and Fuster for the bubble screen is only qualitative. The text states 'good agreement' but also acknowledges 'qualitative differences' in amplitude patterns, and the color scales in Fig. 6 indicate only order-of-magnitude agreement. Since this is the principal demonstration that finite-speed interactions create spatial patterns, a quantitative metric (e.g., relative amplitude error, pattern correlation, or a profile comparison along a centerline) is needed to support the claim of agreement and to allow the reader to judge the significance of the QA-IC differences.","section":"Section VI B, Fig. 6"}],"minor_comments":[{"comment":"The name 'Minneart' is a misspelling; the correct form is 'Minnaert' (as in the cited reference).","section":"Section I and Section VI B"},{"comment":"The axis labels appear as 'R [ m]' and 't [ s]'; the radius axis should be labeled with the proper unit (likely micrometers, 'R [µm]') to avoid ambiguity.","section":"Figures 10 and 11"},{"comment":"The notation N_i for the number of neighbor bubbles is introduced after N has been used for the total number of bubbles; this distinction should be clarified in the text to prevent confusion.","section":"Eq. (26)"},{"comment":"The phrase 'progressive sinusoidal wave' should probably read 'propagating sinusoidal wave'.","section":"Section VI A"},{"comment":"The pruning condition is said to be 'based on Eq. (20)', but Eq. (20) uses the liquid density rho, whereas Eq. (A2) writes rho_0; the notation should be made consistent.","section":"Appendix A, Eq. (A2)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of Physics of Fluids and the underlying formulation is plausible, but the load-bearing numerical results are not yet reproducible because the pruning coefficient C in Eq. (A2) is unreported and case-dependent. I would like to see a sensitivity analysis and reported C values, together with a quantitative comparison for the bubble-screen test case, before recommending acceptance. The self-citation to the authors' prior APECSS framework is not excessive, and the open-source code is a positive feature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real contribution here is the combination: the authors take their existing Lagrangian wave-tracking framework for single bubbles and use it to build a low-dimensional, time-domain model of acoustic interactions in bubble clusters. That is a genuinely useful step. The derivation is internally consistent, the quasi-acoustic reduction to Keller-Miksis is standard but clean, and the model reduces correctly to the bubble-wall boundary conditions. I also give them credit for comparing against several independent results (Haghi, Fan, Ida, Maeda and Colonius) rather than only their own cases, and for shipping the APECSS code and data. No constants are fitted to the target results, which matters.\n\nThe soft spot is real and it is the one the stress test flags. Appendix A introduces a pruning coefficient C, described as case-dependent and \"determined mainly by trial and error,\" with no reported values and no sensitivity analysis. Because the criterion is applied node-by-node before summation, it effectively sets an interaction radius. For the 51x51 bubble screen with D = 400 R0, the order-of-magnitude estimate in the stress test is enough to worry me: at C = 1 the pruning radius is comparable to the nearest-neighbor distance, so the retained interaction range, not the finite propagation speed, could be driving a substantial part of the dramatic quasi-acoustic vs incompressible differences in Fig. 6. That concern applies to the cluster and cavitation-onset cases too. The stress test may be wrong, but the paper as written does not provide the evidence to know. A few lines reporting C values and a sensitivity sweep would do a lot.\n\nA secondary issue: the interaction-averaging in Appendix A, replacing retarded invariants of a neighbor by mean values over emission nodes currently inside the receiving bubble, is not analyzed for error. For dense clusters this approximation might be comparable in size to the finite-speed corrections. I would not call it fatal, but it deserves a paragraph.\n\nThe validation is otherwise adequate in spirit. Some comparisons are only qualitative (Fig. 6 shows order-of-magnitude agreement with Fan et al.), and the authors explain the pattern differences plausibly by time-domain vs frequency-domain solutions, but the explanation is not quantified.\n\nThe core model is believable and the central idea is worth taking seriously. The paper deserves a serious referee, not a desk rejection. My recommendation would be a major revision in which the authors report C values, run sensitivity tests, and ideally demonstrate convergence of the interaction-averaging approximation. This is the kind of paper the cavitation and bubbly-flow community will want to use, and the free parameter needs to be under control before that happens.","headline":"A useful extension of wave tracking to time-delayed bubble interactions, but the free pruning coefficient C needs a sensitivity study before the headline claims are fully credible.","tokens_in":19871,"tokens_out":1268,"would_cite":true,"duration_ms":15570,"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":"Accounting for the finite speed of sound materially changes predicted dynamics of dense cavitation bubble clusters, including the onset of cavitation.","keywords":["cavitation","bubble clusters","acoustic emissions","time-delayed interactions","quasi-acoustic assumption","Keller-Miksis equation","Lagrangian wave tracking","cavitation inception"],"falsifier":"In a dense monodisperse spherical cluster (e.g., 250 bubbles of 2 µm radius in a 232 µm cluster), compare the predicted mean radius evolution and pressure at the cluster center under a 1.75 µs tension pulse against a fully resolved compressible two-phase simulation. If the resolved simulation shows no difference from an incompressible-interaction model in the onset time or peak pressure, the paper's central claim of appreciable time-delay effects in dense clusters is not supported.","tokens_in":18941,"feed_emoji":"🫧","tokens_out":7243,"duration_ms":58697,"temperature":0.7,"pith_summary":"This paper tries to establish that the finite propagation speed of acoustic waves, not just liquid compressibility, must be included to predict how cavitation bubbles interact in clusters. The authors derive a consistent set of quasi-acoustic equations, accurate to first order in the Mach number, for the radial motion of each bubble and for the pressure and velocity fields that the bubbles radiate. Coupling bubbles through retarded potentials carried along outgoing sound characteristics, they show that incompressible-interaction models are unreliable for dense and large bubble systems, overestimating resonance amplitudes in a bubble screen and mispredicting the onset of cavitation near large neighbors. If correct, this provides a low-cost path to more faithful simulations of bubble clusters in medical and engineering applications.","feed_headline":"Finite sound speed alters cavitation in bubble clusters","feed_subtitle":"Time-delayed acoustic interactions shift resonance amplitudes and the onset of cavitation in dense clusters.","key_machinery":"The quasi-acoustic assumption: the Mach number is small, density and sound speed are nearly constant, so the velocity potential satisfies the linear wave equation and both the potential and its time derivative propagate unchanged along outgoing characteristics at the liquid sound speed. These invariants are stored on Lagrangian emission nodes emitted at each bubble wall; superposing the retarded potentials of all neighbor bubbles gives the local velocity and pressure fields, and the resulting driving pressure for each bubble. This yields a Keller-Miksis-style radial equation per bubble that is accurate to first order in the Mach number, with the interaction terms entering through the driving pressure and its time derivative.","core_discovery":"The central claim is that time-delayed, first-order-compressible bubble-bubble interactions materially change the predicted dynamics of dense mono- and polydisperse bubble clusters compared with models that treat interactions as instantaneous. The authors' quasi-acoustic model, built from the Keller-Miksis radial equation and a Lagrangian wave-tracking scheme, reproduces established collective phenomena while revealing strong quantitative differences: at resonance, the incompressible model predicts radial oscillation amplitudes that are only about 2% of the quasi-acoustic value for most bubbles in a bubble screen; in the two-bubble onset-of-cavitation case, the smaller bubble expands more under quasi-acoustic interactions because the pressure rise from the larger neighbor arrives late; and in tension-pulsed spherical clusters, incompressible interactions produce larger, more persistent pressure peaks than the time-delayed model. The paper therefore argues that finite propagation speed is a first-order effect for dense clusters, not a small correction.","pith_inferences":["The Lagrangian emission-node framework already outputs the local velocity and pressure fields, so extending it to track bubble translation under Bjerknes forces should be straightforward and would test the spherical-stationary assumption.","The prediction that time delays lower the pressure felt by a small bubble during expansion suggests that experimentally measured cavitation-inception thresholds near large bubbles could differ from incompressible-model predictions; a calibrated two-bubble experiment measuring the critical negative pressure would be a direct test.","The pruning criterion for obsolete emission nodes relies on a case-dependent coefficient C; a systematic error analysis of this pruning would determine how much it affects the accuracy of the predicted delay effects in polydisperse clusters.","For stronger collapses the constant-density, constant-sound-speed assumption breaks down (errors below 6% up to 50 MPa for water), so coupling this delay framework to a real-fluid equation of state might extend the model to shock-producing collapses."],"forward_implications":["For dense bubble screens near resonance, neglecting the propagation delay can underestimate radial oscillation amplitudes by roughly two orders of magnitude.","The onset of cavitation of a small bubble near a larger one is predicted to occur at less negative applied pressures (or with larger expansion) when time delays are included.","In polydisperse clusters excited by a tension pulse, quasi-acoustic interactions make the system more damped, with faster decay of oscillations and shorter-lived excited states.","Because the model retains first-order compressibility at a computational cost typical of ordinary differential equation systems, it offers a practical alternative to fully resolved simulations for multi-bubble problems.","The differences between interaction models grow with cluster size and density, so finite sound speed matters most in large, dense systems."],"supporting_citations":[{"why":"Supplies the quasi-acoustic assumption and the wave-equation/Bernoulli derivation used to build the model.","marker":"[14]"},{"why":"Provides the Keller-Miksis radial equation that the bubble dynamics section derives and solves.","marker":"[16]"},{"why":"Establishes the time-delay effect on bubbly screens and provides the 51×51 bubble-screen test case and reference patterns.","marker":"[32]"},{"why":"Gives the polydisperse microbubble-cluster dynamics and the incompressible-interaction baseline used for the frequency-response comparison.","marker":"[37]"},{"why":"Provides the Lagrangian wave-tracking algorithm for acoustic emissions that the interaction model is built on.","marker":"[48]"},{"why":"Supplies the bubble-cluster configurations (asymmetric collapse and tension-pulse response) used for validation.","marker":"[52]"},{"why":"Defines the two-bubble cavitation-onset setup and the unstable-equilibrium-radius analysis used for interpreting the onset results.","marker":"[9]"}],"fun_headline_variants":["Finite sound speed shifts cavitation onset in dense clusters","Delayed bubble interactions alter cluster resonance and onset","Time-delayed acoustic effects drive dense bubble clusters","Finite sound speed is a first-order effect for dense clusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes a linear superposition of spherically symmetric wave potentials from all bubbles and neglects wave scattering, the finite size of the receiving bubble, and non-spherical radiation; in the densest clusters these neglected effects may be as large as the time-delay corrections being studied.","fun_headline_variants_meta":{"raw":{"variants":["Finite sound speed shifts cavitation onset in dense clusters","Delayed bubble interactions alter cluster resonance and onset","Time-delayed acoustic effects drive dense bubble clusters","Finite sound speed is a first-order effect for dense clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000595,"raw_usage":{"total_tokens":2767,"prompt_tokens":908,"completion_tokens":1859,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":1804}},"tokens_in":524,"tokens_out":1859,"duration_ms":10691,"temperature":1.0,"reasoning_tokens":1804,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:09:14.531049+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"In a dense monodisperse spherical cluster (e.g., 250 bubbles of 2 µm radius in a 232 µm cluster), compare the predicted mean radius evolution and pressure at the cluster center under a 1.75 µs tension pulse against a fully resolved compressible two-phase simulation. If the resolved simulation shows no difference from an incompressible-interaction model in the onset time or peak pressure, the paper's central claim of appreciable time-delay effects in dense clusters is not supported.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the quasi-acoustic assumption and the wave-equation/Bernoulli derivation used to build the model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Keller-Miksis radial equation that the bubble dynamics section derives and solves."},{"cited_title":"Fan , author H","cited_arxiv_id":null,"evidence_quote":"Establishes the time-delay effect on bubbly screens and provides the 51×51 bubble-screen test case and reference patterns."},{"cited_title":"Maeda \\ and\\ author T","cited_arxiv_id":null,"evidence_quote":"Supplies the bubble-cluster configurations (asymmetric collapse and tension-pulse response) used for validation."},{"cited_title":"Ida ,\\ title en title Multibubble cavitation inception , \\ 10.1063/1.3265547 journal journal Physics of Fluids \\ volume 21 ,\\ pages 113302 ( year 2009 ) NoStop","cited_arxiv_id":null,"evidence_quote":"Defines the two-bubble cavitation-onset setup and the unstable-equilibrium-radius analysis used for interpreting the onset results."}],"review_version":1}