The Forcing Functions

Three north stars.

FrankenSim is driven by three flagship pipelines plus a marquee demo — each a forcing function that must work end-to-end, each producing a certified artifact, not just a number.

The FrankenSim continuum: geometry, physics, optimization, and rendering fused into one memory-safe Rust kernel.
Flagship 01

Ornithoid Aircraft.

A certified Pareto atlas for a bird-like flyer

A multi-inlet, bird-like aircraft optimized jointly for lift-to-drag, stability, and maneuverability across a flight envelope.

Objective
maximize L/D × stability × maneuverability
Methods
BEM + FMM + Kuttavortex-particle wakesDirac-structure couplingSE(3) integratorsKoopman surrogates
The Payoff

delivers a certified Pareto atlas with SOS Lyapunov region-of-attraction proofs.

Delivered by P5 · Aero Stack (Weeks 44–56)Planned

Adjoint-native gradients

differentiate through the solution — not through the iterations

implicit function theorem
design paramsp ∈ ℝ¹²⁸p₀p₁p₂p₃p₄p₁₂₇assembleK(p), f(p)forward solveK·u = fobjectiveJ(u, p)dJ/dp = ∂J/∂p − λᵀ·∂R/∂p — one linear solve, gradient for every pcost to get ∇J over 128 paramsadjoint≈ 2 solvesfinite-diff128 solves
gradient checks gate every merge
adjoint-native gradients · press Play

The adjoint aero solve: solve the flow forward once, then a single adjoint solve returns the full gradient over every shape parameter — differentiating through the converged solution, not the solver's iterations.

Flagship 02

Seismic-Minimal Frame.

Least material, provable fragility

A building frame that uses the minimum material while carrying a certified seismic fragility curve — the anti-paperclip constraint made real.

Objective
minimize material subject to a certified fragility bound
Methods
IGA + Kirchhoff–Love shellsfiber beamsground-structure PDHGKanai–Tajimi + MLMCe-stopping
The Payoff

a fragility curve with anytime-valid stopping — you stop the moment the evidence is decisive.

Delivered by P4 · Structures at Scale (Weeks 34–44)Planned
max e 1.0 · gen 0 · idle
E-VALUE WEALTH · anytime-valid · peek continuously, stop anytime12510203004812generation / evidence stepe-value (log)1/α = 20 · decisive
design A 1.0design B 1.0design C 1.0design D 1.0design E 1.0
e-processes are anytime-valid — no p-hacking penalty for peeking; the first curve to reach 20 wins

Anytime-valid fragility stopping: candidate frames race under a betting e-process, and the Kanai–Tajimi + MLMC campaign halts the instant the seismic evidence is decisive — no fixed-horizon core-hours burned.

Budgets first — Error & Time Ledgers

every op takes an accuracy / time / memory budget; they compose end-to-end

within budget
TIME LEDGERseconds → hours, summing toward the wall budgetwall 2.0h2.6h full-scaleERROR LEDGERrel-error contributions, summing toward the targettarget 2.0e-24.5e-2 full-scalegeom0.2h · 0.3e-2mesh0.3h · 0.6e-2assemble0.3h · 0.3e-2solve0.5h · 0.3e-2adjoint0.1h · 0.3e-2
feasible ✓ · wall 1.4h / 2.0h · rel-err 1.8e-2 / 2.0e-2baseline plan fits both budgets

The Error Ledger and Time Ledger attribute every digit of the fragility bound and every core-second back to the operator that produced it, so a tighter tolerance has a visible price.

Flagship 03

The Spout That Never Dribbles.

A laminar-pour vessel, rendered from the same bytes

A pouring vessel shaped so its stream stays laminar — an Orr–Sommerfeld stability objective validated against a free-surface LBM pour.

Objective
maximize pour stability (minimize spectral growth)
Methods
Orr–Sommerfeld stabilityfree-surface LBMCarreau rheologyspectral rendering
The Payoff

the marketing shot and the physics are the same bytes — a differentiable render of the certified design.

Delivered by P3 · Fluids I (Weeks 24–34)Planned
active 42 · cut 38 · meshing steps: 0
φ < 0φ > 0 · inactivezero body-fitted mesh · outside 208 cells inactive · hover a cut cell
interior (active)cut cell · ghost + Nitscheoutside (inactive)quadrature point

The vessel is a level set — an SDF — so physics runs directly on it with zero body-fitted meshing (shown: cut cells with a ghost penalty and Nitsche BCs). The laminar pour itself is a free-surface LBM solve tuned to an Orr–Sommerfeld stability objective.

The marketing shot is the physics.

For the Spout That Never Dribbles, the hero render is not a separate art asset. Lumen's differentiable renderer draws the certified vessel from the same bytes the Orr–Sommerfeld stability objective and the free-surface LBM pour were computed on. The image you would put on the box is a differentiable function of the design — so you can optimize what the product looks like and what it does at once, and never ship a picture the physics disagrees with.

The Marquee · P2

Topology optimization on a raw SDF — no mesh in the loop.

The forcing function for the whole geometry-physics bridge. A SIMP density field evolves to minimize compliance under a volume fraction, its physics computed by CutFEM directly on the level set. A grey blob resolves into a classic cantilever truss — and every iterate carries a composed error certificate, with the mesh-step counter pinned at zero.

iter 0 · compliance 100.0 ↓ · vol 0.40 · certified ✓ · mesh steps: 0
DESIGN DOMAIN · SIMP density on CutFEM-on-SDFcompliance ↓fixed supportFvolume fraction — held at V/V₀ = 0.40every iterate carries a composed error certificate · body-fitted mesh steps: 0iter 0/40
void 0 → solid 1fixed supportapplied load Fvolume held at 0.40load-carrying members emerge; the mesh never does

The optimizer reshapes the boundary every iteration; compliance falls as the density field condenses into load-bearing structure.

active 42 · cut 38 · meshing steps: 0
φ < 0φ > 0 · inactivezero body-fitted mesh · outside 208 cells inactive · hover a cut cell
interior (active)cut cell · ghost + Nitscheoutside (inactive)quadrature point

Underneath, CutFEM-on-SDF supplies FEM-grade physics on the moving level set — the cut cells are certified, so the optimizer never has to wait on a remesh.

Scope

FEEC elasticity, CutFEM-on-SDF, matrix-free p-MG + AMG, adjoints, SIMP.

Exit Gate

Marquee demo: topology optimization on a raw SDF — no mesh in the loop — with a composed error certificate.

From forcing function to kernel.

Every flagship bottoms out in the same seven-layer continuum. See how the kernel is built, or clone it and run the vertical skeleton yourself.