Lattice dynamics study of microstructures of electrorheological fluids

The lattice dynamics method is used to study the stability of the chain structures formed in electrorheological (ER) fluids. The appearance of the soft modes in the phonon dispersion of the structures indicates that the chains tend to distort and aggregate into thicker columns due to the electrostatic attractive forces and thermal generated forces between them. The results show that the stability of the chains relies on their width and the separation between them. The complete chain structures are more stable than the chains with defects. The results can be used to elucidate the densification phenomenon of the chains in the structuring process of ER fluids in the quiescent state.

Dynamic response times of electrorheological fluids in steady shear

Transient responses of electrorheological fluids to square-wave electric fields in steady shear are investigated by computational simulation method. The structure responses of the fluids to the field with high frequency are found to be very similar to that to the field with very low frequency or the sudden applied direct current field. The stress rise processes are also similar in both cases and can be described by an exponential expression. The characteristic time tau of the stress response is found to decrease with the increase of the shear rate (gamma) over dot and the area fraction of the particles phi(2). The relation between them can be roughly expressed as tau proportional to(gamma) over dot(-3/4)phi(2)(-3/2). The simulation results are compared with experimental measurements. The aggregation kinetics of the particles in steady shear is also discussed according to these results.

Computer simulation study of equilibrium properties of magnetic dipolar fluids

Langevin dynamics simulations are used to investigate the equilibrium magnetization properties and structure of magnetic dipolar fluids. The influence of using different boundary conditions are systematically studied. Simulation results on the initial susceptibility and magnetization curves are compared with theoretical predictions. The effect of particle aggregation is discussed in detail by performing a cluster analysis of the microstructure.

Algebraically special coplanar shear-free perfect fluids in general relativity

We investigate all algebraically special, not conformally flat, shear-free, isentropic (<i>p</i>(<i>w</i>), <i>w</i> + <i>p</i> &ne; 0), perfect fluid solutions of Einstein's field equations. We show, using the GHP formalism, that if the repeated principle null direction of the Weyl tensor is coplanar with the fluid's 4-velocity and vorticity vector (assumed nonzero), then the fluid's expansion must vanish.<br />

Shear- free perfect fluids with solenoidal magnetic curvature and a ��-law equation of state

We show that shear-free perfect fluids obeying an equation of state p = (&gamma; &minus; 1)&mu; are non-rotating or non-expanding under the assumption that the spatial divergence of the magnetic part of the Weyl tensor is zero.<br />

Shear-free perfect fluids with a solenoidal magnetic curvature

We investigate shear-free perfect fluid solutions of the Einstein field equations where the fluid pressure satisfies a barotropic equation of state and the spatial divergence of the magnetic part of the Weyl tensor is zero. We prove, with the exception of certain quite restricted special cases within the class of solutions in which there exists a Killing vector aligned with the vorticity and for which the magnitude of the vorticity &omega; is not a function of the matter density &mu; alone, that such a fluid is either non-rotating or non-expanding. In the restricted cases the equation of state must satisfy an over-determined differential system.<br />

Shear-free perfect fluids with a solenoidal electric curvature

We prove that the vorticity or the expansion vanishes for any shear-free perfect fluid solution of the Einstein field equations where the pressure satisfies a barotropic equation of state and the spatial divergence of the electric part of the Weyl tensor is zero.<br />

Analog Model for Quantum Gravity Effects: Phonons in Random Fluids

Coordena����o de Aperfei��oamento de Pessoal de N��vel Superior (CAPES)

Plane waves in noncommutative fluids

We study the dynamics of the noncommutative fluid in the Snyder space perturbatively at the first order in powers of the noncommutative parameter. The linearized noncommutative fluid dynamics is described by a system of coupled linear partial differential equations in which the variables are the fluid density and the fluid potentials. We show that these equations admit a set of solutions that are monochromatic plane waves for the fluid density and two of the potentials and a linear function for the third potential. The energy-momentum tensor of the plane waves is calculated. �� 2013 Elsevier B.V.

MOLECULAR DYNAMICS STUDIES OF SIMPLE MODEL FLUIDS AND WATER CONFINED IN CARBON NANOTUBE

Molecular Dynamics (MD) simulation is one of the most important computational techniques with broad applications in physics, chemistry, chemical engineering, materials design and biological science. Traditional computational chemistry refers to quantum calculations based on solving Schrodinger equations. Later developed Density Functional Theory (DFT) based on solving Kohn-Sham equations became the more popular ab initio calculation technique which could deal with ~1000 atoms by explicitly considering electron interactions. In contrast, MD simulation based on solving classical mechanics equations of motion is a totally different technique in the field of computational chemistry. Electron interactions were implicitly included in the empirical atom-based potential functions and the system size to be investigated can be extended to ~106 atoms. The thermodynamic properties of model fluids are mainly determined by macroscopic quantities, like temperature, pressure, density. The quantum effects on thermodynamic properties like melting point, surface tension are not dominant. In this work, we mainly investigated the melting point, surface tension (liquid-vapor and liquid-solid) of model fluids including Lennard-Jones model, Stockmayer model and a couple of water models (TIP4P/Ew, TIP5P/Ew) by means of MD simulation. In addition, some new structures of water confined in carbon nanotube were discovered and transport behaviors of water and ions through nano-channels were also revealed.

Realistic fluids as source for dynamically accreting black holes in a cosmological background

We show that a single imperfect fluid can be used as a source to obtain a mass-varying black hole in an expanding universe. This approach generalizes the well-known McVittie spacetime, by allowing the mass to vary thanks to a novel mechanism based on the presence of a temperature gradient. This fully dynamical solution, which does not require phantom fields or fine-tuning, is a step forward in a new direction in the study of systems whose local gravitational attraction is coupled to the expansion of the universe. We present a simple but instructive example for the mass function and briefly discuss the structure of the apparent horizons and the past singularity.

Consistency of the mass variation formula for black holes accreting cosmological fluids

We address the spherical accretion of generic fluids onto black holes. We show that, if the black hole metric satisfies certain conditions, in the presence of a test fluid it is possible to derive a fully relativistic prescription for the black hole mass variation. Although the resulting equation may seem obvious due to a form of it appearing as a step in the derivation of the Schwarzschild metric, this geometrical argument is necessary to fix the added degree of freedom one gets for allowing the mass to vary with time. This result has applications on cosmological accretion models and provides a derivation from first principles to serve as a basis to the accretion equations already in use in the literature.

Computer simulations of colloidal fluids in confinement

Computer-Simulationen von Kolloidalen Fluiden in Beschr��nkten Geometrien Kolloidale Suspensionen, die einen Phasen��bergang aufweisen, zeigen eine Vielfalt an interessanten Effekten, sobald sie auf eine bestimmte Geometrie beschr��nkt werden, wie zum Beispiel auf zylindrische Poren, sph��rische Hohlr��ume oder auf einen Spalt mit ebenen W��nden. Der Einfluss dieser verschiedenen Geometrietypen sowohl auf das Phasenverhalten als auch auf die Dynamik von Kolloid-Polymer-Mischungen wird mit Hilfe von Computer-Simulationen unter Verwendung des Asakura-Oosawa- Modells, f��r welches auf Grund der ���Depletion���-Kr��fte ein Phasen��bergang existiert, untersucht. Im Fall von zylindrischen Poren sieht man ein interessantes Phasenverhalten, welches vom eindimensionalen Charakter des Systems hervorgerufen wird. In einer kurzen Pore findet man im Bereich des Phasendiagramms, in dem das System typischerweise entmischt, entweder eine polymerreiche oder eine kolloidreiche Phase vor. Sobald aber die L��nge der zylindrischen Pore die typische Korrelationsl��nge entlang der Zylinderachse ��berschreitet, bilden sich mehrere quasi-eindimensionale Bereiche der polymerreichen und der kolloidreichen Phase, welche von nun an koexistieren. Diese Untersuchungen helfen das Verhalten von Adsorptionshysteresekurven in entsprechenden Experimenten zu erkl��ren. Wenn das Kolloid-Polymer-Modellsystem auf einen sph��rischen Hohlraum eingeschr��nkt wird, verschiebt sich der Punkt des Phasen��bergangs von der polymerreichen zur kolloidreichen Phase. Es wird gezeigt, dass diese Verschiebung direkt von den Benetzungseigenschaften des Systems abh��ngt, was die Beobachtung von zwei verschiedenen Morphologien bei Phasenkoexistenz erm��glicht ��� Schalenstrukturen und Strukturen des Janustyps. Im Rahmen der Untersuchung von heterogener Keimbildung von Kristallen innerhalb einer Fl��ssigkeit wird eine neue Simulationsmethode zur Berechnung von Freien Energien der Grenzfl��che zwischen Kristall- bzw. Fl��ssigkeitsphase undWand pr��sentiert. Die Resultate f��r ein System von harten Kugeln und ein System einer Kolloid- Polymer-Mischung werden anschlie��end zur Bestimmung von Kontaktwinkeln von Kristallkeimen an W��nden verwendet. Die Dynamik der Phasenseparation eines quasi-zweidimensionalen Systems, welche sich nach einem Quench des Systems aus dem homogenen Zustand in den entmischten Zustand ausbildet, wird mit Hilfe von einer mesoskaligen Simulationsmethode (���Multi Particle Collision Dynamics���) untersucht, die sich f��r eine detaillierte Untersuchung des Einflusses der hydrodynamischen Wechselwirkung eignet. Die Exponenten universeller Potenzgesetze, die das Wachstum der mittleren Dom��nengr����e beschreiben, welche f��r rein zwei- bzw. dreidimensionale Systeme bekannt sind, k��nnen f��r bestimmte Parameterbereiche nachgewiesen werden. Die unterschiedliche Dynamik senkrecht bzw. parallel zu den W��nden sowie der Einfluss der Randbedingungen f��r das L��sungsmittel werden untersucht. Es wird gezeigt, dass die daraus resultierende Abschirmung der hydrodynamischen Wechselwirkungsreichweite starke Auswirkungen auf das Wachstum der mittleren Dom��nengr����e hat.

Strong Dynamical Heterogeneity and Universal Scaling in Driven Granular Fluids

Large-scale simulations of two-dimensional bidisperse granular fluids allow us to determine spatial correlations of slow particles via the four-point structure factor S-4 (q, t). Both cases, elastic (epsilon = 1) and inelastic (epsilon < 1) collisions, are studied. As the fluid approaches structural arrest, i.e., for packing fractions in the range 0.6 <= phi <= 0.805, scaling is shown to hold: S-4 (q, t)/chi(4)(t) = s(q xi(t)). Both the dynamic susceptibility chi(4)(tau(alpha)) and the dynamic correlation length xi(tau(alpha)) evaluated at the alpha relaxation time tau(alpha) can be fitted to a power law divergence at a critical packing fraction. The measured xi(tau(alpha)) widely exceeds the largest one previously observed for three-dimensional (3d) hard sphere fluids. The number of particles in a slow cluster and the correlation length are related by a robust power law, chi(4)(tau(alpha)) approximate to xi(d-p) (tau(alpha)), with an exponent d - p approximate to 1.6. This scaling is remarkably independent of epsilon, even though the strength of the dynamical heterogeneity at constant volume fraction depends strongly on epsilon.