{"id":"230c8bb9-ac5e-465a-aacd-5ed606b1d88f","arxiv_id":"2506.12630","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A vortex formed where a barchan's lee recirculation meets an obstacle's horseshoe vortex determines whether the dune bypasses, passes over, or is trapped by the obstacle.","lead":"Using grain-scale simulations of underwater dunes, this paper shows that the outcome of a dune meeting an obstacle, passing over, going around, or stopping, is governed by a vortex that forms between the dune and the obstacle. The result offers engineers and planetary scientists a concrete mechanism for how dunes behave around buildings, bridge pillars, and Martian crater rims.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Periodic boundary channeling may explain trapped case, undermining vortex mechanism","rationale":"The reader's weakest assumption identifies the spanwise periodic boundary condition as a potential confound for the trapped case. This is indeed the most load-bearing concern because the trapped case is one of the three central outcomes, and the paper's mechanism for it relies on an interaction vortex rather than on large-scale flow confinement. The authors acknowledge the channeling risk but only rebut it by comparison to experiments with a similar wall gap, not by testing an isolated obstacle. The proposed concrete test—repeating the trapped case with a much larger spanwise domain—would directly settle whether channeling or the vortex causes the trapping. This does not change the reader's conditional verdict: the paper remains promising but needs this verification before the mechanism is accepted. I therefore mark agreement_with_reader as 'agree' and recommend no change to the verdict. Other concerns, such as the lack of a fully developed barchan in the trapped case, are secondary because the vortex mechanism could still operate even if the bedform is less mature, whereas the channeling effect is a direct alternative explanation for the grain trajectories.","tokens_in":187,"tokens_out":4350,"duration_ms":61002,"concrete_test":"Run the trapped case with identical parameters but with a spanwise domain of Lz = 0.2 m (obstacle width 70 mm, free gap 65 mm on each side), using either a single obstacle with lateral walls placed far away or a periodic domain with a much larger spacing. Compare grain trajectories, trapping fraction, and vortex structure with the current 30 mm gap case. If grains no longer follow circular paths or trapping substantially decreases, channeling was the cause; if trapping and circular paths persist, the vortex mechanism is supported for a more isolated obstacle.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim for the trapped case is that a strong vortex between the dune lee recirculation and the obstacle horseshoe vortex carries grains around the obstacle. However, the trapped-case obstacle spans 70% of the channel width (70 mm out of 100 mm) with spanwise periodic boundary conditions, creating an infinite side-by-side array with only a 30 mm free gap. The authors concede this may cause channeling and accelerate particles (Section 2.2). If this channeling, rather than the dune-obstacle vortex, produces the observed circular grain paths and trapping, then the trapped case does not demonstrate the claimed mechanism for an isolated obstacle. The authors argue that experiments had a comparable 40 mm wall gap, but that does not test an isolated obstacle and may itself include wall effects. No simulation with a wider span or lateral walls was performed to isolate the effect. Thus, the causal role of the vortex in the trapped case remains unverified, directly threatening the paper's strongest claim for one of its three defining outcomes.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports grain-resolved CFD-DEM simulations of a subaqueous barchan interacting with a prismatic obstacle, in three configurations deliberately chosen from the classification map of Assis et al. (2023): bypass, pass-over, and trapped. The fluid is treated by LES and each grain by DEM, and the outputs are used to present dune morphologies, instantaneous streamlines, Q-criterion vortex visualizations, individual grain trajectories, and space-time averaged force distributions. The central mechanistic claim is that in the bypass and trapped cases a strong vortex formed by the interaction between the dune lee-side recirculation and the obstacle's horseshoe vortex is strong enough to divert the main flow and carry grains around the obstacle, while in the pass-over case this vortex is too weak. The paper concludes that the ratio of flow-disturbance strength to grain inertia sets the outcome, with strength decreasing in the order trapped, bypass, pass-over.","tokens_in":19316,"tokens_out":8294,"duration_ms":93251,"significance":"If substantiated, the paper would add a mechanistic explanation to the previously empirical barchan-obstacle outcome maps and would demonstrate a type of grain-scale force information that experiments in this setting cannot currently provide. The work has clear strengths: it uses open-source CFD-DEM tools, deposits the data and post-processing scripts, resolves the flow at grain scale, and reproduces the three experimentally reported outcome classes. The internal consistency between the flow visualizations, grain trajectories, and force maps in Figures 3, 4, 6, and 7 is a genuine contribution. However, the causal claim that the interaction vortex 'has enough strength' is currently supported only by qualitative visualizations with an arbitrary Q-criterion threshold, and the trapped case is contaminated by a periodic-domain channeling geometry. The central conclusion is therefore not yet established at the level claimed in the abstract.","major_comments":[{"comment":"The spanwise periodic boundary condition turns the 70-mm-wide obstacle in a 100-mm-wide channel into an infinite side-by-side array with only a 30-mm free gap. The authors acknowledge in Section 2.2 that this can create channeling and accelerate particles in the trapped case, but no simulation with a wider span or with lateral walls is reported. Because the trapped outcome is one of the three central results, the observed circular grain paths and trapping in Figure 6c are equally compatible with gap channeling as with the proposed dune-obstacle vortex mechanism. A wider-domain or wall-bounded simulation of the trapped configuration is needed before the vortex can be assigned causal control in that case.","section":"Section 2.2"},{"comment":"The geometric setup is internally inconsistent. The text states that the obstacle is centered in the middle of the bottom wall in both the streamwise and spanwise directions, while Figure 1's caption places the initial pile 3 cm from the CFD inlet and the text says the pile is initially 2R (about 2.9 cm) upstream of the obstacle. If the obstacle were at mid-channel, the pile would start about 17 cm upstream, which contradicts the 2R statement; if the obstacle is 2R downstream, then all three cases have similar development distance. This matters because Section 3.1 states that in the trapped case the pile did not have enough time or distance to develop into a barchan. The geometry, development distance, and the claim about an undeveloped barchan in the trapped case need to be clarified.","section":"Sections 2.2, 3.1"},{"comment":"The key quantitative assertion that the interaction vortex 'has enough strength' to deviate the main flow and carry grains is supported only by the Q=100 iso-surfaces and by visual inspection of streamlines. No measure of vortex strength is reported, such as circulation, peak vorticity, core size, or a dimensionless ratio comparing vortex-induced forces to grain inertia, and the Q-criterion threshold is arbitrary. The ordering in the Conclusions (trapped, bypass, pass-over) is therefore not backed by a quantitative metric. I ask for a quantitative vortex-strength comparison across the three cases, and ideally a direct link between that metric and the grain-force PDFs of Figure 8.","section":"Section 3.2, Figure 4"}],"minor_comments":[{"comment":"The text says y+ remained close to 1 near the wall, but Table 1 lists y+_avg values of about 10.6. Please clarify what the two quantities represent and why the average is an order of magnitude larger.","section":"Section 2.2, Table 1"},{"comment":"The figure callouts in the paragraph describing the pass-over case are confused: the text refers to 'Figure 5b' when discussing the weak vortices of the pass-over case, and then says the pass-over case is plotted at t=100 s and t=200 s in Figures 5a and 5b. The figure numbering and cross-references should be corrected.","section":"Section 3.2, Figures 4 and 5"},{"comment":"The void-fraction smoothing length lambda=3d is imposed without sensitivity analysis; since lambda is a free parameter of the method, at least one sensitivity check or a justification from the earlier validation of Lima et al. (2022) should be given.","section":"Section 2.2, Equation 9"},{"comment":"The counts of grains touching the obstacle are based on data stored every 5000 DEM time steps, as the authors note. The 0.3% versus 20% comparison should be presented explicitly as an order-of-magnitude indicator rather than as a resolved grain count.","section":"Section 3.3"},{"comment":"The comparison with the experiments of Assis et al. (2023) for the trapped case is qualitative ('the number of grains entrained further downstream seems a little higher'). A quantitative comparison, or a statement that only the outcome class is being compared, would be more appropriate.","section":"Section 3.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of JGR: Earth Surface and the simulation effort is substantial. My main concern is not novelty but over-claiming: the mechanistic conclusion rests on qualitative vortex visualization, and the trapped case conflates the obstacle with an infinite array geometry and an undeveloped granular pile. If the authors can supply a wider-domain or wall-bounded test, a quantitative vortex-strength metric, and a precise description of the obstacle location, I would be willing to support publication. I also note for the editor that the validation cases are taken from the same group's experimental classification map, so the agreement is a consistency check with a classification rather than an independent blind prediction; this is not disqualifying but should be reflected in the wording."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Colleague],\n\nQuick take on arXiv:2506.12630 (Lima et al., subaqueous barchans hitting obstacles). This is a competent CFD-DEM study that reproduces three interaction regimes – bypass, pass over, trapped – from Assis et al.'s experimental classification, and offers a mechanism: a strong vortex forms between the dune lee recirculation and the obstacle's horseshoe vortex, deviating grains around the obstacle in the bypass and trapped cases. What's new is the grain-resolved trajectories and resultant-force fields, which are impossible to get experimentally, and the Q-criterion visualizations linking flow structure to outcomes. The data is deposited openly. That part is solid.\n\nThe soft spot is the trapped case. The obstacle spans 70mm of a 100mm-wide channel with spanwise periodic boundary conditions, so the simulation is effectively an infinite side-by-side array of wide obstacles with a 30mm gap. The authors acknowledge this could channel flow and accelerate particles, then argue the experiments had a comparable 40mm gap to the walls. But that doesn't fully answer the charge: a single obstacle with solid side walls, even at 40mm, is not the same as a periodic array, and no simulation with a wider span or lateral walls was run. So the claim that the interaction vortex, rather than gap channeling, causes the trapped outcome is not actually isolated. The vortex is present and strong, and the force maps are consistent, but that's correlation, not causation. The 'enough strength' phrase is never quantified — no vortex circulation or force balance vs. grain inertia. Also, in the trapped case the dune doesn't fully develop before hitting the obstacle, so it's more a pile-obstacle interaction, a point the authors candidly mention.\n\nNone of this sinks the paper. The bypass and pass-over mechanisms are less affected by the boundary issue, and the force/trajectory data are a real resource. The main missing piece is a controlled simulation with a wider domain or lateral walls to test whether the trapped outcome persists for an isolated obstacle. A referee should ask for that, or at least a quantification of vortex strength relative to grain inertia.\n\nMy verdict: deserves a serious referee. It's a legitimate contribution to dune-obstacle dynamics, with honest limitations, and the underlying method is established. I'd want the trapped-case control before fully buying the mechanism, but the paper is clearly worth engaging with.\n\nBest,\n[You]","headline":"A solid grain-scale CFD-DEM study with a plausible vortex mechanism, but the trapped-case causal claim needs a wider-span control simulation before it holds.","tokens_in":19783,"tokens_out":2740,"would_cite":false,"duration_ms":30405,"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":"This paper claims that a strong vortex forming between a subaqueous barchan's lee face and a dune-sized obstacle decides whether grains pass over it, flow around it, or get trapped.","keywords":["subaqueous barchan dunes","dune-obstacle interaction","CFD-DEM","large-eddy simulation","horseshoe vortex","recirculation region","grain-scale force distribution","bedload transport"],"falsifier":"Repeat the trapped-case simulation with a single isolated obstacle, or with side walls far away in the spanwise direction, and check whether the strong dune-obstacle vortex still appears and whether grains still follow circular paths into the recirculation region. Trapping without the vortex would falsify the central mechanism, while trapping only in the periodic side-by-side arrangement would show the result is a channeling effect of the repeated obstacles.","tokens_in":18919,"feed_emoji":"🏜️","tokens_out":10541,"duration_ms":111218,"temperature":0.7,"pith_summary":"This paper asks why a subaqueous barchan dune sometimes climbs over a dune-sized obstacle, sometimes flows around it without touching it, and sometimes is trapped and destroyed. Using simulations that resolve the fluid near the grain scale and track every grain, the authors argue that the answer is a strong vortex that forms between the dune's lee face and the obstacle when the dune's recirculation region meets the obstacle's horseshoe vortex. In the bypass and trapped cases this vortex is strong enough to deflect the main flow and carry grains around the obstacle; in the pass-over case the obstacle is too narrow to generate a strong vortex, so central grains follow nearly straight paths over it. If correct, the result gives a mechanistic explanation for an earlier empirical classification of dune-obstacle encounters and a way to think about dunes approaching hills, crater rims, and human structures.","feed_headline":"Vortex between dune and obstacle steers sand around or traps it","feed_subtitle":"Grain-scale flow simulations show why some barchans pass over hills and others detour or get trapped.","key_machinery":"The central object is the interaction vortex between the dune's lee face and the obstacle, produced when the recirculation region downstream of the barchan crest merges with the horseshoe vortex upstream of the obstacle. A horseshoe vortex is the wrapped, arch-like vortex that forms at the base of a wall-mounted obstacle in a boundary layer. The paper visualizes these structures with the Q-criterion, the second invariant of the velocity gradient tensor, at Q = 100, and shows the resulting vortex is much stronger in the bypass and trapped cases than in the pass-over case. The numerical machinery is an Euler-Lagrange CFD-DEM model with large-eddy simulation: the fluid is solved on a grid whose resolution approaches the grain diameter, each of the 100,000 glass grains is tracked at every time step, and the resultant force on each grain is computed directly from drag, pressure gradient, deviatoric stress, and virtual mass. The strength of this vortex, relative to grain inertia, is what the argument uses to explain all three outcomes.","core_discovery":"The paper claims that the outcome of a subaqueous barchan interacting with a dune-sized obstacle is set by a strong vortex that appears in the gap between the dune's lee face and the obstacle. This vortex forms from the interaction of two familiar flow structures: the recirculation region that exists downstream of a barchan crest and the horseshoe vortex that forms upstream of a wall-mounted prism. In the bypass and trapped cases the vortex is strong enough to deviate the main flow and carry grains around the obstacle, with trapped grains following approximately circular paths into the recirculation region downstream of it; in the pass-over case the obstacle is so narrow that its horseshoe vortex is too weak, and grains in the central region travel in nearly straight lines over the obstacle. Using grain-scale LES-DEM simulations, the paper also shows the distribution of the resultant force on each grain correlates with the fluid streamlines, and reports that in the trapped case 93 percent of grains end up trapped in the recirculation region. The authors present this as the first grain-scale view of the flow and forces during these interactions.","pith_inferences":["A testable extension is to compute the circulation of the interaction vortex and compare it with the grain Stokes number; a threshold in that ratio would turn the empirical classification map into a mechanistic prediction.","If the mechanism holds for isolated obstacles, obstacle shape could be used deliberately: a wide or tall obstacle that triggers a strong horseshoe vortex would deflect sand around itself, while a narrow upright would be more likely to be overrun and buried.","The paper's own caveat suggests an immediate numerical check: run the trapped geometry with one isolated obstacle instead of the periodic side-by-side array; if the circular paths disappear, the reported trapping is partly a channeling artifact.","A saltation-resolving version of the same three geometries would test how far the vortex mechanism reaches into eolian conditions, where grain inertia is much larger than in water."],"forward_implications":["In subaqueous conditions, the outcome of a barchan-obstacle encounter should be predictable from the strength of the vortex formed where the dune's lee recirculation meets the obstacle's horseshoe vortex, relative to grain inertia.","Bypass and trapped cases will show curved or circular grain trajectories steered by that vortex, while pass-over cases will show nearly straight streamwise paths over the obstacle.","Because water-borne grains follow the flow closely, the same vortex mechanism should be much weaker for eolian dunes, where saltating grains carry more inertia; the paper explicitly warns against direct extrapolation to desert and Martian dunes.","Force maps on individual grains provide a grain-scale diagnostic: wide force distributions accompany trapping and dune destruction, while narrow distributions accompany straight pass-over motion."],"supporting_citations":[{"why":"Defined the bypass, pass-over, and trapped outcomes and the classification map that these simulations reproduce.","marker":"Assis et al. (2023)"},{"why":"Earlier obstacle experiments on two-dimensional dunes that established the cross-over versus trapped dichotomy.","marker":"Bacik et al. (2021)"},{"why":"Validated the grain-scale LES-DEM setup on isolated subaqueous barchans, providing the methodological foundation used here.","marker":"Lima et al. (2022)"},{"why":"Developed and tested the grain-scale CFD-DEM approach for subaqueous barchans and for force distributions within them.","marker":"Alvarez and Franklin (2020, 2021)"},{"why":"Characterized the horseshoe vortex around surface-mounted prismatic obstacles, used to interpret the vortex upstream of the obstacle.","marker":"Martinuzzi and Tropea (1993)"},{"why":"Showed that horseshoe-vortex formation depends on obstacle aspect ratio and incoming velocity profile, explaining why the narrow prism has no strong vortex.","marker":"Baines (1963)"},{"why":"Established the recirculation region downstream of barchan crests that merges with the horseshoe vortex.","marker":"Kroy et al. (2002a, 2005)"}],"fun_headline_variants":["Grain-scale flow shows vortex decides if sand passes, bypasses, or traps","Dune-obstacle vortex steers grains: pass, detour, or trap","Simulation reveals vortex carries barchan sand around or into trap","Vortex between dune and obstacle controls sand grain paths","First grain-scale view: vortex routes barchan grains past or into obstacle"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations model the obstacle as one in an endless row of identical obstacles, and in the widest-obstacle case the narrow gaps between copies may channel and accelerate grains; if that channeling causes the circular trapping paths, the vortex mechanism would not hold for a single isolated obstacle in nature.","fun_headline_variants_meta":{"raw":{"variants":["Grain-scale flow shows vortex decides if sand passes, bypasses, or traps","Dune-obstacle vortex steers grains: pass, detour, or trap","Simulation reveals vortex carries barchan sand around or into trap","Vortex between dune and obstacle controls sand grain paths","First grain-scale view: vortex routes barchan grains past or into obstacle"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000673,"raw_usage":{"total_tokens":3079,"prompt_tokens":975,"completion_tokens":2104,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":2006}},"tokens_in":591,"tokens_out":2104,"duration_ms":17674,"temperature":1.0,"reasoning_tokens":2006,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:44:52.734423+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the trapped-case simulation with a single isolated obstacle, or with side walls far away in the spanwise direction, and check whether the strong dune-obstacle vortex still appears and whether grains still follow circular paths into the recirculation region. Trapping without the vortex would falsify the central mechanism, while trapping only in the periodic side-by-side arrangement would show the result is a channeling effect of the repeated obstacles.","supporting_citations":[{"cited_title":", Borges, D S","cited_arxiv_id":null,"evidence_quote":"Defined the bypass, pass-over, and trapped outcomes and the classification map that these simulations reproduce."},{"cited_title":", Canizares, P","cited_arxiv_id":null,"evidence_quote":"Earlier obstacle experiments on two-dimensional dunes that established the cross-over versus trapped dichotomy."},{"cited_title":", Assis, W R","cited_arxiv_id":null,"evidence_quote":"Validated the grain-scale LES-DEM setup on isolated subaqueous barchans, providing the methodological foundation used here."},{"cited_title":"\\ Franklin, E M","cited_arxiv_id":null,"evidence_quote":"Developed and tested the grain-scale CFD-DEM approach for subaqueous barchans and for force distributions within them."},{"cited_title":"APACrefauthors \\ 1963","cited_arxiv_id":null,"evidence_quote":"Showed that horseshoe-vortex formation depends on obstacle aspect ratio and incoming velocity profile, explaining why the narrow prism has no strong vortex."},{"cited_title":", Fischer, S","cited_arxiv_id":null,"evidence_quote":"Established the recirculation region downstream of barchan crests that merges with the horseshoe vortex."}],"review_version":1}