{"id":"8044ce6b-57c0-4562-95e2-746486889423","arxiv_id":"1908.08382","paper_version":3,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A master-slave coupling between a cable's beam centerline and an embedded surface mesh gives energy-conserving cable fluid-structure interaction, verified on aerial refueling and validated on Mars parachute inflation.","lead":"This paper introduces a master-slave kinematic coupling that lets a one-dimensional beam model of a cable drive an embedded surface mesh for fluid-structure interaction, and transfers fluid loads back to the beam with energy conservation. It could make simulations of parachute suspension lines, aerial refueling hoses, and offshore risers more accurate and affordable.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Rigid cross-section kinematics is the least secure link; an ovalization test for the soft hose would settle whether the 10% accuracy claim is robust.","rationale":"The reader identified exactly the assumption that is least secure: the time-independence of d_j^i in Eq. (1), which forces each cable cross-section to remain rigid. I agree with that identification. The virtual-work proof in Eq. (5) is a clean algebraic identity for the load-transfer step, and the numerical comparisons are internally consistent. The concern is not about the proof but about whether the physical kinematics it encodes are adequate for the target applications. The parachute suspension lines are very stiff, so the rigid-section assumption is safe there; the refueling hose is much softer (E=17 MPa) and is verified only against the dressing approach, which shares the same rigid-section modeling assumption. Thus the comparison cannot detect cross-section deformation errors. The proposed 2D ring/shell check is a cheap, targeted way to decide whether the assumption matters: if the hose ovalizes beyond about 1% of its diameter, the slave-surface geometry and the moment transfer (3) would be inaccurate, and the general claim would need qualification. If the check passes, the acceptance stands. Since the parachute validation and the conservation proof are not threatened by this concern, UNCHANGED is the appropriate verdict; the concern is a worthwhile sensitivity check rather than a demonstrated flaw.","tokens_in":18878,"tokens_out":16659,"duration_ms":195964,"concrete_test":"Using the hose parameters in Table 1 and the peak aerodynamic pressure from the master-slave simulation, run a 2D nonlinear ring/shell analysis of a unit-length hose cross-section (E=17 MPa, nu=0.42) to compute ovalization under that load. If the maximum radial deflection exceeds about 1% of D=0.067 m, the time-independent d_j^i assumption is violated and the moment transfer Eq. (3) needs a deformable-section correction; if it is below 1%, the assumption is confirmed for this application.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The weakest link is the rigid-cross-section kinematics encoded in Section 3: the distance vector d_j^i between a slave surface node and its master point is \"assumed to be time-independent\" (just before Eq. (1)). This means every cable cross-section translates and rotates as a rigid plane, so the embedded surface cannot ovalize, warp, or otherwise deform. The conservation proof (Eq. (5)) is correct for the load transfer given those kinematics, but it does not validate the kinematics themselves. For the parachute suspension lines (D=3.175 mm, E=29.5 GPa, braided Technora), the assumption is reasonable. The refueling hose, however, has E=17 MPa and a length-to-diameter ratio of 119; aerodynamic loading at M=0.5 could plausibly ovalize such a soft tube. The dressing-approach comparison in Section 5.1 cannot settle this, because the dressing approach uses the same beam/section model and therefore shares the same rigid-section assumption. If the true cross-section deforms, the slave surface geometry fed to the fluid solver is wrong, the moment in Eq. (3) is miscomputed, and the \"accurate\" cable-driven FSI claim is compromised for exactly the soft-cable regime the method is meant to enable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes an embedded boundary (master-slave kinematic) method for fluid-structure interaction with cable subsystems. The structural dynamics of a cable are modeled by 1D beam/cable elements along its centerline C, while the physical cable surface Sigma_h is embedded in a non-body-fitted CFD mesh. Each surface node is slaved to a master point on the centerline by assuming that the initial distance vector d_j^i is time-independent (rigid cross-section kinematics), and its motion is computed from the interpolated translation and rotation of the beam. Flow loads are computed on Sigma_h using a shifted Gauss-point quadrature and transferred to the beam nodes with a conservative load/moment transfer. The authors prove a global energy-conservation identity for this transfer and integrate the approach into the FIVER Eulerian framework, including a new local ghost-node population strategy for slender surfaces and a doubly-intersected edge AMR criterion. The method is verified against the dressing approach on an aerial refueling hose model and validated against Curiosity Mars parachute flight drag data, with an ablation showing that suspension-line FSI significantly reduces total drag despite the small direct drag of the lines.","tokens_in":19047,"tokens_out":8011,"duration_ms":85643,"significance":"Subject to the issues below, the method addresses a real and under-served problem: two-way coupling between 1D structural cable models and 3D fluids. The global conservation result (Eq. (5)) is cleanly proved and is an algebraic consequence of the chosen kinematics and load-transfer shape functions, which is a useful design property. The validation against independent flight data for the Mars DGB parachute is a strong point and gives credibility to the overall framework; the ablation of suspension-line FSI is informative. The paper also contributes a concrete algorithmic component (local ghost-node population and double-intersection AMR criterion) that is relevant beyond cable FSI. However, the physical fidelity of the rigid-cross-section kinematics is not tested independently, the accuracy order of the shear-load quadrature is not established, and the validation evidence is a single drag time-history without mesh-convergence or sensitivity studies. These gaps are fixable and do not invalidate the central algorithmic idea.","major_comments":[{"comment":"The time-independence of d_j^i imposes a rigid-cross-section kinematics: each cable section translates and rotates as a plane, so the embedded surface cannot ovalize or warp. The energy-conservation proof following Eq. (5) is correct for this kinematics but does not validate it. The comparison with the dressing approach in Section 5.1 cannot validate it either, because the dressing superelement uses rigid massless beams to connect centerline nodes to phantom surface nodes, i.e., it embodies the same rigid-section assumption. The flight-data validation in Section 5.2 concerns Technora suspension lines (E = 29.5 GPa) for which the assumption is reasonable, but the soft-hose regime (E = 17 MPa, Table 1) that motivates the cable-FSI capability remains unvalidated. Please add a test that exercises cross-section deformation (e.g., an ovalization-prone soft cylinder) or provide an asymptotic argument for when the rigid-section assumption is valid, and otherwise explicitly scope the accuracy claims to cables whose sections remain rigid.","section":"Section 3, Eq. (1); Sections 5.1-5.2"},{"comment":"The paper claims that the shear-stress contribution to the shifted-Gauss-point quadrature (12) is a second-order approximation of its counterpart in (11). The Taylor argument is given for pressure using partial p / partial n approximately 0, but no analogous boundary condition is stated for the viscous stress tensor. For a generic smooth shear-stress field, tau(G'_k) = tau(G_k) + O(h), so the shear contribution to (12) is only first-order accurate with respect to (11) unless an additional assumption such as partial tau / partial n approximately 0 is justified, which does not follow from the no-slip wall condition. Please either prove the second-order claim with the appropriate boundary condition or revise the accuracy statement; this is the load-computation step on which the master-slave transfer and the reported results depend.","section":"Section 4.2, Eq. (12) and following Taylor expansion"},{"comment":"The validation of the parachute simulation rests on a single comparison of total drag time-history with Curiosity flight data, with a claimed relative error of less than 10%. The paper does not define the error metric (time interval, norm, inclusion of the inflation transient), does not report a mesh-convergence study for the parachute mesh, and does not test sensitivity to the new doubly-intersected edge AMR criterion or to the stated mesh sizes (3 mm near suspension lines, 5 cm near the canopy). Since the central claim is that the proposed approach accurately resolves cable-driven FSI, this evidence is thinner than the claim requires. Please add at least one coarser/finer resolution comparison or a sensitivity study, and give a precise definition of the reported relative error.","section":"Section 5.2, Figure 11"}],"minor_comments":[{"comment":"The phrase 'in the the case' should be corrected to 'in the case'.","section":"Section 5.2"},{"comment":"The caption labels the top panel as the pinned end while the text describes the cross-section at the free end; please reconcile this inconsistency.","section":"Figure 7"},{"comment":"The doubly-intersected edge criterion is described only heuristically; please specify the geometric tolerance for 'intersected twice' and how the criterion is combined with the existing distance/Hessian AMR criteria in the implementation.","section":"Section 4.3"},{"comment":"The statement that the proposed approach is 'more comprehensive and user-friendly' than the dressing approach could be supported by quantitative measures (e.g., wall-clock time, condition number, DOF counts) in addition to the shown matching time histories.","section":"Section 5.1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for the journal and the core algorithmic idea is sound. The main concern is that the validation section needs strengthening before the accuracy claims can be fully supported. The rigid-section issue in particular should be addressed either with a numerical experiment or an explicit scope caveat. I do not see a citation or novelty concern: the self-citations are appropriate given the FIVER lineage."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this is a real method paper, not a repackaging. What is new is replacing the dressing approach's phantom elements and rigid-link network with direct master-slave kinematics between a 1D beam centerline and an embedded surface mesh, plus an energy-conserving load transfer that follows from a virtual-work identity. The local ghost-node update for slender embedded surfaces and the doubly-intersected-edge AMR criterion are smaller but real additions. I checked the algebra in Eq. (5); it goes through cleanly. Credit where due: the refueling comparison against the dressing approach is a sensible verification, and the Curiosity DGB parachute comparison against flight-measured total drag within 10% is meaningful external validation.\n\nSoft spots, in proportion. The paper does not fully close off the rigid cross-section assumption. Equation (1) treats the distance vector d_j^i as time-independent, so each cable section translates and rotates as a rigid disk. For the braided Technora suspension lines (D ~3 mm, E ~29.5 GPa) that is clearly reasonable. For the refueling hose (E = 17 MPa, L/D = 119) it is less obviously safe; the dressing comparison shares the same beam/section model, so it cannot detect ovalization or section deformation. That is a limitation, not a refutation—the load-transfer identity does not depend on section rigidity, and slender cables are the stated scope. Still, a sensitivity test or an explicit note on when the assumption breaks would complete the story.\n\nSecond, the Mars case validates only total drag time history. No mesh convergence or sensitivity study on the parachute problem is reported, and no error bars on the flight data. That keeps this at solid-not-definitive method validation. The paper also self-cites the FIVER machinery heavily, but that is not a flaw here; the coupling contribution is distinct and the conservation result is not fitted to the data.\n\nWho this is for: computational FSI researchers working on parachutes, refueling hoses, risers, tethers—anyone needing two-way coupling for slender cable subsystems. Worth a serious referee. My recommendation: accept for peer review, but ask for a short section on cross-section deformation limits rather than a rewrite. I would bring it to a reading group if the group cares about embedded boundary methods.","headline":"A genuine method contribution for cable-driven FSI with a clean conservation proof and a meaningful flight-data validation, though the rigid-section kinematics and single-case validation keep it from being definitive.","tokens_in":19602,"tokens_out":1891,"would_cite":true,"duration_ms":21635,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A cable's flow can be captured by slaving its embedded surface to a beam centerline, with exact energy transfer.","keywords":["cable dynamics","embedded boundary method","fluid-structure interaction","immersed boundary","master-slave kinematics","beam centerline","parachute inflation","adaptive mesh refinement"],"falsifier":"A body-fitted simulation of a short, flexible cable segment in cross-flow that resolves the cable with 3D solid elements would settle this: if the resolved cross-section measurably ovalizes, the master-slave reconstruction of the surface from centerline kinematics fails. A water-tunnel experiment measuring cross-sectional deformation of a flexible cable would test Eq. (1) directly.","tokens_in":18641,"feed_emoji":"🪂","tokens_out":9807,"duration_ms":96090,"temperature":0.7,"pith_summary":"The paper argues that cable–fluid interaction can be resolved without a three-dimensional structural mesh by representing the cable as a one-dimensional beam for its dynamics and an embedded discrete surface for its fluid boundary. This matters because cables such as parachute suspension lines and refueling hoses are long, slender, and flexible, so their structural dynamics are naturally modeled by beam elements, yet their effect on the surrounding flow can be strong. The proposed master-slave kinematics slave every surface node to the beam centerline and transfer the computed surface loads back to the beam through a virtual-work-conserving identity. The approach reproduces the dressing approach on an aerial refueling hose and, for a supersonic disk-gap-band parachute inflation, brings the predicted total drag within 10 percent relative error of the flight-measured drag.","feed_headline":"Slaved cable surfaces match Mars parachute drag within 10%","feed_subtitle":"It predicts total drag of a supersonic Mars parachute from suspension-line flow coupling.","key_machinery":"At the heart is the master-slave kinematic relation between the beam centerline C and the embedded surface Sigma_h. Each surface node S_j^i is paired with a master point M_i on the beam, and its motion is computed from the beam's displacement and rotation via u_{S_j^i} = u_{M_i} + R(theta_{M_i}) d_j^i - d_j^i and dot u_{S_j^i} = dot u_{M_i} + omega_{M_i} x R(theta_{M_i}) d_j^i. Loads are transferred back through Eq. (3), which accumulates surface forces and moments at the master point, and Eq. (4), which distributes them to beam nodes with the beam's shape functions. The identity that carries the argument is the virtual-work equality -delta W_F = delta W_S, which proves the transfer conserves energy globally. Supporting machinery includes a shifted Gauss-point quadrature for wall loads that avoids spurious oscillations on slender embedded surfaces, and a doubly-intersected edge criterion for adaptive refinement around cables.","core_discovery":"The central claim is that a cable can be represented by two coupled objects of different dimensions: a topologically 1D finite-element beam for its dynamics, and an embedded discrete surface for its fluid boundary. The surface nodes are slaved to the beam through Eq. (1), using the beam's interpolated displacement and rotation to update surface positions and velocities, with the initial distance vector assumed constant. Flow loads computed on the surface are accumulated at master points and then distributed to beam nodes through Eqs. (3)-(4). The paper proves that this transfer is globally conservative, because the virtual work of fluid tractions on the surface equals the virtual work of the generalized forces on the beam nodes. Numerical tests show that the approach matches the dressing method on a refueling hose and, for the parachute inflation, accounts for the suspension lines' disturbance of the bow shock so that the predicted total drag agrees with flight data to within 10 percent.","pith_inferences":["Because the conservation proof is algebraic, the same master-slave transfer could be dropped into other embedded-boundary solvers and should remain conserving as long as surface loads are evaluated consistently on the embedded surface.","The rigid-cross-section assumption suggests a natural extension: for cables whose cross-sections ovalize, such as very flexible hoses or braided cords, a curvature- or pressure-based correction could be added while keeping the master-slave structure of the method.","The doubly-intersected edge refinement criterion is stated for cables but is really a generic rule for any slender body whose diameter is smaller than the local mesh size; it could apply to tethers, towed arrays, or guidewires in biomedical flows."],"forward_implications":["Cable-driven FSI can be computed on a non-body-fitted fluid mesh, so the cost of resolving the flow around a slender cable is controlled by adaptive refinement rather than by maintaining a body-fitted boundary-layer mesh.","The method removes the need for massless rigid and phantom elements, so structural time integration does not have to handle singular mass matrices or differential-algebraic constraints.","For the supersonic parachute test case, including the suspension lines in the FSI changes the predicted total drag enough to bring it into the measured range; omitting them overpredicts drag.","The virtual-work identity shows that the load transfer cannot add or remove energy between the fluid and structure despite the reduced 1D structural idealization.","The mesh-adaptation criterion deliberately does not fully resolve the cable boundary layer; it resolves the cable's effect on the flow, which is what the flight-data comparison validates."],"supporting_citations":[{"why":"Introduces the dressing approach with massless phantom and rigid elements that the master-slave method is designed to replace and is used as the reference in the refueling test.","marker":"[15]"},{"why":"Supplies the load and motion transfer framework with momentum and energy conservation that Equations (3)-(4) build on.","marker":"[31]"},{"why":"Provides the embedded-boundary finite volume method with exact two-material Riemann problems in which the proposed approach is implemented.","marker":"[18]"},{"why":"Introduces the shifted Gauss-point evaluation for wall loads that removes spurious oscillations on the slender embedded surface.","marker":"[26]"},{"why":"Provides the adaptive mesh refinement machinery that the doubly-intersected edge criterion extends for slender cables.","marker":"[27]"},{"why":"Supplies the implicit-explicit staggered time-integrators used for the parachute inflation simulation.","marker":"[44]"},{"why":"Provides the flight-measured total drag of the Mars parachute used for validation.","marker":"[45]"},{"why":"Supply the aerial refueling hose model and parameters used to compare the master-slave and dressing approaches.","marker":"[11, 12, 13]"}],"fun_headline_variants":["Cable dynamics via embedded surfaces solve parachute FSI","Beam-slaved surface predicts Mars parachute drag to 10%","Embedded surfaces for cables: Mars parachute drag within 10%","Cable-fluid coupling via slaved surface nodes matches flight data","Master-slave transfer yields 10% Mars parachute drag error"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that each cable cross-section stays rigid and attached to the beam centerline, so the surface can be reconstructed from centerline kinematics alone; if the physical cross-section deforms or ovalizes, the slave surface motion and the transferred moment are wrong.","fun_headline_variants_meta":{"raw":{"variants":["Cable dynamics via embedded surfaces solve parachute FSI","Beam-slaved surface predicts Mars parachute drag to 10%","Embedded surfaces for cables: Mars parachute drag within 10%","Cable-fluid coupling via slaved surface nodes matches flight data","Master-slave transfer yields 10% Mars parachute drag error"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000524,"raw_usage":{"total_tokens":2540,"prompt_tokens":961,"completion_tokens":1579,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":1489}},"tokens_in":577,"tokens_out":1579,"duration_ms":10346,"temperature":1.0,"reasoning_tokens":1489,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:57:05.900888+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A body-fitted simulation of a short, flexible cable segment in cross-flow that resolves the cable with 3D solid elements would settle this: if the resolved cross-section measurably ovalizes, the master-slave reconstruction of the surface from centerline kinematics fails. A water-tunnel experiment measuring cross-sectional deformation of a flexible cable would test Eq. (1) directly.","supporting_citations":[{"cited_title":"Large-scale nonlinear aeroelastic computations: flutter, LCO and buffet investigations","cited_arxiv_id":null,"evidence_quote":"Introduces the dressing approach with massless phantom and rigid elements that the master-slave method is designed to replace and is used as the reference in the refueling test."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the load and motion transfer framework with momentum and energy conservation that Equations (3)-(4) build on."},{"cited_title":"FIVER: A finite volume method based on exact two-phase Riemann problems and sparse grids for multi-material flows with large density jumps","cited_arxiv_id":null,"evidence_quote":"Provides the embedded-boundary finite volume method with exact two-material Riemann problems in which the proposed approach is implemented."},{"cited_title":"A family of position-and orientation-independent embedded boundary methods for viscous flow and fluid–structure interaction problems","cited_arxiv_id":null,"evidence_quote":"Introduces the shifted Gauss-point evaluation for wall loads that removes spurious oscillations on the slender embedded surface."},{"cited_title":"Mesh adaptation framework for embedded boundary methods for computational fluid dynamics and fluid-structure interaction","cited_arxiv_id":null,"evidence_quote":"Provides the adaptive mesh refinement machinery that the doubly-intersected edge criterion extends for slender cables."},{"cited_title":"Robust and provably second-order explicit–explicit and implicit–explicit staggered time-integrators for highly non-linear compressible fluid–structure interaction problems","cited_arxiv_id":null,"evidence_quote":"Supplies the implicit-explicit staggered time-integrators used for the parachute inflation simulation."},{"cited_title":"Reconstruction of the Mars science laboratory parachute performance","cited_arxiv_id":null,"evidence_quote":"Provides the flight-measured total drag of the Mars parachute used for validation."}],"review_version":1}