Abstract
This paper details the mathematical and computational architecture of a multi-domain physics engine designed to validate interrelated aerospace systems, specifically bi-ionic vector thrust and solid-state thermoelectric energy recycling. To eliminate computational hallucination and enforce strict physical boundaries, the engine is constructed utilizing highly deterministic Python libraries: SciPy for immutable, universally accepted physical constants; SymPy for strict, symbolic equation rendering and mathematical proofs; and NumPy for continuous spatial array generation. This framework interlinks discrete physical phenomena (magnetic containment, plasma gyration, kinetic wave damping, and thermal energy harvesting) into a singular, codependent mathematical ecosystem. The novel contribution is the Master Resonance Equation (Omega), which forces established institutional physics equations to act as strict boundary conditions for the system as a whole.
Foundational Mathematical Lineage The computational logic of this system relies entirely on accepted, peer-reviewed physics formulas. The engine mathematically proves that the hardware parameters do not violate the following established scientific models:
-
Magnetic Induction (Biot-Savart Law)
-
Institutional Baseline: MIT OpenCourseWare (Physics II: Electricity and Magnetism).
-
Integration: The core magnetic field (B) generated by the Bi-Ionic Golden Ratio coil is derived using the standard approximation for the center of a current loop: B = (mu_0 * N * I) / (2 * R_0)
-
System Function: Determines the strict baseline containment strength for the subsequent plasma injection.
-
Plasma Gyration & Containment (Lorentz Force Law)
-
Institutional Baseline: HyperPhysics, Georgia State University (Department of Physics and Astronomy).
-
Integration: The Larmor radius (r_L) of the injected ions dictates the physical width of the plasma orbit: r_L = (m_p * v) / (e * B)
-
System Function: This equation is cross-referenced against the physical coil radius (R_0). If the math dictates that r_L is greater than or equal to R_0, the system registers a catastrophic physical containment breach.
-
Vector Thrust (Newtonian Mass Flow)
-
Institutional Baseline: NASA Glenn Research Center (General Thrust Equation).
-
Integration: Continuous electric/plasma thrust (F_thrust) is calculated as the product of the propellant mass flow rate (m_dot) and the perpendicular exhaust velocity (v): F_thrust = m_dot * v
-
System Function: Proves the exact mechanical force output of the stabilized plasma confinement structure.
-
Metric Wave Stabilization (D'Alembert Wave Operator)
-
Institutional Baseline: LibreTexts Physics (Damped Harmonic Oscillators and Wave Mechanics).
-
Integration: The energy bled from the kinetic system is modeled as a damped metric wave, driven by the damping ratio (Zeta): Zeta = gamma / (2 * c)
-
System Function: Identifies the precise rate of energy dissipation. An overdamped system (Zeta >= 1.0) indicates instant wave collapse, while an underdamped state preserves kinetic ringing.
-
Solid-State Thermoelectric Harvesting (Seebeck Effect)
-
Institutional Baseline: ScienceDirect / Elsevier (Materials Science Database).
-
Integration: The kinetic and thermal waste mapped by the damping ratio is harvested via Bismuth Telluride thermoelectric structures, measured by the dimensionless Figure of Merit (ZT): ZT = ((S^2 * sigma) / kappa) * T_mean
-
System Function: Proves thermodynamic viability. The system enforces a hard physical limit; calculations yielding ZT > 3.0 trigger thermodynamic failure protocols based on current absolute material science limits.
Master System Resonance and Limit Enforcement The equations detailed above are standard, disparate elements of physics. The engine's novel advancement is the unification of these formulas into a singular diagnostic fitness function, designated as the Master Resonance Equation (Omega_system):
Omega_system = (1 - (r_L / R_0)) * (ZT / Zeta) * ln(F_thrust + 1)
This equation acts as the final mathematical gatekeeper by forcing the accepted science to self-regulate the system:
-
Containment Limit [1 - (r_L / R_0)]: If the HyperPhysics Lorentz calculation (r_L) exceeds the MIT Biot-Savart boundaries (R_0), the multiplier drops to zero or becomes negative, resulting in a mathematical system failure.
-
Recycling Limit [ZT / Zeta]: Demands that the Elsevier thermodynamic efficiency (ZT) outpaces the LibreTexts energy damping ratio (Zeta).
-
Thrust Scaling [ln(F_thrust + 1)]: Applies a natural logarithm to the NASA mass flow output, preventing massive Newtonian thrust numbers from masking or overriding failing containment or thermal metrics.
By mathematically binding these institutional equations together, this solver guarantees that any positive Resonance (Omega > 0) strictly adheres to the accepted laws of thermodynamics and electrodynamics. Use and application of these interrelated mathematical models are subject to strict field of use licensing parameters.
Creative Commons License

This work is licensed under a Creative Commons Attribution-Noncommercial-No Derivative Works 3.0 License.
Recommended Citation
Schramm, Daniel, "Deterministic Multi-Domain Computational Framework for Bi-Ionic Propulsion and Thermoelectric Resonance (STALWART v18.1)", Technical Disclosure Commons, ()
https://www.tdcommons.org/dpubs_series/11538
Screenshot 2026-08-29 213108.png (125 kB)
Screenshot 2026-08-29 213116.png (19 kB)