Researcher profile

Pep Español

Pep Español contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
10works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2016arXiv

Assignment of fields from particles to mesh

In Computational Fluid Dynamics there have been many attempts to combine the power of a fixed mesh on which to carry out spatial calculations with that of a set of particles that moves following the velocity field. These ideas indeed go back to Particle-in-Cell methods, proposed about 60 years ago. Of course, some procedure is needed to transfer field information between particles and mesh. There are many possible choices for this "assignment", or "projection". Several requirements may guide this choice. Two well-known ones are conservativity and stability, which apply to volume integrals of the fields. An additional one is here considered: preservation of information. This means that mesh interpolation, followed by mesh assignment, should leave the field values invariant. The resulting methods are termed "mass" assignments due to their strong similarities with the Finite Element Method. We test several procedures, including the well-known FLIP, on three scenarios: simple 1D convection, 2D convection of Zalesak's disk, and a CFD simulation of the Taylor-Green periodic vortex sheet. The most symmetric mass assignment is seen to be clearly superior to other methods.

preprint2015arXiv

Coupling a nano-particle with isothermal fluctuating hydrodynamics: Coarse-graining from microscopic to mesoscopic dynamics

We derive a coarse-grained description of the dynamics of a nanoparticle immersed in an isothermal simple fluid by performing a systematic coarse graining of the underlying microscopic dynamics. As coarse-grained or relevant variables we select the position of the nanoparticle and the \emph{total} mass and momentum density field of the fluid, which are locally conserved slow variables because they are defined to include the contribution of the nanoparticle. The theory of coarse graining based on the Zwanzing projection operator leads us to a system of stochastic \emph{ordinary} differential equations (SODEs) that are closed in the relevant variables. We demonstrate that our discrete coarse-grained equations are consistent with a Petrov-Galerkin finite-element discretization of a system of formal stochastic \emph{partial} differential equations (SPDEs) which resemble previously-used phenomenological models based on fluctuating hydrodynamics. Under suitable approximations we obtain \emph{closed} approximations of the coarse-grained dynamics in a manner which gives them a clear physical interpretation, and provides \emph{explicit} microscopic expressions for all of the coefficients appearing in the closure. Our work leads to a model for dilute nanocolloidal suspensions that can be simulated effectively using feasibly short molecular dynamics simulations as input to a FEM fluctuating hydrodynamic solver.

preprint2015arXiv

Finite element discretization of non-linear diffusion equations with thermal fluctuations

We present a finite element discretization of a non-linear diffusion equation used in the field of critical phenomena and, more recently, in the context of Dynamic Density Functional Theory. The discretized equation preserves the structure of the continuum equation. Specifically, it conserves the total number of particles and fulfills an H-theorem as the original partial differential equation. Guided by the Theory of Coarse-Graining, we discuss the inclusion of thermal fluctuations in the non-linear diffusion equation. This sheds light on the meaning of such a fluctuating hydrodynamics equation and to the limitations of the approximations involved. The methodology proposed for the introduction of thermal fluctuations in finite element methods is general and valid for both regular and irregular grids in arbitrary dimensions. We focus here on simulations of the Ginzburg-Landau free energy functional using both regular and irregular 1D grids. Convergence of the numerical results is obtained for the static and dynamic structure factors as the resolution of the grid is increased.

preprint2014arXiv

Statistical Mechanics of Hamiltonian Adaptive Resolution Simulations

The Adaptive Resolution Scheme (AdResS) is a hybrid scheme that allows one to treat a molecular system with different levels of resolution depending on the location of the molecules. The construction of a Hamiltonian based on the this idea (H-AdResS) allows one to formulate the usual tools of ensembles and statistical mechanics. We present a number of exact and approximate results that provide a statistical mechanics foundation for this simulation method. We also present simulation results that illustrate the theory.

preprint2013arXiv

Comment on "Logarithmic Oscillators: Ideal Hamiltonian Thermostats" [arXiv 1203.5968]

Campisi, Zhan, Talkner and Hänggi have recently proposed a novel Hamiltonian thermostat which they claim may be used both in simulations and experiments [arXiv:1203.5968v4]. We show, however, that this is not possible due to the length and time scales involved, which depend exponentially on the total energy of the system. The implementation suggested by Campisi et alii implies equilibration times greater than the age of the universe for systems with more than a few dozen particles.

preprint2013arXiv

Monte Carlo adaptive resolution simulation of multicomponent molecular liquids

Complex soft matter systems can be efficiently studied with the help of adaptive resolution simulation methods, concurrently employing two levels of resolution in different regions of the simulation domain. The non-matching properties of high- and low-resolution models, however, lead to thermodynamic imbalances between the system's subdomains. Such inhomogeneities can be healed by appropriate compensation forces, whose calculation requires nontrivial iterative procedures. In this work we employ the recently developed Hamiltonian Adaptive Resolution Simulation method to perform Monte Carlo simulations of a binary mixture, and propose an efficient scheme, based on Kirkwood Thermodynamic Integration, to regulate the thermodynamic balance of multi-component systems.

preprint2001arXiv

Thermohydrodynamics for a van der Waals fluid

Following a cell-method due to van Kampen for the calculation of a coarse-grained free energy functional for the van der Waals gas, we compute a corresponding entropy functional from microscopic principles. This entropy functional is one of the building blocks of the recently developed GENERIC framework. This framework allows to obtain in a thermodynamically consistent way the continuum hydrodynamic equations for a fluid able to display liquid-vapor coexistence. Surface tension appears naturally and the resulting model describes interfaces as diffuse regions, much in the same spirit as the gradient theory for equilibrium situations. We suggest that using interfacial forces in the integral form obtained in the microscopic derivation instead of third order derivatives of the density field might represent an advantage from a computational point of view.

preprint1997arXiv

Dissipative Particle Dynamics with energy conservation

Dissipative particle dynamics (DPD) does not conserve energy and this precludes its use in the study of thermal processes in complex fluids. We present here a generalization of DPD that incorporates an internal energy and a temperature variable for each particle. The dissipation induced by the dissipative forces between particles is invested in raising the internal energy of the particles. Thermal conduction occurs by means of (inverse) temperature differences. The model can be viewed as a simplified solver of the fluctuating hydrodynamic equations and opens up the possibility of studying thermal processes in complex fluids with a mesoscopic simulation technique.

preprint1997arXiv

Stochastic differential equations for non-linear hydrodynamics

We formulate the stochastic differential equations for non-linear hydrodynamic fluctuations. The equations incorporate the random forces through a random stress tensor and random heat flux as in the Landau and Lifshitz theory. However, the equations are non-linear and the random forces are non-Gaussian. We provide explicit expressions for these random quantities in terms of the well-defined increments of the Wienner process.