Predicting the performance of reverse osmosis membranes using optimization, mathematical modeling and solving the model with the help of numerical methods

Number of pages: 123 File Format: word File Code: 31773
Year: 2011 University Degree: Master's degree Category: Chemical - Petrochemical Engineering
  • Part of the Content
  • Contents & Resources
  • Summary of Predicting the performance of reverse osmosis membranes using optimization, mathematical modeling and solving the model with the help of numerical methods

    Academic Thesis for receiving a master's degree in environmental chemical engineering

    Abstract

    Among the models that are able to make successful predictions for the performance of reverse osmosis membranes, that is, the separation rate and flow flux in a wide range of physical conditions (membrane and solution) and laboratory conditions (pressure and concentration), the modification model surface force-cavity flow model and its more complete state i.e. modified generalized model of surface force-cavity flow.  In this research, efforts have been made to find the potential function and the friction function, which are two very important and effective functions in the mentioned models as well as other models of porous membranes, and valuable and new results have been obtained in this field. Examining the results showed that the best model for predicting the laboratory data is the proposed model with potential and friction functions, i.e. the Ex-P3-F3 model with the value of the target function equal to 0.0335. On the other hand, the investigation of the effect of the potential function on the results of the modified generalized model of surface force - pore flow, showed that in these conditions, replacing the proposed potential function in the mentioned model, the Ex-P4-F1 model is obtained, which, having the lowest value of the target function equal to 0.0337, predicts the laboratory data well.  On the other hand, today, due to the progress of mathematics and the invention of new numerical solution methods such as the finite difference method and the finite volume method, it has become possible to solve complex and non-linear equations. In this research, using these methods, the non-linear equations of speed and concentration governing the aforementioned models have been solved and finally, the performance of reverse osmosis membranes has been predicted based on these models.

    Key words: reverse osmosis membranes, modified generalized model of surface force - pore flow, potential function, friction function, numerical method.

    -1- Introduction

    In the late 1950s, Reed and Britton [1] at the University of Florida discovered that the membrane made of Cellulose acetate has the ability to remove salt. The flow of water passing through this dense membrane [2] was so low that it was impossible to use it in practice [1 and 2].

    The turning point of membrane phenomena can be considered to be around 1958, when the first asymmetric cellulose acetate membrane was made at the University of Los Angeles California (UCLA).

    Leob and Surirajan[3] was able to produce an asymmetric membrane of cellulose acetate by making a thin layer (0.1-1.0 µm thick) on a porous structure behind it]1 [.

    Cellulose acetate membranes have two major limitations. First, they are prone to attack by biological substances, in which case the incoming water must be chlorinated. Secondly, under different acidic and basic conditions, Cellulose acetate is hydrolyzed to acetate, so the pH of the system must be controlled between 5.4 and 5.7.

    One ??of the ways to develop asymmetric membranes was to produce the shell layer and the porous protective layer separately and then somehow put them on top of each other. In a new type of membrane, called thin film composite [4] (TFC), this process is performed in two steps. First, the thick layer is covered with a layer consisting of very small holes[5] and then a thin layer is covered over it. One of the merits of this method is that it is possible to optimize the performance of each layer separately. Nowadays, these types of membranes (TFC) are economically very affordable and have high performance [4].

     

    Polyamide membranes were made in 1977. The purpose of making these types of membranes, which are made by the TFC method, was to develop the chemical, biological and mechanical resistance of the surface and increase the useful life of the membrane. One of the limitations of this type of membrane compared to cellulose acetate membranes is the difficulty of their preparation, which incurs a higher cost [3 and 5].

     

     

     

    1-2- Introduction of the topic

    There are different points of view in membrane research. The science of modeling tries to describe and predict the performance of membranes at a higher level, and in this way, both to expand the knowledge related to the transfer mechanism and to use this knowledge in the design of membranes and reverse osmosis devices. Materials science, using physical and chemical criteria, strives to accelerate membrane process technology by preparing and synthesizing newer membranes. Process design researches in the field of general design and optimization of reverse osmosis systems on an industrial scale [6, 7].

    Among this, modeling and predicting the performance of reverse osmosis membranes [6] is particularly important. In order to adequately describe the performance of reverse osmosis membranes, very accurate and suitable mathematical models are needed. This need has led to the development of several transport models.

    The general purpose of transport models is to relate membrane performance, which is usually expressed in terms of separation (the fraction of solute that is separated from the feed) and flux, to operating conditions (such as pressure and feed concentration) or driving force (which are usually concentration and pressure gradients), by coefficients (known as transport coefficients and including model parameters). to give Coefficients or parameters should be determined from laboratory data. The complete success of a model means the power of that model to mathematically describe data with coefficients (or parameters) that remain constant within the range of practical conditions (parameters independent of driving forces). Finally, the model with specific transfer coefficients (or parameters) can predict the performance of a membrane in a wide range of operating conditions. This performance prediction power is the real power of the transmission model. The transfer model can sometimes eliminate high experimental costs, or it can be coupled with a research program in membrane manufacturing to lead to better design criteria, or it can be used together with a process design program and lead to a logical industrial scale for the reverse osmosis system [1]. So far, many models have been presented to express the performance and also predict the performance of reverse osmosis membranes, and some of them have had relatively good success. Over time, these models have become more complete and provided better results by clarifying the new facts about the separation mechanism and the relationships governing the model [8]. Among the models that have provided good results for describing reverse osmosis processes with very high accuracy, there is a model called the modified surface force-pore flow model [7] or the MD-SF-PF model and its more complete model, i.e. Ex-MD-SF-PF model [8], both of which were presented by Mahdizadeh and his colleagues] 10-9[.

    The beauty of these two models is that they have been able to predict the performance of membranes in different operating conditions with high accuracy based on a physical mechanism, in addition to determining the relationships governing reverse osmosis processes.

    The main problem in using these models Solving their nonlinear differential equations. To solve the equations of these two models, the numerical method of orthogonal superposition [9] has been used [11].

    In this research, first, several friction and potential functions were evaluated and selected from among the various functions, and after replacing them in the Ex-MD-SF-PF model with the help of programming in Matlab software and optimization, the ability and power of the new models in predicting the performance of reverse osmosis membranes has been examined, which can be seen later in the results section. It was found that not only the new models are able to predict the performance of reverse osmosis membranes, but sometimes some of the presented models have better results than the MD-SF-PF and Ex-MD-SF-PF models.

    On the other hand, the nonlinear equations governing the model were solved using finite difference [10] and finite volume [11] numerical methods, and then the results obtained for the membrane performance were compared with the laboratory results.

     

    Abstract:

     

    In order to explain and predict the behavior of reverse osmosis membranes and also design of reverse osmosis operation units, a good understanding of principals of transport in membranes is inevitable.

    This means that a powerful and efficient transport model has to be prepared.

  • Contents & References of Predicting the performance of reverse osmosis membranes using optimization, mathematical modeling and solving the model with the help of numerical methods

    List:

    Chapter One - Introduction. 1

    1-1- Introduction. 2

    1-2- Introduction of the topic. 3

    1-3- Objectives. 4

    1-4- Membranes and membrane processes. 5

    Chapter Two - Research theory and review of past research. 7

    2-1- Osmosis, osmotic pressure, reverse osmosis and views. 8

    2-2- Principles of membrane mass transfer. 9

    2-2-1- Membrane function. 9

    2-2-2- concentration polarization. 10

    2-2-3- driving forces for reverse osmosis. 12

    2-3- transmission mechanisms. 13

    2-3-1- sieve filtration mechanism. 13

    2-3-2- wet surface mechanism. 13

    2-3-3- Mechanism of dissolution-penetration. 14

    2-3-4- Preferential absorption mechanism - cavity flow. 14

    2-4- Transmission models. 15

    2-4-1- Mechanism-independent models or phenomenological models of transmission. 16

    2-4-1-1- Irreversible thermodynamics- phenomenological relations of transfer. 16

    2-4-1-1- A- Codem model - Katchalski. 17

    2-4-1-1- B- Codem model - Spiegler. 18

    2-4-2- Mechanism dependent models. 19

    2-4-2-1- non-porous transmission models. 19

    2-4-2-1- A- Dissolution-permeation model. 19

    2-4-2-1- B- Dissolution-permeation-cavity model. 20

    Page Title

    2-4-2-1- C. Developed Model of Dissolution- Infiltration. 21

    2-4-2-2- Transmission models based on porosity. 22

    2-4-2-2- A- Kimura-Surirajan model. 22

    2-4-2-2-B- Model of small holes. 24

    2-4-2-2-c- modified model- small holes. 25

    2-4-2-2- d- The surface force model - cavity flow. 27

    2-5- Modified model of surface force- cavity flow. 30

    2-5-1- determination of concentration distribution. 31

    2-5-2- Determination of speed distribution. 32

    2-5-3- Separation and passing fluxes of solvent and dissolved material inside the membrane. 33

    2-5-4- potential function. 34

    2-5-5- Friction function. 35

    2-6- Modified generalized model of surface force-hole flow. 35

    2-6-1- Radial component of dissolved substance flux. 35

    2-6-2- Axial component of dissolved substance flux. 37

    2-6-3- Determination of speed distribution. 39

    2-6-4- Separation and flux of dissolved substance and solvent passing through the membrane. 39

    2-6-5- potential function. 40

    2-6-6- Friction function. 40

    Chapter 3- Numerical methods of solving differential and nonlinear equations. 41

    3-1- Discretization theory. 42

    3-2- finite difference method. 43

    3-3- Finite element method. 44

    3-4- Limited volume method. 45

    3-4-1- central difference scheme. 47

    3-4-2- Upper hand dispute plan. 48

     

    Title                                                                                                                                                 . 49

    Chapter Four - Mathematical modeling and optimization. 50

    4-1- Introduction. 51

    4-2- The general method of solving MD-SF-PF and Ex-MD-SF-PF model equations. 51

    4-3- discretization of equations. 52

    4-3-1- Discretization of the speed equation. 52

    4-3-2- discretization of concentration equation. 54

    4-4- Solving algebraic equations. 59

    4-5- Optimization. 59

    Chapter Five - Results. 64

    5-1- Results of numerical solution and optimization of MD-SF-PF model. 65

    5-2- The results of the optimization of the Ex-MD-SF-PF model and other proposed models. 68

    5-3- Examining the results of Ex-MD-SF-PF and New-Ex-MD-SF-PF models. 74

    5-3-1- Comparing the results of the numerical solution of the Ex-MD-SF-PF model and the New-Ex-MD-SF-PF model with laboratory data 74

    5-3-2- Examining the concentration distribution of the Ex-MD-SF-PF model and the New-Ex-MD-SF-PF model. 81

    5-3-3- Speed ??distribution comparison of Ex-MD-SF-PF model and New-Ex-MD-SF-PF model. 83

    5-3-4- Examining and comparing the trend of changes in the potential function of the Ex-MD-SF-PF model and the New-Ex-MD-SF-PF model. 84

    5-4- Examining the effect of the potential function on Ex-MD-SF-PF model. 88

    5-4-1- Examining the speed distribution of the proposed model with functions (Ex-P4-F1). 93

    5-4-2- Examining the change process of the potential function of the proposed model with functions (Ex-P4-F1)           94

    Sixth chapter- Conclusion and suggestions. 96

    6-1- Conclusion. 97

    6-2- Suggestions.98

    References. 100

    Source:

    [1] Moradi.A., "Determining performance curves of reverse osmosis membranes using MD-SF-PF model",

    Master's thesis, Tarbiat Modares University, Department of Chemical Engineering (1371).

    [2] Reid, C.E. and Breton, E.J., "Water and Ion Flow Across Cellulosic Membranes", J.Appl. Polym. Sci., 1,113-143(1959).

    [3] Applegate, L.E., “Membrane Separation Processes”, Chem.Eng., June 11, 46-89 (1984).

     

    [4] Cadotte, J.E., and Peterson, R.J., “A New Thin-Film Composite Sea Water Reverse Osmosis Membrane”, Desalination, 1980; 32, 25-31. [5] Cadotte, J.E. "Evaluation of Composite Reverse Osmosis Membranes in Material Science of Synthetic Membranes", D.R. Lloyd, Ed. ACS Symposium series 269, Washington, D.C. (1985)

     

    [6] Kesting, R.E., “Synthetic Polymeric Membranes: A Structural Perspective”, 2nd ed. John Wiley and Sons, New York, NY (1985)

     

    [7] Sourirajan, S., “Reverse Osmosis and Synthetic Membranes: Theory-Technology-Engineering”, National Research Council Canada, Ottawa, Canada (1977).

     

    [8] Dickson, J.M., and Mehdizadeh, H., “Overview of Reverse Osmosis for Chemical Engineers, Part I: Fundamentals of Membranes Mass Transfer”, J.Eng. Islamic Rep. Iran, 1(4), 163-179 (1988).

    [9] Mehdizadeh, H. and Dickson, J.M., "Theoretical Modification of the Surface Force-Pore Flow Model for Reverse Osmosis Transport", J, Memb.Sci, 1989; 42, 119-145.

     

    [10] Mehdizadeh, H. and Dickson, J.M., "Evaluation of Surface Force-Pore Flow and Modified Surface Force-Pore Flow Modes for Reverse Osmosis Transport", Desalination, 1989; 24, 814.

     

    [11] Mehdizadeh, H., "Modeling of Transport Phenomena in Reverse Osmosis Membranes", PH.D. Thesis, McMaster University, Canada, 1990.

     

    [12] Taylor, J.S., Duranceau, S.J., Barrett, W.M., and Goigle, J.F., “Assessment of Potable Water Membrane Applications and Development Needs”, AWWA Res. Foun. And Amer. Water works Assoc. December, 1989.

     

    [13] Soltanieh, M., and Gill, W.N., "Review of Reverse Osmosis Membranes and Transport Models", Chem. Eng. Commun. 1981; 12, 279-363.

    [16] Lonsedale, H.K., Merten, U., and Reily, R.L., “Transport Properties of Cellulose Acetate Osmosis Membranes”, J.Appl. Poly.Sci, 1965; 9, 1341-1362.

     

    [17] Banks, W., and Sharples, A., "Studies on Desalination by Reverse Osmosis: Mechanism of Solute Rejection", J.Appl. Chem., 1966; 16, 153-158.

     

    [18] Sourirajan, s., “Reverse Osmosis”, Academic Press, New York, NY, 1970.

     

    [19] Spiegler, K.S., and Kedem, O., “Thermodynamics of Hyperfiltration (Reverse Osmosis): Criteria for Efficient Membranes”, Desalination, Vol. 1966; 1, 311-326.

     

    [20] Sherwood, T.K., Brian, P. L. T., and Fisher, R. E., “Desalination by Reverse Osmosis”, I & EC Fund, 1967; 6(1), 2-12.

    [21] Jonsson, G., and Boesen, C.E., "Water and Solute Transport through Cellulose Acetate Reverse Osmosis Membranes", Desalination, 1975; 17, 145-165.

     

    [22] Kimura, S., and Sourirajan, S., "Analysis of Data in Reverse Osmosis With Porous Cellulose Acetate Membranes Used". , ALCHE J. 1967; Vol. 13, 497.

     

    [23] Sourirajan, S., and Matsuura, T., “Reverse Osmosis Ultrafiltration Process Principles”, NRC, Ottawa, Canada, 1985.

     

    [24] Mehdizadeh, H. and Dickson, J. M., “Theoretical Modification of the Finely Porous Model for Reverse Osmosis Transport”, J. Appl. Polym. Sci., 1991; 42, 1143-1154.

    [25] Matsuura, T., and sourirajan, S.

Predicting the performance of reverse osmosis membranes using optimization, mathematical modeling and solving the model with the help of numerical methods