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.

Convective heat dissipation with lattice-frame materials

This paper presents experimental results on heat transfer and pressure drop for a compact heat sink made of fully triangulated, lightweight (porosity∼0.938), aluminum lattice-frame materials (LFMs). Due to the inherent structural anisotropy of the LFMs, two mutually perpendicular orientations were selected for the measurements. Constant heat flux was applied to the heat sink under steady state conditions, and dissipated by forced air convection. The experimental data were compared with those predicted from an analytical model based on fin analogy. The experimental results revealed that pressure drop is strongly dependent upon the orientation of the structure, due mainly to the flow blockage effect. For heat transfer measurements, typical local temperature distributions on the substrate under constant heat flux conditions were captured with infrared camera. The thermal behavior of LFMs was found to follow closely that of cylinder banks, with early transition Reynolds number (based on strut diameter) equal to about 300. The Nusselt number prediction from the fin-analogy correlates well with experimental measurements, except at low Reynolds numbers where a slightly underestimation is observed. Comparisons with empty channels and commonly used heat exchanger media show that the present LFM heat sink can remove heat approximately seven times more efficient than an empty channel and as efficient as a bank of cylinders at the same porosity level. The aluminum LFMs are extremely stiff and strong, making them ideal candidates for multifunctional structures requiring both heat dissipation and mechanical load carrying capabilities. © 2003 Elsevier Ltd. All rights reserved.

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.