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.

Development of a novel pyrolysis-gas chromatography/mass spectrometry method for the analysis of poly(lactic acid) thermal degradation products

Thermal degradation of PLA is a complex process since it comprises many simultaneous reactions. The use of analytical techniques, such as differential scanning calorimetry (DSC) and thermogravimetry (TGA), yields useful information but a more sensitive analytical technique would be necessary to identify and quantify the PLA degradation products. In this work the thermal degradation of PLA at high temperatures was studied by using a pyrolyzer coupled to a gas chromatograph with mass spectrometry detection (Py-GC/MS). Pyrolysis conditions (temperature and time) were optimized in order to obtain an adequate chromatographic separation of the compounds formed during heating. The best resolution of chromatographic peaks was obtained by pyrolyzing the material from room temperature to 600 °C during 0.5 s. These conditions allowed identifying and quantifying the major compounds produced during the PLA thermal degradation in inert atmosphere. The strategy followed to select these operation parameters was by using sequential pyrolysis based on the adaptation of mathematical models. By application of this strategy it was demonstrated that PLA is degraded at high temperatures by following a non-linear behaviour. The application of logistic and Boltzmann models leads to good fittings to the experimental results, despite the Boltzmann model provided the best approach to calculate the time at which 50% of PLA was degraded. In conclusion, the Boltzmann method can be applied as a tool for simulating the PLA thermal degradation.

Finite difference calculations of permeability in large domains in a wide porosity range.

Determining effective hydraulic, thermal, mechanical and electrical properties of porous materials by means of classical physical experiments is often time-consuming and expensive. Thus, accurate numerical calculations of material properties are of increasing interest in geophysical, manufacturing, bio-mechanical and environmental applications, among other fields. Characteristic material properties (e.g. intrinsic permeability, thermal conductivity and elastic moduli) depend on morphological details on the porescale such as shape and size of pores and pore throats or cracks. To obtain reliable predictions of these properties it is necessary to perform numerical analyses of sufficiently large unit cells. Such representative volume elements require optimized numerical simulation techniques. Current state-of-the-art simulation tools to calculate effective permeabilities of porous materials are based on various methods, e.g. lattice Boltzmann, finite volumes or explicit jump Stokes methods. All approaches still have limitations in the maximum size of the simulation domain. In response to these deficits of the well-established methods we propose an efficient and reliable numerical method which allows to calculate intrinsic permeabilities directly from voxel-based data obtained from 3D imaging techniques like X-ray microtomography. We present a modelling framework based on a parallel finite differences solver, allowing the calculation of large domains with relative low computing requirements (i.e. desktop computers). The presented method is validated in a diverse selection of materials, obtaining accurate results for a large range of porosities, wider than the ranges previously reported. Ongoing work includes the estimation of other effective properties of porous media.

Multiscale modeling of lamellar mesophases

The mesoscale simulation of a lamellar mesophase based on a free energy functional is examined with the objective of determining the relationship between the parameters in the model and molecular parameters. Attention is restricted to a symmetric lamellar phase with equal volumes of hydrophilic and hydrophobic components. Apart from the lamellar spacing, there are two parameters in the free energy functional. One of the parameters, r, determines the sharpness of the interface, and it is shown how this parameter can be obtained from the interface profile in a molecular simulation. The other parameter, A, provides an energy scale. Analytical expressions are derived to relate these parameters to r and A to the bending and compression moduli and the permeation constant in the macroscopic equation to the Onsager coefficient in the concentration diffusion equation. The linear hydrodynamic response predicted by the theory is verified by carrying out a mesoscale simulation using the lattice-Boltzmann technique and verifying that the analytical predictions are in agreement with simulation results. A macroscale model based on the layer thickness field and the layer normal field is proposed, and the relationship between the parameters in the macroscale model from the parameters in the mesoscale free energy functional is obtained.

Armchair or Zigzag? A tool for characterizing graphene edge

Electronic, magnetic, and structural properties of graphene flakes depend sensitively upon the type of edge atoms. We present a simple software tool for determining the type of edge atoms in a honeycomb lattice. The algorithm is based on nearest neighbor counting. Whether an edge atom is of armchair or zigzag type is decided by the unique pattern of its nearest neighbors. Particular attention is paid to the practical aspects of using the tool, as additional features such as extracting out the edges from the lattice could help in analyzing images from transmission microscopy or other experimental probes. Ultimately, the tool in combination with density-functional theory or tight-binding method can also be helpful in correlating the properties of graphene flakes with the different armchair-to-zigzag ratios. Program summary Program title: edgecount Catalogue identifier: AEIA_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEIA_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 66685 No. of bytes in distributed program, including test data, etc.: 485 381 Distribution format: tar.gz Programming language: FORTRAN 90/95 Computer: Most UNIX-based platforms Operating system: Linux, Mac OS Classification: 16.1, 7.8 Nature of problem: Detection and classification of edge atoms in a finite patch of honeycomb lattice. Solution method: Build nearest neighbor (NN) list; assign types to edge atoms on the basis of their NN pattern. Running time: Typically similar to second(s) for all examples. (C) 2010 Elsevier B.V. All rights reserved.

Shear alignment of a disordered lamellar mesophase

The shear alignment of an initially disordered lamellar phase is examined using lattice Boltzmann simulations of a mesoscopic model based on a free-energy functional for the concentration modulation. For a small shear cell of width 8 lambda, the qualitative features of the alignment process are strongly dependent on the Schmidt number Sc = nu/D (ratio of kinematic viscosity and mass diffusion coefficient). Here, lambda is the wavelength of the concentration modulation. At low Schmidt number, it is found that there is a significant initial increase in the viscosity, coinciding with the alignment of layers along the extensional axis, followed by a decrease at long times due to the alignment along the flow direction. At high Schmidt number, alignment takes place due to the breakage and reformation of layers because diffusion is slow compared to shear deformation; this results in faster alignment. The system size has a strong effect on the alignment process; perfect alignment takes place for a small systems of width 8 lambda and 16 lambda, while a larger system of width 32 lambda does not align completely even at long times. In the larger system, there appears to be a dynamical steady state in which the layers are not perfectly aligned-where there is a balance between the annealing of defects due to shear and the creation due to an instability of the aligned lamellar phase under shear. We observe two types of defect creation mechanisms: the buckling instability under dilation, which was reported earlier, as well as a second mechanism due to layer compression.

High thermopower and ultra low thermal conductivity in Cd-based Zintl phase compounds

By combining first principles density functional theory and electronic as well as lattice Boltzmann transport calculations, we unravel the excellent thermoelectric properties of Zintl phase compounds ACd(2)Sb(2) (where, A = Ca, Ba, Sr). The calculated electronic structures of these compounds show charge carrier pockets and heavy light bands near the band edge, which lead to a large power factor. Furthermore, we report large Gruneisen parameters and low phonon group velocity indicating essential strong anharmonicity in these compounds, which resulted in low lattice thermal conductivity. The combination of low thermal conductivity and the excellent transport properties give a high ZT value of similar to 1.4-1.9 in CaCd2Sb2 and BaCd2Sb2 at moderate p and n-type doping. Our results indicate that well optimized Cd-based Zintl phase compounds have the potential to match the performance of conventional thermoelectric materials.

Role of the Molar Volume on Estimated Diffusion Coefficients

The role of the molar volume on the estimated diffusion parameters has been speculated for decades. The Matano-Boltzmann method was the first to be developed for the estimation of the variation of the interdiffusion coefficients with composition. However, this could be used only when the molar volume varies ideally or remains constant. Although there are no such systems, this method is still being used to consider the ideal variation. More efficient methods were developed by Sauer-Freise, Den Broeder, and Wagner to tackle this problem. However, there is a lack of research indicating the most efficient method. We have shown that Wagner's method is the most suitable one when the molar volume deviates from the ideal value. Similarly, there are two methods for the estimation of the ratio of intrinsic diffusion coefficients at the Kirkendall marker plane proposed by Heumann and van Loo. The Heumann method, like the Matano-Boltzmann method, is suitable to use only when the molar volume varies more or less ideally or remains constant. In most of the real systems, where molar volume deviates from the ideality, it is safe to use the van Loo method. We have shown that the Heumann method introduces large errors even for a very small deviation of the molar volume from the ideal value. On the other hand, the van Loo method is relatively less sensitive to it. Overall, the estimation of the intrinsic diffusion coefficient is more sensitive than the interdiffusion coefficient.

模拟MEMS中气体流动的工具——信息保存法

MEMS中气体流动因特征尺度小而是稀薄气体的领域,本文首先介绍处理低速稀薄气体流动的一些方法:线化Boltzmann方程,Lattice Boltzmann方法,加滑移边界的Navier-Stoks方程,以及DSMC方法,并讨论它们模拟MEMS中过渡领域低速流动所遇到的困难。信息保存法克服了流速低使得信息噪声比小引起统计模拟的困难,已成功模拟了一些一维和二维问题。MEMS中流速低和大的长宽比的特点还引起出入口边界条件相互影响需要协调的问题,通过微槽道流动的算例,在模拟中采用守恒形式的质量守恒方程和超松弛法成功地解决了这一问题。处理有温度变化的MEMS流动问题和跨越领域的混合算法是重要的问题,本文用信息保存法也进行了有益的尝试。

Quantum Simulator for Transport Phenomena in Fluid Flows

Transport phenomena still stand as one of the most challenging problems in computational physics. By exploiting the analogies between Dirac and lattice Boltzmann equations, we develop a quantum simulator based on pseudospin-boson quantum systems, which is suitable for encoding fluid dynamics transport phenomena within a lattice kinetic formalism. It is shown that both the streaming and collision processes of lattice Boltzmann dynamics can be implemented with controlled quantum operations, using a heralded quantum protocol to encode non-unitary scattering processes. The proposed simulator is amenable to realization in controlled quantum platforms, such as ion-trap quantum computers or circuit quantum electrodynamics processors.

Mixing and internal dynamics of droplets impacting and coalescing on a solid surface.

The coalescence and mixing of a sessile and an impacting liquid droplet on a solid surface are studied experimentally and numerically in terms of lateral separation and droplet speed. Two droplet generators are used to produce differently colored droplets. Two high-speed imaging systems are used to investigate the impact and coalescence of the droplets in color from a side view with a simultaneous gray-scale view from below. Millimeter-sized droplets were used with dynamical conditions, based on the Reynolds and Weber numbers, relevant to microfluidics and commercial inkjet printing. Experimental measurements of advancing and receding static contact angles are used to calibrate a contact angle hysteresis model within a lattice Boltzmann framework, which is shown to capture the observed dynamics qualitatively and the final droplet configuration quantitatively. Our results show that no detectable mixing occurs during impact and coalescence of similar-sized droplets, but when the sessile droplet is sufficiently larger than the impacting droplet vortex ring generation can be observed. Finally we show how a gradient of wettability on the substrate can potentially enhance mixing.

Theoretical investigation on the phase stability of Nd2Fe14B and site preference of V, Cr, Mn, Zr and Nb

The stability of the excellent permanent magnetic compound Nd2Fe14B and substitution of Fe in the compound by V, Cr, Mn, Zr and Nb are investigated by using interatomic pair potentials which are converted from lattice-inversion method. Calculation shows that the substitution always makes the cell volume larger, and the increase of the volume is almost linear with substituent concentration. The calculated cohesive energy shows that the preferential order of substitution of Fe is Nb, V, Cr, Mn, Zr. Nevertheless, all the five substituting elements should most preferentially replace Fe in the j(2)' site, which has the greatest space among all six Fe sites. (C) 2005 Elsevier B.V. All rights reserved.

The behavior of the sextic coupling for the scalar field at the intermediate and strong coupling regime

We study the behavior of the renormalized sextic coupling at the intermediate and strong coupling regime for the phi(4) theory defined in d = 2 dimensions. We found a good agreement with the results obtained by the field-theoretical renormalization-group in the Ising limit. In this work we use the lattice regularization method.

Noise and ultraviolet divergences in the dynamics of the chiral condensate in QCD

The time evolution of the matter produced in high energy heavy-ion collisions seems to be well described by relativistic viscous hydrodynamics. In addition to the hydrodynamic degrees of freedom related to energy-momentum conservation, degrees of freedom associated with order parameters of broken continuous symmetries must be considered because they are all coupled to each other. of particular interest is the coupling of degrees of freedom associated with the chiral symmetry of QCD. Quantum and thermal fluctuations of the chiral fields act as noise sources in the classical equations of motion, turning them into stochastic differential equations in the form of Ginzburg-Landau-Langevin (GLL) equations. Analytic solutions of GLL equations are attainable only in very special circumstances and extensive numerical simulations are necessary, usually by discretizing the equations on a spatial lattice. However, a not much appreciated issue in the numerical simulations of GLL equations is that ultraviolet divergences in the form of lattice-spacing dependence plague the solutions. The divergences are related to the well-known Rayleigh-Jeans catastrophe in classical field theory. In the present communication we present a systematic lattice renormalization method to control the catastrophe. We discuss the implementation of the method for a GLL equation derived in the context of a model for the QCD chiral phase transition and consider the nonequilibrium evolution of the chiral condensate during the hydrodynamic flow of the quark-gluon plasma.

Determinação de temperatura de chamas por espectroscopia de espécies quimioluminescentes

In a combustion process involving fossil fuels, there is the formation of species Chemiluminescent, especially CH*, C2* and OH*, whose spontaneous emission can be used as a diagnostic tool. In the present work, mapping and determination of the rotational temperature of the species CH* produced in flames on a burner fueled by Liquefied Petroleum Gas (LPG) was carried out. This study is part of a project involving the characterization of supersonic combustion in scramjets engines, whose study has been conducted in the hypersonic shock tunnel IEAv laboratories. The technique used was the natural emission spectroscopy, which has as main advantage of being non-intrusive. The rotational temperature determination was made using the Boltzmann method, whose principle is to relate the emission intensity of the species to the temperature by means of spectroscopic constants established.The temperature values were determined from the analysis of electronic bands AX and BX of the radical CH*. In order to confirm the results of flame temperatures obtained by the natural emission technique, was also used the technique of line reversal sodium. The results of both techniques showed that the temperature of the flames investigated is about 2500K a 2700K