Crack Tip Fracture Simulator

A crack runs through a plate while the stress at its tip spikes toward infinity. The same run is drawn twice: the raw stress, and a view that stays bounded.

Loading fracture simulator

LIVEILLUSTRATIVE BROWSER DEMO

One damage run near a 1/sqrt(r) stress singularity, drawn as the stress and as its bounded r^0.5-weighted view

Live

Measured in Browser

1/sqrt(r)

Singularity Type

k = -0.5

Diagnostic weight r^0.5 (matches the exponent)

None

Speedup Claim

The Physics

In fracture mechanics, the stress field near a crack tip follows:

sigma = K / sqrt(2*pi*r)

As distance r approaches 0, stress approaches infinity. This is the famous stress singularity that makes crack simulation challenging.

The Bounded Diagnostic

The k=-0.5 view weights the stress by r^0.5; its log-slope is bounded on r > 0:

d/dr[1/sqrt(r)] = -1/(2*r^1.5) (diverges)
L_BG[1/sqrt(r)] = -1/2; D_BG = e^-1/2 (bounded diagnostic)

The stress field itself remains singular. The weighting changes what is plotted, not the damage, the timestep or the stall rule: both surfaces come from the same solver run.

Why This Matters

Engineering Context

Predicting when cracks will propagate matters for aircraft, bridges and pressure vessels. Fracture mechanics already works with the bounded quantity: the stress intensity factor K = lim sigma * sqrt(2*pi*r).

Real-Time Simulation

This browser demo reports the live solver dt, its local error and the peak of each view. It makes no speedup claim and does not reuse the CASCADE benchmark numbers.

Browser Check

This simplified damage model runs entirely in your browser. It is an illustration, not a validated fracture solver.

Related CASCADE Results

This demo is inspired by the pyfrac_crack_tip case in the local synthetic CASCADE benchmark (a toy distance-to-singular-point metric, not citable). The browser panel above is a separate illustration, not a reproduction of that benchmark.See full results →