Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.

On the unstable stacking criterion for ideal and cracked copper-crystals

The unstable stacking criteria for an ideal copper crystal under homogeneous shearing and for a cracked copper crystal under pure mode II loading are analysed. For the ideal crystal under homogeneous shearing, the unstable stacking energy gamma(us) defined by Rice in 1992 results from shear with no relaxation in the direction normal to the slip plane. For the relaxed shear configuration, the critical condition for unstable stacking does not correspond to the relative displacement Delta = b(p)/2, where b(p) is the Burgers vector magnitude of the Shockley partial dislocation, but to the maximum shear stress. Based on this result, the unstable stacking energy Gamma(us) is defined for the relaxed lattice. For the cracked crystal under pure mode II loading, the dislocation configuration corresponding to Delta = b(p)/2 is a stable state and no instability occurs during the process of dislocation nucleation. The instability takes place at approximately Delta = 3b(p)/4. An unstable stacking energy Pi(us) is defined which corresponds to the unstable stacking state at which the dislocation emission takes place. A molecular dynamics method is applied to study this in an atomistic model and the results verify the analysis above.

A lattice dynamical treatment for the total potential energy of single-walled carbon nanotubes and its applications: relaxed equilibrium structure, elastic properties, and vibrational modes of ultra-narrow tubes

In this paper, we propose a lattice dynamic treatment for the total potential energy of single-walled carbon nanotubes (SWCNTs) which is, apart from a parameter for the nonlinear effects, extracted from the vibrational energy of the planar graphene sheet. The energetics, elasticity and lattice dynamics are treated in terms of the same set of force constants, independently of the tube structures. Based upon this proposal, we have investigated systematically the relaxed lattice configuration for narrow SWCNTs, the strain energy, the Young's modulus and Poisson ratio, and the lattice vibrational properties with respect to the relaxed equilibrium tubule structure. Our calculated results for various physical quantities are nicely in consistency with existing experimental measurements. In particular, we verified that the relaxation effect makes the bond length longer and the frequencies of various optical vibrational modes softer. Our calculation provides evidence that the Young's modulus of an armchair tube exceeds that of the planar graphene sheet, and that the large diameter limits of the Young's modulus and Poisson ratio are in agreement with the experimental values of graphite; the calculated radial breathing modes for ultra-narrow tubes with diameters ranging between 2 and 5 angstrom coincide with the experimental results and the existing ab initio calculations with satisfaction. For narrow tubes with a diameter of 20 angstrom, the calculated frequencies of optical modes in the tubule's tangential plane, as well as those of radial breathing modes, are also in good agreement with the experimental measurements. In addition, our calculation shows that various physical quantities of relaxed SWCNTs can actually be expanded in terms of the chiral angle defined for the corresponding ideal SWCNTs.

Simple lattice Boltzmann model for traffic flows

A lattice Boltzmann model with 5-bit lattice for traffic flows is proposed. Using the Chapman-Enskog expansion and multi-scale technique, we obtain the higher-order moments of equilibrium distribution function. A simple traffic light problem is simulated by using the present lattice Boltzmann model, and the result agrees well with analytical solution.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Lattice Boltzmann方法模拟多孔介质Klinkenberg效应

格子Boltzmann数值模拟方法是研究复杂的多孔介质结构特别是Klinkenberg效应的有效方法之一，对处理复杂边值问题尤其有效。用格子Boltzmann方法研究了气流穿越多孔介质问题，并将数值计算结果与实验结果进行了比较，结果表明格子Boltzmann方法是数值模拟气流穿越多孔介质问题的有效方法之一。

A lattice Boltzmann method for KDV equation

We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.

Lattice Boltzmann method for compressible flows with high Mach numbers

In this paper we present a lattice Boltzmann model to simulate compressible flows by introducing an attractive force. This scheme has two main advantages: one is to soften sound speed effectively, which greatly raises the Mach number (up to 5); another is its relative simple procedure. Simulations of the March cone and the comparison between theoretical expectations and simulations demonstrate that the scheme is effective in the simulation of compressible flows with high Mach numbers, which would create many new applications.

A lattice Boltzmann equation for waves

We propose a lattice Boltzmann model for the wave equation. Using a lattice Boltzmann equation and the Chapman-Enskog expansion, we get 1D and 2D wave equations with truncation error of order two. The numerical tests show the method can be used to simulate the wave motions.

Simulation of the Load-Unload Response Ratio and Critical Sensitivity in the Lattice Solid Model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

A Lagrangian lattice Boltzmann method for Euler equations

A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.

Recovery of the solitons using a lattice Boltzmann model

We formulate a lattice Boltzmann model which simulates Korteweg-de Vries equation by using a method of higher moments of lattice Boltzmann equation. Using a series of lattice Boltzmann equations in different time scales and the conservation law in time scale to, we obtain equilibrium distribution function. The numerical examples show that the method can be used to simulate soliton.

Heteroclinic cycles in lattice dynamical systems

A criterion of spatial chaos occurring in lattice dynamical systems-heteroclinic cycle-is discussed. It is proved that if the system has asymptotically stable heteroclinic cycle, then it has asymptotically stable homoclinic point which implies spatial chaos.

Modified Generalized Beam Lattice Model Associated With Fracture Of Reinforced Fiber/Particle Composites

Lattice-type model can simulate in a straightforward manner heterogeneous brittle media. Mohr-Coulomb failure criterion has recently been involved into the generalized beam (GB) lattice model, and as a result, numerical experiments on concrete under various loading conditions can be conducted. The GB lattice model is further used to investigate the reinforced fiber/particle composites instead of only particle composites as the model did before. Numerical examples are given to show the effectiveness of the modeling procedure, and influences of inclusions (particle, fiber and rebar) on the fracture processes are also discussed. (c) 2008 Elsevier Ltd. All rights reserved.

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.