This manuscript presents a finite-time nonlinear numerical study of a two-sided Ellis wormhole in Einstein–æther gravity with direct matter–æther coupling. A genuinely nonspherical perturbation is evolved using a 53-component first-order constrained system on two overlapping stereographic patches. The study emphasizes numerical verification, refinement behavior, constraint control, and the response of both the matter sector and the frame–æther sector.
The main result is that the large quadratic response generated in the frame–æther sector decreases substantially with spatial refinement, whereas the response of the matter sector remains comparatively small, and the sampled wormhole throat shows no discernible large-scale geometric motion over the simulated time interval. The calculations do not establish nonlinear stability or continuum convergence; rather, they provide a finite-time benchmark for directly coupled matter–preferred-frame dynamics in modified gravity and identify the higher-resolution tests required for a continuum-level interpretation.
The record includes the manuscript, numerical tables, figures, machine-readable audit summaries, and original computational source code in Python, MATLAB, and C++.