Single curved fiber sedimentation under gravity

Dynamics of single curved fiber sedimentation under gravity are simulated by using the lattice Boltzmann method. The results of migration and rotation of the curved fiber at different Reynolds numbers are reported. The results show that the rotation and migration processes are sensitive to the curvature of the fiber. (c) 2007 Elsevier Ltd. All rights reserved.

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

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

Effect of chain-length dependence of interaction parameter on spinodals for polydisperse polymer blends

The chain-length dependence of the Flory-Huggins (FH) interaction parameter is introduced into the FH lattice theory for polydisperse polymer-blend systems. The spinodals are calculated for the model polymer blends with different chain lengths and distributions. It is found that all the related variables r(n), r(w), r(z), and chain-length distribution, have effects on the spinodals for polydisperse polymer blends.

Theoretical representation of shear and pressure influence on the phase separation behavior of PVME/PS mixture

In the framework of lattice fluid model, the Gibbs energy and equation of state are derived by introducing the energy (E-s) stored during flow for polymer blends under shear. From the calculation of the spinodal of poly(vinyl methyl ether) (PVME) and polystyrene (PS) mixtures, we have found the influence of E., an equation of state in pure component is inappreciable, but it is appreciable in the mixture. However, the effect of E, on phase separation behavior is extremely striking. In the calculation of spinodal for the PVME/PS system, a thin, long and banana miscibility gap generated by shear is seen beside the miscibility gap with lower critical solution temperature. Meanwhile, a binodal coalescence of upper and lower miscibility gaps is occurred. The three points of the three-phase equilibrium are forecasted. The shear rate dependence of cloud point temperature at a certain composition is discussed. The calculated results are acceptable compared with the experiment values obtained by Higgins et at. However, the maximum positive shift and the minimum negative shift of cloud point temperature guessed by Higgins are not obtained, Furthermore, the combining effects of pressure and shear on spinodal shift are predicted.

STATISTICAL THERMODYNAMICS OF POLYDISPERSE POLYMER BLENDS WITH STRONG INTERACTION

A statistical thermodynamics theory of polydisperse polymer mixtures with strong interaction between dissimilar components based on a lattice fluid model is formulated. Expressions for the free energy, equation of state, phase stability and spinodal for a polydisperse, binary polymer mixture with strong interaction are derived.

Lattice type of fracture model for concrete

Concrete is usually described as a three-phase material, where matrix, aggregate and interface zones are distinguished. The beam lattice model has been applied widely by many investigators to simulate fracture processes in concrete. Due to the extremely large computational effort, however, the beam lattice model faces practical difficulties. In our investigation, a new lattice called generalized beam (GB) lattice is developed to reduce computational effort. Numerical experiments conducted on a panel subjected to uniaxial tension show that the GB lattice model can reproduce the load-displacement curves and crack patterns in agreement to what are observed in tests. Moreover, the effects of the particle overlay on the fracture process are discussed in detail. (C) 2007 Elsevier Ltd. All rights reserved.

Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

The study of banded spherulite patterns by coupled logistic map lattice model

Banded spherulite patterns are simulated in two dimensions by means of a coupled logistic map lattice model. Both target pattern and spiral pattern which have been proved to be existent experimentally in banded spherulite are obtained by choosing suitable parameters in the model. The simulation results also indicate that the band spacing is decreased with the increase of parameter mu in the logistic map and increased with the increase of the coupling parameter epsilon, which is quite similar to the results in some experiments. Moreover, the relationship between the parameters and the corresponding patterns is obtained, and the target patterns and spiral patterns are distinguished for a given group of initial values, which may guide the study of banded spherulite.

Simulation of Banded Spherulite Patterns by Coupled Logistic Map Lattice Model in Three Dimensions

Banded spherulite patterns are simulated in three dimensions by means of a Coupled Logistic map lattice model. The patterns obtained by numerical calculation are consistent with those in experiments. The simulation results also indicate that the hand spacing is decreased with the increase of parameter mu in the Logistic map and increased with the increase of the coupling parameter e for cube lattices, and increased with the increase of the thickness of the lattice for polymer film, which is quite similar to the results in some experiments. Spiral pattern in three dimensions is also shown in this paper, which helps us understand the form of banded spherulite in polymers.

A novel self-consistent-field lattice model for block copolymers

We develop a self-consistent-field lattice model for block copolymers and propose a novel and general method to solve the self-consistent-field equations. The approach involves describing the polymer chains in a lattice and employing a two-stage relaxation procedure to evolve a system as rapidly as possible to a free-energy minimum. In order to test the validity of this approach, we use the method to study the microphases of rod-coil diblock copolymers. In addition to the lamellar and cylindrical morphologies, micellar, perforated lamellar, gyroid, and zigzag structures have been identified without any prior assumption of the microphase symmetry. Furthermore, this approach can also give the possible orientation of the rods in different structures.

Application of the Sanchez-Lacombe lattice fluid theory to the system pvme/ps and model calculations

Cloud point curves reported in the literature for five representatives of the system poly(vinyl methyl ether)/polystyrene were evaluated theoretically by means of the Sanchez-Lacombe lattice fluid theory. The measured phase separation behavior can be described within experimental error using only one adjustable parameter (quantifying the interaction between the unlike mers). The Flory-Huggins interaction parameters calculated from this theoretical description depend in good approximation linearly on composition (volume fractions) and on the inverse temperature. An evaluation of these data yields a maximum heat effect which is almost one order of magnitude less (ca. -0.25 J/cm(3)) than obtained via Hess's cycle (dissolution of the components and of the blend) from calorimetric measurements. Model calculations on the basis of the present theory demonstrate that the critical points shift to a different extent upon a certain relative change in the molar mass of the blend components. The sensitivity of the calculated phase diagrams against changes in the scaling parameter decreases in the following order: interaction energies between unlike mers, differences in the scaling temperatures, pressures and densities.

Lattice instability at a fast moving crack tip

A molecular dynamics method is used to analyze the dynamic propagation of an atomistic crack tip. The simulation shows that the crack propagates at a relatively constant global velocity which is well below the Rayleigh wave velocity. However the local propagation velocity oscillates violently, and it is limited by the longitudinal wave velocity. The crack velocity oscillation is caused by a repeated process of crack tip blunting and sharpening. When the crack tip opening displacement exceeds a certain critical value, a lattice instability takes place and results in dislocation emissions from the crack tip. Based on this concept, a criterion for dislocation emission from a moving crack tip is proposed. The simulation also identifies the emitted dislocation as a source for microcrack nucleation. A simple method is used to examine this nucleation process. (C) 1996 American Institute of Physics.

Phase relaxation model for 2-phase flows

A phase relaxation model (PRM) for 2-phase flows is presented in this paper on the basis of three principal assumptions. The basic equations for PRM arc derived from the Boltzmann equations for gas-partlcle mixture, The general characteristics and solving process of the PRM's basic equations are also presented and discussed. Many terms in the PRM's basic equations contain a factor ε= ρgρp/ρg+ρp2 which is an intrinsic small parameter for 2-phase mixture, with ρg and ρp being respectively the densities of gas and particle phases.This makes it possible to simplify the computation of the PRM's basic equations. The model is applied to for example, studying file steady propagation of shock waves in gas-particle mixture. The analysis shows that with an increase of shock wave strength the relaxation process behind a gasdynamics shock front becomes a kind of dynamics relaxation instead of the standard exponential relaxation process. A method of determining experimentally the velocity and tem...更多perature relaxation rates (or times) of gas-particle flows is suggested and analyzed.

Algorithm for simulating fracture processes in concrete by lattice modeling

To simulate fracture behaviors in concrete more realistically, a theoretical analysis on the potential question in the quasi-static method is presented, then a novel algorithm is proposed which takes into account the inertia effect due to unstable crack propagation and meanwhile requests much lower computational efforts than purely dynamic method. The inertia effect due to load increasing becomes less important and can be ignored with the loading rate decreasing, but the inertia effect due to unstable crack propagation remains considerable no matter how low the loading rate is. Therefore, results may become questionable if a fracture process including unstable cracking is simulated by the quasi-static procedure excluding completely inertia effects. However, it requires much higher computational effort to simulate experiments with not very high loading rates by the dynamic method. In this investigation which can be taken as a natural continuation, the potential question of quasi-static method is analyzed based on the dynamic equations of motion. One solution to this question is the new algorithm mentioned above. Numerical examples are provided by the generalized beam (GB) lattice model to show both fracture processes under different loading rates and capability of the new algorithm.

A simple method to simulate shrinkage-induced cracking in cement-based composites by lattice-type modeling

A new numerical procedure is proposed to investigate cracking behaviors induced by mismatch between the matrix phase and aggregates due to matrix shrinkage in cement-based composites. This kind of failure processes is simplified in this investigation as a purely spontaneous mechanical problem, therefore, one main difficulty during simulating the phenomenon lies that no explicit external load serves as the drive to propel development of this physical process. As a result, it is different from classical mechanical problems and seems hard to be solved by using directly the classical finite element method (FEM), a typical kind of "load -> medium -> response" procedures. As a solution, the actual mismatch deformation field is decomposed into two virtual fields, both of which can be obtained by the classical FEM. Then the actual response is obtained by adding together the two virtual displacement fields based on the principle of superposition. Then, critical elements are detected successively by the event-by-event technique. The micro-structure of composites is implemented by employing the generalized beam (GB) lattice model. Numerical examples are given to show the effectiveness of the method, and detailed discussions are conducted on influences of material properties.