V-MMS-02 — ALE free-stream preservation under prescribed mesh motion

Tags: mms, ale, gcl, free_stream

References:

  • Thomas & Lombard, AIAA J. 17 (1979) — geometric conservation law
  • docs/VerificationSuitePlan_2026-07-06.md Tier 3 (ALE-MMS deliverable)

Problem statement

The sharpest possible discrete-GCL assert: spatially uniform T, ξ with no reactions and zero-flux manufactured BCs, while the mesh is compressed to 70 % of its initial thickness through the Workspace.w_override hook with the uniform-stretch velocity

w*(z, t) = (ṡ/s)·z, s(t) = 1 − 0.15·(1 − cos(2πt/20)), z_bottom pinned

(node trajectories zi(t) = zi(0)·s(t) exactly). A uniform state must stay uniform to machine precision on the moving mesh: the ALE advection of a constant field must vanish discretely (limiter slope ≡ 0), the dilation correction is zero (θ reaction-driven, no reactions), and the conductive/ species fluxes are zero. Any drift is a moving-mesh transport defect that smooth-field order studies (V-MMS-03) can hide inside their truncation error.

The node-motion ODE żi = w*(zi,t) is integrated by the same stiff integrator, so the top-cell-center trajectory doubles as an exactness check on the prescribed-motion path (temporal tolerance only).

Jacobian: colored/dense FD (jacobian = :fd) — the structured geometry Jacobian assumes reaction-driven mesh velocity and does not model the override (documented on Workspace.w_override).

Quantities of interest (n = 20)

QoIvalueexacterrortolerancewithin tolprovenance
free-stream T drift (all cells, all snapshots)5.684e-1405.684e-141e-10yesuniform field must survive 30% mesh compression
free-stream ξ_solid relative drift1.895e-1601.895e-161e-12yesALE advection + dilation must vanish for a uniform state
free-stream ξ_gas relative drift2.22e-1502.22e-151e-12yesALE advection + dilation must vanish for a uniform state
top cell-center trajectory vs z₀·s(t_end)0.0068250.0068254.039e-071.5e-06yesnode ODE ż = (ṡ/s)z integrated at solver tolerance

V-MMS-02 thickness_history

V-MMS-02 thickness_history_error

Solver configuration

settingvalue
integratorKenCarp4 (default)
abstol1.0e-10
reltol1.0e-8