{"id":"bd37dac8-d5c3-47d2-97a9-839bfb979f55","arxiv_id":"2412.19953","paper_version":1,"verdict":"REJECT","confidence":"LOW","novelty_score":1.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"An analysis of isotropic DNS and LES data finds direction-independent fractal dimensions, a result built into the isotropic datasets and reported without quantitative evidence.","lead":"This paper applies a box-counting fractal analysis to turbulence data from the JHU Turbulence Database and an OpenFOAM large-eddy simulation, then concludes that the fractal dimension is the same in all directions. The finding is predetermined by the isotropic setup, and the paper reports no quantitative dimension values, so it offers little beyond a method demonstration.","discovery_kind":"incremental","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The isotropy claim is not about energy dissipation: Eq. (6) substitutes ∂u1/∂t for ε, so the box-counting result is disconnected from Meneveau's dissipation multifractal.","rationale":"The paper's aim is to test whether DNS/LES data reproduce Meneveau's multifractal dissipation analysis; the strongest conclusion is that box-counting dimensions are direction-independent in those simulations. For that conclusion to mean anything, the analyzed scalar field must be the energy dissipation rate. Eq. (6) instead uses ∂u1/∂t, a local acceleration component with different units and no established relation to the multifractal geometry of ε. This is the weakest load-bearing point: even if every other step is internally consistent, the result describes a different field. The reader's verdict of REJECT is therefore supported, though the stated reason can be sharpened: the unsupported surrogate, not simply undersampling, is what severs the link to Meneveau's dissipation multifractal. A recomputation using ε = 2ν S_ij S_ij from the JHU database would settle whether the proxy is innocent or the source of the claimed isotropy. Until that is done, the central claim remains unverified, and the manuscript's lack of code, data, and quantitative dimension values further prevents independent confirmation.","tokens_in":5470,"tokens_out":6878,"duration_ms":70109,"concrete_test":"Re-analyze the JHU isotropic1024 velocity data using the actual dissipation rate ε = 2ν S_ij S_ij, evaluated from velocity gradients fetched from JHTDB over the same 32×32×32 mesh and time window; apply the same thresholding/box-counting procedure in the x, y, and z directions; report the log(N)-versus-log(s) curves, fitted dimensions, and fit ranges for both ε and ε′. If the ε-based dimensions differ from the ∂u1/∂t-based ones (or fail to show a plateau over available box sizes), the paper's isotropy conclusion cannot be attributed to the energy dissipation field.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's key conclusion—that simulated DNS and LES dissipation fields have direction-independent box-counting dimensions—depends entirely on the field defined in Eq. (6). There, the 'energy dissipation' surrogate is ε′ = ∂u1/∂t, normalized by a time average of velocity components. This quantity is not the dissipation rate; true local dissipation is ε = 2ν S_ij S_ij (units L^2/T^3), whereas ∂u1/∂t has units L/T^2. No derivation, empirical correlation, or resolution analysis is given to show that the multifractal geometry of ∂u1/∂t matches that of ε. Meneveau and Sreenivasan's results concern the dissipation field, so equal box-counting dimensions for ∂u1/∂t do not constitute agreement with their multifractal cascade. The issue is load-bearing because if the surrogate is rejected the central claim is about an unrelated acceleration component, not turbulent energy dissipation. The 32×32×32 sampling further limits the box-counting scaling range to roughly five box sizes (2^1…2^5), and no fitted dimensions or standard errors are reported, making the claimed direction-independence quantitatively unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript attempts to test whether the box-counting fractal dimension of the energy dissipation field is isotropic in simulated turbulence, following Meneveau and Sreenivasan's multifractal analyses. The authors use velocity data from the JHU 'isotropic1024' DNS database and from in-house OpenFOAM LES snapshots, define a dissipation surrogate in Eq. (6), compute box-counting dimensions via Eq. (9), and conclude in Section 5 that all simulated datasets exhibit isotropic characteristics with direction-independent box-counting fractal dimensions, in agreement with Meneveau's findings.","tokens_in":5769,"tokens_out":2158,"duration_ms":21921,"significance":"If the central claim were correct, the paper would offer a useful numerical check of multifractal dissipation isotropy in DNS and LES and would connect simulation practice to Meneveau-Sreenivasan experimental results. The manuscript also identifies an important problem: simulating anisotropic, realistic atmospheric cascades requires going beyond idealized isotropic datasets. However, the current presentation does not deliver on this promise. No quantitative box-counting dimensions or uncertainty estimates are reported, no code or data are made available, and the central dissipation proxy is dimensionally inconsistent with the energy dissipation rate, so the conclusion is not supported. The paper's positive contribution is limited to a statement of intent and a qualitative observation that all studied datasets are isotropic, which is largely by construction.","major_comments":[{"comment":"The quantity defined as ε′ = ∂u1/∂t is not an energy dissipation rate. Its units are L/T^2, whereas the true dissipation rate ε = 2νSijSij has units L^2/T^3. No derivation, empirical calibration, or scaling argument is provided to show that the multifractal geometry of ∂u1/∂t matches that of the dissipation field. Since Meneveau and Sreenivasan's results concern the dissipation field, the conclusion that the box-counting dimension is direction-independent for these datasets is a statement about an acceleration component, not about turbulent energy dissipation. This is load-bearing because all subsequent box-counting analyses in Sections 3.5 and 4 use this surrogate.","section":"Section 3.3, Eq. (6)"},{"comment":"The central claim that 'the box-counting fractal dimension remains consistent across all directions' is never backed by any reported numerical value, fitted slope, regression statistic, or error bar. The 32×32×32 sampling mesh described in Section 3.3 has a dynamic range of only about five box sizes, and no convergence or resolution study is presented to show that Eq. (9) yields a well-defined limit. Without reporting the actual dimensions and their uncertainties, the direction-independence conclusion is quantitatively unsupported.","section":"Section 5 and Results (Section 4)"},{"comment":"The conclusion that all simulated datasets are isotropic is partially forced by the dataset selection. The JHU dataset is explicitly the 'isotropic1024' dataset, and Section 3.4 states that the OpenFOAM simulation is a cube 'where the turbulence dissipates isotropically.' Running an isotropic dataset through a box-counting analysis and then 'discovering' isotropy is circular. The paper does not compare the DNS and LES results against any anisotropic dataset, so it cannot test Meneveau's anisotropic multifractal framework, which is the stated motivation of the work.","section":"Section 3.4 and Section 5"},{"comment":"Equations (3)-(5) are presented as 'Meneveau's multiplicative process,' but they are not the equations from Meneveau and Sreenivasan (1991). The notation is undefined (e.g., C(q), φ(q,x), χ(i,q), ψ(i)), the sums and products are not connected to the cascade construction, and these equations are never used in the subsequent analysis. This mischaracterizes the methodological foundation of the paper and should be corrected or removed.","section":"Section 3.2, Eqs. (3)-(5)"}],"minor_comments":[{"comment":"The prose is often unclear; for example, 'the one we have used applies to a more exact representation of turbulence where people use the multifractal representation' and 'vertices' instead of 'vortices' should be corrected. A careful language edit is needed throughout.","section":"Abstract and Introduction"},{"comment":"The LES equations are written without any filtering notation and without the subgrid-scale stress term, so they are formally identical to the incompressible Navier-Stokes equations. If LES is actually used, the filtered equations and the SGS model should be stated explicitly.","section":"Section 3.4, Eqs. (7)-(8)"},{"comment":"Figure captions are incomplete or refer to placeholder labels such as 'Fig. 5' and 'Fig. 6', and Figures 7-11 are not described quantitatively in the text. The reader cannot infer from the captions what is being plotted or what conclusion to draw.","section":"Figures"},{"comment":"The statement in Section 4 that 'we see fractal agreement and consistency at 0.01s' and 'at 0.40s' is not supported by any definition of 'fractal agreement' or by any error metric; please specify what quantity is being compared and how the agreement is quantified.","section":"Results and Conclusions"}],"recommendation":"reject","confidential_remarks":"The manuscript is not ready for publication in its current form. The principal quantitative claim rests on a dimensionally incorrect dissipation surrogate and is presented without any reported fractal dimensions, so the central finding cannot be checked. The dataset selection further makes the isotropy conclusion tautological. A resubmission would need to define the dissipation field properly, report actual box-counting dimensions with uncertainties over an adequate scaling range, and include at least one anisotropic comparison case. I would also note that the manuscript cites a large number of the authors' own conference papers without substantive discussion, which should be reconsidered in any future submission."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the central result doesn't hold up. The authors claim box-counting fractal dimensions are the same in all directions for isotropic turbulence simulations, in agreement with Meneveau. But the quantity they analyze, epsilon' = du1/dt in Eq. (6), is an acceleration component, not the energy dissipation rate epsilon = 2*nu*Sij*Sij. Meneveau and Sreenivasan's multifractal results are about epsilon. So even if their box-counting numbers are correct, they are measuring something else entirely, and the comparison to Meneveau is void. The stress-test note is right on the money. What's actually good: the idea of using public DNS and LES data to test directional dependence of multifractal scaling is a reasonable exercise, and the paper is transparent about its limitations, explicitly calling for experimental data in the conclusion. It also correctly notes that common simulation datasets are isotropic by construction. Soft spots: Eq. (6) is dimensionally inconsistent and no justification is given for why du1/dt should share the multifractal geometry of epsilon. Eqs. (3)-(5) misrepresent Meneveau's multiplicative process and are never used afterward. The results section reports no fractal dimension values, no error bars, and no scaling-range analysis; a 32x32x32 mesh at 0.01 spacing gives only a handful of box sizes, far too few for reliable box counting. And the isotropy conclusion is baked into the dataset selection: they used the JHU \"isotropic1024\" dataset and an OpenFOAM setup that explicitly dissipates isotropically. The paper essentially discovers what it put in. So the central claim is unsupported, and the reasoning is circular. The writing suggests the authors are still learning the tools, which is fine, but this is not ready for peer review. If they redo the analysis using a proper dissipation field, compute dimensions with uncertainties, and compare against a genuinely anisotropic dataset, it might be worth another look. My verdict: desk reject. I would not cite this and wouldn't bring it to reading group. It's a work-in-progress lab note, not a paper.","headline":"The paper's central isotropy claim is undermined by a dimensionally wrong dissipation surrogate and circular dataset selection; not ready for review.","tokens_in":664,"tokens_out":1260,"would_cite":false,"duration_ms":31981,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulated turbulence dissipation shows the same box-counting fractal dimension along every direction, so anisotropy must come from experiments.","keywords":["multifractal cascading","direct numerical simulations","large eddy simulation","atmospheric turbulence","box-counting fractal dimension","energy dissipation surrogate","isotropy","inertial-range cascade"],"falsifier":"Compute the box-counting dimension using the full dissipation tensor $\\varepsilon=2\\nu S_{ij}S_{ij}$ instead of $\\partial u_1/\\partial t$ on the same DNS field; if the direction-independence disappears, the isotropy is an artifact of the surrogate. Also repeat with $64^3$ and $128^3$ meshes: if the dimension shifts with resolution, the scaling range was not resolved.","tokens_in":5273,"feed_emoji":"🌀","tokens_out":9510,"duration_ms":85765,"temperature":0.7,"pith_summary":"This paper asks whether numerical turbulence simulations reproduce the direction-dependent, multifractal structure of energy dissipation seen in atmospheric turbulence. The authors estimate the dissipation field from the time derivative of one velocity component, then measure its box-counting fractal dimension along each spatial direction in both a direct numerical simulation and large-eddy simulations. They find that every simulated dataset is isotropic: the fractal dimension is the same along all directions. If that is right, the isotropic assumption underlying many atmospheric models is consistent with current simulations, but anisotropic multifractal corrections for real turbulence will have to come from experiments rather than from DNS or LES data alone. This matters for remote-sensing and weather-prediction metrics that depend on how dissipation structures the atmosphere.","feed_headline":"Simulated turbulence is isotropic by fractal dimension","feed_subtitle":"Fractal dimension is identical along every axis in simulated data, so multifractal anisotropy must come from experiments.","key_machinery":"The argument is carried by the box-counting fractal dimension $D=-\\lim_{s\\to0} \\log N/\\log s$ (Eq. 9), applied to the normalized dissipation surrogate $\\varepsilon'/\\langle\\varepsilon\\rangle$, where $\\varepsilon'=\\partial u_1/\\partial t$ (Eq. 6). A $32\\times32\\times32$ sampling mesh with $0.01$ spacing supplies the point set, and the multiplicative-process multifractal equations (3)-(5) frame the expected scaling behavior. The central comparison is the value of $D$ along the $x$, $y$, and $z$ directions: identical values mean the dissipation field is declared isotropic.","core_discovery":"The central discovery is that, in every dataset examined, the box-counting fractal dimension of the energy-dissipation field is the same along all three spatial directions. The authors estimate dissipation from the time derivative of one velocity component, normalize it, and apply box-counting to DNS and LES fields; both types of simulation return isotropic fractal signatures. The paper reads this as agreement with the earlier multifractal dissipation work and as a limitation of isotropic numerical datasets: anisotropic multifractal behavior must be sought in experimental data.","pith_inferences":["If the single-component time-derivative surrogate is trusted, the observed isotropy may still be a grid-symmetry effect of the $32^3$ Cartesian mesh; sampling on a rotated or non-axis-aligned mesh would test whether direction-independence is intrinsic.","A natural extension is to apply the same box-counting pipeline to experimental hot-wire or particle-image-velocimetry dissipation fields; the paper's framing implies those would show direction-dependent dimensions that the numerical datasets lack.","The coarse $0.01$ spacing likely probes only a narrow band of the inertial range, so the reported isotropic dimension could change at higher Reynolds numbers; finer DNS data would show whether isotropy persists across the full cascade."],"forward_implications":["Current DNS and LES dissipation fields, as sampled here, behave isotropically at the fractal level, so they cannot by themselves validate direction-dependent multifractal cascade models.","The box-counting dimension of the dissipation surrogate can serve as a quick diagnostic for isotropy in numerical turbulence datasets.","Future multifractal turbulence simulations will need to incorporate experimental data to capture anisotropic dissipation, as the paper explicitly concludes.","Atmospheric remote-sensing metrics built on isotropic dissipation assumptions remain consistent with these numerical datasets over the sampled scales."],"supporting_citations":[{"why":"Supplies the multifractal model of turbulent energy dissipation and the experimental methodology the numerical analysis is designed to replicate.","marker":"[8]"},{"why":"Provides the simple multiplicative cascade model whose equations the paper uses to frame the scaling of dissipation.","marker":"[9]"},{"why":"Grounds the multifractal-measure formalism in a geophysical context relevant to atmospheric cascades.","marker":"[10]"},{"why":"Supports applying multifractal analysis to fully developed turbulence and chaotic systems.","marker":"[11]"},{"why":"Earlier work in this line found the same fractal dimension in all directions, the direct precursor of the central isotropic claim.","marker":"[12]"},{"why":"Provides the box-counting fractal dimension method and equation used for the direction-wise analysis.","marker":"[13]"}],"fun_headline_variants":["Fractal dimension identical along all axes in simulated turbulence","Simulated turbulence shows isotropic fractal dimension, not anisotropic","DNS and LES yield isotropic fractal signatures, anisotropy is experimental","Multifractal anisotropy absent in simulations, must come from experiments","Box-counting fractal dimension isotropic in all simulated fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result stands or falls on treating the time derivative of one velocity component as a faithful measure of the true energy dissipation rate, and on a $32\\times32\\times32$ grid with $0.01$ spacing resolving the scales that box-counting needs; if either fails, the measured fractal dimensions could be artifacts of the method rather than properties of the turbulence.","fun_headline_variants_meta":{"raw":{"variants":["Fractal dimension identical along all axes in simulated turbulence","Simulated turbulence shows isotropic fractal dimension, not anisotropic","DNS and LES yield isotropic fractal signatures, anisotropy is experimental","Multifractal anisotropy absent in simulations, must come from experiments","Box-counting fractal dimension isotropic in all simulated fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000306,"raw_usage":{"total_tokens":1714,"prompt_tokens":863,"completion_tokens":851,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":479,"completion_tokens_details":{"reasoning_tokens":770}},"tokens_in":479,"tokens_out":851,"duration_ms":8370,"temperature":1.0,"reasoning_tokens":770,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T23:44:58.266218+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the box-counting dimension using the full dissipation tensor $\\varepsilon=2\\nu S_{ij}S_{ij}$ instead of $\\partial u_1/\\partial t$ on the same DNS field; if the direction-independence disappears, the isotropy is an artifact of the surrogate. Also repeat with $64^3$ and $128^3$ meshes: if the dimension shifts with resolution, the scaling range was not resolved.","supporting_citations":[{"cited_title":"The Multifractal Nature of Turbulent Energy Dissipation,","cited_arxiv_id":null,"evidence_quote":"Supplies the multifractal model of turbulent energy dissipation and the experimental methodology the numerical analysis is designed to replicate."},{"cited_title":"Simple Multifractal Cascade Model for Fully Developed Turbulence,","cited_arxiv_id":null,"evidence_quote":"Provides the simple multiplicative cascade model whose equations the paper uses to frame the scaling of dissipation."},{"cited_title":"Multifractal Measures, Especially for the Geophysi- cist,","cited_arxiv_id":null,"evidence_quote":"Grounds the multifractal-measure formalism in a geophysical context relevant to atmospheric cascades."},{"cited_title":"Multifractal Nature of Fully Developed Turbulence and Chaotic Systems.,","cited_arxiv_id":null,"evidence_quote":"Supports applying multifractal analysis to fully developed turbulence and chaotic systems."},{"cited_title":"Fractal and Convolutional Analysis for Deep Atmospheric Turbulence Using Machine Learning,","cited_arxiv_id":null,"evidence_quote":"Earlier work in this line found the same fractal dimension in all directions, the direct precursor of the central isotropic claim."},{"cited_title":"High-Reynolds-Number Fractal Signature of Nascent Turbulence during Transition,","cited_arxiv_id":null,"evidence_quote":"Provides the box-counting fractal dimension method and equation used for the direction-wise analysis."}],"review_version":1}