Complex fluids at complex surfaces simply complicated

We study wetting and filling of patterned surfaces by a nematic liquid crystal. We focus on three important classes of periodic surfaces: triangular, sinusoidal and rectangular. The results highlight the similarities and differences of nematic wetting of these surfaces and wetting by simple fluids. The interplay of geometry, surface and elastic energies can lead to the suppression of either filling or wetting. The periodic rectangular surface displays re-entrant transitions, with a sequence dry-filled-wet-filled, in the relevant region of parameter space.

Re-entrant phase behaviour of network fluids: A patchy particle model with temperature-dependent valence

We study a model consisting of particles with dissimilar bonding sites ("patches"), which exhibits self-assembly into chains connected by Y-junctions, and investigate its phase behaviour by both simulations and theory. We show that, as the energy cost epsilon(j) of forming Y-junctions increases, the extent of the liquid-vapour coexistence region at lower temperatures and densities is reduced. The phase diagram thus acquires a characteristic "pinched" shape in which the liquid branch density decreases as the temperature is lowered. To our knowledge, this is the first model in which the predicted topological phase transition between a fluid composed of short chains and a fluid rich in Y-junctions is actually observed. Above a certain threshold for epsilon(j), condensation ceases to exist because the entropy gain of forming Y-junctions can no longer offset their energy cost. We also show that the properties of these phase diagrams can be understood in terms of a temperature-dependent effective valence of the patchy particles. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3605703]

Reentrant Phase Diagram of Network Fluids

We introduce a microscopic model for particles with dissimilar patches which displays an unconventional "pinched'' phase diagram, similar to the one predicted by Tlusty and Safran in the context of dipolar fluids [Science 290, 1328 (2000)]. The model-based on two types of patch interactions, which account, respectively, for chaining and branching of the self-assembled networks-is studied both numerically via Monte Carlo simulations and theoretically via first-order perturbation theory. The dense phase is rich in junctions, while the less-dense phase is rich in chain ends. The model provides a reference system for a deep understanding of the competition between condensation and self-assembly into equilibrium-polymer chains.

Three-dimensional patchy lattice model for empty fluids

The phase diagram of a simple model with two patches of type A and ten patches of type B (2A10B) on the face centred cubic lattice has been calculated by simulations and theory. Assuming that there is no interaction between the B patches the behavior of the system can be described in terms of the ratio of the AB and AA interactions, r. Our results show that, similarly to what happens for related off-lattice and two-dimensional lattice models, the liquid-vapor phase equilibria exhibit reentrant behavior for some values of the interaction parameters. However, for the model studied here the liquid-vapor phase equilibria occur for values of r lower than 1/3, a threshold value which was previously thought to be universal for 2AnB models. In addition, the theory predicts that below r = 1/3 (and above a new condensation threshold which is < 1/3) the reentrant liquid-vapor equilibria are so extreme that it exhibits a closed loop with a lower critical point, a very unusual behavior in single-component systems. An order-disorder transition is also observed at higher densities than the liquid-vapor equilibria, which shows that the liquid-vapor reentrancy occurs in an equilibrium region of the phase diagram. These findings may have implications in the understanding of the condensation of dipolar hard spheres given the analogy between that system and the 2AnB models considered here. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4771591]

Three-Dimensional patchy lattice model for empty fluids

The phase diagram of a simple model with two patches of type A and ten patches of type B (2A10B) on the face centred cubic lattice has been calculated by simulations and theory. Assuming that there is no interaction between the B patches the behavior of the system can be described in terms of the ratio of the AB and AA interactions, r. Our results show that, similarly to what happens for related off-lattice and two-dimensional lattice models, the liquid-vapor phase equilibria exhibit reentrant behavior for some values of the interaction parameters. However, for the model studied here the liquid-vapor phase equilibria occur for values of r lower than 1/3, a threshold value which was previously thought to be universal for 2AnB models. In addition, the theory predicts that below r = 1/3 (and above a new condensation threshold which is < 1/3) the reentrant liquid-vapor equilibria are so extreme that it exhibits a closed loop with a lower critical point, a very unusual behavior in single-component systems. An order-disorder transition is also observed at higher densities than the liquid-vapor equilibria, which shows that the liquid-vapor reentrancy occurs in an equilibrium region of the phase diagram. These findings may have implications in the understanding of the condensation of dipolar hard spheres given the analogy between that system and the 2AnB models considered here.

An electro-rheological study of the nematic liquid crystal 4-n-heptyl-4'-cyanobiphenyl

An experimental and theoretical study of the electro-rheological effects observed in the nematic phase of 4-n-heptyl-4'-cyanobiphenyl has been conducted. This liquid crystal appears to be a model system, in which the observed rheological behaviour can be interpreted by the Leslie-Ericksen continuum theory for low molecular weight liquid crystals. Flow curves are illustrated at different temperatures and under the influence of an external electric field ranging from 0 to 3 kV mm-1, applied perpendicular to the direction of flow. Also presented is the apparent viscosity as a function of temperature, over similar values of electric field, obtained at different shear rates. A master flow curve has been constructed for each temperature by dividing the shear rate by the square of the electric field and multiplying by the square of a reference value of electric field. In a log-log plot, two Newtonian plateaux are found to appear at low and high shear rates, connected by a shear-thinning region. We have applied the Leslie-Ericksen continuum theory, in which the director alignment angle is a function of the electric field and the flow field boundary conditions are neglected, to determine viscoelastic parameters and the dielectric anisotropy.

Equilibrium self-assembly of colloids with distinct interaction sites: Thermodynamics, percolation, and cluster distribution functions

We calculate the equilibrium thermodynamic properties, percolation threshold, and cluster distribution functions for a model of associating colloids, which consists of hard spherical particles having on their surfaces three short-ranged attractive sites (sticky spots) of two different types, A and B. The thermodynamic properties are calculated using Wertheim's perturbation theory of associating fluids. This also allows us to find the onset of self-assembly, which can be quantified by the maxima of the specific heat at constant volume. The percolation threshold is derived, under the no-loop assumption, for the correlated bond model: In all cases it is two percolated phases that become identical at a critical point, when one exists. Finally, the cluster size distributions are calculated by mapping the model onto an effective model, characterized by a-state-dependent-functionality (f) over bar and unique bonding probability (p) over bar. The mapping is based on the asymptotic limit of the cluster distributions functions of the generic model and the effective parameters are defined through the requirement that the equilibrium cluster distributions of the true and effective models have the same number-averaged and weight-averaged sizes at all densities and temperatures. We also study the model numerically in the case where BB interactions are missing. In this limit, AB bonds either provide branching between A-chains (Y-junctions) if epsilon(AB)/epsilon(AA) is small, or drive the formation of a hyperbranched polymer if epsilon(AB)/epsilon(AA) is large. We find that the theoretical predictions describe quite accurately the numerical data, especially in the region where Y-junctions are present. There is fairly good agreement between theoretical and numerical results both for the thermodynamic (number of bonds and phase coexistence) and the connectivity properties of the model (cluster size distributions and percolation locus).

Percolation of colloids with distinct interaction sites

We generalize the Flory-Stockmayer theory of percolation to a model of associating (patchy) colloids, which consists of hard spherical particles, having on their surfaces f short-ranged-attractive sites of m different types. These sites can form bonds between particles and thus promote self-assembly. It is shown that the percolation threshold is given in terms of the eigenvalues of a m x m matrix, which describes the recursive relations for the number of bonded particles on the ith level of a cluster with no loops; percolation occurs when the largest of these eigenvalues equals unity. Expressions for the probability that a particle is not bonded to the giant cluster, for the average cluster size and the average size of a cluster to which a randomly chosen particle belongs, are also derived. Explicit results for these quantities are computed for the case f = 3 and m = 2. We show how these structural properties are related to the thermodynamics of the associating system by regarding bond formation as a (equilibrium) chemical reaction. This solution of the percolation problem, combined with Wertheim's thermodynamic first-order perturbation theory, allows the investigation of the interplay between phase behavior and cluster formation for general models of patchy colloids.

Projecto de entreposto frigorífico a R-744/R-717 para produtos congelados

Pretende-se com este trabalho de Projecto de Mestrado conceber uma câmara frigorífica de 20000 m3 para armazenar produtos congelados e estudar um sistema frigorífico com dois fluidos frigorigénios que serão o CO2 (R-744) e o NH3 (R-717) O trabalho inicia-se com a definição dos objectivos principais para o projecto de um sistema frigorífico deste tipo. Após a definição dos objectivos, o projecto propõe um estudo termodinâmico do CO2 (R-744) como fluido frigorigénio, definindo se o seu historial de utilização, características principais, diagrama pressão-entalpia com a distinção das diversas fases do fluído, comparação em diversos parâmetros com outros fluidos, o porquê da utilização deste fluído, problemas comuns em sistemas com a presença deste fluído, entre outros parâmetros de estudo. De seguida será feito o dimensionamento de uma câmara frigorífica de 20000 m3 para armazenar produtos congelados paletizados através de um balanço térmico manual e um balanço térmico informático através do programa da Centauro comparando no final o resultado dos dois métodos. Será feito uma descrição das características básicas da câmara frigorífica com um esquema simples do edifício. Será abordado de seguida o sistema frigorífico a CO2 (R-744) e a NH3 (R-717) nas suas diversas características, nomeadamente no cálculo de caudais, diagrama pressão-entalpia, dimensionamento de tubagens e finalmente selecção de equipamento. Procede-se assim também ao estudo um sistema apenas a NH3 (R-717) com as suas diversas características nomeadamente no cálculo de caudais, diagrama pressão-entalpia, dimensionamento de tubagens e finalmente selecção dos equipamentos mais importantes, para que desta forma se proceda à comparação dos dois sistemas a nível energético, consumo eléctrico e manutenção de equipamentos. Finalmente proceder-se-á interpretação de resultados com o objectivo final de escolher a melhor solução nos vários parâmetros de comparação, para o esquema em questão.

Nematic wetting and filling of crenellated surfaces

We investigate nematic wetting and filling transitions of crenellated surfaces (rectangular gratings) by numerical minimization of the Landau-de Gennes free energy as a function of the anchoring strength, for a wide range of the surface geometrical parameters: depth, width, and separation of the crenels. We have found a rich phase behavior that depends in detail on the combination of the surface parameters. By comparison to simple fluids, which undergo a continuous filling or unbending transition, where the surface changes from a dry to a filled state, followed by a wetting or unbinding transition, where the thickness of the adsorbed fluid becomes macroscopic and the interface unbinds from the surface, nematics at crenellated surfaces reveal an intriguingly rich behavior: in shallow crenels only wetting is observed, while in deep crenels, only filling transitions occur; for intermediate surface geometrical parameters, a new class of filled states is found, characterized by bent isotropic-nematic interfaces, which persist for surfaces structured on large scales, compared to the nematic correlation length. The global phase diagram displays two wet and four filled states, all separated by first-order transitions. For crenels in the intermediate regime re-entrant filling transitions driven by the anchoring strength are observed.

The condensation and ordering of models of empty liquids

We consider a simple model consisting of particles with four bonding sites ("patches"), two of type A and two of type B, on the square lattice, and investigate its global phase behavior by simulations and theory. We set the interaction between B patches to zero and calculate the phase diagram as the ratio between the AB and the AA interactions, epsilon(AB)*, varies. In line with previous work, on three-dimensional off-lattice models, we show that the liquid-vapor phase diagram exhibits a re-entrant or "pinched" shape for the same range of epsilon(AB)*, suggesting that the ratio of the energy scales - and the corresponding empty fluid regime - is independent of the dimensionality of the system and of the lattice structure. In addition, the model exhibits an order-disorder transition that is ferromagnetic in the re-entrant regime. The use of low-dimensional lattice models allows the simulation of sufficiently large systems to establish the nature of the liquid-vapor critical points and to describe the structure of the liquid phase in the empty fluid regime, where the size of the "voids" increases as the temperature decreases. We have found that the liquid-vapor critical point is in the 2D Ising universality class, with a scaling region that decreases rapidly as the temperature decreases. The results of simulations and theoretical analysis suggest that the line of order-disorder transitions intersects the condensation line at a multi-critical point at zero temperature and density, for patchy particle models with a re-entrant, empty fluid, regime. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3657406]

Phase diagrams of particles with dissimilar patches: X-junctions and Y-junctions

We use Wertheim's first-order perturbation theory to investigate the phase behaviour and the structure of coexisting fluid phases for a model of patchy particles with dissimilar patches (two patches of type A and f(B) patches of type B). A patch of type alpha = {A, B} can bond to a patch of type beta = {A, B} in a volume nu(alpha beta), thereby decreasing the internal energy by epsilon(alpha beta). We analyse the range of model parameters where AB bonds, or Y-junctions, are energetically disfavoured (epsilon(AB) < epsilon(AA)/2) but entropically favoured (nu(AB) >> nu(alpha alpha)), and BB bonds, or X-junctions, are energetically favoured (epsilon(BB) > 0). We show that, for low values of epsilon(BB)/epsilon(AA), the phase diagram has three different regions: (i) close to the critical temperature a low-density liquid composed of long chains and rich in Y-junctions coexists with a vapour of chains; (ii) at intermediate temperatures there is coexistence between a vapour of short chains and a liquid of very long chains with X-and Y-junctions; (iii) at low temperatures an ideal gas coexists with a high-density liquid with all possible AA and BB bonds formed. It is also shown that in region (i) the liquid binodal is reentrant (its density decreases with decreasing temperature) for the lower values of epsilon(BB)/epsilon(AA). The existence of these three regions is a consequence of the competition between the formation of X- and Y-junctions: X-junctions are energetically favoured and thus dominate at low temperatures, whereas Y-junctions are entropically favoured and dominate at higher temperatures.

Filling and wetting transitions of nematic liquid crystals on sinusoidal substrates

Close to sinusoidal substrates, simple fluids may undergo a filling transition, in which the fluid passes from a dry to a filled state, where the interface remains unbent but bound to the substrate. Increasing the surface field, the interface unbinds and a wetting transition occurs. We show that this double-transition sequence may be strongly modified in the case of ordered fluids, such as nematic liquid crystals. Depending on the preferred orientation of the nematic molecules at the structured substrate and at the isotropic-nematic interface, the filling transition may not exist, and the fluid passes directly from a dry to a complete-wet state, with the interface far from the substrate. More interestingly, in other situations, the complete wetting transition may be prevented, and the fluid passes from a dry to a filled state, and remains in this configuration, with the interface always attached to the substrate, even for very large surface fields. Both transitions are observed only for a same substrate in a narrow range of amplitudes.

On the solubility of three disperse anthraquinone dyes in supercritical carbon dioxide: new experimental data and correlation

Solubility measurements of quinizarin. (1,4-dihydroxyanthraquinone), disperse red 9 (1-(methylamino) anthraquinone), and disperse blue 14 (1,4-bis(methylamino)anthraquinone) in supercritical carbon dioxide (SC CO2) were carried out in a flow type apparatus, at a temperature range from (333.2 to 393.2) K and at pressures from (12.0 to 40.0) MPa. Mole fraction solubility of the three dyes decreases in the order quinizarin (2.9 x 10(-6) to 2.9.10(-4)), red 9 (1.4 x 10(-6) to 3.2 x 10(-4)), and blue 14 (7.8 x 10(-8) to 2.2 x 10(-5)). Four semiempirical density based models were used to correlatethe solubility of the dyes in the SC CO2. From the correlation results, the total heat of reaction, heat of vaporization plus the heat of solvation of the solute, were calculated and compared with the results presented in the literature. The solubilities of the three dyes were correlated also applying the Soave-Redlich-Kwong cubic equation of state (SRK CEoS) with classical mixing rules, and the physical properties required for the modeling were estimated and reported.

Patching up dipoles: can dipolar particles be viewed as patchy colloids?

We investigate whether the liquid-vapour phase transition of strongly dipolar fluids can be understood using a model of patchy colloids. These consist of hard spherical particles with three short-ranged attractive sites (patches) on their surfaces. Two of the patches are of type A and one is of type B. Patches A on a particle may bond either to a patch A or to a patch B on another particle. Formation of an AA (AB) bond lowers the energy by epsilon AA (epsilon AB). In the limit [image omitted], this patchy model exhibits condensation driven by AB-bonds (Y-junctions). Y-junctions are also present in low-density, strongly dipolar fluids, and have been conjectured to play a key role in determining their critical behaviour. We map the dipolar Yukawa hard-sphere (DYHS) fluid onto this 2A + 1B patchy model by requiring that the latter reproduce the correct DYHS critical point as a function of the isotropic interaction strength epsilon Y. This is achieved for sensible values of epsilon AB and the bond volumes. Results for the internal energy and the particle coordination number are in qualitative agreement with simulations of DYHSs. Finally, by taking the limit [image omitted], we arrive at a new estimate for the critical point of the dipolar hard-sphere fluid, which agrees with extrapolations from simulation.