FP 20-Layer Unidirectional CF Layup (dim=980)

PARITYS2 · dim 980

No 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.

Materials & composites

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

#SolverSurvivalSCDnfevWallScore
1SolvSRK
100%
11.330,08013.72 s0.889
2SciPy RadauSciPy
100%
10.58,7744.07 s0.870
3SciPy RK23SciPy
100%
9.4476,12613.57 s0.843
4CVODE BDFexternal
100%
7.97,2771.19 s0.807
5SciPy LSODASciPy
100%
7.828,495887 ms0.804
6SciPy BDFSciPy
100%
7.52,710954 ms0.798
7CVODE Adamsexternal
100%
6.836,9209.18 s0.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 →

SolvScout · free

Profile your problem for free

This page shows one published benchmark cell. SolvScout fingerprints your ODE, compares it to the full corpus, and recommends a solver with the same survival / precision / speed ranking you see here — including when a SciPy arm wins.

SolvSRK · 30-day trial

Run the winner on your machine

SolvSRK is the stiffness-adaptive integrator behind the SolvSRK column in these tables. Create an account, activate a machine, and take a 30-day trial — same binary you'd ship after purchase.

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://resonix.tech/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