This paper considers the problem of identifying unknown parameters for a mathematical model of an ion-exchange filter via measurement at the outlet of the filter. The proposed mathematical model consists of a material balance equation, an equation describing the kinetics of ion-exchange for the nonequilibrium case, and an equation for the ion-exchange isotherm. The material balance equation includes a nonlinear term that depends on the kinetics of ion-exchange and several parameters. First, a numerical solution of the direct problem, the calculation of the impurities concentration at the outlet of the filter, is provided. Then, the inverse problem, finding the parameters of the ion-exchange process in nonequilibrium conditions, is formulated. A method for determining the approximate values of these parameters from the impurities concentration measured at the outlet of the filter is proposed.
1. Introduction and Background
One of the key technologies for the preparation of water at thermal power plants and nuclear power plants is the use of ion-exchange materials. Water treatment using ion-exchange methods is based on the passing of source water through a filter bed of ion-exchange material. That material is substantially insoluble in water but capable of ion-exchange with ions contained in the treated water. Water is a good solvent of various salts that have a negative solubility coefficient; that is, when the water temperature rises, the amount of dissolved salts is reduced, and the excess amounts of these salts are precipitated in heated tubes, forming scale with a very low coefficient of thermal conductivity. Power-generating units of thermal power plants (TPP) operate under severely corrosive conditions: high temperature (515–530°C) and pressure (15 MPa) hot steam [1–3]. The main reason for failures of power equipment is failure of the most thermally loaded parts, namely, pipes in which the coolant water flows. Scale restricts flow of heat energy to the coolant (water) and does not allow for the cooling of pipe surfaces; this causes the surface metal temperature to increase, reaching the critical melting point. The result is an emergency stop of the boiler [3]. Another problem is that the rate of corrosion of the metal under scale on the surface of the tubes sharply increases; that is, scale is a corrosion catalyst [4–7]. Reducing metal corrosion requires not only using materials with high corrosion resistance but also taking all measures to reduce aggressiveness of the coolant medium, water [4, 7, 8]. The primary source of harmful impurities is water from the water treatment plant, which is used to compensate for the loss of water coolant in the power plant cycle. The necessary water for the plant production process is usually taken from rivers and lakes.
This technology requires periodic regeneration [9, 10] of the filters. Regeneration of ion-exchange filters results in highly mineralized, acidic, and alkaline waste water [10]. After circulating through the energy plant, it is a highly mineralized aqueous solution. The high volume of water used produces an increase in the salinity of rivers and lakes when it is returned and causes harmful environmental impact due to the deterioration of the quality of the water. Approximately 25–30% of the water purifier output is for water treatment facility needs. For example, given a thermal power plant with 1200 MW of power, the water treatment plant capacity will be approximately 250 m3/hour; that is, an average of approximately 62.5 m3/hour or 547,500 m3/year of water is used for filter recovery. The switching from filtering (water treatment) to backwash, regeneration, and washing modes is carried out either by a time-scheduled mode or based on the volume of passed water and solution. Any of these procedures leads to the consumption of more water and reagents than needed. Thus, it is necessary to find alternative filter operation rules for the aim of reducing the water and reagents consumption. The problem of water treatment for thermal power plants using the ion-exchange method is a multiparametric task. The existing techniques for designing the operation of the water purification plant exclude the possibility of reaching an optimal mode of operation [9, 11, 12].
In this paper we propose a new method based on analysis of the water quality at the filter outlet. That is, we model the dynamics of the ion-exchange process at the filter, and then we solve an inverse problem of finding values for the process parameters. In other words, one can make decisions based on the real work conditions of the equipment. The operation of the water treatment plant under those operating rules can increase the operating cycle of the ion-exchange filter. The amount of waste water and water consumed for its own needs can be reduced. The method is easy to implement for water treatment equipment, either in an existing TPP or at the stage of designing a water treatment plant’s equipment. That allows one to predict the stopping time of the filtration process and the starting time for a regeneration process.
Section 2 of the paper contains a description of the mathematical model of the ion-exchange filtration process. In Section 3, a description of the method to solve the direct problem is provided. Section 4 describes the statement and the main idea of solving the inverse problem. In Section 5, the method for the inverse problem with one unknown parameter is developed. In Section 6, results of the numerical solution of the inverse problem with two unknown parameters of the model are provided. The final section contains a brief summary of the research conducted.
2. Mathematical Modeling of the Ion-Exchange Filtration Process
Several mathematical models have been proposed to describe the dynamics of ion-exchange processes [13–16]. In general, these models consist of three equations: mass balance equation, kinetic equation, and isotherm equation. The ion-exchange kinetics equation determines the speed of delivery of ions to the surface of the ion-exchanger. Since the system comprises a kinetic equation, this system of equations models an ion-exchange filter operating in nonequilibrium conditions. In the following, we describe each of them as presented in [11].
2.1. Mass Balance Equation
Sorption of an ions mixture by an ion-exchange material is subject to the law of mass conservation. Thus, the mathematical model for the ion-exchange in the filter also includes a mass balance equation:(1)wap∂C∂x+ε∂C∂t+∂q∂t=0,x∈0,L,t>0,where
wap is the flow velocity in the free section of the filter in m/sec;
ε is the porosity of the ion-exchange material in m3/m3;
∂C/∂x is the rate of change of the ion concentration in the solution, C(x,t), along the filter axis (x) in (mEq/kg)/m;
∂C/∂t is the rate of change of the ion concentration in the solution, C(x,t), with respect to time (t) in (mEq/kg)/sec;
∂q/∂t is the rate of change of the ion concentration in the ion-exchange material, q(x,t), at the point x with respect to time (t) in (mEq/kg)/sec;
L is the thickness of the filter material in m.
2.2. Kinetic Equation
Processes of the anion exchange and cation exchange of fresh natural waters, as well as processes of ion-exchange of condensed water, run in accordance with the mechanism of the external diffusion kinetics. One of the equations that is customarily used to describe the rate of ion-exchange in the external diffusion kinetics is the Nernst kinetics equation [9]. However, the Nernst theory is not always convenient for describing the diffusion kinetics because of the large number of experimental variables in the equation. To overcome this difficulty, a simpler alternative is proposed in [9]. An equation from the general theory of mass transfer is used, which expresses the hypothesis that the rate of ion-exchange is proportional to the deviation of the ion concentration in the sorbent from the corresponding equilibrium value. The equation is mathematically written in the following form:(2)dqdt=ηC-C¯,where
η is the mass transfer coefficient (proportionality) in m3/m3;
q is the concentration of exchanging ions in the sorbent in mEq/kg;
C is the concentration of the exchanging ions in the solution in mEq/kg;
C¯ is the equilibrium ion concentration with respect to the value q in the solution in mEq/kg, which borders with the surface of the ion-exchange material, in the filtration case C>C¯.
2.3. Exchange Isotherm Equation
The adsorption isotherm equation is used to determine the capacity of each sorbent for various sorbates. In the ion-exchange theory, isotherms of Langmuir [17], Freundlich [18], and Redlich-Peterson [19] are commonly used. These studies have shown that, for a mathematical description of the ion-exchange process in nonequilibrium conditions, the Langmuir isotherm equation described in [9] may be used. The corresponding mathematical expression can be written as follows:(3)φ¯=1-φaθk+1-kθ.Here, we get that
φ¯=C¯/C0 is the relative concentration of exchanging ions in the solution phase in mEq/kg, where C0 is the maximal possible concentration of exchanging ions in the original solution and k is the constant of the exchange isotherm (k>1 for a filtration process and k<1 in a regeneration case);
φa is the relative concentration of additional ions in the ion-exchange material. These ions form a weakly dissociated compound. φa=Ca/C0 where Ca is the concentration of ions in the original solution in mEq/kg that form a weakly dissociated compound;
θ is the relative concentration of exchanging ions in the ion-exchange phase in mEq/kg. θ=q/q0, where q0 is the total exchange capacity of the ion-exchange material.
Thus, taking into account (1)–(3), the mathematical model of the ion-exchange filtration process for nonequilibrium conditions can be written as follows:(4)ε∂C∂t+wap∂C∂x+∂q∂t=0,∂q∂t=ηC-C¯,φ¯=1-φaθk+1-kθ.Because exchange isotherm equation (3) is expressed in terms of relative concentrations of exchanging ions in the solution phase and the ion-exchange phase, we can write the equation of balance and kinetics in terms of relative concentrations. By dividing the equation of (1) by the value of C0 and the equation of (2) by the value of q0, we rewrite (4) in terms of the functions φ(x,t)=C(x,t)/C0 and θ=q(x,t)/q0, given that C¯=C0φ¯:(5)∂φ∂t+wapε∂φ∂x+q0C0ε∂θ∂t=0;dθ∂t=ηC0q0φ-φ¯;φ¯=1-φaθk+1-kθ.Let us introduce the following notations for the dimensionless constants appearing in the system of equations:(6)c=wapε;β=q0C0ε;γ=ηC0q0;m=1-φa.0<m<1 is a constant dealing with possible chemical reactions with impurities of other ions in solution.
Given the notations in (6), the system of equations in (5) can be rewritten as(7)∂φ∂t+c∂φ∂x+β∂θ∂t=0,∂θ∂t=γφ-fθ,t∈0,T,x∈0,L,φ¯=fθ≡mθk+1-kθ.Here, we use the function f(θ) to denote the relationship between the relative concentration of exchanging ions in the solution phase φ¯ and the relative concentration of exchanging ions in the ion-exchanger phase θ:(8)φ¯θ=fθ≡mθk+1-kθ.The result is the system of equations in (7) that describes the model of the ion-exchange filter operating under nonequilibrium conditions. The following are boundary conditions for system (7):(9)φ0,t=φ0t,0≤φ0t<1,φx,0=0,θx,0=0.The physical meaning of these terms is as follows. The condition φ0,t=φ0t indicates that, over the ion-exchanger layer (x=0), at any time t, there is a solute with constant concentration φ0(t). θ(x,0)=0 shows that initially (t=0) the absorbed ion is absent at any level of layer x of the ion-exchanger. The equation φ(x,0)=0 indicates that, at the initial time (t=0), the concentration of ions absorbed from the solution in the filter material is negligibly small.
Calculations showed that a change in porosity of the ion-exchange material ε virtually had no impact on the period of time before reaching the permissible concentration of the ion to be removed at the output of the filter, so the average value of 0,325 is used for calculation.
3. The Numerical Solution of the Direct Problem
Direct problem (7), (9) with known constants k and m is nonlinear; however, the nonlinear function f(θ) satisfies the Lipschitz condition, so it can be shown that the solution of problem (7), (9) exists and is unique. The proof of the uniqueness and existence of the solution to the direct problem is out of the scope of this study. For the numerical solution of problem (7), (9) we apply the finite-difference method. According to this method, the partial derivatives in the system of equations of (7), (9) are replaced with their difference analogues:(10)φin+1-φinτ+cφin+1-φi-1n+1h+βγφin-fθin=0,θin+1-θinτ=γφin-fθin,t∈0,T,x∈0,L,fθin=mθink+1-kθin,(11)φi0=0,θi0=0,i=1,2,…,N,φ0n=φ0=const,n=0,1,…,NT.Here, τ=T/NT is the time step, and h=L/N is the step in space on a uniform grid.
The first equation of the direct problem is approximated by the implicit difference scheme of point-to-point computation, since when using the calculation parameters characteristic of the technological process the explicit scheme appeared to require a very small step equal to, in terms of time, τ = 0,000014 for N=128. Due to the nonlinear nature of the problem, in order to ensure consistency of the calculation, the time step to be used for the point-to-point computation must be equal to τ = 0,000125 for N=128. The deviation in the results of calculations using the explicit scheme and the point-to-point computation scheme was calculated using the following formula:(12)δ≜maxjφ1imtj-φ1extjavgφ1imtj,where φ1im(tj) is the relative concentration of Na+ cation to be removed at the outlet of the filter when calculating using the implicit finite-difference scheme; φ1extj is the relative concentration of Na+ cation to be removed at the outlet of the filter when calculating using the explicit finite-difference scheme.
The deviation was equal to 0,075%. To investigate the influence of parameters k and m on the output data the grids of parameters k and m are given at their characteristic ranges of variation. Then, for each pair of discrete values of k and m, the concentration function φ1(t) is calculated at the output of the filter on a grid of time values.
We compared the results with data from a real case. The data were collected at the 1200 MW thermal power plant located in the Pavlodar Region of the Republic of Kazakhstan. At the thermal power plant, a 3-stage scheme of deep water desalination is provided, and the capacity of the main water treatment plant is 300 t/h.
Calculations are made for the values wap=5.55×10-3 m/sec, ε=0.325 m3/m3, β=2.5×103, L=2.5 m, and C0/q0=2.5×10-3, which are typical for the hydrogen-cation exchange filters used at thermal power plants in the first stage of water treatment.
In Figure 1, we show the characteristic graph of φ depending on t at filter layer x = 2,5 m.
The dependence of φ on t.
Figure 1 shows changes in the relative concentration of Na+ cations in time during operation of the filter in the water filtration mode at the output of the filter according to experimental and calculation data. According to the calculations, the relative concentration of Na+ cations during the first 40 hours of the filter operation is extremely low, ranging from 0.000178 to 0.000477; then breakthrough of Na+ cations is observed which leads to an increase in the relative concentration of Na+ cations at the output from the filter. The filtering process ends when the maximum of the relative concentration of Na+ cations at the output of the filter is achieved equal to 0.297, with the filtration duration being 55 hours. According to the experiment data, the maximum relative concentration of cation removed is achieved at the moment in time being 55 hours. Comparison of the graphs shows that the mathematical model approximation (7), (9) using implicit finite-difference scheme (10) provides satisfactory results.
In Figure 2, we show the characteristic graph of φ depending on x.
The dependence of φ on x.
Figure 2 shows changes in the relative concentration of Na+ cations depending on the height of the layer of the ion-exchange material in the filter based on calculation data 20 hours after the beginning of the filtration process. When filtering water through the ion-exchanger layer, the ion-exchange process flows layer by layer. The total height of the ion-exchange layer is 2.5 m. The graph shows three zones as follows: 0 to 0.3 m high layers of depleted ion-exchanger, which cease to be part of the ion-exchange processes (I); 0.3 to 1.3 m high working ion-exchanger layers in which the ion-exchange processes take place (II); and 1.3 to 2.5 m high layers of fresh ion-exchanger which are not yet involved in the ion-exchange processes (III). As the depletion progresses the upper layers of the ion-exchanger cease to be part of the ion-exchange processes. Instead of them fresh ion-exchanger layers are involved in operation which lie under the working layer, and the ion-exchange zone progressively moves as the water treatment process continues.
4. Formulation of the Inverse Problem
During the operation process, the measurement of the impurities concentration at the outlet of the filter is available:(13)φL,t=φ1t,0<t<T.Here, T is the horizon observation time.
Due to the complexity of the sorption process and heterogeneity of the mixture to the filter, the parameters γ in model (7) and k and m in (8) cannot be determined accurately. In practice, we can operate with approximate values of k, m, and γ. In some cases, when the main ions of the filtered mixture are known, we can assume that the value of k is given. Then, we still have two uncertain model parameters, the exchange isotherm parameter m and parameter γ. Having the measurements of (13), the problem is actually reduced to an inverse problem, which can be stated as follows:
Find the model parameters m and γ of problem (7), (9) via boundary measurements (13).
Although there are examples of formulations and solutions of inverse problems for filtration processes in the literature [20], this statement of the inverse problem is new. The standard method to solve inverse problems is quasisolution [21], according to which the unknown parameters of the model are determined by minimizing the residual function. Earlier, in [22], we investigated the behavior of the residual function of the following form: (14)Jk,m=12∫0TφL,t;k,m-φL,t;k0,m02dtwith the fixed value of the mass transfer coefficient γ. On the basis of the direct numerical simulation described in [23], we have shown that function (14) is not convex. Additionally, the low sensitivity of function (14) with respect to the model parameters k and m is established. This means that the routine use of gradient methods is inefficient. As a result, we propose an alternative approach to the solution of the formulated inverse problem. The main idea of the method is the following; as we have two unknown parameters of the model, m and γ, it is necessary to find at least two integral characteristics of measured data (13), which clearly define these parameters.
5. Solution of the Inverse Problem with One Unknown Parameter, m, of the Model
For solving the inverse problem of finding two unknown parameters m and γ first consider the solution for a simple problem. We solve the problem with one unknown parameter m of the model. Assume first that the parameters k and γ in model (7), (9) are known and state the problem of determining the value of m. The natural measured value is the total amount of impurity Q passed through during the ion observation period:(15)Qk,m,γ=∫0Tφ1tdt.Due to the extremely large scatter of Q(k,γ,m), it is more convenient to use its logarithm, I1=logQ. It is established numerically that the function φ1(t) and the value I1 depend on the ratio s=k/m. This fact indicates that, using measured data (13), it is possible to recover the ratio s=k/m, and it might be difficult to recover the variables k and m simultaneously. It follows that the models with the same value of the function s(k,m)=k/m are equivalent with respect to measured data (13) and observed value of I1=logQ. In Figures 3 and 4, we show the characteristic graphs of logQ depending on s and m.
The dependence of logQ on the parameter s.
The dependence of logQ on the parameter m.
Thus, the following equivalence principle is obtained for the problem of filtration (7)–(9): mathematical models with different m and k and the same s=k/m are equivalent with respect to the observed value I1 and measured data (13).
This result seems to be natural if we expand the function f(θ) with respect to θ as follows:(16)fθ=mθk+1-kθ=mkθ+mk-1k2θ2+⋯.One can see that the function f(θ) can be approximated as f(θ)=θ/s and depends on the ratio s=k/m for small values of θ. Note that work parameters of the filter use small values of θ, which explains why we obtain the same output data with the same value of the ratio s for different k and m.
Calculations show monotonic dependence I1(s) of s in the practically acceptable range of k and m. This implies that if we know the value of k, then, on the basis of the observed value of I1, the value of s=k/m is uniquely determined, and, through it, in turn, the value of m=k/s can be calculated. Thus, we have two uncertain parameters of the model, the parameters γ and s=k/m, which is associated with exchange isotherm (3). For practical applications of this result, we have generated a table of values of I1=logQ for a range of values of s in the interval [1.0,5.0]. It is enough to interpolate the data of the generated table to recover the desired value of s0 using the measured value of logQ0. Table 1 shows the values of the exact (sexact) and restored (srest) parameter s values for hydrogen-cation exchange filters used in the first stage of make-up water treatment.
The exact (sexact) and restored (srest) parameter s values.
Parameters
Parameter values
sexact
1.45
1.7
2.0
2.45
3.1
3.5
4.0
4.45
4.7
srest
1.453
1.71
2.14
2.453
3.12
3.52
4.14
4.451
4.71
δs~
0.2
0.58
7
0.12
0.65
0.57
3.5
0.02
0.21
δs~: relative calculation error s, %.
The test results showed satisfactory accuracy of s variable restoration.
6. Restoring Parameters m and γ
Suppose now that we know parameter k and do not know two parameters, m and γ. To calculate the range of change of γ, we can use definition (6) of γ:(17)γ=η·C0q0.In practice, parameter η is calculated according to [9, 24] using the empirical formula(18)η≈0.009αwap0.53d1.47,where wap is the rate of flow of fluid in the unit in m/s, d is the ion-exchange grain size in m, and α is an empirical coefficient.
For the calculation, we use the parameter values that are typical for the hydrogen-cation exchange unit of the first stage of water treatment, wap = 0,0056 m/s, and the value of d for the ion-exchange material varies from 0.315 to 1.250 mm. The coefficient α for the exchange of H+ ions for Na+ is 1. Based on this data, the values of γ in the processes of interest are in the range of 1.523×103≤γ≤3.165×103. Suppose that another observable quantity that is sensitive to the values of the unknown parameters of the model is the logarithm of the impurities concentration at the outlet of the filter at time t=T:(19)I2k,m,γ=logφ1T.By solving the direct problem in the range of variation of parameters s and γ, we have generated tables of values of data for I1s,γ and I2(s,γ). Based on these two tables, we can calculate the appropriate values of s0 and γ0 from the measured values of I1,0 and I2,0. Below is described a numerical method for calculating s0 and γ0. Thus, suppose that we have a table of values for the two measured variables at the filter outlet, I1(si,γj) and I2(si,γj). Suppose that I1,0 and I2,0 are measured values. Using the tables, we can determine a pair of values, sk, γl, for which the minimum deviation is reached:(20)I1sk,γl-I1,02+I2sk,γl-I2,02=mini,jI1si,γj-I1,02+I2si,γj-I2,02.Assuming that the quantities I1(s,γ) and I2(s,γ) are smooth functions of their arguments, we can write approximate equations at a point of sk, γl:(21)I1,0=I1s0,γ0=I1sk,γl+∂I1∂ssk,γls0-sk+∂I1∂γsk,γlγ0-γl,I2,0=I2s0,γ0=I2sk,γl+dI2dssk,γls0-sk+∂I2∂γsk,γlγ0-γl.The partial derivatives ∂I1sk,γl/∂s, ∂I1sk,γl/∂γ, ∂I2sk,γl/∂s, and ∂I2sk,γl/∂γ are approximately calculated using the tables via formulas of finite differences. The relations in (21) represent a linear system of equations with indeterminate s0, γ0. This system is solved numerically and defines approximate values of the unknown parameters. Table 2 shows the values of the exact (γexact, sexact) and restored (γrest, srest) parameters γ and s for the hydrogen-cation exchange filters used in the first stage of make-up water treatment.
The values of the exact and recovered parameters γ and m.
The test results show satisfactory accuracy of m and γ variable recovery.
The developed method of recovering the parameters was verified with errors ±1%, ±3%, ±5%, and ±7% entered in the exact values Q, φ1(T). The mean relative errors of s and γ parameters recovery are shown in Table 3.
The values of the mean relative errors δs and δγ depending on different measurement errors ε.
S_exact
γ_exact, 10−3
ε=1%
ε=3%
ε=5%
ε=7%
δs (%)
δγ (%)
δs (%)
δγ (%)
δs (%)
δγ (%)
δs (%)
δγ (%)
1.6
1.6
0.0
0.6
0.0
0.8
0.0
1.9
0.0
2.3
1.6
2.0
0.0
0.8
0.0
1.6
0.2
2.9
0.0
3.2
2.0
1.6
0.5
2.1
1.0
3.1
1.0
3.1
2.8
3.6
2.0
2.0
0.3
1.8
1.0
1.9
2.3
13.5
2.8
13.5
2.6
2.0
0.4
14.6
1.3
14.6
2.2
12.2
3.0
6.7
7. Conclusion and Discussions
In the paper a mathematical model of the ion-exchange filter operation mode is considered (7), (9). The direct problem for the ion-exchange filter model was numerically solved. For calculation it is necessary to know the values of m and γ, which in practice are not available for measurement. In the paper a method for the determination of these coefficients based on the output values of ion absorbed at the outlet of the filter was proposed; the inverse problem was formulated. The method is based on the definition of the integral characteristics of the input data that accurately describe these parameters. The integral characteristics were determined. The integral characteristics are the logarithm of the number of impurity ions passing through during the observation period and the increase of the logarithm of the concentration of impurities at the outlet of the filter. The inverse problem was numerically solved given the unknown m and γ parameters. Resetting the parameters using the proposed method yielded satisfactory results.
The developed method can be used in simulation of the ion-exchange unit operation mode. This will allow determining the optimal mode of operation facilitating reduction of the amount of regeneration and washout water.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work is supported by the Ministry of Education and Science of the Republic of Kazakhstan under Grants nos. 0112RK02067 and 0113RK00442.
FortuneH. J.Boiler feed-water treatment for corrosion control—part 219629615415610.1108/eb019971MesherskiyN. A.Moscow, RussiaEnergoatomizdat1984 (Russian)SidelkovskyL. N.YurenevV. N.19883rdEnergoatomizdat(Russian)AvlonitisS. A.PappasM.MoutesidisK.AvlonitisD.KouroumbasK.VlachakisN.PC based SCADA system and additional safety measures for small desalination plants200416516517610.1016/j.desal.2004.06.0192-s2.0-8644290765ChoudhuryM. R.HsiehM.-K.VidicR. D.DzombakD. A.Corrosion management in power plant cooling systems using tertiary-treated municipal wastewater as makeup water20126123124110.1016/j.corsci.2012.04.0422-s2.0-84861924967PhongphiphatA.RyuC.YangY. B.FinneyK. N.LeylandA.SharifiV. N.SwithenbankJ.Investigation into high-temperature corrosion in a large-scale municipal waste-to-energy plant201052123861387410.1016/j.corsci.2010.07.0322-s2.0-77957665384SlobodanV.NikolaV.ŽeljkoD.Power electronics solution to dust emissions from thermal power plants20107223125210.2298/sjee1002231vPrisyazhniukV. A.Alternative trends in development of thermal power plants2008282-319019410.1016/j.applthermaleng.2007.03.0252-s2.0-35348900042GromoglasovA. A.KopylovA. C.Pil'shikovA. V.1990Energoatomizdat(Russian)GlazyrinA. I.MuzykaL. P.KabdualievaM. M.Almaty, KazakhstanRepublican Publishing Office1997 (Russian)KostrikinY. M.MisherskyN. A.KorovinO. V.Moscow, RussiaEnergoatomizdat1990 (Russian)LivshicO. V.Moscow, RussiaEnergy1976 (Russian)BogatyryovV. L.Clathrate-forming ion exchangers199816122326510.1080/073662998089345302-s2.0-0031699177BolančaT.Cerjan-StefanovićŠ.Optimization strategies in ion chromatography2007305–779180610.1080/108260707011910782-s2.0-33947116332DavisJ. A.KentD. B.Surface complexation modeling in aqueous geochemistry1990231177260WangZ.GiammarD. E.Mass action expressions for bidentate adsorption in surface complexation modeling: theory and practice20134793982399610.1021/es305180e2-s2.0-84877598404LangmuirI.The constitution and fundamental properties of solid and liquids—part 1 solids19184013611403FreundlichH.Uber die adsorption in Losungen190657385471JossensL.PrausnitzJ. M.FritzW.SchlünderE. U.MyersA. L.Thermodynamics of multi-solute adsorption from dilute aqueous solutions19783381097110610.1016/0009-2509(78)85015-52-s2.0-0017918391TokmachevM. G.2008(Russian)TikhonovA.ArseninV.1977New York, NY, USAJohn Wiley & SonsMR0455365MukanovaB. G.GlazyrinaN. S.Numerical modeling of the ion-exchange filter work201412026KamkeE.1944LeipzigZagrayY. M.KoganovskyA. M.KulskiyL. A.Kiev, UkrainePublishing Ukr NIITI1966 (Russian)