{"id":"92c1b158-89f9-4894-8a54-57c9bcb592d8","arxiv_id":"1908.04981","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For a submerged cylindrical point absorber, a fully resolved CFD model shows that a linear Cummins potential-flow model over-predicts heave and surge, misses slow drift, and that linear-theory PTO damping is suboptimal for moderate and steep waves.","lead":"The paper compares a linear potential-flow (Cummins) model with a fully resolved Navier-Stokes CFD model for a submerged cylindrical wave-energy converter with three degrees of freedom. It shows the linear model over-predicts heave and surge motions, misses slow drift, and that linear-theory damping is not optimal for steeper waves.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"CFD drift claim lacks convergence/validation: the slow-drift result that separates the models may be a numerical tank artifact.","rationale":"The reader's verdict is CONDITIONAL, and the identified wave-order mismatch is real but secondary. The central comparison between Cummins and CFD for heave/surge/pitch is credible and supported by figures and a heave grid-convergence study, so I do not recommend rejecting the paper. However, the load-bearing step for the strongest claim is not the order of the wave model (Airy vs Stokes) but the assertion that CFD 'reliably' captures slow drift. The slow drift is the one physically new result that motivates fully resolved CFD for mooring design, and it is neither validated nor shown to be grid-converged. A small low-frequency mean displacement in a 2D numerical wave tank with a close top boundary and finite damping zone can easily be contaminated by numerical reflections and boundary artifacts; the paper's own Sec. 4.4 acknowledges top-boundary artifacts. This concern does not invalidate the paper's other results—over-prediction of heave/surge amplitudes, insignificant potential-flow pitch, and the qualitative PTO efficiency trends—but it means the headline drift conclusion should be stated conditionally until a convergence and domain-sensitivity check (or an experimental or benchmark comparison) is supplied. Hence the reader's CONDITIONAL verdict stands unchanged, but for a somewhat different reason than the wave-order mismatch. The wave-order mismatch should also be noted; the proposed NWT-configuration test would help expose whether the Stokes vs Airy input contributes to the drift discrepancy.","tokens_in":24002,"tokens_out":9495,"duration_ms":109476,"concrete_test":"Re-run the 2-DOF case of Sec. 4.2 at coarse, medium, and fine grids from Table 2 (adding a finer level if feasible) and in at least two NWT configurations: (i) 8λ working zone + 2λ damping as in the paper, and (ii) 12λ working zone + 4λ damping with the top boundary raised by at least 1 m. Measure the time-averaged surge drift velocity over t = 64–70 s in each run. If drift velocity varies by more than 20% between medium and fine grids or between tank configurations, the drift claim is unreliable; if it stays within a few percent, the concern is resolved. An optional second check is comparison against a body-fitted URANS or experimental measurement of a submerged cylinder under the same waves.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim includes that the Cummins model under-estimates slow surge drift while the CFD model 'is able to capture it reliably' (abstract; Sec. 4.2, Fig. 12c). This is the key differentiator for the mooring-design conclusions. The only wave-case validation for the CFD model is heave-amplitude grid convergence (Sec. 3.6.2, Fig. 10); surge drift, pitch, and PTO power are not grid-converged or independently benchmarked. The drift shown in Fig. 12c is a small, low-frequency displacement (order mm over 6 s), precisely the kind of signal most sensitive to numerical wave-tank effects. The paper itself admits vortex structures reach the zero-pressure top boundary (Sec. 4.4), a no-slip bottom boundary layer is present, and the damping zone is finite. Prior references validate the wave tank generally, but none validates the drift of a tethered submerged body. Invoking Chakrabarti (Sec. 4.2) to explain zero potential-flow drift does not establish that the CFD drift is physical. Without drift convergence or an independent benchmark, 'reliably' is asserted, and the conclusion that linear tools are non-conservative for mooring design is not yet established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript compares a linear potential-flow model based on the Cummins equation with a fully resolved CFD model based on fictitious-domain Brinkman penalization (FD/BP) for a submerged cylindrical point absorber with heave, surge, and pitch degrees of freedom. The comparison is performed for one, two, and three DOFs, and the paper reports that the potential model over-predicts heave and surge amplitudes, produces negligible pitch, and misses a slow surge drift that the CFD model captures. The paper then uses the CFD model to vary PTO stiffness and damping, buoy density, and wave height, concluding that PTO coefficients from linear theory are suboptimal for moderate-to-steep waves, that low-density buoys increase PTO and mooring tension, and that absorption efficiency decreases with wave height. The authors present this as the first application of the FD/BP method to wave-energy devices.","tokens_in":24260,"tokens_out":8172,"duration_ms":85379,"significance":"If fully substantiated, the results would be practically significant for wave-energy converter design: linear BEM/Cummins models would be non-conservative for motion amplitude, drift, and mooring loads, and would misguide PTO tuning in steeper waves. The paper has genuine strengths: the FD/BP formulation is described in detail, the solver is based on the open-source IBAMR library, the damped-oscillator tests provide analytical validation, and a grid-convergence study is included for heave. However, the evidence for the most consequential claims is incomplete, and one ingredient of the potential-flow model is unreported, so the main comparison is not yet reproducible.","major_comments":[{"comment":"The values of the viscous drag coefficients Cx, Cy, and Cθ in the Cummins model are never given, even though these coefficients are part of the potential-flow model that is compared against CFD. Without them the comparison in Figs. 11–13 cannot be reproduced, and the claimed over-prediction of heave and surge amplitudes may depend on how these empirical coefficients were chosen. Please report the values, cite their source, and show the sensitivity of the Cummins results to them.","section":"§2.1, Eqs. (12)–(14)"},{"comment":"The statement that the CFD model 'is able to capture' slow surge drift reliably is not supported by the evidence presented. The only wave-driven grid-convergence test is for heave amplitude (Sec. 3.6.2, Fig. 10); the drift signal in Fig. 12(c) is of order millimeters over a few seconds and is exactly the kind of low-frequency quantity sensitive to numerical tank effects, including the top-boundary artifacts and bottom no-slip layer that the authors acknowledge in Sec. 4.4. In addition, the CFD is forced by fifth-order Stokes waves while the Cummins model uses first-order Airy excitation; because drift is a nonlinear wave property, part of the difference could reflect wave-model order rather than viscous versus potential modeling. Please provide drift convergence with spatial resolution, tank length, and damping-zone length, or an independent experimental or benchmark confirmation, before claiming reliable drift capture.","section":"§4.2, Fig. 12(c)"},{"comment":"The paper calls bPTO = 80.64 N·s/m the 'optimal control damping' and uses it to claim that linear-theory PTO coefficients are sub-optimal, but no optimization sweep is shown; the text states only that further increases 'did not enhance the performance significantly (data not shown).' This makes the optimality claim effectively in-sample: the value was selected from the same CFD efficiency curves used to demonstrate the improvement. At most the data support the weaker statement that a fourfold larger damping improves efficiency for the tested frequencies. Please show the bPTO sweep or an optimization procedure, or revise the claim accordingly.","section":"§5.1, Fig. 19"}],"minor_comments":[{"comment":"The unit of bPTO is written as 'N·m/s' both in the optimal-control bullet list and in the paragraph following Fig. 19; the correct unit for a translational PTO damper is N·s/m, and the text should be corrected.","section":"§5.1"},{"comment":"The caption lists 'T = 0.909 m' for the wave period; this should be seconds, not meters.","section":"Fig. 16 caption"},{"comment":"The grid-convergence statement would be more convincing if it reported the converged heave amplitude values and a convergence rate rather than a visual comparison; the visual agreement in Fig. 10 is suggestive but not quantified.","section":"§3.6.2"},{"comment":"The data availability statement points to the IBAMR repository, which contains code rather than the results behind Figs. 11–21; please either deposit the actual simulation data or clarify that code, not data, is shared.","section":"Data availability statement"},{"comment":"The high-pass filter used to remove slow drift from the CFD surge signal is not described in terms of cutoff frequency and filter order; without this information the 'drift-free' surge amplitude comparison in Fig. 12(d) is not fully reproducible.","section":"§4.2, Fig. 12(d)"}],"recommendation":"major_revision","confidential_remarks":"The claims are interesting and likely useful for the wave-energy community, but the paper currently asserts more than it demonstrates. The required fixes—reporting Cx, Cy, and Cθ, providing a drift convergence or benchmark study, and showing the PTO damping sweep—are within the scope of a revision. I have no concerns about scope fit or citation practices."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Honest take: this paper is worth a serious referee but needs a revision before I'd trust its headline claims. The central comparison—Cummins/potential flow vs FD/BP CFD for a 3-DOF submerged point absorber—is genuinely new in its specifics: first FD/BP application to a WEC, with PTO forces and torques in all three modes, and a parametric sweep of PTO damping, buoy density, and wave height. The heave and surge over-prediction by potential theory is visually supported, the damped-oscillation validation is solid, and the heave grid-convergence study (though only for heave) is a good sign. The density/mooring-tension trend and the resonance-period shift with density are useful design insights.\n\nThe soft spots are real but not evenly distributed. The biggest one is the slow-drift claim. The abstract says CFD 'captures it reliably,' but the only grid-convergence shown is heave. The drift signal is millimeters over a few seconds, exactly the kind of low-frequency response that numerical wave-tank artifacts can contaminate, and the paper itself admits vortex structures reach the top boundary and a no-slip bottom layer exists. The stress-test note says this; I agree. That section needs either drift-specific convergence or an independent benchmark before 'reliably' is supportable. The mooring-design conclusion is therefore not yet established.\n\nThe 'optimal' PTO damping bPTO = 80.64 N·s/m also bothers me. It is selected from the same CFD efficiency data used to show it outperforms the reactive-control value, and the trials supporting it are 'data not shown.' A simple efficiency-vs-bPTO plot would fix this; without it, the optimality claim is in-sample. The viscous drag coefficients Cx, Cy, Cθ in the Cummins model are never given, which makes the potential-flow model unreproducible. I'd also like the authors to note that the CFD waves are fifth-order Stokes while the potential model uses Airy waves; at H=0.01 m this is likely minor, but the drift discrepancy is itself a nonlinear property, so attributing it purely to viscosity is an unverified assumption.\n\nThese are all fixable. The paper's central comparison holds up qualitatively, and the parametric trends are worth having. I'd send it to peer review with a request for major revision—show the damping sweep, specify the drag coefficients, add drift convergence or an independent check, and soften 'reliably.' A WEC or CFD person would get value from this after that.","headline":"Useful but uneven CFD-vs-potential-flow study; the slow-drift and 'optimal' damping claims need more support before the headline conclusions hold.","tokens_in":24819,"tokens_out":3518,"would_cite":true,"duration_ms":34642,"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":"Potential-flow models for submerged wave-energy buoys overpredict heave and surge amplitudes, predict almost no pitch, and miss the slow drift that fully resolved CFD captures.","keywords":["submerged point absorber","wave energy converter","Cummins equation","potential flow theory","Brinkman penalization","fictitious domain method","slow drift","PTO damping tuning"],"falsifier":"Force both models with identical wave kinematics: feed the fifth-order Stokes free-surface elevation and velocity field into the Cummins model, or run a fully nonlinear potential-flow solver, and compare the mean surge drift and the heave and surge amplitudes. If the potential model then reproduces the CFD drift, the claimed viscous origin of the drift is wrong; if the drift gap persists, the claim is supported. A physical wave-tank experiment with the same buoy, tether stiffness, and damping, measuring mean drift over many wave periods, would settle the same question.","tokens_in":23815,"feed_emoji":"🌊","tokens_out":10100,"duration_ms":97162,"temperature":0.7,"pith_summary":"This paper tests how well a standard linear potential-flow model predicts the motion of a fully submerged, tethered cylindrical wave-energy converter that can move in heave, surge, and pitch. It claims that the potential model overpredicts the amplitudes of heave and surge, predicts essentially no pitch for the axisymmetric body, and misses the slow horizontal drift of the buoy, while a fully resolved viscous CFD model captures all three. The consequence matters for design: if true, cheap linear models are not conservative for sizing PTO hardware or mooring lines, and reliable tuning for moderate or steep waves needs viscous simulation. The paper goes on to show with CFD that the PTO damping recommended by linear resonance theory is too low for steeper waves, that lighter buoys create higher permanent tether tension, and that efficiency falls as wave height grows.","feed_headline":"Linear wave models overpredict buoy motion and miss slow drift","feed_subtitle":"Viscous CFD shows the potential-flow recipe for PTO damping is too weak for steeper waves.","key_machinery":"The argument is carried by a comparison between two motion models of the same tethered buoy. The potential-flow side is the Cummins equation, a time-domain integro-differential equation in which the fluid memory of radiated waves enters as a convolution of a radiation impulse response function with the body velocity; the paper approximates that convolution in state-space form. The CFD side is the fictitious domain Brinkman penalization method, a fully Eulerian technique that treats the solid as a region of vanishing permeability inside the fluid mesh, which lets the same PTO stiffness, damping, and torque act directly on the body. A third piece of machinery is the PTO and mooring model itself: a linear spring-damper tether that resists heave, surge, and pitch and converts the absorbed mechanical power to a damping term $b_{PTO}(d\\Delta l/dt)^2$. The comparison is completed by forcing the two models with different wave inputs, fifth-order Stokes waves in the CFD tank and first-order Airy excitation in the potential model, which is also the point where the two models are not driven identically.","core_discovery":"On the paper's own terms, the discovery is that the time-domain Cummins potential-flow model, despite being the standard fast tool for wave-energy converter design, gives a systematically different picture of a fully submerged cylindrical point absorber than fully resolved Navier-Stokes simulation. For the same buoy, tether, and wave, the potential model overpredicts heave and surge oscillation amplitudes and predicts an insignificant pitch angle, because the pressure forces on an axisymmetric body pass through its center of mass and linear theory has weak cross-mode coupling. The CFD model additionally shows a slow surge drift that the potential model under-predicts at the small wave-steepness parameter $ka$ of the test, and this drift is important for mooring design. In the CFD parametric study, the paper finds that reactive-control PTO damping $B_{33}(\\omega)$ from linear theory is sub-optimal for moderate and high wave steepness; raising the damping coefficient improves absorption efficiency at all tested frequencies. It also finds that lower buoy density raises permanent PTO and mooring tension and shifts the resonance period range, and that absorption efficiency decreases as wave height increases.","pith_inferences":["The paper does not isolate whether the slow-drift difference is caused by viscosity alone or by the different wave-order inputs, fifth-order Stokes in the CFD tank versus first-order Airy in the potential model; forcing the Cummins model with the fifth-order Stokes kinematics would settle that.","The insignificant-pitch conclusion is specific to an axisymmetric body under pressure forces that pass through the center of mass; for non-axisymmetric hulls or PTO geometries that apply torque, the pitch mode could matter even at small steepness.","An obvious testable extension is to repeat the parametric study in three dimensions or in irregular waves, where the drift and vortex-shedding patterns found here in two dimensions may change quantitatively.","The efficiency gain from raising PTO damping above $B_{33}(\\omega)$ is probably device- and sea-state-specific; a useful next step would be a response-surface sweep over stiffness and damping for each wave height to map how far the optimum moves."],"forward_implications":["If the central claim is right, linear BEM/Cummins models overestimate heave and surge amplitudes and the associated PTO power, so motion and load estimates from these models are not conservative for engineering design.","Mooring systems designed from potential-flow drift predictions would be under-sized, because the slow surge drift that sets the tether's mean load is largely invisible to linear theory.","Reactive-control PTO damping should be treated as a lower bound for moderate and steep waves; the CFD results show that a damping coefficient several times larger than $B_{33}(\\omega)$ improves absorption efficiency across the tested frequency range.","Buoy density is a design variable with structural consequences: lower density increases permanent tether tension and shifts the resonance period range, so density and PTO stiffness must be chosen together.","Absorption efficiency degrades with increasing wave height, implying either larger PTO and device sizing or deliberate efficiency-cost trade-offs for energetic sea states."],"supporting_citations":[{"why":"Supplies the Cummins equation, the time-domain potential-flow model whose predictions are compared against CFD.","marker":"[45]"},{"why":"Supplies the state-space approximation of the radiation convolution used to solve the Cummins model.","marker":"[54, 55]"},{"why":"Supplies the fifth-order Stokes wave theory used to generate waves in the CFD numerical wave tank.","marker":"[81]"},{"why":"Provides the wave classification used to argue that the test waves fall in the linear Airy regime for the potential model.","marker":"[82]"},{"why":"Provides the steady drift-force analysis used to explain why potential theory under-predicts slow surge drift at small $ka$.","marker":"[84]"},{"why":"Supplies the reactive-control tuning rule $k_{PTO}=\\omega^2(M+A_{33}(\\omega))$, $b_{PTO}=B_{33}(\\omega)$ and the multi-DOF efficiency theory the CFD results test.","marker":"[83]"},{"why":"Supplies the quadratic drag model added to the Cummins equation for the submerged cylinder.","marker":"[19]"},{"why":"Prior CFD study of a submerged cylindrical converter that the paper cites for the same finding that linear-theory PTO coefficients are sub-optimal.","marker":"[25]"}],"fun_headline_variants":["CFD reveals potential flow's blind spots for wave buoys","Linear theory misses slow drift and optimal damping for buoys","Steep waves need more damping than linear theory suggests","Lower buoy density raises tether tension, CFD finds","Wave absorption efficiency drops as wave height rises"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The comparison assumes the two models see the same wave, but the CFD model is driven by a fifth-order Stokes wave while the potential model is driven by a simpler linear wave, so part of the difference could come from the wave description rather than from viscosity.","fun_headline_variants_meta":{"raw":{"variants":["CFD reveals potential flow's blind spots for wave buoys","Linear theory misses slow drift and optimal damping for buoys","Steep waves need more damping than linear theory suggests","Lower buoy density raises tether tension, CFD finds","Wave absorption efficiency drops as wave height rises"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000703,"raw_usage":{"total_tokens":3244,"prompt_tokens":1093,"completion_tokens":2151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":2088}},"tokens_in":709,"tokens_out":2151,"duration_ms":18615,"temperature":1.0,"reasoning_tokens":2088,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:27:02.240130+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Force both models with identical wave kinematics: feed the fifth-order Stokes free-surface elevation and velocity field into the Cummins model, or run a fully nonlinear potential-flow solver, and compare the mean surge drift and the heave and surge amplitudes. If the potential model then reproduces the CFD drift, the claimed viscous origin of the drift is wrong; if the drift gap persists, the claim is supported. A physical wave-tank experiment with the same buoy, tether stiffness, and damping, measuring mean drift over many wave periods, would settle the same question.","supporting_citations":[{"cited_title":"Cummins, The impulse response function and ship motions, Tech","cited_arxiv_id":null,"evidence_quote":"Supplies the Cummins equation, the time-domain potential-flow model whose predictions are compared against CFD."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fifth-order Stokes wave theory used to generate waves in the CFD numerical wave tank."},{"cited_title":"Le M´ ehaut´ e, An introduction to hydrodynamics and water waves, Springer Science & Business Media, 2013","cited_arxiv_id":null,"evidence_quote":"Provides the wave classification used to argue that the test waves fall in the linear Airy regime for the potential model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the steady drift-force analysis used to explain why potential theory under-predicts slow surge drift at small $ka$."},{"cited_title":"Falnes, Ocean waves and oscillating systems: linear interactions including wave-energy extraction, Cambridge university press, 2002","cited_arxiv_id":null,"evidence_quote":"Supplies the reactive-control tuning rule $k_{PTO}=\\omega^2(M+A_{33}(\\omega))$, $b_{PTO}=B_{33}(\\omega)$ and the multi-DOF efficiency theory the CFD results test."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the quadratic drag model added to the Cummins equation for the submerged cylinder."},{"cited_title":"Anbarsooz, M","cited_arxiv_id":null,"evidence_quote":"Prior CFD study of a submerged cylindrical converter that the paper cites for the same finding that linear-theory PTO coefficients are sub-optimal."}],"review_version":1}