arxiv: v1 [physics.comp-ph] 17 Dec 2018

Transcription

1 Pressures nsde a nano-porous medum. The case of a snge phase fud arxv: v1 [physcs.comp-ph] 17 Dec 2018 Oav Gateand, Dck Bedeaux, and Sgne Kjestrup PoreLab, Department of Chemstry, Norwegan Unversty of Scence and Technoogy, NTNU December 18, 2018 Abstract We defne the pressure of a porous medum n terms of the grand potenta, and compute ts vaue n a nano-confned or nano-porous medum, meanng a medum where thermodynamc equatons need be adjusted for smaness. On the nano-scae, the pressure depends n a cruca way on the sze and shape of the pores. Accordng to H [1], two pressures are needed to characterze ths stuaton; the ntegra pressure and the dfferenta pressure. Usng H s formasm for a nano-porous medum, we derve an expresson for the dfference between the ntegra and the dfferenta pressures n a spherca phase α of radus R, ˆp α p α = γ/r. We recover the aw of Young-Lapace for the dfferenta pressure dfference across the same curved surface. We dscuss the defnton of a representatve voume eement for the nano-porous medum and show that the smaest REV s haf a unt ce n the drecton of the pore n the fcc attce. We aso show, for the frst tme, how the pressure profe through a nano-porous medum can be defned and computed away from equbrum. 1 Introducton The descrpton of transport processes n porous meda poses many chaenges, we descrbed n the terature, see e.g. [2 6]. There s, for nstance, no consensus, nether on the defnton, nor on the measurement or the cacuaton of the pressure n a porous medum wth fow of mmscbe fuds. The probem wth the defned mcroscopc pressure tensor [5, 7] s accentuated n a heterogeneous system wth nterfaces between sods and fuds. In a homogeneous fud phase one may defne and cacuate a pressure and a pressure gradent from the equaton of state. In a porous medum the presence of curved surfaces and fud confnements makes t dffcut to appy accepted methods for cacuaton of the mcroscopc pressure tensor and, consequenty, the pressure gradent as drvng force for fud fow. The scae at whch we choose to work w be decsve for the answer. Moreover, the scae that the hydrodynamc equatons of transport refer to, remans to be gven for nano-porous as we as mcro-porous meda. A centra eement n the dervaton of the equatons of transport on the macro-scae s the defnton of a representatve voume eement (REV), see e.g. [8, 9]. The sze of the REV shoud be arge compared to the pore-sze and sma compared to sze of the porous medum. It shoud contan a statstcay representatve coecton of pores. We have recenty dscussed [10] a new scheme to defne as bass set of addtve varabes: the nterna energy, entropy, and masses of a the components of the REV. These varabes are addtve n the sense that they are sums of contrbutons of a phases, nterfaces and contact nes wthn the REV. Usng Euer homogenety 1

2 of the frst knd, we were abe to derve the Gbbs equaton for the REV. Ths equaton defnes the temperature, pressure and chemca potentas of the REV as parta dervatves of the nterna energy of the REV [10]. As dscussed n [11] the grand potenta, Υ, of the REV s gven by mnus k B T tmes the ogarthm of the grand partton functon, Z g, where k B s Botzmann s constant and T s the temperature. The grand potenta s equa to mnus the contrbuton to the nterna energy from the pressure-voume term, k B T n Z g = Υ = pv, whch we w from now on refer to as the compressona energy. For a snge fud f n a porous medum r, the resut was [10, 11] pv = p f V f + p r V r γ fr Ω fr, (1) where p and V are the pressure and the voume of the REV. Furthermore p f and V f are the the pressure and the voume of the fud n the REV, p r and V r are the the pressure and the voume n the grans n the REV, and γ fr and Ω fr are the surface tenson and the surface area between the fud and the gran. The assumpton behnd the expresson was the addtve nature of the grand potenta. Ths defnton of the REV, and the expresson for the grand potenta opens up a possbty to defne the pressure on the hydrodynamc scae. The am of ths work s to expore ths possbty. We sha fnd that t w work very we for fow of a snge fud n a porous medum. As a non-mtng ustratve exampe, we use grans postoned n a fcc attce. The work can be seen as a contnuaton of our earer works [10, 11]. The work so far consdered transport processes n mcro-porous, not nano-porous meda. In mcro-porous meda, the pressure of any phase (the surface tenson of any nterface) s ndependent of the voume of the phase (the area between the phases). Ths was cruca for vadty of equaton 1. For nano-porous systems, we need to step away from equaton 1. Foowng H s procedure for sma systems thermodynamcs [1], we generaze equaton 1 to provde an expresson for the thermodynamc pressure n a nano-porous medum. We sha see that not ony one, but two pressures are needed to hande the addtona compcatons that arse at the nano-scae; the mpact of confnement and of rad of curvature of the nterfaces. In the thermodynamc mt, the approach presented for the nano-scae must smpfy to the one for the macro-scae. We sha see that ths s so. In order to work wth controed condtons, we w frst nvestgate the pressure of a fud around a snge sod nano-scae gran and next around a attce of sod nano-scae grans. The new expresson, whch we propose as a defnton of the pressure n a nano-porous medum, w be nvestgated for vabty and vadty for ths case. The present work can be seen as a frst step n the drecton towards a defnton and use of pressure and pressure gradents n rea porous meda. The pressure s not unquey defned at moecuar scae. Ths ack of unqueness becomes apparent n moecuar dynamcs (MD) smuatons, for whch the compuatona agorthm has to be carefuy desgned [7]. The predomnant method for pressure cacuatons n partcuar systems s usng the Irvng-Krkwood contour for the force between two partces [12]. Ths agorthm works for homogeneous systems, but speca care must be taken for heterogeneous systems [5, 6]. However, f the contro voume (REV) used for pressure cacuaton s arge compared wth the heterogenety ength scae, one may argue that the the agorthm for homogeneous systems gves a good approxmaton to the true resut. We are nterested n the sotropc pressure averaged over the REV, on a scae where the porous medum can be consdered to be homogeneous. The paper s organzed as foows. In secton 2 we derve the pressure of a REV for one sod gran surrounded by fud partces (Case I) and for a three-dmensona face-centered cubc (fcc) attce of sod grans (Case II). Secton 3 descrbes the moecuar dynamcs smuaton technque when the system s n equbrum and n a pressure gradent. In secton 4 we use the theory to nterpret resuts of equbrum moecuar dynamcs smuatons for one sod gran and for an array of sod grans n a fud. Fnay we appy the resuts to descrbe the system under 2

3 a pressure gradent. We concude n the ast secton that the expressons and the procedure deveoped provde a vabe defnton of the pressures and pressure gradents n nano-porous meda. 2 The pressure of a nano-porous medum Equaton 1 appes to a mcro-porous medum, a medum where the pore-sze s n the mcrometer range or arger [10,11]. For a nano-porous medum we need to appy the thermodynamcs of sma systems [1]. In nano-porous meda, ths technque s therefore we suted for the nvestgaton. The thermodynamc propertes ke nterna energy, entropy and masses of components of a sma system are not proportona to the system s voume. As H expaned, ths eads to the defnton of two dfferent pressures, for whch he ntroduced the names ntegra and dfferenta pressure, ˆp and p, respectvey. For a system wth a voume V, these pressures are reated by p(v ) = (ˆp (V ) V ) V = ˆp (V ) + V (ˆp (V )). (2) V The symbo p (the dfferenta pressure) s gven to the varabe that we normay understand as the pressure on the macroscopc eve. It s ony when ˆp depends on V, that the two pressures are dfferent. For arge systems, ˆp does not depend on V and the two pressures are the same. The ntegra and dfferenta pressures connect to dfferent types of mechanca work on the ensembe of sma systems. The dfferenta pressure tmes the change of the sma system voume s the work done on the surroundngs by ths voume change. The name dfferenta derves from the use of a dfferenta voume. Ths work s the same, whether the system s arge or sma. The ntegra pressure tmes the voume per repca s, however, the work done by addng one sma system of constant voume to the remanng ones, keepng the temperature constant. Ths work s speca for sma systems. It derves from the ensembe vew, but s equay we measurabe. The word ntegra derves from the addton of a sma system. From statstca mechancs for macro-scae systems, we know that pv = ˆpV equas k B T tmes the natura ogarthm of the grand-canonca partton functon. For a sma (nano-szed) system, H ( [1], equaton 1-17), showed that ths ogarthm gves ˆpV. In nano-porous meda ths product s dfferent from pv, cf. equaton 2. Energes are st addtve and the tota compressona energy wthn the sma system s smar to equaton 1. We repace equaton 1 by: ˆpV = ˆp f V f + ˆp r V r ˆγ fr Ω fr, (3) where ˆp f, ˆp r are ntegra pressures of the sub-voumes V f and V r, and ˆγ fr s the ntegra surface tenson. We consder here a nano-porous medum, so ntegra pressures and ntegra surface tensons appy. The ntegra pressure and ntegra surface tenson normay depend on the system sze. In the porous medum there are two characterstc szes: the sze of a gran and the dstance between the surfaces of two grans. 1 The quanttes ˆp, ˆp f, ˆp r and ˆγ fr may depend on both. We sha here examne a system (cf. secton 3) of spherca, monodsperse grans, for whch the radus R s a good measure of the sze. The voume of the grans may be a good aternatve measure, whch we w aso use. The dependence on the gran sze and on the dstance between the surfaces of the grans w be studed n an effort to estabsh equatons 1 and 3. In the foowng, we consder a snge spherca gran confned by a snge phase fud (Case I) and a face-centered cubc (fcc) attce of spherca grans confned by a snge phase fud (Case 1 Another vad characterstc sze s the sze of the pores between the grans, but ths foows from the two we have chosen. 3

4 Fgure 1: A partce n a confned system (A) n equbrum wth a buk fud phase (B). II). The sze of the REV does not need to be arge, and we w show n secton 4.2 that the smaest REV s haf a unt ce n the drecton of the pore n the fcc attce. 2.1 Case I. Snge spherca gran Consder the spherca ncuson of a gran r n a box wth a fud phase f. Ths s system A of fgure 1. Phase f has the voume V f and phase r has the voume V r. The tota voume s V = V f + V r. The surface area between phases f and r s Ω fr. System A s n equbrum wth system B, whch has the same voume V, and contans fud f ony. We now make the assumptons that ˆp r depends ony on V r, whe ˆp, ˆp f and ˆγ fr do not depend on V r or any other sze varabe. Ths means that the spherca gran r s consdered to be sma n H s sense, whe the fud phase s a arge thermodynamc system. When the surface tenson depends on the curvature, there s a dependence of ˆγ fr on Ω fr [13, 14]. Ths nterestng effect, whch we w not consder here, becomes reevant as the gran sze decreases. The assumptons w be used to fnd the contrbutons to equaton 3. The assumptons mean that Equaton 3 for system A n fgure 1 reduces to ˆp r p r, ˆp = p, ˆp f = p f and ˆγ fr = γ fr. (4) pv = p f V f + ˆp r V r γ fr Ω fr. (5) An nterestng mpcaton of ˆp f = p f s that we can use the standard MD agorthms as mpemented for arge systems to cacuate p f n the fud phase. It aso mpes that n equbrum, p fa = p fb as p f has to be contnuous on the border between A and B. To proceed, we examne the state where system A n fgure 1 s n equbrum wth system B. The phase f of system B has the pressure p f. The equbrum condton that makes the pressure of system A equa to the pressure of system B, s p = p f. (6) By ntroducng ths nto equaton 5, we obtan the ntegra pressure nsde V r at equbrum, ˆp r = p f fr Ωfr + γ V r 4 = pf + 3γfr R. (7)

5 where we used that V r = 4πR 3 /3 and Ω fr = 4πR 2 for a sphere of radus R. From the defnton of H [1] the dfferenta pressure s, cf. equaton 2, where we used that p r = (ˆpr V r ) V r = (pf V r ) fr Ωfr V r + γ V r = p f fr Ωfr + γ V r Ω fr V r = Ωfr R ( V r = pf + 2γfr R, (8) R ) 1 = 2 R. (9) We see that the dfferenta pressure nsde the gran satsfes Young-Lapace s aw for pressure dfferences across curved nterfaces. The fact that we recover ths we-known aw, supports the vadty of the assumptons made. By subtractng equaton 8 from equaton 7, we obtan an nterestng new reaton ˆp r p r = γfr R The expresson reates the ntegra and dfferenta pressure for a spherca phase r of radus R. It s cear that ths pressure dfference s amost equay senstve to the radus of curvature as s the pressure dfference n Young-Lapace s aw. We aso see from ths exampe how the ntegra pressure enters n the descrpton of sma systems. The ntegra pressure s not equa to our norma buk pressure, caed the dfferenta pressure by H, ˆp r p r. Whe two dfferenta pressures satsfy Young-Lapace s aw, the ntegra pressure does not! The ntegra pressure has the property that when averaged over system A usng equaton 5, t s the same as n system B, cf. equaton 6. Ths anayss shows that system A s a possbe, or as we sha see proper, choce of a REV that contans the sod gran, whe system B s a possbe choce of a REV that contans ony fud. An aternatve way to derve equaton 6 s gven n the Appendx. 2.2 Case II. Lattce of spherca grans The above expanaton concerned a snge spherca gran, and was a frst step n the deveopment of a procedure to determne the pressure of a nano-porous medum. To create a more reastc mode, we ntroduce now a attce of spherca grans. The ntegra pressure of a REV contanng many grans s gven by an extenson of equaton 3 (10) ˆpV = ˆp f V f + n ˆp r V r n ˆγ fr Ω fr, (11) where n s the number of grans n the voume. For each gran one may foow the same dervaton for the ntegra and dfferenta pressure as for the snge gran. The equbrum condton s By usng equaton 7, we obtan ˆp = ˆp f = p f. (12) ˆp r = p f + γ fr Ω fr V r = p f + 3γfr, (13) R 5

6 Fgure 2: A sce of the smuaton box n case II. The box has sde engths L x, L y, L z, and propertes are cacuated aong the x-axs n ayers of wdth x. Bue partces are gran r and red partces are fud f. The A s the attce constant of the fcc attce, d 1 and d 2 are the two shortest surface-to-surface dstances. where the ast dentty appes to spherca grans ony. The dfferenta pressure of the grans s gven by a generazaton of equaton 8 p r = (ˆpr V r ) V r = p f + γ fr = (pf V r ) V r Ω fr V r + γ fr Ω fr V r = p f + 2γfr, (14) R where the ast dentty s ony for spherca grans. The dfferenta pressures agan satsfy Young- Lapace s aw. When a grans are dentca spheres and postoned on a fcc attce, a propery chosen ayer coverng haf the unt ce can be a proper choce of the REV. We sha see how ths can be understood n more deta from the moecuar dynamcs smuatons beow. 3 Moecuar dynamcs smuatons Cases I and II were smuated at equbrum, whe Case II was smuated aso away from equbrum. Fgures 4 and 6 ustrate the equbrum smuatons of the two cases. 3.1 Systems The smuaton box was three-dmensona wth sde engths L x, L y, L z, where t was eongated n the x-drecton, L x > L y = L z. Perodc boundary condtons were used n a drectons n the equbrum smuatons. In the non-equbrum smuaton refectng partce boundares [15] were apped to the x-drecton, cf. secton 3.5. Aong the x-axs, the smuaton box was dvded nto n rectanguar cubods (caed ayers) of sze x, L y, L z, where x = L x /n. The voume of each ayer s V = xl y L z. There are two regons A and B n the smuaton box. Regon 6

7 Fgure 3: The rada dstrbuton functon of fud partces around a gran, as shown n fgure 4. Resuts are shown for denstes that vary between ρ = 0.1 and ρ = 0.7. A contans fud (red partces) and gran (bue partces) and regon B contans ony fud, see fgure 2. The regons, B = B 1 + B 2 and A do not have the same sze, but the ayers have the same thckness, x. The compressona energy of the fud n one ayer s, ˆp f V f = p V f. The smuaton was carred out wth LAMMPS [16] n the canonca ensembe usng the Nosé- Hoover thermostat [17], at constant temperature T = 2.0 (n Lennard-Jones unts). The crtca temperature for the Lennard-Jones/spne potenta (LJ/s) s approxmatey Tc 0.9. Fud denstes range from ρ = 0.01 to ρ = Case studes In Case I the snge spherca gran was paced n the center of the box. A perodc mage of the spherca gran s a dstance L x, L y and L z away n the x, y and z-drectons, see fgure 4a. The surface to surface dstance of the spherca grans s d = L α 2R, where R s the radus of the gran, and α = y, z. In Case I, each spherca gran has four nearest neghbours n the perodc attce that s but when we use perodc boundary condtons. We consdered two nearest neghbour dstances; d = 4σ 0 and d = 11σ 0, where σ 0 s the dameter of the fud partces. In case II, the spherca grans were paced n a fcc attce wth attce constant A. The two shortest dstances between the surfaces were characterzed by d 1 = 1 2 ( 2a 4R) and d 2 = a 2R, see fgure 2, where d 1 < d 2. We used d 1 = 4.14σ 0 and d 1 = 11.21σ 0, whch s amost the same as the dstances consdered n Case I. The correspondng other dstances were d 2 = 10σ 0 and d 2 = 20σ 0. Each gran has 12 nearest neghbours at a dstance d 1. The rada rock-fud par correaton functon was computed for varous fud denstes, see fgure 3. The radus of the partce s R = 5.4σ 0. In a cases we computed the voume of the grans V r, the surface area Ωfr and the compressona energy of each ayer,, n the x-drecton. 7

8 3.3 Partce nteracton potenta The partces nteract wth the Lennard-Jones/spne potenta, [ f r < R j ( ) 12 ( ) ] 6 σj R 4ɛ j σj R j r R u j (r) = j j r R j f R j < r < r s,j. (15) a j (r r c,j ) 2 + b j (r r c,j ) 3 f r s < r < r c,j 0 f r > r c,j Each partce type has a hard-core dameter R and a soft-core dameter σ. There were two types of partces, sma partces wth σ ff = σ 0, R ff = 0 and arge partces wth σ rr = 10σ 0, R rr = 9σ 0. The sma partces are the fud (f), and the arge partces are the gran (r). The hard-core and soft-core dameters for fud-gran pars are gven by the Lorentz mxng rue R fr = 1 2 (R ff + R rr ) and σ fr = 1 2 (σ ff + σ rr ). (16) We defne the radus of the gran partces as R (σ ff + σ rr )/2 = 5.5σ 0, whch s the dstance from the gran center where the potenta energy s zero. Ths mpes that fud partces can occupy a poston coser to the gran than ths, especay at hgh pressures, ths s shown n fgure 3. The nteracton strength ɛ j was set to ɛ 0 for a partce-partce pars. The potenta and ts dervatve are contnuous n r = r c,j. The parameters a j, b j and r s,j are determned so that the potenta and the dervatve of the potenta (the force) are contnuous at r = r s,j. 3.4 Pressure computatons The contrbuton of the fud to the grand potenta of ayer s [12] p f V f = 1 m (v v ) 1 N (r j f j ), (17) 3 6 j=1 where p f s the fud ntegra pressure, V f the fud voume, m and v are the mass and veocty of fud partce. The frst two sums are over a fud partces ayer, whe the second sum s over a other partces j. r j r r j s the vector connectng partce and j, and f j = u j / r j s the force between them. The means a nner product of the vectors. The computaton gves ˆp f whch s the contrbuton to the ntegra pressure n ayer from the fud partces, accountng for ther nteracton wth the gran partces. 3.5 The porous medum n a pressure gradent We used the refectng partce boundary method deveoped by L et a. [15] to generate a pressure dfference across the system aong the x-axs. Where partces pass the perodc boundary at x = 0 and x = L x wth probabty (1 α p ) and refected wth probabty α p. A arge α p gves a hgh pressure dfference and a ow α p gves a ow pressure dfference. 4 Resuts and dscusson The resuts of the moecuar dynamcs smuatons are shown n fgures 4-8 (equbrum) and fgures 9, 10 (away from equbrum). The porous medum structure was characterzed by ts 8

9 Fgure 4: a) Iustraton of case I, a snge spherca gran surrounded by a fud phase wth d = 4σ0. b) Voume of gran, V r, c) surface area Ωf r and d) compressona energy pv as a functon of the x-axs of the smuaton box. par correaton functon, cf. fgure 3. The compressona energy was computed accordng to equaton 5 n case I wth a snge spherca gran and case II wth a attce of spherca grans. We computed the compressona energy, p V, n the buk qud (regon B) and n the nanoporous medum (regon A). In the buk qud we computed the pressure drecty from the compressona energy, because p V = pf Vf (not shown). Fgures 4 and 6 show the contrbutons from the sod phase p r Vr and fud-sod surface γf r Ωf r, cf. equaton 5. The gran partces were dentca and the system was n equbrum, so the ntegra pressure n the grans was everywhere the same, p r = p r. Smary, the surface tenson was everywhere the same, γf r = γ f r. In fgures 5, 7 and 8 the gran pressure p r and surface tenson γ f r are potted as a functon of the fud pressure pf. The resuts for Case II were next used n fgures 9 and 10 to determne the pressure gradent across the sequence of REVs n the porous medum. 4.1 Case I. Snge spherca gran. Equbrum The snge sphere case s ustrated n fgure 4a. Fgures 4b and 4c show the varaton n the voume of the porous medum (rock), Vr, and the surface area between the rock and the fud, Ωf r, aong the x-axs of the smuaton box. The two quanttes were determned for a ayers, 9

10 Fgure 5: Ftted gran pressure ˆp r and surface tenson γ fr as a functon of pressure p for a sphere (characterstc ength d = 4σ 0 and d = 11σ 0 )., and these resuts were used n the pots of fgure 4b and 4c. To be representatve, the REV must ncude the sod sphere wth boundares eft and rght of the sphere. In order to obtan p REV V REV we summed p V over a the ayers n the REV. At equbrum, p REV = p, where p s the pressure n the fud n regon B. For the REV we then have pv REV = p f V f + ˆp r V r γ fr Ω fr, (18) REV REV REV where we used that ˆp r = ˆp r and γ fr = γ fr. We know the vaues of a the eements n ths equaton, except ˆp r and γ fr. The vaues of ˆp r and γ fr are next chosen such that equaton 18 hods. Wth these ftted vaues avaabe, we cacuated p V of each ayer from p V = p f V f + ˆp r V r γ fr Ω fr. (19) The contrbutons to the compressona energy n ths equaton for Case I are shown n the bottom fgure 4d. We see the contrbuton from (1) the buk fud p f V f, (2) the buk fud and gran p f V f + ˆp r V r and (3) the tota compressona energy, p V = p f V f + ˆp r V r γ fr Ω fr, whch gves the pressure of the REV when summed and dvded wth the voume of the REV. Fgure 4d shows ceary, that the buk pressure energy gves the argest contrbuton, as one woud expect. It s aso cear that the surface energy s sgnfcant. As the surface to voume rato ncreases, the buk contrbutons may become smaer than the surface contrbuton. In the present case, ths w happen when the radus of the sphere s 2.25σ 0. For our grans wth R = 5.5σ 0 ths does not happen. The pots of ˆp r and γ fr as functons of p n regon B are shown n fgure 5. The vaues for d = 4σ 0 and d = 11σ 0 are gven n the same pots. We see that the pots fa on top of each other. Ths shows that the ntegra pressure and the surface tenson are ndependent of the dstance d n the nterva consdered, and therefore above the ower mt 4σ 0 < d. If confnement effects were essenta, we woud expect that ˆp r and γ fr were functons of the dstance d between the surfaces of the spheres. When the vaue of d decreases beow 4σ 0, devatons may arse, for nstance due to contrbutons from the dsjonng pressure. Such a contrbuton s expected to vary wth the surface area, and ncrease as the dstance between nterfaces become shorter. In pots ke fgure 5, we may see ths as a decrease n the surface tenson. 4.2 Case II. Lattce of spherca grans. Equbrum Consder next the attce of spherca grans, ustrated n fgure 6a. Fgure 6b and 6c gve the varaton n the voume of the porous medum V r and surface area, Ω fr, aong the x-axs. 10

11 Fgure 6: a) Iustraton of case II, a attce spherca gran surrounded by a fud phase wth a = 20σ0. b) Voume of gran, V r, c) surface area Ωf r and d) compressona energy pv as a functon of the x-axs of the smuaton box. The smaest REV s haf a unt ce. 11

12 When the REV n regon A s propery chosen, we know that p REV = p. In equbrum, the pressure of the REV s constant n the buk qud phases, n regons B 1 or B 2, where p s the pressure of the fud n regon B. In order to obtan pv REV n regon A, we sum p V over a the ayers that make up the REV, and obtan pv REV = REV p f V f + ˆp r REV V r γ fr REV Ω fr, (20) To proceed, we fnd frst the vaues of a the eements n ths equaton, except ˆp r and γ fr. The vaues of ˆp r and γ fr are then chosen such that equaton 20 hods. Usng these ftted vaues, we next cacuated ˆp V of each ayer usng p V = p f V f + ˆp r V r γ fr Ω fr (21) The contrbutons to the compressona energy n ths equaton are shown n three stages n fgure 6d: (1) buk fud contrbuton p f V f, (2) buk fud and gran contrbuton pf V f + ˆp r V r and (3) the tota compressona energy, p V = p f V f + ˆp r V r γ fr Ω fr. Fgure 6d shows ceary that the buk contrbuton s argest, as s expected. However, the surface energy s sgnfcant. From fgure 6b t foows that a proper choce of the REV s haf a unt ce, because a REVs are then dentca, (except the REVs at the boundares). The ntegra over p V n these REVs s the same and equa to pv REV. The ayers are smaer than the REV and as a consequence ˆp V w fuctuate, a fuctuaton that s seen n fgure 6d. The vaues for ˆp r and γ fr are shown as functon of p for Case II n fgure 7 for d 1 = 4.14σ 0 and d 1 = 11.21σ 0. We see now a systematc dfference between the vaues of ˆp r and γ fr n the two cases. The ntergra pressure and the surface tenson ncreases as the dstance between the grans decreases. The dfference n one set can be estmated from the other. Say, for a dfference n surface tenson γ fr we obtan for the same fud pressure from equaton 10, a dfference n ntegra pressure of ˆp r = 3 γ fr /R. Ths s neary what we fnd n by comparng the nes n fgure 6, the nes can be predcted from one another usng R = 6.5σ 0 whe the vaue n fgure 3 s R = 5.5σ 0. The dfference must be due to the dsjonng pressure. Its dstrbuton s not sphercay symmetrc, whch may expan the dfference between 6.5σ 0 and 5.5σ 0. The resuts shoud be the same as for Case I for the arger dstance, and ndeed that s found, cf. fgure 8. As the dstance between the gran surfaces ncreases, we expect the dependence of confnement to dsappear, and ths documented by fgure 8 where the two Cases are shown wth dstances d = 11σ 0 and d 1 = 11.21σ 0, respectvey. The curves for the snge gran and attce of grans overap. Our anayss therefore shows that the pressure nsde grans n a fcc attce and the surface tenson, depends n partcuar on the dstances between the surfaces of the spheres, ncudng on ther perodc repcas. A procedure has been deveoped to fnd the pressure of a REV, from nformaton of the (equbrum) vaues of ˆp r and γ fr as a functon of p. It has been documented for partcuar nano-porous medum, but s key to hod for other attces. 4.3 The representatve eementary voume. The knowedge ganed above on the varous pressures at equbrum, can be used to frst construct the REV and next the pressure varaton away from equbrum. Foowng Kjestrup et a. [10], we assume then that there s oca equbrum n the REV. The sze of the REV was estabshed above, as the mnmum sze that cover the compete range of potenta nteractons avaabe n the system. To fnd a REV-property, we need to sampe the whoe space of possbe nteractons. 12

13 Fgure 7: Ftted gran pressure ˆp r and surface tenson γ fr as a functon of pressure p for the attce of spheres (characterstc ength d 1 = 4.14σ 0 and d 1 = 11.21σ 0 ). Fgure 8: Ftted gran pressure ˆp r and surface tenson γ fr as a functon of pressure p for the sphere (characterstc ength d = 11σ 0 ) and a attce of spheres (characterstc ength d 1 = 11.21σ 0 ). 13

14 Fgure 9: a) Iustraton of case II n a pressure gradent. b) Compressona energy pv varaton across the system. Fgure 10: Compressona energy pv varaton across the system smoothed over the representary eementary voume. The thckness of the REV s therefore arger than the ayer thckness used n the smuatons. To show how a REV-property s determned from the ayer-property, we consder agan the compressona energes of each ayer, but now for a system n a pressure gradent. In the anayss we used the fcc attce wth attce parameter a = 20σ0. The voume of the gran, V r, and the surface area, Ωf r, vared of course n the exact same way as n fgure 6b and c. The pressure gradent was generated as expaned n secton 3.5. The pressure dfference between the externa reservors B1 and B2 was arge, gvng a gradent wth order of magntude 1012 bar/m. The fud on the eft hand sde s qud-ke, whe the fud on the rght-hand sde s gas-ke. Fgure 9 ustrates the system n the pressure gradent, where subfgure 9b shows the compressona energy, pv, aong the x-axs. The dp n the pressure cose to x = 0 s caused by the refectng partce boundary, cf. secton 3.5. The refectng partce boundary ntroduces a surface between the hgh fud on the eft sde and the ow pressure on the rght sde. The smaest REV as obtaned n the anayss at equbrum s ndcated n the fgure. Wthn a REV, the varabes fuctuate, as ceary shown n the equbrum resuts, fgure 6. In order to compute a REV varabe, we foow the procedure descrbed by Kjestrup et a [10] and choose a ayer as a reference pont. We then compute the movng average usng fve ayers, 14

15 two to the eft, two to the rght and the centra ayer. Movng one ayer down the gradent, we repeat the procedure, and n ths manner we obtan the property varaton on the REV scae. The resuts of the smuaton gave, for each ndvdua ayer, p V, as potted n fgure 9b. The profe created by the REV-centers s shown n fgure 10. We see a smooth near profe (centra curve) as one woud expect from the boundary condtons that are mposed on the system. Some traces of oscaton are st eft n the separate contrbutons to the tota compressona energy. The tota energy s constant n regons A, as expected. 4.4 Case II. Lattce of spherca grans. Non-equbrum The resuts shown n fgure 10 mean that the ntegra as we as the dfferenta pressures vary through the system. The reaton between p and ˆp r, γ fr n the REV s the same as n fgure 8 at oca equbrum. We can therefore cacuate ˆp r and γ fr for each ayer from ths fgure for each REV. The compressona energy p V, as determned from p V = p f V f + ˆp r V r γ fr Ω fr (22) gves the tota compressona energy p V, and fnay the REV compressona energy. The ast property shows an essentay near varaton from the vaue n one buk fud to the vaue n the other buk fud, cf. fgure 10. The fuctuatons around the mean vaue were dscussed above. The REV shoud be arge enough to emnate these. We have found that a nano-porous medum has an ncreased number of varabes. We need consder the pressures n the fud and the sod phases, as we as the surface tenson between the fud and the sod. When one reduces the sze of a thermodynamc system to the nano-meter sze, the pressures and the surface tensons become dependent on the sze of the system. An mportant observaton s then that there are two reevant pressures rather than one. H [1] caed them the ntegra and the dfferenta pressure, respectvey. It s maybe surprsng that the smpe vra expresson works so we for a pressure cacuatons n a fud, but we have found that t can be used. Ths means that we w next be abe to study transport processes, where the externa pressure dfference s a drvng force. The method, to compute the mechanca force ntrnsc to the porous medum may open nterestng new possbtes to study the effects that are characterstc for porous meda. In a macro-scae descrpton the so-caed representatve eementary voume (REV) s essenta. Ths makes t possbe to obtan thermodynamc varabes on ths scae. We have here dscussed how the fact, that the macro-scae pressure s constant n equbrum, makes t possbe to obtan the ntegra pressure n the sod, as we as the surface tenson of the REV. An observaton whch confrms the soundness of the procedure s that we recover that the resutng dfferenta pressure n the sod satsfes Young-Lapace s aw. The exstence of a REV for systems on the nano-scae suppements the REV that can be defned for pores of mcrometer dmenson [10]. There s no confct between the eves of descrpton as they merge n the thermodynamc mt. The REV, as defned n the present work, may aow us to deveop a non-equbrum thermodynamc theory for the nano-scae. 5 Concusons The foowng concusons can be drawn from the above studes We have obtaned the frst support for a new way to compute the pressure n a nano-porous medum. The ntegra pressure s defned by the grand potenta. The defnton appes to 15

16 the thermodynamc mt, as we as to systems whch are sma accordng to the defnton of H [1]. It foows that nano-porous meda need two pressures n ther descrpton, the ntegra and the dfferenta pressure. Ths s new knowedge n the context of nano-porous meda For a spherca rock partce of radus R, we derve a reaton between the ntegra and the dfferenta pressure n terms of the surface tenson, ˆp r p r = γ/r. Ther dfference s non-neggbe n the cases where Young-Lapace s aw appes. To ustrate the cacuaton, we have constructed a system wth a snge fud. The rock pressure and the surface tenson are constant throughout the porous medum at equbrum. The assumptons were confrmed for a porosty change from φ = 0.74 to 0.92, for a REV wth mnmum sze of haf a unt ce. From the assumpton of oca equbrum, we can fnd the pressure nterna to a REV of the porous medum, under non-equbrum condtons, and a contnuous varaton n the pressure on a macro-scae. The procedure to use a vrtua (sster) path to fnd the varabe n queston, may appy aso to other varabes than the pressure [10]. To obtan these concusons, we have used moecuar dynamcs smuatons of a snge spherca gran n a pore and then for face-centered attce of spherca grans n a pore. Ths too s rrepaceabe n ts abty to test assumptons made n the theory. The smuatons were used here to compute the ntegra rock pressure and the surface tenson, as we as the pressure of the representatve voume, and through ths to deveop a procedure for porous meda pressure cacuatons. Ony one fud has been studed here. The stuaton s expected to be much more compcated wth two-phase fow and an amorphous medum. Nevertheess, we beeve that ths frst step has gven usefu nformaton for the work to foow. We sha contnue to use the grand potenta for the more compcated cases, n work towards a non-equbrum thermodynamc theory for the nano-scae. 6 Acknowedgement The cacuaton power was granted by The Norwegan Metacenter of Computatona Scence (NOTUR). Thanks to the Research Counc of Norway through ts Centres of Exceence fundng scheme, project number , PoreLab. Appendx Ths appendx descrbes an aternatve way to derve equaton 6 n secton 2.1 and an ntegra form of the Young-Lapace aw. Insertng equaton 3 nto equaton 2 gves p(v ) = (ˆp f V ) f + (ˆpr V r ) (ˆγ fr Ω fr) (23) V V V 16

17 From H [1] we can express the dfferenta pressures and the nterfaca tenson as p f = (ˆp f V ) f V f p r = (ˆp r V ) f V r γ fr = (ˆγ fr Ω fr) Ω fr The porous medum s characterzed by a porosty, (24a) (24b) (24c) φ = V f V If the change n voume mped by the dfferentaton n equaton 23 s performed at constant porosty (proportona sweng of the two phases) 2, we get and V = φ V f = (1 φ) ) p = φp f + (1 φ) (p r 2γfr R (25) V r = (1 φ) 2 R Ω fr (26) In the mt φ 1 p p f as one shoud expect. However, Equaton 27 nvoves ony dfferenta pressures. We can then appy Young-Lapace s aw, p r = p f + 2γ fr /R, to the ast term and fnd that p = φp f + (1 φ)p f = p f, (28) whch s equaton 6. Usng the concept of porosty, we can rearrange equaton 3: (27) ˆp = φˆp f + (1 φ)ˆp r (1 φ) 3ˆγfr R Ths suggests an aternatve to Young-Lapace s aw on the ntegra eve: ˆp r = ˆp f + 3ˆγfr R (29) (30) References [1] Terre L. H, Thermodynamcs of Sma Systems, New York: Dover, [2] W. G. Gray and C. T. Mer, Thermodynamcay constraned averagng theory approach for modeng fow and transport phenomena n porous medum systems: 8. nterface and common curve dynamcs, Advances n Water Resources, vo. 33, no. 12, pp , [3] L. S. Bennethum and T. Wensten, Three pressures n porous meda, Transport n Porous Meda, vo. 54, no. 1, pp. 1 34, Ths may be the ony way to do the dfferentaton consstent wth Euer homogenety 17

18 [4] J. Magda, M. Trre, and H. Davs, Moecuar dynamcs of narrow, qud-fed pores, The Journa of Chemca Physcs, vo. 83, no. 4, pp , [5] B. Todd, D. J. Evans, and P. J. Davs, Pressure tensor for nhomogeneous fuds, Physca Revew E, vo. 52, no. 2, p. 1627, [6] T. Ikeshoj, B. Hafskjod, and H. Furuhot, Moecuar-eve cacuaton scheme for pressure n nhomogeneous systems of fat and spherca ayers, Moecuar Smuaton, vo. 29, no. 2, pp , [7] B. Hafskjod and T. Ikeshoj, Mcroscopc pressure tensor for hard-sphere fuds, Physca Revew E - Statstca, Nonnear, and Soft Matter Physcs, vo. 66, no. 1, pp. 1 4, [8] S. M. Hassanzadeh and W. G. Gray, Mechancs and thermodynamcs of mutphase fow n porous meda ncudng nterphase boundares, Advances n water resources, vo. 13, no. 4, pp , [9] W. G. Gray and S. M. Hassanzadeh, Macroscae contnuum mechancs for mutphase porous-meda fow ncudng phases, nterfaces, Advances n Water Resources, vo. 21, pp , [10] S. Kjestrup, D. Bedeaux, A. Hansen, B. Hafskjod, and O. Gateand, Non-sotherma transport of mut-phase fuds n porous meda. the entropy producton, Fronters n Physcs, vo. 6, p. 126, [11] S. Kjestrup, D. Bedeaux, A. Hansen, B. Hafskjod, and O. Gateand, Non-sotherma transport of mmscbe fuds n porous meda. II. Consttutve equatons, Accepted, Fronters n Physcs. [12] J. H. Irvng and J. G. Krkwood, The Statstca Mechanca Theory of Transport Processes. IV. The Equatons of Hydrodynamcs, The Journa of Chemca Physcs, vo. 18, no. 6, pp , [13] R. C. Toman, The effect of dropet sze on surface tenson, The Journa of Chemca Physcs, vo. 17, no. 3, pp , [14] W. Hefrch, Eastc propertes of pd bayers: theory and possbe experments, Zetschrft für Naturforschung C, vo. 28, no , pp , [15] J. L, D. Lao, and S. Yp, Coupng contnuum to moecuar-dynamcs smuaton: Refectng partce method and the fed estmator, Physca Revew E, vo. 57, no. 6, pp , [16] S. Pmpton, Fast Parae Agorthms for Short Range Moecuar Dynamcs, Journa of Computatona Physcs, vo. 117, no. June 1994, pp. 1 19, [17] W. G. Hoover and B. L. Hoan, Knetc moments method for the canonca ensembe dstrbuton, Physcs Letters-Secton A, vo. 211, no. 5, pp ,