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.

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.

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.