FP 20-Layer Unidirectional CF Layup (dim=980)
PARITYS2 · dim 980No clear winner. The survival gap is under 10 percentage points and the balanced-score gap is under 0.05, so neither SolvSRK nor the best baseline clears the win threshold. Either works - choose on cost, licensing, or integration effort. All verdicts →
20-layer unidirectional carbon-fiber laminate with through-thickness 1D thermal diffusion coupled to Kamal-Sourour autocatalytic cure kinetics. 490 spatial nodes, state = [T(490), alpha(490)]. Bottom face Dirichlet at 453 K, top face convective BC. Stiffness from fine spatial resolution (dx ~ 20 um) vs slow cure timescale.
Problem definition
Kamal & Sourour (1973); Pojman (2012)
Canonical RHS excerpt from the registered callable used for this benchmark cell. Expand it to verify the state equations; it is not a standalone runnable fixture.
Show canonical RHS excerpt
def _fp_layup_20_kamal_rhs(t, y):
n = _FL20_N_NODES
T = np.clip(y[:n], _T_FLOOR, _T_CEIL)
alpha = np.clip(y[n:], 0.0, 1.0)
inv_RT = 1.0 / (_R_GAS * T)
arg1 = np.clip(_FL20_E1 * inv_RT, 0.0, _EXP_ARG_MAX)
arg2 = np.clip(_FL20_E2 * inv_RT, 0.0, _EXP_ARG_MAX)
k1 = _FL20_A1 * np.exp(-arg1)
k2 = _FL20_A2 * np.exp(-arg2)
dadt = (k1 + k2 * np.power(alpha, _FL20_M)) * np.power(1.0 - alpha, _FL20_N_ORD)
dT = np.empty(n)
# Interior nodes — vectorised Laplacian
dT[1:-1] = (
_FL20_DIFF * _FL20_INV_DX2 * (T[:-2] - 2.0 * T[1:-1] + T[2:])
+ _FL20_SRC * dadt[1:-1]
)
# Bottom node: Dirichlet (held at cure temperature)
dT[0] = 0.0
# Top node: convective BC → ghost = T[-1] + (h·dx/κ)·(T_amb - T[-1])
T_ghost = T[-1] + (_FL20_H_CONV * _FL20_DX / _FL20_KAPPA) * (_T_AMBIENT - T[-1])
dT[-1] = (
_FL20_DIFF * _FL20_INV_DX2 * (T[-2] - 2.0 * T[-1] + T_ghost)
+ _FL20_SRC * dadt[-1]
)
dy = np.empty(_FL20_DIM)
dy[:n] = dT
dy[n:] = dadt
return dy- Parameters
- _EXP_ARG_MAX = 500
- _FL20_A1 = 20000
- _FL20_A2 = 1e+07
- _FL20_DIFF = 1.38889e-07
- _FL20_DIM = 980
- _FL20_DX = 2.04499e-05
- _FL20_E1 = 60000
- _FL20_E2 = 80000
- _FL20_H_CONV = 10
- _FL20_INV_DX2 = 2.39121e+09
- _FL20_KAPPA = 0.2
- _FL20_M = 0.5
- _FL20_N_NODES = 490
- _FL20_N_ORD = 1.5
- _FL20_SRC = 291.666666667
- _R_GAS = 8.314
- _T_AMBIENT = 298
- _T_CEIL = 5000
- _T_FLOOR = 200
- Initial condition
- y(0) = [453.15, 298, 298, 298, 298, 298, …] [shape=(980,), min=0.001, max=453.15]
- Horizon
- t ∈ [0, 300]
Canonical RHS excerpt captured from the same registered callable used for the published benchmark. Frozen closure values are summarized below; helper imports and solver settings are intentionally omitted.
Fingerprint
Spread: extreme
Default noise: low
Recommendation snapshot
Clean best: SolvSRK
Noisy best: SciPy BDF
Coverage
14 solver arms · clean + 5 noise levels
Ranked on survival, precision, and speed
Versions & freeze
Methodology →- Freeze
- 2026-08-13
- libsolvsrk
- 2.3.0
- SciPy
- 1.14
- SUNDIALS
- CVODE (bundled backend)
20 seeds/cell default · 14 arms · TRL 4–5 · simulation-lab validated · this page: FP 20-Layer Unidirectional CF Layup (dim=980) (fp-20-layer-unidirectional-cf-layup-dim-980)
Governed SolvTune benchmark freeze; per-arm medians only. RHS definitions and raw trial rows are not published.
Self-reported by Resonix Labs · not independently verified
Results matrix
Pick an objective and a noise level to rank all arms on survival, median SCD, median nfev, and median wall time. Medians across seeds.
Objective
Best overall trade-off of survival, precision, and speed.
Noise level
| # | Solver | Survival | SCD | nfev | Wall | Score |
|---|---|---|---|---|---|---|
| 1 | SolvSRK | 100% | 11.3 | 30,080 | 13.72 s | 0.889 |
| 2 | SciPy RadauSciPy | 100% | 10.5 | 8,774 | 4.07 s | 0.870 |
| 3 | SciPy RK23SciPy | 100% | 9.4 | 476,126 | 13.57 s | 0.843 |
| 4 | CVODE BDFexternal | 100% | 7.9 | 7,277 | 1.19 s | 0.807 |
| 5 | SciPy LSODASciPy | 100% | 7.8 | 28,495 | 887 ms | 0.804 |
| 6 | SciPy BDFSciPy | 100% | 7.5 | 2,710 | 954 ms | 0.798 |
| 7 | CVODE Adamsexternal | 100% | 6.8 | 36,920 | 9.18 s | 0.781 |
| - | SciPy RK45SciPy | 0% | - | - | - | - |
| - | SciPy DOP853SciPy | 0% | - | - | - | - |
| - | Tsit5external | 0% | - | - | - | - |
At Clean, best balanced arm is SolvSRK.
Values are medians across seeds, measured by Resonix Labs on Resonix hardware and not independently verified; nfev and wall are on reference lab hardware (indicative). Under injected noise only SolvSRK and the SciPy arms are run. How we measure accuracy → · Verification status →
Cite this page
Replace the access date. Pin the freeze ID and library versions when comparing against a later export. Cite it as what it is - a self-reported vendor benchmark, not an independently verified result. The note field says so; please keep it.
@misc{resonix_evidence_fp_20_layer_unidirectional_cf_layup_dim_980_2026,
title = {Resonix Evidence Portal: FP 20-Layer Unidirectional CF Layup (dim=980)},
author = {{Resonix Labs (Canada) Inc.}},
year = {2026},
howpublished = {\url{https://resonixusa.com/evidence/problems/fp-20-layer-unidirectional-cf-layup-dim-980}},
note = {Self-reported vendor benchmark; internally generated by Resonix Labs and not independently verified. Accessed YYYY-MM-DD. Freeze 2026-08-13; libsolvsrk 2.3.0; SciPy 1.14.}
}Related
TRL 4–5 · simulation-lab validated · 398 problems · 14 solver arms · clean + 5 noise levels
Freeze: 2026-08-13 · scipy 1.14 · libsolvsrk 2.3.0 · Methodology
Self-reported by Resonix Labs · not independently verified · Verification status