Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Cm Math Balances

1 · Prerequisites

2 · Summary

The finite-balance and continuation layer of the non-relativistic classical mechanics collection, all of it mathematics: the algebraic identities behind centre-of-mass and kinetic-energy splittings, momentum, angular-momentum and work-energy balances for finitely many labelled curves, the compact-continuation criterion for locally Lipschitz ordinary differential equations, and the distributional and algebraic machinery of an isolated impulse. Every item on this page carries domain: mathematics and depends only on mathematical suppliers; nothing here asserts that the adopted balances describe nature, and the physical reading of masses, forces, impulses and energies remains with the consuming classical-mechanics pages.

The differentiation layer is Differentiation of sums, scalar multiples, inner products, cross products and squared norms of differentiable curves in R3, which proves the sum, scalar-multiple, inner-product, cross-product and squared-norm rules for differentiable curves in R3. On it rest the balance identities: Centre of mass, the identity ∑imiρi=0, and the splitting of the kinetic sum ∑imi∥xi′∥22 into centre-of-mass and relative parts (with the two-body reduced-mass case) introduces the centre of mass X, proves ∑imiρi=0 and the splitting ∑imi∥xi′∥22=M∥X′∥22+∑imi∥ρi′∥22, with the two-body reduced-mass identity ∑imi∥ρi′∥22=μ∥r′∥22; Total momentum and total angular momentum balances for a finite system with antisymmetric pair forces; pairwise central forces make the internal torque vanish derives P′=∑ifi and J′=∑ixi×fi+∑i<j(xi−xj)×fij from antisymmetric pair forces and shows that pairwise central forces kill the internal torque term; and Work-energy balance ddt(K+U(t,x(t)))=∂tU(t,x(t))+∑i⟨di(t),xi′(t)⟩ for forces mixi′′=−∇xiU+di proves ddt(K+U(t,x(t)))=∂tU(t,x(t))+∑i⟨di,xi′⟩ for forces mixi′′=−∇xiU+di, with the autonomous conserved case. The potential has a C1 extension to an open time neighborhood; relative derivatives retain included endpoints of closed and half-open time intervals. The variable-mass material accounting is Variable-mass momentum balance: from m′=rin−rout and (mv)′=F+rinwin−routwout to mv′=F+rin(win−v)−rout(wout−v), whose two adopted balance identities imply mv′=F+rin(win−v)−rout(wout−v) and exhibit why (mv)′=F alone is not the open-system law.

The continuation layer is the pair Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it and A non-extendable solution leaves every compact subset of the domain near its finite endpoint; an autonomous solution trapped in a compact subset of its state domain cannot break down there: a solution of x′=F(t,x) whose graph stays in a compact subset of the open domain near a finite endpoint b extends past b (and the autonomous trapped-state case is clause 2), while a non-extendable solution must leave every compact subset of the domain as the endpoint is approached. The criterion is about compact trapping, not about pointwise finiteness or boundedness of F.

The impulse layer begins with The distributional derivative of a piecewise C1 function with a jump: (up)′=upreg′+[p] δt∗, which computes the distributional derivative of a piecewise C1 function with a jump, (up)′=upreg′+[p]δt∗, and declares the Axiom of Countable Choice used by the Lebesgue-integration suppliers it cites. On it rests Isolated reciprocal impulse with full-vector restitution: momentum jumps ±J, conservation of total momentum, the restitution solution and the kinetic-energy change −1−e22μ∥v2−−v1−∥22: under the isolated reciprocal impulse model, the momentum jumps are ±J, total momentum is preserved, the restitution law determines the post-impulse velocities uniquely, and the kinetic-energy change is −1−e22μ∥v2−−v1−∥22. The adopted impulse model is stated as a mathematical hypothesis about pairings against smooth test functions; no contact or restitution law is derived from smooth forces.

The companion page cm-math-balances-examples records explicit arithmetic checks of the impulse formulas, a rocket computation, and worked continuation boundary cases.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

Differentiation of sums, scalar multiples, inner products, cross products and squared norms of differentiable curves in R3

Statement

Let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let t0∈I, let x,y:I→R3 be differentiable at t0 (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral) and let λ:I→R be differentiable at t0 (The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set). Write x′(t0),y′(t0)∈R3 for the derivatives, ⟨⋅,⋅⟩ for the Euclidean inner product and ∥⋅∥2 for the Euclidean norm (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), and × for the cross product (The cross product in R3). Then:

  1. x+y and λx are differentiable at t0, with (x+y)′(t0)=x′(t0)+y′(t0) and (λx)′(t0)=λ′(t0)x(t0)+λ(t0)x′(t0);
  2. t↦⟨x(t),y(t)⟩ is differentiable at t0, with derivative ⟨x′(t0),y(t0)⟩+⟨x(t0),y′(t0)⟩;
  3. t↦x(t)×y(t) is differentiable at t0, with derivative x′(t0)×y(t0)+x(t0)×y′(t0);
  4. t↦∥x(t)∥22 is differentiable at t0, with derivative 2⟨x′(t0),x(t0)⟩.

Facts & Assumptions

Given: An interval I⊆R with at least two points, a point t0∈I, curves x,y:I→R3 and a function λ:I→R, with x,y,λ differentiable at t0.

[L1]

A function f:I→Rm is differentiable at t0 if and only if every component fk is differentiable at t0 as a real function on I, and then f′(t0)k=fk′(t0); moreover sums and scalar multiples of differentiable curves are differentiable with (f+g)′(t0)=f′(t0)+g′(t0) and (μf)′(t0)=μf′(t0) for μ∈R (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[L3]

The Euclidean inner product and norm are ⟨u,w⟩=∑k<3ukwk and ∥u∥2=⟨u,u⟩; the inner product is bilinear and symmetric (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn, Laws of finite sums and finite products).

[L4]

For all u,v∈R3 the cross product is bilinear-alternating: (u,v)↦u×v is bilinear, u×u=0, and u×v=−(v×u); in coordinates (u×v)0=u1v2−u2v1, (u×v)1=u2v0−u0v2, (u×v)2=u0v1−u1v0 (The cross product in R3, The cross product is bilinear, alternating, and orthogonal to both factors).

Proof

technique · direct
1.1givenL1L2

By [L1] the components xk,yk:I→R are differentiable at t0 with xk′(t0)=x′(t0)k and yk′(t0)=y′(t0)k, and by [L2] each of the functions xk+yk, λxk, xkyk is differentiable at t0 with the sum, scalar-multiple and product rules; this proves clause 1 componentwise, since (x+y)k=xk+yk and (λx)k=λxk.

2.1step 1.1L2L3

The scalar function t↦⟨x(t),y(t)⟩=∑k<3xk(t)yk(t) is a finite sum of products of functions differentiable at t0, hence by [L2] differentiable at t0 with derivative ∑k<3(xk′(t0)yk(t0)+xk(t0)yk′(t0)); by [L3] this is ⟨x′(t0),y(t0)⟩+⟨x(t0),y′(t0)⟩, which is clause 2.

2.2step 1.1L1L2L4

Each component of t↦x(t)×y(t) is, by [L4], one of the differences x1y2−x2y1, x2y0−x0y2, x0y1−x1y0 of products of functions differentiable at t0, so [L2] makes it differentiable at t0 with derivative (x1′y2+x1y2′)−(x2′y1+x2y1′), and by [L4] and bilinearity this is the corresponding component of x′(t0)×y(t0)+x(t0)×y′(t0); reassembling the three components with [L1] gives clause 3.

3.1step 2.1L3algebra∎

Applying clause 2 at y:=x gives differentiability of t↦⟨x(t),x(t)⟩=∥x(t)∥22 with derivative ⟨x′(t0),x(t0)⟩+⟨x(t0),x′(t0)⟩=2⟨x′(t0),x(t0)⟩, the two terms being equal by symmetry of the inner product [L3]; this is clause 4.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Centre of mass, the identity ∑imiρi=0, and the splitting of the kinetic sum ∑imi∥xi′∥22 into centre-of-mass and relative parts (with the two-body reduced-mass case)

Statement

Let N∈N, N≥1, let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let m1,…,mN>0 and let xi:I→R3 be differentiable at every point of I (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral). All sums over i range over 1,…,N: for a real family a1,…,aN, define ∑iai:=∑k<Nbk, where b:N→R has bk=ak+1 for k<N and bk=0 for k≥N, using the finite sums of Laws of finite sums and finite products; vector sums are defined coordinatewise, and sums of curves pointwise. Put M:=∑imi>0, X:=M−1∑imixi and ρi:=xi−X. The symbols ⟨⋅,⋅⟩ and ∥⋅∥2 denote the Euclidean inner product and norm on R3 (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn). Then:

  1. X and every ρi are differentiable on I, with X′=M−1∑imixi′ and ρi′=xi′−X′, and the identities ∑imiρi=0 and ∑imiρi′=0 hold on I;
  2. for every t∈I, ∑imi∥xi′(t)∥22=M∥X′(t)∥22+∑imi∥ρi′(t)∥22;
  3. if N=2, then with r:=x1−x2 and μ:=m1m2/M one has ∑imi∥ρi′(t)∥22=μ∥r′(t)∥22 and hence ∑imi∥xi′(t)∥22=M∥X′(t)∥22+μ∥r′(t)∥22 for every t∈I.

Facts & Assumptions

Given: A natural N≥1, an interval I⊆R with at least two points, positive reals m1,…,mN and curves xi:I→R3 differentiable at every point of I.

[L1]

If f,g:I→Rm are differentiable at t and μ∈R, then f+g and μf are differentiable at t with (f+g)′(t)=f′(t)+g′(t) and (μf)′(t)=μf′(t); a vector-valued function is differentiable exactly when all its components are (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[L2]

The Euclidean inner product is bilinear and symmetric, ⟨u,w⟩=∑k<3ukwk, and ∥u∥22=⟨u,u⟩ (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

[L3]

Under the statement's shifted indexing convention, additivity and scaling in Laws of finite sums and finite products clauses 1–2 give ∑i(αai+βbi)=α∑iai+β∑ibi for real families and, by applying those scalar laws to each coordinate, vector families. Clause 4 also gives ∑iai≥aj for every 1≤j≤N when all ai≥0.

Proof

technique · direct
1.1givenL1L3

Since N≥1 and all masses are positive, [L3] gives M≥m1>0, so M−1 is defined. Each xi is differentiable on I, so repeated application of [L1] makes the finite sum ∑imixi and its scalar multiple X=M−1∑imixi differentiable on I with X′=M−1∑imixi′, and then every ρi=xi−X is differentiable on I with ρi′=xi′−X′.

2.1step 1.1L1L3algebra

Using the definition of X and [L3], ∑imiρi=∑imixi−MX=∑imixi−∑imixi=0 on I; differentiating this identity through the sum formula of step 1.1 gives ∑imiρi′=0 on I.

3.1step 2.1L2L3algebra

For every t∈I, substituting xi′=X′+ρi′ and expanding with bilinearity and symmetry [L2], ∑imi∥xi′∥22=∑imi(∥X′∥22+2⟨X′,ρi′⟩+∥ρi′∥22)=M∥X′∥22+2⟨X′,∑imiρi′⟩+∑imi∥ρi′∥22, and the middle term vanishes by step 2.1; this is clause 2.

4.1step 3.1L1L2L3algebra∎

Let N=2. From the definitions, X=m1x1+m2x2M and r=x1−x2, so ρ1=x1−X=Mx1−m1x1−m2x2M=m2Mr and ρ2=x2−X=Mx2−m1x1−m2x2M=−m1Mr; differentiating by [L1] and expanding squared norms by [L2] gives ∑imi∥ρi′∥22=(m1m22+m2m12M2)∥r′∥22=m1m2(m1+m2)M2∥r′∥22=μ∥r′∥22, because M=m1+m2. Substitution into clause 2 gives clause 3.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Total momentum and total angular momentum balances for a finite system with antisymmetric pair forces; pairwise central forces make the internal torque vanish

Statement

Let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let N≥1, let m1,…,mN>0, and let xi:I→R3 be twice differentiable, i.e. differentiable with derivative xi′ differentiable (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral). Write pi:=mixi′ and suppose there are maps fi:I→R3 and, for every ordered pair i≠j in {1,…,N}, maps fij:I→R3 such that

fij(t)=−fji(t),pi′(t)=fi(t)+∑j≠ifij(t)(t∈I).

Write P:=∑ipi, J:=∑ixi×pi (cross product in R3, The cross product in R3). Then:

  1. P is differentiable on I with P′(t)=∑ifi(t) for all t∈I;
  2. J is differentiable on I with J′(t)=∑ixi(t)×fi(t)+∑i<j(xi(t)−xj(t))×fij(t)(t∈I);
  3. if in addition, for every i<j, fij is parallel to xi−xj — that is, there is a real function λij on I with fij=λij(xi−xj) — then J′(t)=∑ixi(t)×fi(t) for all t∈I.

Facts & Assumptions

Given: An interval I⊆R with at least two points, N≥1, positive reals m1,…,mN, twice differentiable curves xi:I→R3, velocities vi:=xi′, momenta pi=mivi, and maps fi,fij with fij=−fji and pi′=fi+∑j≠ifij on I.

[L1]

Sums and scalar multiples of curves differentiable at t are differentiable at t with the expected derivatives (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

[L2]

Finite sums over finite sets may be reindexed along bijections and are additive over disjoint unions (Finite commutative-monoid sums are invariant under bijective reindexing, split over disjoint unions, and satisfy the finite Fubini rule clauses 1 and 2); consequently, for vectors aij indexed by the ordered pairs (i,j) with i≠j, the set of such pairs is the disjoint union over the unordered pairs i<j of the two-element sets {(i,j),(j,i)}, so ∑i∑j≠iaij=∑i<j(aij+aji). Scalar factors distribute over finite sums (Laws of finite sums and finite products).

[L3]

For differentiable curves u,v:I→R3 and t∈I, the cross product t↦u(t)×v(t) is differentiable at t with derivative u′(t)×v(t)+u(t)×v′(t) (Differentiation of sums, scalar multiples, inner products, cross products and squared norms of differentiable curves in R3).

[L4]

The cross product is bilinear and alternating, so u×u=0, u×v=−(v×u) and u×(v+w)=u×v+u×w for u,v,w∈R3 (The cross product in R3, The cross product is bilinear, alternating, and orthogonal to both factors).

Proof

technique · direct
1.1givenL1L2algebra

Since P is the finite sum of the differentiable curves pi, [L1] gives P′=∑ipi′; substituting the balance equation and using fij=−fji with the reindexing rule [L2], ∑i∑j≠ifij=∑i<j(fij+fji)=0, so P′=∑ifi; this is clause 1.

2.1step 1.1givenL1L3

Each term xi×pi is differentiable on I with derivative xi′×pi+xi×pi′ by [L3], and [L1] makes their finite sum J differentiable with J′=∑ixi′×pi+∑ixi×pi′=∑ixi′×pi+∑ixi×fi+∑i∑j≠ixi×fij.

3.1step 2.1L4algebra

The first sum in step 2.1 vanishes termwise: xi′×pi=mixi′×xi′=0 by alternation [L4].

4.1step 2.1step 3.1L2L4algebra

The double sum in step 2.1 equals ∑i<j(xi×fij+xj×fji)=∑i<j(xi−xj)×fij: the first equality reindexes the ordered pairs by the unordered ones [L2], and the second uses fji=−fij, bilinearity and u×(−w)=−(u×w) [L4]. Combining with steps 2.1 and 3.1 gives clause 2.

5.1step 4.1givenL4algebra∎

Under the parallelism hypothesis, (xi−xj)×fij=λij(xi−xj)×(xi−xj)=0 for every i<j by alternation [L4]; hence the double sum in clause 2 vanishes and J′=∑ixi×fi, which is clause 3.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Work-energy balance ddt(K+U(t,x(t)))=∂tU(t,x(t))+∑i⟨di(t),xi′(t)⟩ for forces mixi′′=−∇xiU+di

Statement

Let N≥1, let I⊆R be an interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let Q⊆R3N be nonempty and open, and let U:I×Q→R be the restriction of a function U^:J×Q→R of class C1, where J⊆R is an open interval containing I (Ck Euclidean maps and diffeomorphisms). Write points of Q as q=(q1,…,qN) with qi∈R3; for t∈I, q∈Q and each i let ∇xiU(t,q)∈R3 be the gradient of u↦U(t,q1,…,qi−1,u,qi+1,…,qN) at u=qi (The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case), and let ∂tU and the spatial derivatives on I×Q denote the derivatives of U^ restricted there. At an included endpoint of I, derivatives of curves are limits relative to I (one-sided there); this preserves closed and half-open time intervals. Let m1,…,mN>0, let x:I→Q be twice differentiable (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral), put K(t):=12∑imi∥xi′(t)∥22 (Euclidean norm, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), let di:I→R3 be maps, and suppose

mixi′′(t)=−∇xiU(t,x(t))+di(t)(t∈I, 1≤i≤N).

Then:

  1. K is differentiable on I with K′(t)=∑i⟨mixi′′(t),xi′(t)⟩;
  2. t↦U(t,x(t)) is differentiable on I with ddtU(t,x(t))=∂tU(t,x(t))+∑i⟨∇xiU(t,x(t)),xi′(t)⟩;
  3. consequently ddt(K(t)+U(t,x(t)))=∂tU(t,x(t))+∑i⟨di(t),xi′(t)⟩(t∈I); in particular, if di≡0 for all i then ddt(K+U(t,x(t)))=∂tU(t,x(t)), and if in addition ∂tU≡0 then t↦K(t)+U(t,x(t)) is constant on I.

Facts & Assumptions

Given: The data of the statement: N≥1, an interval I with at least two points, open Q⊆R3N, a function U:I×Q→R with the stated C1 extension U^ to J×Q, positive masses mi, a twice differentiable x:I→Q, and maps di:I→R3 with mixi′′=−∇xiU(t,x(t))+di on I.

[L3]

The total derivative computes directional derivatives: Dvf(a)=Df(a)v (A total derivative computes every directional derivative, and its matrix is the Jacobian).

[L4]

The Euclidean inner product is bilinear and symmetric, ⟨u,w⟩=∑k<3ukwk (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), and finite sums are linear (Laws of finite sums and finite products).

Proof

technique · direct
1.1givenL1L4algebra

Each xi′ is differentiable on I because x is twice differentiable, so [L1] gives differentiability of t↦∥xi′(t)∥22 with derivative 2⟨xi′′(t),xi′(t)⟩; by the sum and scalar rules of [L1] and linearity of the finite sum [L4], K is differentiable with K′=12∑imi⋅2⟨xi′′,xi′⟩=∑i⟨mixi′′,xi′⟩, which is clause 1.

1.2givenL2L3L4algebra

The curve γ:I→R×Q, γ(t):=(t,x(t)), is differentiable (its components are the identity and the components of x; The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral), with γ′(t)=(1,x′(t)) under the identification R×R3N=R1+3N. Since U^ is C1 on the open neighborhood J×Q of γ(I), [L2] gives the usual chain rule at interior points. At an included endpoint, γ(t+h)=γ(t)+hγ′(t)+o(h) as h→0 relative to I, and total differentiability gives U^(γ(t)+v)=U^(γ(t))+DU^(γ(t))v+o(∣v∣). Substituting the first expansion into the second proves the same relative derivative there. Thus (U∘γ)′(t)=DU^(γ(t))γ′(t), and [L3] evaluates the right side as the matrix product of the Jacobian of U at γ(t) with (1,x′(t)), that is ddtU(t,x(t))=∂tU(t,x(t))+∑i∑k<3∂i,kU(t,x(t)) xi,k′(t)=∂tU(t,x(t))+∑i⟨∇xiU(t,x(t)),xi′(t)⟩, the last step regrouping the coordinates of Q into the N blocks R3 and using the coordinate formula for the inner product [L4]; this is clause 2.

2.1step 1.1step 1.2givenL4L5algebra∎

Adding the identities of clauses 1 and 2 and substituting mixi′′=−∇xiU(t,x(t))+di, the two occurrences of ∑i⟨∇xiU(t,x(t)),xi′(t)⟩ cancel, leaving clause 3. If di≡0 the derivative is ∂tU(t,x(t)), and if also ∂tU≡0 then the derivative of K+U(t,x(t)) vanishes identically on the interval I, so that function is constant by [L5].

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Variable-mass momentum balance: from m′=rin−rout and (mv)′=F+rinwin−routwout to mv′=F+rin(win−v)−rout(wout−v)

Statement

Let I⊆R be an interval with nonempty interior (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let m:I→(0,∞) and v:I→R3 be differentiable (The derivative f′(c)=lim⁡x→cf(x)−f(c)x−c of f:A→R at a point c∈A that is a limit point of A, and differentiability on a set, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral), and let F:I→R3, rin,rout:I→[0,∞) and win,wout:I→R3 be maps. Suppose the two material balance identities

m′(t)=rin(t)−rout(t),(mv)′(t)=F(t)+rin(t)win(t)−rout(t)wout(t)

hold for every t∈I, where mv denotes the pointwise scalar multiple. Then

m(t)v′(t)=F(t)+rin(t)(win(t)−v(t))−rout(t)(wout(t)−v(t))(t∈I).

In particular, if F≡0, rin≡0 and rout≡r for a constant r≥0, then mv′=−r(wout−v): the acceleration is governed by the momentum carried away per unit time relative to the body, not by (mv)′=0.

The two balance identities are the adopted content of the material momentum accounting — the mathematical rendering of "gains rin dt of matter with velocity win, loses rout dt with velocity wout, and feels the external force F". They are not derived here from any closed-system law; the deduction below is their product-rule consequence.

Facts & Assumptions

Given: An interval I with nonempty interior, a differentiable positive real function m and vector curve v, maps F,win,wout into R3, nonnegative functions rin,rout, and the two material balance identities on I.

Proof

technique · direct
1.1givenL1

By [L1] the product rule applies coordinatewise to mv: (mv)′=m′v+mv′, since (mv)k=mvk and (mvk)′=m′vk+mvk′ for each coordinate k.

2.1step 1.1givenL1algebra∎

Substituting the first balance identity into the second and using step 1.1, mv′=(mv)′−m′v=F+rinwin−routwout−(rin−rout)v=F+rin(win−v)−rout(wout−v), the last step grouping the rin-terms and the rout-terms and using bilinearity of scalar multiplication. This is the displayed conclusion; the special case F≡0, rin≡0, rout≡r is its immediate substitution, in that case the assumed momentum balance is (mv)′=−rwout, so (mv)′=0 at a time t exactly when rwout(t)=0.

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it

Statement

Let n≥1, let D⊆R×Rn be open, and let F:D→Rn be continuous and locally Lipschitz in the state variable, locally uniformly in time, in the sense of Local Lipschitz continuity in the state variable, locally uniform in time and parameters. Let t0<b be reals, let x0∈Rn with (t0,x0)∈D, and let x:[t0,b)→Rn be a solution of the initial value problem

x′(t)=F(t,x(t)),x(t0)=x0

on [t0,b) (First-order systems, initial value problems, and solutions on intervals). Suppose there is a compact set K⊆D and a real δ>0 such that

(t,x(t))∈Kfor every t∈(b−δ,b)∩[t0,b).

Then:

  1. there is ε>0 and a solution x~:[t0,b+ε)→Rn of the same initial value problem with x~(t)=x(t) for all t∈[t0,b); in particular x~(b)=lim⁡t↑bx(t) exists and lies in the state projection of K;
  2. if moreover D=R×U for an open U⊆Rn and F is autonomous, F(t,z)=G(z) with G:U→Rn continuous and locally Lipschitz in the state variable (Local Lipschitz continuity in the state variable, locally uniform in time and parameters), and if there is a compact K0⊆U with x(t)∈K0 for every t∈(b−δ,b)∩[t0,b), then the extension conclusion of clause 1 holds as well.

Facts & Assumptions

Given: The data of the statement: open D⊆R×Rn, a continuous F that is locally Lipschitz in the state variable in the sense of [L2], a solution x of the initial value problem on [t0,b), and a compact K⊆D with (t,x(t))∈K for all t∈(b−δ,b)∩[t0,b).

[L1]

At every point (t∗,y∗)∈D, the initial value problem has a unique local solution on some interval around t∗; any two solutions through the same initial data agree on every common subinterval containing the initial time (Picard-Lindelöf local existence and uniqueness for first-order systems, First-order systems, initial value problems, and solutions on intervals).

[L2]

In the cylinder terminology of Local Lipschitz continuity in the state variable, locally uniform in time and parameters, this item uses C=[a,c]×B‾(z,r)⊆D, where a,c∈R, a≤c, z∈Rn, r>0, and B‾(z,r)={w∈Rn:∥w−z∥2≤r} (Open ball, closed ball and sphere in a metric space). F is locally Lipschitz in the state variable when every such cylinder carries a finite L with ∥F(t,u)−F(t,w)∥2≤L∥u−w∥2 for (t,u),(t,w)∈C (Local Lipschitz continuity in the state variable, locally uniform in time and parameters).

[L4]

For s<t and an integrable curve g with values in Rn one has ∥∫stg∥2≤∫st∥g∥2, and for a real integrable u≤M on [s,t] one has ∫stu≤M(t−s) (For a≤b and f:[a,b]→Rm integrable when a<b, ∥∫abf∥2≤∫ab∥f∥2; for a<b, ∥f∥2 is integrable, If f≤g on [a,b] and both are integrable then ∫abf≤∫abg; and m(b−a)≤∫abf≤M(b−a)).

[L5]

The Euclidean norm is a norm on Rn, so it satisfies the triangle inequality ∥u+w∥2≤∥u∥2+∥w∥2 (Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation) and the reverse triangle inequality ∣∥u∥2−∥w∥2∣≤∥u−w∥2 (The reverse triangle inequality in a normed space); moreover (Rn,d2) is complete: every Cauchy sequence in Rn converges (R and Rn for n≥1 with the Euclidean metric are complete, componentwise from the Cauchy criterion in R, Complete metric space: every Cauchy sequence converges in the space, Vector-valued functions f:A→Rm, their limits and continuity, with the dictionary to the metric notions).

[L8]

A continuous real-valued function on a closed interval [s,t], s<t, is integrable; hence a continuous curve into Rn on [s,t] is integrable componentwise (A continuous function on [a,b] is Riemann integrable, by Heine-Cantor and Riemann's criterion, The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral).

Proof

technique · direct
1.1givenL5L6L7

If δ>b−t0 replace δ by b−t0; then the interval (b−δ,b)∩[t0,b) is nonempty. Choose any t in it; the hypothesis gives (t,x(t))∈K, so K is nonempty. The norm function z↦∥F(z)∥2 is continuous on the compact metric subspace K: given z0∈K and ε>0, continuity of F at z0 [L7] gives δ0>0 with ∥F(z)−F(z0)∥2<ε whenever d2(z,z0)<δ0, and the reverse triangle inequality [L5] then gives ∣∥F(z)∥2−∥F(z0)∥2∣≤∥F(z)−F(z0)∥2<ε. Hence [L6] gives M:=max⁡{∥F(z)∥2:z∈K}<+∞.

2.1step 1.1givenL3L4L7L8

For s,t with b−δ<s<t<b: since x solves the initial value problem, [L3] gives x(t)=x0+∫t0tF(u,x(u)) du and x(s)=x0+∫t0sF(u,x(u)) du, hence x(t)−x(s)=∫stF(u,x(u)) du; the integrand is continuous in u on [s,t] because u↦F(u,x(u)) is a composition of continuous maps (the components of x are continuous by differentiability, [L7]) and (u,x(u))∈K on this interval, so its components are integrable on [s,t] by [L8]. By [L4], ∥x(t)−x(s)∥2≤∫st∥F(u,x(u))∥2 du≤M(t−s).

3.1step 2.1L5L7algebra

The limit xb:=lim⁡t↑bx(t) exists. Indeed, choose k0≥1 with 1/k0<min⁡{δ,b−t0} [L7] and put yk:=x(b−1/k) for k≥k0; by step 2.1, ∥yk−yl∥2≤M ∣1/k−1/l∣ for all k,l≥k0, so (yk) is a Cauchy sequence in Rn and converges, by [L5], to some xb∈Rn. Given η>0, choose k≥k0 with 1/k<η/(4(M+1)) and ∥yk−xb∥2<η/2; then for every t∈(max⁡{b−δ,b−η/(4(M+1))},b), step 2.1 gives ∥x(t)−xb∥2≤∥x(t)−yk∥2+∥yk−xb∥2≤M(∣t−b∣+1/k)+η/2<M⋅η/(2(M+1))+η/2<η. Hence lim⁡t↑bx(t)=xb.

4.1step 3.1L6

The point (b,xb) lies in K, hence in D: for t∈(b−δ,b)∩[t0,b) the points (t,x(t)) lie in K and converge to (b,xb) as t↑b (in the product metric, the first coordinate trivially and the second by step 3.1); since K is closed [L6], its complement is open, so a point outside K would have a ball around it disjoint from K, which is impossible for (b,xb) because every such ball contains some (t,x(t)).

5.1step 2.1step 3.1step 4.1givenL3L4L6L8algebra

Extend x to the continuous curve xˉ:[t0,b]→Rn by xˉ(t):=x(t) for t<b and xˉ(b):=xb; it is continuous by step 3.1. Let 0<h0<δ. For s∈(b−h0,b) and t∈(s,b), step 2.1 and additivity of the integral [L3] give (x(t)−x(s))−∫sbF(u,xˉ(u)) du=−∫tbF(u,xˉ(u)) du; since (u,xˉ(u))∈K for every u∈[s,b] (for u<b by hypothesis, for u=b by step 4.1), one has ∥F(u,xˉ(u))∥2≤M there, and the norm inequality [L4] gives ∥(x(t)−x(s))−∫sbF(u,xˉ(u)) du∥2≤∫tb∥F(u,xˉ(u))∥2 du≤M(b−t)→0 as t↑b. With x(t)→xb (step 3.1) this yields xb−x(s)=∫sbF(u,xˉ(u)) du. Hence the left difference quotient satisfies ∥xˉ(b)−xˉ(s)b−s−F(b,xb)∥2≤sup⁡s≤u≤b∥F(u,xˉ(u))−F(b,xb)∥2⟶0(s↑b) because for every η>0 continuity of xˉ at b and of F at (b,xb) makes the displayed integrand smaller than η for all u∈[s,b] once s is sufficiently close to b; [L6] ensures each supremum is finite. So xˉ is differentiable at every point of [b−h0,b] as a function on that interval, with xˉ′(t)=F(t,xˉ(t)) for t<b and left derivative F(b,xb) at t=b; it is therefore, for every such h0, a solution of the initial value problem z′=F(t,z), z(b)=xb on [b−h0,b] in the sense of [L1] and [L3].

5.2step 4.1L1L2

Apply the local existence conclusion of Picard-Lindelöf [L1] at (b,xb)∈D, whose hypotheses on F are the assumed continuity and local Lipschitz continuity in the state variable [L2]. It gives a solution through (b,xb) on some interval around b; choose h>0 so [b−h,b+h] is contained in that interval and denote the restricted solution by y. Then y is a solution on [b−h,b+h] through (b,xb).

6.1step 5.2step 5.1L1algebra

Let h′:=min⁡{h,δ/2}>0, where h is the Picard-Lindelöf half-length of step 5.2; by step 5.1 the curve xˉ is a solution on [b−h′,b] through (b,xb), and the restriction of y is a solution on [b−h,b]⊇[b−h′,b] through the same point; by the uniqueness clause of [L1] they agree on [b−h′,b]. Define x~:[t0,b+h)→Rn by x~(t):=xˉ(t) for t≤b and x~(t):=y(t) for b<t<b+h. Then x~ is continuous, and it is differentiable with x~′(t)=F(t,x~(t)) at every t∈[t0,b+h): for t<b and for t>b this is the respective solution property, while at t=b the left difference quotient tends to F(b,xb) by step 5.1 and the right difference quotient is the two-sided derivative y′(b)=F(b,xb) of the solution y; the two agree, so the derivative of x~ at b exists and equals F(b,xb)=F(b,x~(b)). With ε:=h this is an extension as asserted in clause 1.

7.1step 1.1step 2.1step 3.1step 4.1step 5.2step 5.1step 6.1givenL6algebra∎

Under the hypotheses of clause 2, the argument of steps 1.1-6.1 applies with K replaced by the given compact K0⊆U: the bound of step 1.1 is M:=max⁡{∥G(z)∥2:z∈K0}<+∞ by [L6], the Lipschitz estimate of step 2.1 becomes ∥x(t)−x(s)∥2≤M(t−s) because F(u,x(u))=G(x(u)) and x(u)∈K0 for u∈(b−δ,b)∩[t0,b), the limit point of step 4.1 lies in K0 since K0 is compact hence closed [L6], and (b,xb)∈R×U=D because xb∈K0⊆U; the local solution of step 5.2 and the gluing of step 6.1 are then carried out for the same field F(t,z)=G(z) at (b,xb). This proves clause 2 with the same ε=h.

CorollaryStatement: AI-adaptedProof: AI-adaptedOpen item page →

A non-extendable solution leaves every compact subset of the domain near its finite endpoint; an autonomous solution trapped in a compact subset of its state domain cannot break down there

Statement

Let n≥1, let D⊆R×Rn be open, and let F:D→Rn be continuous and locally Lipschitz in the state variable in the sense of Local Lipschitz continuity in the state variable, locally uniform in time and parameters. Let t0<b be reals, let x0∈Rn with (t0,x0)∈D, and let x:[t0,b)→Rn be a solution of the initial value problem x′=F(t,x), x(t0)=x0 (First-order systems, initial value problems, and solutions on intervals). Suppose x admits no extension past b: there do not exist ε>0 and a solution x~:[t0,b+ε)→Rn of the same initial value problem with x~(t)=x(t) for every t∈[t0,b). Then:

  1. for every compact K⊆D and every η>0 there is t∈(b−η,b)∩[t0,b) with (t,x(t))∉K; that is, the graph of x leaves every compact subset of D that could contain it arbitrarily close to b;
  2. if D=R×U for an open U⊆Rn, F is autonomous, F(t,z)=G(z) with G:U→Rn continuous and locally Lipschitz in the state variable, and there are a compact K0⊆U and δ>0 with x(t)∈K0 for all t∈(b−δ,b)∩[t0,b), then the no-extension hypothesis is false: x does extend past b.

Facts & Assumptions

Given: The data of the statement and the no-extension hypothesis on x.

[L1]

Compact continuation: under the hypotheses of the statement without the no-extension hypothesis, if there are a compact K⊆D and δ>0 with (t,x(t))∈K for all t∈(b−δ,b)∩[t0,b), then x extends past b; and the same holds for an autonomous field F(t,z)=G(z) on D=R×U when x(t)∈K0 for a compact K0⊆U and all t∈(b−δ,b)∩[t0,b) (Compact continuation for locally Lipschitz ODEs: a solution whose graph stays in a compact subset of the domain near a finite endpoint extends past it).

Proof

technique · contrapositive
1.1givenL1

Clause 1 is the contrapositive of the compact continuation statement in [L1]: if for some compact K⊆D and some η>0 every t∈(b−η,b)∩[t0,b) had (t,x(t))∈K, then [L1] would produce an extension of x past b, contradicting the no-extension hypothesis; hence for every compact K⊆D and every η>0 such a t exists.

2.1step 1.1givenL1∎

Clause 2 likewise contradicts the no-extension hypothesis: the autonomy and trapping hypotheses are exactly those of the second part of [L1], which would give an extension of x past b; therefore no such trapping compact set K0 can exist under the no-extension hypothesis, and a fortiori every autonomous solution that is trapped in such a K0 near b does extend past b.

LemmaStatement: AI-adaptedProof: AI-adaptedOpen item page →

The distributional derivative of a piecewise C1 function with a jump: (up)′=upreg′+[p] δt∗

Statement

Assume the Axiom of Countable Choice (ACω, The Axiom of Countable Choice (ACω)), used only through the Borel-measurability, Lebesgue-integration and distribution suppliers cited below. Let I⊆R be an open interval (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let t∗∈I, and let p:I→R satisfy:

  1. p is continuous on I∖{t∗} and of class C1 on each of the two open intervals I∩(−∞,t∗) and I∩(t∗,∞);
  2. the one-sided limits p(t∗−):=lim⁡t↑t∗p(t), p(t∗+):=lim⁡t↓t∗p(t) and the one-sided limits of the derivative p′(t∗−):=lim⁡t↑t∗p′(t), p′(t∗+):=lim⁡t↓t∗p′(t) all exist in R.

Fix a value p(t∗)∈R arbitrarily (for instance p(t∗):=0) and define preg′:I→R by preg′(t):=p′(t) for t≠t∗ and preg′(t∗):=0. Put [p]:=p(t∗+)−p(t∗−) and let up,upreg′∈D′(I) be the associated regular distributions (Regular distribution from a locally integrable function, A locally integrable function on Rn): up(φ)=∫Ip(t) φ(t) dλ(t),upreg′(φ)=∫Ipreg′(t) φ(t) dλ(t).

Then p and preg′ are locally integrable on I, and

(up)′=upreg′+[p] δt∗

in D′(I) (Distributional derivative, Dirac delta and its derivatives): that is, for every test function φ∈D(I),

−∫Ip(t) φ′(t) dλ(t)=∫Ipreg′(t) φ(t) dλ(t)+[p] φ(t∗).

Both integrals are over the Lebesgue measure λ on R, and the pairing convention of the two sides is complex-linear in φ; it therefore suffices that the identity hold for real-valued test functions.

Facts & Assumptions

Given: ACω, an open interval I, a point t∗∈I, and a real function p satisfying the piecewise-C1 and one-sided-limit hypotheses above.

[L1]

The test space D(I) consists of the compactly supported C∞ functions on I; for φ∈D(I) the pairing of a distribution u is written ⟨u,φ⟩, and (u′)(φ)=−u(φ′); the Dirac distribution satisfies δt∗(φ)=φ(t∗); ∫I and ∫E denote the Lebesgue integral, with ∫Ef dλ=∫fχE dλ; a locally integrable f defines the regular distribution uf(φ)=∫Ifφ dλ (Distribution, Distributional derivative, Dirac delta and its derivatives, Test function space d of an open set, Regular distribution from a locally integrable function, A locally integrable function on Rn, Integral over a measurable subset).

[L2]

Complex Lebesgue integration by parts: for complex C1 functions u,v on a closed interval [a,b], a<b, one has ∫abu′v=u(b)v(b)−u(a)v(a)−∫abuv′ and ∫abu′=u(b)−u(a) (Complex integration by parts on intervals and decaying lines).

[L3]

Functions agreeing almost everywhere have the same integral over every measurable set, and the Lebesgue integral is linear on L1 (Two integrable functions are equal almost everywhere exactly when all of their indefinite integrals agree, The Lebesgue integral is linear on L1(μ)); a bounded real function on a closed bounded interval that is continuous except at finitely many points is Riemann integrable, and a bounded Borel Riemann integrable function on a box is Lebesgue integrable with the same integral (A bounded function on [a,b] that is continuous except at finitely many points is Riemann integrable, Riemann–Lebesgue comparison for distribution test integrands).

[L4]

A real function continuous on [a,b] and differentiable on (a,b) has f(b)−f(a)=f′(c)(b−a) for some c∈(a,b). (The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a))

[L5]

A compact subset of the metric space R is closed, hence Borel, and under ACω every singleton in R is Lebesgue-null (A compact subset of a metric space is closed and bounded, Every at most countable subset of Rn is Lebesgue null; in particular λ1(Q)=0).

[L6]

For nonnegative measurable functions, pointwise order implies the same order of their Lebesgue integrals; with the integral-over-a-set convention, E⊆F gives ∫Eh dλ≤∫Fh dλ (Monotonicity and nonnegative homogeneity of the nonnegative integral, Integral over a measurable subset).

[L7]

A continuous real function on a nonempty compact metric space is bounded and attains its bounds (A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value).

Proof

technique · direct
1.1givenL1L3L5L6L7

Local integrability: put I−:=I∩(−∞,t∗) and I+:=I∩(t∗,∞). Both are open, and p and preg′ are continuous on each. For any open O⊆R, the preimage on either half is open relative to that open half, hence open in R; the preimage at t∗ is either empty or {t∗}. Thus both functions are Borel on I. Let K⊆I be compact. If K=∅, its integral is zero. Otherwise choose c<t∗<d with K∪{t∗}⊂(c,d) and [c,d]⊂I. Extending each one-sided reading of p and preg′ to t∗ by its finite one-sided limit gives continuous functions on [c,t∗] and [t∗,d]; [L7] bounds them there, and the fixed finite values at t∗ preserve boundedness. Hence p and preg′ are bounded on [c,d], and their absolute values are Borel and continuous except possibly at t∗. By [L3] they are Riemann integrable and Lebesgue integrable on [c,d]. The set K is Borel by [L5], and [L6] gives ∫K∣p∣ dλ≤∫[c,d]∣p∣ dλ<∞ and ∫K∣preg′∣ dλ≤∫[c,d]∣preg′∣ dλ<∞. Thus both functions are locally integrable and define the regular distributions in [L1].

1.2givenL1L2L3L5algebra

Let φ∈D(I) be real-valued. Choose a<t∗<b with supp⁡φ∪{t∗}⊂(a,b) and [a,b]⊂I; compactness of the support and openness of I allow this, without requiring a symmetric interval around t∗. Since φ vanishes outside [a,b], [L1] and [L3] give ∫Ipφ′ dλ=∫[a,b]pφ′ dλ=∫[a,t∗]pφ′ dλ+∫[t∗,b]pφ′ dλ, the second equality splitting the integral over the two overlapping half-intervals by linearity and the identity χ[a,b]=χ[a,t∗]+χ[t∗,b]−χ{t∗} (the subtracted term contributes 0 because the singleton {t∗} has measure zero).

2.1step 1.2givenL2L3L4L5algebra

On the left half-interval define pˉ−(t):=p(t) for t<t∗ and pˉ−(t∗):=p(t∗−); on [a,t∗] this extension is continuous. For t<t∗, the mean value theorem [L4] on [t,t∗] gives (pˉ−(t∗)−pˉ−(t))/(t∗−t)=p′(ξt) for some t<ξt<t∗. As t↑t∗ this tends to p′(t∗−), proving the endpoint derivative. Its derivative is continuous up to t∗; extend beyond t∗ by the affine function p(t∗−)+p′(t∗−)(t−t∗) to obtain a C1 function on an open neighborhood of [a,t∗]. On [a,t∗], p and pˉ− agree except at t∗, so [L3] and [L5] give ∫[a,t∗]pφ′ dλ=∫at∗pˉ−φ′ dλ. Lebesgue integration by parts [L2] for the C1 pair (pˉ−,φ) gives ∫at∗pˉ−φ′=pˉ−(t∗)φ(t∗)−pˉ−(a)φ(a)−∫at∗pˉ−′φ=p(t∗−)φ(t∗)−∫[a,t∗]preg′φ dλ, because φ(a)=0 and pˉ−′=preg′ on (a,t∗); the last integrals agree by [L3] and [L5].

2.2step 1.2givenL2L3L4L5algebra

On the right half-interval define pˉ+(t):=p(t) for t>t∗ and pˉ+(t∗):=p(t∗+). The mean value theorem [L4] on [t∗,t] shows that the right derivative of pˉ+ at t∗ is p′(t∗+); extending affinely to the left gives a C1 function on a neighborhood of [t∗,b]. Since p and pˉ+, and preg′ and pˉ+′, agree except at t∗, [L3] and [L5] give equality of their integrals against the tests on [t∗,b]. Integration by parts [L2], together with φ(b)=0, now gives ∫[t∗,b]pφ′ dλ=∫t∗bpˉ+φ′ dλ=pˉ+(b)φ(b)−pˉ+(t∗)φ(t∗)−∫t∗bpˉ+′φ=−p(t∗+)φ(t∗)−∫[t∗,b]preg′φ dλ.

3.1step 2.1step 2.2L1L3algebra∎

Adding the two identities of steps 2.1 and 2.2 and substituting into step 1.2 gives ∫Ipφ′ dλ=−(p(t∗+)−p(t∗−))φ(t∗)−∫Ipreg′φ dλ=−[p] φ(t∗)−∫Ipreg′φ dλ, that is, ⟨(up)′,φ⟩=−∫Ipφ′ dλ=∫Ipreg′φ dλ+[p]φ(t∗)=⟨upreg′+[p]δt∗,φ⟩ for every real-valued φ∈D(I) by [L1]. Since both sides are complex-linear in φ, the identity extends to all φ∈D(I) by applying it to Re⁡φ and Im⁡φ; hence (up)′=upreg′+[p]δt∗ in D′(I).

TheoremStatement: AI-adaptedProof: AI-adaptedOpen item page →

Isolated reciprocal impulse with full-vector restitution: momentum jumps ±J, conservation of total momentum, the restitution solution and the kinetic-energy change −1−e22μ∥v2−−v1−∥22

Statement

Assume the Axiom of Countable Choice (ACω, The Axiom of Countable Choice (ACω)), inherited from the jump lemma below. Let I⊆R be an open interval with at least two points (Intervals of R: the nine order-convex forms, nondegeneracy, and length), let t∗∈I, let m1,m2>0, and let x1,x2:I→R3 be continuous on I and twice continuously differentiable on each of the two open intervals I∩(−∞,t∗) and I∩(t∗,∞) (The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral, Ck Euclidean maps and diffeomorphisms). Assume the one-sided limits vi−:=lim⁡t↑t∗xi′(t),vi+:=lim⁡t↓t∗xi′(t)(i=1,2) exist in R3. Define pi:I→R3 by pi(t):=mixi′(t) for t≠t∗ and fix the representative value pi(t∗):=0 (Vector-valued functions f:A→Rm, their limits and continuity, with the dictionary to the metric notions). Let F1,F2:I→R3 be continuous and let J∈R3. The isolated reciprocal impulse model at t∗ consists of the two requirements:

  1. regular dynamics off the impulse: mixi′′(t)=Fi(t) for every t∈I∖{t∗} and i=1,2;
  2. impulse pairing: for every Φ∈Cc∞(I;R3) — meaning every map whose three components lie in Cc∞(I) (Test function space d of an open set, Ck Euclidean maps and diffeomorphisms) — −∫I⟨p1(t),Φ′(t)⟩ dλ(t)=∫I⟨F1(t),Φ(t)⟩ dλ(t)+⟨J,Φ(t∗)⟩, −∫I⟨p2(t),Φ′(t)⟩ dλ(t)=∫I⟨F2(t),Φ(t)⟩ dλ(t)−⟨J,Φ(t∗)⟩, with integrals against Lebesgue measure λ and the Euclidean inner product ⟨⋅,⋅⟩ (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn), using the fixed representative pi(t∗)=0 above.

Under this model, with M:=m1+m2 and μ:=m1m2/M:

a. m1(v1+−v1−)=J and m2(v2+−v2−)=−J; consequently the total momentum is continuous across the impulse, p1(t∗+)+p2(t∗+)=p1(t∗−)+p2(t∗−).

b. Let e∈[0,1] satisfy the restitution law v2+−v1+=−e(v2−−v1−). Then the post-impulse velocities are uniquely determined: with V:=m1v1−+m2v2−M=m1v1++m2v2+M, v1+=V+m2Me (v2−−v1−),v2+=V−m1Me (v2−−v1−).

c. The kinetic-energy change across the impulse is ΔK:=m1∥v1+∥22+m2∥v2+∥222−m1∥v1−∥22+m2∥v2−∥222=−1−e22 μ ∥v2−−v1−∥22≤0; in particular ΔK=0 for e=1 and ΔK=−12μ∥v2−−v1−∥22 for e=0.

Facts & Assumptions

Given: ACω; an open interval I, a point t∗∈I, masses m1,m2>0; position curves x1,x2 continuous on I and twice differentiable off t∗ with one-sided velocity limits vi±; continuous F1,F2 and a vector J; the isolated reciprocal impulse model of the statement; and, for claims b and c, a number e∈[0,1] satisfying the restitution law.

[L1]

Jump lemma: for a real function p on I that is continuous and C1 off t∗ and whose one-sided limits p(t∗±) and p′(t∗±) exist, the regular distributions of p and of its off-t∗ derivative satisfy (up)′=upreg′+[p]δt∗ with [p]=p(t∗+)−p(t∗−) (The distributional derivative of a piecewise C1 function with a jump: (up)′=upreg′+[p] δt∗, The Axiom of Countable Choice (ACω)).

[L2]

Two-body kinetic identity: for all vectors u,w∈R3 and V:=m1u+m2wM, m1∥u∥22+m2∥w∥22=M∥V∥22+μ∥u−w∥22; this is the two-body case of the kinetic splitting identity, whose proof is the algebraic expansion of the right-hand side (Centre of mass, the identity ∑imiρi=0, and the splitting of the kinetic sum ∑imi∥xi′∥22 into centre-of-mass and relative parts (with the two-body reduced-mass case), The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn). The Euclidean norm is a norm on R3, hence homogeneous: ∥cw∥2=∣c∣ ∥w∥2 for real c (Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation).

[L3]

Test-function cutoffs: for compact K⊆Ω with Ω⊆Rm open there is a smooth compactly supported χ∈Cc∞(Ω) with 0≤χ≤1 and χ=1 on a neighbourhood of K; in particular some χ∈Cc∞(I) has χ(t∗)=1 (Test function cutoffs and euclidean localization, Test function space d of an open set, Ck Euclidean maps and diffeomorphisms). The Euclidean inner product is bilinear and symmetric (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn).

Proof

technique · direct
1.1givenL3algebra

For k∈{0,1,2}, let ek∈R3 be the vector with kth coordinate 1 and other coordinates 0, and let δkl=1 if k=l and 0 otherwise. Fix a real-valued φ∈Cc∞(I); the vector test function Φ:=φ ek has components (φek)l=δklφ and lies in Cc∞(I;R3). Reading the impulse pairing of the model for this Φ gives −∫Ip1,k(t)φ′(t) dλ(t)=∫IF1,k(t)φ(t) dλ(t)+Jkφ(t∗),−∫Ip2,k(t)φ′(t) dλ(t)=∫IF2,k(t)φ(t) dλ(t)−Jkφ(t∗).

2.1step 1.1givenL1

Each component pi,k satisfies the hypotheses of the jump lemma [L1]: it is continuous and C1 off t∗, its one-sided limits exist, and its one-sided derivative limits exist too, because off t∗ one has pi,k′=mixi,k′′=Fi,k by the regular-dynamics requirement of the model and Fi,k is continuous on I. Moreover the off-t∗ derivative of pi,k is Fi,k, so (pi,k)reg′=Fi,k except possibly at t∗, which does not change the regular distribution or its integral. Hence [L1] gives, for every real φ∈Cc∞(I), −∫Ipi,kφ′ dλ=∫IFi,kφ dλ+[pi,k] φ(t∗).

3.1step 2.1givenL3algebra

Comparing the identity of step 2.1 with the model identity of step 1.1 for the same φ, the terms ∫IFi,kφ dλ cancel and leave ([p1,k]−Jk)φ(t∗)=0 and ([p2,k]+Jk)φ(t∗)=0 for every φ∈Cc∞(I); choosing by [L3] a φ with φ(t∗)=1 gives [p1,k]=Jk and [p2,k]=−Jk. As k is arbitrary, m1(v1+−v1−)=J and m2(v2+−v2−)=−J, and adding these two vector equations gives p1(t∗+)+p2(t∗+)=p1(t∗−)+p2(t∗−); this is claim a.

4.1step 3.1givenalgebra

For claim b, let V:=m1v1−+m2v2−M and d:=v2−−v1−. By step 3.1, V=m1v1++m2v2+M as well. The restitution law reads v2+−v1+=−ed. Writing d+:=v2+−v1+ and solving the pair of linear equations (v2+=v1++d+ and m1v1++m2v2+=MV) gives v1+=V−m2Md+=V+m2Med and v2+=V+m1Md+=V−m1Med; conversely these displayed vectors satisfy both equations, so they are the unique solution and claim b holds.

5.1step 4.1givenL2algebra∎

For claim c, apply the kinetic identity [L2] to the pairs (v1−,v2−) and (v1+,v2+), whose centre-of-mass velocities are both V: m1∥v1±∥22+m2∥v2±∥222=M∥V∥222+μ∥v2±−v1±∥222. Subtracting the minus identity from the plus identity cancels the centre-of-mass term, and the restitution law gives ∥v2+−v1+∥22=e2∥v2−−v1−∥22 by homogeneity of the norm [L2], so ΔK=μ2(e2−1)∥v2−−v1−∥22=−1−e22μ∥v2−−v1−∥22≤0 because 0≤e≤1; the cases e=1 and e=0 are the stated substitutions. This is claim c.

5 · Examples, experiments and counterexamples

None yet.

Sources