A Threshold Model for Micrometeoroid Atmospheric Entry: Filippov Dynamics, Survival Estimates, and Survivor-Only Inverse Limits

Md Shahrier Islam Arham; Min Heo; Prasun Panthi

arxiv: 2603.28785 · v2 · pith:MTT7VD2Lnew · submitted 2026-03-19 · ⚛️ physics.ao-ph · astro-ph.EP· math.DS

A Threshold Model for Micrometeoroid Atmospheric Entry: Filippov Dynamics, Survival Estimates, and Survivor-Only Inverse Limits

Md Shahrier Islam Arham , Prasun Panthi , Min Heo This is my paper

Pith reviewed 2026-05-15 07:54 UTC · model grok-4.3

classification ⚛️ physics.ao-ph astro-ph.EPmath.DS

keywords micrometeoroidsatmospheric entryablationsurvival thresholdFilippov dynamicsinverse problemAllen-Eggers trajectory

0 comments

The pith

A reduced threshold model for micrometeoroid entry recovers the classical survival scaling of critical radius as the inverse cube of entry velocity.

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

The paper establishes a compact analytical model for the thermal survival boundary of spherical micrometeoroids by treating the melting temperature as a Filippov switching surface where ablation begins only when local heating exceeds radiative cooling. This framework combines free molecular drag, an exponential atmosphere, and a surplus-heat ablation rule to determine which particles reach the surface intact. Under standard simplifying assumptions it reproduces the known approximate relation in which the smallest surviving initial radius scales as the inverse cube of entry speed. The model also sets up an inverse reconstruction from observed survivor distributions back to pre-atmospheric entry conditions and identifies entire classes of entry events that leave no survivors and therefore cannot be recovered from ground data alone. A reader would care because the approach supplies closed-form estimates of survival fractions and quantifies the information loss inherent in any survivor-only sample.

Core claim

Under the Allen-Eggers assumptions of constant radius, constant entry angle, negligible gravity during the main heating interval, and constant transport coefficients, the threshold model yields an exact radius-loss identity along the prescribed trajectory and the classical approximate survival scaling r_0^crit ~ v_0^{-3}. The switching surface T = T_m is handled as a Filippov complementarity condition, with sustained melting occurring precisely when the local heating-to-radiation ratio exceeds unity. The inverse problem is formulated via a transfer matrix from pre-atmospheric entry bins to observed survivor bins, revealing survivor-only null spaces that require external information for full

What carries the argument

The Filippov complementarity surface at the melting temperature T = T_m, which enforces the switch between heating-only and melting-plus-ablation regimes once the heating-to-radiation ratio exceeds one.

If this is right

The critical initial radius for survival decreases with increasing entry velocity according to the inverse-cube relation.
An exact analytic identity for total radius loss holds along any prescribed Allen-Eggers trajectory.
Perturbative stability estimates bound the error when the reduced model approximates the full dynamic equations.
The transfer matrix from entry conditions to survivors is rank-deficient, so some pre-atmospheric populations remain invisible in survivor-only data.

Where Pith is reading between the lines

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

The same threshold construction could be applied to other planetary atmospheres by simply changing the density profile and radiation parameters.
Combining the survivor-only inverse limits with independent constraints from radar or satellite detection of ablation trails would allow statistical recovery of the full entry population.
The framework suggests that laboratory hypervelocity experiments with controlled initial sizes and velocities could directly test the predicted radius-loss identity.

Load-bearing premise

The assumptions of constant particle radius, constant entry angle, negligible gravity during the main heating phase, and constant transport coefficients must hold to obtain both the exact radius-loss identity and the classical v_0^{-3} scaling.

What would settle it

Direct measurement of the size-frequency distribution of recovered micrometeoroids across a range of known entry velocities would falsify the model if the observed critical radius deviates significantly from the predicted inverse-cube dependence.

Figures

Figures reproduced from arXiv: 2603.28785 by Md Shahrier Islam Arham, Min Heo, Prasun Panthi.

**Figure 2.** Figure 2: FIG. 2. Fate map for silicate micrometeoroids as a function of entry speed [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Fate map for iron micrometeoroids. The upper-right quadrant (red) shows complete ablation for [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Singular value spectrum of [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. L-curve for Tikhonov regularisation. The optimal regularisation parameter [PITH_FULL_IMAGE:figures/full_fig_p013_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. Synthetic inversion test. Left [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗

read the original abstract

Micrometeoroids enter Earth's atmosphere at hypervelocity speeds and experience rapid coupling between drag, heating, radiation, melting, ablation, and deceleration. This paper develops a reduced threshold model for the thermal survival boundary of spherical micrometeoroids. The model uses free molecular drag, an exponential atmosphere, projected-area heating, full-sphere radiative cooling, and a surplus-heat ablation rule at the melting temperature. The switching surface $T=T_m$ is treated as a Filippov/complementarity surface. Sustained melting occurs when the local heating-to-radiation ratio exceeds unity. Under the additional Allen--Eggers assumptions of constant radius, constant entry angle, negligible gravity during the main heating interval, and constant transport coefficients, this threshold yields the classical approximate survival scaling $r_0^{\rm crit}\sim v_0^{-3}$. An exact radius-loss identity is obtained along the prescribed Allen--Eggers trajectory, and a perturbative stability estimate explains when this expression approximates the full reduced model. The inverse problem is formulated through a transfer matrix from pre-atmospheric entry bins to observed survivor bins. Entry bins with zero survival probability lie in the survivor-only null space and require external information for reconstruction. The framework gives a compact analytical description of threshold entry survival and identifies the information lost when only surviving particles are observed.

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.

Desk Editor's Note private letter to a colleague

This paper wraps the classical micrometeoroid survival scaling in a Filippov threshold model and adds a survivor-only inverse formulation with null-space analysis.

read the letter

The paper introduces Filippov dynamics to handle the melting threshold for micrometeoroids and derives the standard r_0^crit ~ v_0^{-3} scaling under Allen-Eggers assumptions. It also provides an exact radius-loss identity and a perturbative stability estimate. On the inverse side, it sets up a transfer matrix from entry to survivor bins and identifies the null space where survival probability is zero. This approach does a good job of giving a compact analytical description. The complementarity surface treatment for sustained melting when heating exceeds radiation is a neat technical step. The point about information lost when observing only survivors is practical for anyone trying to reconstruct pre-atmospheric populations. The soft spots are the reliance on strong assumptions like constant radius and negligible gravity, which are necessary for the exact identity but restrict the model's reach. The abstract does not show error analysis or comparisons to full simulations, so the approximation quality remains to be checked in the full text. Since the scaling itself is classical, the novelty sits mainly in the formulation rather than the result. This work is for researchers in atmospheric entry physics and meteoroid dynamics who want an analytical handle on thresholds and inverse problems. It is coherent on its own terms and shows clear thinking about the limitations of survivor data. It deserves a serious referee because the framework is new enough to warrant scrutiny and potential refinement. I would recommend sending it to peer review, focusing on whether the Filippov model adds predictive power beyond the assumptions.

Referee Report

2 major / 3 minor

Summary. The paper develops a reduced threshold model for the thermal survival boundary of spherical micrometeoroids during hypervelocity atmospheric entry. It combines free-molecular drag in an exponential atmosphere, projected-area heating, full-sphere radiative cooling, and a surplus-heat ablation rule at the melting temperature T_m, treating the T=T_m surface via Filippov/complementarity dynamics. Sustained melting is conditioned on the local heating-to-radiation ratio exceeding unity. Under the Allen-Eggers assumptions of constant radius, constant entry angle, negligible gravity during the main heating phase, and constant transport coefficients, the threshold recovers the classical scaling r_0^crit ∼ v_0^{-3} via an exact radius-loss identity along the prescribed trajectory. A perturbative stability estimate is supplied, and the inverse problem is cast as a transfer matrix from pre-atmospheric entry bins to observed survivor bins, with zero-survival bins forming a survivor-only null space that requires external information for reconstruction.

Significance. If the central derivations hold, the work supplies a compact, analytically tractable framework that recovers established scaling laws while incorporating Filippov switching and quantifying information loss in survivor-only data. The exact radius-loss identity and the explicit identification of the null space are clear strengths for inverse reconstruction of entry populations. The parameter-free character under the stated assumptions and the scoping to classical limits are positive features that could support more rigorous population estimates once validated.

major comments (2)

[§3] §3 (Allen-Eggers limit and radius-loss identity): the exact radius-loss identity and the perturbative stability estimate are load-bearing for the claim that the threshold recovers the classical scaling; however, the manuscript provides neither the explicit algebraic form of the identity nor quantitative error bounds or direct numerical comparisons against integration of the full reduced ODE system, limiting verification of the approximation's domain of validity.
[§5] §5 (inverse transfer matrix): the construction of the transfer matrix from the forward model is central to the survivor-only null-space claim, yet the paper does not specify the numerical or analytical procedure for populating the matrix elements or discuss conditioning/regularization, which is required to assess practical utility of the inverse formulation.

minor comments (3)

[Abstract] The abstract and §1 use the phrase 'survivor-only inverse limits'; this should be aligned with the later terminology of 'null space' or 'kernel' for notational consistency.
Notation for the heating-to-radiation ratio and the Filippov switching condition should be defined explicitly with an equation number at first use to aid readability.
A brief table or paragraph contrasting the model's assumptions (free-molecular regime, constant coefficients, etc.) with those in prior micrometeoroid entry literature would help situate the contribution.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation and constructive comments on our manuscript. We address each major comment below and will revise the paper accordingly to enhance clarity and verifiability.

read point-by-point responses

Referee: [§3] §3 (Allen-Eggers limit and radius-loss identity): the exact radius-loss identity and the perturbative stability estimate are load-bearing for the claim that the threshold recovers the classical scaling; however, the manuscript provides neither the explicit algebraic form of the identity nor quantitative error bounds or direct numerical comparisons against integration of the full reduced ODE system, limiting verification of the approximation's domain of validity.

Authors: We agree that the explicit algebraic form of the radius-loss identity and supporting quantitative validation were not presented in sufficient detail. The identity follows directly from integrating the surplus-heat ablation rate along the prescribed Allen-Eggers trajectory under constant-radius and constant-coefficient assumptions. In the revised manuscript we will include the full algebraic derivation, the explicit perturbative stability estimate, and direct numerical comparisons of the threshold model against full reduced-ODE integrations across a range of entry velocities and radii, together with quantitative error bounds that delineate the approximation's domain of validity. revision: yes
Referee: [§5] §5 (inverse transfer matrix): the construction of the transfer matrix from the forward model is central to the survivor-only null-space claim, yet the paper does not specify the numerical or analytical procedure for populating the matrix elements or discuss conditioning/regularization, which is required to assess practical utility of the inverse formulation.

Authors: We acknowledge that the numerical procedure for constructing the transfer-matrix elements and any discussion of conditioning or regularization were omitted. The matrix is populated by discretizing the pre-atmospheric (r0, v0) parameter space into bins and computing the forward survival probability for each bin via the threshold model. In the revision we will explicitly describe this binning and integration procedure, report the resulting matrix conditioning, and outline regularization approaches (e.g., Tikhonov) suitable for handling the survivor-only null space in practical inverse reconstructions. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained under stated assumptions

full rationale

The paper constructs the threshold model from first-principles equations (free-molecular drag, exponential atmosphere, projected-area heating, full-sphere radiation, surplus-heat ablation at T_m) and treats the T=T_m surface via Filippov dynamics. Under the explicitly scoped Allen-Eggers assumptions (constant radius, constant entry angle, negligible gravity, constant transport coefficients), the heating-to-radiation >1 criterion produces an exact radius-loss identity that recovers the known r_0^crit ~ v_0^{-3} scaling by direct integration along the prescribed trajectory. The inverse transfer matrix is built forward from the same model, and the survivor-only null space follows immediately from bins with zero survival probability. No step reduces by construction to a fitted parameter, self-definition, or load-bearing self-citation; all central results are independent consequences of the stated equations and assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 5 axioms · 0 invented entities

The central claim relies on these standard domain assumptions from atmospheric physics and dynamical systems theory, with the Filippov framework providing the mathematical structure for the threshold.

axioms (5)

domain assumption Free molecular drag regime applies
Used for the drag force calculation.
domain assumption Exponential atmosphere model
For density variation with altitude.
domain assumption Projected-area heating and full-sphere radiative cooling
For the energy balance.
domain assumption Surplus-heat ablation rule at melting temperature
For material loss when T = Tm.
domain assumption Allen-Eggers assumptions: constant radius, constant entry angle, negligible gravity, constant transport coefficients
To derive the classical scaling and exact identity.

pith-pipeline@v0.9.0 · 5556 in / 1712 out tokens · 98343 ms · 2026-05-15T07:54:02.831425+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

[1]

Absolute and relative tolerances are set to 10 −10 and 10 −8 respectively

ODE Solver The system (8) is integrated using SciPy’ssolve ivp with the RK45 adaptive method. Absolute and relative tolerances are set to 10 −10 and 10 −8 respectively. The ablation switching is regularised with a sigmoid function Hδ(T) = [1 + exp(−(T−T melt)/δ)]−1 withδ= 10 K to avoid numerical stiffness. Terminal events: •Ground impact:z= 0 •Complete ab...

work page
[2]

Each trajectory is clas- sified into one of four outcome classes based onTmax, time aboveT melt, and final mass fraction

Fate Map Computation Fate maps are computed on a grid of 40×35 points in (r0, v0) space, withr 0 ∈[1,1000]µm (log-spaced) and v0 ∈[11.2,72] km/s (log-spaced). Each trajectory is clas- sified into one of four outcome classes based onTmax, time aboveT melt, and final mass fraction. Parallel computa- tion with 8 processes completes the 2800 trajectories in a...

work page
[3]

Base sample sizeN= 8192 givesN×(2n+ 2) = 114 688 model evaluations for n= 6 parameters

Sobol Analysis Global Sobol indices are computed using the SALib library with Saltelli sampling. Base sample sizeN= 8192 givesN×(2n+ 2) = 114 688 model evaluations for n= 6 parameters. First-order indicesS i and total-order indicesS T i are estimated using the method of Saltelli et al. [8]

work page
[4]

For each of theN e = 300 entry bins, NMC = 30 trajectories are launched with parameters jit- tered uniformly within the bin bounds

K-Matrix Construction The transfer matrixKis constructed via Monte Carlo trajectory sampling. For each of theN e = 300 entry bins, NMC = 30 trajectories are launched with parameters jit- tered uniformly within the bin bounds. The matrix ele- mentK ji is the fraction of trajectories from binilanding in surface binj. Total computation: 9000 trajectories. Ap...

work page
[5]

S. G. Love and D. E. Brownlee, Icarus89, 26 (1991)

work page 1991
[6]

J. M. C. Plane, Chemical Society Reviews41, 6507 (2012)

work page 2012
[7]

V. A. Bronshten,Physics of Meteoric Phenomena, Astro- physics and Space Science Library, Vol. 104 (D. Reidel Publishing Company, Dordrecht, Netherlands, 1983)

work page 1983
[8]

C. B. Henderson, AIAA Journal14, 707 (1976)

work page 1976
[9]

A. F. Filippov,Differential Equations with Discontinu- ous Right-Hand Sides, Mathematics and Its Applications (Kluwer Academic Publishers, Dordrecht, Netherlands, 1988)

work page 1988
[10]

H. J. Allen and A. J. Eggers,A study of the motion and aerodynamic heating of ballistic missiles entering the Earth’s atmosphere, Technical Report 1381 (NACA, 1958)

work page 1958
[11]

P. C. Hansen, SIAM Review34, 561 (1992)

work page 1992
[12]

Saltelli, M

A. Saltelli, M. Ratto, T. Andres, F. Campolongo, J. Cari- boni, D. Gatelli, M. Saisana, and S. Tarantola,Global Sen- sitivity Analysis: The Primer(John Wiley & Sons, Chich- ester, UK, 2008)

work page 2008

[1] [1]

Absolute and relative tolerances are set to 10 −10 and 10 −8 respectively

ODE Solver The system (8) is integrated using SciPy’ssolve ivp with the RK45 adaptive method. Absolute and relative tolerances are set to 10 −10 and 10 −8 respectively. The ablation switching is regularised with a sigmoid function Hδ(T) = [1 + exp(−(T−T melt)/δ)]−1 withδ= 10 K to avoid numerical stiffness. Terminal events: •Ground impact:z= 0 •Complete ab...

work page

[2] [2]

Each trajectory is clas- sified into one of four outcome classes based onTmax, time aboveT melt, and final mass fraction

Fate Map Computation Fate maps are computed on a grid of 40×35 points in (r0, v0) space, withr 0 ∈[1,1000]µm (log-spaced) and v0 ∈[11.2,72] km/s (log-spaced). Each trajectory is clas- sified into one of four outcome classes based onTmax, time aboveT melt, and final mass fraction. Parallel computa- tion with 8 processes completes the 2800 trajectories in a...

work page

[3] [3]

Base sample sizeN= 8192 givesN×(2n+ 2) = 114 688 model evaluations for n= 6 parameters

Sobol Analysis Global Sobol indices are computed using the SALib library with Saltelli sampling. Base sample sizeN= 8192 givesN×(2n+ 2) = 114 688 model evaluations for n= 6 parameters. First-order indicesS i and total-order indicesS T i are estimated using the method of Saltelli et al. [8]

work page

[4] [4]

For each of theN e = 300 entry bins, NMC = 30 trajectories are launched with parameters jit- tered uniformly within the bin bounds

K-Matrix Construction The transfer matrixKis constructed via Monte Carlo trajectory sampling. For each of theN e = 300 entry bins, NMC = 30 trajectories are launched with parameters jit- tered uniformly within the bin bounds. The matrix ele- mentK ji is the fraction of trajectories from binilanding in surface binj. Total computation: 9000 trajectories. Ap...

work page

[5] [5]

S. G. Love and D. E. Brownlee, Icarus89, 26 (1991)

work page 1991

[6] [6]

J. M. C. Plane, Chemical Society Reviews41, 6507 (2012)

work page 2012

[7] [7]

V. A. Bronshten,Physics of Meteoric Phenomena, Astro- physics and Space Science Library, Vol. 104 (D. Reidel Publishing Company, Dordrecht, Netherlands, 1983)

work page 1983

[8] [8]

C. B. Henderson, AIAA Journal14, 707 (1976)

work page 1976

[9] [9]

A. F. Filippov,Differential Equations with Discontinu- ous Right-Hand Sides, Mathematics and Its Applications (Kluwer Academic Publishers, Dordrecht, Netherlands, 1988)

work page 1988

[10] [10]

H. J. Allen and A. J. Eggers,A study of the motion and aerodynamic heating of ballistic missiles entering the Earth’s atmosphere, Technical Report 1381 (NACA, 1958)

work page 1958

[11] [11]

P. C. Hansen, SIAM Review34, 561 (1992)

work page 1992

[12] [12]

Saltelli, M

A. Saltelli, M. Ratto, T. Andres, F. Campolongo, J. Cari- boni, D. Gatelli, M. Saisana, and S. Tarantola,Global Sen- sitivity Analysis: The Primer(John Wiley & Sons, Chich- ester, UK, 2008)

work page 2008