MaximaLabs
All guides
Analytical ultracentrifuge, bioseparations laboratory

How to simulate differential sedimentation — stokes settling, size selectivity, svedberg

Four centrifuges run side by side on the same medium so the **d² law is visible as a result rather than asserted**: 20, 50 and 100 nm protein particles at 200,000 × g, plus a 10.24 nm particle at 250,000 × g. `centrifuge` is sized by **Sigma theory**, `Σ = Q/(2 v_g)`, where `v_g` is the g-amplified Stokes settling velocity `v = d²(ρ_p − ρ_m)g/(18η)`. So the reported `sigma_m2` carries the settling velocity, recoverable as `v = Q_liquid/(2Σ)`. Against the closed form the solver agrees to machine precision (relative difference 0 to 4e-16), which makes this a check of the tool and not just a demonstration of it. What it shows: **v(20 nm) = 1.53e-3 cm/s**; the 100 nm particle settles **exactly 4×** faster than the 50 nm one, since velocity goes as the square of diameter and every other term cancels; and the 250,000 × g branch is sized so its velocity is 5e-4 cm/s, giving a sedimentation coefficient **s = v/ω²r = 20.4 S** — the Svedberg range real proteins occupy (catalase 11.3 S, ribosome ~70 S). **Bounded, and the bound matters.** This is terminal Stokes velocity at constant field: no wall, no concentration gradient, no Boycott effect, and no hindered settling. Run the 50 and 100 nm particles for 30 minutes and the arithmetic says they separate by 51 cm, which no rotor can deliver — both pellet against the tube bottom first, the 100 nm one in about four minutes. The *ratio* is robust; the distance is what the formula says rather than what a centrifuge does. It is also the reference wiring for a solids flowsheet: the particles arrive as their own `phase: "solid"` feed with a flat `solids` payload, wired **directly** to the centrifuge. Declaring them `liquid` or `mixed`, or routing them through a mixer first, leaves the unit with no solid-phase inlet and the solve fails.

LIQ P21
SOL P21
solid
liquor
cake
filt
CF P21
PEL P21
SUP P21
LIQ P22a
SOL P22a
solid
liquor
cake
filt
CF P22a
PEL P22a
SUP P22a
LIQ P22b
SOL P22b
solid
liquor
cake
filt
CF P22b
PEL P22b
SUP P22b
LIQ P23
SOL P23
solid
liquor
cake
filt
CF P23
PEL P23
SUP P23
  1. 1
    Open the ready-made model

    Open the "Differential sedimentation — Stokes settling, size selectivity, Svedberg" model in the MaximaLabs workspace — no install, no license. It loads live on the canvas, ready to edit and run.

  2. 2
    Confirm the thermodynamics

    This process is modeled with the PENG-ROBINSON property package over water — already selected, so the phase equilibrium and enthalpy are physically consistent from the first run.

  3. 3
    Review the flowsheet

    The flowsheet chains 4× CF P23. Every block is a real, solvable unit op you can reconfigure on the canvas.

  4. 4
    Run the simulation

    Click Run. The deterministic solver converges the material and energy balances (recycles included) and fills the live stream table — the AI never invents a number.

  5. 5
    Read the results and iterate

    Inspect the converged streams, tweak a spec, and re-run — or ask the AI copilot to explain a result or diagnose a failed solve in plain English.

What you'll build
Thermodynamics
PENG-ROBINSON
Components
water
Unit operations
4× CF P23
Open this model in the workspace

Opens live on the canvas — free, no install.

Explore the model & flowsheet

Frequently asked questions

What does the Differential sedimentation — Stokes settling, size selectivity, Svedberg model simulate?
Four centrifuges run side by side on the same medium so the **d² law is visible as a result rather than asserted**: 20, 50 and 100 nm protein particles at 200,000 × g, plus a 10.24 nm particle at 250,000 × g. `centrifuge` is sized by **Sigma theory**, `Σ = Q/(2 v_g)`, where `v_g` is the g-amplified Stokes settling velocity `v = d²(ρ_p − ρ_m)g/(18η)`. So the reported `sigma_m2` carries the settling velocity, recoverable as `v = Q_liquid/(2Σ)`. Against the closed form the solver agrees to machine precision (relative difference 0 to 4e-16), which makes this a check of the tool and not just a demonstration of it. What it shows: **v(20 nm) = 1.53e-3 cm/s**; the 100 nm particle settles **exactly 4×** faster than the 50 nm one, since velocity goes as the square of diameter and every other term cancels; and the 250,000 × g branch is sized so its velocity is 5e-4 cm/s, giving a sedimentation coefficient **s = v/ω²r = 20.4 S** — the Svedberg range real proteins occupy (catalase 11.3 S, ribosome ~70 S). **Bounded, and the bound matters.** This is terminal Stokes velocity at constant field: no wall, no concentration gradient, no Boycott effect, and no hindered settling. Run the 50 and 100 nm particles for 30 minutes and the arithmetic says they separate by 51 cm, which no rotor can deliver — both pellet against the tube bottom first, the 100 nm one in about four minutes. The *ratio* is robust; the distance is what the formula says rather than what a centrifuge does. It is also the reference wiring for a solids flowsheet: the particles arrive as their own `phase: "solid"` feed with a flat `solids` payload, wired **directly** to the centrifuge. Declaring them `liquid` or `mixed`, or routing them through a mixer first, leaves the unit with no solid-phase inlet and the solve fails.
Which thermodynamic method does it use?
The PENG-ROBINSON property package, over water — already selected. You can switch the method on the canvas before running.
Which unit operations are in the flowsheet?
It chains 4× CF P23. Every block is a real, solvable unit operation you can reconfigure, add to, or remove.
Do I need to install software or buy a license?
No. Differential sedimentation — Stokes settling, size selectivity, Svedberg runs entirely in your browser on MaximaLabs — free, no install, no license. Open the model to load it live and run the deterministic solver.

More guides like this

Ludwigshafen, Germany

Polymerization reactor comparison: CSTR cascade vs. tubular (PFR)

The same monomer feed, split 50/50 into two trains at equal total reactor volume: a 3-stage CSTR cascade (each stage blending its own new polymer with the MWD already carried by the stream from the stage before it) versus a single tubular (PFR) reactor. Both converge to a real, different Mn/Mw/PDI — the PFR reaches higher conversion at the same total volume, but its cumulative polydispersity is measurably broader than the cascade's, because chain length drifts continuously with conversion along the tube while each CSTR stage makes its own local production under one fixed monomer concentration. See the 'Polymerization Reactors' docs page for the full reasoning.

Amarillo, Texas, USA

Helium liquefaction (Linde-Hampson cycle)

A Joule-Thomson liquefaction loop: helium is compressed, aftercooled, precooled well below its ~40-51 K JT-inversion temperature, then throttled through a valve into a flash drum that draws off liquid helium while the cold vapor recycles back to the compressor suction — the same JT-cycle physics as the platform's LNG/refrigeration loops, run on a fluid whose inversion curve makes precooling mandatory before the throttle can cool it at all.

Brevik, Norway (Heidelberg Materials Norcem — first full-scale cement CCS)

Cement kiln calcination + CO2 liquefaction

Preheated limestone (CaCO₃) enters an indirectly-heated rotary kiln that calcines it toward CaO, releasing a CO₂-rich off-gas; the gas passes a knockout drum, then is compressed and cooled to liquefy the CO₂ for transport/storage.

Mont Belvieu NGL fractionation complex, Texas, USA

NGL fractionation: single-shell Petlyuk dividing-wall column

A natural-gas-liquids cut (ethane / propane / n-butane) split into three on-spec products by ONE thermally-coupled dividing-wall column instead of two columns in series. This uses the native Petlyuk unit op: a prefractionator whose reflux and boilup are supplied by the main column (the single condenser + single reboiler), solved by an outer Wegstein loop that converges the bidirectional vapor/liquid coupling to self-consistency — the coupling legacy sequential-modular solvers approximate with two connected columns whose recycle fails to converge.

Pekin, Illinois, USA

Ethanol–water distillation

An 8-stage column concentrating ethanol overhead toward the azeotrope (the headline demo).

Tray efficiency — real trays vs ideal stages

The same ethanol–water column solved with a Murphree vapor tray efficiency of 0.7 instead of ideal equilibrium stages. A real sieve/valve tray never reaches full vapor-liquid equilibrium — the vapor leaving it only partly approaches the equilibrium composition with the tray liquid, mixing in un-equilibrated vapor from the tray below: y = E·K·x + (1−E)·y_below (Murphree 1925). At E = 0.7 each of these 12 trays does 70% of an ideal stage's work, so the overhead ethanol is lower than an equilibrium column of the same tray count would predict — which is exactly why a real column needs more trays than a shortcut (ideal-stage) calculation says. Both HYSYS and Aspen RadFrac expose this per-tray efficiency; set murphree_efficiency back to 1.0 to recover the ideal-stage column. The efficiency auto-selects the component-flow Naphtali-Sandholm solver (the reduced-form solvers carry no explicit per-tray VLE row to apply an efficiency to).

Stop fighting legacy software. Build your first flowsheet in 60 seconds.