A Molecular Dynamics Study on Rotational Nanofluid and Its Application to Desalination

In this work, we systematically study a rotational nanofluidic device for reverse osmosis (RO) desalination by using large scale molecular dynamics modeling and simulation. Moreover, we have compared Molecular Dynamics simulation with fluid mechanics modeling. We have found that the pressure generated by the centrifugal motion of nanofluids can counterbalance the osmosis pressure developed from the concentration gradient, and hence provide a driving force to filtrate fresh water from salt water. Molecular Dynamics modeling of two different types of designs are performed and compared. Results indicate that this novel nanofluidic device is not only able to alleviate the fouling problem significantly, but it is also capable of maintaining high membrane permeability and energy efficiency. The angular velocity of the nanofluids within the device is investigated, and the critical angular velocity needed for the fluids to overcome the osmotic pressure is derived. Meanwhile, a maximal angular velocity value is also identified to avoid Taylor-Couette instability. The MD simulation results agree well with continuum modeling results obtained from fluid hydrodynamics theory, which provides a theoretical foundation for scaling up the proposed rotational osmosis device. Successful fabrication of such rotational RO membrane centrifuge may potentially revolutionize the membrane desalination technology by providing a fundamental solution to the water resource problem.


Introduction
Clean water resources are essential to the survival of human life and human civilization. Developments of new techniques for water treatment and desalination have been a focal point for both sustainability research and environmental engineering [1]. In the past several decades, the membrane-based reverse osmosis (RO) technologies have been more prevalent and have accounted for more than 40 percent of the total desalination capacity [2], due to its energy efficiency compared to other desalination methods [3], such as thermal-based multi-stage flash method [4] and electrodialysis method [5]. Among all membrane materials, nanomaterials have many advantageous physicochemical properties that make them particularly attractive as separation media for water purification [6][7][8][9][10][11].
Recently, exciting progress has been made in the developments of carbon-based membrane materials and structures (such as Graphene, Graphene Oxide (GO) and Carbon Nanotube (CNT) [12][13][14][15][16]), which are identified to be excellent candidates for RO desalination because of their porous microstructures and excellent ion rejection abilities. Graphene membranes have excellent impermeability to molecules as small as helium [12,13] and a high water flux rate because of the property of hydrophobicity [14]. Surwade et al. manufactured a nanoporous graphene (NPG) membrane and observed that water molecules could pass through the membrane quickly with a very high salt rejection rate by controlling the radius of the nanopores on the NPG membrane [17]. Grossman theoretically investigated graphene membrane desalination and found NPG can not only enhance the rejection of ions, but also improve the permeation of water [18,19]. Graphene Oxide (GO) membranes are believed to be suitable alternative materials to overcome problems existing in pristine Graphene membranes. The existence of functional groups (such as hydroxyl, epoxy and carbonyl groups) [20,21] makes the GO membrane more hydrophilic and more accessible for synthesis. Hummer et al. studied water molecules flowing through CNT membranes and observed a high flow rate [15], which was verified by MD calculation showing that water flux rate inside the CNTs can be several magnitudes higher than that in the open environment [16].
However, several inherent drawbacks from the materials mentioned above hinder their practical applications in real desalination systems. First, self-contamination (fouling) of the membrane material can significantly degrade the efficiency of desalination. The fouling problem commonly exists in almost all current RO membrane methods, due to the high salt concentrations around nanopores and complicated physicochemical interactions between membrane materials and solutions [22,23]. One of the main reasons for fouling is initiated from the applied high hydrostatic pressure in the conventional desalination procedures [24]. Although this high hydraulic pressure on the porous membrane is necessary to balance osmotic pressure and maintain high freshwater flux [18], it also facilitates the blockage of pores on the membrane by pushing more ions inside it. Second, the high hydrostatic pressure on top of the membrane can decrease the interlayer distance between a layered membrane material. The narrower channels inside the membrane make it harder for water molecules to flow through the membrane and therefore decrease the filtration efficiency. Meanwhile, the high pressure can potentially fracture the membrane, which means the selected membrane material has to satisfy mechanical stability. Finally, although CNT membranes are natural porous materials, it is problematic to select suitable diameters of CNTs so that only water molecules can pass through while other ions are blocked. It is also challenging to align and assemble such CNTs into membranes and use them in real experiments.
In order to minimize the fouling problem, we propose a nanofluidic device-a rotating system with a nanoporous membrane, motivated by a centrifuge. A scaling-up strategy to fabricate such a nanofluidic device at the macroscale has previously been proposed in [25]. By utilizing the advantages of Couette flow around the nanoporous wall of the device, we theoretically proved that the issues mentioned above can be resolved and a higher desalination efficiency can be achieved. As shown in Figure 1 and in Appendix B, holes with mesoscale sizes are created on the wall of the macroscale size container, which separates the fluid into two reservoirs. The first holding saltwater inside the container and the second holding the pure water outside the container. The mesoscale holes are on the order of micron, which are filled with holes on the micro-scale that are in turn filled with nanopores.
Nanoporous membrane materials (such as Graphene or GO) are filled in the mesoscale holes with sub-nanometer pores on the porous material. An osmotic pressure P os pointing inward of the container is developed at the interface of fresh water and salt water, because of the concentration difference of ions between them. Naturally, water molecules will flow into the container because of this osmosis pressure; however, when a rotator inside the centrifuge rotates with an angular velocity ω, a centrifugal force will be developed among the ions and the water molecules. The centrifugal force grows larger with larger ω. Reverse osmosis will happen if the angular velocity ω is large enough so that the developed "centrifugal pressure" P ext (centrifugal force per unit area) is greater than the osmotic pressure: ∆P = P ext − P os > 0. The net outward pressure ∆P will drive the water solution inside the container to flow outward, but the finite sized pores will block the ions from exiting.
In this work, we use the classic molecular dynamics (MD) method to model and simulate such a rotating nanoporous structure at the submicron scale. Its scaled-up macroscale model is also analyzed by applying conventional fluid dynamics theory in a comparison study. Both the nanoscale and the macroscale structure are studied and compared. This paper is structured as followed: In the Results section, velocity profiles, slipping properties, and Taylor instability are investigated with data from the MD simulations. The relationship between angular velocity, equivalent pressure, and volumetric flow is discussed. The spiral motion of salt ions is provided as evidence for preventing fouling with the current model. Theoretical analysis based on Couette flow theory is employed to link the ion transportation through the nanoporous membrane with pressure profiles and volumetric flow rates. In the Discussion section, key findings from the results are summarized and the possibility of scaling up the model to a macroscale system is discussed, which can be used in real engineering desalination or water purification practice.

Model Description
The molecular models are shown in Figure 2. Two different models are built to study the effect of two different boundary conditions on the rotational fluid fields. In both models, the outer-layer large tubes are a nanoporous membrane (such as Graphene or GO) with different sizes and numbers of pores. The rotator that generates a rotational fluid field and develops the centrifugal force is modeled as a small tube in Model I, and a blade in Model II. Notably, for simplicity, hydrophobic materials (Carbon Nanotube for Model I and Graphene for Model II) are used here in construction and modeling of rotator. In the real application, hydrophilic materials are preferable for the rotator in order to increase the friction between the fluid and the rotator, while hydrophobic materials are suitable for the porous membrane. However, for comparison purposes, we build a rotating blade in Model II in order to generate as large a rotating fluid as possible, which functions similar to the behavior of a rotator made of hydrophilic material.
In all MD calculations, outer-layer tubes have the same radius R 0 = 4.2 nm and same length L 0 = 9 nm. On the wall of the outer-layer tubes, twelve pores with radius r p = 0.42 nm are created. (Only six pores are shown in Figure 2, and other pores may be seen when rotating the centrifuge).
The size of pores is selected according to the suggestion in the literature e.g., [26,27], so that water molecules can flow through the membrane freely while ions are blocked. The inner-layer tube in Model I and the blade in Model II have the same radius R i = 2 nm with varying angular velocity ω, to generate different fluid motions. However, the mechanisms of the generated rotating field are different. In Model I, the water swirling motion is generated by the viscous force of the rotation of inner-layer tube, whereas in Model II the rotation is generated by a direct propelling force exerted by the inner blade. Fluid velocity in Model I will be smaller than that in Model II under the same angular velocity applied on the two types of rotators. Models are solvated into water boxes with size 10 nm × 10 nm × 9 nm. The salinity of seawater is around 35 grams/liter (corresponds to a NaCl concentration of ∆C = 0.6 M), which can be equivalently converted into 66 pairs of Na + and Cl − ions with a total number of 9770 water molecules in Model I, and 96 pairs of Na + and Cl − ions with a total number of 14,110 water molecules in Model II.

Motion of Water Molecules and Ions
Snapshots from a typical run of Model II are depicted in Figure 3. Water molecules that are initially outside the nanoporous membrane are not shown here to highlight the flow of water molecules from inside to outside of the membrane. A constant angular velocity ω is applied on the blade so that all molecules and ions rotate with the rotator. As time progresses, water molecules gradually flow out from inside the membrane to the outside through pores on the wall while Na + and Cl − ions are blocked inside because of the sub-nanometer pores. An appreciable distance of water molecules and ions to the membrane (∼0.2 nm) is observed, which is because of the hydrophobic property of the membrane material. Notably, the hydrophobicity of the membrane can be beneficial in two aspects. First, less friction of the membrane to water molecules enables a larger centrifugal pressure P ω . Second, there is little friction force between the hydrophobic graphene membrane and hydrated ions, so that it enables a fast boundary slip motion of the ions, which reduces the fouling of the membrane by the hydrated ions. The water density profile computed is shown in Figure 4a, in which the blue area (3.9 nm < r < 4.4 nm) indicates zero water density, corresponding to the impermeable wall occupied by carbon atoms. The three channels correspond to the three nanopores on the membrane. The layered structure of water molecules along the radial direction is observed here, which was also reported in other works e.g., [28]. Among different mechanisms, two of the possible causes for the layered structure of water molecule distribution are the molecular attraction of the membrane and the centrifugal force during rotational motion. The first peak is located near the membrane at 3.6 nm < r < 3.9 nm and 4.4 m < r < 4.7 nm, with magnitude ∼1.5 times larger than bulk water. The second peak is located at 2.8 nm < r <3.1 nm with magnitude ∼1.2 times larger than bulk. The three valleys between three peaks are located at 2.4 nm < r < 2.8 nm, 3.1 nm < r < 3.6 nm and 4.7 nm < r < 5 nm with magnitude ∼0.8. The peaks and valleys in the water density profile are noted to better describe the behavior of the nanofluid in this rotational device, because they show where the water molecules are most likely to be found inside the centrifuge. Notably, the density of water molecules increases in the radial direction due to the centrifugal force. Figure 4b,c display the trajectory of Na + and Cl − ions from the center to the nanoporous membrane. The spiral shape trajectories indicate that ions are affected by the rotational field. A centrifugal force is developed on these ions and water molecules because of the rotational fluid field which facilitates their radial motion. Some water molecules near pores might flow out of the container under this centrifugal force while all ions will continue rotating near the inner wall of the nanoporous membrane. As mentioned in the previous section, one of the main advantages of this system is the self-assisted cleaning mechanism by preventing the blockage of nanopores from ions. Although the centrifugal force exerted on ions near the entrance of pores can push them into the filtration channels and cause blockage, which is the main reason for fouling in classic RO methods, the unique tangential force in the rotating fluid field, which is mainly attributed to the Coriolis force, may flush these ions away from the surface of the membrane and hence the nano channels. This dynamic motion of ions will ensure the clearance of most nanopores and the continuous transport of water molecules from inside to outside through those active pores.

Azimuthal Velocity and Boundary Slipping
Snapshots of velocity vectors of water molecules for Model II are depicted in Figure 5a,b. The average distributions of the azimuthal velocity v θ along the radial direction for all water molecules of Model I and Model II are shown in Figure 5c,d respectively. During MD calculations, high angular velocities ω ranging from 0.1 rad/ns to 100 rad/ns are applied to rotators, which can hardly be achieved in real macro-scale experiments. However, with high angular velocity ω, MD calculations can speed up and more data can be collected for analysis to reduce statistic errors. A similar approach has been employed and verified by many other similar MD calculations to save computational resource [29,30]. In addition, this approach can be justified by our data from MD calculations and derivations from fluid dynamics that water flux scales linearly with net pressure. This means relations obtained from this work can be used to extrapolate results under lower angular velocities used in the real system.
Although rotators of Model I and Model II have the same tangential velocity at the inner radius r = R i under the same angular velocity ω with value ωR i , the rotating fields of water generated by them are entirely different, as shown in Figure 5c,d. Within the range of 0 < r < 2 nm, the azimuthal velocity v θ of Model I is zero (0 < r < 2 nm excluded from plot in Figure 5c) because nothing exists inside the cylinder rotator, while v θ of Model II increases almost linearly (Figure 5d) because any water molecules and ions within this range (0 < r < 2 nm) can be accelerated by ω of the blade rotator. Within the range of 2 nm < r < 4 nm, water molecules in both models have almost constant v θ , but different magnitudes under different angular velocities. In general, the azimuthal velocity v θ of Model II is much larger than that of Model I, and both are smaller than their tangential velocity at the inner radius, ωR i . Note that the azimuthal velocity near the CNT walls increases with larger angular velocity ω. Multiple rotor velocities were analyzed to observe this behavior. It is also notable that only the most extreme ω in Figure 5d leads to an azimuthal velocity close to 150 m/s, while in a Maxwell distribution the mean molecular velocity of H 2 O at 20 degrees centigrade is about 595 m/s [31].
A slipping parameter α i can be defined here to describe the rotating of fluid caused by the rotator: is the azimuthal velocity of fluids at the inner radius r = R i . From Figure 5c,d, we can get the slipping parameter α i for both Model I and Model II. Averaging data in range 3 nm < r < 4 nm for each curve in Figure 5c,d and plotting as a function of the angular velocity of rotators (as shown in Figure 6a), we find that α i = 0.16 for Model I and α i = 0.82 for Model II. The small value of α i in Model I means friction between water molecules and the surface of the carbon nanotube (served as a rotator in Model I) are very small, which agrees well with other studies on the slip properties of Graphene/CNT/GO membranes [7,32]. However, while most of these studies focused on the unimpeded water transport in basal CNT nanochannels, we discovered that a slight friction still exists between the porous graphene membrane and water molecules during rotational motions. As discussed later, a no-slip boundary condition is preferred (α i → 1) for fluid near the rotator, therefore, Model II is a better choice than Model I. Notably, changing the material of rotator in Model I to more hydrophilic material can also increase α i , such as GO and MoS 2 which show much lower slipping property than carbon nanotube.

Derivation of Centrifugal Pressure
Following the classic fluid dynamics theory, we can derive the centrifugal pressure P ω required to balance the osmotic pressure, as shown in the Method section. Since boundary conditions in Model I and Model II corresponds to the maximal and minimal slipping situation respectively, we can obtain the lower and upper limits of the centrifugal pressure P ω by substituting the slipping parameter α i into Equations (5) and (6), The upper and lower limits of the centrifugal pressure correspond to the two idealized boundary conditions from fluid dynamics theory, which is further discussed later. The results from MD calculations can be compared to the upper and lower theoretical limits of centrifugal pressure to better understand how fluid dynamics theory relates to the nanoscale. Note the water density ρ used here is 1.5 g/cm 3 , based on results in Figure 4a. Meanwhile, the centrifugal pressure P ω can be obtained from MD calculations by summing up the radial force (F ir ) of all water molecules within 0.3 nm (thickness of the first water layer) of the nanoporous membrane and then divided by the area of the membrane (A m ). Centrifugal pressure obtained from both fluid dynamics and MD simulation are plotted in Figure 6b with α i = 0.82 from Model II.
The black and red lines in Figure 6b correspond to the lower and upper limits of the centrifugal pressure P ω derived from classic fluid dynamics. MD results are plotted in blue dots. From Figure 5d, the value of v θ is almost a constant from the inner radius r = R i to the outer radius r = R o under the same angular velocity ω, which means the boundary condition for fluid near the nanoporous membrane is almost ideal-slip (with outer slipping parameter α o ≈ α i ). Therefore, it is expected that the centrifugal pressure P ω acquired from MD calculations (notated as P MD ω ) would be close to the upper limit P  Figure 6b, we discover MD results agree well with the upper limit at relatively smaller angular velocity ω but start to deviate when ω becomes larger. The deviation is probably from the Taylor instability (such as turbulence) which will be discussed later. With these results, we conclude that the modified fluid dynamics theory (by incorporating slipping property of boundaries) is still valid at the nanoscale before the fluid becomes unstable under large ω.
From Figure 6b, we can also obtain the critical angular velocity ω cr with value 32.55 rad/ns (the vertical dashed line) needed to balance the osmotic pressure with value Π = 3 MPa (the horizontal dashed line), which is the threshold value for water filtration. Water molecules can flow out from pores on the nanoporous membrane only when the angular velocity ω of the rotator in Model II is higher than ω cr .

Taylor Instability Due to High Rotational Flux
Higher water flux can be reached by increasing the angular velocity ω of the rotator. However, we may encounter other problems as ω becomes too large. First, a larger centrifugal force might be developed. This force can push ions more in-depth into the nano-channels causing blockage or even damage to the membrane material. It is also possible that ions are "squeezed out" from the inside leading to degradation of the rejection rate. Second, the Taylor-Couette instability can be developed when Reynolds' number becomes significant. We note that once the angular velocity of the rotator exceeds a specific value, hydrodynamic instability and turbulence will occur, which may affect the efficiency of the desalination process, as shown in Figure 7.
Two typical vortexes symmetric to the central plane can be seen inside the model when ω ≥ 70 rad/ns. Figure 7a shows the shape of the vortex and Figure 7b shows the density of water molecules along radial and longitudinal directions, with the azimuthal direction averaged out. The dark blue regions correspond to the two vortexes in Figure 7a. A more detailed analysis of water density is shown in Figure 7c, a quadratic shape is observed in the vortex area when both azimuthal and longitudinal directions are averaged out. Outside the vortex area, water density increases linearly to around 1.2 g/cm 3 , and then follows the layered structure when approaching the membrane, as discussed in Figure 4a. However, the width of the layered water structure here is much smaller than that in Figure 4a, because of the squeezing effect from large angular velocity ω. Further study is needed to investigate the phenomenon of Taylor instability at the nanoscale.

Discussion
One of the most significant features of the present system regards the angular velocity dependence of water flux. Figure 8a shows the total volume of water molecules flowing out from inside of the container as a function of simulation time (assuming one water molecule occupies a volume ≈3 × 10 −2 nm 3 ). Linear interpolating of the three curves under three angular velocity ω is shown in dashed lines for the calculation of water flux. When ω < ω cr , the centrifugal pressure P ω is smaller than the osmotic pressure Π, leading to an inward net flow of water molecules (corresponds to the negative slope of black curves). When ω > ω cr , the net water flow becomes outward (positive slopes of red and blue curves) and has a larger slope when ω is larger. Water flux (water flow per time) are the slopes of dashed lines in Figure 8a, which are plotted as a function of angular velocities ω in Figure 8b. Water flux predicted from MD calculations is a quadratic function of angular velocity ω, which agrees well with Equation (9) derived from fluid dynamics. Recall that pressure is also quadratically proportional to ω (Figure 6b), therefore we can conclude that water flux is linearly proportional to the overall centrifugal pressure as reported in other studies [32]. The centrifugal pressure in their works is from hydrostatic pressing while in our work it is from the centrifugal motion of the rotational fluid. Since we already proved fluid dynamics theory is still valid at the nanoscale, MD data in Figure 8b are fitted with Equation (9) to obtain the permeability of this system, which is 60 (L/cm 2 /day/MPa). A comparison between the permeability of the proposed device and the experimentally measured values from Pendergast et al. [33] and Grossman et al. [18,19,26] is shown below in Figure 9a. The water permeability of the proposed system is much better than three of them and comparable to the best two materials, which means that the desalination efficiency of this system could be 3-4 orders of magnitude higher than that of commercially available RO devices. Energy efficiency of this system is also estimated with the relation: η = Energy Input Filtrated Water . Energy Input is defined as the kinetic energy change of the system within a time period ∆t, and Filtrated Water is defined as the total number of water molecules passing through the nanoporous walls within the same time ∆t. A longer duration MD simulation (>30 ns) is conducted with applied angular velocity ω = 35 rad/ns. Figure 9b shows the value of energy efficiency as a function of simulation time. With the given angular velocity, Energy input = 2.045 × 10 −18 J. The Filtrated Water at a specific time t 0 can be obtained by integrating the water flux shown in Figure 9a from initial time to t 0 . For example, Filtrated Water = 296.13 nm 3 at the end of 30 ns. With relation: 1 J = 2.78 ×10 −7 kWh, the energy efficiency is: 1.92 kWh/m 3 , which falls below the lower limit of traditional reverse osmosis RO (2∼5 kWh/m 3 ) [34].
Our results indicate that this new nanopore membrane structure can act as a high-permeability and low energy consumption membrane for water filtration. We also provided guidelines for the design of macroscale nanoporous membrane co-axial centrifuge. Due to computational resource limitations, we only present a MD simulation of a nanoscale membrane centrifuge of 8 nm in diameter. However, the radius of the centrifuge may not affect the simulation result qualitatively, as long as the nanopores are created on the co-axial centrifuge wall, because we have no obvious size effect for such nanofluidic devices. From the theoretical perspective, if the size of the nanoporous membrane centrifuge reaches to macroscale, the critical angular velocity needed will be much smaller, which indicates that it is possible to fabricate such desalination centrifuge at engineering scale. Moreover, we can increase the water permeability by increasing the porosity of the membrane centrifuge, i.e., increasing the pore density on the membrane surface. Therefore, we propose a macroscale centrifuge that is decorated with nanoporous membrane to solve the scaling up problem existed in present nanomaterial based desalination technology.

MD Simulation Details
The open source molecular dynamics software package, Gromacs [35], is used for all the MD calculations. The OPLS-aa force field [36] is used to describe bonded and nonbonded interactions between atoms, ions and molecules. The particle mesh Ewald method [37] is used for calculating long range electrostatic interactions, with cut-off distance of 0.8 nm, as is used in other similar MD simulations [25,27,38]. The SPC/E water model was used due to its advantage of simplicity, computational efficiency, and speed [39]. This water model is rigid and non-polarizable consisting of three interaction sites, located on the oxygen and hydrogen atomic centers, which shows the best overall agreement with experiments [40]. Periodic boundary conditions in all directions are employed. A total simulation time of 6 ns is conducted with data in the last 2 ns collected for data analysis. The integration time step is set to be 1 fs with a leap-frog algorithm [41]. The neighbor list is updated every time step to avoid intrinsic errors [42]. The MD system is first relaxed to minimize its free energy in a duration of 50 ps, followed by an NVT calculation by using the Nose-Hoover thermostat [43] to maintain a temperature at 298 K with damping coefficient of 0.1 ps. Then the non-equilibrium molecular dynamics simulation (NEMD) [29,44] featured with the rotating blade is then performed with the same Nose-Hoover thermostat still applied to the system. The rotating blade is prescribed to rotate around the central axis of the centrifuge rotator with a given angular velocity, and a tangential force perpendicular to the radial direction of the rotating blade may be exerted to each atom inside the centrifuge. Ten different values of the angular velocities ranging from 0.1 to 100 rad/ns are selected in the simulations. Even though the high angular velocities used in the simulations are almost practically impossible, however, since MD calculations are limited by their time scale, large angular velocities were required in order to satisfy the critical angular velocity condition to generate a sufficient pressure to balance the osmosis pressure. For practical desalination applications of the centrifuge, the readers are referred to the discussions on how to scale-up such desalination centrifuge.

Feed Flow Motion Solution Derived from Fluid Dynamics Theory
From a continuum point of view, rotational fluid motion inside the model in Figure 2 can be considered as slipped Taylor-Couette flow in a confined nano-size space. The complex interactions of water molecules and hydrated ions near pores of the nanoporous membrane can be viewed as competitions of different forces: the centrifugal force due to internal rotation of the rotator, the osmosis force due to the solute concentration difference, and the Coriolis force due to the inertia effect. In addition to MD calculations, the critical centrifugal pressure, P c , that balances the osmotic pressure can be derived from classical fluid dynamics.
According to water transport theory in membranes, the net water flux J (in unit cm/s) of a porous membrane may be expressed as: Here, N p is the number of nano-sized pores at the wall of the nanoporous membrane. The first term on the right side is the water flux generated by centrifugal pressure (notated as the centrifugal pressure P ω since it depends on the angular velocity ω in our model), with L p (in unit cm 3 /N·s) defined as hydraulic permeability. The second term is the water flux generated by osmotic pressure (notated as Π), with ∆C defined as the solute concentration gradient and A c the osmotic permeability (in unit cm/s). V H 2 O is the molar volume of pure water. For seawater with salt concentration around 35 g/L at room temperature, ∆C = 0.6 M and V H 2 O = 18.91 cm 3 /mol.
Since the salt solution in the present study is relatively dilute, the osmotic pressure is related to ∆C through the linear equation: Π = 2R v T∆C = 3 MPa (R v is the gas constant and T = 300 K is the environmental temperature. The double value comes from pairs of Na + and Cl − ions in saltwater). The relation between osmotic permeability A c and hydraulic permeability L p can be derived according to the equilibrium condition (zero net water flux without centrifugal pressure: Substituting this relation into the previous water flux equation, we can obtain the relation of hydraulic permeability as a function of net freshwater flux: In this work, the centrifugal pressure P ω near pores of the membrane is generated by a rotational field, which can be derived by incorporating a classic Newtonian Couette flow under an angular velocity ω from the rotator. Consider the Navier-Stokes equation in fluid dynamics: where ρ is the density and µ is dynamic viscosity. v and P are the velocity vector field and pressure field of the fluid respectively. Since our main focus in the current study is a first order approximation of pressure near the membrane wall, it is reasonable to assume that fluid motions in longitudinal and azimuthal directions are homogeneous without gradients of both velocities and pressure, which means ∂v ∂z = ∂v ∂θ = 0, and ∂P ∂z In addition, since the angular velocity ω is large at the nanoscale, velocities in the longitudinal (v z ) and radial (v r ) directions are negligible in comparison with the velocity in azimuthal direction (v θ ), which is usually two to three orders of magnitude larger. Therefore, the velocity and pressure can be approximated as: v = [v z = v r = 0, v θ (r)], P = P(r). With these assumptions and considering only the steady state solution, Equation (2) can be reduced to two ODEs: The general solutions for the azimuthal velocity v θ and pressure P of the fluid are: where c 1 and c 2 are constants depending on boundary conditions (BCs). R i is the radius of inner-layer rotator where angular velocity ω is applied. P(R i ) is set to be the atmospheric pressure (P(R i ) = P 0 ≈ 1 Bar because we are assuming the centrifuge is in equilibrium with the external environment). Although the azimuthal velocities of the rotator at r = R i and the nanoporous membrane at r = R o are ωR i and 0 respectively, these values are not necessarily equal to velocities of fluid near to the two boundaries. The slipping phenomenon for materials at nanoscale has been studied extensively [45], and different materials show different viscous and slipping properties [32,46]. Therefore, a detailed investigation is needed in MD calculation in order to quantify the slipping effect.
Here a generalized situation is discussed: Assuming the azimuthal velocities of fluid near the inner radius r = R i and the outer radius Where α i and α 0 are the slipping parameters of the inner and outer radius of the device, respectively. Note α o ≤ α i since no external force is applied on fluid between R i and R o and v θ (r) will decay from R i to R o because of viscosity. Obviously, α i = 1, α o = 0 correspond to the classic no-slip BCs and α i = α o = 0 correspond to the ideal-slip BCs. Therefore, the range of the two slipping parameters are: For energy efficiency, a non-slip BC for fluids near r = R i (α i = 1) and an ideal-slip BC for fluids near r = R o (α o = α i ) are preferred for best performance of the system, since v θ (R i ) reaches maximal by the friction of the rotator and v θ (R o ) is not affected by the friction of the static nanoporous membrane. Although neither the non-slip nor the ideal-slip conditions can be achieved for real materials, different materials can be selected to approach a limiting value.
The solutions of v θ and P can be obtained as functions of α i and α o by substituting v θ (R i ) and v θ (R o ) into Equation (4). Here two types of BCs are listed explicitly for further discussion: Boundary conditions case 1: where The pressure at the boundary of the cylinder P(r = R o ) can be solved because P ω = P(R o ), after obtaining α i from MD calculation. Here boundary conditions Case 1 and Case 2 correspond to the lower and upper bounds of v θ and P. Therefore, these two types of BCs define the limiting performance that a given system can achieve. Moreover, a comparison between solutions from Equation (4) and MD results can quantify the derived analytical results from fluid dynamics for the nanofluidics systems.

Comparisons between MD and Fluid Dynamics Theory
Within fluid dynamics theory, the expression of pressure in Equations (5) and (6) can be further simplified as is a function of geometry or the dimension (size) of the centrifuge that varies depending on different boundary condition, and it has the unit of L 2 . In the case of the first boundary condition, we have and in the case of the second boundary condition, we have Thus, the parameter B is proportional to the radius of the centrifuge, i.e., B ∼ R 2 o . The centrifugal pressure P ω can be substituted into Equation (1) to obtain the expression of net water flux J as a function of angular velocity ω shown below. Meanwhile, a minimum or critical angular velocity ω cr can be derived by setting J = 0, i.e., Within MD calculations, slipping parameter α i can be obtained by investigating the velocity profiles of water molecules. The net water flux J can be obtained by simply counting the number of water molecules (N H 2 O ) passing through the nanoporous membrane within the simulation time, where N A is Avogadro's number; A m is the area of the membrane, and ∆t is the simulation time). Stable J can be obtained for each ω when ∆t is big enough. The relationship between J and angular velocity ω can be plotted and directly compared with Equation (9). The critical angular velocity, ω cr , derived from Equation (9), can also be compared with the applied ω in the MD calculation when J = 0. Velocities obtained from MD calculations can also be compared to the analytical velocity result, v θ (r), derived from Equation (4). However, the comparison of the centrifugal pressure between MD calculations and Equations (5) and (6) needs further analysis, because the pressure defined in the MD simulation is for the whole system, instead of the value at r = R o . Nonetheless, the original definition of pressure (force per area) can be used to derive the centrifugal pressure P ω based on MD theory, with the expression where F ir is the radial component of the force vector exerted on water molecule i. Here, the pressure calculated from the sum of the total force experienced on each water molecule in the direction of the membrane divided by the area of the membrane is purely based on Newton's laws of motion. By comparing the centrifugal pressure P ω from MD calculations with values obtained from Equations (5) and (6), it is possible to investigate whether or not P ω at micro-scale is still within the bounds defined at macro-scale.
In general, by comparing results from two independent theories, it would be possible to analyze the size effect and discuss how the performance of the system varies if the size of the system is scaled up from micro-scale to macro-scale.
In the ensuing derivation, the following assumptions are made, as discussed in the main text of the paper, 3. Velocity in θ direction is much higher than z and r directions, i.e., ν θ >> ν z , ν r , therefore it is assumed that ν z = ν r ≈ 0.
Considering these assumptions, from Equation (A1), we obtain, Considering these assumptions in Equation (A2), we then obtain, Equations (A3) and (A4) can be cast as Equation (A5), which is the same as Equation (3) in the main text,

. Derivation of General Solutions (4) from N-S Equation
The general solution of ν θ in Equation (A5) is, where the constants c 1 and c 2 are defined from the boundary conditions. Substituting ν θ into equation, dP dr = P ν 2 θ r , and integrating over the whole domain leads to, Appendix A.3. Solution of Equation (4) under Two Special Boundary Conditions Equations (5) and (6) I. From the boundary conditions: ν θ (R i ) = α i R i ω i and ν θ (R 0 ) = 0, we have Then c 1 and c 2 can be solved as, where Then ν θ (r) can be rewritten as, When r = R 0 , the pressure P ω is given as, (A11) II. For the boundary conditions: ν θ (R i ) = ν θ (R 0 ) = α i R i ω, we have, Then c 1 and c 2 can be solved as, where Then ν θ (r) can be rewritten as, When r = R 0 , the pressure P ω is given as, Appendix A.4. Derivation of the Critical Angular Velocity, ω cr : Consider the fact that When J(ω cr ) = 0, it leads to, Finally, ω cr can be solved as follows, where B = B(R i , R o ) is a function of dimension of the centrifuge that varies under different boundary conditions, and it has units of L 2 . For the first type of boundary conditions, we can derive that and the second set of boundary conditions, we have, where 0 < r = R i /R o ≤ 1. Therefore, in general, B ∼ C(α i , β i , r)R 2 o , and C is a function of the inner radius of the centrifuge and the boundary condition of the centrifuge.

Appendix B. Scale-Up Idea to Fabricate a Macroscale Nanoporous Rotating Centrifuge
In Appendix B, we display a graphic illustration on how to scale up the proposed rotational nanofluidic device with nanoporous membrane wall. The illustration demonstrate both the scale-up concept as well as the detailed fabrication steps in Figure A1. The key to fabricate such scale up rotating nano-porous graphene membrane centrifuge is to create a multiscale porous structure. The first step is to drill micron scale tine holes on the metal wall or other structural material wall that can provide structure support for a macroscale system with its dimension R o > 1 m. The second step is to paste small nano-porous graphene patches to cover those micron scale holes as illustrated in Figure A1. By doing so, we can achieve: 1. having a nano-porous membrane on the wall of the centrifuge; 2. having a strong wall structure to support macroscale loads such as overall pressure under working conditions; 3. having a realistic critical angular velocity of a few hundreds of radian per second to generate a centrifugal pressure P ω that is greater than the osmosis pressure.