pith. sign in

arxiv: 2607.02128 · v1 · pith:HSQYTZK2new · submitted 2026-07-02 · 📡 eess.SY · cs.SY· physics.space-ph

Reachability-Based Safe-Start Regions for Approach to a Tumbling Target with Rotating LOS Constraints

Omer Burak Iskender , Keck Voon Ling , Wee Seng Lim , Erick Lansard This is my paper

Pith reviewed 2026-07-03 07:46 UTC · model grok-4.3

classification 📡 eess.SY cs.SYphysics.space-ph
keywords safe-start regiontumbling targetrotating LOS constraintsreachabilityMPC guidanceCWH dynamicssynchronization boundinner certificate
0
0 comments X

The pith

Two closed-form criteria define a conservative safe-start region for approaching a tumbling target under rotating line-of-sight constraints.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows that directional per-constraint erosion combined with a synchronization range bound produces an inner approximation to the reachable set for safe initiation of approach. This approximation supports real-time go/no-go decisions because full reachability solvers cannot run onboard at the required speed. The resulting region is validated by comparing its predictions against closed-loop Monte Carlo trials and Hamilton-Jacobi level sets on hundreds of cases. If the criteria hold, the chaser can reach the hold point, cancel apparent rotational velocity, and remain inside the time-varying LOS corridor without violating thrust limits. The work therefore supplies an explicit, computationally cheap certificate rather than relying on numerical propagation at every decision point.

Core claim

The central claim is that the synchronization set defined by the range bound r less than 2 a_max over omega_t squared together with directional erosion margins forms a sound inner certificate: any starting state inside this set yields a feasible closed-loop trajectory under the MPC law with CWH dynamics and explicit LOS corridor constraints, even though the set is a strict subset of the full positional reachable set.

What carries the argument

The pair of conservative criteria consisting of directional per-constraint erosion (margin consumed by rotation-induced drift) and the synchronization range bound r less than 2 a_max over omega_t squared.

If this is right

  • The analytical region supports onboard go/no-go decisions without invoking Hamilton-Jacobi or polytopic reachability at runtime.
  • Prediction accuracy reaches precision 0.80 and recall 0.91 against 500 closed-loop simulations.
  • The gap between the synchronization set and the full positional reachable set widens as tumble rate increases.
  • The three-regime tracking law transitions the chaser from inertial approach through co-rotation to synchronized hold while respecting the time-varying corridor.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same erosion-plus-synchronization construction could be applied to other rotating polyhedral constraints beyond LOS corridors.
  • Hardware experiments with realistic sensor noise would reveal whether the 6 percent false-positive margin remains safe under state estimation error.
  • Replacing the CWH linearization with full nonlinear relative dynamics would quantify how much the inner certificate shrinks.

Load-bearing premise

The two criteria remain conservative enough that their false-positive rate stays acceptable for the intended docking application.

What would settle it

A Monte Carlo trial in which a state inside the analytical region produces a closed-loop trajectory that violates the LOS corridor or exceeds thrust limits.

Figures

Figures reproduced from arXiv: 2607.02128 by Erick Lansard, Keck Voon Ling, Omer Burak Iskender, Wee Seng Lim.

Figure 1
Figure 1. Figure 1: Nominal, stochastic, and robust backward reachable sets on the canonical double-integrator ( [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Planar CWH (𝑛 = 1.1 × 10−3 rad/s, blue solid hull) vs. a same-order double integrator (𝑛 = 0, grey dashed hull) backward reachable set for the three BRS variants; target hold [0, 30, 0, 0] ⊤, 𝑁 = 20, 𝑑𝑡 = 2 s, 𝑎max = 0.10 m/s2 . Orbital dynamics contract the radial (𝑥𝐵) reach to ∼81% of the double integrator’s (see text). 1. Far approach (𝑟 > 𝑟sync): LVLH-frame PD track￾ing of the spiral reference with vel… view at source ↗
Figure 3
Figure 3. Figure 3: Combined safe-start region overview for six tumble rates ( [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
read the original abstract

This paper presents a reachability-aware guidance architecture for autonomous approach to a tumbling, uncooperative target under a rotating line-of-sight (LOS) docking corridor. The LOS admissible set rotates with the target body frame, producing time-varying polyhedral constraints in the chaser's relative coordinates. A safe-start region is constructed via two conservative criteria: (i) directional per-constraint erosion, the margin consumed by rotation-induced drift before thrust can arrest it, and (ii) a synchronization range bound $r < 2a_{\max}/\omega_t^2$ ensuring the chaser can cancel the apparent rotational velocity without overshooting the hold point. Closed-loop guidance uses a receding-horizon MPC controller with Clohessy-Wiltshire-Hill (CWH) prediction dynamics and explicit LOS corridor constraints in the quadratic program. Truth propagation uses the exact discrete CWH state-transition matrix with sub-stepping, so feasibility claims are physically honest: no reference blending or state projection is applied. A three-regime tracking law manages the transition from long-range inertial approach to body-frame co-rotation and synchronized hold. The analytical safe-start region is benchmarked against four standard reachability engines (backward and forward polytopic reachable sets, Hamilton-Jacobi level sets, and closed-loop Monte Carlo): the closed-form criteria are 250x faster than Hamilton-Jacobi reachability while predicting closed-loop feasibility with precision 0.80 and recall 0.91 on a 500-case sweep. The residual 6% false-positive rate and the IoU gap against Hamilton-Jacobi quantify a structural property: the synchronization set (reach and co-rotate) is a strict subset of the positional reachable set, the gap widening with tumble rate. The analytical bound is thus a sound inner certificate for onboard go/no-go decisions where Hamilton-Jacobi is prohibitively expensive.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper presents a reachability-aware guidance architecture for autonomous approach to a tumbling uncooperative target under rotating LOS docking constraints. Safe-start regions are constructed from two conservative criteria—directional per-constraint erosion and the synchronization range bound r < 2 a_max / ω_t²—followed by receding-horizon MPC with explicit CWH dynamics and LOS constraints. The analytical criteria are benchmarked against backward/forward polytopic sets, Hamilton-Jacobi level sets, and closed-loop Monte Carlo, claiming 250× speedup over HJ while achieving precision 0.80 and recall 0.91 on a 500-case sweep, with the result positioned as a sound inner certificate for onboard go/no-go decisions.

Significance. If the inner-certificate claim holds, the work supplies a fast, closed-form onboard safety filter for a practically relevant class of rendezvous problems where full numerical reachability is prohibitive. The explicit multi-engine validation, the reported speed-up factor, the quantified conservatism (IoU gap widening with tumble rate), and the honest use of exact discrete CWH propagation without projection are concrete strengths that would be useful to the community.

major comments (2)
  1. [Safe-start region construction (criteria (i) and (ii))] The synchronization range bound r < 2 a_max / ω_t² is obtained from a simplified 1D stopping-distance argument (v = r ω, distance = v² / (2 a) < r) that does not incorporate the time-varying polyhedral LOS constraints or the directional erosion term during the three-regime transition. Under CWH dynamics the conjunction of the two criteria may therefore fail to guarantee that the trajectory remains inside the rotating admissible set, which directly affects the soundness of the inner-certificate claim.
  2. [Benchmark results section] Table or figure reporting the 500-case sweep: the 6 % false-positive rate and precision 0.80 indicate that the criteria are not perfectly sound. Without a breakdown of the false-positive trajectories (e.g., whether they violate during co-rotation because the LOS coupling was omitted from the bound), it is impossible to judge whether the residual error is acceptable for safety-critical use.
minor comments (2)
  1. The three-regime tracking law is mentioned in the abstract but the transition logic, switching conditions, and how the MPC cost and constraints change across regimes should be stated explicitly for reproducibility.
  2. Notation: define a_max, ω_t, and the erosion margins at first use and ensure they appear consistently in all equations and figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback and for recognizing the practical value of the analytical safe-start criteria. We address each major comment below and outline targeted revisions to strengthen the presentation of conservatism and validation.

read point-by-point responses

  1. Referee: [Safe-start region construction (criteria (i) and (ii))] The synchronization range bound r < 2 a_max / ω_t² is obtained from a simplified 1D stopping-distance argument (v = r ω, distance = v² / (2 a) < r) that does not incorporate the time-varying polyhedral LOS constraints or the directional erosion term during the three-regime transition. Under CWH dynamics the conjunction of the two criteria may therefore fail to guarantee that the trajectory remains inside the rotating admissible set, which directly affects the soundness of the inner-certificate claim.

    Authors: We agree that the synchronization bound originates from a simplified 1D stopping-distance argument and does not explicitly fold in the polyhedral LOS constraints or couple directly with the directional erosion term. The two criteria are applied sequentially as conservative filters, with directional erosion intended to pre-compensate rotation-induced drift on each facet. No formal invariance proof for the conjunction under full time-varying CWH dynamics is supplied; the claim is positioned as an inner certificate whose conservatism is quantified by the reported IoU gap and Monte Carlo results. We will revise the manuscript to (a) state more explicitly that the region is a conservative inner approximation rather than a complete invariant set, and (b) add a short discussion of the omitted LOS-coupling effects as a source of additional conservatism. revision: partial

  2. Referee: [Benchmark results section] Table or figure reporting the 500-case sweep: the 6 % false-positive rate and precision 0.80 indicate that the criteria are not perfectly sound. Without a breakdown of the false-positive trajectories (e.g., whether they violate during co-rotation because the LOS coupling was omitted from the bound), it is impossible to judge whether the residual error is acceptable for safety-critical use.

    Authors: We will augment the benchmark section with a breakdown of the false-positive trajectories from the 500-case Monte Carlo sweep. The added analysis will classify failures by regime (long-range approach, co-rotation transition, synchronized hold), report associated tumble rates and initial conditions, and note whether violations occur inside the rotating LOS corridor. This will allow readers to assess whether the 6 % residual error remains acceptable for the intended onboard go/no-go filter use case. revision: yes

Circularity Check

0 steps flagged

No circularity detected; analytical criteria validated against independent engines

full rationale

The paper constructs two conservative analytical criteria (directional erosion and the synchronization bound r < 2a_max/ω_t²) from simplified dynamics arguments, then benchmarks the resulting safe-start region against four external reachability engines (backward/forward polytopic sets, Hamilton-Jacobi level sets, closed-loop Monte Carlo) on a 500-case sweep, reporting precision 0.80 and recall 0.91. No load-bearing derivation step reduces by construction to a fitted parameter, self-definition, or self-citation chain; the central claim is an inner approximation whose soundness is assessed via independent numerical truth models. The 6% false-positive rate and IoU gap are explicitly quantified as structural properties of the subset relation, not hidden by redefinition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method rests on standard orbital relative-motion assumptions and conservative inner approximations whose tightness is quantified only numerically.

axioms (1)
  • domain assumption Clohessy-Wiltshire-Hill equations accurately describe relative motion for the approach distances considered
    Used both for MPC prediction and for truth propagation with the exact discrete state-transition matrix.

pith-pipeline@v0.9.1-grok · 5893 in / 1264 out tokens · 20794 ms · 2026-07-03T07:46:00.802700+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

14 extracted references · 1 canonical work pages

  1. [1]

    Researchesinthelunartheory.American Journal of Mathematics, 1(1):5–26, 1878

    G.W.Hill. Researchesinthelunartheory.American Journal of Mathematics, 1(1):5–26, 1878

  2. [2]

    W. H. Clohessy and R. S. Wiltshire. Terminal guid- ance system for satellite rendezvous.Journal of the Aerospace Sciences, 27(9):653–658, 1960

    1960

  3. [3]

    Borrelli, A

    F. Borrelli, A. Bemporad, and M. Morari.Predictive Control for Linear and Hybrid Systems. Cambridge University Press, 2017

    2017

  4. [4]

    Introduction to Model Predictive Control

    F.Borrelli,M.Morari,andC.N.Jones. Introduction to Model Predictive Control. Lecture notes, ETH Zürich / UC Berkeley, 2015

    2015

  5. [5]

    Bansal, M

    S. Bansal, M. Chen, S. Herbert, and C. J. Tomlin. Hamilton-Jacobi reachability: A brief overview and recent advances. InProceedings of the 56th IEEE Conference on Decision and Control (CDC), pages 2242–2253, 2017

    2017

  6. [6]

    AIAA Education Series, 4th edition, 2018

    H.SchaubandJ.L.Junkins.AnalyticalMechanicsof Space Systems. AIAA Education Series, 4th edition, 2018

    2018

  7. [7]

    D. Q. Mayne, M. M. Seron, and S. V. Raković. Robust model predictive control of constrained lin- ear systems with bounded disturbances.Automatica, 41(2):219–224, 2005

    2005

  8. [8]

    O. B. Iskender, K. V. Ling, and V. Dubanchet. Con- straintstighteningapproachtowardsmodelpredictive control based rendezvous and docking with uncoop- erativetargets. InProceedingsofthe2018European Control Conference (ECC), pages 380–385, 2018. 9

    2018

  9. [9]

    O. B. Iskender.Model Predictive Control for Space- craft Rendezvous and Docking with Uncooperative Targets. PhD thesis, Nanyang Technological Univer- sity, Singapore, 2020

    2020

  10. [10]

    O. B. Iskender and Y. Ling. Horizon-dependent tube MPC for elliptical-orbit rendezvous under mass un- certainty. InProceedings of the 77th International AstronauticalCongress(IAC2026),2026. toappear

    2026

  11. [11]

    Sloshing awaretubeMPCwithpartialtighteningforspacecraft proximity operations

    O.B.IskenderandJ.L.RedondoGutiérrez. Sloshing awaretubeMPCwithpartialtighteningforspacecraft proximity operations. InProceedings of the 77th International Astronautical Congress (IAC 2026),

    2026

  12. [12]

    Simulation-basedLQR-trees with input and state constraints

    P.ReistandR.Tedrake. Simulation-basedLQR-trees with input and state constraints. InProceedings of the IEEE International Conference on Robotics and Automation (ICRA), pages 5504–5510, 2010

    2010

  13. [13]

    Regionofattractionestimationforfree-floatingLQR controllers

    L.Maywald,S.Vyas,M.Ben-Larbi,andM.Janković. Regionofattractionestimationforfree-floatingLQR controllers. InProceedingsoftheGermanAerospace Congress (GVW-II), 2021

    2021

  14. [14]

    Shala, S

    E. Shala, S. Vyas, M. Ben-Larbi, M. Kumar, and E. Stoll. Region of attraction estimation for free- floating systems under time-varying LQR control. arXiv:2405.06726, 2024. 10

    work page arXiv 2024