Adsorption Equilibrium, Fixed Beds and Ion Exchange

Engineering context

A packed bed treats a feed with intermittent contaminant spikes. A batch isotherm shows high equilibrium capacity, yet the effluent exceeds the specification sooner than expected. This is compatible with slow intraparticle diffusion, channeling or competing constituents rather than evidence the isotherm is wrong. Bed depth, flow, temperature and regeneration history should be documented. A service period must be scheduled to a leakage criterion, not the nominal equilibrium loading of every gram of media.

Unit learning outcomes

1. Develop and check adsorption equilibrium and selectivity using q_e=q_max K C/(1+K C); identify the principal assumptions and explain an experimental validation path.
2. Develop and check adsorption kinetics and film versus pore transport using dq/dt=k_LDF(q_star-q); identify the principal assumptions and explain an experimental validation path.
3. Develop and check fixed-bed breakthrough and regeneration using C_out/C_in=1/[1+exp(k(t_50-t))]; identify the principal assumptions and explain an experimental validation path.
4. Develop and check ion-exchange capacity and exhaustion using Q_use=m_resin q_eq eta_use; identify the principal assumptions and explain an experimental validation path.

Concept map

1. Adsorption equilibrium and selectivity

Langmuir’s isotherm posits equivalent adsorption sites and a monolayer-like saturation capacity. It is often a useful fitting equation even when its physical assumptions are imperfect. When concentration is low, the isotherm is approximately linear; at high concentration it approaches a maximum loading. An equilibrium fit does not predict uptake rate or packed-bed breakthrough. Separate the loading units, such as mg/g sorbent, from the equilibrium concentration units, such as mg/L solution.

Central relation: q_e=q_max K C/(1+K C). Operating boundary: These coefficients represent a specified adsorbent–solute–temperature combination; do not transfer them to another chemistry without data.

2. Adsorption kinetics and film versus pore transport

A linear-driving-force approximation replaces detailed diffusion in irregular porous particles with an effective rate constant. It can describe a measured uptake transient but does not uniquely identify the actual resistance. If the external film changes with mixing intensity, the fitted effective constant may change; if pore diffusion dominates, particle radius can have a stronger effect. Always inspect the early-time and late-time portions of kinetic data before extrapolating.

Central relation: dq/dt=k_LDF(q_star-q). Operating boundary: Single-exponential uptake assumes constant equilibrium concentration; coupled batch depletion needs a liquid-phase balance.

3. Fixed-bed breakthrough and regeneration

The logistic curve is a convenient empirical representation of a breakthrough front. It should not be mistaken for a universal solution to axial dispersion and adsorption kinetics. The front steepness is affected by the feed concentration, flow, bed length, pellet size, and isotherm nonlinearity. Breakthrough criteria are generally lower than complete exhaustion. Design around an allowable outlet concentration and a credible safety factor; regeneration efficiency and cycle scheduling determine effective capacity.

Central relation: C_out/C_in=1/[1+exp(k(t_50-t))]. Operating boundary: A logistic front does not resolve multicomponent displacement, irreversible fouling or changing adsorption temperature.

4. Ion-exchange capacity and exhaustion

Ion-exchange resins are rated on charge equivalents rather than a universal mass removal per gram. For hardness removal, divalent ion moles consume twice as many equivalents as monovalent ions. Working capacity is smaller than theoretical capacity because of leakage, incomplete regeneration and nonuniform utilization. Resin swelling, osmotic shock, fouling and chemical compatibility set additional constraints. Mass balances must include counterions leaving the resin, not just the target ion being removed.

Central relation: Q_use=m_resin q_eq eta_use. Operating boundary: The model ignores mass-transfer zone leakage and competitive species; verify bed performance experimentally.

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 08: Langmuir adsorption equilibrium](../labs/lab-08.html)
• [Lab 09: Linear-driving-force uptake](../labs/lab-09.html)
• [Lab 10: Fixed-bed breakthrough front](../labs/lab-10.html)
• [Lab 11: Ion-exchange working capacity](../labs/lab-11.html)

Further work

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

Unit 3, Lesson 1 — Adsorption equilibrium and selectivity

Purpose and operating question

Adsorption equilibria depend on temperature, solvent, concentration and surface chemistry. Why this matters: A packed bed treats a feed with intermittent contaminant spikes. A batch isotherm shows high equilibrium capacity, yet the effluent exceeds the specification sooner than expected. This is compatible with slow intraparticle diffusion, channeling or competing constituents rather than evidence the isotherm is wrong. Bed depth, flow, temperature and regeneration history should be documented. A service period must be scheduled to a leakage criterion, not the nominal equilibrium loading of every gram of media.

Physical interpretation

Langmuir’s isotherm posits equivalent adsorption sites and a monolayer-like saturation capacity. It is often a useful fitting equation even when its physical assumptions are imperfect. When concentration is low, the isotherm is approximately linear; at high concentration it approaches a maximum loading. An equilibrium fit does not predict uptake rate or packed-bed breakthrough. Separate the loading units, such as mg/g sorbent, from the equilibrium concentration units, such as mg/L solution.

The principal mathematical relationship used in this lesson is

q_e=q_max K C/(1+K C)

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: q_e=q_max K C/(1+K C). 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.

These coefficients represent a specified adsorbent–solute–temperature combination; do not transfer them to another chemistry without data.

Worked investigation

For qmax=125 mg/g and K=0.08 L/mg, calculate equilibrium loading at C=5, 20 and 100 mg/L. Determine the solute mass held by 1.5 kg of clean adsorbent at C=20 mg/L if every gram reaches the stated equilibrium. Explain why a finite batch may not have enough dissolved solute to reach this loading.

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 q_e=q_max K C/(1+K C) 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

The Langmuir term KC is 0.4, 1.6 and 8.0 at C=5, 20 and 100 mg/L respectively. Thus q=125(0.4/1.4)=35.71 mg/g, q=125(1.6/2.6)=76.92 mg/g and q=125(8/9)=111.11 mg/g. At C=20 mg/L and a mass of 1.5 kg, nominal equilibrium capacity is 1500 g×76.92 mg/g≈115.4 g. That mass is only possible if the contact solution supplies the required solute; an actual batch equilibrium must solve adsorbent loading and liquid depletion simultaneously.

Questions for independent work

1. Explain the physical meaning and units of every symbol in q_e=q_max K C/(1+K C) 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 3, Lesson 2 — Adsorption kinetics and film versus pore transport

Purpose and operating question

External film transfer and intraparticle diffusion often control the time to reach equilibrium. Why this matters: A packed bed treats a feed with intermittent contaminant spikes. A batch isotherm shows high equilibrium capacity, yet the effluent exceeds the specification sooner than expected. This is compatible with slow intraparticle diffusion, channeling or competing constituents rather than evidence the isotherm is wrong. Bed depth, flow, temperature and regeneration history should be documented. A service period must be scheduled to a leakage criterion, not the nominal equilibrium loading of every gram of media.

Physical interpretation

A linear-driving-force approximation replaces detailed diffusion in irregular porous particles with an effective rate constant. It can describe a measured uptake transient but does not uniquely identify the actual resistance. If the external film changes with mixing intensity, the fitted effective constant may change; if pore diffusion dominates, particle radius can have a stronger effect. Always inspect the early-time and late-time portions of kinetic data before extrapolating.

The principal mathematical relationship used in this lesson is

dq/dt=k_LDF(q_star-q)

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: dq/dt=k_LDF(q_star-q). 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.

Single-exponential uptake assumes constant equilibrium concentration; coupled batch depletion needs a liquid-phase balance.

Worked investigation

For qstar=48 mg/g, initial loading zero and LDF rate constant 0.035 min^-1, calculate q after 5, 20 and 60 min. Compute the half-time and show why increasing a batch contact time from 60 to 120 min has diminishing returns under this assumption. Identify two process observations that would contradict first-order uptake.

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 dq/dt=k_LDF(q_star-q) 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

For q=48 mg/g and k=0.035 min⁻¹, q(5)=48(1−e⁻⁰·¹⁷⁵)=7.71 mg/g; q(20)=24.16 mg/g; and q(60)=42.12 mg/g. Half-equilibrium time is ln2/0.035=19.80 min. At 120 min, loading approaches 47.28 mg/g, adding only about 5.16 mg/g beyond the one-hour point. If the observed uptake is initially linear but later follows a much slower long tail, a single constant-k model may be inadequate. Concentration depletion can also make a nominal q time-dependent.

Questions for independent work

1. Explain the physical meaning and units of every symbol in dq/dt=k_LDF(q_star-q) 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 3, Lesson 3 — Fixed-bed breakthrough and regeneration

Purpose and operating question

A moving mass-transfer zone produces a finite outlet breakthrough curve. Why this matters: A packed bed treats a feed with intermittent contaminant spikes. A batch isotherm shows high equilibrium capacity, yet the effluent exceeds the specification sooner than expected. This is compatible with slow intraparticle diffusion, channeling or competing constituents rather than evidence the isotherm is wrong. Bed depth, flow, temperature and regeneration history should be documented. A service period must be scheduled to a leakage criterion, not the nominal equilibrium loading of every gram of media.

Physical interpretation

The logistic curve is a convenient empirical representation of a breakthrough front. It should not be mistaken for a universal solution to axial dispersion and adsorption kinetics. The front steepness is affected by the feed concentration, flow, bed length, pellet size, and isotherm nonlinearity. Breakthrough criteria are generally lower than complete exhaustion. Design around an allowable outlet concentration and a credible safety factor; regeneration efficiency and cycle scheduling determine effective capacity.

The principal mathematical relationship used in this lesson is

C_out/C_in=1/[1+exp(k(t_50-t))]

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: C_out/C_in=1/[1+exp(k(t_50-t))]. 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.

A logistic front does not resolve multicomponent displacement, irreversible fouling or changing adsorption temperature.

Worked investigation

A measured bed follows a logistic approximation with t50=210 min and steepness k=0.032 min^-1. Determine C/C0 at 150, 210 and 260 min. Find when the outlet reaches 10% of feed concentration; compare that breakthrough time to t50. Discuss whether increasing the bed by 20% must increase t50 by exactly 20%.

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 C_out/C_in=1/[1+exp(k(t_50-t))] 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

At t=150 min, k(t50−t)=0.032×60=1.92, so C/C0=1/(1+e¹·⁹²)=0.1279. At t=210 min, C/C0=0.5. At 260 min, C/C0≈0.832. For the 10% leakage criterion, solve 0.1=1/[1+exp(k(t50−t))], yielding t=t50−ln9/k≈141.3 min. A design using 210 min as the service life would allow 50% of inlet contaminant to leave. Fitted logistic curves must not be extrapolated indefinitely beyond measured conditions.

Questions for independent work

1. Explain the physical meaning and units of every symbol in C_out/C_in=1/[1+exp(k(t_50-t))] 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 3, Lesson 4 — Ion-exchange capacity and exhaustion

Purpose and operating question

Equivalent charge and resin capacity constrain the amount of ions a column can exchange. Why this matters: A packed bed treats a feed with intermittent contaminant spikes. A batch isotherm shows high equilibrium capacity, yet the effluent exceeds the specification sooner than expected. This is compatible with slow intraparticle diffusion, channeling or competing constituents rather than evidence the isotherm is wrong. Bed depth, flow, temperature and regeneration history should be documented. A service period must be scheduled to a leakage criterion, not the nominal equilibrium loading of every gram of media.

Physical interpretation

Ion-exchange resins are rated on charge equivalents rather than a universal mass removal per gram. For hardness removal, divalent ion moles consume twice as many equivalents as monovalent ions. Working capacity is smaller than theoretical capacity because of leakage, incomplete regeneration and nonuniform utilization. Resin swelling, osmotic shock, fouling and chemical compatibility set additional constraints. Mass balances must include counterions leaving the resin, not just the target ion being removed.

The principal mathematical relationship used in this lesson is

Q_use=m_resin q_eq eta_use

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: Q_use=m_resin q_eq eta_use. 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.

The model ignores mass-transfer zone leakage and competitive species; verify bed performance experimentally.

Worked investigation

A bed contains 45 kg dry resin with rated capacity 3.2 eq/kg and 65% usable fraction. Calculate working capacity in equivalents. For inlet total hardness 4 meq/L and a desired treated-water volume, estimate exhaustion throughput in liters. Repeat at 6 meq/L; quantify the resulting reduction in service time at fixed flow.

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 Q_use=m_resin q_eq eta_use 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

Working capacity is 45 kg×3.2 eq/kg×0.65=93.6 equivalents. A feed concentration of 4 meq/L equals 0.004 eq/L, so the ideal service volume is 93.6/0.004=23,400 L. At 6 meq/L it is 15,600 L, one third lower. Convert to runtime by dividing the service volume by actual water flow in L/h. An ideal exhaustion balance is not a breakthrough guarantee, especially when competitive ions or a finite mass-transfer zone are significant.

Questions for independent work

1. Explain the physical meaning and units of every symbol in Q_use=m_resin q_eq eta_use 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 3: engineering project — Adsorption, packed beds and ion exchange

Design brief

A packed bed treats a feed with intermittent contaminant spikes. A batch isotherm shows high equilibrium capacity, yet the effluent exceeds the specification sooner than expected. This is compatible with slow intraparticle diffusion, channeling or competing constituents rather than evidence the isotherm is wrong. Bed depth, flow, temperature and regeneration history should be documented. A service period must be scheduled to a leakage criterion, not the nominal equilibrium loading of every gram of media.

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

• Adsorption equilibrium and selectivity — employ q_e=q_max K C/(1+K C) only within its assumptions. These coefficients represent a specified adsorbent–solute–temperature combination; do not transfer them to another chemistry without data.
• Adsorption kinetics and film versus pore transport — employ dq/dt=k_LDF(q_star-q) only within its assumptions. Single-exponential uptake assumes constant equilibrium concentration; coupled batch depletion needs a liquid-phase balance.
• Fixed-bed breakthrough and regeneration — employ C_out/C_in=1/[1+exp(k(t_50-t))] only within its assumptions. A logistic front does not resolve multicomponent displacement, irreversible fouling or changing adsorption temperature.
• Ion-exchange capacity and exhaustion — employ Q_use=m_resin q_eq eta_use only within its assumptions. The model ignores mass-transfer zone leakage and competitive species; verify bed performance experimentally.

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 05 — Nanofiltration rejection and concentrate

Unit: 2. Model boundaries: Observed, fixed rejection and constant volumetric densities; neglects charge/solvent coupling.

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

Feed solute concentration (mg/L), Observed rejection (%), Water recovery (%).

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 Water recovery (%) 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

Observed, fixed rejection and constant volumetric densities; neglects charge/solvent coupling. 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 06 — Gas permeation stage cut

Unit: 2. Model boundaries: Assumed permeate composition; the overall balance imposes feasibility limits.

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

Feed A fraction, Permeate A fraction, Stage cut.

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 Stage cut 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

Assumed permeate composition; the overall balance imposes feasibility limits. 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 3: 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: Adsorption equilibrium and selectivity

Q1.1 — Interpretation: Starting from q_e=q_max K C/(1+K C), identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q1.2 — Calculation: For qmax=125 mg/g and K=0.08 L/mg, calculate equilibrium loading at C=5, 20 and 100 mg/L. Determine the solute mass held by 1.5 kg of clean adsorbent at C=20 mg/L if every gram reaches the stated equilibrium. Explain why a finite batch may not have enough dissolved solute to reach this loading.

Q1.3 — Model validation: These coefficients represent a specified adsorbent–solute–temperature combination; do not transfer them to another chemistry without data. 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 q_e=q_max K C/(1+K C), uses consistent concentration or mass units, explicitly honors the limit (“These coefficients represent a specified adsorbent–solute–temperature combination; do not transfer them to another chemistry without data.”), 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: Adsorption kinetics and film versus pore transport

Q2.1 — Interpretation: Starting from dq/dt=k_LDF(q_star-q), identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q2.2 — Calculation: For qstar=48 mg/g, initial loading zero and LDF rate constant 0.035 min^-1, calculate q after 5, 20 and 60 min. Compute the half-time and show why increasing a batch contact time from 60 to 120 min has diminishing returns under this assumption. Identify two process observations that would contradict first-order uptake.

Q2.3 — Model validation: Single-exponential uptake assumes constant equilibrium concentration; coupled batch depletion needs a liquid-phase balance. 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 dq/dt=k_LDF(q_star-q), uses consistent concentration or mass units, explicitly honors the limit (“Single-exponential uptake assumes constant equilibrium concentration; coupled batch depletion needs a liquid-phase balance.”), 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: Fixed-bed breakthrough and regeneration

Q3.1 — Interpretation: Starting from C_out/C_in=1/[1+exp(k(t_50-t))], identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q3.2 — Calculation: A measured bed follows a logistic approximation with t50=210 min and steepness k=0.032 min^-1. Determine C/C0 at 150, 210 and 260 min. Find when the outlet reaches 10% of feed concentration; compare that breakthrough time to t50. Discuss whether increasing the bed by 20% must increase t50 by exactly 20%.

Q3.3 — Model validation: A logistic front does not resolve multicomponent displacement, irreversible fouling or changing adsorption temperature. 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 C_out/C_in=1/[1+exp(k(t_50-t))], uses consistent concentration or mass units, explicitly honors the limit (“A logistic front does not resolve multicomponent displacement, irreversible fouling or changing adsorption temperature.”), 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: Ion-exchange capacity and exhaustion

Q4.1 — Interpretation: Starting from Q_use=m_resin q_eq eta_use, identify the independent driving force and the dependent response. Give their dimensions and explain which parameter must be measured.

Q4.2 — Calculation: A bed contains 45 kg dry resin with rated capacity 3.2 eq/kg and 65% usable fraction. Calculate working capacity in equivalents. For inlet total hardness 4 meq/L and a desired treated-water volume, estimate exhaustion throughput in liters. Repeat at 6 meq/L; quantify the resulting reduction in service time at fixed flow.

Q4.3 — Model validation: The model ignores mass-transfer zone leakage and competitive species; verify bed performance experimentally. 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 Q_use=m_resin q_eq eta_use, uses consistent concentration or mass units, explicitly honors the limit (“The model ignores mass-transfer zone leakage and competitive species; verify bed performance experimentally.”), 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.