Source author record

A. Santos

A. Santos appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2020arXiv

One-dimensional Janus fluids. Exact solution and mapping from the quenched to the annealed system

The equilibrium properties of a Janus fluid confined to a one-dimensional channel are exactly derived. The fluid is made of particles with two faces (active and passive), so that the pair interaction is that of hard spheres, except if the two active faces are in front of each other, in which case the interaction has a square-well attractive tail. Our exact solution refers to quenched systems (i.e., each particle has a fixed face orientation), but we argue by means of statistical-mechanical tools that the results also apply to annealed systems (i.e., each particle can flip its orientation) in the thermodynamic limit. Comparison between theoretical results and Monte Carlo simulations for quenched and annealed systems, respectively, shows an excellent agreement.

preprint2020arXiv

Structural properties of additive binary hard-sphere mixtures

An approach to obtain the structural properties of additive binary hard-sphere mixtures is presented. Such an approach, which is a nontrivial generalization of the one recently used for monocomponent hard-sphere fluids [S. Pieprzyk, A. C. Brańka, and D. M. Heyes, Phys. Rev. E 95, 062104 (2017)], combines accurate molecular-dynamics simulation data, the pole structure representation of the total correlation functions, and the Ornstein-Zernike equation. A comparison of the direct correlation functions obtained with the present scheme with those derived from theoretical results stemming from the Percus-Yevick (PY) closure and the so-called rational-function approximation (RFA) is performed. The density dependence of the leading poles of the Fourier transforms of the total correlation functions and the decay of the pair correlation functions of the mixtures are also addressed and compared to the predictions of the two theoretical approximations. A very good overall agreement between the results of the present scheme and those of the RFA is found, thus suggesting that the latter (which is an improvement over the PY approximation) can safely be used to predict reasonably well the long-range behavior, including the structural crossover, of the correlation functions of additive binary hard-sphere mixtures.

preprint2016arXiv

Theoretical approaches to the structural properties of the square-shoulder fluid

A comparison of simulation results with the prediction of the structural properties of square-shoulder fluids is carried out to assess the performance of three theories: Tang--Lu's first-order mean spherical approximation, the simplified exponential approximation of the latter and the rational-function approximation. These three theoretical developments share the characteristic of being analytical in Laplace space and of reducing in the proper limit to the Percus--Yevick result for the hard-sphere fluid. Overall, the best agreement with the simulation data is obtained with the simplified exponential approximation.

preprint2015arXiv

Equation of state of a seven-dimensional hard-sphere fluid. Percus-Yevick theory and molecular dynamics simulations

Following the work of Leutheusser [Physica A 127, 667 (1984)], the solution to the Percus-Yevick equation for a seven-dimensional hard-sphere fluid is explicitly found. This allows the derivation of the equation of state for the fluid taking both the virial and the compressibility routes. An analysis of the virial coefficients and the determination of the radius of convergence of the virial series are carried out. Molecular dynamics simulations of the same system are also performed and a comparison between the simulation results for the compressibility factor and theoretical expressions for the same quantity is presented.

preprint2014arXiv

Role of roughness on the hydrodynamic homogeneous base state of inelastic spheres

A gas of inelastic rough spheres admits a spatially homogeneous base state which turns into a hydrodynamic state after a finite relaxation time. We show that this relaxation time is hardly dependent on the degree of inelasticity but increases dramatically with decreasing roughness. An accurate description of translational-rotational velocity correlations at all times is also provided. At a given inelasticity, the roughness parameter can be tuned to produce a huge distortion from the Maxwellian distribution function. The results are obtained from a Grad-like solution of the Boltzmann-Enskog equation complemented by Monte Carlo and molecular dynamics simulations.

preprint2013arXiv

Study of SEY degradation of amorphous carbon coatings

Deposition of low secondary electron yield (SEY) carbon coatings by magnetron sputtering onto the inner walls of the accelerator seems to be the most promising solution for suppressing the electron cloud problem. However, these coatings change their electron emission properties during long term exposure to air. The ageing process of carbon coated samples with initial SEY of about 0.9 received from CERN is studied as a function of exposure to different environments. It is shown that samples having the same initial SEY may age with different rates. The SEY increase can be correlated with the surface concentration of oxygen. Annealing of samples in air at 100-200 °C reduces the ageing rate and even recovers previously degraded samples. The result of annealing is reduction of the hydrogen content in the coatings by triggering its surface segregation followed by desorption.

preprint2012arXiv

Rational-function approximation for fluids interacting via piece-wise constant potentials

The structural properties of fluids whose molecules interact via potentials with a hard-core plus n piece-wise constant sections of different widths and heights are derived using a (semi-analytical) rational-function approximation method. The results are illustrated for the cases of a square-shoulder plus square-well potential and a shifted square-well potential and compared both with simulation data and with those that follow from the (numerical) solutions of the Percus-Yevick integral equation.

preprint2012arXiv

Steady base states for non-Newtonian granular hydrodynamics

We study in this work steady laminar flows in a low density granular gas modelled as a system of identical smooth hard spheres that collide inelastically. The system is excited by shear and temperature sources at the boundaries, which consist of two infinite parallel walls. Thus, the geometry of the system is the same that yields the planar Fourier and Couette flows in standard gases. We show that it is possible to describe the steady granular flows in this system, even at large inelasticities, by means of a (non-Newtonian) hydrodynamic approach. All five types of Couette-Fourier granular flows are systematically described, identifying the different types of hydrodynamic profiles. Excellent agreement is found between our classification of flows and simulation results. Also, we obtain the corresponding non-linear transport coefficients by following three independent and complementary methods: (1) an analytical solution obtained from Grad's 13-moment method applied to the inelastic Boltzmann equation, (2) a numerical solution of the inelastic Boltzmann equation obtained by means of the direct simulation Monte Carlo method and (3) event-driven molecular dynamics simulations. We find that, while Grad's theory does not describe quantitatively well all transport coefficients, the three procedures yield the same general classification of planar Couette-Fourier flows for the granular gas

preprint2012arXiv

The penetrable-sphere fluid in the high-temperature, high-density limit

We consider a fluid of $d$-dimensional spherical particles interacting via a pair potential $ϕ(r)$ which takes a finite value $ε$ if the two spheres are overlapped ($r<σ$) and 0 otherwise. This penetrable-sphere model has been proposed to describe the effective interaction of micelles in a solvent. We derive the structural and thermodynamic functions in the limit where the reduced temperature $k_BT/ε$ and density $ρσ^d$ tend to infinity, their ratio being kept finite. The fluid exhibits a spinodal instability at a certain maximum scaled density where the correlation length diverges and a crystalline phase appears, even in the one-dimensional model. By using a simple free-volume theory for the solid phase of the model, the fluid-solid phase transition is located.

preprint2011arXiv

Hydrodynamics of inelastic Maxwell models

An overview of recent results pertaining to the hydrodynamic description (both Newtonian and non-Newtonian) of granular gases described by the Boltzmann equation for inelastic Maxwell models is presented. The use of this mathematical model allows us to get exact results for different problems. First, the Navier--Stokes constitutive equations with explicit expressions for the corresponding transport coefficients are derived by applying the Chapman--Enskog method to inelastic gases. Second, the non-Newtonian rheological properties in the uniform shear flow (USF) are obtained in the steady state as well as in the transient unsteady regime. Next, an exact solution for a special class of Couette flows characterized by a uniform heat flux is worked out. This solution shares the same rheological properties as the USF and, additionally, two generalized transport coefficients associated with the heat flux vector can be identified. Finally, the problem of small spatial perturbations of the USF is analyzed with a Chapman--Enskog-like method and generalized (tensorial) transport coefficients are obtained.

preprint2011arXiv

Structure of the square-shoulder fluid

The structural properties of square-shoulder fluids are derived from the use of the rational function approximation method. The computation of both the radial distribution function and the static structure factor involves mostly analytical steps, requiring only the numerical solution of a single transcendental equation. The comparison with available simulation data and with numerical solutions of the Percus-Yevick and hypernetted-chain integral equations shows that the present approximation represents an improvement over the Percus-Yevick theory for this system and a reasonable compromise between accuracy and simplicity.

preprint2010arXiv

A Numerical Test of a High-Penetrability Approximation for the One-Dimensional Penetrable-Square-Well Model

The one-dimensional penetrable-square-well fluid is studied using both analytical tools and specialized Monte Carlo simulations. The model consists of a penetrable core characterized by a finite repulsive energy combined with a short-range attractive well. This is a many-body one-dimensional problem, lacking an exact analytical solution, for which the usual van Hove theorem on the absence of phase transition does not apply. We determine a high-penetrability approximation complementing a similar low-penetrability approximation presented in previous work. This is shown to be equivalent to the usual Debye-Hückel theory for simple charged fluids for which the virial and energy routes are identical. The internal thermodynamic consistency with the compressibility route and the validity of the approximation in describing the radial distribution function is assessed by a comparison against numerical simulations. The Fisher-Widom line separating the oscillatory and monotonic large-distance behavior of the radial distribution function is computed within the high-penetrability approximation and compared with the opposite regime, thus providing a strong indication of the location of the line in all possible regimes. The high-penetrability approximation predicts the existence of a critical point and a spinodal line, but this occurs outside the applicability domain of the theory. We investigate the possibility of a fluid-fluid transition by Gibbs ensemble Monte Carlo techniques, not finding any evidence of such a transition. Additional analytical arguments are given to support this claim. Finally, we find a clustering transition when Ruelle's stability criterion is not fulfilled. The consequences of these findings on the three-dimensional phase diagrams are also discussed.

preprint2007arXiv

Third and fourth degree collisional moments for inelastic Maxwell models

The third and fourth degree collisional moments for $d$-dimensional inelastic Maxwell models are exactly evaluated in terms of the velocity moments, with explicit expressions for the associated eigenvalues and cross coefficients as functions of the coefficient of normal restitution. The results are applied to the analysis of the time evolution of the moments (scaled with the thermal speed) in the free cooling problem. It is observed that the characteristic relaxation time toward the homogeneous cooling state decreases as the anisotropy of the corresponding moment increases. In particular, in contrast to what happens in the one-dimensional case, all the anisotropic moments of degree equal to or less than four vanish in the homogeneous cooling state for $d\geq 2$.

preprint2004arXiv

DSMC evaluation of the Navier-Stokes shear viscosity of a granular fluid

A method based on the simple shear flow modified by the introduction of a deterministic non-conservative force and a stochastic process is proposed to measure the Navier-Stokes shear viscosity in a granular fluid described by the Enskog equation. The method is implemented in DSMC simulations for a wide range of values of dissipation and density. It is observed that, after a certain transient period, the system reaches a hydrodynamic stage which tends to the Navier-Stokes regime for long times. The results are compared with theoretical predictions obtained from the Chapman-Enskog method in the leading Sonine approximation, showing quite a good agreement, even for strong dissipation.

preprint2003arXiv

Exact steady state solution of the Boltzmann equation: A driven 1-D inelastic Maxwell gas

The exact nonequilibrium steady state solution of the nonlinear Boltzmann equation for a driven inelastic Maxwell model was obtained by Ben-Naim and Krapivsky [Phys. Rev. E 61, R5 (2000)] in the form of an infinite product for the Fourier transform of the distribution function $f(c)$. In this paper we have inverted the Fourier transform to express $f(c)$ in the form of an infinite series of exponentially decaying terms. The dominant high energy tail is exponential, $f(c)\simeq A_0\exp(-a|c|)$, where $a\equiv 2/\sqrt{1-α^2}$ and the amplitude $A_0$ is given in terms of a converging sum. This is explicitly shown in the totally inelastic limit ($α\to 0$) and in the quasi-elastic limit ($α\to 1$). In the latter case, the distribution is dominated by a Maxwellian for a very wide range of velocities, but a crossover from a Maxwellian to an exponential high energy tail exists for velocities $|c-c_0|\sim 1/\sqrt{q}$ around a crossover velocity $c_0\simeq \ln q^{-1}/\sqrt{q}$, where $q\equiv (1-α)/2\ll 1$. In this crossover region the distribution function is extremely small, $\ln f(c_0)\simeq q^{-1}\ln q$.

preprint2001arXiv

Comment on "Theory and computer simulation for the equation of state of additive hard-disk fluid mixtures"

A flaw in the comparison between two different theoretical equations of state for a binary mixture of additive hard disks and Monte Carlo results, as recently reported in C. Barrio and J. R. Solana, Phys. Rev. E 63, 011201 (2001), is pointed out. It is found that both proposals, which require the equation of state of the single component system as input, lead to comparable accuracy but the one advocated by us [A. Santos, S. B. Yuste, and M. López de Haro, Mol. Phys. 96, 1 (1999)] is simpler and complies with the exact limit in which the small disks are point particles.

preprint2000arXiv

Computer simulation of uniformly heated granular fluids

Direct Monte Carlo simulations of the Enskog-Boltzmann equation for a spatially uniform system of smooth inelastic spheres are performed. In order to reach a steady state, the particles are assumed to be under the action of an external driving force which does work to compensate for the collisional loss of energy. Three different types of external driving are considered: (a) a stochastic force, (b) a deterministic force proportional to the particle velocity and (c) a deterministic force parallel to the particle velocity but constant in magnitude. The Enskog-Boltzmann equation in case (b) is fully equivalent to that of the homogeneous cooling state (where the thermal velocity monotonically decreases with time) when expressed in terms of the particle velocity relative to the thermal velocity. Comparison of the simulation results for the fourth cumulant and the high energy tail with theoretical predictions derived in cases (a) and (b) [T. P. C. van Noije and M. H. Ernst, Gran. Matt. 1, 57 (1998)] shows a good agreement. In contrast to these two cases, the deviation from the Maxwell-Boltzmann distribution is not well represented by Sonine polynomials in case (c), even for low dissipation. In addition, the high energy tail exhibits an underpopulation effect in this case.

preprint2000arXiv

Demixing in binary mixtures of hard hyperspheres

The phase behavior of binary fluid mixtures of hard hyperspheres in four and five dimensions is investigated. Spinodal instability is found by using a recent and accurate prescription for the equation of state of the mixture that requires the equation of state of the single component fluid as input. The role played by the dimensionality on the possible metastability of the demixing transition with respect to a fluid-solid transition is discussed. The binodal curves in the pressure--chemical potential representation are seen to lie on a common line, independent of the size ratio

preprint2000arXiv

Nonlinear Couette flow in a low density granular gas

A model kinetic equation is solved exactly for a special stationary state describing nonlinear Couette flow in a low density system of inelastic spheres. The hydrodynamic fields, heat and momentum fluxes, and the phase space distribution function are determined explicitly. The results apply for conditions such that viscous heating dominates collisional cooling, including large gradients far from the reference homogeneous cooling state. Explicit expressions for the generalized transport coefficients (e.g., viscosity and thermal conductivity) are obtained as nonlinear functions of the coefficient of normal restitution and the shear rate.These exact results for the model kinetic equation are also shown to be good approximations to the corresponding state for the Boltzmann equation via comparison with direct Monte Carlo simulation for the latter

preprint2000arXiv

On the influence of gravity on the thermal conductivity

In this paper we evaluate the corrections to the Navier-Stokes constitutive equations induced by the action of a gravitational field in a gas subjected to a thermal gradient parallel to the field with no convection. The analysis is performed from an exact perturbation solution of the BGK kinetic model for Maxwell molecules through sixth order in the field. The reference state (zeroth order approximation) corresponds to the exact solution in the pure planar Fourier flow, which holds for arbitrary values of the thermal gradient. The results show that the pressure tensor becomes anisotropic, so that the momentum flux along the field direction is enhanced. In addition, the heat flux increases (decreases) with respect to its Navier-Stokes value when the gas is heated from above (below).