Physlib

PhyslibAlpha.ClassicalMechanics.NortonDome.PhysicalSpace

The Norton dome in physical space

i. Overview

`NortonDome.Basic` writes the dynamics on the arc length `r`, with kinetic energy `½ m ṙ²` and potential energy `-m g h(r)`. This module identifies that chart dynamics with a point mass constrained to the surface of the dome, as `SimplePendulum.Geometric.PhysicalSpace` does for the pendulum.

The motion is radial, so it suffices to work in the vertical plane through the apex, with the apex at the origin and the vertical coordinate measured upwards. The profile is the arc-length parametrised curve `s ↦ (x(s), -h(s))`, where `x'(s)² + h'(s)² = 1` and `h'(s) = √s / g` give `x'(s) = √(1 - s/g²)` and `x(s) = (2g²/3) (1 - (1 - s/g²)^{3/2})`. The profile exists for `0 ≤ s ≤ g²`; at `s = g²` it is vertical and Norton's model stops. While the arc length is in that range the chart kinetic energy is `½ m ‖v‖²` by unit speed, and the chart potential energy is the gravitational potential `m g y` at height `y`, so the chart Lagrangian is the constrained Lagrangian.

ii. Key results

- `NortonDome.horizontal` is the horizontal distance `x(s)` of the profile from the apex, with its derivative `hasDerivAt_horizontal`; `NortonDome.unit_speed` is the identity `x'(s)² + h'(s)² = 1`. - `NortonDome.profile` is the profile curve `s ↦ (x(s), -h(s))` in the vertical plane, and `NortonDome.dome` its image on `0 ≤ s ≤ g²`. - `NortonDome.spaceTrajectory` is the position of the particle in the plane along a chart trajectory; it lies on the dome, `spaceTrajectory_mem_dome`, while the arc length is in range. - `NortonDome.deriv_spaceTrajectory` is the velocity of the particle in the plane, with the square of its speed `norm_sq_deriv_spaceTrajectory`. - `NortonDome.kineticEnergy_eq_space`, `NortonDome.potentialEnergy_eq_height`, `NortonDome.lagrangian_eq_space` and `NortonDome.energy_eq_space` identify the chart energies and Lagrangian with those of the point mass in physical space.

iii. Table of contents

- A. The profile of the dome in the plane - A.1. The horizontal coordinate and the unit-speed property - A.2. The profile curve and the dome - B. The particle's position along a trajectory - C. The particle's velocity and speed - D. The energies and the Lagrangian in physical space

iv. References

- Norton, J. D., *The dome: an unexpectedly simple failure of determinism*, Philosophy of Science 75 (2008), 786–798. - `Physlib.ClassicalMechanics.Pendulum.SimplePendulum.Geometric.PhysicalSpace` (the same identification for the pendulum).

A. The profile of the dome in the plane

A.1. The horizontal coordinate and the unit-speed property

The horizontal distance from the apex has derivative `√(1 - s/g²)`; with the slope `√s / g` of the depth this gives unit speed.

A.2. The profile curve and the dome

B. The particle's position along a trajectory

Along a chart trajectory `r` the particle is at the point of the profile at arc length `r t 0`, on the dome while that is in range.

C. The particle's velocity and speed

The velocity in the plane is the unit tangent `(x'(r), -h'(r))` times `ṙ`, so by unit speed the square of the speed is `ṙ²` while the arc length is in range.

D. The energies and the Lagrangian in physical space

The chart kinetic energy is that of the particle in the plane while the arc length is in range, and the chart potential energy is `m g y` at its height `y` at all times; the chart Lagrangian is the constrained Lagrangian `T - V`.

23 declarations

definition

Horizontal distance x(s)x(s) of the Norton dome profile

For a Norton dome characterized by a parameter gg (representing gravitational acceleration), the horizontal distance x(s)x(s) from the apex at a given arc length ss is defined by the function: x(s)=2g23(1−(1−sg2)3/2) x(s) = \frac{2g^2}{3} \left( 1 - \left( 1 - \frac{s}{g^2} \right)^{3/2} \right) The profile of the dome is valid for arc lengths 0≤s≤g20 \le s \le g^2.

theorem

Explicit formula for the horizontal distance x(s)x(s) of the Norton dome profile

For a Norton dome with gravitational acceleration gg, the horizontal distance x(s)x(s) from the apex at arc length ss is given by the formula: x(s)=2g23(1−(1−sg2)3/2) x(s) = \frac{2g^2}{3} \left( 1 - \left( 1 - \frac{s}{g^2} \right)^{3/2} \right)

theorem

x(0)=0x(0) = 0 for the Norton dome profile

For a Norton dome, the horizontal distance x(s)x(s) from the apex evaluated at arc length s=0s = 0 is 00.

theorem

Derivative of the horizontal distance dxds=1−sg2\frac{dx}{ds} = \sqrt{1 - \frac{s}{g^2}}

Let x(s)x(s) be the horizontal distance of the Norton dome profile from the apex as a function of the arc length ss, and let gg denote the characteristic gravitational parameter of the dome. For any s∈Rs \in \mathbb{R}, the derivative of the horizontal distance with respect to the arc length is given by: dxds=1−sg2 \frac{dx}{ds} = \sqrt{1 - \frac{s}{g^2}}

theorem

The horizontal distance function x(s)x(s) is differentiable

Let x(s)x(s) be the horizontal distance of the Norton dome profile from the apex, considered as a function of the arc length s∈Rs \in \mathbb{R}. The function xx is differentiable on R\mathbb{R}.

theorem

Unit Speed Identity x′(s)2+h′(s)2=1x'(s)^2 + h'(s)^2 = 1 for the Norton Dome Profile

For any real number ss representing the arc length such that 0≤s≤g20 \le s \le g^2, where gg is the characteristic parameter of the Norton dome, the following identity holds: (1−sg2)2+(sg)2=1\left( \sqrt{1 - \frac{s}{g^2}} \right)^2 + \left( \frac{\sqrt{s}}{g} \right)^2 = 1 In the context of the dome's geometry, this identity represents the unit-speed property x′(s)2+h′(s)2=1x'(s)^2 + h'(s)^2 = 1, where x′(s)=1−s/g2x'(s) = \sqrt{1 - s/g^2} is the derivative of the horizontal distance and h′(s)=s/gh'(s) = \sqrt{s}/g is the derivative of the vertical depth with respect to the arc length ss.

definition

Profile curve s↦(x(s),−h(s))s \mapsto (x(s), -h(s)) of the Norton dome

The profile of the Norton dome is a curve in R2\mathbb{R}^2 (representing the vertical plane through the apex) parameterized by the arc length s∈Rs \in \mathbb{R}. The position at arc length ss is given by the vector (x(s),−h(s))(x(s), -h(s)), where x(s)x(s) is the horizontal distance from the apex and h(s)h(s) is the depth of the dome surface below the apex.

theorem

The first component of the Norton dome profile is the horizontal coordinate x(s)x(s)

For any real number ss representing the arc length, the first component (index 0) of the Norton dome's profile curve P(s)P(s) is equal to the horizontal distance x(s)x(s) from the apex. Specifically, if the profile is given by P(s)=(x(s),−h(s))P(s) = (x(s), -h(s)), then P(s)0=x(s)P(s)_0 = x(s).

theorem

The vertical coordinate of the Norton dome profile is −h(s)-h(s)

For any arc length s∈Rs \in \mathbb{R}, the vertical coordinate (the second component) of the Norton dome's profile curve at ss is equal to −h(s)-h(s), where h(s)h(s) is the depth of the dome surface below the apex.

theorem

profile(0)=0\text{profile}(0) = 0

The profile curve profile(s)\text{profile}(s) of the Norton dome evaluated at arc length s=0s = 0 is the origin 00. This identifies the apex of the dome with the origin of the coordinate system.

definition

The Norton dome as the image of the profile curve on [0,g2][0, g^2]

The Norton dome is defined as the subset of the vertical plane R2\mathbb{R}^2 through its apex formed by the image of the profile curve profile(s)\text{profile}(s) for arc lengths ss in the closed interval [0,g2][0, g^2]. That is, dome={profile(s)∣0≤s≤g2}\text{dome} = \{ \text{profile}(s) \mid 0 \le s \le g^2 \}, where gg is the gravitational acceleration parameter.

theorem

profile(s)∈dome\text{profile}(s) \in \text{dome} for 0≤s≤g20 \le s \le g^2

For any real number ss representing the arc length such that 0≤s≤g20 \le s \le g^2, where gg is the gravitational acceleration, the point profile(s)\text{profile}(s) on the profile curve of the Norton dome is an element of the set dome\text{dome}.

definition

Physical space trajectory from arc length trajectory on the Norton dome

The function maps a chart trajectory r(t)r(t), which represents the arc length of a particle on the Norton dome as a function of time, to its corresponding trajectory in 2-dimensional physical space. At each time tt, the position in the vertical plane is given by the profile curve evaluated at the current arc length r(t)0r(t)_0, i.e., t↦profile(r(t)0)t \mapsto \text{profile}(r(t)_0).

theorem

The horizontal component of the physical trajectory is x(r(t)0)x(r(t)_0)

For a chart trajectory r(t)r(t) representing the arc length s=r(t)0s = r(t)_0 of a particle on a Norton dome at time tt, the horizontal component (the 0-th coordinate) of the particle's physical position vector in 2D space is equal to the horizontal distance function x(s)x(s) evaluated at the current arc length. That is, (spaceTrajectory r)(t)0=x(r(t)0) (\text{spaceTrajectory } r)(t)_0 = x(r(t)_0) where x(s)=2g23(1−(1−sg2)3/2)x(s) = \frac{2g^2}{3} \left( 1 - \left( 1 - \frac{s}{g^2} \right)^{3/2} \right) is the horizontal coordinate of the dome's profile.

theorem

Vertical component of physical trajectory is −h(r)-h(r)

For a trajectory r(t)r(t) representing the arc length of a particle on the Norton dome at time tt, the vertical coordinate (the second component) of its position in 2D physical space is the negative of the depth of the dome at that arc length: (spaceTrajectory(r,t))1=−h(r(t)0)(\text{spaceTrajectory}(r, t))_1 = -h(r(t)_0) where hh is the depth function of the Norton dome and r(t)0r(t)_0 is the scalar value of the arc length at time tt.

theorem

The physical trajectory lies on the dome for arc lengths in [0,g2][0, g^2]

Let r(t)r(t) be a trajectory representing the arc length of a particle as a function of time tt. If at time tt the arc length r(t)r(t) satisfies 0≤r(t)≤g20 \le r(t) \le g^2 (where gg is the gravitational acceleration), then the particle's position in 2-dimensional physical space, given by the mapping spaceTrajectory(r,t)\text{spaceTrajectory}(r, t), is an element of the set defining the Norton dome.

theorem

Differentiability of the Physical Space Trajectory on the Norton Dome

Let r:Time→Rr: \text{Time} \to \mathbb{R} be an arc length trajectory on the Norton dome. If rr is differentiable with respect to time, then the corresponding physical trajectory in 2-dimensional space, given by t↦spaceTrajectory(r,t)t \mapsto \text{spaceTrajectory}(r, t), is also differentiable with respect to time.

theorem

Velocity of the physical trajectory on the Norton dome **End of informal name**

Let r:Time→Rr: \text{Time} \to \mathbb{R} be a differentiable arc length trajectory of a particle on a Norton dome, and let gg denote the gravitational acceleration. The velocity of the particle in the vertical 2D plane at time tt is the derivative of its physical space trajectory, given by: ddtspaceTrajectory(r,t)=(1−r(t)g2r˙(t),−r(t)gr˙(t))\frac{d}{dt} \text{spaceTrajectory}(r, t) = \left( \sqrt{1 - \frac{r(t)}{g^2}} \dot{r}(t), -\frac{\sqrt{r(t)}}{g} \dot{r}(t) \right) where r˙(t)\dot{r}(t) is the rate of change of the arc length with respect to time. This velocity vector corresponds to the unit tangent vector to the dome's profile curve multiplied by the scalar speed r˙\dot{r}.

theorem

The Squared Physical Speed Equals the Squared Derivative of the Arc Length for the Norton Dome

Let r(t)r(t) be a differentiable trajectory representing the arc length of a particle on the Norton dome. If at time tt the arc length is within the valid range 0≤r(t)≤g20 \le r(t) \le g^2, where gg is the acceleration due to gravity, then the square of the speed of the particle in 2D physical space (the norm of the time derivative of the physical trajectory) is equal to the square of the time derivative of the arc length: ∥ddtx(t)∥2=(r˙(t))2\|\frac{d}{dt} \mathbf{x}(t)\|^2 = (\dot{r}(t))^2 where x(t)\mathbf{x}(t) is the position of the particle in the vertical plane (the `spaceTrajectory`).

theorem

Chart Kinetic Energy Equals Physical Kinetic Energy for the Norton Dome

Let r:Time→R1r: \text{Time} \to \mathbb{R}^1 be a differentiable trajectory representing the arc length of a particle of mass mm on a Norton dome with gravitational acceleration gg. For any time tt such that the arc length s=r(t)0s = r(t)_0 is within the valid range 0≤s≤g20 \le s \le g^2, the chart kinetic energy is equal to the kinetic energy in 2D physical space: T(t)=12m∥ddtx(t)∥2 T(t) = \frac{1}{2} m \left\| \frac{d}{dt} \mathbf{x}(t) \right\|^2 where x(t)\mathbf{x}(t) (the `spaceTrajectory`) is the position of the particle in the vertical plane.

theorem

Chart Potential Energy V=mgyV = m g y for the Norton Dome Trajectory

For any trajectory r(t)r(t) representing the arc length of a particle on the Norton dome, and for any time tt, the chart potential energy V(r(t))V(r(t)) is equal to the gravitational potential energy in physical space, mgy(t)m g y(t), where mm is the mass, gg is the gravitational acceleration, and y(t)y(t) is the vertical coordinate (height) of the particle in 2D physical space at time tt.

theorem

Chart Lagrangian of the Norton Dome Equals the Physical Space Lagrangian

Consider a particle of mass mm on a Norton dome with gravitational acceleration gg. Let r(t)r(t) be a differentiable trajectory representing the arc length of the particle from the apex. For any time tt such that the arc length s=r(t)0s = r(t)_0 is within the range 0≤s≤g20 \le s \le g^2, the chart Lagrangian LL evaluated at the current state (t,s,s˙)(t, s, \dot{s}) is equal to the constrained Lagrangian of the particle in physical space: L(t,s,s˙)=12m∥v∥2−mgyL(t, s, \dot{s}) = \frac{1}{2} m \|\mathbf{v}\|^2 - m g y where v\mathbf{v} is the velocity vector ddtq(t)\frac{d}{dt} \mathbf{q}(t) of the particle in 2D physical space and yy is its vertical coordinate (height).

theorem

Chart Energy Equals Physical Space Energy for the Norton Dome

Consider a particle of mass mm on a Norton dome under gravitational acceleration gg. Let r(t)r(t) be a differentiable trajectory where r(t)0r(t)_0 represents the arc length from the apex at time tt. For any time tt such that the arc length s=r(t)0s = r(t)_0 is within the valid range 0≤s≤g20 \le s \le g^2, the total chart energy E(t)E(t) is equal to the total mechanical energy of the particle in 2D physical space: E(t)=12m∥v(t)∥2+mgy(t)E(t) = \frac{1}{2} m \|\mathbf{v}(t)\|^2 + m g y(t) where v(t)\mathbf{v}(t) is the velocity vector of the particle in physical space (the derivative of the spatial trajectory) and y(t)y(t) is its vertical coordinate (height).