Viscoelastic stress past the conformation limit: a positive-definite convected step and a log-conformation history

[lmoresi]

Stacked on #789. Lets the viscoelastic Stokes solve reach the high-Weissenberg regime of the benchmark flows: with it the creeping UCM cross-slot shows its purely elastic pitchfork, where before every run from De 0.5 stalled or hung.

What lost the conformation, measured on the cross-slot at De 0.8:

• The step. The upper-convected BDF-1 step is F c* Fᵀ (F = I + dt L) less its second-order term dt² L c* Lᵀ, so it goes indefinite once dt|L| is of order one. A smaller timestep did not help: at a corner |L| grows with resolution.
• The store. Interpolation, projection and per-cell fits undershoot near a corner singularity, and the forward store's linear fit extrapolates to the cell edges; the stress read back was indefinite although every stored value was positive.

Two options on ViscoElasticPlasticFlowModel (upper-convected, first order):

• convected_step="deformation" keeps the dropped term: positive-definite for any step, still first order. The exponential step completes the relaxation target's stretch to a product the same way.
• stress_history="log_conformation" (2-D) stores ψ = log c and reads σ* = G(exp ψ* − I); it implies the deformation step. The flavours are untouched: each transports its stored tensor by pure advection, and the step stays on c, where relaxation is linear, so neither integrator changes.

The model owns the encoding. The default NavierStokes decodes its history flux; SNES_NavierStokes and the multi-material model refuse the log store. stokes.tau projects the decoded stress; the health check reports where the logarithm's floor acted.

Result, cross-slot, full resolution, forward history, seeded: growth rate of the asymmetry −0.77, −0.24, +0.06, +0.27 at De 0.5, 0.65, 0.72, 0.8 (De = λU/H, H the half-width), onset near 0.71, against 0.62 expected from Poole et al. (full-width 0.31). The conformation stays above 0.33 everywhere.

Tests (test_1065, and log cases in test_1064): the one-step extension closed form for both integrators and both representations, plus the linear step's indefinite −0.818; the shear start-up across four flavours and both integrators against the discrete recurrence in σ_xy and N1; the units model. The stress-history, swarm and advection files pass (316); the style gate is clean.

Limits: first order; the log history 2-D only; not for stresses that are not conformations (geodynamic VEP with no objective rate).

Underworld development team with AI support from Claude Code

🤖 Generated with Claude Code

https://claude.ai/code/session_017kSkAq7oJ5J3XisuvLBovo

[lmoresi]

Adversarial review before the push (two reviewers: the units change and this branch, briefed for over-complexity as well as defects). Applied here:

• One option, not a matrix. The log history implies the deformation step and refuses the linear one; the linear step could hand it an indefinite conformation, which the floor then hid.
• One owner of the encoding. The model's encode_history and its carried-stress read; one solver helper stores the record for both the flux and the commit (the fallback chain was written three ways); the history gets the encoding once, for the inflow datum, read at one place. Removed: the private maps pushed into the history, the _inflow_value property and its __dict__ access, three of the five encode/decode methods.
• Readers that took ψ as a stress. The default NavierStokes flux (decodes now), stokes.tau (projects the decoded stress, not the next step's predictor), the multi-material model (refuses).
• The exponential term at the elastic limit. Written in x = dt/λ; the integrator's α is clamped to one there and the old form divided by zero.
• The floor. 1e-12 (round-off only after a positive-definite step), not 1e-3, which clipped legitimate conformations such as 1/(1+2Wi) in extension; the health check now reports where it acted, since a decoded conformation is positive by construction.
• max_elastic_timestep reports no limit for the deformation step; the order guard sits where the order is set.
• Tests. The shear box checked only σ_xy, where every new term vanishes; it now follows the discrete recurrence in N1 too, and the exponential deformation term has its closed form.

Recorded: set_initial_history takes stored values (documented), while the inflow datum is given as a stress; a checkpoint does not record which representation it holds.