The development of membranes and membrane-based separation processes should be accompanied by a stan- dardization of the protocols applied for membrane characterization and for data analysis. Here, streamlined equations for the estimation of the water flux and of the observed salt permeability coefficient in pressure- driven processes deploying dense membranes are presented. Also, a protocol for the experimental character- ization of the transport properties of dense membranes is presented and the results are validated against the proposed equations. The proposed water flux equation is algebraic, whereas the ordinary equation needs to be solved iteratively. Moreover, in contrast to the traditional expression for the solute transport coefficient, which requires estimation of the concentration polarization, the respective equation proposed in this study only requires bulk parameters. Dimensionless variables for water flux, driving pressure, and mass transfer are introduced, and a filtration efficiency is defined, a useful parameter in terms of process design. Keywords: algebraic water flux equation; membrane characterization; dimensionless water flux; external concentration polarization; process design & optimization; reverse osmosis 1. Introduction The design and the development of the next-generation membranes for reverse osmosis (RO), nanofiltration (NF), forward osmosis, and other separation processes based on high selectivity between water and solutes or among solutes, cannot do without robust membrane characterization protocols and transport modeling tools1-4. In turn, the deployment of current and future membranes in high-value applications require the ability to predict system performance, chiefly membrane flux, in the presence of transport-limiting phenom- ena, such as concentration polarization5-7. While significant efforts are made to synthesize membranes with materials previously unimaginable, these research endeavors are often accompanied by unclear and highly differentiated characterization approaches, which limit the fair comparison between membranes, impair their adoption, and not so rarely thwart the community’s confidence in their applicability. In the ideal situation, a standardized, straightforward but robust evaluation procedure including both an experimental protocol and a data modeling strategy would be applied and would allow the evaluation of clear, univocal results by all parties, from materials scientists to the final stakeholders. In processes utilizing an applied pressure on the feed side of asymmetric, e.g., thin-film composite, membranes (i.e., RO, NF), when the objective of the research is the characterization of a membrane, the approach should provide values for the transport coefficients of the main species, namely, that related to water (also known as water permeance), A, and those related to each solute of interest, B(here referred to as observed salt permeance or as observed salt permeability coefficient), measured under relevant conditions8-10. According to the solution-diffusion transport model, these parameters are solely related to the membrane properties and to the interaction between the membrane and the relative species in the feed solution, and they do not depend ∗Corresponding author Email address: fynnja@math.aau.dk (Fynn Jerome Aschmoneit) 1 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . on the experimental conditions11. More recent models, such as the solution-friction model, have improved our understanding of membrane transport and describe instead a dependence of the observed salt permeability coefficient on feed salt concentration and applied pressure, through elimination of the assumption of a constant hydraulic pressure across the membrane and by accounting for interactions among water, salt ions, and the membrane12-14. Therefore, according to the solution-diffusion model, a uniform Bvalue would be used to compare different membranes. On the other hand, according to the more accurate solution-friction model, attention should be paid that observed salt permeability values are compared only when obtained under the same pressure and salt concentration conditions. Regardless, Bremains a powerful and simple parameter to compare different membranes deployed under similar conditions and/or applications. Unfortunately, reports on membrane characterization usually include permeate fluxes and observed rejection rates which, while interesting for practical purposes, are not as comparable and significant as permeance values15,16. Also, the related tests are rarely performed under reliably representative conditions for real applications, such as pressure values, hydrodynamics conditions, recovery rates, feed composition. On the other hand, when the transport characteristics of a membrane are known, adequate predictions of the water flux become possible under a variety of engineering conditions or applications, also with the goal to support the design of such systems4,17. One of the main obstacles for the correct predictions of water flux or for the calculation of the transport parameters, is that concentration polarization must be taken into account18. Accounting for concentration polarization means knowing or being able to predict the exact value for the osmotic pressure at the feed membrane surface, πm, whether in laboratory setups or in full-scale modules19-20. The fluid dynamics in the feed channel are governed by the convection-diffusion equation and the separation characteristics of the membrane. In particular, the solute transport in the membrane feed channel can be quantified locally or broadly, by means of the mass transfer coefficient kd. If the mass transfer coefficient is known, the concentration at the feed membrane cmcan be estimated from solving the convection-diffusion equation in the membrane boundary layer, yielding cm=cfexp(jw/kd). The usually unknown value of cmis therefore expressed through the more accessible bulk feed concentration cf, the water fluxjw, and the mass transfer coefficient. Under the assumptions of ideal thermodynamics, i.e., the osmotic pressure is linearly related to concentration, and of the intrinsic membrane Avalue being independent of salt concentration and other operating conditions, the osmotic pressures in the bulk feed and permeate, πf, πp, and the hydraulic feed pressure pfdefine the well-known ordinary water flux equation21: jw=A pf+πp−πfexpjw kd (1) This model accurately predicts the water flux under the influence of concentration polarization, when the mass transfer coefficient is known. However, it has one central disadvantage: It is a so-called transcendental equation, because the term jwcannot be isolated. As such, it has to be solved iteratively or numerically, and this necessity somewhat complicates the prediction of water flux under different conditions when salt is present in the feed solution. Such intricacy and the reliance on manual calculations often invite errors, dubious approximations, inconsistency, and significant wasted time. An important initiative has recently started to make membrane performance data and membrane evaluation results more easily findable, accessible, inter-operable, and reusable (FAIR). One such attempt is the Open Membrane Database (OMD), a web-based interface that collects data about membranes worldwide and ”allow the easy exploration and comparison of membrane performance, physicochemical properties, and synthesis conditions”22. The OMD website also includes effective explanations and calculations tools for concentration polarization and membrane performance evaluation. Another parallel project is related to the development of the so-called ”membrane-toolkit” (https://rkingsbury.github.io/membrane-toolkit/), a software serving as a library of validated calculators, with thorough documentation and high test coverage with the following goals: (i) automate routine tasks around membrane investigation to save time and reduce human error, (ii) promote standardization of membrane characterization, (iii) facilitate the creation and curation of large membrane data sets. This work fits within these ongoing efforts and its main aims are to: (a) propose and assess a robust ex- perimental protocol and a simplified equation to estimate Bfrom experimental data only based on bulk parameters, as well as (b) propose and assess a simplified non-transcendental, algebraic, equations used to 2 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . reliably estimate water fluxes in the presence of concentration polarization. The main hypothesis is that the ordinary water flux equation 1 can be expressed in algebraic form without loss of significant informa- tion following rational approximation. The validity of the simplified algebraic water flux equation is thus evaluated under an ample range of working conditions and tested against the results of the experimental characterization of several membranes following the proposed protocol. 2. Methods The algebraic water flux framework consists of three fundamental elements: (i) Dimensionless process variables for the characterization of membrane processes allow for a better com- parability of processes and a phenomenological perspective on membrane processes. (ii) The central piece of the framework is the algebraic water flux equation providing a simple method for evaluating the filtration efficiency from dimensionless bulk variables. (iii) Related characterization equations for process characterization and design optimization, using dimen- sionless bulk variables. 2.1.Dimensionless variables for filtration efficiency, pressure modulus and transportiveness Dimensionless variables allow for a better comparability of membrane processes, as will be shown later. The filtration efficiency Jis defined as the ratio of the water flux and the ideal, maximum, water flux that would be obtained without concentration polarization. Depending on the magnitude of the concentration polarization, the filtration efficiency assumes a value between 0% and 100% and therefore poses a quantity for assessing the efficiency of the filtration process. The pressure modulus, P, is defined as the ratio of the net driving pressure and the feed osmotic pressure minus the observed rejection rate, and it is positive for pressure-driven processes. Finally, the transportiveness Kis a measure for effectiveness of mixing in the feed channel. It is defined as the mass transfer coefficient, divided by the theoretical counter flow. A large transportiveness indicates good solute mixing. Only in the case of very saline solutes and at slow cross-flow rates, may the transportiveness be smaller than unity. The dimensionless flow variables are: J=jw A(pf−Rπf), P =pf πf−R, K =kd Aπf, (2) where R= 1−cp/cfis the observed rejection. Note that these three variables and the equations proposed in this study are applicable to all membrane configurations, e.g., cross-flow or dead-end systems. This is true also forK, as long as the mass transfer coefficient in the feed channel is known. Note that in the case of unsteady flow, e.g. in dead-end systems, the mass transfer coefficient is time-dependent. Also, note that the three dimensionless variables are calculated using macroscopic experimentally-observed parameters, with the only exception of A. One of the main assumptions of this study is that Ais independent of operating conditions, and its definition strictly follows the solution-diffusion model. However, as explained by previous studies discussing the solution-friction model, the friction between water and the membrane, which is independent of salt concentration, dominates the hydraulic pressure drop across the membrane, thus resulting in near stable water permeability under different operating conditions12-14. Therefore, the three variables and the equations proposed in this study are applicable under practical conditions from the point of view of different models used to describe membrane transport. 2.2.Equations for process design and membrane characterization: water flux The dimensionless process variables above are used to derive simpler expressions for the water flux and related quantities. The equations presented in this section are based on mathematical derivations, which can be found in the Supplementary Information. Equation 3 represents the dimensionless version of the ordinary water flux equation 1. This equation is in itself not simpler, nor more expressive than the full water flux equation, since it remains a transcendental 3 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . function of the filtration efficiency, J. However, it states that the efficiency is one minus the very right- hand term, which can therefore be directly linked to the effect of concentration polarization. In the limit of perfect solute mixing, K→ ∞ , the concentration polarization term becomes zero and the filtration efficiency equation reduces to J= 1. That relates to jw=A(pf+πp−πf), i.e., a perfect filtration efficiency with no concentration polarization11. J= 1−1 P expJP K −1 | {z } concentration polarization(3) The fundamental problem with this equation is its transcendental character, where the filtration efficiency Jis both on the left-hand side and inside the exponential function on the right-hand side. It is shown in the Supplementary Information (Note 1) how the filtration efficiency Jcan be written as a function of the pressure modulus Pand the transportiveness K, where the Jstands isolated on the left-hand side: J= 1−1 1 +K−PK 2(1 + K)3 | {z } concentration polarization(4) In contrast with the ordinary water flux equation, this equation is algebraic, and can be conveniently solved without the need for iterative solvers. It is therefore referred to as the algebraic water flux equation, consti- tuting one of the central items of this article. It must be emphasized that the algebraic water flux equation is an approximation to the ordinary water flux equation. It is therefore also based on the classical convection- diffusion model. Similarly to the ordinary equation, it is seen that the concentration polarization term disappears, under the assumption of perfect mixing, K→ ∞ . On the other hand, when the cross-flow is stagnating and K→0, the water flux will converge to zero as a consequence of overwhelming external concentration polarization (ECP). In this form, the algebraic water flux equation expresses the filtration efficiency as ’one minus the effect of concentration polarization’, that makes it especially convenient for the analysis of membrane processes. For example, suppose that a lab experiment or that a system is designed with ( P, K) = (4 ,6), yielding a filtration efficiency of J≈82%. That directly implies that 18% efficiency is lost to concentration polarization. These values for P and K could exemplarily relate to a brackish water process with A=4 LMH/bar, R=98%, pf=12 bar πf=4 bar, kd=96 LMH, or a wastewater process with A=10 LMH/bar, R=95%, pf=2.3 bar, πf=0.75 bar, kd=45 LMH. For completeness, the water flux is also presented in absolute terms in equation 5. It follows from substituting the dimensionless variables in equation 4. This equation is the algebraic approximation to equation 1. jw=A(pf−Rπf) 1−Aπf Aπf+kd−A2(pf−πfR)πfkd 2(Aπf+kd)3 | {z } =J(5) Comparing equations 4 and 5 demonstrates the benefit of using dimensionless variables as they not only exemplify physical meaning, but also improve readability and conciseness. Equation 5 shall remain here as a stand-alone result, while the remaining of this article is concerned with the dimensionless water flux equation 4. As the algebraic water flux equation 4 is an approximation, it must be analyzed how much it deviates from the ordinary water flux equation. Error quantification:. Figure 1a presents contour maps of the filtration efficiency and Figure 1b shows the deviations introduced through the approximations in the algebraic water flux equation. Both contour plots deploy the same, representative, ranges of pressure modulus Pand transportiveness K. The dashed and solid lines in Figure 1a show the filtration efficiency for the ordinary water flux equation 3 and the algebraic water flux equation 4, respectively. Qualitatively, it is seen that both equations predict that the filtration efficiency 4 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . increases significantly for increasing Kvalues and moderately for increasing Pvalues. Furthermore, it is seen seen that the approximated, algebraic water flux equation deviates from the ordinary equation for increasing Pvalues. For filtration efficiencies J >60%, the algebraic equation underestimates the water flux, while it overestimates the water flux where J <60%. Figure 1b illustrates a quantification of how severe the model deviations are. The continuous contour lines indicate the relative error of calculating Jwith equation 4 (algebraic) vs. calculating Jwith equation 3 (ordinary). The derivation of the algebraic water flux equation makes two central assumptions; see Supplementary Information. The contour lines of relative error illustrate nicely how these underlying assumptions create two regions of inaccuracies (the region ( P > 5, K > 3) relates to assumption 1 and the region ( K < 3) relates to assumption 2; see Supplementary Information). While inaccuracies in the top-right region are negligibly small, the region for diminishing Kvalues indicates close error contours and hence, significant inaccuracies of the algebraic water flux equation. The region of great inaccuracies can conveniently be related to the empirically found relation 4 P > K (1+K)2, which corresponds to the red area in Figure 1b. The black dashed line indicates the contour of J= 50% efficiency. It is seen that for P > 3 the region of increasing inaccuracies corresponds to J <50%. This allows to define a practical condition for the validity of the algebraic water flux equation as 4 P < K (1 +K)2, or, as a rule of thumb, J >50%, where P > 3. Error discussion:. In practical terms, the validity condition basically renders the algebraic water flux equation invalid for very small mass transport coefficients in the feed chamber/channel only, where the concentration polarization is extreme. It is shown in the following that this limitation is of little practical significance. Table 1 presents process exemplifications of three typical medium and high-pressure membrane applications, namely, seawater, brackish water, and wastewater desalination. For each of these applications, the table presents typical operating conditions and the flow variables at the feed inlet and outlet of a membrane. The values of the mass transfer coefficients were chosen conservatively (i.e., lower than typical values), causing significant concentration polarization in each of the processes. From the operation conditions follow the (P, K) operation points at the inlet and outlet, using the definitions in equation 2. The third-last row in Table 1 indicates that all ( P, K) operation points pass the validity condition, see also Figure 1b for locations of these operation points in the map of relative errors. The filtration efficiency Jwill therefore be very well- approximated with the algebraic water flux equation 4. It is seen in the second-last row that the filtration efficiency is much greater at the inlets, while it reduces to only 30 −35% at the outlet, which is due to the build-up of concentration polarization caused, in turn, by the high concentration of the retentate stream (and a small mass transfer coefficient in cross-flow systems due to low flow rates in the feed channel at the outlets). In conclusion, only if concentration polarization is impractically high, the algebraic equation is not accurate. However, such conditions are virtually never found in real or in laboratory applications23. The algebraic water flux equation therefore accurately reproduces the prediction of water flux of the ordinary equation for a wide range of realistic KandP. Its accuracy is impaired only for conditions that are rarely found in laboratory or full-scale applications. 2.3.Equations for process design and membrane characterization: observed salt permeability coefficient Focusing the attention to the simplification of an expression to calculate the observed solute permeance of the membrane, B. Note that in this study, Bis regarded as an experimental based parameter and does not necessarily refer to an intrinsic membrane property. Bis defined as the value of salt permeance observed under specific and fixed operating conditions. It should be reminded that this parameter requires accounting for concentration polarization, since it depends on cm, the feed concentration at the membrane interface, which is affected by the concentration polarization. This is commonly accomplished through indirect experimental- based estimation of the solute concentration at the feed/membrane interface and/or using the mass transfer coefficient, kd. The mass transfer coefficient kdis a complex quantity in the sense that its value depends on the operation conditions through, e.g., the cross-flow velocity, the feed concentration and the solute’s diffusion coefficient, but also on the membrane geometry through, e.g., channel widths, spacers, fouling24. However, a reasonable estimate for the mass transfer can be obtained as a function of the other process variables. Therefore, within the framework of the dimensionless process variables, the transportiveness K can be expressed as a function of the filtration efficiency Jand the pressure modulus Pas: 5 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . K=JP ln( 1 + P(1−J)|{z} CPmod)(6) This expression for the transportiveness follows directly from the dimensionless water flux equation 3, see derivation in the Supplementary Information. It is therefore not an approximation and no errors are intro- duced. The argument in the logarithm function of equation 6 is identified as the concentration polarization modulus πm/πf, which yields this very simple expression for the concentration polarization modulus: CPmod= 1 + P(1−J) (7) Note that, quite significantly, in this way CPmodandKare determined without the need to know cmorπm, i.e., without the need to know or estimate the mass transfer coefficient kd. On the contrary, equation 7 can be used to estimate πmandBin a simple manner and only considering bulk values of the parameters, by applying the following equation: B=jw1−R CPmod−1 +R(8) The derivation of the above equation can be found in the Supplementary Information. It does not include any approximations and will therefore not introduce inaccuracies. Importantly, equation 8 deploys Jand P, which are based on readily available bulk parameters, only. Additionally, as already mentioned above, this equation allows calculation of an observed salt permeability coefficient, not necessarily an intrinsic membrane parameter. No assumptions are made in this study regarding the behavior of Bas a function of operating conditions. When conditions change, e.g., the salt concentration is higher in the feed solution, the solution-friction model predicts a different, e.g., larger, value for B. In fact, this behavior would produce condition-specific experimentally accessible values of jw,CPmod, and R, which are applied in equation 8 and translate into a respective condition-specific value of the observed salt permeability coefficient, B. In this sense, this simple equation based on experimentally available bulk parameters is applicable from both the point of view of the solution-diffusion and the solution-friction model. It should be noted that, according to the solution-friction model, the observed salt permeance, here referred to as B, approaches the value of an intrinsic membrane salt permeance when membranes are not charged or when the feed salt concentration is high compared to the membrane charge density. 2.4.A robust protocol for the experimental characterization of membrane transport In this section, a protocol aimed at characterizing the transport properties of a membrane comprising a dense active layer is proposed. Such a protocol should be easily reproducible and practical so as to be used as a standard practice, should be as simple as possible while being as detailed as necessary, and it should be reliable and consistent in terms of outcomes and for that reason involve overdetermined data. These characteristics are collectively referred to as ’robust’. Setup and operating conditions: the membrane characterization rig should include systems to control tem- perature, pressure, and, when relevant, cross-flow. When the membrane cross-flow cell is small, each of this parameter may be measured at the feed inlet or concentrate outlet, only. However, for large cross-flow mem- brane housings, both inlet and outlet values should be obtained. Experimentalists should use deionized water as a feed solution to determine Aand should prepare a stock solution of concentrated salt or salt mixture to be dissolved and diluted into the deionized water solution to determine rejection rates and B. Experimentalists should be able to maintain steady-state when needed, for example by running a cross-flow rig in closed-loop, i.e., concentrate and permeate streams recirculated into the feed tank; however, any configuration that allows keeping conditions (e.g., composition of the feed tank) constant in time are suitable. Probes or analytical instrumentation for pH and solute measurements are necessary to ensure the desired water composition and to measure rejection rates. The conditions deployed for membrane characterization, including feed pressures 6 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . as well as the nature and the concentration of the solute(s) in the feed solution, should be chosen in the range that is relevant or representative of the specific application for which a membrane has been fabricated. An important parameter related to water composition is pH, which should be checked and adjusted to the desired value. With regard to the sampling of the feed and permeate solutions, attention shall be paid that the total volume of all samples collected during the test is negligible compared to the initial feed solution volume. Execution of the test: a first compaction step should be performed, using deionized water as feed solution (phase 0). This step should be run at a higher feed pressure compared to the pressure values subsequently used for characterization, and until steady-state in flux is achieved. At this point, water flux should be measured at varying pressure, each time at steady-state (phase 1). The authors suggest taking measurements with at least three different feed pressures. Subsequently (phase 2), the concentrated stock solution containing solute(s) should be added into the feed tank and pH adjusted to obtain the desired feed composition. This step should be done while letting the system run. For the characterization phase of membrane selectivity, it is advised to obtained different values of water flux and solute concentrations at various combinations of feed applied pressure and mass transfer coefficient in the feed channel/chamber (obtained, e.g., varying the cross-flow velocity in cross-flow configurations). The authors suggest using at least three combinations. Based on the current understanding of membrane transport, water flux should always increase with feed pressure at a constant value of the mass transfer coefficient, due to an increase in bulk driving force in the feed channel/chamber, and it should also increase at increasing value of the mass transfer coefficient if the feed pressure remains constant, due to lesser concentration polarization. For analogous reasons and according to currently available transport models, observed rejection should increase if either feed pressure or mass transfer coefficient are increased. These trends may be used to verify whether the experiment is running in accordance to expectations. It is imperative that all values of water flux and solute concentrations used for subsequent analyses are recorded when the system is at steady-state. Therefore, a sufficient amount of time should be allowed upon each change of conditions, also to make sure that all tubings/pipings are well flushed with the solutions relative to the new conditions obtained after changing the parameters. Analysis: The water flux values obtained with deionized water in phase 1 should be fit with a line and this line should pass through the point of zero flux for zero applied pressure (origin of the flux vs. feed pressure graph). If the intercept at zero pressure is significantly distant from the zero value, determination of Amay not be accurate and this result may be due to insufficient compaction of the membrane, membrane defects, or experimental deviations occurred during the test. The Avalue is the slope of the line that best fits the water flux data and that passes through the zero-zero point. For the determination of membrane selectivity from the data collected in phase 2, rejection rates of each solute under each condition are calculated using the concentrations of that solute determined in the feed sample and in the permeate sample, both related to the same time of sampling. Note that when a proxy parameter is used in place of solute concentration, for example, electric conductivity in place of salinity, experimentalists should consider that the proxy parameter may not correlate linearly with the actual parameter in the entire spectrum of values relevant for the test. Therefore, direct substitution of one with the other in the equation used to calculate rejection rate may provide inaccurate results. A calibration curve providing the exact correlation between the proxy and the actual parameter should be determined in advance and then applied to translate the proxy parameter into the actual concentration. A value of Bcan be thus obtained for each solute and for each operating condition, by applying equation 8 proposed in this study. It is important to highlight here, once again, that according to the solution-diffusion model, all Bvalues obtained for the same solute in phase 2 are supposed to be equal, regardless of the conditions: thus, according to this model, an average value may be presented. However, more recent and accurate models, e.g., the solution-friction model, suggest that observed Bvalues are not intrinsic to the membrane but they depend also on the operating conditions: thus, presenting an average value is not relevant or correct, and separate observed Bvalues obtained with equation 8 shall be presented, together with the conditions under which they were measured. 2.5.Proposed protocol in action: Experimental conditions and analyses applied in this study The transport properties of various polyamide membranes characterized by active layers of different densities were evaluated using a laboratory-scale cross-flow unit25. The unit comprises a high-pressure pump, a feed tank, a flat membrane housing cell, and a chiller with heat exchanger coils immersed in the feed tank for 7 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . temperature control. The effective membrane active area was 20.1 cm2and the temperature was constant at 23±0.5°C. Membranes suitable for processes classifiable as seawater reverse osmosis (SW), brackish water reverse osmosis (BW), and nanofiltration (NF) were deployed. Prior to each experiment, the membrane sample was immersed in water overnight. The filtration tests consisted of two different phases: initially, deionized water (resistivity >107Ohms) was used as feed solution to evaluate the water permeance of the membrane, A; subsequently, an appropriate volume of NaCl stock solution (stock solution concentration = 5 mol/L) or of MgSO 4stock solution (1 mol/L) was directly added into the feed tank to evaluate rejection rates and the solute transport coefficienc, B. The pH was fixed at 8.0 by addition of a minimal amount of buffer compound (NaHCO 3) and via adjustment with NaOH. The solute concentrations were consistent with those typically utilized by membrane manufacturer for standard membrane testing and commonly reported in the specification sheets. In the first phase, the applied feed pressure, pf, was changed to obtain different values of the water flux as a function of pfwith a feed solution of deionized water. In the second phase, both pfand the cross-flow velocity (cfv) were changed to obtain different measurements of the permeate flux and of the solute rejection in the presence of solutes in the feed solution. Specifically, three cfv values were investigated referred to as high, medium, and low cfv. These values were chosen because consistent with those typically encountered in spiral- wound elements installed in reverse osmosis and nanofiltration plants and consistent with diminishing flow rates along the elements within the plant. The values of pfand solute concentration were chosen according to the density of the membranes, higher for the membranes with denser active layers and lower for the ones with lower expected rejection. All the testing conditions can be found in Table 2. According to currently available models, the observed rejection, R, is a function of applied feed pressure and feed salt concentration, and is also affected by concentration polarization, which, in turn, is influenced by hydrodynamics conditions, e.g., cross-flow velocity26. In the beginning of each test, the membrane sample was compacted with deionized water as feed solution at the highest value of applied pressure until the permeate flux reached a steady-state (generally 2-3 h)27. In this first phase involving deionized water as feed solution, pfwas then lowered in a step-wise fashion. In each step, the pressure was changed gradually to avoid shocks in the system and to the membrane, and the water volume passing through the membrane was then measured by means of a computer-interface balance; see also Figure 2a for an example of experimental data related to membrane ”SW-1”. The pure water flux was calculated by dividing the volumetric permeate rate, obtained at steady-state (reached typically after only a few minutes), by the membrane active area. Awas determined as the slope of the best fitting line for the water flux data as a function of pf, with the line passing through the origin (Figure 2b). In the second phase, after addition of solute in the feed solution, five different steps were performed under operating conditions consisting of different combinations of the same values of pfinvestigated in the first phase and of three cfv values (Table 2). Solute concentrations in the feed and permeate streams were calculated from conductivity values measured using a conductivity meter (Oakton CON 450), calibrated for each salt. The permeate flux, jw, was calculated by dividing the volumetric permeate rate by the membrane area. R, was then computed from the concentrations determined in bulk feed, cf, and in the permeate stream, cp, as R= 1−cp cf(9) The observed rejection rates and the permeate fluxes were always measured at steady state (reached typically after 10-15 minutes after each change of condition but probed after roughly 20-25 min from each change of conditions). Therefore, within each step, i.e., for each combination of pfand cfv, the values of these parameters were always constant in time, within experimental error. Two separate measurements were performed, distanced 10-20 min from each other and the values were averaged. Except for collection periods, both the concentrate and the permeate streams were recirculated back into the feed tank. In the end, B was computed from experimentally available bulk data: first, JandPwere calculated from the input or measured data of the experiment. Then, equation 7 was applied to calculate the concentration polarization modulus. Finally, Bwas obtained using equation 8. 8 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . 3. Results and Discussion 3.1.Results of experimental membrane characterizations Figure 2 presents the experimental results obtained with SW-1, as representative membrane. Figure 2a shows the water flux data measured as a function of time in the various phases and steps of the experiment. The average water flux data obtained at steady state in the initial phase of the test is thus reported in Figure 2b as a function of pf, where the best fitting line passing through the origin is shown for the calculation ofA. Figure 2c presents the values of CPmodandR(NaCl) evaluated in the five steps of the second testing phase as a function of pf, in the presence of 32 g/L NaCl in the feed solution (pH 8.0). As expected from theoretical considerations, the permeate flux increased with increasing feed pressure. Consequently, the observed NaCl rejection and the CP modulus also increased. More interestingly, flux and rejection data increased slightly but significantly with increasing cfv at a given value of pf, thus the value of CPmod decreased. Higher cfv increased the mixing in the feed channel, reducing the thickness of the unmixed boundary layer and reducing the magnitude of ECP28. This phenomenon translated into a lower solute concentration at the feed-membrane interface, c f,m, which in turns allows a higher effective driving force and lower salt passage across the membrane active layer. These observations were consistently achieved for all membranes, suggesting the reliability of the experimental protocol and the accordance between experimental results and conceptual understanding of the phenomena underlying mass transport across the active layer of asymmetric membranes21. Figure 3 summarizes the results in terms of Aand average Bfor all membrane types. The data are plotted for the six membranes, from the least permeable to the most permeable one from left to right. Just for simplicity and conciseness, the observed salt permeances are here presented as average values for each of the membrane. The individual values of B, obtain under each operating condition, can be found in the Supplementary Information. As expected, the highest productivity is achievable with NF membranes, followed by BW and SW membranes, respectively. The values of average Bcorrelate well with those of the parameter A, except for SW-2, which displayed both better productivity and rejection rate than SW-1. Note that the value of average Bestimated for NaCl with the NF membrane is significantly higher than that estimated for MgSO 4, since the latter solute includes a divalent cation, hence associated with better rejection29. More importantly, note that the standard deviations for the parameter Bare relatively small, i.e., low coefficient of variation, despite the fact that these are average of the five different steps conducted at varying feed pressure and cfv combinations. 3.2.Analysis of the filtration efficiency using the algebraic water flux equation Figure 4 shows the experimental data in the framework of the dimensionless variables, while the unprocessed experimental data are found in the Supplementary Information (see Tables). The curves are contours of the pressure modulus in a J−Kmap, calculated with the algebraic water flux equation. Apart from the SW-3 data, all data reside in regions where the accuracy of the algebraic water flux equation is at least 99% (97% for the BW-1 membrane). Hence, the figure indicates how well the experimental data adhere to the convection-diffusion model of polarization, as well as how robust the data are in terms of experimental estimation of the hydrodynamics parameters. A first take-home message from the graphs is that the experimental data are much more in line with the currently available transport models as the density of the membrane active layer increases. Note that the scale of the y-axis is different for the various graphs, with SW-1, SW-2, SW-3 utilizing a smaller range ofJ. This result is consistent with theoretical expectations, since mechanisms of partition of the solvent and of the solutes in the membrane and their diffusion across the active layer become relatively less im- portant compared to other mechanisms of transport, e.g., Donnan exclusion, as the ratio between species and membrane pores is reduced30. Even more importantly, in this work the model was computed assuming that the reflection coefficient is equal to 1, i.e., impermeable solute, which is only a fair approximation for high-rejection membranes31. Note that the experimental data almost always sits above the theoretical curves for all membranes, i.e., higher Jvalues, and that for the NF membrane the consistency of the data with the theoretical curves improves for MgSO 4compared to NaCl. Both these observations indicate that higher rejection rates undoubtedly allow to safely neglecting the reflection coefficient. When considering the effect ofK, namely, of hydrodynamics, note that the width of the horizontal error bars imply a certain uncertainty 9 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . in accurately estimating the value of the mass transport coefficient, one of the main obstacles of membrane characterization, also highlighted above. Additional noteworthy conclusions can be draw by assessing the absolute values of J, which may be thought as a filtration efficiency or, in other words, how much of the nominal driving force actually goes into producing a water flux. It is important here to underline the difference between filtration efficiency, which refers to the inner workings of a membrane-based system, and absolute productivity, which refers to the amount of product water obtained per unit time from the same system. While the two variables are connected, the latter is what water utilities are most concerned with and it represents a design objective of the plants, while the former is associated to the means to achieving these objectives. Looking at the behavior of three seawater membranes (SW-1, SW-2, SW-3), all tested at the same value applied feed pressures, the filtration efficiency dropped dramatically from the least permeable to the most permeable membrane, despite the fact that the water fluxes were obviously higher with the latter. This observation implies that attempting to increase flux above a certain range by applying a high applied feed pressure, produces only marginal returns, as the increase in flux brings about a sustained concentration polarization that in turn limits the flux increase itself. Therefore, the energy expense associated with higher feed pressures is not entirely justified, implying that the driving force should be adjusted for each membrane water permeance to be within a certain range, if the goal is to improve efficiency. In real applications, the system productivity is often set by the needs of an industry or a community and the degree of freedom in adjusting filtration efficiency may be lower. However, the results of this study suggest that an increase in membrane area may be more advantageous than that of applied feed pressure, to maintain overall productivity while also increasing efficiency. Indeed, economic considerations are outside the scope of this study, and they must be taken into consideration for real-scale operation. 3.3.Application of the algebraic equation in process design and membrane characterization Figure 5 presents sensitivity maps from the implementation of equations presented above in terms of various interdependencies among K,P,J, and CPmod. Note that such maps do not attempt to aid in overall system-scale design, which should be performed using appropriate modeling tools and should account for the flux and concentration profiles along the entire plant (also considering different stages, passes, recirculations, and bypasses). However, the maps may be used to estimate individual water flux values related to specific combinations of operating conditions (hence, if desired, every water flux value in its changing profile along a membrane plant). Also, the maps nicely illustrate the underlying principles and implications of the algebraic equation and the usefulness, and indeed limitations, of the dimensionless parameters introduced in this study. Specifically, Figure 5a is a map of equation 4, Figure 5d is a map of equation 7, while Figures 5b,c are alternative representations of Figures 5d,a, respectively. Several conclusions may be drawn about the strong or weak dependency of the variables on each other. For example, at a fixed value of P, the filtration efficiency can only be increased by increasing the mass transfer coefficient in the feed channel, hence higher K. On the other hand, at a fixed value of K, the filtration efficiency can partly be also increased by reducing P, that is, by working with a smaller driving force and a smaller overall productivity. To exemplify this discussion, a hypothetical optimization strategy may be assumed, in which an initial process with 80% filtration efficiency should be modified to reach a value of efficiency equal to 85%. If the process is characterized by P= 6, K= 5.9 (see starting point for the two arrows in Figure 5a), the ECP can be reduced by moving into different directions. The red vertical arrow relates to the case where the feed pressure is held constant and Kis increased, i.e., the cross-flow velocity. The blue horizontal arrow indicates the case whereby ECP is reduced by lowering the feed pressure at constant cross-flow. Therefore, Figure 5a suggests that the effect of Kis more significant in influencing Jthan that of P. A strategy for optimizing a process in terms of filtration efficiency would thus favor adjusting the cross-flow velocity rather than the pressure. Similar to Figure 1a, the orange region indicates ranges in PandK, for which the water flux equation is invalid. Figures 5b-d indicate how Jand the CP modulus change in the same optimization process. The initial process with 80% efficiency has a CPmodof 2.2. When reducing the feed pressure at constant cross-flow (blue leftward arrow), the CPmoddrops to 1.1 as a filtration efficiency of 85% is reached. In the case of increasing cross- flow at constant feed pressure, the CPmodis instead reduced from 2.2 to 1.9. This observation implies that reducing the CPmodcan be very efficiently done by lowering the feed pressure, while increasing the cross-flow only has a limited effect. Note that the two outcomes in Figure 5a and Figure 5d are not in contradiction, 10 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w . but they actually suggest something less than trivial and related to the definitions of filtration efficiency and CPmod.CPmodindicates how much concentration exist at the membrane-feed interface, but it does not necessarily indicate how much water flux is lost with respect to ideality. On the other hand, Jdoes not indicate what the concentration is at the membrane-feed interface, but rather how much the water flux will be reduced because of it. Figure 5 may therefore be used to help design a filtration process. If the goal is maintaining a high filtration efficiency, thus allowing the correct exploitation of a certain driving force and membrane transport properties, the combinations of PandKcan be determined from the maps for a certain target value of J. Based on the membrane properties and on the needed system productivity, one can then calculate the required values of the absolute design variables, pfandkd. Or alternatively, if the goal is to help choosing an appropriate membrane, the required value of A, that is, the most appropriate membrane for a certain application, can be estimated to achieve a certain fixed productivity or filtration efficiency, known or hypothesized the operating conditions of a system. In summary, the availability of a robust standard protocol for membrane characterization and a simple way to estimate the transport parameters from experimental data would incentivize the adoption of common practices in the membrane field, with positive implications for membrane development and for the progress of science through clear and shared gathering, curation, and interpretation of data. Furthermore, an equation that allows for the straightforward estimation of the water flux across the active layer of dense membranes, without the need of numerical methods, would allow for the streamlined exploration of the productivity of a system under a wide range of operating conditions. It would also promote understanding of the functioning of different membranes characterized by diverse transport parameters and comparison between membranes and materials. The equations proposed in this study include dimensionless parameters with physical meaning, all based on bulk values, and they are conceived so that their terms are strongly correlated to the efficiency of the process. Indeed, the highlight of the equation terms on system efficiency and the possibility to easily estimate the magnitude of concentration polarization allow for a better understanding of the performance of a system and of a membrane, beyond sole assessment of productivity. Acknowledgements M.Mo. would like to thank Eni S.p.A. for funding her Ph.D. scholarship. Competing Interests The Authors declare no Competing Financial or Non-Financial Interests. Data availability The data generated during the current study are found in the Supplementary Information. Author contribution A.T.: conceptualizing protocol, data analysis, writing of manuscript M.Ma.: generating experimental data, writing of manuscript M.Mo.: generating experimental data, writing of manuscript M.G.: data analysis, writing of manuscript F.J.A.: conceptualizing & derivation of equations, writing of manuscript 11 Accepted version. This article appeared in npj Clean Water 6, 58 (2023) and may be found at https://doi.org/10.1038/s41545-023-00270-w .