m physical modeling / IICSM

PHYSICS / ENGINEERING / COMPUTING

Spallation calculation for satellite survivability

Connect stress-wave failure to spacecraft debris containment, test-based damage criteria and mission-level reliability.

Subject library · 51 guides · derivations & worked examples

Fundamental physics

Space and flight

Materials and transport

Systems and mathematics

Marine and road vehicles

Electrical and computing

Companion site

Matter pathway: atom → solid → liquid → gas → plasma. Quantum mechanics and quantum field theory provide foundations across the pathway; they are not additional phases. This is a connected modeling route, not a universal heating curve. Actual phases depend on pressure, composition, and kinetics.

A satellite protection workflow

  1. Define the orbit, exposure duration and meteoroid/orbital-debris environment. NASA MEM and the orbital-debris environment are distinct inputs.
  2. Define failure for each component: detached rear-wall material, perforation, pressure leakage, electrical interruption or optical damage.
  3. Select a shield family and an experimentally supported ballistic-limit model within its stated material, velocity, angle and geometry range.
  4. Use dynamic material characterization and hydrocodes to investigate mechanisms; compare against witness plates, free-surface histories and recovered damage.
  5. Propagate uncertainty and inspect shared-fragment and common-cause failures before crediting redundancy.

Backing fabrics and catchers can absorb or redistribute debris, but their performance depends on construction, temperature, vacuum exposure, aging and mounting. Generic anti-spall marketing or a coating name is not a spacecraft qualification. This guide does not supply personnel-incapacitation coefficients or weapon-specific behind-armor models.

NASA impact physics and damage criteria · NASA meteoroid/debris shielding report · NASA meteoroid environment

1. Stress waves and tensile failure

Definitions & inputs. σ tensile-positive stress,ρ density,c wave speed,Z acoustic impedance,σsp measured dynamic spall strength.

  1. Linear acoustics links particle-velocity change and stress amplitude.

    Z=ρc,∣Δσ∣≃Z∣Δv∣Z=\rho c,\quad |\Delta\sigma|\simeq Z|\Delta v|

  2. A free surface reflects a compressive stress increment as a release/tensile increment; total traction at the surface remains zero.

    Rσ=(Z2−Z1)/(Z2+Z1)→−1(Z2→0)R_\sigma=(Z_2-Z_1)/(Z_2+Z_1)\rightarrow-1\quad(Z_2\rightarrow0)

  3. Interacting release histories can place interior material in tension, initiating voiding or cracking.

    σtensile,max>σspis a screening criterion\sigma_{tensile,max}>\sigma_{sp}\quad\text{is a screening criterion}

Interpretation. A threshold alone does not predict detached mass or debris speed. Spall strength depends on microstructure, pulse duration, strain rate, temperature and measurement interpretation; not every free-surface reflection creates spall.

↑ Return to definitions and contents

2. Damage models and hydrocode consistency

Definitions & inputs. D scalar damage0to1,σeff effective undamaged stress,εp plastic strain,α internal variables.

  1. Scalar degradation shows how damage can reduce load-carrying capacity; tension/compression treatment needs a specified law.

    σ=(1−D)σeff,0≤D≤1\sigma=(1-D)\sigma_{eff},\quad0\le D\le1

  2. Evolution must be calibrated to dynamic tests rather than inferred from static tensile strength.

    D˙=F(σ,ϵ˙,T,α)\dot D=\mathcal F(\sigma,\dot\epsilon,T,\alpha)

  3. A hydrocode evolves mass, momentum, energy and constitutive variables consistently.

    ∂tU+∇⋅F(U)=S(U,α)\partial_tU+\nabla\cdot F(U)=S(U,\alpha)

Interpretation. Grady–Kipp families connect flaw populations, activation and fragmentation; Cochran–Banner-type models use tensile deformation and damage evolution. Names do not define a unique implementation. Mesh regularization, energy dissipation, EOS, strength and failure parameters must be documented and verified against tests.

↑ Return to definitions and contents

3. Debris statistics and containment

Definitions & inputs. m fragment mass,N greater cumulative count,μ scale,n shape,N0 normalization,A exposed area,Φ environmental flux.

  1. A stretched-exponential cumulative model can summarize fragment counts above a threshold.

    N>m=N0exp⁡[−(m/μ)n]N_{>m}=N_0\exp[-(m/\mu)^n]

  2. μ is a scale, not generally the mean: for n1 it equals the mean; for n1/2 the mean is2μ.

    ⟨m⟩=μ Γ(1+1/n)\langle m\rangle=\mu\,\Gamma(1+1/n)

  3. An integrated directional environment may be simplified to this flux–area–time count only with consistent projected-area conventions.

    Nimpact=ΦAΔtN_{impact}=\Phi A\Delta t

Interpretation. Cone shapes, fragment velocity and damage cannot be transferred from unrelated armor tests or assigned a universal fraction of sound speed. Use witness plates, velocity diagnostics, recovered fragments and validated spacecraft-specific debris models.

↑ Return to definitions and contents

4. From shield damage to mission survival

Definitions & inputs. λfail expected mission failure count,pf conditional component failure given an event,λevent expected impact count.

  1. Integrate over size, velocity, direction, configuration and time using a calibrated failure criterion.

    λfail=∫Φ(z)A(z)pf(z) dz dt\lambda_{fail}=\int\Phi(z)A(z)p_f(z)\,dz\,dt

  2. Poisson thinning produces a survival probability only for the defined event process.

    P(no modeled failure)=e−λfailP(\text{no modeled failure})=e^{-\lambda_{fail}}

  3. A simplified bookkeeping example separates environment frequency from conditional damage.

    λfail=λeventpffor constant pf\lambda_{fail}=\lambda_{event}p_f\quad\text{for constant }p_f

Interpretation. Distinguish cratering, attached damage, detached spall, perforation, leakage and loss of function. Whipple and stuffed-Whipple shields, internal catchers, equipment placement and redundancy mitigate different outcomes; thermal blankets or coatings are not automatically qualified spall liners.

↑ Return to definitions and contents

Graphical worked example

Poisson mission survival for a defined component failure process; input rate needs validated environment and damage models. X axis: Expected modeled failure count (dimensionless). Y axis: Probability of no modeled failure.
Poisson mission survival for a defined component failure process; input rate needs validated environment and damage models. Related worked calculation · Download SVG · Plot data

Twenty worked examples

Open a problem to see its defined inputs, assumptions, equation, numerical substitution, result, and interpretation. Values are illustrative analytical exercises.

Example 01. One-way wave time

Definitions & inputs. L.01m,c5000m/s.

  1. Choose the governing model and isolate the requested quantity.

    t=L/ct=L/c

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    .01/5000.01/5000

  3. Evaluate the expression; the result uses the units shown.

    Result=2×10−6 s\mathrm{Result}=2\times10^{-6}\ \mathrm s

Interpretation. Elastic-wave timing only.

↑ Return to definitions and contents
Example 02. Round-trip time

Definitions & inputs. Same path.

  1. Choose the governing model and isolate the requested quantity.

    t=2L/ct=2L/c

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    2(.01)/50002(.01)/5000

  3. Evaluate the expression; the result uses the units shown.

    Result=4×10−6 s\mathrm{Result}=4\times10^{-6}\ \mathrm s

Interpretation. Reflection timing benchmark.

↑ Return to definitions and contents
Example 03. Acoustic impedance

Definitions & inputs. ρ2700kg/m³,c5000m/s.

  1. Choose the governing model and isolate the requested quantity.

    Z=ρcZ=\rho c

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    2700(5000)2700(5000)

  3. Evaluate the expression; the result uses the units shown.

    Result=1.35×107 Pa s/m\mathrm{Result}=1.35\times10^{7}\ \mathrm{Pa\,s/m}

Interpretation. Illustrative material.

↑ Return to definitions and contents
Example 04. Free-surface stress reflection

Definitions & inputs. Z2=0,Z1positive.

  1. Choose the governing model and isolate the requested quantity.

    R=(0−Z1)/(0+Z1)R=(0-Z_1)/(0+Z_1)

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    −1-1

  3. Evaluate the expression; the result uses the units shown.

    Result=−1 \mathrm{Result}=-1\ {}

Interpretation. Stress increment reverses sign.

↑ Return to definitions and contents
Example 05. Matched impedance

Definitions & inputs. Z2=Z1.

  1. Choose the governing model and isolate the requested quantity.

    R=(Z1−Z1)/(2Z1)R=(Z_1-Z_1)/(2Z_1)

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    00

  3. Evaluate the expression; the result uses the units shown.

    Result=0 \mathrm{Result}=0\ {}

Interpretation. Ideal interface has no reflection.

↑ Return to definitions and contents
Example 06. Higher-impedance reflection

Definitions & inputs. Z2=2Z1.

  1. Choose the governing model and isolate the requested quantity.

    R=(2−1)/(2+1)R=(2-1)/(2+1)

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    1/31/3

  3. Evaluate the expression; the result uses the units shown.

    Result=0.3333333 \mathrm{Result}=0.3333333\ {}

Interpretation. Stress reflection coefficient.

↑ Return to definitions and contents
Example 07. Damage stiffness fraction

Definitions & inputs. D.2.

  1. Choose the governing model and isolate the requested quantity.

    Eeff/E=1−DE_{eff}/E=1-D

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    1−.21-.2

  3. Evaluate the expression; the result uses the units shown.

    Result=0.8 \mathrm{Result}=0.8\ {}

Interpretation. Illustrative degradation law only.

↑ Return to definitions and contents
Example 08. Residual stiffness

Definitions & inputs. E70GPa,D.2.

  1. Choose the governing model and isolate the requested quantity.

    Eeff=(1−D)EE_{eff}=(1-D)E

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    .8(70).8(70)

  3. Evaluate the expression; the result uses the units shown.

    Result=56 GPa\mathrm{Result}=56\ \mathrm{GPa}

Interpretation. No compressive closure model included.

↑ Return to definitions and contents
Example 09. Fragment threshold fraction

Definitions & inputs. m=μ,n1.

  1. Choose the governing model and isolate the requested quantity.

    N>m/N0=e−1N_{>m}/N_0=e^{-1}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    e−1e^{-1}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.3678794 \mathrm{Result}=0.3678794\ {}

Interpretation. Distribution identity,no predictive material calibration.

↑ Return to definitions and contents
Example 10. Two-scale threshold fraction

Definitions & inputs. m2μ,n1.

  1. Choose the governing model and isolate the requested quantity.

    S=e−2S=e^{-2}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    e−2e^{-2}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.1353353 \mathrm{Result}=0.1353353\ {}

Interpretation. Exponential tail illustration.

↑ Return to definitions and contents
Example 11. Half-shape threshold

Definitions & inputs. m4μ,n.5.

  1. Choose the governing model and isolate the requested quantity.

    S=e−4S=e^{-\sqrt4}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    e−2e^{-2}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.1353353 \mathrm{Result}=0.1353353\ {}

Interpretation. Same formula,different parameter meaning.

↑ Return to definitions and contents
Example 12. Mean-to-scale ratio

Definitions & inputs. Shape n.5.

  1. Choose the governing model and isolate the requested quantity.

    ⟨m⟩/μ=Γ(3)\langle m\rangle/\mu=\Gamma(3)

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    22

  3. Evaluate the expression; the result uses the units shown.

    Result=2 \mathrm{Result}=2\ {}

Interpretation. μis not mean in this case.

↑ Return to definitions and contents
Example 13. Exponential median-to-scale

Definitions & inputs. n1,S.5.

  1. Choose the governing model and isolate the requested quantity.

    m50/μ=ln⁡2m_{50}/\mu=\ln2

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    ln⁡2\ln2

  3. Evaluate the expression; the result uses the units shown.

    Result=0.6931472 \mathrm{Result}=0.6931472\ {}

Interpretation. Cumulative survival convention.

↑ Return to definitions and contents
Example 14. Exposure count

Definitions & inputs. Flux.01events/(m²yr),area2m²,duration5yr.

  1. Choose the governing model and isolate the requested quantity.

    λ=ΦAt\lambda=\Phi At

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    .01(2)(5).01(2)(5)

  3. Evaluate the expression; the result uses the units shown.

    Result=0.1 \mathrm{Result}=0.1\ {}

Interpretation. Assumed direction-integrated flux.

↑ Return to definitions and contents
Example 15. Zero-event probability

Definitions & inputs. λ.1.

  1. Choose the governing model and isolate the requested quantity.

    P0=e−λP_0=e^{-\lambda}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    e−.1e^{-.1}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.9048374 \mathrm{Result}=0.9048374\ {}

Interpretation. Independent Poisson model.

↑ Return to definitions and contents
Example 16. At-least-one event

Definitions & inputs. λ.1.

  1. Choose the governing model and isolate the requested quantity.

    P≥1=1−e−λP_{\ge1}=1-e^{-\lambda}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    1−e−.11-e^{-.1}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.09516258 \mathrm{Result}=0.09516258\ {}

Interpretation. Event is not automatically failure.

↑ Return to definitions and contents
Example 17. Expected component failures

Definitions & inputs. λevent.1,pf.02from assumed test model.

  1. Choose the governing model and isolate the requested quantity.

    λf=λepf\lambda_f=\lambda_ep_f

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    .1(.02).1(.02)

  3. Evaluate the expression; the result uses the units shown.

    Result=0.002 \mathrm{Result}=0.002\ {}

Interpretation. Hypothetical conditional probability.

↑ Return to definitions and contents
Example 18. Component survival

Definitions & inputs. λfailure.002.

  1. Choose the governing model and isolate the requested quantity.

    Ps=e−λfP_s=e^{-\lambda_f}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    e−.002e^{-.002}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.998002 \mathrm{Result}=0.998002\ {}

Interpretation. Limited to defined Poisson failure process.

↑ Return to definitions and contents
Example 19. Mass of protective layer

Definitions & inputs. ρ1000kg/m³,A2m²,t.001m.

  1. Choose the governing model and isolate the requested quantity.

    m=ρAtm=\rho At

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    1000(2)(.001)1000(2)(.001)

  3. Evaluate the expression; the result uses the units shown.

    Result=2 kg\mathrm{Result}=2\ \mathrm{kg}

Interpretation. Mass budget does not establish impact performance.

↑ Return to definitions and contents
Example 20. Zero-failure test bound

Definitions & inputs. n30independent identical tests,95percent one-sided.

  1. Choose the governing model and isolate the requested quantity.

    pU=1−.051/np_U=1-.05^{1/n}

  2. Insert the stated inputs in consistent units or the explicitly defined normalized units.

    1−.051/301-.05^{1/30}

  3. Evaluate the expression; the result uses the units shown.

    Result=0.09503385 \mathrm{Result}=0.09503385\ {}

Interpretation. Zero failures do not prove zero failure probability;tests represent only their conditions.

↑ Return to definitions and contents

Symbols and units

Each derivation and problem defines its own symbols and inputs. Symbols may be reused with different meanings in other subjects. Keep units consistent, retain sufficient precision during calculation, and apply the stated validity limits.