1 · Spiral-wound element
2 · Solvent & solute
3 · CentrivaLab cartridge
estimated collection time only — preparation and stabilisation not included
expected flux — LMH
wall concentration cm — g/L in both
expected cp — g/L
Test programme
—
| Run | target cm (g/L) | ΔP (bar) | rpm | feed C0 (g/L) | Jv (LMH) | Vp (mL) | time (min) | expected cp (g/L) | Robs in cartridge |
|---|
Element rejection = 1 − cp / C0,element — divide the measured permeate concentration by the element feed, not by the cartridge feed.
Inside the element
| Position | u (m/s) | Re | k (LMH) | cb (g/L) | Γ | cm (g/L) | Jv (LMH) |
|---|
Model & validity
Element. 1-D plug flow along the feed channel, 60 cells per element. Mass transfer from Schock & Miquel (1987), Sh = 0.065 Re0.875 Sc0.25 with Re on the interstitial velocity and the spacer hydraulic diameter; channel friction λ = 6.23 Re−0.3. Film model cm = cp + (cb − cp) eJv/k with cp = (1 − Rint) cm. Osmotic pressure by van't Hoff (dilute solutions; for concentrated aqueous electrolytes use the Pitzer-based CP Calculator). Feed-channel cross-section Ach = A·h / 2L from the active area, element length and spacer thickness.
Cartridge. Concentration polarisation from the buoyancy-driven model of Geraldes et al. (2026), k = 0.28 U* 𝒢1/4 Sc−1/4, at the run-mean pressure and bulk concentration. The planner solves for the cartridge feed concentration whose run-mean wall concentration cm equals the area-mean cm of the element at the same ΔP. The same membrane at the same ΔP and cm gives the same flux, so the cartridge flux is taken equal to the element mean flux.
Properties. Solvent density and viscosity from handbook data (Andrade fit, 10–50 °C). Solute diffusivity by Stokes–Einstein with the molecular radius from molar mass and solute density; density contrast β = (1 − ρsolvent/ρsolute)/ρsolvent (volume-additive mixture). Both can be overridden by measured values.
Validity. The O-CMF correlation is dimensionless: it applies to any Newtonian solvent whenever the buoyancy parameter 𝒢, the Schmidt number and Jv/k lie in the ranges over which it was established (𝒢 7×10−9–5×10−6, Sc 700–2200, Jv/k ≤ 1.5). The planner checks each run and flags any group outside its range; Sc enters only as Sc−1/4, so a modest excursion in Sc is a short extrapolation. The property inputs — β (volume-additive mixture) and D (Stokes–Einstein) — are estimates: replace them by measured values when you have them. A quick in-system check: two runs at the same speed with feed concentrations differing by a factor of 3 should show the observed rejection rising with C0 as k ∝ C01/4 predicts. Solvents denser than the solute (DCM, chloroform, DMSO with most solutes) give β ≤ 0: the retained layer is stably stratified, the buoyancy model does not apply and the planner stops.
Rejection. A membrane's intrinsic rejection can depend on cm. That is why the planner reproduces cm instead of transferring a rejection measured elsewhere. When the element spans a wide range of cm (high recovery, several elements in series), the programme adds inlet, mid-length and outlet runs so the membrane is tested across that range.