The bond-graph method was applied to study the hybrid membrane extraction,

which occurs in the system with an ion-exchange polymer membrane placed

between an aqueous and organic liquid phase. The numerical calculations were

carried out to describe the time evolution of the transport systems as dependent on

initial conditions for pH and the aqueous solution concentration as well as the

carrier concentration in a liquid membrane. The results of the numerical

calculations were compared with the experimental results

Szczepańska, Grażyna

Szczepański, Piotr

Wódzki, Romuald

English

Repository of Nicolaus Copernicus University

 XXV ARS SEPARATORIA – Toruń, Poland 2010  254 BOND-GRAPH DESCRIPTION AND SIMULATION  OF MEMBRANE PROCESSES:  HYBRID MEMBRANE EXTRACTION  Grażyna SZCZEPAŃSKA, Piotr SZCZEPAŃSKI, Romuald WÓDZKI  Nicolaus Copernicus University, Faculty of Chemistry,  87-100 Toruń, Gagarina 7, Poland  Abstract The bond-graph method was applied to study the hybrid membrane extraction, which occurs in the system with an ion-exchange polymer membrane placed between an aqueous and organic liquid phase. The numerical calculations were carried out to describe the time evolution of the transport systems as dependent on initial conditions for pH and the aqueous solution concentration as well as the carrier concentration in a liquid membrane. The results of the numerical calculations were compared with the experimental results.  1. INTRODUCTION  The classic, common methods of dispersive liquid-liquid extraction are affected by various inconvenient phenomena due to water uptake, loss of active components, instability, short time of operation, etc. These drawbacks can be overcome by applying the non-dispersive extraction with the interface immobilized by a porous membrane [1,2]. In the case of ionic substances, this type of membrane can be replaced by dense ion-exchange membrane which leads to the system of the hybrid membrane extraction (HME) [3,4]. Ion-exchange membranes (IEM) are expected to mediate the extraction occuring as the ion-exchange processes between the external aqueous solution and the reactive extractant (liquid ion-exchanger-specific reagent dissolved in the organic phase).  The hybrid membrane extraction can be theoretically studied by  a quantitative description of the involved processes by thermodynamic network analysis (TNA) with the formalism of the bond-graph method [5]. On this basis, the mathematical model in the form of a set of ordinary differential equations was proposed and used for calculating the evolution of local concentrations and fluxes. Having assumed model parameters and experimentally estimated transport conditions, the mathematical model was numerically solved using the Berkeley Madonna (ver.8.0.1.) program. The numerical calculations were carried out to describe the time evolution of the transport systems as dependent on the initial conditions for pH and the  XXV ARS SEPARATORIA – Toruń, Poland 2010  255 aqueous solution concentration as well as the carrier concentration in a liquid membrane. To verify the model correctness and applicability, the calculations were compared with the experimental results.  2. EXPERIMENTAL  The experiments were performed in the system with the Nafion-117 (17 cm2) cation-exchange polymer membrane (CEM) separating the aqueous solution of KNO3 (1000 cm3) and the organic phase (50 cm3) composed of D2EHPA (di(2-ethylhexyl)phosphoric acid) dissolved in kerosene. In the system, extracted cations are transported by cation exchange groups in the polymer membrane and then desorbed into the organic phase by a liquid exchanger (extractant).  3. RESULTS AND DISCUSSION  3.1. EXPERIMENTAL VERIFICATION OF THE NUMERICAL CALCULATIONS The comparison of theoretically calculated and typical experimental dependence of the concentration [K]oc of extracted K+ cations in the organic phase (in the form of 3( )CK CH⋅  complex) vs. time of transport is presented in Fig.1a-b. The results correspond with the organic containing D2EHPA of the concentration equal to 0.01 or 0.325 mol/dm3, respectively. The first dependence (Fig.1a) exhibits a maximum of experimental concentrations which can be ascribed to the system changes caused by the osmotic permeation of water. In the second case (Fig. 1b) a regular, typical extraction curve with the concentration [K]oc attaining an equilibrium was observed.             Fig.1. Concentration of K+ - D2EHPA complex in organic phase vs. time  of extraction; concentration of the aqueous phase: 0.01 mol KNO3/dm3, pH=5.95, concentration of D2EHPA in the organic phase: (a). 0.1 and (b). 0.325 mol/dm3 Continuous lines were calculated by respective TNA model. Time [h]0 5 10 15 20 25 30[K]oc  [mol/dm3]0.05.0e-61.0e-51.5e-52.0e-52.5e-53.0e-5[K]oc -  experimental[K]oc  - calculationsTime [h]0 5 10 15 20 25[K]oc  [mol/dm3]0.05.0e-51.0e-41.5e-42.0e-4[K]oc -  experimental[K]oc  - calculationsa) b)  XXV ARS SEPARATORIA – Toruń, Poland 2010  256 Time [h]0 5 10 15 20 25[K]oc  [mol/dm3]01e-52e-53e-54e-55e-5Concentration of the feed solution [mol/dm3]0.00 0.02 0.04 0.06 0.08 0.10Extraction rate  [mol/cm2s]02e-114e-116e-118e-111e-10experimentalcalculationsIn order to verify the hypothesis concerning negative influence of water uptaken into the organic phase on the HME efficiency, an additional experiment with the same system but with water pervaporation (PV) from the organic phase was carried out. The results in Fig.2 indicate that the HME-PV system exhibits a regular shape typical for the extraction kinetics curves. This means that osmotic permeation of water can diminish equilibrium concentration of extracted species, although the maximum of extraction rates in the HME and HME-PV systems are comparable, i.e. 1.45 and 1.38 mol/cm2s, respectively.            3.2. INFLUENCE OF FEED CONCENTRATION  The properties of the HME system were studied as dependent on the feed concentration with KNO3 conc. changing in the range: 0.005-0.1 mol/dm3 (with pHf, t=0 = 5.95). The results of the investigation are presented in Fig.3 as a dependence of the maximum extraction rate (mol/cm2s) on the initial feed KNO3 concentration.              It can be concluded from the plots in Fig.3 that the model calculations and experimental results coincide in the region of relatively low concentrations, i.e. up to 0.02 mol/dm3. Qualitatively however, the both dependencies are similar and characteristic for consecutive processes.  Fig.2. Extraction of K+ into organic phase with water pervaporation: aqueous feed solution: 0.01 mol KNO3 /dm3; pHf,t=0 =5.95,  organic phase: 0.1 mol D2EHPA/dm3, cation exchange membrane: Nafion-117, pervaporation membrane: Nafion-120. Fig.3 Experimental and theoretical (from TNA) dependence of maximum extraction rate on KNO3 concentration in the feed solution; organic phase: 0.1 mol D2EHPA/dm3 of kerosene; pHf, t=0 =5.95.   XXV ARS SEPARATORIA – Toruń, Poland 2010  257 pH of the feed solution2 3 4 5 6Extraction rate [mol/cm2s]0.03.0e-126.0e-129.0e-121.2e-111.5e-11experimentalcalculationsConcentration of D2EHPA [mol/dm3]0.0 0.2 0.4 0.6 0.8Extraction rate  [mol/cm2s]0.07.0e-111.4e-102.1e-102.8e-103.5e-10experimentalcalculations3.3. INFLUENCE OF FEED pH One of most frequently studied extraction property of D2EHPA is its dependence on the pH of the feed solution. Therefore, a similar investigation was carried out for the HME system using the KNO3 feed solution of 0.01 mol/dm3 concentration and pH ranging from 2 to 6. The corresponding results are presented and compared with the results of model calculations in Fig.4.              The results of experiments and model investigations confirm that the extraction in the HME system is fast for pHf ranging from 3.5 to 6. Below this region extraction rate significantly decreases, and for pHf>6 the formation of emulsion in the organic phase was observed.  3.4. INFLUENCE OF D2EHPA CONCENTRATION ON EXTRACTION RATE The extraction rates in the HME system as dependent on the extractant concentration were studied. The experimental and calculated results are presented in Fig.5.             From the results presented in Fig.5, one can conclude that the extraction rate in the HME system is strongly dependent on D2EHPA concentration above Fig.4. Dependence of maximum extraction rate on pH of feed solution: 0.01 mol KNO3 /dm3, organic phase 0.1 mol D2EHPA/dm3 of kerosene.  Fig.5 Dependence of maximum extraction rate on extractant concentration; feed solution: 0.01 mol KNO3/dm3, pHf, t=0 =5.  XXV ARS SEPARATORIA – Toruń, Poland 2010  258 pH0 1 2 3 4 5 6 7log βo|f -3.8-3.6-3.4-3.2-3.0-2.8-2.6-2.4    ( )f btpH →∞|( )m o rtpH →∞( 0)f btpH =0.4 mol/dm3. For very high concentrations (0.7 mol/dm3), an increase of the organic phase viscosity can diminish the extraction rate which is observed as convex curvature of the studied dependence.   3.5. DEPENDENCE OF DISTRIBUTION COEFFICIENTS ON REACTION CONDITION The extraction properties of D2EHPA are known to depend on the composition of a reaction medium. However, the HME system differs from the classic one because the organic phase is not immediately contacted with the aqueous solution of the known composition but with the internal solution of the cation exchange membrane. From this reason, some calculations were carried out to correlate pH of the bulk feed solution with the local pH at the organic phase and cation exchange membrane (reaction sub-layer). The calculations were carried out for ( 0)f btpH = =2.4÷7 with other model constants and parameters assumed to be the same. The results of calculations are presented in Fig.6.             The calculated linear dependence of |logo fβ  on equilibrium ( )f btpH →∞  and |( )m o rtpH →∞  complies with a general theory and D2EHPA properties. Also, the results of calculations indicate that an increase of ( 0)f btpH =  above 4 does not change significantly its equilibrium value ( )f btpH →∞  and respective distribution coefficients |o fβ . This effect seems to be caused by a stabilizing function of the acidic cation exchange membrane. The plots in Fig. 6 allow predicting the lower distribution coefficients for the HME in comparison to the classic L-L extraction because of lower pH in the reaction medium.   4. CONCLUSION  Fig.6. Dependence of K+ distribution coefficient o fβ  on ( 0)f btpH =  and equilibrum ( )f btpH →∞  of feed solution and equilibrium |( )m o rtpH →∞  at the interface cation exchange membrane | organic solution m|o (model calculations).  XXV ARS SEPARATORIA – Toruń, Poland 2010  259 According to the results of this study, thermodynamic network analysis can be successfully exploited for the qualitative or quantitative description and prediction of the hybrid membrane extraction. For example, the results of calculations simulating the HME system properties comply with experimentally observed dependencies of extraction rates on quantities such as feed or stripping solution concentration, its pH, and the concentration of an extractant in the organic phase.   REFERENCES [1] R. Prasad, K.K. Sirkar, „Membrane-based solvent extraction” in: W.S.W. Ho,  K.K. Sirkar (Eds), Membrane Handbook, Chapman and Hall, New York 1992, 727-763. [2] S. Schlosser, E. Sabolova, R. Kertesz, L. Kubisova, J. Sep. Sci., 24 (2001) 509-518. [3] R. Wódzki, P. Szczepański, J. Membr. Sci., 197 (2002) 297-308. [4] R-S. Juang, H. Huang, J. Membr. Sci., 156 (1999) 179-186. [5] R. Wódzki, G. Szczepańska, P. Szczepański, Sep. Purif. Technol., 36 (2004) 1-16.  

Bond-graph description and simulation of membrane processes: hybrid membrane extraction

https://repozytorium.umk.pl/bitstream/handle/item/4114/P31%20Ars%20Separatoria%202010.pdf?sequence=1

Bond-graph description and simulation of membrane processes: hybrid membrane extraction

Abstract

Similar works

Full text

Available Versions

Repository of Nicolaus Copernicus University