Chromatography and Preparative Resolution

Engineering context

An analytical method resolves two protein peaks in a small injection. During scale-up the second peak tails and a collection cut sacrifices yield. At preparative loading, adsorption-site saturation and gradients can shift retention and broaden peaks. A separation decision must combine recovery and purity for the specific fraction window, and use a validated calibration curve where detector response factors differ. Pressure drop and cleaning compatibility also influence the acceptable loading.

Unit learning outcomes

1. Develop and check retention, selectivity and capacity factor using k_prime=(t_R-t_0)/t_0; identify the principal assumptions and explain an experimental validation path.
2. Develop and check plate theory and band broadening using H=A+B/u+C u; identify the principal assumptions and explain an experimental validation path.
3. Develop and check peak resolution, detection and fraction cuts using R_s=2(t_R2-t_R1)/(w_1+w_2); identify the principal assumptions and explain an experimental validation path.
4. Develop and check gradient elution, scaling and preparative recovery using ln k(phi)=ln k_0-S phi; identify the principal assumptions and explain an experimental validation path.

Concept map

1. Retention, selectivity and capacity factor

The retention factor compares residence time beyond an unretained tracer with the void time. It is neither a partition coefficient nor a resolution metric. Two compounds can have strongly different retention but still overlap if their bands are sufficiently broad. When the mobile-phase composition changes during a gradient, a single constant k prime can be misleading. A reported retention time must be interpreted with column void time and flow rate.

Central relation: k_prime=(t_R-t_0)/t_0. Operating boundary: Retention factors are meaningful only for a specified mobile phase, column, temperature and measurement convention.

2. Plate theory and band broadening

The van Deemter relationship offers a useful view of why extremely slow and extremely fast mobile-phase velocities can both give poor efficiency. The B/u term is associated with longitudinal diffusion, while the C u term represents finite equilibration rates. Coefficients vary with packing, chemistry and operating conditions. More theoretical plates generally narrow peaks, but high efficiency alone does not guarantee selectivity for a difficult compound pair.

Central relation: H=A+B/u+C u. Operating boundary: The equation is a semi-empirical teaching correlation; values differ with solute and mobile phase.

3. Peak resolution, detection and fraction cuts

Chromatograms are distributions of signal in elution time, not discrete compound arrival events. Peak widths are usually defined by an agreed convention, often baseline widths or widths at half height, and formulas must match that convention. A collection window can increase target purity by discarding shoulders, but recovery then decreases. If detector signals have different response factors, peak areas cannot be equated to mass without calibration.

Central relation: R_s=2(t_R2-t_R1)/(w_1+w_2). Operating boundary: Gaussian peaks omit tailing from adsorption heterogeneity and overloaded injections.

4. Gradient elution, scaling and preparative recovery

Solvent-strength models capture the exponential decrease in retention with a stronger mobile phase for certain chromatography modes. A preparative process must consider mass loading, pressure drop and detector linearity in addition to retention. A gradient that separates an analytical microgram injection may fail on a kilogram process because mass-transfer and saturation behavior change. Fraction cuts should be justified by purity and yield measurements.

Central relation: ln k(phi)=ln k_0-S phi. Operating boundary: Only the local solvent-strength correlation is represented; time-dependent gradient propagation is beyond this model.

What to test experimentally

Establish the basis. A useful calculation begins with a consistent basis: mass rate, molar rate, solution volume, dry adsorbent mass or solvent mass. For equipment sizing a rate must not be confused with a cumulative inventory. Record every conversion as a separate line. Does the result represent an instantaneous rate, a time-integrated amount, a concentration, or a dimensionless ratio?

Identify equilibrium and rate assumptions. Equilibrium constrains an attainable composition or loading; kinetics determines how quickly a process approaches that constraint. Where an empirical coefficient is used, report its calibration range, units and experimental provenance. Which measured observation would show that an assumed equilibrium relation is not applicable?

Connect the model to a unit operation. A local model should be embedded within balances, operating constraints and product targets. Report useful recovery, purity and waste flows on comparable bases rather than optimizing a single local response. Which outlet becomes the limiting waste stream when product recovery increases?

Test uncertainty and reversal. If two alternatives are close, measurement error may reverse their apparent ranking. Estimate the sensitivity of performance to the largest uncertain inputs, and check that the conclusion survives a credible parameter range. Could an operating choice preferred at the reference feed become inferior when feed quality changes?

Laboratories supporting this unit

• [Lab 12: Chromatographic peak resolution](../labs/lab-12.html)
• [Lab 13: van Deemter optimal velocity](../labs/lab-13.html)
• [Lab 14: Solvent-strength retention](../labs/lab-14.html)

Further work

• [Unit 4 assessment](../assessments/unit-04.html)
• [Unit 4 engineering project](../projects/unit-04.html)

Unit 4, Lesson 1 — Retention, selectivity and capacity factor

Purpose and operating question

Migration speed depends on mobile- and stationary-phase partitioning. Why this matters: An analytical method resolves two protein peaks in a small injection. During scale-up the second peak tails and a collection cut sacrifices yield. At preparative loading, adsorption-site saturation and gradients can shift retention and broaden peaks. A separation decision must combine recovery and purity for the specific fraction window, and use a validated calibration curve where detector response factors differ. Pressure drop and cleaning compatibility also influence the acceptable loading.

Physical interpretation

The retention factor compares residence time beyond an unretained tracer with the void time. It is neither a partition coefficient nor a resolution metric. Two compounds can have strongly different retention but still overlap if their bands are sufficiently broad. When the mobile-phase composition changes during a gradient, a single constant k prime can be misleading. A reported retention time must be interpreted with column void time and flow rate.

The principal mathematical relationship used in this lesson is

k_prime=(t_R-t_0)/t_0

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: k_prime=(t_R-t_0)/t_0. Define the control volume or interface. Separate a measurable inlet/outlet quantity from an unobservable local state; do not fit a transport parameter to a mixed-stream observation without stating the model.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

Retention factors are meaningful only for a specified mobile phase, column, temperature and measurement convention.

Worked investigation

An unretained tracer emerges at 1.3 min and solutes A and B at 4.7 and 5.4 min. Compute each retention factor and the selectivity alpha=kB/kA for the more strongly retained component. Describe why changing flow rate can alter all three times without changing the equilibrium selectivity.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation k_prime=(t_R-t_0)/t_0 as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

Measure a full material balance over a representative interval and identify sensors whose uncertainty can invert the claimed trend. Measure a full material balance over a representative interval and identify sensors whose uncertainty can invert the claimed trend.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

For t0=1.3 min, kA=(4.7−1.3)/1.3=2.615 and kB=(5.4−1.3)/1.3=3.154. Selectivity α=kB/kA=1.206. This is relatively modest selectivity. Resolution also depends on peak width, so retention factors alone do not indicate whether preparative cuts are feasible. Changing flow can shorten both retention times, but column efficiency may change due to longitudinal diffusion and mass-transfer resistance. An empirical method should report the void-time determination and peak-width convention.

Questions for independent work

1. Explain the physical meaning and units of every symbol in k_prime=(t_R-t_0)/t_0 and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 4, Lesson 2 — Plate theory and band broadening

Purpose and operating question

Axial dispersion, eddy diffusion and mass transfer combine to set column efficiency. Why this matters: An analytical method resolves two protein peaks in a small injection. During scale-up the second peak tails and a collection cut sacrifices yield. At preparative loading, adsorption-site saturation and gradients can shift retention and broaden peaks. A separation decision must combine recovery and purity for the specific fraction window, and use a validated calibration curve where detector response factors differ. Pressure drop and cleaning compatibility also influence the acceptable loading.

Physical interpretation

The van Deemter relationship offers a useful view of why extremely slow and extremely fast mobile-phase velocities can both give poor efficiency. The B/u term is associated with longitudinal diffusion, while the C u term represents finite equilibration rates. Coefficients vary with packing, chemistry and operating conditions. More theoretical plates generally narrow peaks, but high efficiency alone does not guarantee selectivity for a difficult compound pair.

The principal mathematical relationship used in this lesson is

H=A+B/u+C u

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: H=A+B/u+C u. Choose a driving-force definition that matches the constitutive equation. A concentration gradient, pressure difference and electrochemical potential difference are not interchangeable.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

The equation is a semi-empirical teaching correlation; values differ with solute and mobile phase.

Worked investigation

Given A=0.014 cm, B=0.009 cm^2/s and C=0.055 s, calculate H at velocities 0.2, 0.4 and 0.8 cm/s. Differentiate H with respect to u to obtain u_opt=sqrt(B/C), and determine the plate count for a 30 cm column at the optimum. Explain which experimental measurements would test this model.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation H=A+B/u+C u as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

Hold all parameters but the candidate driving force constant; then document whether the observed output changes linearly, monotonically or nonmonotonically. Hold all parameters but the candidate driving force constant; then document whether the observed output changes linearly, monotonically or nonmonotonically.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

Differentiate H=A+B/u+C u to get dH/du=−B/u²+C. Setting this to zero gives uopt=√(B/C)=√(0.009/0.055)=0.4045 cm/s. Minimum H=A+2√(BC)=0.014+2√(0.000495)=0.0585 cm. A 30-cm column would then have N≈30/0.0585=513 theoretical plates. At u=0.2, H=0.070 cm; at u=0.8, H≈0.06925 cm. Similar non-optimal heights on both sides of the minimum illustrate competing diffusion and mass-transfer mechanisms.

Questions for independent work

1. Explain the physical meaning and units of every symbol in H=A+B/u+C u and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 4, Lesson 3 — Peak resolution, detection and fraction cuts

Purpose and operating question

Quantitative recovery depends on peak separation and selected collection windows. Why this matters: An analytical method resolves two protein peaks in a small injection. During scale-up the second peak tails and a collection cut sacrifices yield. At preparative loading, adsorption-site saturation and gradients can shift retention and broaden peaks. A separation decision must combine recovery and purity for the specific fraction window, and use a validated calibration curve where detector response factors differ. Pressure drop and cleaning compatibility also influence the acceptable loading.

Physical interpretation

Chromatograms are distributions of signal in elution time, not discrete compound arrival events. Peak widths are usually defined by an agreed convention, often baseline widths or widths at half height, and formulas must match that convention. A collection window can increase target purity by discarding shoulders, but recovery then decreases. If detector signals have different response factors, peak areas cannot be equated to mass without calibration.

The principal mathematical relationship used in this lesson is

R_s=2(t_R2-t_R1)/(w_1+w_2)

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: R_s=2(t_R2-t_R1)/(w_1+w_2). Estimate uncertainty and parameter sensitivity before investing in a more elaborate model. Distinguish uncertainty in measured properties from uncertainty in model structure.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

Gaussian peaks omit tailing from adsorption heterogeneity and overloaded injections.

Worked investigation

Two peaks have baseline retention times 8.0 and 9.4 min and widths 0.9 and 1.1 min. Compute Rs. Compare a second operating condition with retention times 8.3 and 9.8 min but wider peaks 1.2 and 1.4 min. Use Gaussian curves to estimate visually which fraction window loses more target solute.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation R_s=2(t_R2-t_R1)/(w_1+w_2) as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

At the reference operating point, perturb two high-impact inputs separately by about ten percent and record the response; interpret direction and magnitude. At the reference operating point, perturb two high-impact inputs separately by about ten percent and record the response; interpret direction and magnitude.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

Using baseline peak widths, resolution is Rs=2×(9.4−8.0)/(0.9+1.1)=1.40. Under the second operating condition, Rs=2×(9.8−8.3)/(1.2+1.4)≈1.154. Although the second chromatogram has a larger retention-time spacing (1.5 versus 1.4 min), its broader peaks reduce resolution. With Gaussian peaks, a fraction window selected to exclude impurity shoulders will sacrifice a different fraction of target signal under the two conditions. Actual collected mass requires detector-response calibration and careful use of baseline-width definitions.

Questions for independent work

1. Explain the physical meaning and units of every symbol in R_s=2(t_R2-t_R1)/(w_1+w_2) and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 4, Lesson 4 — Gradient elution, scaling and preparative recovery

Purpose and operating question

Mobile-phase strength changes retention and fractionation throughout a gradient run. Why this matters: An analytical method resolves two protein peaks in a small injection. During scale-up the second peak tails and a collection cut sacrifices yield. At preparative loading, adsorption-site saturation and gradients can shift retention and broaden peaks. A separation decision must combine recovery and purity for the specific fraction window, and use a validated calibration curve where detector response factors differ. Pressure drop and cleaning compatibility also influence the acceptable loading.

Physical interpretation

Solvent-strength models capture the exponential decrease in retention with a stronger mobile phase for certain chromatography modes. A preparative process must consider mass loading, pressure drop and detector linearity in addition to retention. A gradient that separates an analytical microgram injection may fail on a kilogram process because mass-transfer and saturation behavior change. Fraction cuts should be justified by purity and yield measurements.

The principal mathematical relationship used in this lesson is

ln k(phi)=ln k_0-S phi

The equation is a compact model rather than an instruction to ignore transport resistance, property dependence or other equipment conditions. Before using it, translate every mathematical term into a measurable quantity. Identify whether the driving force is a local value at an interface, an average over equipment, or a measured inlet-to-outlet difference. Units and operating conditions must accompany any quoted property coefficient.

Derivation strategy and limitations

Start from the model: ln k(phi)=ln k_0-S phi. Associate every ideal parameter with a realistic location in a piece of equipment: a membrane surface, packed bed, crystallizer slurry or filter cake.

The physical assumptions matter more than a successful algebraic manipulation. If the property varies with temperature or composition, calculate locally or integrate over the changing state rather than using a single value without justification. For time-dependent problems, distinguish a true transient solution from a sequence of pseudo-steady snapshots. For multiphase systems, check that a concentration or composition basis has not silently changed from mole to mass fraction.

Only the local solvent-strength correlation is represented; time-dependent gradient propagation is beyond this model.

Worked investigation

Consider k0=20, S=4 and solvent fraction phi=0.35. Estimate retention factor and tR for a void time of 1.5 min using a locally constant-composition approximation. Repeat at phi=0.50. Explain why these results cannot be used directly to forecast a continuously programmed gradient.

Solution approach: First list the information provided, and convert all values to consistent units. Substitute only the values that are supported by the statement; identify any additional properties or conditions that would need measurement for a full numerical result. Use the central relation ln k(phi)=ln k_0-S phi as the mathematical starting point, then evaluate a limiting or reference case. The final conclusion should distinguish an explicitly calculated response from a hypothetical plant-level extrapolation. Do not invent the unspecified flow, area, equilibrium or mass-transfer parameters to make a single definitive answer appear possible.

Engineering decision and diagnostic testing

Write an equipment-level feed/product/residue table. Do not claim whole-plant recovery from a local equilibrium or point-transport relation. Write an equipment-level feed/product/residue table. Do not claim whole-plant recovery from a local equilibrium or point-transport relation.

A meaningful experiment uses paired samples or simultaneous inlet/outlet measurements at a steady operating condition. Record the sampling interval and the instrument calibration or reference standard. If the observed response changes as throughput rises, test whether the change originates in the material property, the boundary layer, equilibrium constraints or residence-time distribution. When comparing equipment, express the outcome as recovery, product purity, energy or material consumption on the same basis.

Failure mechanisms and how to investigate them

• Inconsistent basis: Recalculate using one mass or molar basis and check whether every outflow is physically nonnegative. In a composition balance, sum the species to unity after solving.
• A parameter outside its fitted range: Compare the new operating condition with the measurements that established the coefficient. An apparently precise result may be unreliable under a different feed chemistry.
• An unmodeled resistance: Determine whether external transfer, internal diffusion, hydraulic pressure losses or phase equilibrium changes could constrain the response.
• An ambiguous measurement: Confirm the sample location and timing. An apparent difference in purity may be caused by mixing product samples across a transition.
• A neglected waste or utility duty: Draw the full separation train and account for regenerant, wash solvent, backwash, concentrate or mechanical power as applicable.

Sensitivity study

Retain the reference case from the worked investigation, then change the principal driving force by +10% while holding the model's other independent parameters fixed. Repeat for the strongest resistance or equilibrium coefficient. If the response is strongly nonlinear, a percentage change in the input does not imply the same percentage change in product. Where the result is dimensionless, verify that it remains within physical bounds. Where the result is a cumulative amount, confirm the time and inventory basis.

Fully worked numerical demonstration

Using k=k0 exp(−Sφ), at φ=0.35 the capacity factor is 20 exp(−1.4)=4.932, giving tR=1.5(1+4.932)=8.90 min. At φ=0.50, k=20 exp(−2)=2.707, so tR≈5.56 min. Increasing solvent strength reduces retention for this example. In a true gradient run, fluid of different composition propagates through the column; the analyte experiences a time-dependent solvent history, so neither isocratic estimate is an exact gradient retention time.

Questions for independent work

1. Explain the physical meaning and units of every symbol in ln k(phi)=ln k_0-S phi and identify one property that would need experimental calibration.
2. Solve the worked investigation with a justified perturbation to its largest uncertain input, then report the direction of the result and a dimensional check.
3. Describe a realistic counterexample where the relation would fail even though all arithmetic was correctly implemented.
4. Develop a plant-level material balance that includes both recovered product and rejected or retained material.
5. Compare the model's predicted effect of increasing the driving force with an actual equipment constraint and state which measurement would reveal the discrepancy.

Further study

Read the related unit chapter and complete the unit assessment. Use the relevant simulation in the [laboratory directory](../index.html#labs) to investigate a parameter change. Keep the model limitations visible when drawing an engineering conclusion.

Unit 4: engineering project — Chromatography and high-resolution separations

Design brief

An analytical method resolves two protein peaks in a small injection. During scale-up the second peak tails and a collection cut sacrifices yield. At preparative loading, adsorption-site saturation and gradients can shift retention and broaden peaks. A separation decision must combine recovery and purity for the specific fraction window, and use a validated calibration curve where detector response factors differ. Pressure drop and cleaning compatibility also influence the acceptable loading.

Deliverables and evidence

1. A stream table with all independent feeds, products, residues and recycle flows; define a common mass or molar basis.
2. A labeled process diagram naming the physical location of the driving force and the location of the main measurement.
3. A one-page model statement with every coefficient, unit, validity range and inferred rather than measured property.
4. A nominal calculation and at least four parameter-variation cases, including a counterexample outside the simplifying assumptions.
5. A safety and environmental note describing rejected material, regenerant, solvent or concentrate management.
6. A short conclusion that states what the simulation can support, what it cannot, and the next experimental validation needed.

Technical focus

• Retention, selectivity and capacity factor — employ k_prime=(t_R-t_0)/t_0 only within its assumptions. Retention factors are meaningful only for a specified mobile phase, column, temperature and measurement convention.
• Plate theory and band broadening — employ H=A+B/u+C u only within its assumptions. The equation is a semi-empirical teaching correlation; values differ with solute and mobile phase.
• Peak resolution, detection and fraction cuts — employ R_s=2(t_R2-t_R1)/(w_1+w_2) only within its assumptions. Gaussian peaks omit tailing from adsorption heterogeneity and overloaded injections.
• Gradient elution, scaling and preparative recovery — employ ln k(phi)=ln k_0-S phi only within its assumptions. Only the local solvent-strength correlation is represented; time-dependent gradient propagation is beyond this model.

Assessment rubric

• Material and energy accounting: 25%.
• Correct governing equations and units: 25%.
• Experimental test strategy: 20%.
• Sensitivity, uncertainty and limits: 20%.
• Clear engineering communication: 10%.

Safety boundary: No classroom model in this volume substitutes for validated pressure-equipment, chemical-handling or environmental-compliance design.

Laboratory 07 — Electrodialysis Faraday balance

Unit: 2. Model boundaries: Monovalent equivalent removal with fixed current efficiency; does not calculate stack energy.

Learning objective

Identify the physical driving force, calculate the principal response, and compare the effect of two independent inputs. Explain the difference between an equation that fits measurements and an equipment design that accounts for full mass and energy balances.

Adjustable variables

Stack current (A), Run time (min), Current efficiency (%), Diluate volume (L).

Experiment instructions

1. Record the reference-case output, its units and the assumed basis; confirm that the plotted point matches the numerical readout.
2. Change Stack current (A) to its lower and upper tested values. Export the CSV response, inspect monotonicity and identify any warnings or nonphysical outcomes.
3. Change one of the other parameters while holding the x-axis variable fixed. Explain the response in terms of the governing model and its dominant resistance or equilibrium contribution.
4. Estimate a practical equipment limitation omitted by this model, and describe a measurement that could reveal it. Do not claim equipment safety, product specification or regulatory compliance from this simulator alone.

Interpretation and model limits

Monovalent equivalent removal with fixed current efficiency; does not calculate stack energy. This is a transparent educational approximation with explicit parameter inputs. It is not a validated design program. The UI may omit temperature changes, feed composition dynamics, axial dispersion, fouling, energy duties or multicomponent coupling depending on the selected phenomenon. Use the experimental log to identify whether those omitted effects would change the conclusion.

Three questions

1. Which manipulated parameter most strongly changes the displayed response over the chosen operating range, and why?
2. Where can the response be nonphysical or the model cease to represent the real process? State a specific boundary.
3. How would you validate one point of the predicted curve using laboratory or vendor measurements?

Laboratory 08 — Langmuir adsorption equilibrium

Unit: 3. Model boundaries: Single-site capacity model on a fixed dry-sorbent mass basis.

Learning objective

Identify the physical driving force, calculate the principal response, and compare the effect of two independent inputs. Explain the difference between an equation that fits measurements and an equipment design that accounts for full mass and energy balances.

Adjustable variables

Maximum capacity (mg/g), Langmuir coefficient (L/mg), Equilibrium concentration (mg/L).

Experiment instructions

1. Record the reference-case output, its units and the assumed basis; confirm that the plotted point matches the numerical readout.
2. Change Equilibrium concentration (mg/L) to its lower and upper tested values. Export the CSV response, inspect monotonicity and identify any warnings or nonphysical outcomes.
3. Change one of the other parameters while holding the x-axis variable fixed. Explain the response in terms of the governing model and its dominant resistance or equilibrium contribution.
4. Estimate a practical equipment limitation omitted by this model, and describe a measurement that could reveal it. Do not claim equipment safety, product specification or regulatory compliance from this simulator alone.

Interpretation and model limits

Single-site capacity model on a fixed dry-sorbent mass basis. This is a transparent educational approximation with explicit parameter inputs. It is not a validated design program. The UI may omit temperature changes, feed composition dynamics, axial dispersion, fouling, energy duties or multicomponent coupling depending on the selected phenomenon. Use the experimental log to identify whether those omitted effects would change the conclusion.

Three questions

1. Which manipulated parameter most strongly changes the displayed response over the chosen operating range, and why?
2. Where can the response be nonphysical or the model cease to represent the real process? State a specific boundary.
3. How would you validate one point of the predicted curve using laboratory or vendor measurements?

Unit 4: assessment and answer guidance

These questions require a declared basis, explicit assumptions and units; no safety-critical design decision is authorized by classroom examples.

Topic 1: Retention, selectivity and capacity factor

Q1.1 — Interpretation: Starting from k_prime=(t_R-t_0)/t_0, identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q1.2 — Calculation: An unretained tracer emerges at 1.3 min and solutes A and B at 4.7 and 5.4 min. Compute each retention factor and the selectivity alpha=kB/kA for the more strongly retained component. Describe why changing flow rate can alter all three times without changing the equilibrium selectivity.

Q1.3 — Model validation: Retention factors are meaningful only for a specified mobile phase, column, temperature and measurement convention. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q1.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q1.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites k_prime=(t_R-t_0)/t_0, uses consistent concentration or mass units, explicitly honors the limit (“Retention factors are meaningful only for a specified mobile phase, column, temperature and measurement convention.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.

Topic 2: Plate theory and band broadening

Q2.1 — Interpretation: Starting from H=A+B/u+C u, identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q2.2 — Calculation: Given A=0.014 cm, B=0.009 cm^2/s and C=0.055 s, calculate H at velocities 0.2, 0.4 and 0.8 cm/s. Differentiate H with respect to u to obtain u_opt=sqrt(B/C), and determine the plate count for a 30 cm column at the optimum. Explain which experimental measurements would test this model.

Q2.3 — Model validation: The equation is a semi-empirical teaching correlation; values differ with solute and mobile phase. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q2.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q2.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites H=A+B/u+C u, uses consistent concentration or mass units, explicitly honors the limit (“The equation is a semi-empirical teaching correlation; values differ with solute and mobile phase.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.

Topic 3: Peak resolution, detection and fraction cuts

Q3.1 — Interpretation: Starting from R_s=2(t_R2-t_R1)/(w_1+w_2), identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q3.2 — Calculation: Two peaks have baseline retention times 8.0 and 9.4 min and widths 0.9 and 1.1 min. Compute Rs. Compare a second operating condition with retention times 8.3 and 9.8 min but wider peaks 1.2 and 1.4 min. Use Gaussian curves to estimate visually which fraction window loses more target solute.

Q3.3 — Model validation: Gaussian peaks omit tailing from adsorption heterogeneity and overloaded injections. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q3.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q3.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites R_s=2(t_R2-t_R1)/(w_1+w_2), uses consistent concentration or mass units, explicitly honors the limit (“Gaussian peaks omit tailing from adsorption heterogeneity and overloaded injections.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.

Topic 4: Gradient elution, scaling and preparative recovery

Q4.1 — Interpretation: Starting from ln k(phi)=ln k_0-S phi, identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q4.2 — Calculation: Consider k0=20, S=4 and solvent fraction phi=0.35. Estimate retention factor and tR for a void time of 1.5 min using a locally constant-composition approximation. Repeat at phi=0.50. Explain why these results cannot be used directly to forecast a continuously programmed gradient.

Q4.3 — Model validation: Only the local solvent-strength correlation is represented; time-dependent gradient propagation is beyond this model. Explain an experimental plan to identify whether this limitation is important at the proposed operating point.

Q4.4 — Decision: Compare a higher-throughput operating condition with a lower-throughput but more selective or stable condition. Identify the missing measurements needed to rank them for a plant.

Q4.5 — Fault diagnosis: Suppose an observed result deviates by 25% from the relation. Distinguish data-unit mistakes, instrument bias, a boundary condition change and genuinely invalid model assumptions.

Answer-check framework: A correct solution cites ln k(phi)=ln k_0-S phi, uses consistent concentration or mass units, explicitly honors the limit (“Only the local solvent-strength correlation is represented; time-dependent gradient propagation is beyond this model.”), and reports a physically sensible direction of change. Numerical calculations must show substitutions and a final unit. Alternative answers can be valid when a different, clearly documented operating assumption is made.