1887
banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Nonequilibrium molecular dynamics simulation of water transport through carbon nanotube membranes at low pressurea)
a)This paper was presented, in part, at NAMS’11, Las Vegas, NV, June 2011.
Rent:
Rent this article for
USD
10.1063/1.4734484
/content/aip/journal/jcp/137/4/10.1063/1.4734484
http://aip.metastore.ingenta.com/content/aip/journal/jcp/137/4/10.1063/1.4734484

Figures

Image of FIG. 1.
FIG. 1.

Orthographic snapshots of the CNT membrane model produced by the VMD package for the MD simulation: (a) Top view down through the CNT membrane model and (b) side view of the CNT membrane model. Two graphene sheets are located at the ends of the CNT and each sheet has an approximated hole in the center to match the pore size of the CNT.

Image of FIG. 2.
FIG. 2.

Perspective snapshots of the simulation systems produced by VMD package at the beginning status of MD simulations. Shown are the two graphene sheets which act as the movable walls, the CNT membrane model shown in Fig. 1, and the two water reservoirs connected by the membrane. Forces f t and f b act by the movable walls on the fluid. Carbon atoms in green, hydrogen atoms in white, and oxygen atoms in red.

Image of FIG. 3.
FIG. 3.

Pressure produced in each water reservoir as a function of the simulation time (10 ns as an example shown here, real simulation time is longer and dependent on each case). Each pressure value is averaged in 10 ps after the pressure difference stabilizes in the NEMD simulation. The error about the mean value of the pressure is less than 0.2% for the high pressure on the top reservoir and is less than 10% for the low pressure on the bottom reservoir. The standard deviation for the pressure differences over all seven simulations is about 0.2 MPa.

Image of FIG. 4.
FIG. 4.

Number flow rate of #4 CNT membranes in 30 ns simulations, shown as the slope of the linear trend in the net number of water molecules passing through the membrane as a function of the simulation time. The data shown are for membrane simulations of CNT #4 (Table I): (12, 12) CNT; L z = 6 nm. ΔP = 5.0 MPa; T = 300 K. The standard error of the slope is ±0.05 molecules/nm.

Image of FIG. 5.
FIG. 5.

Flux (rate per membrane area), shown both as molecular flux and conventional flux (in m/s), as a function of the membrane pore size as represented by the pore area. The points are based on membrane simulations on CNT #1–4 (Table I): #1 (8, 8), #2 (10, 10), #3 (11, 11), and #4 (12, 12) CNT, respectively; L z = 6 nm. ΔP = 5.0 MPa; T = 300 K. H-P equation refers to the Hagen-Poiseuille equation (Eq. (1)), and solid line is the trend based on the simulation results. The statistical error for each data point is evaluated and less than 0.5% (see Fig. 4 for the case of the #4 membrane).

Image of FIG. 6.
FIG. 6.

Effect of pore size as represented by the pore radius on the density distributions along radial direction for the four membrane models: density as a function of r; values are averaged over a thin annular section in the pore at r (excluding the entrance and exit regions). The points are based on membrane simulations on CNT #1–4 (Table I): #1 (8, 8), #2 (10, 10), #3 (11, 11), and #4 (12, 12) CNT, respectively; L z = 6 nm. ΔP = 5.0 MPa; T = 300 K. Solid lines are the trend curves based on the simulation results. The statistical error for each data point is less than 1%.

Image of FIG. 7.
FIG. 7.

Orthographic snapshots of water molecules in the different CNT membranes. Atoms are shown by balls with the van der Waals volumes: carbon atoms in green, hydrogen atoms in white, and oxygen atoms in red.

Image of FIG. 8.
FIG. 8.

Effect of pore size as represented by the pore radius on the distributions of z component velocity along the radial direction for the four membrane models: velocity as a function of r; values are averaged over a thin annular section in the pore at r (excluding the entrance and exit regions). The points are based on membrane simulations on CNT #1–4 (Table I): #1 (8, 8), #2 (10, 10), #3 (11, 11), and #4 (12, 12) CNT, respectively; L z = 6 nm. ΔP = 5.0 MPa; T = 300 K. N-S equation refers to the Navier-Stokes equation (Eq. (2)). Dashed lines are the trends based on the Navier-Stokes equation, and solid lines are the trend curves based on the simulation results. The statistical errors for most data points are less than 10%.

Image of FIG. 9.
FIG. 9.

Effect of pore size on the density distributions along the flow direction for the four membrane models: density as a function of z; values are averaged over a thin section at z (cylindrical sections in the pore and cuboid sections in the water reservoirs). The points are based on membrane simulations on CNT #1–4 (Table I): #1 (8, 8), #2 (10, 10), #3 (11, 11), and #4 (12, 12) CNT, respectively; L z = 6 nm. ΔP = 5.0 MPa; T = 300 K. Dashed lines are the membrane boundaries along z direction separating the system in three parts from left to right: the water reservoir at high pressure side, the pore of the CNT membrane, and the water reservoir at low pressure side. Solid lines are the trend curves based on the simulation results. The statistical error for each data point is less than 1%.

Image of FIG. 10.
FIG. 10.

Effect of pore size on the velocity (in z direction) distributions along the flow direction for the four membrane models: velocity as a function of z; values are averaged over a thin section at z (cylindrical sections in the pore and square cuboids sections in the water reservoirs). The points are based on membrane simulations on CNT #1–4 (Table I): #1 (8, 8), #2 (10, 10), #3 (11, 11), and #4 (12, 12) CNT, respectively; L z = 6 nm. ΔP = 5.0 MPa; T = 300 K. Dashed lines are the membrane boundaries along z direction separating the system in three parts from left to right: the water reservoir at high pressure side, the pore of the CNT membrane, and the water reservoir at low pressure side. Solid lines are the trend curves based on the simulation results. The statistical error for each data point is less than 15%.

Image of FIG. 11.
FIG. 11.

Flux (rate per membrane area), shown both as molecular flux and conventional flux (in m/s), as a function of the inverse membrane thickness. The points are based on membrane simulations on CNT #5, #4, #6, and #7 (Table I) corresponding to different lengths of 3 nm, 6 nm, 9 nm, and 15 nm, respectively, and r p = 0.643 nm. ΔP = 5.0 MPa; T = 300 K. H-P equation refers to the Hagen-Poiseuille equation (Eq. (1)), and the solid line is the trend based on the simulation results. The statistical error for each data point is less than 0.5% (see Fig. 4 for the case of #4 membrane).

Image of FIG. 12.
FIG. 12.

Effect of membrane thickness on the distributions along the radial direction for the four membrane models; values are averaged over a thin annular section in the pore at r (excluding the entrance and exit regions): (a) Density as a function of r; (b) Velocity, in z direction, as a function of r. The points are based on membrane simulations on CNT #5, #4, #6, and #7 (Table I) corresponding to (12, 12) CNT with different lengths of 3 nm, 6 nm, 9 nm, and 15 nm, respectively, and r p = 0.643 nm. ΔP = 5.0 MPa; T = 300 K. N-S equation refers to the Navier-Stokes equation (Eq. (2)). Dashed lines are the trends based on the Navier-Stokes equation, and solid lines are the trend curves based on the simulation results. The statistical error for each density data point is less than 1% and the statistical errors for most data points are less than 10%.

Image of FIG. 13.
FIG. 13.

Effect of membrane thickness on the distributions along flow direction for the four membrane models; values are averaged over a thin section at z (cylindrical sections in the pore and square cuboids sections in the water reservoirs): (a) density distributions as a function of z and (b) velocity, in z direction, as a function of z. The points are based on membrane simulations on CNT #5, #4, #6, and #7 (Table I) corresponding to (12, 12) CNT with different length of 3 nm, 6 nm, 9 nm, and 15 nm, respectively, and r p = 0.643 nm. ΔP = 5.0 MPa; T = 300 K. Dashed lines are the membrane boundaries along the z direction separating the system in three parts, from left to right: the water reservoir at high pressure side, the pore of the CNT membrane, and the water reservoir at low pressure side. Solid lines are the trend curves based on the simulation results. The statistical error for each density data point and each velocity data point are less than 1%, and less than 15%, respectively.

Tables

Generic image for table
Table I.

Specifications of different CNT membrane physical parameters used in the simulations.

Generic image for table
Table II.

Comparison of the enhancement factorsa between our NEMD simulation and the experimental data of Qin.45

Generic image for table
Table III.

Enhancement factorsa of the #4–7 CNT membranes at 5.0 MPa pressure difference.

Loading

Article metrics loading...

/content/aip/journal/jcp/137/4/10.1063/1.4734484
2012-07-23
2014-04-25
Loading

Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Nonequilibrium molecular dynamics simulation of water transport through carbon nanotube membranes at low pressurea)
http://aip.metastore.ingenta.com/content/aip/journal/jcp/137/4/10.1063/1.4734484
10.1063/1.4734484
SEARCH_EXPAND_ITEM