Lattice Boltzmann simulations of drops colliding with solid surfaces

Video images of drops colliding with solid surfaces shown by Rioboo et al. (2002) reveal that, for large drop velocities, the drops flatten and form a ring structure before receding and, in some cases, rebounding from the surface. They described the sequence of events in terms of four distinct regimes. During the initial kinematic phase, the dimensionless wetting radius of the drop follows a universal form if the drop Weber and Reynolds numbers are sufficiently large. In the second phase, the drop becomes highly flattened and the values of the Weber and Reynolds numbers influence the time evolution of the dimensionless wetting radius and its maximum value. This is followed by a third phase in which the wetting radius begins to decrease with time and the wettability of the surface influences the dynamics. This paper presents simulation results for the early stages of drop impact and spreading on a partially wetting solid surface. The simulations were performed with a modified version of the lattice Boltzmann method (LBM) developed by Inamuro et al. (2004) for a liquid-gas density ratio of 1000. The Inamuro et al. version of the LBM was modified by incorporating rigid, no-slip boundary conditions and incorporating a boundary condition on the normal derivative of the order parameter to impose the desired equilibrium contact angle.     

Monte Carlo Simulation of Hard Hyperspheres in Six,Seven and Eight Dimensions for Low to Moderate Densities

The pair correlation function of hard hyperspheres in six, seven and eight dimensions is obtained from Monte Carlo simulations. The value of the pair correlation function at contact is compared with the results from molecular dynamics calculations and a variety of theoretical approaches. Remarkably good agreement is found with the simple, closed-form equations of Y. Song, E. A. Mason and R. M. Stratt, J. Phys. Chem., 93:6916–6919 (1989). The Monte Carlo results for the equation of state are compared with the theoretical expressions of M. Baus and J. L. Colot, Phys. Rev. A, 36:3912 (1987), M. Luban and J. P. J. Michels, Phys. Rev A, 41:6796 (1990), and high order virial expansions. In addition, in seven dimensions, comparisons are made with the exact PY solution provided by M. Robles, M. L. de Haro and A. Santos, J. Chem. Phys., 120:9113 (2004). Very good agreement was observed between theory and computer simulation in all dimensions.     

Angular Distributions for

[1]G.T.Bodwin,E.Braaten,and G.P.Lepage,Phys.Rev.D 51 (1995) 1125;[Erratum-ibid.D 55 (1997) 5853][arXiv:hep-ph/9407339]; J.Boltz,P.Kroll,and G.A.Schulre,Phys.Lett.B 392 (1997) 198; J.Boltz,P.Kroll,and G.A.Schulre,Phys.J.C 2 (1998) 705. [2]S.M.Wong,Nucl.Phys.A 674 (2000) 185; S.M.Wong,Eur.Phys.J.C 14 (2000) 643. [3]J.Z.Bai,Y.Ban,J.G.Bian,et al.,Phys.Rev.D 67 (2003)112001. [4]M.Jacob and G.C.Wick,Ann.Phys.7 (1959) 404. [5]S.U.Chung,Phys.Rev.D 48 (1993) 1225; S.U.Chung,Phys.Rev.D 57 (1998) 431; B.S.Zou and D.V.Bugg,Eur.Phys.J.A 16 (2003) 537. [6]Particle Data Group,Phys.Lett.B 592 (2004) pp.924-966. [7]M.A.Doncheski,et al.,Phys.Rev.D 42 (1990) 2293; E.Eichten,et al.,Phys.Rev.D 21 (1980) 203; K.J.Sebastian,Phys.Rev.D 26 (1982) 2295; G.Hardekopf and J.Sucher,Phys.Rev.D 25 (1982) 2938; R.McClary and N.Byers,Phys.Rev.D 28 (1983) 1692; P.Moxhay and J.L.Rosner,Phys.Rev.D 28 (1983) 1132. [8]B.S.Zou and F.Hussain,Phys.Rev.C 67 (2003) 015204.     

Universal Scaling Law for Atomic Diffusion and Viscosity in Liquid Metals

The recently proposed scaling law relating the diffusion coefficient and the excess entropy of liquid [Samanta A et al. 2004 Phys. Reu. Lett. 92 145901; Dzugutov M 1996 Nature 381 137], and a quasi-universal relationship between the transport coefficients and excess entropy of dense fluids [Rosenfeld Y 1977 Phys. Rev. A 15 2545],are tested for diverse liquid metals using molecular dynamics simulations. Interatomic potentials derived from the glue potential and second-moment approximation of tight-binding scheme are used to study liquid metals.Our simulation results give sound support to the above-mentioned universal scaling laws. Following Dzugutov,we have also reached a new universal scaling relationship between the viscosity coefficient and excess entropy.The simulation results suggest that the reduced transport coefficients can be expressed approximately in terms of the corresponding packing density.     

Sigma Terms and Strangeness Contents of Baryon Octet in Modified Chiral Perturbation Theory

[1]J.Gasser,H.Leutwyler,and M.E.Sainio,Phys.Lett.B 253 (1991) 252. [2]John Ellis,Eur.Phys.J.A 24S2 (2005) 3,[arXive:hepph/0411369]. [3]T.Inoue,V.E.Lyubovitskij,Th.Gutsche,and Amand Faessler,Phys.Rev.C 69 (2004) 035207,[arXive:hepph/0311275]. [4]M.M.Pavan,I.I.Strakovsky,R.L.Workman,and R.A.Arndt,PiN Newslett.16 (2002) 110,[arXive:hepph/0111066]. [5]V.E.Lyubovitskij,Th.Gutsche,Amand Faessler,and E.G.Drukarev,Phys.Rev.D 63 (2001) 054026,[arXive:hep-ph/0009341]. [6]S.D.Bass,Phys.Lett.B 329 (1994) 358,[arXive:hepph/9404294]. [7]Marc Knecht,PiN Newslett.15 (1999) 108,[arXive:hepph/9912443]. [8]P.Schweitzer,Phys.Rev.D 69 (2004) 034003. [9]B.C.Lehnhart,J.Gegelia,and S.Scherer,J.Phys.G 31(2005) 89,[arXive:hep-ph/0412092]. [10]P.J.Ellis and K.Torikoshi,Phys.Rev.C 61 (1999)015205. [11]Gerald E.Hite,William B.Kaufmann,and Richard J.Jacob,Phys.Rev.C 71 (2005) 065201. [12]S.Weinberg,Physica A 96 (1979) 327. [13]J.Gasser and H.Leutwyler,Nucl.Phys.B 250 (1985)465. [14]J.Gasser,M.E.Sainio,and A.Svarc,Nucl.Phys.B 307(1988) 779. [15]P.Papazoglou,D.Zschiesche,S.Schramm,J.SchaffnerBielich,H.St(o)cker,and W.Greiner,Phys.Rev.C 59(1999) 411. [16]T.Fuchs and J.Gegelia,Phys.Rev.D 68 (2003) 056005.     

Re-study of Nucleon Pole Contribution in J/

[1]V.D.Burkert,Phys.Lett.B 72 (1997) 109. [2]S.Capstick and W.Roberts,Prog.Part.Nucl.Phys.45 (2000) S241,and references therein. [3]B.S.Zou,Nucl.Phys.A 675 (2000) 167c; B.S.Zou,Nucl.Phys.A 684 (2001) 330; BES Collaboration (J.Z.Bai,et al.) Phys.Lett.B 510 (2001) 75; BES Collaboration (M.Ablikim,et al.),hep-ex/0405030. [4]R.Sinha and Susumu Okubo,Phys.Rev.D 30 (1984)2333. [5]W.H.Liang,P.N.Shen,B.S.Zou,and A.Faessler,Euro.Phys.J A 21 (2004) 487. [6]Particle Data Group,Euro.Phys.J.C 15 (2000) 1. [7]K.Tsushima,A.Sibrtsev,and A.W.Thomas,Phys.Lett.B 390 (1997) 29. [8]J.Kogut,Rev.Mod.Phys.51 (1979) 659; Rev.Mod.Phys.55 (1983) 775. [9]Q.Haider and L.C.Liu,J.Phys.G 22 (1996) 1187; L.C.Liu and W.X.Ma,J.Phys.G 26 (2000) L59. [10]V.G.J.Stoks,R.A.M.Klomp,C.P.F.Terheggen,and J.J.de Swart,Phys.Rev.C 49 (1994) 2950. [11]H.Haberzettl,C.Bennhold,T.Mart,and T.Feuster,Phys.Rev.C 58 (1998) R40. [12]Y.Oh,A.I.Titov,and T.-S.H.Lee,Phys.Rev.C 63(2001) 25201.     

Top Quark Flavor Changing Decay t

[1]R. Casalbuoani, A. Deandrea, and M. Oertel, JHEP 032(2004) 0402. [2]G. Hooft, In Search of the Ultimate Building Blocks, Cambridge University Press, Cambridge (1997). [3]J. Belazey, Searches for New Physics at Hadron Coliders,Northern Illinois University (2005). [4]N. Arkani-hamed, A.G. Cohen, and H. Georgi, Phys. Lett.B 513 (2001) 232 [hep-ph/0105239]. [5]I. Low, W. Skiba, and D. Smith, Phys. Rev. D 66 (2002)072001 [hep-ph/0207243]. [6]N. Arkani-hamed, A.G. Cohen, E. Katz, and A.E. Nelson,JHEP 0207 (2002) 304 [hep-ph/0206021]. [7]N. Arkani-hamed, A.G. Cohen, E. Katz, A.E. Nelson, T.Gregoire, and J. G. Wacker, JHEP 0208 (2002) 021 [hepph/0206020]. [8]T. Gregoire and J.G. Wacker, JHEP 0208 (2002) 019[hep-ph/0206023]. [9]For a recent review, see e.g., M. Schmaltz, Nucl. Phys. B (Proc. Suppl.) 117 (2003) 40. [10]N. Arkani-hamed, A.G. Cohen, T. Gregoire, and J.G.Jacker, JHEP 0208 (2002) 020 [hep-ph/0202089]. [11]or a recent review, see e.g., M. Schmaltz, Nucl. Phys.Proc. Suppl. 117 (2003) 40 [hep-ph/0210415]. [12]E. Katz, J. Lee, A.E. Nelson, and D.G. Walker, hepph/0312287. [13]M. Beneke, I. Efthymiopoulos, M.L. Mangano, et al., hepph/0003033. [14]D.O. Carlson and C.-P. Yuan, hep-ph/9211289. [15]R. Frey, D. Gerdes, and J. Jaros, hep-ph/9704243. [16]G. Eilam, J.L. Hewett, and A. Soni, Phys. Rev. D 44(1991) 1473; W.S. Hou, Phys. Lett. B 296 (1992) 179; K.Agashe and M. Graesser, Phys. Rev. D 54 (1996) 4445;M. Hosch, K. Whisnant, and B.L. Young, Phys. Rev. D56 (1997) 5725. [17]C.S. Li, R.J. Oakes, and J.M. Yang, Phys. Rev. D 49(1994) 293, Erratum-ibid. D 56 (1997) 3156; G. Couture,C. Hamzaoui, and H. Koenig, Phys. Rev. D 52 (1995)1713; G. Couture, M. Frank, and H. Koenig, Phys. Rev.D 56 (1997) 4213; G.M. de Divitiis, et al., Nucl. Phys. B 504 (1997) 45. [18]B. Mele, S. Petrarca, and A. Soddu, Phys. Lett. B 435(1998) 401. [19]B. Mele, hep-ph/0003064. [20]J.M. Yang and C.S. Li, Phys. Rev. D 49 (1994) 3412,Erratum, ibid. D 51 (1995) 3974; J.G. Inglada, hepph/9906517. [21]L.R. Xing, W.G. Ma, R.Y. Zhang, Y.B. Sun, and H.S.Hou, Commun. Theor. Phys. (Beijing, China) 41 (2004)241. [22]L.R. Xing, W.G. Ma, R.Y. Zhang, Y.B. Sun, and H.S.Hou, Commun. Theor. Phys. (Beijing, China) 40 (2003)171. [23]T. Han, H.E. Logan, B. McElrath, and L.T. Wang, Phys.Rev. D 67 (2003) 095004. [24]I. Low, W. Skiba, and D. Smith, Phys. Rev. D 66 (2002)072001. [25]T. Han, H.E. Logan, B. McElrath, and L.T. Wang, hepph/0302188. [26]A.J. Buras, A. Poschenrieder, and S. Uhlig, hepph/0410309. [27]S. Eidelman, et al., Phys. Lett. B 592 (2004) 1. [28]F. Legerlehner, DESY 01-029, hep-ph/0105283.     

Mixing,ergodicity and slow relaxation phenomena

Investigations on diffusion in systems with memory [I.V.L. Costa, R. Morgado, M.V.B.T. Lima, F.A. Oliveira, Europhys. Lett. 63 (2003) 173] have established a hierarchical connection between mixing, ergodicity, and the fluctuation–dissipation theorem (FDT). This hierarchy means that ergodicity is a necessary condition for the validity of the FDT, and mixing is a necessary condition for ergodicity. In this work, we compare those results with recent investigations using the Lee recurrence relations method [M.H. Lee, Phys. Rev. B 26 (1982) 2547; M.H. Lee, Phys. Rev. Lett. 87 (2001) 250601; M.H. Lee, J. Phys. A: Math. Gen. 39 (2006) 4651]. Lee shows that ergodicity is violated in the dynamics of the electron gas [M.H. Lee, J. Phys. A: Math. Gen. 39 (2006) 4651]. This reinforces both works and implies that the results of [I.V.L. Costa, R. Morgado, M.V.B.T. Lima, F.A. Oliveira, Europhys. Lett. 63 (2003) 173] are more general than the framework in which they were obtained. Some applications to slow relaxation phenomena are discussed.     

Lattice Boltzmann simulations of micron-scale drop impact on dry surfaces

A lattice Boltzmann equation (LBE) method for incompressible binary fluids is proposed to model the contact line dynamics on partially wetting surfaces. Intermolecular interactions between a wall and fluids are represented by the inclusion of the cubic wall energy in the expression of the total free energy. The proposed boundary conditions eliminate the parasitic currents in the vicinity of the contact line. The LBE method is applied to micron-scale drop impact on dry surfaces, which is commonly encountered in drop-on-demand inkjet applications. For comparison with the existing experimental results [H. Dong, W.W. Carr, D.G. Bucknall, J.F. Morris, Temporally-resolved inkjet drop impaction on surfaces, AIChE J. 53 (2007) 2606–2617], computations are performed in the range of equilibrium contact angles from 31° to 107° for a fixed density ratio of 842, viscosity ratio of 51, Ohnesorge number (Oh) of 0.015, and two Weber numbers (We) of 13 and 103.     

Dynamics of Two-Photon Lasers with

[1]C.O.Weiss and R.Vilaseca,Dynamics of Lasers,VCH,Weinheim (1991); Instabilities and Chaos in Quantum Optics,eds.F.T.Arecchi and R.G.Harrison,Springer-Verlag,Berlin (1987). [2]H.Haken,Phys.Lett.A 53 (1975) 77. [3]Ju Rui,Huang Hong-Bin,Yang Peng,Xie Xia,and Zhao Huan,Commun.Theor.Phys.(Beijing,China) 44 (2005) 65; Ju Rui,Zhang Ya-Jun,Huang Hong-Bin,and Zhao Huan,Acta Phys.Sin.53 (2004) 2191 (in Chinese). [4]C.Z.Ning and H.Haken,Z.Phys.B 77 (1989) 247; B 77 (1989) 157; B 77 (1989) 163; J.Zakrenwski and M.Lewenstein,Phys.Rev.A 45 (1992) 2057. [5]G.J.deValearcel,E.Roldan,and R.Vilaseca,Phys.Rev.A 45 (1992) R2674; Phys.Rev.A 49 (1994) 1243. [6]X.Xie,H.B.Huang,F.Qian,Y.J.Zhang,P.Yang,and G.X.Qi,Commun.Theor.Phys.(Beijing,China) 46 (2006) 1042. [7]X.L.Deng,H.Q.Ma,B.D.Chen,and H.B.Huang,Phys.Lett.A 290 (2001) 77. [8]C.Benkert,and M.O.Scully,Phys.Rev.A 42 (1990) 2817. [9]M.O.Scully and M.S.Zubairy,Quantum Optics,Cambridge University Press,Cambridge (1997).     

Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer-Meshkov Instability

The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time (MRT) Lattice Boltzmann (LB) model proposed in [Europhys. Lett. 90 (2010) 54003]. We discuss the reason why the Gram-Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 (2004) 035701(R)] is only valid in subsonic flows.The von Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer-Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.     

The debate on the dependence of apparent contact angles on drop contact area or three-phase contact line: A review

A sessile drop is an isolated drop which has been deposited on a solid substrate where the wetted area is limited by the three-phase contact line and characterized by contact angle, contact radius and drop height. Although, wetting has been studied using contact angles of drops on solids for more than 200 years, the question remains unanswered: Is wetting of a rough and chemically heterogeneous surface controlled by the interactions within the solid/liquid contact area beneath the droplet or only at the three-phase contact line? After the publications of Pease in 1945, Extrand in 1997, 2003 and Gao and McCarthy in 2007 and 2009, it was proposed that advancing, receding contact angles, and contact angle hysteresis of rough and chemically heterogeneous surfaces are determined by interactions of the liquid and the solid at the three-phase contact line alone and the interfacial area within the contact perimeter is irrelevant. As a consequence of this statement, the well-known Wenzel (1934) and Cassie (1945) equations which were derived using the contact area approach are proposed to be invalid and should be abandoned. A hot debate started in the field of surface science after 2007, between the three-phase contact line and interfacial contact area approach defenders. This paper presents a review of the published articles on contact angles and summarizes the views of the both sides. After presenting a brief history of the contact angles and their measurement methods, we discussed the basic contact angle theory and applications of contact angles on the characterization of flat, rough and micropatterned superhydrophobic surfaces. The weak and strong sides of both three-phase contact line and contact area approaches were discussed in detail and some practical conclusions were drawn.     

Scale-separation scheme for simulating superfluid turbulence: Kelvin-wave cascade

A Kolmogorov-type cascade of Kelvin waves-the distortion waves on vortex lines-plays a key part in the relaxation of superfluid turbulence at low temperatures. We propose an efficient numeric scheme for simulating the Kelvin-wave cascade on a single vortex line. This idea is likely to be generalizable for a full-scale simulation of different regimes of superfluid turbulence. With the new scheme, we are able to unambiguously resolve the cascade spectrum exponent, and thus to settle the controversy between recent simulations of Vinen, Tsubota, and Mitani [Phys. Rev. Lett. 91, 135301 (2003)]] and recently developed analytic theory [Phys. Rev. Lett. 92, 035301 (2004)]].     

Quasi-equilibrium lattice Boltzmann method

A general lattice Boltzmann method for simulation of fluids with tailored transport coefficients is presented. It is based on the recently introduced quasi-equilibrium kinetic models, and a general lattice Boltzmann implementation is developed. Lattice Boltzmann models for isothermal binary mixtures with a given Schmidt number, and for a weakly compressible flow with a given Prandtl number are derived and validated.     

基于六粒子纠缠态和Bell态测量的量子信息分离

通过介绍六粒子纠缠态的新应用研究,提出了一个二粒子任意态的信息分离方案.在这个方案中,发送者Alice、控制者Charlie和接受者Bob共享一个六粒子纠缠态,发送者先执行两次Bell基测量|然后控制者执行一次Bell基测量|最后接受者根据发送者和控制者的测量结果,对自己拥有的粒子做适当的幺正变换,从而能够重建要发送的二粒子任意态.这个信息分离方案是决定性的,即成功概率为100%.与使用相同的量子信道进行二粒子任意态的信息分离方案相比,本文提出的方案只需要进行Bell基测量而不需要执行多粒子的联合测量,从而使得这个方案更简单、更容易,并且在目前的实验室技术条件下是能够实现的.     

Pair-Breaking Critical Current Density of Two-Band Superconductor MgB2

[1]J. Nagamatsu, N. Nakagava, T. Muranaka, Y. Zenitani,and J. Akimitsu, Nature 410 (2001) 63. [2]C. Buzea and T. Yamashita, Supercond. Sci. Techn. 14(2001) R115. [3]S. Budko, G. Lapertot, C. Petrovic, C.E. Gunningham, N.Anderson, and P.C. Canfield, Phys. Rev. Lett. 86 (2001)1877. [4]H. Kotegawa, K. Ishida, Y. Kitaoka, T. Muranaka, and J. Akimitsu, Phys. Rev. Lett. 87 (2001) 127001. [5]J. Kortus, I.I. Mazin, K.D. Belashchenko, V.P. Antropov,and L.L. Boyer, Phys. Rev. Lett. 87 (2001) 4656. [6]A. Liu, I.I. Mazin, and J. Kortus, Phys. Rev. Lett. 87(2001) 087005. [7]X.K. Chen, M.J. Konstantinovich, J.C. Irwin, D.D.Lawrie, and J.P. Frank, Phys. Rev. Lett. 87 (2001)157002. [8]H. Giublio, D. Roditchev, W. Sacks, R. Lamy, D.X.Thanh, J. Kleins, S. Miraglia, D. Fruchart, J. Markus,and P. Monod, Phys. Rev. Lett. 87 (2001) 177008. [9]F. Bouquet, R.A. Fisher, N.E. Phillips, D.G. Hinks, and J.D. Jorgensen, Phys. Rev. Lett. 87 (2001) 04700. [10]S.V. Shulga, S.-L. Drechsler, H. Echrig, H. Rosner, and W. Pickett, Cond-mat/0103154 (2001). [11]A.A. Golubov, J. Kortus, O.V. Dolgov, O. Jepsen, Y.Kong, O.K. Andersen, B.J. Gibson, K. Ahn, and R.K.Kremer, J. Phys. Condens. Matter 14 (2002) 1353. [12]H. Doh, M. Sigrist, B.K. Chao, and Sung-Ik Lee, Phys.Rev. Lett. 85 (1999) 5350. [13]I.N. Askerzade, N. Guclu, and A. Gencer, Supercond. Sci.Techn. 15 (2002) L13. [14]I.N. Askerzade, N. Guclu, A. Gencer, and A. Kiliq, Supercond. Sci. Techn. 15 (2002) L17. [15]I.N. Askerzade and A. Gencer, J. Phys. Soc. Jpn. 71(2002) 1637. [16]I.N. Askerzade, Physica C 397 (2003) 99. [17]V.V. Anshukova, B.M. Bulychev, A.I. Golovashkin, L.I.Ivanova, A.A. Minakov, and A.P. Rusakov, Phys. Solid State 45 (2003) 1207. [18]A.A. Abrikosov, Fundamentals of the Theory of Metals,North-Holland, Amsterdam (1988). [19]M.N. Kunchur, S.I. Lee, and W.N. Kang, Phys. Rev. B 68 (2003) 064516.     

Field-induced vortices in weakly anisotropic ferroelectrics

In microscale and nanoscale ferroelectric samples, the formation and growth of domains are the usual stages of the polarization switching mechanism. By assuming weak polarization anisotropy and by solving the Ginzburg–Landau–Khalatnikov equation we have explored an alternative mechanism which consists in ferroelectric switching induced by vortex formation. We have studied the polarization dynamics inside a ferroelectric circular capacitor where switching leads to the formation of a metastable vortex state with a rotational motion of the polarization. Our results are consistent with recent first-principle simulations [I.I. Naumov, H.X. Fu, Phys. Rev. Lett. 98, 077603 (2007)] and with experiments on PbZr_0.2Ti_0.8O₃ [A. Gruverman, D. Wu, H.J. Fan, I. Vrejoiu, M. Alexe, R.J. Harrison, J.F. Scott, J. Phys. Condens. Matter 20 342201(2008)] and demonstrate that vortex-induced polarization switching can be an effective mechanism for circular nanocapacitors.     

Shape Coexistence for 179Hg in Relativistic Mean-Field Theory

[1]J.H. Hamilton,A. VRamayya, W.T. Pinkston, et al.,Phys. Rev. Lett. 32 (1974) 239. [2]R. Julin, K. Helariutta, and M. Muikku, J. Phys. G 27(2001) R109. [3]J.H. Hamilton, Nukleonika 24 (1979) 561. [4]W.C. Ma, et al., Phys. Lett. B 139 (1984) 276. [5]R. Bengtsson, et al., Phys. Lett. B 183 (1987) 1. [6]S. Yoshida and N. Takigawa, Phys. Rev. C 55 (1996)1255. [7]T. Niksic, D. Vretenar, P. Ring, et al., Phys. Rev. C 65(2002) 054320. [8]F.G. Condev, M.P. Carpenter, R.V.F. Janssens, et al.,Phys. Lett. B 528 (2002) 221. [9]D.G. Jenkins, A.N. Andreyev, R.D. Page, et al., Phys.Rev. C 66 (2002) 011301(R). [10]B.D. Serot and J.D. Walecka, Adv. Nuc]. Phys. 16 (1986)1. [11]P. Ring, Prog. Part. Nucl. Phys. 37 (1996) 193. [12]J. Meng and P. Ring, Phys. Rev. Lett. 77 (1996) 3963. [13]J. Meng and P. Ring, Phys. Rev. Lett. 80 (1998) 460. [14]S.K. Patra, S. Yoshida, N. Takigawa, and C.R. Praharaj,Phys. Rev. C 50 (1994) 1924. [15]S. Yoshida, S.K. Patra, N. Takigawa, and C.R. Praharaj,Phys. Rev. C 50 (1994) 1938. [16]G.A. Lalazissis and P. Ring, Phys. Lett. B 427 (1998)225. [17]Jun-Qing Li, Zhong-Yu Ma, Bao-Qiu Chen, and Yong Zhou, Phys. Rev. C 65 (2002) 064305. [18]G. Audi and A.H. Wapstra, Nucl. Phys. A 565 (1993) 1. [19]G. Audi and A.H. Wapstra, Nucl. Phys. A 595 (1995)409. [20]G. Audi and A.H. Wapstra, Nucl. Phys. A 624 (1997) 1. [21]P. MOller and J.R. Nix, Atom. Data and Nucl. Data Table 59 (1995) 307.     

润湿梯度作用下液滴在倾斜表面向上运动的格子Boltzmann模拟

采用基于Shan-Chen伪势模型的格子Boltzmann方法,对液滴在存在润湿梯度的倾斜表面上克服重力、自下而上运动的过程进行模拟。探究润湿梯度、液滴尺寸、Bond数以及表面倾斜角度对液滴运动的影响。计算结果表明：液滴在运动过程中,内部会出现沿斜面向上的速度矢量,润湿梯度越大,液滴运动速度越快,润湿长度也越长,且动态接触角减小速率越快。液滴尺寸和Bond数对液滴运动的影响较小,但存在临界Bond数,超过该临界Bond数时,液滴将沿梯度润湿表面向下运动。表面倾角对液滴运动有显著影响,倾角增大,液滴运动速度和润湿长度都明显减小。