Evaluation of Transport Properties and Energy Conversion of Bacterial Cellulose Membrane Using Peusner Network Thermodynamics

We evaluated the transport properties of a bacterial cellulose (BC) membrane for aqueous ethanol solutions. Using the Rr version of the Kedem–Katchalsky–Peusner formalism (KKP) for the concentration polarization (CP) conditions of solutions, the osmotic and diffusion fluxes as well as the membrane transport parameters were determined, such as the hydraulic permeability (Lp), reflection (σ), and solute permeability (ω). We used these parameters and the Peusner (Rijr) coefficients resulting from the KKP equations to assess the transport properties of the membrane based on the calculated dependence of the concentration coefficients: the resistance, coupling, and energy conversion efficiency for aqueous ethanol solutions. The transport properties of the membrane depended on the hydrodynamic conditions of the osmotic diffusion transport. The resistance coefficients R11r, R22r, and Rdetr were positive and higher, and the R12r coefficient was negative and lower under CP conditions (higher in convective than nonconvective states). The energy conversion was evaluated and fluxes were calculated for the U-, F-, and S-energy. It was found that the energy conversion was greater and the S-energy and F-energy were lower under CP conditions. The convection effect was negative, which means that convection movements were directed vertically upwards. Understanding the membrane transport properties and mechanisms could help to develop and improve the membrane technologies and techniques used in medicine and in water and wastewater treatment processes.


Introduction
Membranes exert multiple functions, including protective, regulatory, and coordinating functions [1]. The protective and regulatory functions are based on the membrane selectivity of a barrier that regulates the transport between the interior of a system and its surroundings [2]. The separation properties of synthetic polymer membranes enable their application in many fields of science, technology, and medicine, such as food production, water treatment, hemodialysis, wastewater treatment, and membrane dressings [1,3]. The coordinating function seems to apply only to biological systems, where the membrane simultaneously plays the role of a receiver, regulator, and coordinator of environmental signals, which are the driving forces of membrane transport [4]. These driving forces cause various types of physical fields, such as concentration, pressure, temperature, or electric potentials, which participate in shaping the field character of nature [5,6].
Selective permeability is one of the basic properties of porous media, including polymeric membrane-forming materials. It is required for large-scale models of fluid flow and mass transport. These models operate within the framework of nonequilibrium thermodynamics, hydrodynamics, and statistical physics [7]. The last 20 years have seen a significant increase in the use of modeling to study multiphase flow and transport in porous media. Starting with models of fluid systems in single pores, calculations of the relative permeability, interfacial area, dissolution rate, and many other physical properties have been carried out [8]. One of the more interesting methods for pore-scale numerical studies is direct hydrodynamic simulation (DHD) technology [9], which uses a description of the dispersed interface and is applied to various fluid-rock or fluid-fluid interactions for equal rheological conditions.
In the paper [10], a stochastic method based on simulated annealing and X-ray microtomography was used to study the pore structures of various porous solids that differ in pore space morphology and topology. In the process of verifying the developed models, many interesting simulation and experimental results were obtained, which confirmed the pore space models. In addition, it was shown that predictions based on tomographic pore space models were more effective than stochastic models and that the time-dependent effective diffusivity is particularly sensitive to small morphological deviations between the actual and reconstructed pore structures. It was also shown that the combined prediction of the effective permeability, effective pore size, geometric coefficient, and time-dependent effective water diffusivity is needed to reliably evaluate pore space reconstruction.
The paper [11] presented a description of volume membrane transport using the Kedem-Katchalsky equations of homogeneous aqueous solutions of ethanol and glucose. The flows generated by the hydrostatic pressure differential, the osmotic pressure differential, and the simultaneous action of these two thermodynamic drives were analyzed independently. In addition, a formula for the membrane filtration coefficient was presented, taking into account the density and viscosity of ethanol, and the corresponding calculations were made. In this way, it was shown that the membrane filtration coefficient depends on both the membrane properties and the flowing fluid. In turn, the paper [12] presented formulas for calculating the viscosity and diffusion coefficients of binary aqueous nonelectrolyte solutions as a function of the solution concentration under isothermal conditions. In the process of verifying the obtained formulas, the dependence of the diffusion coefficient, filtration coefficient, and dynamic viscosity coefficient for these solutions on the solution concentration was calculated. Based on the obtained formulas, the results of calculations of the diffusion coefficient, membrane filtration coefficient, and dynamic viscosity of aqueous solutions of ethanol and glucose were presented. Processes such as diffusion or osmosis can modify physical fields, including concentration fields. In addition, a concentration field can be modified by the concentration polarization (CP) as a consequence of the creation of the concentration boundary layers (CBLs) l r h and l r l on both sides of the membrane [5,13]. The thickness of CBL l r h is δ r h , and the thickness of CBL l r l is δ r l . As a consequence of the CBL formation, the concentration difference decreases from the value of C h -C l to the value of C r h − C r l , where C r h > C r l , C h > C r h , and C r l > C l , and the density difference increases from ρ h − ρ l to the value of ρ r h − ρ r l , where ρ r h > ρ r l , ρ h > ρ r h , and ρ r l > ρ l . When a lower density solution is placed in the compartment under the membrane and a higher density solution is placed in the compartment above the membrane, the complex l r h /M/l r l loses its hydrodynamic stability. Hydrodynamic instability is manifested by natural convection in near-membrane areas [13][14][15].
Then, the concentration Rayleigh number (R C ), which controls the process of the appearance of gravitational convection, exceeds its critical value, and hydrodynamic instabilities appear in the near-membrane areas [5,14,[16][17][18]. Over time, the destructive effect of gravitational convection limits the growth of δ r h and δ r l and accelerates the diffusion of substances beyond the layers, which extends the effect of convection to the entire volume of the solution. Under certain conditions, even self-organization of the liquid may occur, which is manifested in the "plum structure" [19]. The creation of CBLs can be visualized by a Mach-Zehnder laser interferometer [14]. The consequence of CP is a significant reduction in concentration gradients, as evidenced by the minimization of the osmotic and diffusion fluxes of dissolved substances and the membrane potentials [6,13,15]. Under certain conditions, depending on the compositions of the solutions and the orientation of the artificial biomembrane in relation to the gravity vector, concentration gradients can be reconstructed by gravitational convection [6,13,20]. The basic research tools for describing membrane transport are the Kedem-Katchalsky (KK) equations, both the classical version [3,21] and the modified forms [13,16,22]. For the concentration polarization conditions of the solutions, the equations have the form where L P , σ, and ω are the hydraulic permeability, reflection, and solute permeability coefficients; ∆P = P h − P l and ∆π = RT∆C are the hydrostatic and osmotic pressure differences (RT is the product of the gas constant and temperature, ∆C = C h − C l , C h and C l are the solutes concentrations, and ∆C is the solution concentration difference); J r v and J r s are the volume and solute fluxes; σ is the reflection coefficient; ω is the permeability coefficient of the solute; is the average concentration of the solutes; and ζ r p , ζ r v , ζ r s , and ζ r a are, respectively, the hydraulic, osmotic, diffusive, and advective coefficients of the CP [23]. For dilute nonelectrolyte solutions, σ v = σ s . In contrast, for nondilute solutions, σ v = σ s [16].
From Equations (1) and (2), the phenomenological coefficients of homogeneous solutions (ζ r p = ζ r v = ζ r s = ζ r a = 1) are defined as follows: The tetrad of membrane transport parameters (L P , σ v , σ s , and ω) play the roles of proportionality coefficients. The coefficients ζ r p , ζ r v , ζ r s , and ζ r a play similar roles. The products L P ζ r p , σζ r v , ωζ r s , and σζ r a determine the transport properties of the membrane complex and the concentration boundary layers.
The values of the L P , σ, and ω coefficients, for isotropic and electrically neutral artificial membranes and for dilute solutions, are constant. Examples of such membranes are those made of regenerated cellulose (Nephprophan and Cuprophan) and bacterial cellulose (Biofill) [24][25][26]. The values of these coefficients for compound and ion-exchange membranes (Nafion and Textus bioactiv) are concentration-dependent [27][28][29][30].
The Kedem-Katchalsky-Peusner (KKP) equations are the network forms of the KK equations proposed by L. Peusner, which are obtained by means of the symmetrical or hybrid transformations proposed by Peusner network thermodynamics [25]. These equations contain the Peusner coefficients (R ij ), which for the conditions of homogeneity of solutions are a combination of the phenomenological coefficients of the membrane (L p , σ, and ω) and the average concentration of the solutions (C).
In previous papers [24,31], we showed descriptions of the membrane transport of binary solutions of nonelectrolytes under conditions of heterogeneity of solutions by introducing the R r forms of the KKP equations and the R versions of the KKP equations for binary solutions of nonelectrolytes for the conditions of solution homogeneity. Here, we evaluated the transport properties of a membrane for aqueous ethanol solutions and the conditions of CP using network KKP equations. We experimentally determined the time and concentration characteristics of the volume (J r v ) and solute (J r s ) fluxes for conditions of homogeneity and CP. Next, we calculated the time and concentration dependencies of the CP coefficients (ζ r v and ζ r s ) and resistance coefficients (R r ij and R r det ; i, j ∈ {1, 2}, r = A, B). We used J r v , J r s , and R r ij to calculate the energy conversion efficiency coefficients ((e r ij ) R ) and the flux of dissipated energy (S-energy) ((Φ r S ) R ). Then, we used (e r ij ) R and (Φ r S ) R to calculate the flux of free energy (F-energy) ((Φ r F ) R ) and the flux of internal energy (U-energy) ((Φ r U ) R ).

Membrane System
The system used as a model to study membrane transport, illustrated schematically in Figure 1, consisted of a membrane (M) situated in the horizontal plane and separating two aqueous solutions of ethanol with concentrations at the initial moment of C h and C l = constant (C h ≥ C l ). The density of the solutions with concentrations of C h and C l fulfilled the condition ρ h ≤ ρ l = constant. In configuration A, a solution with the concentration C l was located in the compartment above the membrane and a solution with the concentration C h was in the compartment under the membrane. In configuration B, the solutions with the concentrations of C l and C h were swapped. The study on volume ( ) and solute ( ) fluxes was carried out using the measuring set described in the paper [32] and is presented in Figure 2. It consisted of two cylindrical measuring vessels with volumes of 200 cm 3 each containing aqueous ethanol solutions, one with a concentration in the range of 1-501 mol m −3 and the other with a constant concentration of 1 mol⋅m −3 . The solutions in the vessels were separated by a previously described bacterial cellulose (BC) membrane called Bioprocess ® (Biofill Produtos Biotechnologicos S.A., Curitiba, Brasile) [33][34][35][36] positioned in a horizontal plane with an area of A = 3.36 cm 2 . The BC membrane was produced in flat sheets, and its structure was made of The study on volume (J r v ) and solute (J r s ) fluxes was carried out using the measuring set described in the paper [32] and is presented in Figure 2. It consisted of two cylindrical measuring vessels with volumes of 200 cm 3 each containing aqueous ethanol solutions, one with a concentration in the range of 1-501 mol m −3 and the other with a constant concentration of 1 mol·m −3 . The solutions in the vessels were separated by a previously described bacterial cellulose (BC) membrane called Bioprocess ® (Biofill Produtos Biotechnologicos S.A., Curitiba, Brasile) [33][34][35][36] positioned in a horizontal plane with an area of A = 3.36 cm 2 . The BC membrane was produced in flat sheets, and its structure was made of microcellulose fibers produced by Acetobacter Xylinum [8,37]. is the increase in the total concentration of the solutions. The was measured by a Rayleigh interferometer based on previously calculated feature curves, i.e., the experimental dependence of the shift of the interference bars (Δn) as a function of the ethanol concentration (C) [38]. The study was carried out at = 295 K. A laser interferometry method can also be used to determine [39][40][41]. We measured Δ and under intense mechanical stirring of the solutions at 500 rpm, and when steady-state flows were obtained, the stirring of the solutions was turned off. In the second step, the increments of Δ and were measured until steady-state flows were obtained. The volume flux was from a vessel with a lower concentration of solutions to a vessel with a higher concentration of solutions, and the solute flux was in the opposite direction. The where V u is the volume of the measuring vessel and dC r s is the increase in the total concentration of the solutions. The dC r s was measured by a Rayleigh interferometer based on previously calculated feature curves, i.e., the experimental dependence of the shift of the interference bars (∆n) as a function of the ethanol concentration (C) [38]. The study was carried out at T = 295 K. A laser interferometry method can also be used to determine dC r s [39][40][41]. We measured ∆V r and dC r s under intense mechanical stirring of the solutions at 500 rpm, and when steady-state flows were obtained, the stirring of the solutions was turned off. In the second step, the increments of ∆V r and dC r s were measured until steadystate flows were obtained. The volume flux was from a vessel with a lower concentration of solutions to a vessel with a higher concentration of solutions, and the solute flux was in the opposite direction. The ∆V r and dC r s were measured in a series of independent experiments. From the measurements of ∆V r and dC r s , the characteristics J r v = f (t) and J r s = f (t) were determined for different concentrations of ethanol solutions. For each characteristic, three independent experiments were performed. The relative error in the determination of J r v = f (t) and J r s = f (t) was no greater than 10%. Based on the time characteristics of J r v and J r s for the steady state, we calculated the concentration characteristics of J r v and J r s . Next, we used the characteristics J v = f (t) (for the homogeneity of solutions) and J r v = f (t) (for the conditions of CP) to calculate the dependence ζ r v = f (t). Similarly, we used the characteristics of J s = f (t) (for the homogeneity of solutions) and J r s = (t) (for CP conditions) to calculate the dependence ζ r s = f (t). Additionally, the dependence ζ r v = f (∆C) was determined based on the characteristics J v = f (∆C) (for the homogeneity conditions of solutions and J r v = f (∆C) for CP conditions. Similarly, the dependence ζ r s = f (∆C) was determined based on the characteristics J s = f (∆C) for the conditions of CP and J r s = f (∆C) for CP conditions. Moreover, the dependences R r ij = f (∆C),

The R r Form of Kedem-Katchalsky Equations for Binary Nonelectrolyte Solutions
For the interpretation of the obtained results, we used the R r form of the KKP equations, which can be obtained using simple algebraic transformations presented in the paper [24,31]: The above equations can be written in matrix form: where [R r ] is the matrix of resistance coefficients given by From the above equation, it follows that R r 12 = R r 21 and that the matrix determinant [R r ] is equal to To write Equations (3)- (7) for the homogeneity conditions of the solutions, it is enough to leave the 'r' index and assume ζ r p = ζ r v = ζ r s = ζ r a = 1. Then, we have R r ij = R ij and R r det = R det . The first part of the right-hand side of Equation (4) has the sense of the membrane Peclét number [42]. The classical definition of this number has the form Pe = (1 − σ)J v ℘ −1 = R 21 J v and appears in Equation (4) for conditions of homogeneity of solutions (ζ r s = ζ r a = 1, J r v = J v , and σ s = σ). In this equation (1 − σ)℘ −1 ≡ α is the Peclét coefficient and ℘ is the solute permeability coefficient. For conditions of concentration polarization, this number can be written in the form where ωRT = ℘ v and α r v is expressed in s m −1 . The second part of the right side of Equation (3) is similar to the membrane Peclét number. However, in this case, the definition of this number has the form Pe = (1 − σ)J s ℘ −1 = R 12 J s and appears in Equation (3) for conditions of homogeneity of solutions (ζ r s = ζ r v = 1, J r s = J s , and σ v = σ). For conditions of concentration polarization, this number can be written in the form where ωRTC = ℘ s and α r s is expressed in m 2 s mol −1 . It follows that the coefficients R 21 = a v RT, R 12 = a s RTC, R r 21 = α r v RT, and R r 12 = α r s RTC (r = A, B) are related to Peclét's coefficients, which are known from the literature [24,42]. Using the coefficients R r ij , R ij , R r det , and R det , it is possible to define the coefficients (φ r ij ) R and φ r det R , which are measures of the CP effect, and the coefficients (ϕ ij ) R and (ϕ det ) R , which are measures of the effect of gravitational convection in osmotic and diffusive membrane transport. The definitions of these coefficients can be written as The coefficients φ r ij and φ r det are measures of the distance of the membrane system from the CP state, and the coefficients ϕ ij and ϕ det are measures of the distance of the membrane system from the unstable state.
In thermodynamic systems, including membrane systems, U-energy can be converted into F-energy and S-energy (TS) [3,22]. If the solutions contain a solvent and one solute, the flux of S-energy for the CP conditions (Φ r S ) is described by the equation [5] ( where Φ r S R is the global S-energy for CP conditions, Φ r S R J r v is the S-energy produced by J r v , and Φ r S R J r s is the S-energy produced by J r s . Taking into account Equations (3) and (4), in Equation (14) we obtain An explicit form of the coefficients R r 11 , R r 12 , R r 21 , and R r 22 appearing in the above equation is given in Equation (6). To obtain a global S-energy for the conditions of homogeneity of solutions (Φ S ) in Equations (14) and (15), one should assume the condition of ζ r p = ζ r s = ζ r v = ζ r a = 1. In turn, from Equation (6), the coefficients (e r ij ) R , r r ij , and Q r R can be expressed using the coefficients R r ij . Using the definition proposed by Kedem and Caplan [43] and Peusner [30], we present the definitions of the energy conversion efficiency coefficients for CP conditions: is the coupling coefficient [30,43]. To obtain the expressions for r ij and r ji , it is enough to omit the superscripts "r" due to the fact that R r ij ≈ R r ji , R ij = R ji , r r ij = r r ji , and r ij = r ji . In turn, to obtain the expressions for (e ij ) R and (e ji ) R , it is enough to omit the superscripts "r" due to the fact that r r ij ≈ r r ji , r ij = r ji , (e r ij ) R ≈ (e r ji ) R , and (e ij ) R = (e ji ) R .
According to the first law of thermodynamics, for isothermal isochoric processes, the following equation is correct: where (Φ r S ) R = A −1 Td i S r /dt is the flux of dissipated energy (S-energy), (Φ r F ) R = A −1 dF r /dt is the flux of free energy (F-energy), and (Φ r U ) R = A −1 dU r /dt is the flux of internal energy (U-energy). All of these fluxes are expressed in Wm −2 .
We calculate the fluxes (Φ r F ) R and (Φ r S ) R using the expression below: Transforming this expression, we obtain The transport properties of the BC membrane were determined by the hydraulic permeability (L p ), reflection (σ), and solute permeability (ω) coefficients. The values of these coefficients, determined in a series of independent experiments that were carried out according to a previously described procedure [14], were L p = (62.8 ± 0.5) × 10 −12 m 3 N −1 s −1 , σ = (0.23 ± 0.01) × 10 −2 , and ω = (15.3 ± 0.5) × 10 −10 mol N −1 s −1 .  Figure 3c,d. The shapes of these graphs indicate that both J r v and J r s reached steady states relatively quickly and that in the steady states J A v > J B v and J A s > J B s . This dependence was a consequence of the emergence of gravitational convection, which is destructive to CBLs. This means that, in this case, CP and gravitational convection were antagonistic processes.

Results and Discussion
The characteristics of Figure 3a-d, illustrated by plots 2B, are typical for solutions whose densities decrease with increasing concentrations and CP conditions. Examples include aqueous solutions of ethanol, methanol, or ammonia [13]. If an aqueous solution of such a substance is placed in the compartment above the membrane (configuration B), a stable system of CBLs is formed, which causes a reduction in the value of the osmotic pressure difference, which results in the J B v and J B s fluxes. In configuration A, which refers to the situation when an aqueous solution of such a substance is placed in a compartment under the membrane, natural convection occurs, which decreases the reduction in the value of the osmotic pressure difference and causes fluxes J A v and J A s .   Figure 3a-d, illustrated by plots 2B, are typical for solutions whose densities decrease with increasing concentrations and CP conditions. Examples include aqueous solutions of ethanol, methanol, or ammonia [13]. If an aqueous solution of such a substance is placed in the compartment above the membrane (configuration B), a stable system of CBLs is formed, which causes a reduction in the value of the osmotic pressure difference, which results in the and fluxes. In configuration A, which refers to the situation when an aqueous solution of such a substance is placed in a compartment under the membrane, natural convection occurs, which decreases the reduction in the value of the osmotic pressure difference and causes fluxes and .

The Time and Concentration Dependencies of and
The time dependencies of and were calculated based on the results illustrated in Figure 4a As shown in Figure 4a, ζ r v and ζ r s take values from the intervals 0.03 ≤ ζ r v ≤ 1, 0.26 ≤ ζ r v ≤ 1, 0.06 ≤ ζ r s ≤ 1, and 0.29 ≤ ζ r s ≤ 1. Based on these time dependencies, the concentration dependencies ζ r v = f (∆C) and ζ r s = f (∆C) were determined for the steady states (Figure 4b). The coefficients ζ r v and ζ r s take their values from the . As shown in Figure 4b, ζ r v and ζ r s take their values in the range between 0.03 ≤ ζ r v ≤ 0.26 and 0.06 ≤ ζ r s ≤ 0.29. Therefore, the coefficients ζ r v and ζ r s are a measure of the CP in both convection and nonconvection states. tervals 0.03 ≤ ≤ 1, 0.26 ≤ ≤ 1, 0.06 ≤ ≤ 1, and 0.29 ≤ ≤ 1. Based on these time dependencies, the concentration dependencies = f(∆ ) and = f(∆ ) were determined for the steady states (Figure 4b). The coefficients and take their values from the intervals ( ) . ≤ ≤ ( ) . and ( ) . ≤ ≤ ( ) . . As shown in Figure 4b, and take their values in the range between 0.03 ≤ ≤ 0.26 and 0.06 ≤ ≤ 0.29. Therefore, the coefficients and are a measure of the CP in both convection and nonconvection states.  The transition from nonconvective to convective states is controlled by the Rayleigh concentration number ( ). The critical value of this number can be calculated from Equation [18]: where is the gravitational acceleration, D is the diffusion coefficient in the solution, 0 is the mass density, 0 is the kinematic viscosity of the solution, is the solute permeability coefficient through the membrane, / is the variation in density with concentration, and ζ is the concentration polarization coefficient.  (19), we obtain = −1155.07. The minus sign indicates that the convective currents are directed vertically upwards. In contrast, for aqueous glucose solutions, studied previously, convective currents were directed vertically downwards, and therefore had a positive sign [5]. The obtained critical value of is consistent with the values presented in the papers [44,45].

Concentration Dependencies of the Resistance Coefficients and
The concentration dependencies of the resistance coefficients 11 , 12 , 21 , 22 , and , calculated based on Equations (6) and (7), are shown in Figure 5a-d for homogenous solutions (curves 1A and 1B) and CP conditions (curves 2A and 2B). For all studied dependences of the resistance coefficients, curves 1A and 1B were symmetrical, whereas curves 2A and 2B were asymmetrical with respect to the point Δ = 0. The transition from nonconvective to convective states is controlled by the Rayleigh concentration number (R C ). The critical value of this number can be calculated from Equation (18): where g is the gravitational acceleration, D is the diffusion coefficient in the solution, ρ 0 is the mass density, ν 0 is the kinematic viscosity of the solution, ω is the solute permeability coefficient through the membrane, ∂ρ/∂C is the variation in density with concentration, and ζ is the concentration polarization coefficient.  (19), we obtain R C = −1155.07. The minus sign indicates that the convective currents are directed vertically upwards. In contrast, for aqueous glucose solutions, studied previously, convective currents were directed vertically downwards, and therefore R C had a positive sign [5]. The obtained critical value of R C is consistent with the values presented in the papers [44,45].

Concentration Dependencies of the Resistance Coefficients R r ij and R r det
The concentration dependencies of the resistance coefficients R r 11 , R r 12 , R r 21 , R r 22 , and R r det , calculated based on Equations (6) and (7), are shown in Figure 5a-d for homogenous solutions (curves 1A and 1B) and CP conditions (curves 2A and 2B). For all studied dependences of the resistance coefficients, curves 1A and 1B were symmetrical, whereas curves 2A and 2B were asymmetrical with respect to the point ∆C = 0.
The comparison of the dependencies R r 11 = f (∆C) (Figure 5a) for homogenous and CP conditions indicates that for |−∆C| = ∆C the condition R A 11 = R B 11 = R 11 was fulfilled. In turn, for the same values (−∆C), the condition R A 11 > R 11 was fulfilled, and for the same values of ∆C the condition R B 11 > R 11 . Moreover, for |−∆C| = ∆C the condition R A 11 < R B 11 was fulfilled, and for |−∆C| = ∆C = 62.5 mol m −3 the condition R A 11 = R B 11 was fulfilled. The dependencies R r ij = f (∆C) and R r ji = f (∆C), presented in Figure 5b,   The comparison of the dependencies 11 = f(Δ ) (Figure 5a) for homogenous and CP conditions indicates that for |−∆ | = ∆ the condition 11 = 11 = 11 was fulfilled. In turn, for the same values (−∆ ), the condition 11 > 11 was fulfilled, and for the same values of ∆ the condition 11 > 11 . Moreover, for |−∆ | = ∆ the condition 11 < 11 was fulfilled, and for |−∆ | = ∆ = 62.5 mol m −3 the condition 11 = 11 was fulfilled.
The dependencies = f(ΔC) and = f(Δ ), presented in Figure 5b, indicate that for |−∆ | = ∆ , the conditions were 12 = 21 , 12 = 21 , and 12 = 21 . In turn, for the same values (−∆ ), the condition 12 = 21 < 12 = 21 was fulfilled, and for the same values of ∆ the condition 12 = 21 < 12 = 21 was fulfilled. Moreover, for |−∆ | = ∆ the condition was 12 = 21 > 12 = 21 , and for |−∆ | = ∆ = 62.5 mol m −3 the condition was 12 = 21 = 12 = 21 . From Equation   Figure 5d, indicate that for |−∆C| = ∆C the condition R A det = R B det = R det was satisfied. In turn, for the same values (−∆C), the condition R A det > R det was fulfilled, and for the same values of ∆C-the condition R B det > R det . Moreover, for |−∆C| = ∆C the condition R A det < R B det was fulfilled, and for |−∆C| = ∆C = 62.5 mol m −3 the condition R A det = R B det was fulfilled. The coefficients R r 11 , R r 22 , and R r det were positive, and the coefficients R r 12 and R r 21 were negative and dependent on ∆C for the CP conditions. Compared to the conditions of homogeneity of solutions, the CP increased the value of the coefficients R r 11 , R r 22 , and R r det and reduced the value of the coefficients R r 12 and R r 21 . For the same ∆C, the values of the coefficients R r 11 , R r 22 , and R r det were smaller, and the coefficients R r 12 and R r 21 were higher for the nonconvective state.
Using the dependencies shown in Figures 3c,d and 5a 3.4. Concentration Dependencies (φ r ij ) R and (φ r det ) R Figure 6a,b shows the concentration dependencies of the CP effects (φ r ij ) R and (φ r det ) R , calculated based on Equations (10) and (11). Figure 6a shows that for the same indices ij the dependencies (φ r ij )  Figure 6b shows the concentration dependencies of the effect of gravitational convection, (ϕ ij ) R = f (∆C) and (ϕ det ) R = f (∆C), calculated based on Equations (12) and (13).
were negative and dependent on ∆ for the CP conditions. Compared to the conditions of homogeneity of solutions, the CP increased the value of the coefficients 11 , 22 , and and reduced the value of the coefficients 12 and 21 . For the same ∆ , the values of the coefficients 11 , 22 , and were smaller, and the coefficients 12 and 21 were higher for the nonconvective state.
Using the dependencies shown in Figures 3c,d and 5a Figures 6a and 6b shows the concentration dependencies of the CP effects ( ) and ( ) , calculated based on Equations (10) and (11). Figure 6a shows that for the same indices ij the dependencies ( ) = f(-∆ ) were asymmetric to the dependence ( ) = f(∆ ). The plots 2A, 3A, and 4A overlap for negative ∆ values, whereas plots 2B, 3B, and 4B overlap for positive ∆ values. This means that ) , and ( 11 ) were negative. According to the convention adopted in [5], the convection currents in a system containing aqueous ethanol solutions, which are a consequence of the hydrodynamic instabilities of CBLs, are directed vertically upwards. In contrast, the convection currents in a system containing aqueous glucose solutions are directed vertically downward. The results from Equations (12) and (13)  As shown in Figure 6b, for the same ∆C values of the dependency, (ϕ 12 ) R ≈ (ϕ 21 ) R ≈ (ϕ 22 ) R ≈ (ϕ det ) R < (ϕ 11 ) R . Moreover, the values of the coefficients (ϕ 12 ) R , (ϕ 21 ) R , (ϕ 22 ) R , (ϕ det ) R , and (ϕ 11 ) R were negative. According to the convention adopted in [5], the convection currents in a system containing aqueous ethanol solutions, which are a consequence of the hydrodynamic instabilities of CBLs, are directed vertically upwards. In contrast, the convection currents in a system containing aqueous glucose solutions are directed vertically downward. The results from Equations (12) and (13) state that for the nonconvective state (ϕ ij ) R = 0 and (ϕ det ) R = 0 and for the convective state (ϕ ij ) R < 0 and (ϕ det ) R < 0. Therefore, from Figure 6b we can conclude that (ϕ ij ) R = 0 and (ϕ det ) R = 0 for ∆C = 80 mol m −3 .

Concentration Dependencies ( ) and ( )
The coefficients (φ r ij ) R and (φ r det ) R were positive, and the coefficients (ϕ ij ) R and (ϕ det ) R were negative. The negative convection effect shows that the convection movements were directed vertically upwards. Moreover, for the same ∆C, the values of these coefficients were higher for the nonconvective state. . On the other hand, when comparing curves 1A and 1B to curves 2A and 2B, it can be seen that for the same values of (−∆C) and (∆C) the values of (Φ A S ) R and (Φ B S ) R for homogeneous solution conditions were greater than for CP conditions. Compared to the conditions of homogeneity of solutions, for the same ∆C, CP reduced the value of the source of entropy (Φ r S ) R . For the concentration polarization conditions, the values of (Φ r S ) R were positive and depended on both the value and the sign of ∆C. CP reduced the flux of dissipated energy (S-energy).

Concentration
Dependencies of ( ) , ( ) , ( ) , and ( ) As shown in Figure 7a, the dependencies ( ) = f(∆ ) (r = A, B) calculated based on Equation (15) were nonlinear and symmetrical with respect the point ΔC = 0 for homogenous conditions (curves 1A and 1B), whereas they were complex and asymmetrical with respect to the vertical axis passing through the point ∆ = 0 for CP conditions (curves 2A and 2B). The comparison of curves 2A and 2B shows that for the conditions of concentration polarization and the same values (−∆ ) and (∆ ) the values of ( ) (for the convective state) were greater than ( ) (for the nonconvection state). On the other hand, when comparing curves 1A and 1B to curves 2A and 2B, it can be seen that for the same values of (−∆ ) and (∆ ) the values of ( ) and ( ) for homogeneous solution conditions were greater than for CP conditions. Compared to the conditions of homogeneity of solutions, for the same ∆ , CP reduced the value of the source of entropy ( ) . For the concentration polarization conditions, the values of ( ) were positive and depended on both the value and the sign of ∆ . CP reduced the flux of dissipated energy (S-energy).   Figure 7a,b and Equations (19) and (20), the dependencies (Φ r F ) R = f (C) and (Φ r U ) R = f (∆C) were calculated. As shown in Figure 8a, the dependencies (Φ r F ) R = f (∆C) (r = A, B) calculated based on Equation (19) were nonlinear and symmetrical with respect to the point ∆C = 0 for homogenous conditions (curves 1A and 1B), whereas they were complex and asymmetrical with respect to the vertical axis passing through the point ∆C = 0 for CP conditions (curves 2A and 2B). The comparison of curves 2A and 2B shows that for the conditions of CP and the same values of (−∆C) and (∆C) the values of (Φ A F ) R (for the convective state) were greater than (Φ B F ) R (for nonconvection conditions). On the other hand, when comparing curves 1A and 1B to curves 2A and 2B, it can be seen that for the same values of (−∆C) and (∆C) the values of (Φ A F ) R and (Φ B F ) R for homogeneous solution conditions were greater than for CP conditions.
to the vertical axis passing through the point ∆ = 0 for CP conditions (curves 2A and 2B). The comparison of curves 2A and 2B shows that for the conditions of CP and the same values of (-∆ ) and (∆ ) the values of ( ) (for the convective state) were greater than ( ) (for nonconvection conditions). On the other hand, when comparing curves 1A and 1B to curves 2A and 2B, it can be seen that for the same values of (-∆ ) and (∆ ) the values of ( ) and ( ) for homogeneous solution conditions were greater than for CP conditions. This tendency was maintained for dependence ( ) = f(∆ ). Figure 8a shows that the dependencies ( ) = f(∆ ), calculated based on Equation (20), were nonlinear and symmetrical with respect the point ΔC = 0 for homogenous conditions (curves 1A and 1B), whereas they were complex and asymmetrical with respect to the vertical axis passing through the point ∆ = 0 for CP conditions (curves 2A and 2B). The comparison of curves 2A and 2B shows that for the conditions of CP and the same values of (-∆ ) and (∆ ) the values of ( ) (for the convective state) were greater than ( ) (for nonconvection conditions). On the other hand, when comparing curves 1A and 1B to curves 2A and 2B, it can be seen that for the same values of (−∆ ) and (∆ ) the values of ( ) and ( ) for homogeneous solution conditions were greater than for CP conditions.
As shown in Figures 7a and 8a The ( ) is a measure of the flux of F-energy, that is, that part of the ( ) that can be converted into useful work. The conversion efficiency of U-energy to F-energy for the same values of ∆ ranged from 12 to 36% (for conditions of homogeneity of solutions separated by the membrane). For diffusion conditions (configuration B), the efficiency of This tendency was maintained for dependence (Φ r U ) R = f (∆C). Figure 8a shows that the dependencies (Φ r U ) R = f (∆C), calculated based on Equation (20), were nonlinear and symmetrical with respect the point ∆C = 0 for homogenous conditions (curves 1A and 1B), whereas they were complex and asymmetrical with respect to the vertical axis passing through the point ∆C = 0 for CP conditions (curves 2A and 2B). The comparison of curves 2A and 2B shows that for the conditions of CP and the same values of (-∆C) and (∆C) the values of (Φ A U ) R (for the convective state) were greater than (Φ B U ) R (for nonconvection conditions). On the other hand, when comparing curves 1A and 1B to curves 2A and 2B, it can be seen that for the same values of (−∆C) and (∆C) the values of (Φ A U ) R and (Φ B U ) R for homogeneous solution conditions were greater than for CP conditions.
As shown in Figures 7a and 8a , and (Φ B U ) R < (Φ U ) R . The (Φ r F ) R is a measure of the flux of F-energy, that is, that part of the (Φ r U ) R that can be converted into useful work. The conversion efficiency of U-energy to F-energy for the same values of ∆C ranged from 12 to 36% (for conditions of homogeneity of solutions separated by the membrane). For diffusion conditions (configuration B), the efficiency of the U-energy to F-energy conversion was contained in the range from 50 to 57%, and for diffusion-convection conditions (configuration B) the efficiency ranged from 50 to 79%. Although this efficiency was relatively high, the amount of F-energy produced was small. The amount of this energy can be important in biological microsystems.
The procedure presented in this paper for evaluating membrane transport properties and energy conversion in a membrane system could be useful for any biological or artificial membranes [46][47][48][49][50].

Conclusions
This study presents the following results: • Developed within the framework of the Kedem-Katchalsky-Peusner formalism, the procedure using the Peusner coefficients R r ij (i = j ∈ {1, 2}, r = A, B) and R r det is suitable for evaluating the transport properties of polymer membranes and assessing the conversion of internal energy (U-energy) to useful energy (F-energy) and degraded energy (S-energy). • Peusner coefficients R r 12 and R r 21 are related to the membrane Peclét coefficients α r s and α r v .

•
The procedure developed in this paper to evaluate the conversion of internal energy (U-energy) to useful energy (F-energy) and degraded energy (S-energy) requires the calculation of the value of the flux of S-energy (Φ r S ) R and efficiency factors (e r 12 ) R and (e r 21 ) R , followed by the fluxes of F-energy ( Φ r F ) R and U-energy ((Φ r U ) R ).

•
The procedure proposed in the paper can be applied to membranes for which the coefficients L P , σ v , σ s , and ω can be determined experimentally. Funding: This research received no external funding.

Data Availability Statement:
The datasets for this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.