FEM methods & multiphysics · 108

Shared FEM machinery: H1/H(div) elements, transfers, PML, solvers, and coupled multiphysics contracts. Method shown on hover; tolerance is the pass gate. back to QA

CaseTierReferenceGate
fem.h1-curved-p3-transfer.v1physics-gatesame-mesh complete global-state identitymaximum absolute global-DOF transfer error <= 2e-10
fem.adaptive-lineage-transfer.v1physics-gatelinear H1 and H(curl) fields, zero RT0/RT1 states, stale-history rejectionH1 error <= 2e-12; H(curl) error <= 2e-10; two typed rejections
fem.adaptive-lineage-chain.v1physics-gatelinear H1/H(curl) fields and finite zero RT0/RT1 statesH1 error <= 2e-12; H(curl) error <= 2e-10; two typed rejections
fem.adaptive-mpc-transfer.v1physics-gateconstant H1 state, zero vector states, complex-phase and constraint-residual rejectionH1 error <= 2e-10; at least two typed MPC rejections
fem.spatial-source-locator.v1physics-gateanalytic linear scalar transfer, design transfer finiteness, outside/non-finite rejectionscalar transfer error <= 2e-12; at least two typed rejection checks
fem.complex-dirichlet-elimination.v1physics-gateindependent dense complex system with constant prescribed boundary stateboundary and interior errors <= 1e-12; at least two rejection checks
fem.cross-space-assembly.v1physics-gateintegral of each affine P1 basis against an element-local constant statelocal mass error <= 2e-12; at least one typed rejection
fem.hdiv-curved-cross-space-assembly.v1physics-gatefinite/nonzero curved RT mass, swapped-order transpose parity, malformed-topology rejectioncurved transpose error <= 2e-10; nonzero vector block; at least one typed rejection
fem.hdiv-divergence-cross-space-assembly.v1physics-gatefinite/nonzero physical divergence block, swapped-order transpose parity, malformed-topology rejectiontranspose error <= 1e-12; nonzero physical block; at least one typed rejection
fem.hdiv-curved-divergence-cross-space-assembly.v1physics-gateindependent published Piola reference, finite/nonzero curved blocks, swapped-order transpose parity, malformed-topology rejectiontranspose error <= 2e-10 in both dimensions; nonzero physical blocks; typed rejection in both dimensions
fem.mixed-full-assembly.v1physics-gateindependent four weighted mass forms with analytic rectangular block dimensionsblock/monolithic apply error <= 2e-12; lower-right block nonempty; at least one typed rejection
fem.h1-mixed-order-transfer.v1physics-gateanalytic quadratic and cubic polynomial fieldsboth transfer errors <= 2e-12; at least two typed rejections
fem.h1-curved-p2-transfer.v1physics-gatesame-mesh complete global-state identity with curved edge midpointcurved transfer error <= 2e-11; at least two typed rejections
fem.cross-source-vector-transfer.v1physics-gateconstant H(curl) parity, finite RT0/RT1 transfer, typed outside-domain rejectionH(curl) error <= 2e-10; at least one typed rejection
fem.hdiv-rt1-refinement-transfer.v1physics-gatefinite multi-element refinement, 3-D identity parity, typed outside/order rejection3-D identity error <= 2e-8; two typed rejections
fem.pml-open-boundary-reference.v1physics-gateanalytically continued outgoing plane wave with independent attenuation integralphysical error <= 2e-3; full error <= 3e-3; reflection ratio <= 3e-3; residual <= 1e-11
fem.pml-directional-contract.v1physics-gateindependent complex coordinate-stretch coefficient identities on all 2-D sides and XY corner28 finite coefficient/direction checks pass to 1e-12
fem.validation-matrix.v1smokeexplicit maturity/reference/gate recordsexact
fem.element-contracts.v1smokeexact nodal partition/unisolvence, lowest-order DOF counts, edge traces, orientation parity, constant-curl property, and independent curl(grad) finite-difference checkexact counts; 2e-13 identity residual; finite-difference curl(grad) <= 5e-8
fem.dg-p3-contract.v1physics-gateindependent P3 nodal partition, contiguous DG numbering, expected face traces, and cell-bubble exclusionelement_checks >= 4; boundary_trace_checks >= 4; partition residual <= 2e-12
fem.h1-p3-curved-geometry-reference.v1physics-gateindependent cubic nodal coordinates, partition of unity, positive Jacobian, finite gradients, and symmetrypartition residual <= 2e-12; symmetry residual <= 2e-10; positive finite Jacobian
fem.operator-invariants.v1physics-gateconstant-field diffusion conservation, transpose symmetry, unit-domain mass, positive quadratic energy5e-13 residuals; mass total = 1 within 5e-13; energy > 0
fem.high-order-orientation-patch.v1smokeanalytic quadratic polynomial reproduction and canonical Nedelec edge parityP2 value/gradient residuals <= 2e-12; Nedelec parity residual <= 2e-13
fem.p3-tetra-trace-orientation.v1physics-gateindependent cubic polynomial and shared-face continuity referencetwo reference checks pass below the focused trace threshold
fem.p3-tetra-reproduction.v1smokeanalytic cubic polynomial on affine tetrahedral elementscubic interpolation residual <= 5e-12
fem.p3-h1-triangle-quad-reproduction.v1physics-gateindependent cubic polynomial on affine P3 triangles and quadrilateralscubic interpolation residual <= 2e-11
fem.vector-laplacian-convergence.v1physics-gateindependent sine manufactured vector field with analytic vector-Laplacian sourceL2 refinement rate >= 1.5 and H1-seminorm refinement rate >= 0.75
fem.waveholtz-direct-reference.v1physics-gateseparate assembled direct Helmholtz solvetemporal order >= 1.8 and all relative reference errors <= 1e-3
fem.nonlinear-tangent.v1physics-gatequadratic residual Taylor remainder under deterministic state/directionorder >= 1.8 across four steps
fem.mesh-transfer-adjoint.v1physics-gateindependent weighted primal/adjoint pairing identity1e-11
fem.mesh-transfer-design.v1physics-gatesource element design values and analytic area-weighted integral1e-11
fem.helmholtz-plane-wave-pollution.v1physics-gateindependent exact oblique complex plane wavefocused test convergence and finite residual/phase gates pass
fem.helmholtz-mms-convergence.v1physics-gateclosed-form plane wave, MMS, and scattering fieldsnonzero free unknowns on every leg; quadrature L2 <= 3e-3; far-field solved L2 <= 4e-2; decay mismatch <= 0.3; convergence mean rate >= 1.5
fem.biharmonic-3d.v1contract-gateanalytical mixed block structureexact block structure
fem.electrochemistry-3d.v1contract-gatepartition of unity and prescribed-coordinate bilinear energy identities2e-10
fem.semiconductor-3d.v1contract-gatepartition of unity and prescribed-coordinate bilinear energy identities2e-10
fem.thermoelectric-3d.v1contract-gatepartition of unity and prescribed-coordinate bilinear energy identities2e-10
fem.biharmonic-mms-3d.v1physics-gateu=sin(pi*x)sin(pi*y)sin(pi*z), q=3*pi^2*u, f=9*pi^4*uu <= 0.02; q <= 0.01; both u rates >= 1.7; both q rates >= 1.7; residual <= 1e-10
fem.aeroacoustics-drp-reference.v1reference-onlyanalytical characteristic phase/group dispersion and Fourier-mode propagationconvergence rate >= 3.5; Fourier-mode relative L2 <= 2e-5
coupling.conservative-power-contracts.v1smokeconstant preservation, weighted adjoint identity, analytic power accounting, typed invalid-input rejectiontrace and power residuals <= 1e-12; two invalid inputs rejected
coupling.body-fitted-monolithic.v1physics-gatefinite driven states, two-driver residual blocks, and boundary-spanning brace patch areaall residual blocks <= 1e-7; body-fitted mesh and brace patch non-empty
coupling.body-fitted-ported-monolithic.v1physics-gateproduction body-fitted FEM/BEM monolithic port continuity and driver blocksall residual blocks and port pressure jump <= 1e-7 across 60/90/120 Hz; finite non-empty port flow
coupling.scheduler-preconditioner-contracts.v1smokescheduler quota/refund counters, CPU completion bridge, saddle-point identity application, typed rejection pathsexact counters and rejection counts; identity error <= 1e-12
fem.schwarz-pml-reflection.v1physics-gatedirect sparse Helmholtz solve on the identical mesh and RHSrelative error < 1e-3 for target reflection 1e-2, 1e-4, and 1e-6
fem.hdiv-rt0-reference.v1physics-gateindependent affine RT0 moments, Piola metric, and unit-flux divergence identitymass <= 2e-11; divergence/orientation <= 2e-10; mixed conservation <= 2e-11
fem.hdiv-rt0-3d-reference.v1physics-gateindependent tetrahedron/simplex and tensor-product hexahedron RT0 moment integralsmass <= 5e-11; divergence/permutation <= 5e-10; mixed conservation <= 5e-11
fem.hdiv-rt1-triangle-reference.v1physics-gateindependent edge/cell moment quadrature, finite-difference divergence, topology, orientation, and admission checksfive independent RT1 reference checks pass
fem.hdiv-rt1-quadrilateral-reference.v1physics-gateindependent edge/cell moment quadrature, finite-difference divergence, topology, and shared-edge orientationfour independent RT1 quadrilateral reference checks pass
fem.hdiv-rt1-tetra-reference.v1physics-gateindependent tetrahedral face/cell quadrature, finite-difference divergence, topology and Piola checksthree independent RT1 tetrahedron reference checks pass
fem.hdiv-rt1-hexa-reference.v1physics-gateindependent tensor-product face/cell quadrature, finite-difference divergence, topology Piola checksthree independent RT1 hexahedron reference checks pass
fem.hdiv-curved-p2-reference.v1physics-gateindependent P2 isoparametric physical-coordinate reconstruction and physical-divergence finite differencetwo independent curved H(div) reference checks pass
fem.piola-consumer-contract.v1physics-gateproduction checked tabulation, topology-owned shared-edge orientation, and typed singular/non-finite Jacobian rejectionall accepted values finite; exact typed rejection variants; minimum checks 4/2/4
fem.natural-boundary-contract.v1physics-gateindependent unit-square/unit-cube boundary resultants and P1 Robin mass/RHS integralsresultants/integrals <= 1e-12; two typed rejection paths leave systems unchanged
fem.mpc-system-contract.v1physics-gateindependent P^H A P/P^H b fixtures, recovered periodic residual, and ordering/error checksanalytic values <= 1e-12; recovery residual <= 1e-14; three typed rejection checks
fem.hanging-node-contract.v1physics-gateindependent affine edge parameters, conforming-mesh no-op, typed malformed/order/tolerance rejectiongeometric admission tolerance 1e-10; interpolation errors <= 1e-14; four typed rejection checks
fem.higher-order-hanging-node-contract.v1physics-gateindependent coarse-edge/triangular-face trace partition, polynomial recovery, typed rejectiongeometric admission tolerance 1e-10; interpolation errors <= 1e-10; five typed rejection checks
fem.higher-order-natural-boundary-contract.v1physics-gateindependent unit-square/unit-cube constant resultants, P2 Robin mass/RHS, P2/P3 curved trace finitenessresultants/integrals <= 1e-10; finite geometry; transactional rollback
fem.hdiv-rt1-mixed-conservation.v1physics-gateindependent quadrature of RT1 Piola divergence against constant P1/Q1 fieldfour domain checks constant-field conservation contract pass
fem.solver-failure-contract.v1smokeshared SolverStatus and AccuracyStatus mapping for FEM solver errorsall four status/accuracy mappings match the shared contract
fem.prebaked-checked-contract.v1physics-gateindependent malformed-triplet admission and valid CSR merge contractvalid merge succeeds; invalid triplets reject
fem.p2-poisson-convergence-3d.v1physics-gateindependent quartic manufactured Poisson solution and refinement ratesfinite residuals and expected L2/H1 refinement rates
fem.curved-p2-h1-geometry.v1physics-gateindependent P2 nodal interpolation and reference-coordinate gradient finite differencepositive Jacobian and finite-difference Hessian error
fem.mixed-hdiv-conservation.v1physics-gateindependently integrated RT0 divergence and P1 partition of unitymaximum conservation error <= 2e-12 at refinements 1, 2, 4
fem.curved-p2-assembly-reference.v1physics-gateindependent curved-P2 stiffness symmetry conservation energy checkssymmetry <= 2e-12; constant-field row sums <= 2e-11; positive energy
fem.hdiv-curved-rt1-mixed-physical-reference.v1physics-gateindependent Q2 physical map and RT1/P1 mixed divergence quadratureall 288 checks pass; max block error <= 2e-10; minimum Jacobian determinant >= 0.05
fem.vector-natural-boundary-contract.v1physics-gateindependent unit-square/unit-cube constant traction resultants, midpoint loading, and rollbackconstant resultants <= 1e-10; finite callbacks; transactional rollback
fem.higher-order-hexa-hanging-node-contract.v1physics-gateindependent four-patch coarse-face partition, constant reproduction, polynomial recovery, typed rejectiongeometric admission tolerance 1e-10; interpolation and recovery errors <= 1e-9
fem.vector-h1-hanging-node-contract.v1physics-gateindependent scalar/vector constraint-count multiplication, polynomial recovery, zero-component rejectiongeometric admission tolerance 1e-10; imaginary recovery error <= 1e-14
fem.hdiv-rt1-hanging-edge-contract.v1physics-gateindependent RT1 edge traces, signed orientation, MPC recovery, topology/order rejectiongeometric admission tolerance 1e-10; coefficient/recovery error <= 1e-12
fem.hdiv-rt0-hanging-edge-contract.v1physics-gateindependent positive/negative half-edge checks, MPC recovery, and unsupported-input rejectiongeometric admission tolerance 1e-10; signed-fraction errors <= 1e-12; imaginary recovery error <= 1e-14
fem.hdiv-rt1-hanging-tri-face-contract.v1physics-gateindependent owner-face quadrature, signed subface fractions, recovery, unsupported-input rejectiongeometric admission tolerance 1e-10; coefficient/recovery errors <= 1e-11
fem.hdiv-rt0-hanging-tri-face-contract.v1physics-gateindependent one-third subface checks, conforming no-op, recovery, and unsupported-input rejectiongeometric admission tolerance 1e-10; signed-fraction errors <= 1e-12; imaginary recovery error <= 1e-14
fem.hdiv-rt0-periodic-face-contract.v1physics-gatecyclic/reversed face orientation, complex MPC recovery, typed rejectionphase/recovery residual <= 1e-12
fem.hdiv-rt1-periodic-face-contract.v1physics-gatecyclic/reversed face mode change, complex MPC recovery, typed rejectionphase/recovery residual <= 1e-12
fem.hdiv-rt0-refinement-transfer.v1physics-gateindependent affine RT0 field edge moments, 3-D identity recovery, unsupported-input rejection2-D error <= 2e-10; 3-D identity error <= 2e-10; at least 4 rejection checks
fem.hdiv-rt1-curved-p2-transfer.v1physics-gatecurved self-recovery, affine-to-curved finite/complete transfer, typed rejectioncurved identity <= 2e-8; affine-to-curved recovery checks >= 2
fem.hdiv-rt1-periodic-quad-face-contract.v1physics-gateidentity, quarter-turn, reflected square maps, complex MPC recovery, typed rejectionphase/recovery residual <= 1e-12
fem.hdiv-rt1-periodic-mixed-face-contract.v1physics-gatemixed pairing, seven-mode MPC recovery, 7x7 condensation, typed rejectionphase/recovery residual <= 1e-12
fem.hdiv-rt1-multielement-ownership.v1physics-gateunit-cube tetrahedral generator, shared-face reuse, checked topology layoutexact integer topology and numbering checks
fem.eigenphysics-wrappers.v1physics-gatebeam theory, PEC cutoff, and Rayleigh-quotient consistencyEM wavenumber error <= 5e-3; residuals <= 1e-6; all legs converged
fem.biharmonic-validation.v1physics-gatemanufactured field and Navier 0.00406 simply-supported coefficientfine u <= 10%; fine q <= 15%; rates >= 1.0/0.7
fem.electroacoustics-validation.v1physics-gateindependent decoupled block solve, physical scaled reciprocity under swapped forward drives, and power balancereciprocity and power imbalance <= 1e-9
fem.electrostatics-point-charge-image.v1physics-gateindependent Coulomb source minus mirror charge, potential and electric fieldpotential <= 1e-2; field <= 1e-1; rates >= 1 and 0.8; residual/boundary <= 1e-9; source <= 1e-12
fem.electrostatics-two-cylinder.v1physics-gatebipolar potential and field; full capacitance pi epsilon/acosh(d/(2r)); explicit finite-sector truncationfull capacitance <= 2e-2; finite-domain capacitance <= 3e-3; flux <= 4e-2; field <= 6e-2; rates >= 0.7; conservation/residual <= 1e-9
fem.piezoelectric-1d-axial-load.v1physics-gateindependent stress-charge constitutive relations and uniform loaded/driven rod fields; cD=cE+e^2/epsilon for open electrodesall relative errors/residuals and symmetry <= 1e-9; 24 actually solved cases
fem.piezoelectric-1d-thickness-mode.v1physics-gatef=sqrt((cE+e^2/epsilon)/rho)/(2L), u=cos(pi x/L), phi=e/epsilon*(u-u(0)); coupling-off frequency ratiofrequency <= 5e-4; field <= 2e-3; second-order rates >= 1.8; residual/potential/stiffening <= 1e-9
fem.thermoelectric-3d-mms.v1physics-gateindependent strong equations with sigma=2, alpha=0.005, Pi=1.5, k=1.5; Lap(V)=-1.2, Lap(theta)=-80; heat source corrected by sigma|grad(V)|^2field rates >= 1.8; gradient rates >= 0.8; P2 field/gradient errors <= 1e-9; original residual <= 1e-9; missing-term controls require nonzero accuracy errors
fem.acoustics-driven-sweep.v1physics-gateindependent exact one-dimensional Helmholtz pressure; fresh stiffness and mass assembly and direct solve at each sweep frequencypressure 100/1Hz <=2%/1%; GMRES/direct <=1e-6; residual <=1e-8; sweep solution/matrix <=1e-9/1e-12; three swept solves
fem.electroacoustics-cavity-bar.v1physics-gateindependent closed-form damped waves and interface impedance sum; weighted physical reciprocity and fluid/solid dissipationpressure/displacement errors <=0.003 and rates >=1.8; residual/reciprocity/power <=1e-9; omitted coupling error >=0.9
fem.semiconductor-3d-coupled-mms.v1physics-gateindependent quadratic potential/electron/hole strong-form sources, Einstein diffusivities and signed drift; steady native E.grad(c) modelfield/gradient rates >=1.8/0.9; fine errors <=0.005/0.1; P2/residual/charge balance <=1e-9; five omitted-term errors >=1e-3
fem.electrochemistry-3d-pnp-electrode.v1physics-gateindependent mixed-xyz quadratic fields and strong forms; inverse symmetric BV prescribed metal-voltage field, solved concentration activity, electric/ionic fluxes and electrode currentfield/current spatial rates >=1.8; flux rates >=0.9; fine current <=0.01; P2/residual/weak balance <=1e-9; four omitted-term controls >=1e-3
fem.autofrettage-spatial.v1physics-gateincompressible elastic-plastic interface solution and Lame elastic unloadingfine residual stress <= 1%; history stress <= 2%; convergence rate >= 0.7; equilibrium <= 1e-9
fem.axisymmetric-cylinder-harvest.v1physics-gategenerated independent OpenCourant QUAD N2D3D=1 pressure-hold history and Lame fieldsengine displacement/energy <= 0.5%; Lame displacement <= 0.1%; refinement rate >= 1.5; equilibrium <= 1e-8
fem.necking-bar-spatial.v1physics-gateindependent Voce uniaxial reaction turning point and fine integration of the nonuniform equilibrium branchpeak extension <= 1%; peak force <= 0.1%; displacement/reaction <= 0.1%; refinement rate >= 1.5; free equilibrium <= 1e-9
fem.foam-laws-integrity.v1contract-gateoriginal LOIS engine golden seriesfailed_probes = 0; numerical_solve_executed = 0
fem.foam-law51-material-replay.v1physics-gateoriginal engine stress/density/energy histories and independent exact PG adiabatic energystress <= 2e-4; density <= 2e-5; native exact energy <= 2e-4; original golden energy band usage <= 1; frozen-energy control >= 0.05
fem.foam-law33-material-replay.v1physics-gateoriginal independently generated and harvested OpenCourant single-brick historystress/work <= 0.5%; density/geometry <= 1e-5; air-disabled stress gap >= 1%; plastic strain >= 1%
fem.foam-law76-compression-replay.v1physics-gateoriginal independently generated and harvested OpenCourant brick explicitly selecting LAW76stress <= 0.2%; plastic strain/work <= 0.5%; density/geometry <= 1e-5; pressure-blind stress gap >= 1%; plastic strain >= 0.0005
fem.foam-law76-tension-replay.v1physics-gateoriginal independently generated and harvested OpenCourant brick explicitly selecting LAW76stress <= 0.2%; plastic strain/work <= 0.5%; density/geometry <= 1e-5; pressure-blind stress gap >= 1%; plastic strain >= 0.0005
fem.foam-law33-mixed-spatial.v1physics-gateindependent continuum stress/reaction, closed-form compressive inverse, integrated displacement and external workfine displacement/work <=0.5%; spatial rate >=1.5; mixed finest-pair rate reported without bound; separate smooth case gates finest-pair >=1.5; J error <=1e-7; load refinement <=2e-4; free force <=1e-
fem.foam-law33-spatial.v1physics-gateindependent continuum stress/reaction, closed-form compressive inverse, integrated displacement and external workfine displacement/work <=0.5%; spatial rate >=1.5; finest-pair rate >=1.5; J error <=1e-7; load refinement <=2e-4; free force <=1e-9; reaction <=1e-8; stress <=1e-3; elastic gap >=5%
fem.foam-law76-spatial.v1physics-gateindependent continuum stress/reaction, closed-form compressive inverse, integrated displacement and external workfine displacement/work <=0.5%; spatial rate >=1.5; finest-pair rate >=1.5; J error <=1e-7; load refinement <=2e-4; free force <=1e-9; reaction <=1e-8; stress <=1e-3; elastic gap >=5%
fem.crash-box-impact-harvest.v1physics-gateindependently harvested original CRASH box1_42 energy, impulse and displacement historiesnine history errors <= 1%; timestep changes <= 0.3%; energy residual <= 1e-6; elastic response change >= 50%; wall response change >= 10%
fem.opencourant.hertz-sphere.native-parity.v1physics-gateOpenCourant external engine health is recorded separately from native parity; OPENCOURANT_EXE or legacy RADISS_EXE selects the referenceall three full native/reference histories and final same-time probes <=3%; half-CFL changes <=0.3%; energy drift <=1%; final wall-force balance <=1%; original stiffness window 1e8..3e9 N/m; exact nati
fem.opencourant.taylor-bar.native-parity.v1physics-gateOpenCourant external engine health is recorded separately from native parity; OPENCOURANT_EXE or legacy RADISS_EXE selects the referencefull displacement/plastic fields and peak radius/length/plastic strain <=3%; half-CFL fields and peaks <=0.3%; exact native80us stop, complete external history, no erosion, numerical energy drift <=1%
fem.opencourant.rail-crush.native-parity.v1physics-gateOpenCourant external engine health is recorded separately from native parity; OPENCOURANT_EXE or legacy RADISS_EXE selects the referenceexternal analytical probe budgets remain case-specific; native_parity_executed must equal 1 before release success; unavailable optional engine yields structured SKIP, SONIUM_QA_REQUIRE_OPENCOURANT=1
fem.opencourant.bending-3pt.native-parity.v1physics-gateOpenCourant external engine health is recorded separately from native parity; OPENCOURANT_EXE or legacy RADISS_EXE selects the referencenative/OpenCourant displacement <=3%; original Euler-Bernoulli continuum approximation budget 10%; spatial changes <=2% then <=0.5%, contraction <=0.6; timestep change <=0.3%; every-step independent d
fem.opencourant.necking-bar.native-parity.v1physics-gateOpenCourant external engine health is recorded separately from native parity; OPENCOURANT_EXE or legacy RADISS_EXE selects the referenceexternal analytical probe budgets remain case-specific; native_parity_executed must equal 1 before release success; unavailable optional engine yields structured SKIP, SONIUM_QA_REQUIRE_OPENCOURANT=1