Search

Through acoustically levitated ice formation and subsequent release onto a controlled area, this study introduces a third class of ice-countering system beyond de- and anti-icing: ice-impacting. By subjecting stainless steel 316 (SS), epoxy resin-coated (ER), and laser-textured (LT) surfaces with known surface roughness, hardness, and structural characteristics to 40 ice droplet impacts each, the effect on surface properties and their effect on solid-solid interfacial impact dynamics, in turn, was examined using a novel analysis framework based on fundamental conservation laws. For the velocities experienced in this study, the impacts did not affect the surface properties; they were consistent after each impact. Elasticity was the most significant factor in droplet behavior: the ER surface exhibited rebounding for 78% of impacts (important for moving surfaces). Surface roughness also played a role, particularly for droplets with rotational motion, as immobilization occurred for 66% of impacts on the rougher LT surface. However, the nanostructures on that textured surface resulted in droplet redirection perpendicular to the surface directionality (critical for stationary surfaces). In contrast, the other surfaces saw no change or no consistent change in rebound angle. Elasticity also affected momentum retention, where the ER surface had a translational restitution coefficient of 0.32 compared to 0.17 for the two stainless steel surfaces. Surface roughness was the predominant aspect of energy retention: the LT surface had a translational-to-rotational energy transfer coefficient of 0.07 (0.23 for the smoother surfaces), resulting in an overall energy retention coefficient of 0.09 compared to 0.28 for the SS and ER surfaces on average.

(1) Introduction: Nanoparticles have multiple applications, including drug delivery systems, biosensing, and carbon capture. Non-Einstein-like viscosity reduction has been reported in nanoparticle-polymer blends at low nanoparticle concentrations. More recently, a similar non-Einsteinian viscosity reduction effect has been observed in aqueous ultra-low concentration carbon-based nanofluids. (2) Methods: We use a boron nitride nanotube functionalized with hydrophilic groups in rheological experiments to investigate the viscosity reduction in ultra-low concentration nanofluids (0.1-10 ppm). We measure the dynamic viscosity in an air atmosphere and methane (0-5 MPag) at low temperatures (0-10 C). (3) Results: A negligible effect on the temperature dependence of viscosity was found. Ultra-low concentrations of BNNT reduced the viscosity of the nanofluid by up to 29% at 10 ppm in the presence of methane. The results presented here were compared to similar studies on O-GNF and O-MWCNT nanofluids, which also reported significant viscosity reductions. (4) Conclusions: This work identified a non-Einsteinian viscosity reduction in BNNT nanofluids, which was exacerbated by methane dissolved in the nanofluid.

The viscosity of oxygen-functionalized multi-walled carbon nanotube (O-MWCNT) nanofluids was measured for concentrations from 0.1 to 10 ppm under conditions of 0 to 30 MPag pressures and 0 to 10 C temperatures. The presence of O-MWCNTs did not affect the temperature dependence of viscosity but did reduce the effective viscosity of solution due to cumulative hydrogen bond-disrupting surface effects, which overcame internal drag forces. O-MWCNTs added a weak pressure dependence to the viscosity of solution because of their ability to align more with the flow direction as pressure increased. In the liquid to hydrate phase transition, the times to reach the maximum viscosity were faster in O-MWCNT systems compared to the pure water baseline. However, the presence of O-MWCNTs limited the conditions at which hydrates formed as increased nanoparticle collisions in those systems inhibited the formation of critical clusters of hydrate nuclei. The times to viscosity values most relevant to technological applications were minimally 28.02 % (200 mPa s) and 21.08 % (500 mPa s) slower than the baseline, both in the 1 ppm system, even though all systems were faster to the final viscosity. This was attributed to O-MWCNT entanglement, which resulted in a hydrate slurry occurring at lower viscosity values.

Water's viscosity dependence on pressure was also not affected by O-GNFs, except at 10 ppm, where the shuttle effect may have increased the presence of hydrophobic methane bubbles in the solution. Under high pressure, the relative viscosity of the system remained non-Einsteinian at all temperatures except 2C. This may have been because the density anomaly of water was shifted to a colder temperature as the hydrogen bonding network was weaker. The phase transition from liquid to hydrate was identical to that of pure water, indicating that the presence of different stages of growth was not affected by the presence of O-GNF. However, the times to reach a maximum viscosity were faster in O-GNF systems compared to pure water. This said, the hydrate formation limitations inherent to the measurement system were not overcome by the presence of O-GNFs. The times to application-relevant viscosity values were maximized in the 1 ppm system at 49.75 % (200 mPa.s) and 31.93 % (500 mPa.s) faster than the baseline. Therefore, the presence of O-GNFs allowed for shorter times to desired viscosities and at lower driving forces than the baseline, improving the viability of the hydrate technologies to which they can be added.

In this study, using direct digital and infrared imaging techniques, the freezing of up to three simultaneous THF hydrate droplets was investigated for the first time. Nucleation was initiated at the aqueous solution-air interface. Two pseudo-heterogeneous mechanisms created additional nucleation interfaces: one from cavitation effects entraining microbubbles and another from subvisible ice particles, also called hydrate nucleating particles (HNPs), impacting the droplet surface. For systems containing droplets in both the second and third positions, nucleation was statistically simultaneous between all droplets. This effect may have been caused by the high liquid-solid interfacial pressures that developed at nucleation, causing some cracking in the initial hydrate shell around the droplet and releasing additional HNPs (now of hydrate) into the air. During crystallization, the THF hydrate droplets developed a completely white opacity, termed optical clarity loss or OCL. It was suggested that high hydrate growth rates within the droplet resulted in the capture of tiny air bubbles within the solid phase. In turn, light refraction through many smaller bubbles resulted in the OCL. These bubbles created structural inhomogeneities, which may explain how the volumetric expansion of the droplets upon complete solidification was 23.6% compared with 7.4% in pure, stationary THF hydrate systems. Finally, the thermal gradient that developed between the top and bottom of the droplet during melting resulted in a surface tension gradient along the air-liquid interface. In turn, convective cells developed within the droplet, causing it to spin rapidly about the horizontal axis.

The viscosity of methane hydrate slurries with poly(vinylpyrrolidone) (PVP) at 700 and 7000 ppm by weight, molecular weights of 40,000 (PVP40) and 360,000 (PVP360) Da, and shear rates of 400 and 80 1/s, were measured in a high-pressure rheometer with pressures up to 30 MPag and compared to pure water systems. The additives successfully reduced the formation of high-viscosity slurries, but at low concentrations were incapable of delaying hydrate agglomeration at the late growth stage. The average relative time required for PVP40 solutions at 700 ppm to grow to 50 mPa.s was 1.9 times the water reference value, but only 1.2 times to reach 200 mPa.s. Improved inhibition was observed for the higher concentration and higher molecular weight sets, where the relative time to reach 50 mPa.s were 8.2 and 2.6 times the water reference value, respectively. While the additives demonstrated anti-nucleation properties and suppressed crystal growth initially, they accelerated the hydrate clusters agglomeration rate, and potentially weakened the hydrate mechanical properties.

Initially, the acoustic field forced the droplets into an oblate spheroid shape, though the counteracting force of the cooling stream caused them to circularize. Droplet geometry was thus the net result of streaming forces and surface tension at the acoustic boundary layer/air-liquid interface. Nucleation was determined to be neither homogeneous nor heterogeneous but secondary, and thus dependent on the cooling rate and not on the degree of supercooling. It was likely initiated by aerosolized ice particles from the air or from droplets that had already nucleated and broken up. The latter secondary ice production process resulted in multi-drop systems with statistically identical nucleation times. Notably, this meant that the presence of interfacial rupture at an adjacent droplet could influence the crystallization behaviour of another. After the formation of an initial ice shell around the individual droplets, dendritic protrusions grew from the droplet surface, likely seeded by the same ice particles that caused nucleation, but at a quasi-liquid layer. When freezing was complete, it was determined that the frozen core had undergone a volumetric expansion of 30.75%, compared to 9% for pure, sessile water expansion. This significantly greater expansion may have resulted from entrained air bubbles at the inner solid-liquid interface and oscillations at the moving phase boundary caused by changes in local acoustic forces. Soon after melting began, acoustic streaming, the buoyancy of the remaining ice, and convective currents caused by an inner thermal gradient and thermocapillary effects along the air-liquid interface, all contributed to the droplet spinning about the horizontal axis.

(1) Introduction: New technologies that leverage gas hydrates phenomena include carbon capture and sequestrations. These processes are often semi-continuous and require regulation of the system's flow properties for proper operation. Accurate computational models for the viscosity of carbon dioxide hydrate systems at pre-nucleation conditions can be important for process design and control of such technologies. (2) Methods: This work validates the viscosity predictions of molecular dynamics simulations using previously measured experimental data. The TIP4P/Ice force field was used to model water, while the EPM2, TraPPE, and Zhang force fields were used for carbon dioxide. The Green-Kubo and Einstein formulations of viscosity and diffusivity were used in this work. (3) Results: All force fields overpredicted viscosity when compared to experimental data, but EPM2 resulted in lower discrepancies. Additionally, EPM2 was determined to model molecular behavior expected from the macroscopic trends in viscosity with respect to temperature and pressure. (4) Conclusions: The EPM2 force field more accurately predicted the viscosity of carbon dioxide hydrates systems at pre-nucleation conditions relative to TraPPE and Zhang.

The viscosity of methane and carbon dioxide hydrate systems were measured using a high-pressure rheometer up to 30 MPag. Where hydrate formation was not detected, the effect of temperature on the viscosity was one order of magnitude higher than the pressure effect on viscosity in most of the experimental pressure range (-0.048 mPa s/C at 1 MPag and 0.009 mPa s/MPag at 2C). The pressure effect on the viscosity of carbon dioxide systems where no hydrate formation was observed was up to one order of magnitude higher than that of the methane systems, due to carbon dioxide's higher solubility in water. Novel rheological phases diagrams were developed to further characterize the gas hydrate systems. Several systems with high driving forces for hydrate formation (2.07 MPag to 4.1 MPag) did not form gas hydrates. System limitations to the formation of hydrates were categorized as kinetic, mass diffusion, and/or heat of crystallization effects.

A systematic analysis of defect textures in facetted nanoparticles with polygonal configurations embedded in a nematic matrix is performed using the Landau-de Gennes model, homeotropic strong anchoring in a square domain with uniform alignment in the outer boundaries. Defect and textures are analyzed as functions of temperature T, polygon size R, and polygon number N. For nematic nanocomposites, the texture satisfies a defect charge balance equation between bulk and surface (particle corner) charges. Upon decreasing the temperature, the central bulk defects split and together with other -1/2 bulk defects, are absorbed by the nanoparticle's corners. Increasing the lattice size decreases confinement and eliminates bulk defects. Increasing the polygon number increases the central defect charge at high temperature and the number of surface defects at lower temperatures. The excess energy per particle is lower in even than in odd polygons, and it is minimized for a square particle arrangement. These discrete modeling results show for first time that even under strong anchoring, defects are attached to particles as corner defects, leaving behind a low energy homogeneous orientation field that favors nanoparticle ordering in nematic matrices. These new insights are consistent with recent thermodynamic approaches to nematic nanocomposites that predict the existence of novel nematic/crystal phases and can be used to design nanocomposites with orientational and positional order., Comment: Author's version. Published in Langmuir, http://pubs.acs.org/doi/abs/10.1021/la203226g

A new computational strategy is proposed to approximate, with a simple but accurate expression, the Maier- Saupe free energy for nematic order. Instead of the traditional approach of expanding the free energy with a truncated Taylor series, we employ a least-squares fitting to obtain the coefficients of a polynomial expression. Both methods are compared, and the fitting with at most five polynomial terms is shown to provide a satisfactory fitting, and to give much more accurate results than the traditional Taylor expansion. We perform the analysis in terms of the tensor order parameter, so the results are valid in uniaxial and biaxial states., Comment: This is an Author's Original Manuscript of an article Published in Liquid Crustals (copyright Taylor & Francis); available online at http://www.tandfonline.com/doi/abs/10.1080/02678292.2010.539303#.Uo4PVifZ8l0

We present theory and simulation of simultaneous chemical demixing and phase ordering in a polymer-liquid crystal mixture in conditions where isotropic-isotropic phase separation is metastable with respect to isotropic-nematic phase transition. In the case the mechanism is nucleation and growth, it is found that mesophase growth proceeds by a transient metastable phase that surround the ordered phase, and whose lifetime is a function of the ratio of diffusional to orientational mobilities. In the case of spinodal decomposition, different dynamic regimes are observed depending on the mobility ratio: metastable phase separation preceding phase ordering, phase ordering preceding phase separation, or simultaneous phase ordering and phase separation. Not only the overall dynamics but also the final structure of the material can be different for each kinetic regime., Comment: 29 pages, 11 figures. Accepted author{s version, printed version in http://link.springer.com/article/10.1140/epjb/e2011-20375-7. arXiv admin note: substantial text overlap with arXiv:0904.4460

In this work, a theoretical study of phase equilibrium in mixtures of a calamitic nematic liquid crystal and hard spherical nanoparticles is presented. A mean-field thermodynamic model is used, where the interactions are considered to be proportional to the number of contacts, which in turn are proportional to the areas and area fractions of each component. It is shown that, as the radius of the particle is increased, the effect of the particles on the isotropic-nematic transition is less pronounced, and that for large radius the miscibility increases as the particle radius increases., Comment: This is an Author's Accepted Manuscript of an article published in "Molecular Crystals and Liquid Crystals" (see reference), available online at: http://www.tandfonline.com/doi/abs/10.1080/15421406.2011.609447

Liquid crystals are synthetic and biological viscoelastic anisotropic soft matter materials that combine liquid fluidity with crystal anisotropy and find use in optical devices, sensor/actuators, lubrication, super-fibers. Frequently mesogens are mixed with colloidal and nanoparticles, other mesogens, isotropic solvents, thermoplastic polymers, cross-linkable monomers, among others. This comprehensive review present recent progress on meso and macro scale thermodynamic modelling, highlighting the (i) novelties in spinodal and binodal lines in the various phase diagrams, (ii) the growth laws under phase transitions and phase separation, (iii) the ubiquitous role of metastability and its manifestation in complex droplet interfaces, (iv) the various spinodal decompositions due to composition and order fluctuations, (v) the formation of novel material architectures such as colloidal crystals, (vi) the particle rich phase behaviour in liquid crystal nanocomposites, (vii) the use of topological defects to absorb and organize nanoparticles, and (viii) the ability of faceted nanoparticles to link into strings and organize into lattices. Emphasis is given to highlight dominant mechanisms and driving forces, and to link them to specific terms in the free energies of these complex mixtures. The novelties of incorporating mesophases into blends, solutions, dispersions and mixtures is revealed by using theory, modelling, computation, and visualization., Comment: This is an Author's Accepted Manuscript of an article published in "Molecular Simulations" (see reference), available online at: http://www.tandfonline.com/doi/abs/10.1080/08927022.2012.669478

Theory and simulations of simultaneous chemical demixing and phase ordering are performed for a mixed order parameter system with an isotropic-isotropic (I-I) phase separation that is metastable with respect to an isotropic-nematic (I-N) phase ordering transition. Under certain conditions, the disordered phase transforms into an ordered phase via the motion of a double front containing a metastable phase produced by I-I demixing, a thermodynamically driven mechanism not previously reported. Different kinetic regimes are found depending on the location of the initial conditions in the thermodynamic phase diagram and the ratio between diffusional and nematic phase ordering mobilities. For a diffusional process, depending if the temperature is above or below the critical co-dissolution point, an inflection point or a phase separation takes place in the depletion layer. This phase separation leads to the formation of a second interface where the separation of the two metastable isotropic phases grows monotonically with time. The observed deviations from the typical Fickian concentration profiles are associated with strong positive deviations of the mixture from ideality due to couplings between concentration and nematic ordering. Although systems of interest include liquid crystalline nanocomposites, this novel mechanism may apply to any mixture that can undergo an order-disorder transition and demix., Comment: 24 pages, 8 figures, final accepted authors version, published version in http://pre.aps.org/abstract/PRE/v86/i1/e011605

A systematic computational and scaling analysis of defect textures in polygonal arrangement of cylindrical particles embedded in a cholesteric (Ch) liquid crystal matrix is performed using the Landau de Gennes model for chiral self assembly, with strong anchoring at the particles' surface. The defect textures and LC phases observed are investigated as a function of chirality, elastic anisotropy, particle polygonal arrangement and particle size. The presence of a polygonal network made of N circular inclusions results in defect textures of a net charge of -(N-2)/2 per unit polygonal cell. As the chirality increases, the LC matrix shows the following transition sequence: weakly twisted cholesterics, 2D blue phases with non-singular/ singular defect lattices, cholesteric phases with only disclinations, and fingerprint cholesteric textures with disclinations and dislocations. For monomeric mesogens at concentrations far from the I/Ch phase transition and low chirality, for a given symmetry of the LC phase, the particle with weaker (stronger) confinement results in a phase with lower (higher) elastic energy, while at high chirality the elastic energy of a LC phase is proportional to the number of particles that form the polygonal network. Thus, hexagonal (triangular) particle arrangement results in low elastic energy at low (high) chirality. For semiflexible polymeric mesogens (high elastic anisotropy), defect textures with fewer disclinations/ dislocations arise but due to layer distortions we find a higher elastic energy than monomeric mesogens. The defects arising in the simulations and the texture rules established are in agreement with experimental observations in cellulosic liquid crystal analogues such as plant cell wall and helical biological polymeric mesophases made of DNA, PBLG and xanthan.

Experimental and modeling/simulation studies of phase equilibrium and growth morphologies of novel polymer-dispersed liquid crystal (PDLC) mixtures of PS (polystyrene) and liquid crystals that exhibit a direct isotropic/smectic-A (lamellar) mesophase transition were performed for PS/10CB (decyl- cyanobiphenyl) and PS/12CB (dodecyl-cyanobiphenyl). Partial phase diagrams were determined using polarized optical microscopy (POM) and differential scanning calorimetry (DSC) for different compositions of both materials, determining both phase separation (liquid/liquid demixing) and phase ordering (isotropic/smectic-A transition) temperatures. The Flory-Huggins theory of isotropic mixing and Maier-Saupe-McMillan theory for smectic-A liquid crystalline ordering were used to computationally determine phase diagrams for both systems, showing good agreement with the experimental results. In addition to thermodynamic observations, growth morphology relations were found depending on phase transition sequence, quench rate, and material composition. Three stages of liquid crystal-rich domain growth morphology were observed: spherical macroscale domain growth ("stage I"), highly anisotropic domain growth ("stage II"), and sub-micron spheroid domain growth ("stage III"). Nano-scale structure of spheroidal and spherocylindrical morphologies were then determined via two-dimensional simulation of a high-order Landau-de Gennes model. Morphologies observed during stage II growth are typical of di- rect isotropic/smectic-A phase transitions, such as highly anisotropic "batonnets" and filaments. These morphologies, which are found to be persistent in direct isotropic/smectic-A PDLCs, could provide new functionality and applications for these functional materials., Comment: First Revision, 21 pages, 11 figures, submitted to Macromolecules as an article 17JUL2009

An extension to a high-order model for the direct isotropic/smectic-A liquid crystalline phase transition was derived to take into account thermal effects including anisotropic thermal diffusion and latent heat of phase-ordering. Multi-scale multi-transport simulations of the non-isothermal model were compared to isothermal simulation, showing that the presented model extension corrects the standard Landau-de Gennes prediction from constant growth to diffusion-limited growth, under shallow quench/undercooling conditions. Non-isothermal simulations, where meta-stable nematic pre-ordering precedes smectic-A growth, were also conducted and novel non-monotonic phase-transformation kinetics observed., Comment: First revision: 20 pages, 7 figures

We present theory and simulation of simultaneous chemical demixing and phase ordering in a polymer-liquid crystal mixture in conditions where isotropic- isotropic phase separation is metastable with respect to isotropic-nematic phase transition. It is found that mesophase formation proceeds by a transient metastable phase that surround the ordered phase, and whose lifetime is a function of the ratio of diffusional to orientational mobilities. It is shown that kinetic phase ordering in polymer-mesogen mixtures is analogous to kinetic crystallization in polymer solutions., Comment: 17 pages, 5 figures accepted for publication in EPL

Modeling and simulation studies using a high-order Landau-de Gennes type model for the direct isotropic/smectic-A liquid crystal transition are presented showing the existence of meta-stable nematic pre-ordering under certain conditions, in agreement with experimental observations (Tokita, M. et al. Macromolecules 2006, 39, 2021-2023), Comment: submitted for publication in Macromolecules as a Communication to the Editor 11APR2009

A Landau-de Gennes model that integrates the nematic quadrupolar tensor order parameter and complex smectic-A order parameters is used to simulate the two-dimensional growth of an initially homogeneous smectic-A spherulite in an isotropic matrix. These simulations are performed in the shape-dynamic (nano-scale) regime of growth under two material conditions: isotropic nematic elasticity and equal splay-bend nematic elasticity. A comparison of the growth kinetics, spherulite morphology, interfacial/bulk energy landscapes between both cases is made showing that equal nematic splay-bend elasticity is required to reproduce past experimental and theoretical observations. Additionally, a previously unknown undulation instability during spherulite growth is found which, in conjunction with preferred planar anchoring and defect shedding mechanisms at micron length scales, could explain the formation mechanism of focal conic curvature defects and ultimately smectic-A "batonnet" structures observed experimentally., Comment: submitted to Liquid Crystals 09FEB2009

Extensions of a previously presented Landau-de Gennes type liquid crystalline phase transition model for the direct isotropic/smectic-A (lamellar) a phase transition to the direct isotropic/smectic-C (tilted lamellar) transition are studied. Two different proposed extensions to the model are studied both in the context of ideal ordering (point volume) and full three-dimensional uniaxial scalar/vector decomposition. Recommendations based upon the inclusion of essential physics and computational feasibility are made, distinguishing each of the two proposed extensions based upon these criteria. Additionally, it is found that the approach used for the isotropic/smectic-A model to deterministically compute phase diagram data is not possible for either of the two proposed isotropic/smectic-C models., Comment: manuscript submitted for presentation at 8th World Congress of Chemical Engineering August 23-27, 2009, Montreal, Quebec, Canada

A comprehensive modeling and simulation study of the first-order isotropic/smectic-A transition is presented and applied to phase diagram computation and two-dimensional spherulite growth. An approach based on nonlinear optimization, that incorporates experimental data (from 12CB, dodecyl-cyanobiphenyl), is used to determine physically realistic model parameters. These parameters are then used in conjunction with an optimized phase diagram computation method. Additionally, a time-dependent formulation is presented and applied to the study of two-dimensional smectic-A spherulite growth. These results show the growth kinematics and defect dynamics of nanoscale smectic-A spherulite growth in an isotropic phase with an initially radial layer configuration.

Two-dimensional simulations of the coarsening process of the isotropic/smectic-A phase transition are presented using a high-order Landau-de Gennes type free energy model. Defect annihilation laws for smectic disclinations, elementary dislocations, and total dislocation content are determined. The computed evolution of the orientational correlation length and disclination density is found to be in agreement with previous experimental observations showing that disclination interactions dominate the coarsening process. The mechanism of smectic disclination movement, limited by the absorption and emission of elementary dislocations, is found to be facilitated by curvature walls connecting interacting disclinations. At intermediate times in the coarsening process, split-core dislocation formation and interactions displaying an effective disclination quadrupole configuration are observed. This work provides the framework for further understanding of the formation and dynamics of the diverse set of curvature defects observed in smectic liquid crystals and other layered material systems.

A computational study of the growth of two-dimensional nematic spherulites in an isotropic phase was performed using a Landau-de Gennes type quadrupolar ensor order parameter model for the first-order isotropic/nematic transition of 5CB (pentyl-cyanobiphenyl). An energy balance, taking anisotropy into account, was derived and incorporated into the time-dependent model. Growth laws were determined for two different spherulite morphologies of the form tn, with and without the inclusion of thermal effects. Results show that incorporation of the thermal energy balance correctly predicts the transition of the growth law exponent from the volume driven regime (n=1) to the thermally limited regime (approaching n=0.5), agreeing well with experimental observations. An interfacial nemato-dynamic model is used to gain insight into the interactions that result in the progression of different spherulite growth regimes.