Monte Carlo simulation of phase behavior of polymer blends with special interactions

Full Paper: The phase, behavior of A-B-random copolymer/C-homopolymer, blends with special interaction was studied by a. Monte, Carlo simulation in two dimensions. The interaction between I segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. The simulation results showed that the blend became two large co-continuous phase domains at lower segment-B component compositions, indicating that the blend showed spinodal decomposition. With an increase of the segment-B component, the miscibility between the copolymer,and the polymer was gradually improved up to being miscible. In addition, it was found that segment B tended to move to the surface of the copolymer phase in the case of a lower component of segment B. On the other hand, if was observed that the average, end-to-end distances ((h) over bar) for both copolymer and polymer changed slowly with increasing segment-B component of the copolymer up to 40%, thereafter they increased considerably with increasing segment B component. Moreover, it was found that the (h) over bar of the copolymer was obviously shorter than that of the homopolymer for the segment-B composition, region from 0% to 80%. Finally, a, phase diagram showing I phase and - II phase regions under the condition of constant-temperature is presented.

Monte Carlo simulation of phase separation of A/B/A-B ternary mixtures

Monte Carlo simulations were used to model A/B/A-B ternary mixtures with different AB diblock copolymer volume fractions for which both the dispersed and continuous phase volume fractions were kept constant. For concentrations of the diblock copolymer below a critical value, the domain size increment of the dispersed phase decreases linearly with the copolymer concentration. This is in agreement with the predictions of Noolandi and Hong. The dependence of the domain size as a function of the copolymer volume fraction can also be fitted by the equation of Tang and Huang. Our simulations indicate, for the first time, that the micelles form before saturation of the interface occurs. This means that the formation of the micelles is not a result of the saturation of the interface.

Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

Bubble-to-Slug Transition for Low-Re Two-Phase Flows

A critical review on the mechanism and models on the bubble-to-slug transition of two-phase gas-liquid flows are presented in the present paper. It is shown that the most possible mechanism controlled the bubble-to-slug transition is the bubble coalescence. Focusing on the bubble-to-slug transition for the low-Re two-phase flow, a simple Monte Carlo method is used to simulate the influence of the initial bubble size on the bubble-to-slug transition. Some secondary factors, such as the liquid viscosity, the surface tension, and the relative slip between the two phases, are ignored in the present study. It is found that the locus of the dimensionless rate of collision is a universal curve. Based on this curve, it is determined that the bubble initial size can affect the phase distribution and flow pattern when its dimensionless value is in the range from 0.03 to 0.4. A simple relationship between the critical void fraction and the initial bubble size is proposed, which agrees very well with the experimental data.

Information Preservation Modeling of Rayleigh-Benard Transition from Thermal Conduction to Convection

Onset and evolution of the Rayleigh-Benard (R-B) convection are investigated using the Information Preservation (IP) method. The information velocity and temperature are updated using the Octant Flux Splitting (OFS) model developed by Masters & Ye based on the Maxwell transport equation suggested by Sun & Boyd. Statistical noise inherent in particle approaches such as the direct simulation Monte Carlo (DSMC) method is effectively reduced by the IP method, and therefore the evolutions from an initial quiescent fluid to a final steady state are shown clearly. An interesting phenomenon is observed: when the Rayleigh number (Ra) exceeds its critical value, there exists an obvious incubation stage. During the incubation stage, the vortex structure clearly appears and evolves, whereas the Nusselt number (Nu) of the lower plate is close to unity. After the incubation stage, the vortex velocity and Nu rapidly increase, and the flow field quickly reaches a steady, convective state. A relation of Nu to Ra given by IP agrees with those given by DSMC, the classical theory and experimental data.

Rarefied gas dynamics: fundamentals, simulations and micro flows

Kinetic studies of gas flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.

蒙特卡罗模拟四极阱中原子的参变激发

在中性原子的磁囚禁实验中,磁阱线圈的电流噪声会激发磁阱中的原子运动,势必对原子团的温度和寿命产生不可忽视的影响。对于非简谐阱,这种激发具有能量选择特性,它又取决于电流噪声的频谱分布。选择了实验中常用的四极阱为研究对象,用直接模拟蒙特卡罗方法来模拟四极阱中原子运动的参变激发现象,得到了原子温度与原子数损失随激发频率的变化关系,并进一步计算了两个共振峰处原子温度随调制时间和调制深度的变化曲线。此外,还研究了弹性碰撞速率对参变激发过程中原子温度上升的影响。这些结果对四极阱参变激发的实验有较好的参考价值。

光通过纳米颗粒随机散射体透射光强的计算及分析

光通过随机散射体后透射光强的估算,对于提取散射体内部无法直接测量或者无法直接观察的信息至关重要。在分析平行光束通过纳米圆形颗粒随机散射体出射面上光强组成的基础上推导了透射光强的计算公式,并指出多重散射理论、一阶多重散射理论和朗伯比尔定律在一定近似程度上可以相对精确地估算透射光强。这三种方法的估算结果之间会出现四种相对关系。根据这些相对关系,分析了光在散射体内部传输时散射过程的特征以及各种散射过程对出射面光强的贡献大小。

星载激光对水下目标通信可行性研究

The characteristics of media in communication channel are analyzed briefly and the reasonable optical parameters of media are adopted. With certain communication system parameters the temporal and spatial distributions of the received signal from submerged platform are simulated using Monte Carlo method. The upper limit of the ratio of Monte Carlo estimated error to averaged value is about 0.3%. From the simulated results, the optimized sampling timing of receiver and field of view of telescope are obtained. Also the signal-to-noise ratio of the receiver is calculated. Based on this, the error probability of the communication system is deduced from laser pulse position modulation and maximum likelihood detection. The results show that under severe environment robust laser communication from a satellite to a submerged platform can be achieved.

Extended bidirectional reflectance distribution function for polarized light scattering from subsurface defects under a smooth surface

By introducing the scattering probability of a subsurface defect (SSD) and statistical distribution functions of SSD radius, refractive index, and position, we derive an extended bidirectional reflectance distribution function (BRDF) from the Jones scattering matrix. This function is applicable to the calculation for comparison with measurement of polarized light-scattering resulting from a SSD. A numerical calculation of the extended BRDF for the case of p-polarized incident light was performed by means of the Monte Carlo method. Our numerical results indicate that the extended BRDF strongly depends on the light incidence angle, the light scattering angle, and the out-of-plane azimuth angle. We observe a 180 degrees symmetry with respect to the azimuth angle. We further investigate the influence of the SSD density, the substrate refractive index, and the statistical distributions of the SSD radius and refractive index on the extended BRDF. For transparent substrates, we also find the dependence of the extended BRDF on the SSD positions. (c) 2006 Optical Society of America.

The effect of diffusion and overhangs/vacancies on the microstructure of zig-zag thin film

Basing on some growth models of thin film, we have investigated the growth mechanism of glancing angle deposition (GLAD) film. The simulation verifies that the overhangs/vacancies also contribute to the columnar growth as well as the self-shadowing effect for GLAD thin film. Besides, we have studied the effect of the deposition rate, surface and bulk diffusions on the microstructure of thin film using the time-dependent Monte Carlo method. The results show that the surface and bulk diffusions can significantly enhance the packing density of thin film in GLAD growth, and the increase of the deposition rate induce the moderate decrease of the packing density. (c) 2006 Elsevier B.V. All rights reserved.

Bayesian inference of population expansions in domestic bovines

< p > The past population dynamics of four domestic and one wild species of bovine were estimated using Bayesian skyline plots, a coalescent Markov chain Monte Carlo method that does not require an assumed parametric model of demographic history. Four dom

Single-electron multiple-valued memory using ultra-small floating gate metal-oxide-semiconductor field effect transistor

This paper proposes a novel single-electron multiple-valued memory. It is a metal-oxide-semiconductor field effect transistor (MOS)-type memory with multiple separate control gates and floating gate layer, which consists of nano-crystal grains. The electron can tunnel among the grains (floating gates) and between the floating gate layer and the MOS channel. The memory can realize operations of 'write', 'store' and 'erase' of multiple-valued signals exceeding three values by controlling the single electron tunneling behavior. We use Monte Carlo method to simulate the operation of single-electron four-valued memory. The simulation results show that it can operate well at room temperature.