REVIEW 3 major objections 6 minor 2 references
A Modified Moving Reference Frame Method for Propeller Resolution
T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read A smoothed reference frame removed the velocity jump at propeller-domain interfaces in ship CFD.
desk verdict Clean derivation and solid integral validation, but the local-accuracy payoff is calibrated rather than independently tested. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the modified momentum equation (Eq. 19), which adds the term (∇f·u_MRF)u to the classical MRF formulation, where u_MRF = f Ω_MRF × (r−r0) is the spatially varying reference-frame velocity and f is a sigmoidal damping function (Eq. 21) that transitions from 1 near the rotating propeller surface to 0 at the domain interface. This term is the only new addition and it enforces the compatibility condition at the interface.
What would settle it
A direct comparison of the mMRF prediction against a sliding-interface reference in a configuration where the transition band cuts through a strongly non-uniform wake, with the parameters q and k fixed a priori (not fitted to the SI solution), and where the term (∇f·u_MRF)u is not small: if the disc-averaged deviation from SI is not reduced below classical MRF, the claimed benefit is a calibration artifact. In particular, the Taylor-Couette verification cannot falsify the model because u_MRF·∇f vanishes there.
Extended reading notes
Core claim
The central claim is that scaling the MRF rotation rate by a smooth scalar function f(r), which is unity near the propeller and zero at the domain interface, eliminates the velocity discontinuity that classical MRF exhibits at the rotating-stationary boundary. The key new term in the momentum equation is (∇f·u_MRF) u, which arises from the spatially varying rotation and ensures that the interface compatibility condition f Ω_MRF × u = f u_MRF·∇u is satisfied automatically where f=0. In practice, the method reduces the disc-averaged deviation from the sliding-interface reference to 58-83% of the classical MRF values across tested planes and rotation ratios, with no additional computational cost.
Load-bearing premise
The single scalar damping function f, applied to the MRF rotation and accompanied by the extra momentum term (∇f·u_MRF)u, is a physically consistent model of a partially rotating reference frame, with its shape parameters q and k chosen for convenience rather than derived from blade motion.
Editorial extensions
If this is right
- Ship self-propulsion simulations can use larger time steps and higher MRF fractions without suffering the interface artifacts that classical MRF produces, making partially rotating grid methods more reliable for hull-propeller interaction studies.
- The method offers a favorable cost-accuracy trade-off for adjoint-based optimization, since it requires fewer time steps than sliding-interface simulations while retaining better fidelity than steady MRF.
- The mMRF formulation can be applied to configurations with energy-saving devices or pre-swirl ducts, where local flow accuracy at the propeller plane matters.
- At very high MRF fractions, the method still produces some flow distortions, so fully modeled rotation should be avoided when accurate local flow prediction is required.
Reading between the lines
- The same smoothing approach could be adapted to other rotating-machine contexts, such as pumps, turbines, or wind turbines, where an MRF interface cuts through a non-uniform wake and classical MRF produces similar discontinuities.
- The scaling function f could be optimized rather than hand-tuned: the paper fixes q and k by matching the sliding-interface solution, but a calibration-free choice based on local flow gradients or turbulence length scales might be possible.
- The method's benefit likely grows with the non-uniformity of the inflow: in nearly uniform open-water flow, classical MRF already performs well, while the largest gains appear in the ship wake, suggesting the term (∇f·u_MRF)u corrects precisely the missing convection of the velocity jump.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a modified Moving Reference Frame (mMRF) method for propeller-hull interaction simulations. The rotation rate of the reference frame is scaled by a position-dependent damping function f(r) that is intended to be unity near the propeller and zero at the interface to the stationary domain. The governing equations are derived in Sec. 2.2 and Appendix A, leading to an extra momentum term (∇f·u_MRF)u in Eq. (19). The method is implemented in the RANS solver FreSCo+, verified against the analytical Taylor–Couette solution, and applied to open-water and JBC self-propulsion cases. The paper reports that mMRF reproduces integral propulsion quantities as accurately as classical MRF while markedly reducing interface discontinuities and local flow-field deviations from a sliding-interface reference, at comparable computational cost.
Significance. If the central claim holds, the method offers a practical cost-accuracy trade-off for propeller-hull interaction CFD: it retains the efficiency of MRF while reducing the interface artifacts that limit classical MRF at large MRF fractions. The derivation is self-contained, the reduction to standard MRF for f=1 is explicitly shown, and the open-water grid/time-step convergence study in Appendix B is properly documented with monotonic convergence ratios. The application to the JBC benchmark with validation against experimental integral data is a strength. However, the verification strategy leaves a gap between the governing-equation claim and the numerical evidence: the only case that exercises the new term is precisely the case in which the damping-function parameters are calibrated, so the reported local-accuracy improvement may be partly a calibration effect.
major comments (3)
- [§3.1 and §4.3] The Taylor–Couette verification in Sec. 3.1 cannot test the new term in Eq. (19), because u_MRF·∇f=0 there by axisymmetry, as the paper itself notes. The first nontrivial exercise of this term is the JBC self-propulsion case, but Sec. 4.3 selects q=0.8 and k=30 by comparing K_T and K_Q against the sliding-interface solution F01 (Table 5), and the same F01 solution is later used as the reference for the disc-averaged deviation metric in Eq. (35) and Table 8. The reported reductions R_ε≈0.58–0.83 are therefore in-sample measures: they compare a calibrated model against the calibration reference. To make the central local-accuracy claim load-bearing, the authors should either provide an independent test case in which ∇f is active and a known or experimentally measured solution exists, or demonstrate that the qualitative conclusions are insensitive to q and k over a wide range, or evaluate the local deviations against an independent reference (e.g., EFD fields) rather than only the SI solution.
- [§2.2.1, Eq. (21)] With the selected parameters q=0.8 and k=30, Eq. (21) gives f≈(1+e^6)^-1≈2.5×10^-3 at d_int=0, not zero. Hence the compatibility condition (20) is satisfied only approximately at the interface, and the abstract's claim that mMRF "restores velocity and pressure continuity" is stronger than what is implemented. The residual discontinuity is small but nonzero, and its magnitude should be quantified in the manuscript. The authors should also discuss whether forcing f to exactly zero at the interface (e.g., by using a compactly supported damping function) would change the reported improvements.
- [Eq. (19) and Appendix A] The derivation of the absolute-velocity formulation (19) is algebraically plausible, but it is presented in a condensed form that is difficult to audit. In particular, the step leading to Eq. (41) involves a rearrangement of the convective term that is not fully explained. Since the entire method rests on this equation, the authors should either expand the derivation in Appendix A to show each cancellation explicitly or provide a supplementary symbolic algebra check. This is not a request for cosmetic changes; it is needed to rule out a hidden sign or factor error in the extra ∇f term.
minor comments (6)
- [Abstract and Sec. 2.2.1] The phrase "restoring velocity and pressure continuity" should be qualified as approximate, given that f does not vanish at the interface for the chosen sigmoid parameters (see major comment).
- [Fig. 1 caption] The caption refers to the middle panel as "an unmodified simulation with a partially rotating grid"; for clarity it should state that this is classical MRF with n_MRF/n=0.5, consistent with the terminology used elsewhere in the paper.
- [Sec. 4.3, Table 5] The table reports ΔK_T and ΔK_Q as percentages but uses the notation "0.944%" without a plus sign for positive values; using a consistent signed percentage format (e.g., +0.944%) would avoid ambiguity.
- [Appendix B, Tables 9 and 10] The convergence ratios R=0.67 and R=0.40 for the grid and time-step studies are reported without a confidence interval; stating the number of significant digits and the definition of R exactly as in Eq. (44) is fine, but the authors should note that these are single-computation estimates and not based on a Richardson-extrapolation uncertainty calculation.
- [Sec. 4.5 and Figs. 21–24] The EFD data in these figures are stated to be time-averaged experimental measurements, while the numerical fields are instantaneous at a common phase (for the deviation metric) or time-averaged over one revolution (for the contour comparisons). The caption should state clearly which averaging is used in each panel, since the mixing of instantaneous and averaged fields can mislead the reader.
- [References] The reference to Durasević et al. (2022, 2023) is cited in the text as "Durasevi´c et al. (2022; 2023)" but appears in the reference list as "Durasevi´c, S., Gatin, I., Uroi´c, T., and Jasak, H."; please ensure consistent spelling and formatting.
Circularity Check
Local-accuracy claim is not independently supported: the Taylor-Couette verification cannot activate the new ∇f term, and the q,k parameters are calibrated to the same SI reference used for the local deviation metric.
-
fitted input called prediction
[Sec. 4.3 (Table 5) and Sec. 4.5 (Eq. 35, Table 8)]
"Table 5 shows that, among the tested parameter combinations, the combination q=0.8 and k=30 (F10) exhibits the best agreement with the SI solution and is therefore selected for further use."
The free parameters q and k of the sigmoid scaling function f (Eq. 21) are selected by comparing integral KT/KQ values against the sliding-interface solution of the same JBC self-propulsion setup. The later local-accuracy headline is then measured as the disc-averaged deviation from that same SI solution (Eq. 35, Table 8). The reported error reductions (R_ε ≈ 0.58–0.83) are therefore in-sample: q,k were chosen to make the mMRF results agree with SI on integral quantities, and the local comparison uses the same reference solution. This is not a fully forced prediction because the fit target is integral while the evaluation is local, but the local improvement claim is not a blind, parameter-free test.
-
other
[Sec. 3.1 (Taylor-Couette verification)]
"Due to the axial symmetry of the flow field (27), the compatibility relation (13) is maintained and uMRF·∇f=0. Therefore, the sliding interface (SI), moving reference frame (MRF), modified moving reference frame (mMRF) simulations should always produce identical velocity fields."
The only analytical verification explicitly renders the novel term (∇f·u_MRF)u inactive. The mMRF result is therefore identical to the classical MRF result by symmetry, and this exercise cannot validate the new physical contribution of the method. The first configuration in which the ∇f term actually acts is the JBC self-propulsion case, but there q and k are calibrated against the SI reference. Consequently, the central claim that mMRF reduces interface artifacts lacks an independent verification that isolates the new term from parameter calibration.
full rationale
The derivation of the modified governing equation is self-contained algebra from the stated spatially varying rotation-rate ansatz (Eqs. 14–19), so the method is not circular at the level of equation derivation. However, the paper's headline local improvement is not independently tested. The Taylor-Couette verification cannot exercise the new ∇f term, as the paper itself notes. In the JBC application, the sigmoid parameters q and k are chosen as the combination with the best agreement with the SI solution (F10 in Table 5), and the same SI solution is used as the reference for the disc-averaged local deviation (Eq. 35, Table 8). That makes the local-accuracy comparison partly a calibrated, in-sample evaluation rather than a prediction. The effect is tempered because the calibration target is integral (KT, KQ) and the advertised reduction is local, so the reduction is not forced by construction. In addition, Eq. (21) with q=0.8 and k=30 gives f≈(1+e^6)^-1≈2.5×10^-3 at d_int=0, so the implemented scaling function does not strictly vanish at the interface and Eq. (20) is only approximately satisfied; the abstract's 'restoring continuity' overstates the implemented boundary condition. Overall, the integral self-propulsion results are credible and externally anchored to EFD, but the central local artifact-reduction claim should be scored as partially circular: score 4.
Assumptions & free parameters
free parameters (3)
- q =
0.8
- k =
30
- PID gains K_P, K_I, K_D =
0.5e-2, 0.5e-5, 0.5e-4
assumptions (6)
- domain assumption Incompressible RANS with Boussinesq eddy-viscosity closure (SST k-omega) is an adequate model for propeller-hull flow.
- ad hoc to paper The propeller rotation can be decomposed into a grid-resolved part and an MRF part, with the MRF rate scaled by a position-only function f(r)=f_R(r_R).
- ad hoc to paper The extra momentum term (nabla f dot u_MRF) u in Eq. (19) is the correct physical representation of a spatially varying frame rotation.
- domain assumption Compatibility at the interface is sufficient if f=0 there (Eq. 20).
- domain assumption The thrust identity method, plus the experimental form factor and skin friction correction, gives valid propulsion factors.
- domain assumption Wave resistance from the bare-hull two-phase simulation can be added to the single-phase double-body self-propulsion simulation without propeller interaction.
invented entities (1)
-
Spatially varying reference-frame rotation scaling f(r)
Cite this review
Pith. "Pith review of A Modified Moving Reference Frame Method for Propeller Resolution." pith.science (2026). https://pith.science/paper/4XXV6U6J
@misc{pith2026260724630,
author = {Pith},
title = {Pith review of: A Modified Moving Reference Frame Method for Propeller Resolution},
year = {2026},
howpublished = {\url{https://pith.science/paper/4XXV6U6J}},
note = {Machine review of arXiv:2607.24630}
}
abstract
Accurate resolution of propeller-hull interaction is essential for predicting the self-propulsion point in ship CFD, yet motion-resolving methods such as sliding interfaces (SI) are computationally expensive, while the classical Moving Reference Frame (MRF) approach cannot capture unsteady interaction effects. Partially rotating grid methods bridge this gap by splitting the propeller rotation into a grid-resolved and an MRF component, but the abrupt transition between the rotating and stationary domains introduces discontinuities in the velocity field. This work presents a modified MRF (mMRF) formulation in which the reference-frame rotation rate is scaled by a spatially varying function that decays smoothly from unity near the propeller to zero at the domain interface, restoring velocity and pressure continuity across the boundary. The governing equations are derived and implemented in the RANS solver FreSCo$^+$, verified against the analytical Taylor--Couette solution, and applied to open-water propeller and Japan Bulk Carrier self-propulsion simulations at model scale. Both MRF and mMRF reproduce the principal integral propulsion quantities ($n$, $K_{\mathrm{T}}$, $K_{\mathrm{Q}}$, $1-t$, $1-w_{\mathrm{T}}$, $\eta_{\mathrm{R}}$) accurately, but the mMRF markedly reduces interface discontinuities and non-physical artifacts in the local flow field, particularly at large MRF fractions, at essentially the same computational cost.
Figures
Figures from the paper (21 more)
Reference graph
Works this paper leans on
-
[1]
Bakica, A., Gatin, I., Vuk ˇcevi´c, V ., Jasak, H., and Vladimir, N. (2019). Accurate assessment of ship-propulsion characteristics using CFD.Ocean Engineering, 175:149–162. Benek, J. A., Steger, J. L., and Dougherty, F. C. (1983). A flexible grid embedding technique with application to the Euler equations. InProceedings of the 6th AIAA Computational Flui...
work page 2019
-
[2]
1978 ITTC performance prediction method
International Towing Tank Conference (2021a). 1978 ITTC performance prediction method. ITTC Recommended Procedure 7.5-02-03-01.4. International Towing Tank Conference (2021b). Uncertainty analysis in CFD verification and validation, methodology and procedures. ITTC Recommended Procedure 7.5-03-01-01. Kawamura, T., Miyata, H., and Mashimo, K. (1997). Numer...
work page 2021
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.