Search NASA⌕ Search

Engineering topics

Chen, Chau-Chyun

Publications and source records attributed to Chen, Chau-Chyun.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adsorption Thermodynamics for Process Simulation

Adsorption has rapidly evolved in recent decades and is an established separation technology extensively practiced in gas separation industries and others. However, rigorous thermodynamic modeling of multicomponent adsorption equilibrium remains elusive, and industrial practitioners rely heavily on expensive and time-consuming trial-and-error pilot studies to develop adsorption units. Here, this article highlights the need for rigorous adsorption thermodynamic models and the limitations and deficiencies of existing models such as the extended Langmuir isotherm, dual-process Langmuir isotherm, and adsorbed solution theory. It further presents a series of recent advances in the generalization of the classical Langmuir isotherm of single-component adsorption by deriving an activity coefficient model to account for the adsorbed phase adsorbate–adsorbent interactions, substituting adsorbed phase adsorbate and vacant site concentrations with activities, and extending to multicomponent competitive adsorption equilibrium, both monolayer and multilayer. Requiring a minimum set of physically meaningful model parameters, the generalized Langmuir isotherm for monolayer adsorption and the generalized Brunauer–Emmett–Teller isotherm for multilayer adsorption address various thermodynamic modeling challenges including adsorbent surface heterogeneity, isosteric enthalpies of adsorption, BET surface areas, adsorbed phase nonideality, adsorption azeotrope formation, and multilayer adsorption. Also discussed is the importance of quality adsorption data that cover sufficient temperature, pressure, and composition ranges for reliable determination of the model parameters to support adsorption process simulation, design, and optimization.

09 BIOMASS FUELS↗
DOE OSTI2024

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data for Aqueous Solutions of Binary-Acid and Ternary-Acid Mixtures of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-Exchange Resin at Initial pH Levels of ∼3–7 and at 25–55 °C

This work is part of an effort to develop thermophysical property data and models supporting adsorptive process development for organic acid separation from a dilute aqueous solution of fermentation broth. It presents systematic experimental measurements for aqueous-phase adsorption equilibrium, kinetics, and column breakthrough for three binary-acid aqueous mixtures (acetic acid + lactic acid, butyric acid + lactic acid, and acetic acid + butyric acid) and one ternary-acid aqueous mixture (acetic acid + butyric acid + lactic acid) on Amberlite IRN-78 ion-exchange resin. The equilibrium measurements covered broad ranges of the initial acid concentration (100–400 mmol/L), initial pH (∼3–7), and temperature (25–55 °C). The equilibrium data for the binary-acid and ternary-acid aqueous mixtures indicate selective adsorption of lactic acid at initial pH levels of ∼3 and 4, while acetic acid and butyric acid are selectively adsorbed at initial pH levels of 5, 6, and 7. The subsequent kinetics and column breakthrough experiments were performed at 200 mmol/L with equimolar ratios, an initial pH of 6, and 25 °C. In conclusion, the measurements provide essential data sets for rigorous thermodynamic modeling and process simulation of adsorptive separation processes for organic acid separation from fermentation broth.

Adsorption↗
DOE OSTI2024

Adsorption Equilibrium, Kinetics, and Column Breakthrough Data of Acetic Acid, Butyric Acid, and Lactic Acid on IRN-78 Ion-exchange Resin at Initial pH ~3 – 7 and Temperature 25 – 55 °C

This work presents systematic aqueous-phase adsorption equilibrium, kinetics, and column breakthrough measurements with three key biointermediates that are common compounds in many bioprocesses. Adsorption equilibrium experiments were carried out with acetic acid, butyric acid, and lactic acid on a commercial ion-exchange resin, Amberlite IRN-78, at wide ranges of acid concentration (8-500 mmol/L), initial pH (similar to 3-7), and temperature (25-55 degrees C), simulating the effluent characteristics from different fermenter operations. The kinetics and column breakthrough experiments were conducted at an initial pH of 6 and a concentration of 200 mmol/L. The equilibrium study shows a higher loading at the initial pH < pK(a) and a lower loading at the initial pH > pK(a). Overall removal varies between 16 and 99% depending on the initial pH, temperature, and organic acid concentration and type. The study further indicates monolayer adsorption at the equilibrium pH > 10 and multilayer adsorption at the equilibrium pH < 6. The thermodynamic modeling of adsorption isotherm data was carried out using Langmuir and Freundlich isotherms. IRN-78 presents fast adsorption kinetics as the maximum loading was attained in <= 10 min and nearly the same breakthrough time for all three organic acids involved in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Generalized Brunauer–Emmett–Teller Isotherm for Mixed-Gas Multilayer Adsorption Equilibria

Here, this work presents a rigorous thermodynamic framework to calculate mixed-gas multilayer adsorption equilibria from single-component isotherms and vapor–liquid equilibria. Named the generalized Brunauer–Emmett–Teller (gBET) isotherm, the newly formulated isotherm considers adsorbent surface heterogeneity, competitive adsorption on the monolayer, condensation–evaporation on the subsequent layers, and adsorbed phase nonideality for the monolayer and the subsequent layers. The monolayer adsorbed phase nonideality is tracked by using an area-based adsorption nonrandom two-liquid activity coefficient model. The adsorbed phase composition and corresponding nonideality in the subsequent layers are calculated at the dew point condition of the mixed gas with either an equation of state or an activity coefficient model for the vapor–liquid equilibria. The proposed model is validated with six single and three binary multilayer adsorption equilibrium systems, and the model results are compared against those from the classical BET isotherm for single-component adsorption and ideal adsorbed solution theory for mixed-gas adsorption equilibria.

09 BIOMASS FUELS↗
DOE OSTI2023

Thermodynamic modeling of double azeotropy systems

Azeotropes are boiling liquid mixtures which have identical vapor and liquid compositions. Known as polyazeotropy, azeotropic mixtures may exhibit more than one azeotrope at a given temperature or pressure. While the thermodynamic condition for azeotrope formation in vapor-liquid equilibrium is well established, thermodynamic modeling of mixtures exhibiting polyazeotropy remains challenging and prior modeling studies often failed to satisfactorily correlate the polyazeotropic behavior of associating systems. Based on the association NRTL model, this study presents successful correlation and analysis of double azeotrope formation for three diverse systems: strongly associating methanol + diethylamine binary, weakly associating ammonia + pentafluoroethane binary, and non-associating benzene + hexafluorobenzene binary systems. Furthermore, these systems exhibit polyazeotropic behavior as the temperature increases and approaches their Bancroft points.

42 ENGINEERING↗
DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗
DOE OSTI2023

Addressing spreading pressure dependence of real adsorbed solution theory with generalized Langmuir isotherm

This work addresses spreading pressure dependence of Real Adsorbed Solution Theory (RAST) for mixed-gas adsorption equilibria using generalized Langmuir (gL) isotherm. Considering vacant sites as an integral part of competitive multicomponent adsorption on a constant adsorbent surface area, the gL isotherm properly accounts for surface heterogeneity and loading, adsorbate composition, and temperature dependence. We show the spreading pressure dependence of adsorbate activity coefficient expression in the RAST framework can be generated from the gL isotherm. The procedure is illustrated with a spreading pressure dependent adsorption Nonrandom Furthermore, two-Liquid activity coefficient model, and the results are validated for ten binary mixed-gas adsorption equilibria systems including two highly nonideal azeotropic systems.

42 ENGINEERING↗
DOE OSTI2023

Comparison of Heterogeneous Langmuirian Models for Mixed-Gas Adsorption Equilibria

Mixed-gas adsorption remains a very active area of research, where newly proposed adsorption equilibrium models need to be compared to standards in the field. The generalized Langmuir (gL) isotherm, proposed by Hamid et al. to address adsorption of different size molecules on heterogeneous adsorbents without requiring evaluation of the Gibbs isotherm, has thus far been tested mainly against the adsorbed solution theory (AST). The gL isotherm was shown to be thermodynamically consistent and gave either better or similar representation than AST models. Here, the gL isotherm is contrasted with the Dual-Process Langmuir isotherm popularized by Ritter and his co-workers and the Loading Ratio Correlation popularized for simple mixtures, both of which are considered the “standard” models for process simulation. Furthermore, the comparison is based on their theoretical roots, interpretations, and representation of the equilibrium loading and the isosteric enthalpy of adsorption of unary, binary, and ternary mixtures of various gaseous adsorbates on different adsorbents over a range of experimental pressures and temperatures.

42 ENGINEERING↗
DOE OSTI2023

Comparative Energetics of Various Membrane Distillation Configurations and Guidelines for Design and Operation

This paper presents a comparative performance study of single-stage desalination processes with major configurations of membrane distillation (MD) modules. MD modules covered in this study are (a) direct contact MD (DCMD), (b) vacuum MD (VMD), (c) sweeping gas MD (SGMD), and (d) air gap MD (AGMD). MD-based desalination processes are simulated with rigorous theoretical MD models supported by molecular thermodynamic property models for the accurate calculation of performance metrics. The performance metrics considered in MD systems are permeate flux and energy efficiency, i.e., gained output ratio (GOR). A general criterion is established to determine the critical length of these four MDs (at fixed width) for the feasible operation of desalination in a wide range of feed salinities. The length of DCMD and VMD is restricted by the feed salinity and permeate flux, respectively, while relatively large AGMD and SGMD are allowed. The sensitivity of GOR flux with respect to permeate conditions is investigated for different MD configurations. AGMD outperforms other configurations in terms of energy efficiency, while VMD reveals the highest permeate production. With larger MD modules, utilization of thermal energy supplied by the hot feed for evaporation is in the order of VMD > AGMD > SGMD > DCMD. Simulation results highlight that energy efficiency of the overall desalination process relies on the efficient recovery of spent for evaporation, suggesting potential improvement in energy efficiency for VMD-based desalination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Estimation of isosteric heat of adsorption from generalized Langmuir isotherm

Simulation and design of adsorptive separation units demand accurate estimation of thermodynamic properties. Isosteric heat of adsorption as calculated from generalized Langmuir (gL) isotherm coupled with Clausius-Clapeyron expression for pure component and mixed-gas adsorption equilibria is presented in this work. The estimated isosteric heat of adsorption as functions of surface loading and composition is validated against the experimental data for various adsorption systems. Furthermore, the gL results are compared against classical Langmuir (cL) and Toth isotherm for pure components and with Ideal Adsorbed Solution Theory (IAST) for mixed-gas adsorption equilibria. Furthermore, the comparison highlights that gL outperforms cL and Toth for pure component adsorption and IAST for mixed-gas adsorption, and gL reliably captures the loading dependence and the composition dependence for isosteric heat of adsorption.

42 ENGINEERING↗
DOE OSTI2022

Thermodynamic modeling of aqueous lithium salt solutions with association electrolyte nonrandom two-liquid activity coefficient model

The high charge density of lithium ion and the resulting strong association phenomena make thermodynamic modeling of aqueous lithium electrolyte solutions extremely challenging. In this study, the association electrolyte nonrandom two-liquid activity coefficient model of Lin et al. (AIChE J. 2022, 68(2), e17422) is utilized to correlate and predict thermodynamic properties and solubility behavior of aqueous single electrolyte solutions of LiCl, LiBr, LiI, and LiNO 3 , and their mixed electrolyte solutions. Capturing self-association of water, cross-association of ion and water for hydration, and cross-association of cation and anion for ion-pairing, the association model accurately represents the literature experimental data up to saturation concentrations and at the temperature ranging from 263 K to 523 K. Here, this study further investigated the effect of anions of the lithium salts, and re-confirmed that the order of solution non-ideality as LiI > LiBr > LiCl > LiNO 3 because the anions with stronger association strengths are more likely to form ion pairs and thus lower the mean ionic activity coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Process Model-Based Validation of the Intensification of Biomass Fast Pyrolysis in a Fluidized Bed via Autothermal Operation

A model for the fluidized-bed pyrolysis of biomass is extended to enable the simulation of intensified autothermal operation. In this system, partial oxidation of char and pyrolysate species replaces an external heat source to supply the enthalpy of pyrolysis, greatly increasing the throughput of a given fast pyrolysis reactor. Oxidation reactions are compiled from CRECK and other literature sources and extended to cover the species found in pyrolysis products, including a derivation of the catalytic effect of ash on char combustion from experimental studies at Iowa State University (ISU). Results indicate a roughly 3-fold increase in biomass throughput for a given reactor volume at the pilot scale with minimal loss of valuable products, in agreement with published data from the 3 × 10 -2 m 3 reactor operated by ISU. Additionally, a proposal for a 250 ton/day commercial-scale biomass pyrolysis reactor is analyzed, showing that the same size reactor operated by heating the fluidizing gas externally would have only one-tenth the capacity; if operated by heating the sand externally, the sand would need to circulate at an impractical 0.7 bed volume per minute rate to maintain a 250 ton/day biomass throughput.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Thermodynamic modeling of aqueous polyelectrolyte solutions with mixed-valent counterions

Thermodynamic modeling of aqueous polyelectrolyte solutions with salts is of significant interest for many industrial applications. Here this study applies the polyelectrolyte Nonrandom Two-liquid activity coefficient model to aqueous polyelectrolyte solutions with mixed-valent counterions. A modified Delocalized Binding Theory was proposed to determine the polyion condensation fractions of the mixed counterions. This modified theory accounts for the electrostatic binding of the counterions on the polyion, the dissociation entropy of the counterions, and the electrostatic interactions between the uncondensed ionic species. Given the polyion condensation fractions, the critical value of Manning’s parameter ξ and the amounts of uncondensed polyions and counterions can be calculated along with the activity coefficients of mobile ions. The model successfully correlates experimental data for various aqueous polyelectrolyte systems with mixed-valent counterions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Process Modeling of CO 2 Absorption with Monoethanolamine Aqueous Solutions Using Rotating Packed Beds

A first-principle process simulation model is presented for the chemical absorption of carbon dioxide (CO 2 ) with monoethanolamine (MEA) aqueous solutions using rotating packed beds (RPB). Built on a proven rate-based packed bed absorber model, the RPB model rigorously simulates the phase and chemical equilibria at the vapor-liquid interface, the heat and mass transfer across the gas and liquid films, the fast reactions between MEA and CO 2 in the liquid film, and the RPB hydraulics. Here, estimation of the mass transfer rate across the liquid film is central to accurate simulation of the CO 2 absorption process with MEA aqueous solutions. We show that the literature lab-scale RPB data for CO 2 removal efficiency can be satisfactorily correlated by introducing a correction factor for the effective packing surface area predicted by the Onda correlation. Given the validated RPB model, we further show that, among the gas-phase mass transfer coefficient, the liquid-phase mass transfer coefficient, and the reaction rate constant for the reaction between amine and CO 2 , the reaction rate constant is the controlling step with the highest potential to enhance the CO 2 absorption performance in RPB.

42 ENGINEERING↗
DOE OSTI2021

Thermodynamic analysis of hydrogel swelling in aqueous sodium chloride solutions

Here, a comprehensive thermodynamic model is presented for swelling and salt partitioning of poly (n-isopropyl acrylamide) hydrogels in aqueous sodium chloride (NaCl) solutions at 298 K. Together with a Helmholtz energy expression for network elastic energy, a modified electrolyte Nonradom Two-Liquid activity coefficient model is used to represent both the long-range electrostatic interactions and the short-range van der Waals interactions present in the gel phase. The model parameters include one network parameter per hydrogel system plus two binary parameters per interaction pair in the polymer–solvent-salt systems. With semiquantitave agreement with experimental data for both hydrogel swelling and salt partitioning up to salt saturation, the analysis suggests the ion hydration between NaCl and water is the root cause for hydrogel deswelling in aqueous NaCl solution at high NaCl concentrations, i.e., NaCl weight fraction >0.03. On the other hand, the strong repulsive interaction between NaCl and hydrogel polymer is the dominant factor for salt partitioning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗