TNB frames

When a jet fighter or a rocket travels along a complex 3D flight path, analyzing its motion using standard fixed Cartesian coordinates is incredibly inefficient. Instead, engineers map the vehicle's telemetry directly onto its moving TNB (Frenet-Serret) frame.

• Tangential Acceleration points entirely along the Tangent vector (T). This component measures how much the vehicle's speed is changing (e.g., rocket engines firing or air brakes deploying).
• Normal Acceleration points entirely along the Normal vector (N). This component measures how much the vehicle's direction is changing due to local curvature (κ).
• There is zero acceleration along the Binormal vector (B). This mathematical absolute provides immense engineering utility:
  1. The Osculating Plane: The vectors T and N span a plane called the osculating plane. Because acceleration rests completely within this plane, the vehicle's instantaneous turning maneuver is always localized to a 2D slice of space.
  2. Tracking Pilot G-Forces: The Binormal vector (B = T × N) acts as the dynamic axle of the turn. When an aircraft executes a sharp bank, the pilot experiences lateral forces. By evaluating Torsion (τ) along the flight path, aerospace software computes how fast the plane's plane-of-motion is twisting. This is used to ensure a missile or aircraft doesn't execute a maneuver that structurally tears its fuselage apart or causes a pilot to experience a fatal G-induced loss of consciousness (G-LOC, stands for G-force induced loss of consciousness}. It is a dangerous medical state where the brain is temporarily starved of oxygen (cerebral hypoxia) because high acceleration forces drain blood downward into the lower body).

A fiber optic cable consists of a pure silica core that transmits data via Total Internal Reflection (TIR). If a fiber is bent or twisted too aggressively, it suffers from two physical degradation mechanisms:

1. Macrobending Loss: Light leaks out of the core because the angle of incidence drops below the critical angle required for TIR.
2. Polarization Mode Dispersion (PMD): Physical twisting induces mechanical stress inside the glass cylinder, breaking its geometric symmetry. This creates birefringence, where different polarizations of light travel at different speeds, smearing the optical pulses and causing data corruption.

• Curvature (κ) Limits: The radius of curvature R = 1/κ cannot drop below the fiber's structural minimum bending radius. If κ spikes, macrobending losses skyrocket.
• Torsion (τ) Limits: Torsion measures how fast the Osculating Plane (spanned by T and N) rotates around the Tangent vector (T). This structural "twist rate" directly correlates with the torsional shear stress applied to the silica. Engineers use τ to calculate the mechanical twist limits during the high-speed spinning process of manufacturing ribbon cables.

The following Mathematica code models a fiber optic thread spiraling inside a cable that is navigating a sweeping 3D bend. It extracts and displays the physical geometric stressors (κ and τ) along the path.

DNA (deoxyribonucleic acid) is the hereditary molecule that carries the genetic instructions used in the development, functioning, growth, and reproduction of all known living organisms

Inside a biological cell, DNA is not an unconstrained, relaxed ladder; it is packed under immense mechanical tension. When enzymes like RNA Polymerase track along the DNA to transcribe genes, they forcefully unwind the double helix ahead of them. This action introduces severe mechanical strain upstream and downstream, a phenomenon known as DNA Supercoiling.

To quantify this, biophysicists model the DNA backbone as a continuous space curve and apply the Călugăreanu-White-Fuller Theorem, which elegantly maps topological properties to the TNB frame parameters:

• Linking Number (Lk): An invariant topological integer counting how many times the two strands intertwine. It can only change if a cell enzyme physically cuts the DNA backbone.
• Twist (Tw): The local rotation of the DNA strands around the central axis. This is calculated directly by integrating the total torsion (τ) along the curve, combined with the normal vector rotation.
• Writhe (Wr): A measure of the non-local, spatial coiling of the central axis itself (how the DNA loops over itself like an old telephone cord). Writhe is governed strictly by the curvature (κ) profile and non-local geometry.

The following Mathematica code surves to Visualize DNA Supercoiling Energy:

• s is the arc length parameter.
• A is the Bending Stiffness Coefficient (also called flexural rigidity) (for DNA, this relates to persistence length; for glass fiber, it depends on Young’s Modulus).
• B is the Torsional/Twisting Stiffness Coefficient (governed by Shear Modulus).
• τ₀ is the intrinsic torsion of the material when completely relaxed (e.g., the natural twist rate of DNA base pairs).

1. Bending Yield Point (κ_max): Exceeding this causes micro-fractures in glass fibers or ruptures base-pair bonds in DNA.
2. Torsional Buckling Instability (τ_max): When a thread accumulates too much torsion, it undergoes a geometric phase transition called writhing—buckling directly into a tight plectonemic loop to lower its twist energy.

Notice that during integration, the density equations are multiplied by speed[t] (which represents ds/dt). This converts parameter integration directly into physical arc-length integration seamlessly, keeping calculations geometrically precise regardless of how fast or slow the curve's parametrization is pacing.    ■

1. Plot the TNB frame for the spacecurve given by r(t) = [t, t cos(3t), t sin(3t)] with 0 ≤ t ≤ 4π.
2. Plot the TNB frame for the spacecurve given by r(t) = [cos(2 t), sin(3 t), cos(3 t)] with 0 ≤ t ≤ 4π.
3. Plot the TNB frame for the spacecurve given by r(t) = [(2 + cos(2t)) cos(3t), (2 + cos(2t)) sin(3t), sin(4t)] with 0 ≤ t ≤ 3π.
4. The following space curve lies on a torus and is given by r(t) = [(3 + 2 cos(t)) cos(10t), (3 + 2 cos(t)) sin(10t), 2 sin(t)] with 0 ≤ t ≤ 2π. Plot, in the same figure, the following graphs: r(t), the TNB frame for r(t), and the torus, given parametrically by s(u, v) = [(3 + 2 cos(v)) cos(u), (3 + 2 cos(v)) sin(u), 2 sin(v)], for −π ≤ u, v ≤ π.
5. The space curve r(t) = [(3 + 2 cos(3t)) cos(7t), (3 + 2 cos(3t)) sin(7t), 2 sin(3t)] for 0 ≤ t ≤ 2π, also lies on a torus given parametrically by s(u, v) = [(3 + 2 cos(v)) cos(u), (3 + 2 cos(v)) sin(u), 2 sin(v)], for −π ≤ u, v ≤ π. Plot, in the same figure, r(t), the TNB frame for r(t), and the torus.
6. Plot the curvature function κ(t) for the space curve from the previous problklem.
7. Consider the standard circular helix parametrized by r(t) = [𝑎 cos(t), 𝑎 sin(t), bt], where 𝑎 > 0 and b > 0 are constants.
   1. Find the unit tangent vector T(t).
   2. Find the principal normal vector N(t) and show that it always points directly toward the z-axis.
   3. Compute the binormal vector B(t).
   4. Show that both the curvature κ(t) and torsion τ(t) are completely constant for all t.
8. Let r(t) = [t, t², 0] be a parabola constrained to the xy-plane.
   1. Compute the TNB frame vectors at the vertex t = 0.
   2. Calculate the torsion τ(t) for any arbitrary t.
   3. Conceptual question: What must the torsion always be for any curve that lies completely flat within a single plane? Why?
9. A drone is tracking a path modeled by \( \displaystyle \quad \mathbf{r}(t) = \left[ e^t , \ \sqrt{2}\,t, \ e^{-t} \right] . \)
   1. Compute the velocity vector v(t) and the acceleration vector 𝑎(t).
   2. Decompose the acceleration vector into its tangential component (𝑎_T(t)) and normal component (𝑎_N(t)) such that 𝑎 = 𝑎_TT + 𝑎_NN.
   3. If the drone's structural hull fails when the directional bending acceleration ( 𝑎_N) exceeds 2.0 m/s², determine if the drone survives its maneuver at t = 0.
10. An engineer is pulling a fiber optic cable through a curved machine guide. The line's trajectory follows Viviani's curve r(t) = [1 + cost, sint, 2 sin(t/2)]
    1. Find the points along the interval t ∈ [0, 4π] where the torsion τ(t) = 0.
    2. Explain what is happening to the moving osculating plane at these exact coordinates. Is the fiber twisting or momentarily flattening out?
11. Define the (4,3) toroidal knot: \[ \mathbf{r}(t) = \left[ (3 + \cos (4t))\,\cos (3t) , \ (3 + \cos (4t))\,\sin (3t) , \ \sin (4t) \right] . \]
    1. Write a script to calculate and plot the curvature κ(t) and torsion τ(t) side-by-side for t ∈ [0, 2π].
    2. By inspecting your generated graphs, visually identify how many distinct curvature "peaks" (local maxima) occur over one full revolution. What physical geometric property of the torus causes this oscillation?

1. Calladine, C. R., Drew, H. R., Luisi, B. F., & Travers, A. A. (2004). Understanding DNA: The Molecule and How it Works (3rd ed.). Academic Press.
2. Do Carmo, M. P. (2016). Differential Geometry of Curves and Surfaces (2nd ed.). Dover Publications.
3. Ghatak, A., & Thyagarajan, K. (1998). Introduction to Fiber Optics. Cambridge University Press.
4. Greenwood, D. T. (1988). Principles of Dynamics (2nd ed.). Prentice Hall.
5. Stewart, J. (2020). Calculus: Early Transcendentals (9th ed.). Cengage Learning.
6. White, J. H. (1969). Self-linking and the Gauss integral in higher dimensions. American Journal of Mathematics, 91(3), 693-728.