Pot Still ABV Calculator: Run Simulator and Cuts

Simulate ethanol–water pot still ABV over time using Rayleigh batch distillation. Watch the boiler empty and compare heads, hearts, and tails for education and licensed distilling.

#ABV #batch distillation #distillation cuts #ethanol–water #pot still #Rayleigh distillation #VLE

Caution

The cuts model ethanol and water only, not methanol or other contaminants: they do not show that a distillate is safe to drink. Spirits for drinking need a tested, licensed process.

This pot still ABV calculator follows a batch of ethanol and water as ethanol leaves preferentially in the vapour. It is intended for licensed distillers comparing collection plans, chemical engineering students studying separation, and educators explaining why the first condensate is stronger than the liquid remaining in the boiler.

Enter the initial charge volume (20 L = 0.020 m³) and ABV, then choose the boiler strength at which the simulated run should stop. All strengths and liquid volumes refer to 20 °C. The educational composition option accepts ethanol mole fractions instead. The model uses atmospheric equilibrium data at approximately 101 kPa; it does not offer a pressure correction.

Enter the condensate volume collected in a one-minute sample before blending (50 cm³ per minute equals 3 L/h). The clock starts with the first condensate. Heat-up time is excluded. Holding this rate constant gives a useful comparison between runs, but it requires a separate operating assumption: the simulator does not calculate heater power, heat losses, or changing boil-up from an energy balance. The vapour molar rate changes with composition even when this liquid volume rate is fixed.

Press play or scrub the time control to see the boiler inventory fall and the collection jars fill. The outlet result is the instantaneous strength of newly condensed vapour. Each jar has its own final-volume scale, so equal liquid heights do not imply equal amounts. The ABV chart shows the collection windows and the declining boiler strength; the inventory chart shows the remaining liquid and the metered condensate.

Distillation cuts are user-defined windows here. Heads end after your selected collected volume. Hearts continue until the outlet ABV reaches your tails threshold, and tails continue to the selected boiler endpoint. The final table reports each cut’s blended volume, average ABV, and pure ethanol content. Average ABV is calculated from the combined component amounts, rather than by averaging ABV readings. Blending can change volume, so the sum of jar volumes can differ slightly from the metered collection volume.

The underlying Rayleigh distillation equation uses measured ethanol–water vapour–liquid equilibrium and a density interpolation that includes mixing contraction. The model assumes a well-mixed boiler, immediate equilibrium, and total condensation, with no reflux or pipe inventory. Real pot geometry, entrainment, passive reflux, dissolved solids, and rate changes can shift observed results. The dilute end of the equilibrium curve is an approximation; the density data support process estimates rather than certified alcoholometry.

The heads and tails labels do not predict methanol, other congeners, flavour, or product safety. This is a binary separation model, not an operating procedure or a recommended cut recipe. Alcohol production requires applicable authorization: US federal distilled-spirits approval and German alcohol-production permissions must be checked before practical use. Legal framing was checked on 29 October 2026; verify the rules currently in force with the relevant authority.

How it is calculated

The boiler is a closed, well-mixed ethanol–water batch. A differential component balance gives d(nx)=y∗(x) dnd(nx)=y^*(x)\,dn, and hence

ln⁡n0n(x)=∫xx0dξy∗(ξ)−ξ.\ln\frac{n_0}{n(x)}=\int_x^{x_0}\frac{d\xi}{y^*(\xi)-\xi}.

This is the Rayleigh equation presented in the IIT ISM laboratory manual, attached source s2, p.56. The equilibrium curve is a shape-preserving cubic interpolation (pchip) of measured 101 kPa data from Susial Badajoz et al. (2026), Table 5, DOI 10.1021/acsomega.6c04827. Two experimental points are withheld from interpolation for reference checks. A pure-water endpoint closes the dilute interval.

Use s=ln⁡(x0/x)s=\ln(x_0/x), so the integral becomes ∫0sx/(y∗−x) ds\int_0^s x/(y^*-x)\,ds. The tool evaluates a 601-point logarithmic composition grid with cumulative trapezoidal integration. Remaining moles are n=n0e−In=n_0e^{-I}; distilled moles and ethanol moles follow directly from total and ethanol balances.

At each composition, molar mass is M=xME+(1−x)MWM=xM_E+(1-x)M_W and cooled molar volume is v=M/ρ(x)v=M/\rho(x). Density comes from the experimental columns of Khattab et al. (2012), Table 1, attached s1 p.2. Linear interpolation from 293 to 298 K gives 293.15 K (20 °C); pchip interpolates composition. ABV is 100xME/(ρEv)100xM_E/(\rho_E v). An inverse bisection converts ABV inputs into mole fractions.

The metered collection volume is Vm=∫v(y∗) dDV_m=\int v(y^*)\,dD. Elapsed time is t=Vm/Qt=V_m/Q, where Q is the constant liquid take-off rate inferred from the one-minute collection sample. This permits changing vapour molar flow without assuming a heater energy balance. Cut boundaries use the heads metered volume and the outlet-ABV threshold. Each cut's component amounts are subtracted at its two boundaries; its combined composition determines its cooled volume and average ABV. Pure ethanol inventory closes independently across all cuts and the residue.

Assumptions

  • One equilibrium vapour–liquid stage; no reflux, rectification, entrainment, reaction, or noncondensable gas. All vapour is condensed and collected without delay or losses.
  • Fixed pressure approximately 101 kPa; binary ethanol and water only. Sugars, salts, acids, methanol, esters, and other congeners are excluded.
  • The liquid take-off rate is inferred from a user-entered one-minute volume sample and held constant at its 20 °C volume equivalent. The time axis starts at first collection and excludes heat-up, cooling delays, and warm equipment inventory.
  • All ABV and liquid volumes use 20 °C. The boiler drawing uses the cooled equivalent volume, not the thermally expanded hot liquid level. Jar height shows each jar's fraction of its own final volume.
  • Heads end at the selected metered volume. Tails start at the selected outlet ABV, after heads. If the threshold is above the initial outlet strength, there may be no hearts; if below the final outlet strength, there may be no tails. Those cases produce explanatory conditions.
  • The Rayleigh balance, the logarithmic integration coordinate, and cut component balances are derived here; the VLE and density values are measured source data. Molar masses use 46.06844 g/mol for ethanol and 18.01528 g/mol for water.
  • Charge and result volumes use m³ (US switch: ft³); the heads boundary and one-minute sample use cm³ (US switch: in³). The schematic retains explicitly labelled litres for readability.

Where it applies

Use for educational comparisons and process estimates for licensed distilling, not certified alcohol content, excise accounting, product release, or equipment design. Default input range: 1–60% initial ABV, 0.05–20% boiler stop ABV, with stop below charge; educational mole-fraction mode is limited to 0.005–0.6 initially. Runs stop with liquid remaining, not at a dry boiler.

The first measured nonzero VLE point is x=0.011 (roughly 3.4% boiler ABV); lower compositions interpolate toward the pure-water limit. Tail strength and long-run time there have greater model uncertainty. Density measurements are relatively sparse and have systematic differences from reference alcoholometric data. No accuracy claim finer than about one ABV percentage point is made. Numerical resolution is finer than these physical-data limits.

The supplied cut boundaries are illustrative controls, not advice for removing methanol or obtaining a safe or desirable product. Actual cuts depend on congeners, feedstock, equipment, and professional assessment. No safety conclusion follows from a predicted ABV.

For education and appropriately authorized distillers. US TTB guidance requires prior approval for distilled-spirits production; personal beverage production has no federal exemption. Germany's Alkoholsteuergesetz §32 restricts private alcohol production without the required permissions. These jurisdictional statements were checked on 29 October 2026; verify current federal, state, local, and German customs requirements with the responsible authority before practical operation.

References

  • IIT (ISM) Dhanbad, Advanced Chemical Engineering Laboratory Manual, Experiment 7, attached s2 p.56: Rayleigh differential-distillation balance. https://people.iitism.ac.in/~download/lab%20manuals/chemical_engg/Advanced_chemical_Engineering_Lab_Manual.pdf
  • Susial Badajoz, García Montesdeoca & Santiago (2026), A Novel Equipment for VLE Studies: Ethanol and Water System at Moderate Pressures and Data Verification with a Direct Point-to-Point Test, ACS Omega, Table 5. https://doi.org/10.1021/acsomega.6c04827 ; open text https://www.ebi.ac.uk/europepmc/webservices/rest/PMC13491997/fullTextXML
  • Khattab et al. (2012), Density, viscosity, and surface tension of water+ethanol mixtures from 293 to 323 K, Korean Journal of Chemical Engineering 29, 812–817, Table 1, attached s1 p.2. https://doi.org/10.1007/s11814-011-0239-6 ; open copy https://pcprakt.userpage.fu-berlin.de/SKRIPT/T13/DensityH2O_EtOH.pdf
  • LearnChemE, Single-Stage Batch Distillation: Summary. https://learncheme.com/quiz-yourself/interactive-self-study-modules/single-stage-batch-distillation/single-stage-batch-distillation-summary/ ; its separate ethanol–water simulation starts with an equimolar charge and constant heating, whereas this calculator uses a specified liquid collection rate: https://learncheme.com/simulations/separations/simple-batch-distillation-of-an-ethanol-water-mixture/
  • US Alcohol and Tobacco Tax and Trade Bureau, Alcohol FAQs G1 and distilled spirits permits; legal context checked 29 October 2026. https://www.ttb.gov/faqs/alcohol ; https://www.ttb.gov/regulated-commodities/beverage-alcohol/distilled-spirits/permits
  • Germany, Alkoholsteuergesetz §32(2), official consolidated text; legal context checked 29 October 2026. https://www.gesetze-im-internet.de/alkstg/__32.html
  • Perry's Chemical Engineers' Handbook, 7th ed., §13, p.13-96, worked Rayleigh batch example (100 mol, x0=0.18 to xf=0.06, 70.0 mol residue). Institutional excerpt https://schneider.cheme.cmu.edu/Files/Perry%27s%20Chemical%20Engineers%27%20Handbook%207th%20Edition/13.%20Distillation.pdf . Newer VLE data change the residue fraction by about 0.0094; see the comparison case.

Inputs

InputSymbolUnitDefault
Charge and stop compositionbasisABV at 20 °C
Initial liquid volume at 20 °C
20 L = 0.020 m³. The US switch displays cubic feet; strengths still refer to 20 °C.
V_chargem^30.02
Initial ABV (%)
ABV means pure ethanol volume divided by solution volume, both at 20 °C.
A_charge10
Stop at boiler ABV (%)A_stop0.5
Initial ethanol mole fractionx_charge0.18
Stop at boiler mole fractionx_stop0.06
Condensate in a one-minute sample
50 cm³ per minute equals 3 L/h. Assume this measured liquid collection rate stays constant throughout the run.
V_samplecm^350
End heads after collected volume
150 cm³ = 150 mL. User-defined collection boundary, not a recommendation or a methanol removal rule.
V_headscm^3150
Begin tails at outlet ABV (%)
Switch after heads when the instantaneous condensed vapour falls below this ABV. No congeners are calculated.
A_tails35
Run time ▶
Play or scrub from the first condensate to the selected boiler endpoint.
t_playmin30

Results

With the default inputs above, the calc gives:

ResultValue
Current outlet ABV (%)43.1
Boiler liquid at displayed time (20 °C)0.0185 m³
Collection duration to boiler stop144.9 min
Final hearts volume at 20 °C0.00237 m³
Final hearts average ABV (%)43.9
Initial ethanol recovered in hearts (%)52
Ethanol recovered in all cuts (%)96.8

Checked against sources

Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).

Withheld measurement: dilute ethanol–water equilibrium

Source: Independent held-out experimental point, Susial Badajoz et al., ACS Omega 2026, Table 5 (101 kPa): T=361.75 K, x1=0.071, y1=0.384; expanded composition uncertainty 0.003 (k=2). This row is excluded from the interpolation table. https://doi.org/10.1021/acsomega.6c04827

Inputs: Charge and stop composition = Ethanol mole fraction (education); Initial ethanol mole fraction = 0.071; Stop at boiler mole fraction = 0.01; the others as above.

ResultSourceThis calc
First vapour ethanol mole fraction0.384 ± 0.0060.3833

Withheld measurement: stronger ethanol–water equilibrium

Source: Independent held-out experimental point, Susial Badajoz et al., ACS Omega 2026, Table 5 (101 kPa): T=355.55 K, x1=0.199, y1=0.547. This row is excluded from the interpolation table; tolerance covers experimental uncertainty and interpolation. https://doi.org/10.1021/acsomega.6c04827

Inputs: Charge and stop composition = Ethanol mole fraction (education); Initial ethanol mole fraction = 0.199; the others as above.

ResultSourceThis calc
First vapour ethanol mole fraction0.547 ± 0.0060.5489

Measured mixed-liquid density near the charge range

Source: Khattab et al. (2012), Table 1, experimental density at 293 K and ethanol mole fraction 0.117 = 963.0 kg/m³. Model uses the measurements interpolated 0.15 K higher to 20 °C (about −0.16 kg/m³); tolerance 1 kg/m³ allows this stated temperature difference. This case verifies transcription and unit conversion, not an independent density model.

Inputs: Charge and stop composition = Ethanol mole fraction (education); Initial ethanol mole fraction = 0.117; Stop at boiler mole fraction = 0.01; the others as above.

ResultSourceThis calc
Initial liquid density at 20 °C963 ± 1 kg/m³962.8 kg/m³

Published Rayleigh batch example: residue fraction

Source: Perry's Chemical Engineers' Handbook, 7th ed., §13, p.13-96, simple batch distillation worked example: 100 mol initially at x=0.18 gives 70.0 mol residue at x=0.06 (fraction 0.700), using graphical integration and its older VLE table. Institutional excerpt: https://schneider.cheme.cmu.edu/Files/Perry%27s%20Chemical%20Engineers%27%20Handbook%207th%20Edition/13.%20Distillation.pdf . The present measured 2026 VLE curve gives a residue fraction about 0.7094 because its vapour compositions in this range are slightly richer; the 0.015 absolute tolerance allows the documented equilibrium-data and graphical-integration difference. The formulas are not adjusted to reproduce the older example.

Inputs: Charge and stop composition = Ethanol mole fraction (education); the others as above.

ResultSourceThis calc
Fraction of initial moles remaining0.7 ± 0.0150.7094

Close endpoint: honour the requested strength and conserve ethanol

Source: Independent reviewer reproduction supplied with this task (29 October 2026): otherwise default inputs, initial ABV 10%, requested boiler stop ABV 9.999%, zero heads volume. Reviewer gives duration 0.00941076 min; final boiler ABV must equal the requested endpoint and ethanol conservation requires zero relative balance residual. This is regression coverage of an input contract and conservation law, not an additional experimental validation.

Inputs: Stop at boiler ABV (%) = 9.999; End heads after collected volume = 0 cm³; the others as above.

ResultSourceThis calc
Actual boiler ABV at the endpoint (%)9.9999.999
Collection duration to boiler stop0.009411 min0.009411 min
Relative cut ethanol-balance residual (rounded to 12 decimals)00

Responses

Your email address will not be published. Required fields are marked *

L ading...